Research March 2017  Volume 125  Issue 3
Mediation of the Relationship between Maternal Phthalate Exposure and Preterm Birth by Oxidative Stress with Repeated Measurements across Pregnancy
Kelly K. Ferguson,^{1} YinHsiu Chen,^{2} Tyler J. VanderWeele,^{3} Thomas F. McElrath,^{4} John D. Meeker,^{1} and Bhramar Mukherjee^{2}
PDF Version (301 KB)

Background: Mediation analysis is useful for understanding mechanisms and has been used minimally in the study of the environment and disease.
Objective: We examined mediation of the association between phthalate exposure during pregnancy and preterm birth by oxidative stress.
Methods: This nested case–control study of preterm birth (n = 130 cases, 352 controls) included women who delivered in Boston, Massachusestts, from 2006 through 2008. Phthalate metabolites and 8isoprostane, an oxidative stress biomarker, were measured in urine from three visits in pregnancy. We applied four counterfactual mediation methods: method 1, utilizing exposure and mediator averages; method 2, using averages but allowing for an exposure–mediator interaction; method 3, incorporating longitudinal measurements of the exposure and mediator; and method 4, using longitudinal measurements and allowing for an exposure–mediator interaction.
Results: We observed mediation of the associations between phthalate metabolites and all preterm birth by 8isoprostane, with the greatest estimated proportion mediated observed for spontaneous preterm births specifically. Fully utilizing repeated measures of the exposure and mediator improved precision of indirect (i.e., mediated) effect estimates, and including an exposure–mediator interaction increased the estimated proportion mediated. For example, for mono(2ethylcarboxypropyl) phthalate (MECPP), a metabolite of di(2ethylhexyl) phthalate (DEHP), the percent of the total effect mediated by 8isoprostane increased from 47% to 60% with inclusion of an exposure–mediator interaction term, in reference to a total adjusted odds ratio of 1.67 or 1.48, respectively.
Conclusions: This demonstrates mediation of the phthalate–preterm birth relationship by oxidative stress, and the utility of complex regression models in capturing mediated associations when repeated measures of exposure and mediator are available and an exposure–mediator interaction may exist.

Citation: Ferguson KK, Chen YH, VanderWeele TJ, McElrath TF, Meeker JD, Mukherjee B. 2017. Mediation of the relationship between maternal phthalate exposure and preterm birth by oxidative stress with repeated measurements across pregnancy. Environ Health Perspect 125:488–494; http://dx.doi.org/10.1289/EHP282
Address correspondence to B. Mukherjee, University of Michigan School of Public Health, Department of Biostatistics, M4166 SPH II, 1415 Washington Heights, Ann Arbor, MI 481092029 USA. Phone: (734) 7646544. Email: bhramar@umich.edu
We thank K. Kneen, S. Clipper, G. Pace, D. Weller, and J. Bell of NSF International in Ann Arbor, Michigan, for phthalate analysis; E. Hurst and colleagues of Cayman Chemical in Ann Arbor, Michigan, for oxidative stress biomarker analysis; and R. Hauser of Harvard School of Public Health for guidance on study design.
Funding for the study was provided by the National Institute of Environmental Health Sciences, National Institutes of Health (R01ES018872, P42ES017198, P01ES022844, and P30ES017885).
The authors declare they have no actual or potential competing financial interests.
Received: 16 September 2015
Revised: 31 March 2016
Accepted: 27 May 2016
Published: 28 June 2016Note to readers with disabilities: EHP strives to ensure that all journal content is accessible to all readers. However, some figures and Supplemental Material published in EHP articles may not conform to 508 standards due to the complexity of the information being presented. If you need assistance accessing journal content, please contact ehponline@niehs.nih.gov. Our staff will work with you to assess and meet your accessibility needs within 3 working days.
 Supplemental Material PDF (514 KB)
Note to readers with disabilities: EHP has provided a 508conformant table of contents summarizing the Supplemental Material for this article (see below) so readers with disabilities may determine whether they wish to access the full, nonconformant Supplemental Material. If you need assistance accessing journal content, please contact ehponline@niehs.nih.gov. Our staff will work with you to assess and meet your accessibility needs within 3 working days.
Supplemental Table of Contents PDF (116 KB)
Introduction
Understanding causal mechanism in observational studies of environmental exposure and complex disease is challenging. A useful approach may be to screen human populations for biomarkers of exposures as well as mechanistic intermediates and assess relationships with mediation analyses to aid in establishing biological pathways. We recently observed that urinary concentrations of phthalate metabolites were associated with increased odds of preterm birth (Ferguson et al. 2014). Phthalate metabolites are indicative of exposure to phthalate diesters which are found ubiquitously in the environment in plastics, personal care products, and medications (ATSDR 2001, 2002). Exposure to these compounds during pregnancy has been linked to preterm birth in other previous studies as well, although some results are conflicting (Adibi et al. 2009; Meeker et al. 2009; Whyatt et al. 2009; Wolff et al. 2008). We hypothesized that this relationship may be mediated by phthalateinduced maternal oxidative stress. Our previous work circumstantially supports this hypothesis; we have established that urinary phthalate metabolites are associated with an increase in urinary 8isoprostane, a biomarker of oxidative stress (Ferguson et al. 2015b), and that furthermore urinary 8isoprostane is associated with an increased risk of preterm birth (Ferguson et al. 2015a). In the present analysis, we sought to quantify and test the role of 8isoprostane in the relationship between maternal phthalate exposure during pregnancy and prematurity using a mediation analysis within a causal framework.
Mediation analysis has been used in social and epidemiological research for decades to understand causal pathways, biological mechanisms, and to design policy interventions (Pearl 2014). Its development originated in structural equation modeling, and with the work of Baron and Kenny (1986) became widely used in the context of linear models. Application of a counterfactual framework (Rubin 1974) has further provided a strong theoretical basis for causal inference from mediation analysis by precisely defining the necessary assumptions. This framework has also facilitated natural extension of mediation analysis to more complex models that include exposure–mediator interaction, nonlinear terms, and longitudinal data (Valeri and VanderWeele 2013; VanderWeele and Tchetgen Tchetgen 2014; VanderWeele and Vansteelandt 2010).
The use of mediation analysis in environmental and exposure epidemiology has been infrequent and more recent. For example, in the Normative Aging Study partial mediation of the relationship between black carbon particulate exposure and increased fibrinogen levels was observed through epigenetic modifications (Bind et al. 2014). The limited application of mediation analysis is surprising, because many studies in environmental epidemiology measure markers of intermediate biological changes (e.g., hormone levels, epigenetic modifications, and inflammatory cytokines) in addition to examining exposure and disease outcomes.
The goal of the present study was to examine causal mediation in an observational study of preterm birth. Specifically, we investigated mediation of the relationship between exposure to phthalates, as indicated by urinary phthalate metabolites, and preterm birth by oxidative stress, as indicated by urinary 8isoprostane. Within the causal framework, we investigated the natural direct effect (NDE) and the natural indirect (i.e., mediated) effect (NIE), as illustrated in Figure 1. Statistically speaking, the NDE refers to the change in the odds of preterm birth (Y) in association with a defined change in the urinary phthalate concentration (A; e.g., from a to a*) while holding the urinary 8isoprostane concentration (M) at the level it would have naturally been at with A set at the original level (e.g., a). The NIE refers to the change in the odds of preterm birth (Y) when urinary phthalate concentration is held at a specific level (e.g., a) and 8isoprostane is set at what it would have naturally been at for the defined change in urinary phthalate concentration (e.g., a*) (Pearl 2014).
Figure 1. Simple conceptual model for mediation analysis in the context of the present study (adapted from VanderWeele and Tchetgen Tchetgen 2014).
In addition to applying this counterfactual framework to allow for causal interpretation of our results, we wanted to appropriately use the repeated measures available in this data set to fully capture mediated effects. There are very few empirical examples of mediation analysis with longitudinal data in the literature (MacKinnon 2012), and none in a causal framework with timevarying exposure and mediators and a non–timevarying binary end point. Thus, we applied four analytic strategies to this study using recent advancements in the causal inference literature for mediation analysis, using the counterfactual framework. These include method 1, an approach utilizing exposure and mediator averages; method 2, using averages but allowing for an exposure–mediator interaction; method 3, a longitudinal approach utilizing repeated measures of exposure and mediator; and method 4, using repeated measures and allowing for an exposure–mediator interaction.
Methods
Study Population
The study population has been described in detail previously (Ferguson et al. 2014). Briefly, mothers included in this nested case–control study were selected from a prospective birth cohort examining predictors of preeclampsia in women who delivered at the Brigham and Women’s Hospital in Boston, Massachusetts, from 2006 through 2008. The present study, designed with the primary purpose of examining phthalate exposure in relation to preterm birth, included 130 women who delivered preterm as well as 352 randomly selected controls. Participants provided written informed consent, and institutional review board approval was obtained from Brigham and Women’s Hospital as well as the University of Michigan. Preterm birth was defined as delivery before 37 weeks gestation based on last menstrual period with verification by firsttrimester ultrasound (hereafter denoted all preterm birth). We also examined spontaneous preterm births alone (delivery < 37 weeks gestation preceded by spontaneous preterm labor and/or preterm premature rupture of membranes), because we had previously observed strongest associations within this subgroup and we hypothesized that oxidative stress would play a more important role in this pathway based on biologic plausibility (Ferguson et al. 2014).
Urine samples were available for analysis of phthalate and oxidative stress biomarkers from up to four visits per subject (mean, 3.52 visits per subject) across gestation, at median 10, 18, 26, and 35 weeks gestation. For the present analysis, we used measures from visits 1–3 only because many of the preterm cases had delivered by the 35week visit and those measurements may be biased. Demographic characteristics that were included in regression models in previous analyses within this population included maternal age at visit 1, race/ethnicity (white, African American, other), education level (high school, technical school, junior college/some college, college graduate), health insurance provider (private vs. public), and prepregnancy body mass index (BMI; < 25 kg/m^{2}, 25 to < 30 kg/m^{2}, ≥ 30 kg/m^{2}). These covariates were included in each model for the present analysis; categorical variables were modeled with the following reference levels: race/ethnicity (white), education level (high school), health insurance provider (private), and prepregnancy BMI (< 25 kg/m^{2}). In addition, we performed a sensitivity analysis treating the most frequently occurring categorical variables (white, college graduate, private health insurance, prepregnancy BMI < 25 kg/m^{2}) as reference levels for comparison.
Urinary Exposure and Mediator Biomarkers
NSF International (Ann Arbor, MI) measured nine urinary phthalate metabolites by high performance liquid chromatography and tandem mass spectrometry as described previously (Ferguson et al. 2014; Lewis et al. 2013). All metabolites were detected in > 95% of urine samples, and levels below the limit of detection were kept as is if reported and otherwise replaced by the limit of detection divided by the square root of 2. In addition to individual metabolites, a summed measure of di(2ethylhexyl) phthalate (DEHP) metabolites, including mono(2ethylhexyl) phthalate (MEHP), mono(2ethyl5hydroxyhexyl) phthalate (MEHHP), mono(2ethyl5oxohexyl) phthalate (MEOHP), and mono(2ethyl5carboxypentyl) phthalate (MECPP), was created based on nanomolar concentrations. Because DEHP metabolites are highly correlated, and the strongest associations observed with all preterm birth were for MEHP and MECPP, MEHHP and MEOHP were not examined separately in this analysis. Cayman Chemical (Ann Arbor, MI) measured total 8isoprostane in affinity purified urine samples via enzyme immunoassay, and detection was to 3.9 pg/mL. Levels below the limit of detection (4%) for 8isoprostane concentrations were treated the same as phthalate metabolites. In addition, to adjust for urinary dilution, specific gravity was measured in all samples at the time of phthalate analysis with a digital handheld refractometer (Atago Co., Ltd., Tokyo, Japan). We did not correct phthalate metabolite or 8isoprostane concentrations to specific gravity measures, but instead models were adjusted for specific gravity to achieve more precise estimates of urine biomarkers. Urinary phthalate metabolites and 8isoprostane measurements are logarithmtransformed and then standardized with mean of 0 and standard deviation of 1 throughout the mediation analysis for better interpretability and comparability across different phthalates.
Of the overall 482 subjects, 357 had measures of phthalate metabolites and 8isoprostane available from all three study visits. One hundred twentyfive subjects had at least one missing 8isoprostane measurement (n = 8, 61, and 73 for visits 1–3, respectively) and 118 had at least one missing urinary phthalate metabolite measurement (n = 3, 60, 70 for visits 1–3, respectively). For time points where 8isoprostane or phthalate metabolites were missing, we applied single imputation of the geometric average of that subject’s urinary concentrations from other visits. We additionally performed sensitivity analyses to compare results across three methods of treating missing data: a) subjects with all exposure mediator measures available or a complete case analysis; b) data with imputation based on subjectspecific average; and c) data with imputation based on average of all subjects at that visit.
Mediation Methods
Traditional mediation analysis enables the researcher to identify the proportion of a statistical relationship between exposure A and outcome Y that occurs through a change in the mediator M using a sum or product coefficient method. The counterfactual approach to mediation analysis differs from the traditional approaches developed by Baron and Kenny (1986) in that it clarifies the assumptions that allow for causal interpretation of results and enables extension to more sophisticated and general data structures and models (Pearl 2014). The underlying assumptions are key to application of the counterfactual framework. Most important, this includes the assumption of no unmeasured confounding within a) the exposure–outcome relationship; b) the mediator–outcome relationship; and c) the exposure–mediator relationship. Furthermore, the analysis assumes that none of the mediator–outcome confounders are affected by exposure. Our ability to meet these assumptions in the present example is addressed in the Discussion section.
The statistical notation for the models is as follows: Let A(t) represent the logarithm of urinary phthalate metabolites at visit t, A = [A(1), A(2), A(3)], and A̅ represent the average of the logtransformed phthalate metabolite concentrations across the three visits; M(t) represent the logarithm of the mediator 8isoprostane at visit t, M = [M(1), M(2), M(3)], and M̅ represent the average of the logtransformed concentrations of 8isoprostane across the three visits; Y represent the outcome, all or spontaneous preterm birth; V represent the set of timeinvariant baseline covariates listed above; and L(t) represent the only timevarying covariate, urinary specific gravity, at visit t, L = [L(1), L(2), L(3)], and L̅ represent the average across the three visits. A̅ and M̅ have been standardized to have a mean of zero and a standard deviation of one for ease of interpretability. Finally, uppercase letters represent the random variables and the corresponding lowercase letters represent possible realizations of the random variables. All statistical analyses were performed using R (version 3.1.0; R Project for Statistical Computing). A more detailed and complete description of the subsequent statistical models, methods, and corresponding R code are available in Section S1. Steps for Mediation Analysis with Longitudinal Exposure and Mediator with a Binary Outcome.
Method 1: utilizing exposure and mediator averages. We first applied a simple counterfactual approach toward performing a mediation analysis on these data by creating subjectspecific averages of the timevarying exposure, mediator, and timevarying covariate variables. The two statistical models used are as follows:
Model 1: logit[P(Y = 1a, m, v, l)] = β_{y0} + β_{ya}a̅ + β_{ym}m̅ + βT
yvv + β_{yl}l̅
Model 2: E[M̅a, v, l] = β_{m0} + β_{ma}a̅ + βT
mvv + β_{ml}l̅
In model 2 we applied inverse probability weightings to account for the nested case–control study design. The estimates of NDE and NIE are based on the expressions in the appendix of VanderWeele and Vansteelandt (2010), and the standard errors are obtained by bootstrap. The total effect can be expressed as the product of NDE and NIE in the original beta coefficient scale for linear models. On the difference scale for log odds ratios, the total effect can be defined as the sum of NDE and NIE and the proportion mediated is then .
Method 2: utilizing exposure and mediator averages with interaction. Our second approach was to employ the same counterfactual framework, but allow for an interaction between the exposure and mediator averages. The NDE and NIE can be estimated in the same fashion as the previous method from the following two models:
Model 3: logit[P(Y = 1a, m, v, l)] = β_{y}_{0} + β_{ya}a̅ + β_{ym}m̅ + β_{yi}a̅m̅ + βT
yvv + β_{yl}l̅
Model 4: E[M̅a, v, l] = β_{m}_{0} + β_{ma}a̅ + βT
mvv + β_{ml}l̅
The interaction between urinary phthalate metabolite and urinary 8isoprostane is taken into account in model 3. As in the previous method, for model 4 we applied inverse probability weightings to adjust for the case–control study design. Because the estimate of NDE depends on the level of the covariates, we fixed the continuous covariates at their arithmetic means and set the categorical covariates at their reference levels for this and subsequent analysis.
Method 3: longitudinal approach utilizing repeated measures of exposure and mediator. Our third approach was to apply the randomized interventional analogue of natural direct and indirect effects defined by VanderWeele and Tchetgen Tchetgen (2014) in the longitudinal setting. Using longitudinal detail entails additional assumptions. These include a) that no timevarying confounder is affected by prior exposure or mediator; and b) that at each time point the exposure–mediator, exposure–outcome, and mediator–outcome relationships conditional on the past are unconfounded. Under the latter assumption the randomized interventional analogues of natural direct and indirect effects can be derived from the gformula (Robins 1986). When there are no timevarying confounders, this gformula reduces to the longitudinal version of Pearl’s mediation formula (2014). There are no timevarying confounders in our scenario, because specific gravity is not a confounder of the exposure–outcome relationship. Thus, we applied the latter formula to estimate the natural direct and indirect effects using marginal structural models (Robins et al. 2000). The marginal structural models are as follows:
Model 5: logit[P(Y_{am} = 1 / a, m, v, l)] = β_{y}_{0} + β_{ya}cum(a) + β_{ym}cum(m) + βT
yvv + β_{yl}l̅
Model 6: E[M_{a}(t)a(1), …, a(t), v, l(t)] = β_{m}_{0}(t) + β_{ma}(t)avg[a(t)] + βT
mv(t)v + β_{ml}(t)l(t) for t = 1, 2, 3
where cum is the cumulative total of the measures across all time points, and avg[a(t)] is the average of a (exposure), up to and including time point t. The conditional model can be interpreted as a structural model if there is no timedependent confounding. Inverse probability weighting to adjust for the casecontrol study design are applied to each of the three regressions in model 6. The direct and indirect effects are statistically defined in a similar fashion to the second approach, except that change from one exposure history to another exposure history in timevarying sense rather than change from one level of exposure to another level now needs to be specified. The expressions for direct and indirect effects are detailed in Section S1. Steps for Mediation Analysis with Longitudinal Exposure and Mediator with a Binary Outcome; Method 3 – Counterfactual Approach with Longitudinal Exposure and Longitudinal Mediator without Interaction (VanderWeele and Tchetgen Tchetgen 2014). The scenario constructed for interpreting the mediation effect here corresponds to an lnunit increase in exposure across all three time points. The NDE, NIE, and proportion mediated are then expressed accordingly and the corresponding standard errors are obtained by bootstrap for inference purposes.
Method 4: longitudinal approach utilizing repeated measures of exposure and mediator with interaction. Our fourth approach was to extend the third approach to account for an exposure–mediator interaction, in the same way that Method 2 extends Method 1. More detailed descriptions are provided in Section S1. Steps for Mediation Analysis with Longitudinal Exposure and Mediator with a Binary Outcome; Method 4 – Counterfactual Approach with Longitudinal Exposure and Longitudinal Mediator with Interaction (VanderWeele and Tchetgen Tchetgen 2014). Similarly, we applied a gformula approach to estimate the direct and indirect effects using marginal structural models and inverse probability of treatment weighting (Robins et al. 2000) under the same set of assumptions. The marginal structural models are:
Model 7: logit[P(Y_{am} = 1 a, m, v, l)] = β_{y}_{0} + β_{ya}cum(a) + β_{ym}cum(m) + β_{yi}cum(x)cum(m) + βT
yvv + β_{yl}l̅
Model 8:E[M_{a}(t)a(1), …, a(t), v, l(t)] = β_{m}_{0}(t) + β_{ma}(t)avg[a(t)] + βT
mv(t)v + β_{ml}(t)l(t) for t = 1, 2, 3
Error Terms in Model 8 ϵ = [ϵ(1), ϵ(2), ϵ(3)]^{T} ~ MVN(0, σ^{2}∑),
where Σ can be an arbitrary 3 × 3 positivedefinite matrix with 1s along the diagonal. Model 8 is the same as Model 6 except for an additional assumption, that the three errors terms from the three regressions follow a zeromean multivariate normal distribution with common variance and unstructured correlation. The conditional model can be interpreted as a structural model if there is no timedependent confounding. The timevarying coefficients are jointly estimated using generalized least squares. The reason for assuming a joint correlation structure across the three regression residuals is to obtain the estimated variance of Σ3
t = 1 M(t), a necessary quantity for estimating NDE in the scenario with exposure–mediator interaction. The inverse probability weighting to adjust for the case–control study design is applied to model 8 as before. The expressions for direct and indirect effects and their derivations are provided in Section S1. Steps for Mediation Analysis with Longitudinal Exposure and Mediator with a Binary Outcome; Method 4 – Counterfactual Approach with Longitudinal Exposure and Longitudinal Mediator with Interaction (VanderWeele and Tchetgen Tchetgen 2014). We used the same scenario as constructed in Method 3 where we assessed the direct effect, indirect effect, and proportion mediated based on unit change in exposure level at each time point and subsequently drew inference based on the standard errors obtained by bootstrap.
Results
Beta coefficients and standard errors from regression models for each approach are similar to results observed previously (Tables S1–S4; Ferguson et al. 2014, 2015a, 2015b). Urinary phthalate metabolites were positively associated with preterm birth, and associations were stronger for spontaneous preterm birth specifically. Effect estimates were greatest in magnitude, and confidence intervals most precise, for associations with DEHP metabolites or MBP. Also, 8isoprostane was positively associated with preterm and particularly spontaneous preterm birth, and with urinary phthalate metabolites.
Mediation results for Methods 1–4 are shown in Tables 1–4, respectively. Each table shows results for all and spontaneous preterm births separately. Effect estimates represent lnodds ratios in association with a 1 standard deviation change (from mean – 1 to the mean) in exposure average (Methods 1–2) or across the three study visits (Methods 3–4). For all preterm birth, mediation analysis was performed for MEHP, MECPP, ΣDEHP, and mononbutyl phthalate (MBP) only, because associations in regression models were not statistically significant for other metabolites.
Table 1. Effect estimates^{a} (95% confidence intervals) with lnunit increase in exposure and estimated percent mediated, calculated from regression estimates and standard errors generated from models 1–2 (Table S1) under Method 1: Counterfactual approach utilizing exposure and mediator averages.
Table 2. Effect estimates (95% confidence intervals) with lnunit increase in exposure and estimated percent mediated calculated from regression estimates and standard errors generated from models 3–4 (Table S2) under Method 2: Counterfactual approach utilizing exposure and mediator averages with interaction.
Table 3. Effect estimates (95% confidence intervals) with lnunit increase in exposure and estimated percent mediated calculated from regression estimates and standard errors generated from models 5–6 (Table S3) under Method 3: Longitudinal approach utilizing repeated measures of exposure and mediator.
Table 4. Effect estimates (95% confidence intervals) with lnunit increase in exposure and estimated percent mediated calculated from regression estimates and standard errors generated from models 7–8 (Table S4) under Method 4: Longitudinal approach utilizing repeated measures of exposure and mediator with interaction.
We observed significant mediation (i.e., NIE with confidence intervals that did not include zero) for the relationship between urinary phthalate metabolites and all preterm birth as well as spontaneous preterm birth by 8isoprostane. Significant NIE were observed across all Methods, and we consistently estimated a greater proportion mediated for relationships with spontaneous versus all preterm births. In regard to specific phthalate metabolites, the greatest estimated proportion of the preterm birth relationship mediated by 8isoprostane was for MBP across all methods. For spontaneous preterm birth, the greatest estimated proportion mediated for each approach was either for MECPP or MBP.
Through applying the four different mediation approaches, we were able to compare a) the effect of incorporating an exposure–mediator interaction and b) the effect of using longitudinal exposure and mediator trajectories instead of reducing repeated measurements to an average. For comparison of results across mediation methods in the text, we present results for MECPP and MBP only because we observed the greatest mediation by 8isoprostane for models of these metabolites in relation to spontaneous preterm birth (see other results in Tables 1–4). Incorporating interaction into mediation models generally increased the estimated percent mediated. In methods with exposure and mediator averages (Methods 1 and 2), estimated percent mediated increased from 51% to 61% for MECPP and from 45% to 70% for MBP. In methods with exposure and mediator treated longitudinally (Methods 3 and 4), estimated percent mediated increased from 47% to 60% for MECPP and from 39% to 81% for MBP.
Using repeated measures of exposure and mediator variables did not have a clear impact on the estimated proportion mediated. However, it did affect the precision of NIE (i.e., mediated effect) as evidenced by narrowing of the confidence intervals for these estimates. When methods without interaction terms (Methods 1 and 3) were compared, the width of the NIE confidence intervals for spontaneous preterm birth decreased from 0.358 to 0.314 for MECPP and from 0.429 to 0.333 for MBP. When methods with interaction terms (Methods 2 and 4) were compared, the width of the NIE confidence intervals for spontaneous preterm birth decreased from 0.369 to 0.307 for MECPP and from 0.436 to 0.342 for MBP.
In sensitivity analyses, we examined the impact of using imputation by subjectspecific average for exposure and mediator values using the most general Method 4. (Method 4 is considered the most general modeling framework because it includes multiple measures of exposure and mediator and no exposure by mediator interaction. Methods 1–3 can all be viewed as special cases of Method 4, because they are all reduced versions of that method.) As alternatives to imputation with subjectspecific exposure or mediator averages, we estimated NIE and NDE when imputing by a) using only subjects with available measures (complete case analysis); or b) using the population averages of those measures. This comparison showed that the estimated percent mediated was attenuated in models with imputation by population averages and that our approach with single imputation gave point estimates that were most consistent with those observed in models using subjects without missing values (Table S5).
Additionally, in Methods 2 and 4 with interaction we examined the impact of changing reference categories for categorical covariates to those observed most frequently in the study population. The estimated proportion mediated was similar across these choices (data not shown).
Discussion
We previously observed that urinary phthalate metabolites measured at multiple time points during pregnancy were associated with increases in urinary 8isoprostane, a systemic biomarker of oxidative stress, and that both were associated with increased risk of all and especially spontaneous preterm birth. In the present mediation analysis, we statistically demonstrated within a causal framework that the relationship between phthalate exposure and spontaneous preterm birth is mediated in part by phthalateinduced oxidative stress, which has not been elucidated clearly in human or animal studies previously. The higher estimated proportion mediated by 8isoprostane observed in models of spontaneous preterm births alone indicates that the oxidative stress pathway may be particularly relevant for this subtype.
The latter finding has biologic plausibility based on what is known about the role of oxidative stress and prematurity. One of the leading explanations for pathways to spontaneous preterm birth is an increase in inflammation at the maternal–fetal interface, leading to early initiation of parturition pathways (Challis et al. 2009). Oxidative stress is tightly linked to inflammation, and could be the origin of inflammation in this pathway. Additionally, recent evidence suggests that oxidative damage to the membranes surrounding the fetus could result in preterm membrane rupture—a subtype of preterm birth that we included in our “spontaneous” category (Menon 2014). This mediation analysis provides a quantitative assertion of this plausible mechanism.
In addition to our substantive findings, we quantitatively compared mediation analyses using average versus repeated measures of the exposure and mediator, as well as the effect of allowing for an exposure–mediator interaction. Fully using repeated measures of exposure and mediator, rather than condensing them to averages, resulted in narrower confidence intervals for indirect effect estimates for models both with and without interaction. This clearly illustrates the improved precision of these advanced mediation models. Furthermore, longitudinal analysis is a more appropriate approach when there is temporal variation in exposure, autocorrelation in the repeated measures, and potential timevarying confounding. Thus, to maintain power and address these issues, longitudinal models should be considered by researchers who have repeated measures of exposure and mediator variables available.
Allowing for an interaction term between exposure and mediator increased the observed percent mediated for both average and repeated measures models. This finding suggests that a biologic interaction may exist in these relationships. Specifically, elevated exposure to phthalates in combination with higher levels of oxidative stress may result in a greater than cumulative increase in risk of spontaneous preterm birth. This has important connotations for future research in environmental impacts on pregnancy, because many contaminants found in pregnant mothers have demonstrated capacity to cause oxidative stress. Examining exposure to mixtures of oxidative stress inducing compounds in relation to preterm birth will be an important next step in this line of research. Notably, although we observed the largest estimated proportion mediated in the model with interaction, this may not be the case for all analyses. Other investigators using a causal framework for mediation analysis should test the inclusion of an exposure–mediator interaction using the models presented here, and should keep these in the model only if they have substantial impact (i.e., show a large change in the estimated proportion mediated), have biologic plausibility, and meet the assumptions described above.
Despite the fact that we noticed a larger estimated proportion of mediated effects and narrower confidence intervals under Method 4, we do not know the truth in this or any given data set and thus cannot claim Method 4 is the best. However, because Method 4 is the most general framework that most fully utilizes the longitudinal information in the data and allows for potential exposure by mediator interaction, an investigator can start with this approach. The advantage of starting with a more general model is that if the interaction has any impact, this approach will capture that situation, and in absence of interaction, the interaction coefficient will be close to null and can be removed.
There are several limitations to our study and the present mediation analysis. In regard to our study design, we measured circulating biomarkers of phthalate exposure and oxidative stress, despite the fact that measures at the maternal–fetal interface may be more biologically relevant. However, markers in urine may be indicative at least in part of activity in the uterine compartment, and collection of urine is much more feasible (i.e., less invasive) than tissue or fluid samples from the uterus during pregnancy. A second limitation is that 8isoprostane is not a direct measure of reactive oxygen species production but only a proxy. Although we attempted to examine the phthalate–preterm birth mediation by oxidative stress, in reality we only examined the mediation of the relationship by 8isoprostane, or the mediation by oxidative stress detected by 8isoprostane in urine. Other factors may influence the relationship between phthalate exposure and 8isoprostane levels, making 8isoprostane incompletely representative of the oxidative stress that phthalates produce. Thus the mediation by total oxidative stress may be underestimated in this analysis.
In addition to these study design limitations, the mediation analysis within the counterfactual framework makes strong assumptions for “no unmeasured confounding” of the outcome–exposure, outcome–mediator, and mediator–exposure relationships as described in the Methods section. Although we were able to examine a large number of potential confounders in this analysis, there is always the possibility of unmeasured confounding. One example could include changes in diet, which have been linked to phthalate exposure biomarkers, oxidative stress, and prematurity. However, associations between this potential confounder and preterm birth are questionable, and 8isoprostane is known for being less sensitive to diet compared to other oxidative stress biomarkers (Milne et al. 2005). Sensitivity analyses can in some instances be applied to theoretically observe how violations of these “no unmeasured confounding” assumptions would affect the mediation results and are described in detail in the literature (VanderWeele and Arah 2011) but are not applied here.
Finally, a major assumption in this analysis is that phthalates cause oxidative stress (i.e., temporal ordering). If oxidative stress causes an increase in urinary phthalate metabolite excretion, or if the connection between the two is not causal but instead a result of unmeasured confounding, then the interpretation of these results would be different. However, based on animal and cellular studies, there is moderate evidence that phthalate exposure during pregnancy causes an increase in reactive oxygen species release which can be measured by urinary excretion of 8isoprostane (Rusyn et al. 2001; Tetz et al. 2013).
Despite these limitations and assumptions our study has many strengths, namely the ability to examine these associations in a casecontrol population with a large number of subjects and repeatedly measured biomarkers of both exposure and mediator. It is also the first analysis, to our knowledge, that attempts to identify through an epidemiologic study the mediators of relationships between an environmental contaminant exposure and preterm birth. A number of studies examine relationships between environmental chemicals and prematurity, but are limited by the inability to establish causality. Although mediation analysis still does not concretely establish a causal pathway, it provides an additional step that none of these previous studies have been able to take. These methods may be particularly useful for future studies in this realm, but also for the vast number of epidemiologic studies attempting to identify causal pathways using molecular biomarkers of intermediate effects.
In conclusion, our methodological approach of gradually building more complex mediation models incorporating exposure–mediator interactions and repeated measures of the exposure and mediator illustrate the utility of the longitudinal study design and application of these more sophisticated analytical approaches. Additionally, conditional on the underlying assumptions, these findings provide causal evidence for mediation of the previously observed associations between phthalate exposure and preterm birth by oxidative stress.
References
Adibi JJ, Hauser R, Williams PL, Whyatt RM, Calafat AM, Nelson H, et al. 2009. Maternal urinary metabolites of di(2ethylhexyl) phthalate in relation to the timing of labor in a US multicenter pregnancy cohort study. Am J Epidemiol 169:1015–1024.
ATSDR (Agency for Toxic Substances and Disease Registry). 2001. Toxicological Profile for Dinbutyl Phthalate (DBP). http://www.atsdr.cdc.gov/substances/toxsubstance.asp?toxid=167 [accessed 9 September 2015].
ATSDR. 2002. Toxicological Profile for Di(2ethylhexyl) phthalate (DEHP). http://www.atsdr.cdc.gov/substances/toxsubstance.asp?toxid=65 [accessed 9 September 2015].
Baron RM, Kenny DA. 1986. The moderator–mediator variable distinction in social psychological research: conceptual, strategic, and statistical considerations. J Pers Soc Psychol 51:1173–1182.
Bind MA, Lepeule J, Zanobetti A, Gasparrini A, Baccarelli A, Coull BA, et al. 2014. Air pollution and genespecific methylation in the Normative Aging Study: association, effect modification, and mediation analysis. Epigenetics 9:448–458.
Challis JR, Lockwood CJ, Myatt L, Norman JE, Strauss JF III, Petraglia F. 2009. Inflammation and pregnancy. Reprod Sci 16:206–215.
Ferguson KK, McElrath TF, Chen YH, LochCaruso R, Mukherjee B, Meeker JD. 2015a. Repeated measures of urinary oxidative stress biomarkers during pregnancy and preterm birth. Am J Obstet Gynecol 212(208):208.e1–208.e8, doi: 10.1016/j.ajog.2014.08.007.
Ferguson KK, McElrath TF, Chen YH, Mukherjee B, Meeker JD. 2015b. Urinary phthalate metabolites and biomarkers of oxidative stress in pregnant women: a repeated measures analysis. Environ Health Perspect 123:210–216, doi: 10.1289/ehp.1307996.
Ferguson KK, McElrath TF, Meeker JD. 2014. Environmental phthalate exposure and preterm birth. JAMA Pediatr 168:61–67.
Lewis RC, Meeker JD, Peterson KE, Lee JM, Pace GG, Cantoral A, et al. 2013. Predictors of urinary bisphenol A and phthalate metabolite concentrations in Mexican children. Chemosphere 93:2390–2398.
MacKinnon DP. 2012. Introduction to Statistical Mediation Analysis. Hoboken, NJ:Taylor & Francis Group.
Meeker JD, Hu H, Cantonwine DE, LamadridFigueroa H, Calafat AM, Ettinger AS, et al. 2009. Urinary phthalate metabolites in relation to preterm birth in Mexico City. Environ Health Perspect 117:1587–1592, doi: 10.1289/ehp.0800522.
Menon R. 2014. Oxidative stress damage as a detrimental factor in preterm birth pathology. Inflammation 5:567, doi: 10.3389/fimmu.2014.00567.
Milne GL, Musiek ES, Morrow JD. 2005. F_{2}isoprostanes as markers of oxidative stress in vivo: an overview. Biomarkers 10(suppl 1):S10–S23.
Pearl J. 2014. Interpretation and identification of causal mediation. Psychol Methods 19:459–481.
Robins J. 1986. A new approach to causal inference in mortality studies with a sustained exposure period—application to control of the healthy worker survivor effect. Mathematical Modelling 7:1393–1512.
Robins JM, Hernán MA, Brumback B. 2000. Marginal structural models and causal inference in epidemiology. Epidemiology 11:550–560.
Rubin DB. 1974. Estimating causal effects of treatments in randomized and nonrandomized studies. J Educ Psychol 66:688–701.
Rusyn I, Kadiiska MB, Dikalova A, Kono H, Yin M, Tsuchiya K, et al. 2001. Phthalates rapidly increase production of reactive oxygen species in vivo: role of Kupffer cells. Mol Pharmacol 59:744–750.
Tetz LM, Cheng AA, Korte CS, Giese RW, Wang P, Harris C, et al. 2013. Mono2ethylhexyl phthalate induces oxidative stress responses in human placental cells in vitro. Toxicol Appl Pharmacol 268:47–54.
Valeri L, VanderWeele TJ. 2013. Mediation analysis allowing for exposure–mediator interactions and causal interpretation: theoretical assumptions and implementation with SAS and SPSS macros. Psychol Methods 18:137–150.
VanderWeele TJ, Arah OA. 2011. Bias formulas for sensitivity analysis of unmeasured confounding for general outcomes, treatments, and confounders. Epidemiology 22:42–52.
VanderWeele TJ, Tchetgen Tchetgen E. 2014. Mediation Analysis with TimeVarying Exposures and Mediators. Harvard University Biostatistics Working Paper Series. Working Paper 168. http://biostats.bepress.com/harvardbiostat/paper168/ [accessed 21 January 2016].
VanderWeele TJ, Vansteelandt S. 2010. Odds ratios for mediation analysis for a dichotomous outcome. Am J Epidemiol 172:1339–1348.
Whyatt RM, Adibi JJ, Calafat AM, Camann DE, Rauh V, Bhat HK, et al. 2009. Prenatal di(2ethylhexyl)phthalate exposure and length of gestation among an innercity cohort. Pediatrics 124:e1213–e1220.
Wolff MS, Engel SM, Berkowitz GS, Ye X, Silva MJ, Zhu C, et al. 2008. Prenatal phenol and phthalate exposures and birth outcomes. Environ Health Perspect 116:1092–1097, doi: 10.1289/ehp.11007.
Announcements
Come visit EHP EditorinChief Sally Darney at the 29th Annual Scientific Conference of the International Society of Environmental Epidemiology, to be held in Sydney, Australia, 24–28 September 2017. She looks forward to talking with EHP authors, reviewers, and editors. Contact Sally directly to make an appointment with her.
EHP is pleased to announce that “Maternal and Cord Blood Manganese Concentrations and Early Childhood Neurodevelopment among Residents near a MiningImpacted Superfund Site,” published in EHP on 28 June 2017, has been selected by the Children’s Environmental Health Network (CEHN) as its August 2017 Article of the Month. (more…)