Urinary phthalate metabolites and biomarkers of oxidative stress in pregnant women: a repeated measures analysis.
Phthalate diesters are used as plasticizers and solvents in a variety of consumer products, and can readily enter human systems through ingestion, inhalation, and dermal absorption (Agency for Toxic Substances and Disease Registry 2001, 2002). Although diesters are metabolized and excreted quickly, constant contact results in daily exposures for most of the U.S. population. Metabolites are consistently detected in urine of pregnant women in populations worldwide (Cantonwine et al. 2013; Lin et al. 2011; Woodruff et al. 2011; Zeman et al. 2013).
Although phthalates are best known for their action as endocrine disruptors, there is also evidence from in vitro and animal studies that mono(2-ethylhexyl) phthalate (MEHP) may cause oxidative stress by inducing release of reactive oxygen species (ROS) and/or impairing antioxidant defenses (Erkekoglu et al. 2010; Kasahara et al. 2002; Tetz et al. 2013; Zhao et al. 2012). However, few studies have examined this association in humans. Three cross-sectional studies have observed associations between some phthalate metabolites and serum levels of bilirubin, a potent antioxidant, and systemic markers of oxidative stress including serum gammaglutamyl transferase, and urinary malondialdehyde (MDA) and 8-hydroxydeoxyguanosine (8-OHdG) (Ferguson et al. 2011, 2012; Hong et al. 2009). A more recent study in elderly subjects also observed an association between summed di(2-ethylhexyl) phthalate (DEHP) metabolites and MDA (Kim et al. 2013).
To our knowledge, no studies have examined the relationship between urinary phthalate metabolites and biomarkers of oxidative stress during pregnancy, when these possible effects represent particular concern given the gestational vulnerability of the developing fetus. Increases in oxidative stress biomarkers in pregnant women have been associated with pregnancy loss, preeclampsia, preterm birth, and fetal growth restriction (Agarwal et al. 2012; Stein et al. 2008). Additionally, the potential effects of phthalate exposures on oxidative stress are relevant to a number of other outcomes in the general population, such as infertility, various cancers, and type 2 diabetes. In the present study, we examined associations between repeated measures of urinary phthalate metabolites and biomarkers of oxidative stress in pregnant women.
Study population. Pregnant women were recruited before 15 weeks gestation at Brigham and Women's Hospital in Boston, Massachusetts, from 2006 through 2008 as part of a large prospective cohort study and provided informed consent upon enrollment. Participants were included if they had a singleton pregnancy that resulted in a live birth. Women were followed throughout the duration of pregnancy and provided demographic and anthropometric data, urine samples from up to four study visits (targeted for 10, 18, 26, and 35 weeks gestation), as well as birth outcome data at delivery. Gestational age was calculated from last menstrual period and validated with first-trimester ultrasound; if gestational ages calculated by the two methods differed by > 8%, ultrasound dating was used. For the present study, 130 women who delivered preterm and 352 random controls were selected, and their urine samples were extracted from -80[degrees]C storage for laboratory analysis (n = 482 subjects total). Institutional review board approval for this study was obtained from the University of Michigan and Brigham and Women's Hospital.
The nested case-control study was designed with the intention of examining associations between urinary phthalate metabolites across gestation and preterm birth (Ferguson et al. 2014b). The present analysis examining the relationship between urinary phthalate metabolites and biomarkers of oxidative stress was a secondary aim of this study to help inform potential biological mechanisms involved. Our goal was to characterize these associations in a population that would be generalizable to the overall cohort. Therefore, the analyses were weighted using inverse probability weightings calculated based on selection probabilities for cases (90.1%) and controls (33.9%) from the parent cohort population (Jiang et al. 2006).
Urinary phthalate metabolites. All available urine samples (n = 1,693) were assayed for concentrations of nine phthalate metabolites using high-performance liquid chromatography and tandem mass spectrometry by NSF International in Ann Arbor, Michigan (Lewis et al. 2013; Silva et al. 2007). All metabolites were highly detectable (> 95%); levels below the limit of detection (LOD) were replaced with the LOD divided by the square root of 2 (Hornung and Reed 1990). All distributions were right-skewed and natural log-transformed to meet normality assumptions for statistical analysis. To adjust for urine dilution, specific gravity (SG) was measured at the time of phthalate analysis using a handheld refractrometer (Atago Co., Ltd., Tokyo, Japan). For bivariate analysis we created SG-corrected concentrations using the following formula: [P.sub.c] = P[(1.015 - 1)/(SG - 1)], where [P.sub.c] is the corrected phthalate concentration, P is the raw concentration, 1.015 is the median SG for the study population, and SG is the specific gravity for the sample (Meeker et al. 2009). In regression models uncorrected phthalate metabolite levels were used and models were adjusted for SG as a covariate.
In addition to examining individual urinary phthalate metabolites, we created a summed measure of DEHP metabolites ([summation]DEHP). MEHP, mono(2-ethyl-5-hydroxyhexyl) phthalate (MEHHP), mono(2ethyl-5-oxohexyl) phthalate (MEOHP), and mono(2-ethyl-5-carboxypentyl) phthalate (MECPP) concentrations were converted from micrograms per liter to micromoles per liter using molecular weights (278, 294, 292, and 308 g/mol) and summed to create [summation]DEHP.
Oxidative stress biomarkers. All urine samples with sufficient volume remaining following phthalate measurement (n = 1,678) were assayed for levels of 8-OHdG and total 8-isoprostane using enzyme immunoassay by Cayman Chemical (Ann Arbor, MI). For 8-isoprostane measurement only, samples were first hydrolyzed and passed through columns for affinity purification. LODs were 3.9 pg/mL and 10.3 pg/mL for urine 8-isoprostane and 8-OHdG concentrations, respectively. As with phthalate metabolites, undetected oxidative stress measures were replaced with the LOD divided by the square root of 2 (Hornung and Reed 1990). For calculation of distributions overall and by categorical covariates biomarkers were corrected for urinary SG using the formula applied to phthalate measures above. Uncorrected biomarker concentrations were used for multivariate models with SG as a covariate.
Statistical analysis. All statistical analysis was performed in R version 2.15.2 (R Foundation for Statistical Computing, Vienna, Austria) and SAS version 9.2 (SAS Institute Inc., Cary, NC). Unless stated otherwise, all analysis was performed with inverse probability weightings. Overall distributions of urinary phthalate metabolites and oxidative stress biomarkers were assessed using SG-corrected concentrations. Variability in oxidative stress biomarker concentrations across gestation was examined using intraclass correlation coefficients (ICCs), which represent a ratio of within to between individual variability (Rosner 2011). ICCs for urinary phthalate metabolites have been examined previously (unpublished data), and were similar to other studies observing greater reproducibility of mono-benzyl phthalate [MBzP; ICC = 0.61; 95% confidence interval (CI) = 0.56-0.65], mono-n-butyl phthalate (MBP; ICC = 0.57; 95% CI: 0.53, 0.62), and monoisobutyl phthalate (MiBP; ICC = 0.52, 95% CI: 0.48, 0.57) compared with DEHP metabolites (ICC range, 0.19-0.31) across pregnancy (Adibi et al. 2008; Braun et al. 2012). Geometric means and standard deviations of corrected concentrations were created by categorical covariate groups and compared using linear mixed models (LMMs) with subject-specific random intercepts to adjust for intra-individual correlation of measurements at multiple time points (nlme package in R) (Pinheiro et al. 2013). Covariates examined included race/ ethnicity (white, African American, or other), education level (high school, technical school, some college or junior college, college graduate or above), health insurance provider (public or private), body mass index (BMI) at the initial visit (continuous), tobacco and alcohol use during pregnancy (yes or no), parity (nulliparous or parous), sex of fetus (male or female), and use of assisted reproductive technology (yes or no). We also examined time-varying covariates, including BMI (continuous) and time of day of urine sample collection (dichotomized into before versus after 1300 hours based on histograms which displayed a nadir in urinary phthalate metabolite concentrations at that time of day).
Associations (fixed effects) between uncorrected urinary phthalate metabolites and oxidative stress biomarkers were estimated using LMM with subject-specific random intercepts to adjust for intraindividual correlation (random slopes did not improve model fit based on the Akaike information criterion). For all statistical models very concentrated urine samples (SG > 1.04) were excluded because biomarkers measured in those samples may be inaccurate (n = 4) (Boeniger et al. 1993; Braun et al. 2012). One oxidative stress biomarker was regressed on one phthalate metabolite per model. Crude models were adjusted for gestational age and urinary SG. Full models additionally included covariates that were significantly (p < 0.05) associated with one or both oxidative stress biomarkers as well as one or more urinary phthalate metabolites. Final full models were adjusted for urinary specific gravity, gestational age at sample collection, race/ethnicity, education level, health insurance provider, BMI, time of day of urine sample collection, and parity of infant. Tobacco and alcohol use were excluded from adjusted models because of the small number of subjects who used either during pregnancy (n = 31 and 20, respectively). Subjects with missing visit 1 covariates were excluded from LMM models; if time-varying covariates were missing, then individual time points only (not all data for that subject) were excluded. Statistical significance of effect estimates was assessed with an alpha level of 0.05.
Several sensitivity analyses were performed. First we examined associations in a stratified analysis of cases and controls separately. Second, we examined associations after excluding mothers who used alcohol and tobacco during pregnancy, because the number of users was too small to include these as covariates. Third, we created generalized additive mixed models (GAMMs) to investigate the possibility that the relationships between oxidative stress biomarkers and urinary phthalate metabolites were nonlinear. Fourth, to examine whether the relationship between phthalate exposure and oxidative stress differed based on time point in pregnancy, we examined interaction terms between urinary metabolites and either study visit or gestational age at sample collection, also in LMMs with random intercepts. These models were created in preterm cases and controls separately, because of the difference in proportions of cases compared with controls with measurements available at each time point, particularly at visit 4. Finally, we estimated associations based on LMM models adjusted for multiple urinary phthalate metabolites. We selected metabolites for multiple metabolite models based on their correlations with one another as well as their individual associations with the oxidative stress biomarkers.
Population demographics have been presented elsewhere for the case-control study population (Ferguson et al. 2014b) and are presented in Table 1 for the weighted population examined in the present analysis. Demographics for the population with complete covariates for regression models were almost identical and are presented in Supplemental Material, Table S1. Maternal age at visit 1 was 32.0 years on average, and most women were white (59%), well educated (71% with junior college, some college, or above), did not use tobacco (94%) or alcohol (95%) during pregnancy (Table 1), and did not use assisted reproductive technology to get pregnant (91%). More than half of infants were female (55%). Approximately half of the participants had normal to low BMI at visit 1 (53%) (Table 1), and as expected this proportion decreased across pregnancy (visit 2 = 47%, visit 3 = 28%, visit 4 = 15%). After weighting, deliveries in the study population were approximately 12% preterm (< 37 weeks gestation).
Urine samples were collected from study participants at up to four visits (mean 3.52 samples per subject; n = 3 subjects with one urine sample; n = 38 subjects with two samples; n = 148 subjects with three samples; n = 293 subjects with four samples). Samples for visit 1 (median 9.95 weeks gestation) were collected for 474 subjects (range, 4.17-19.1 weeks); for visit 2 (median, 18.0 weeks) samples were collected for 421 subjects (range, 14.9-32.1 weeks); for visit 3 (median, 26.0 weeks) samples were collected for 409 subjects (range, 22.9-36.3 weeks); and for visit 4 (median, 35.2 weeks) samples were collected for 374 subjects (range, 33.1-38.3 weeks). The proportion of cases with samples available was consistent for visits 1-3 (86-100%) but low for visit 4 (51%) because many had already delivered by that time point. Most urine samples were collected before 1300 hours at visits 1 (60%), 2 (66%), 3 (70%), and 4 (69%).
Urinary phthalate metabolite and oxidative stress biomarker distributions are presented in Table 2. Phthalate metabolites were detected in 95-100% of all samples measured, 8-OHdG was detected in all samples, and 8-isoprostane was below the LOD in 67 (4.0%) samples. As reported previously, correlations between phthalate metabolites were strongest within DEHP metabolites, as expected (Spearman r = 0.68-0.91), were moderate between MBzP, MBP, and MiBP (r = 0.46-0.62) and between mono(3-carboxypropyl) phthalate (MCPP) and DEHP metabolites (r = 0.35-0.46), and were weak between all other metabolites (r = 0.01-0.21) (Ferguson et al. 2014a). Correlations between the two oxidative stress markers at each visit were weak but statistically significant (r = 0.10-0.20, p < 0.05). 8-OHdG concentrations were more variable across pregnancy (ICC = 0.32; 95% CI: 0.27, 0.38) compared with 8-isoprostane (ICC = 0.60; 95% CI: 0.56, 0.64).
When oxidative stress biomarker concentrations were compared by categorical covariates, we observed different patterns for each marker (Table 3). Levels of 8-isoprostane were lowest in mothers who were white, had higher levels of education and private health insurance, were of lower BMI at visit 1, did not use tobacco or alcohol, and were parous. Few differences were observed by categorical covariates for 8-OHdG, although mothers with private health insurance had significantly lower levels. Associations with covariates were also different for urinary phthalate metabolites. MEHP concentrations (Table 3) and other DEHP metabolites (see Supplemental Material, Table S2) were higher in African-American compared with white mothers but no other differences were observed. Patterns for MBzP, MiBP, and monoethyl phthalate (MEP) (see Supplemental Material, Table S2) were similar to those for MBP (Table 3); higher concentrations were observed in mothers who were African American or other race/ethnicity compared with white, in mothers with lower education levels, in mothers with public compared with private health insurance, and in mothers with higher BMI at visit 1. No statistically significant differences in oxidative stress biomarkers and few differences in urinary phthalate metabolites were observed by fetus sex (lower MiBP concentrations in mothers of female vs. male fetus) or use of assisted reproductive technology (higher MEHHP concentrations and lower MBzP and MiBP concentrations in mothers who used assisted reproductive technology compared with those who did not) (data not shown).
For time-varying covariates, significantly higher oxidative biomarkers concentrations were observed in samples collected before versus after 1300 hours, and significantly lower urinary phthalate metabolite concentrations were observed in samples collected before versus after 1300 hours for all metabolites except MiBP and MEP, which were slightly higher in the morning (data not shown). When BMI was examined as a time-varying covariate, both oxidative stress biomarkers were positively associated with increasing BMI category (data not shown). MEHP was inversely associated with the highest BMI category, and MBzP was positively associated with the highest BMI category, but otherwise associations with urinary phthalate metabolites were close to the null (data not shown). Because BMI measures at each study visit may more accurately capture confounding by this variable, time-varying BMI was included as a covariate in fully adjusted models.
Effect estimates from adjusted models were similar to those from crude models (data not shown); adjusted results alone are presented in Table 4 (n = 464 subjects with complete data). Fixed-effect results are presented in the form of percent change in oxidative stress biomarker with an inter-quartile range (IQR) increase in untransformed phthalate metabolite. All phthalate metabolites were associated with higher 8-OHdG concentrations; the largest percent changes with an IQR increase in exposure observed were for MBzP (20.7%; 95% CI: 15.6, 26.1%), MBP (18.1%; 95% CI: 13.5, 22.9%), and MiBP (30.3%; 95% CI: 24.4, 36.5%). All metabolites were associated with significantly higher 8-isoprostane concentrations, and coefficients were larger compared with those estimated for 8-OHdG. However, as with 8-OHdG, the largest differences in 8-isoprostane were in association with MBzP (42.7%; 95% CI: 31.8, 54.4%), MBP (42.0%; 95% CI: 32.0, 52.7%), and MiBP (56.4%; 95% CI: 43.9, 69.9%).
Sensitivity analyses. LMM were replicated when stratifying by preterm birth case status. Effect estimates were similar with some exceptions (Table S3). In cases alone, associations between 8-OHdG and MBP (8.92%; 95% CI: 2.03, 16.3%), MiBP (20.3%; 95% CI: 10.2, 31.3%), and MEP (19.2%; 95% CI: 10.3, 28.8%) and also between 8-isoprostane and MBP (30.9%; 95% CI: 15.8, 48.0%) were smaller than corresponding associations estimated using weighted models. Effect estimates from models of controls alone were similar to those from weighted models. Effect estimates from weighted models will be closer to control rather than case estimates from Supplemental Material, Table S3, because there are many more controls than cases in the weighted population. Also, the effect estimates from the weighted models do not necessarily fall between case and control effect estimates because of differences in IQR ranges for cases compared with controls as well as differences in covariate distributions and association with exposure/outcome within case and control groups.
We also examined the effect of removing alcohol and tobacco users from the population; model estimates were similar to those for the overall population (data not shown). GAMM models were created to examine deviation of the oxidative stress-phthalate metabolite relationships from linearity. Models included the same covariates as in LMM. Based on visual inspection of smooth plots, relationships were linear for most metabolites and for those that were not, smooth curves deviated minimally from linearity (data not shown).
To identify potentially sensitive time points for the relationship between phthalate exposure and maternal oxidative stress, we additionally examined interaction terms in fully adjusted LMM between phthalate metabolite and gestational age at sample collection or visit of sample collection. No significant (p < 0.05) interactions with gestational age or study visit were observed for 8-OHdG or 8-isoprostane models (data not shown).
Finally, we examined the effect of including multiple phthalate metabolites in the same model. First we examined [summation]DEHP metabolites with MBP, because these measures were weakly correlated and would not create problems with multicollinearity (Chatterjee et al. 2000; Ferguson et al. 2014a), were strong predictors of oxidative stress in the present analysis, and previous animal and in vitro research on phthalate-induced oxidative stress has focused on the parent compounds of these two metabolites, DEHP (Erkekoglu et al. 2010; Kasahara et al. 2002; Rusyn et al. 2001) and dibutyl phthalate (DBP) (Kim et al. 2002; Shono and Taguchi 2014; Zhou et al. 2010). In the model for 8-OHdG, the effect estimate was lower for [summation]DEHP (4.32%; 95% CI: 0.98, 7.78%) and very similar for MBP (17.0%; 95% CI: 12.4, 21.9%). The same was true in models for 8-isoprostane ([summation]DEHP: 13.7%; 95% CI: 7.55, 20.1%; MBP: 37.3%; 95% CI: 27.6, 47.9%). Second, we examined the effect of including one predictor from each primary parent compound of interest in the same model ([summation]DEHP, MBzP, MBP, MEP, and MCPP). In models of 8-OHdG, effect estimates were diminished compared with single-exposure models (e.g., 13.3%; 95% CI: 7.84, 19.1% for MBzP), and coefficients for [summation]DEHP and MCPP lost statistical significance (see Supplemental Material, Table S4). In models of 8-isoprostane, effect estimates were also smaller in magnitude (e.g., 21.8%; 95% CI: 11.3, 33.2% for MBzP) but remained statistically significant for all phthalate metabolites (see Supplemental Material, Table S4).
We examined the association between repeated measures of urinary phthalate metabolites and 8-OHdG and 8-isoprostane as biomarkers of oxidative stress during pregnancy. We observed that all phthalate metabolites were associated with higher concentrations of both biomarkers. Associations were stronger with 8-isoprostane compared with 8-OHdG, and, among phthalates, MBzP, MBP, and MiBP showed strongest associations with both outcome measures.
Many different biomarkers have been used in environmental and other epidemiologic studies as a proxy of systemic levels of oxidative stress. These can have very different specificities, both in terms of mechanism (i.e., how they are produced) and downstream physiologic effect that the biomarkers themselves can have. The long-time goal of the National Institute of Environmental Health Sciences Biomarkers of Oxidative Stress Study has been to identify sensitive and specific markers of oxidative injury, as well as best methods for measuring these markers in animal and eventually human matrices (Kadiiska et al. 2013; National Institute of Environmental Health Sciences 2012). However, this task remains difficult, because of the numerous available markers and assays for detection, the long list of mechanisms that can cause oxidative stress, and the temporal instability of some biomarkers, among other reasons. We selected 8-OHdG and 8-isoprostane for measurement in this study because of their well-documented usefulness as systemic biomarkers of oxidative stress for establishing association with adverse health outcomes (Il'yasova et al. 2012), but also their representation of different cellular reactions to ROS exposure and the potential downstream effects of the biomarkers themselves.
Consistent with their low correlation in this and other studies (Stein et al. 2008), urinary levels of 8-OHdG and 8-isoprostane represent two distinct cellular processes. 8-OHdG is a DNA adduct formed in the presence of excess ROS (e.g., hydroxyl radicals) (Wu et al. 2004). Via repair mechanisms, oxidized nucleotides are excised from DNA and excreted in the urine, making repair capabilities an important factor in urinary concentrations (Wu et al. 2004). 8-isoprostane is formed in a non-enzymatic reaction between ROS and arachidonic acid and is advantageous because it is very specific to lipid oxidation, yet is not affected by dietary lipid intake, and is highly detectable in urine samples (Roberts and Morrow 2000). Measurement of specific isomers (e.g., the biologically active 8-iso-[PGF.sub.2[alpha]]) with liquid chromatography/ mass spectrometry is preferable to enzyme immunoassay, although multiple isomers can be representative of oxidative lipid damage, and the liquid chromatography/mass spectrometry method in a study of this size is cost prohibitive (Il'yasova et al. 2012; Smith et al. 2011).
In addition to indicating systemic oxidative stress, these different products may be markers of or play a direct role in physiologic changes that have adverse consequences for pregnancy. Oxidative DNA damage, specifically indicated by 8-OHdG, occurring in the intrauterine compartment could result in apoptosis at the maternal-fetal interface (Heazell et al. 2007), which can lead to poor vascularization of the placenta and consequently preeclampsia and/ or intrauterine growth restriction (Potdar et al. 2009). On the other hand, increased levels of prostaglandins such as 8-isoprostane may be particularly dangerous later in pregnancy because of their direct involvement in the preterm parturition pathway (Challis et al. 2009). This study illustrates that phthalates are associated with increases in both oxidative stress biomarkers in pregnant women, which suggests that exposure to phthalates could play a role in downstream pregnancy outcomes via multiple mechanisms.
In a number of in vitro studies, phthalates have been shown to cause increases in ROS and various markers of oxidative stress, potentially via activation of peroxisome proliferator-activated receptors or by increasing permeability of mitochondrial membranes (Hurst and Waxman 2003; Rosado-Berrios et al. 2011). These studies have been performed using DEHP and/or MEHP in a number of cell types, including placental cells (Tetz et al. 2013), Leydig cells (Erkekoglu et al. 2010; Zhou et al. 2013), neutrophils (Vetrano et al. 2010), and Kupffer cells (Rusyn et al. 2001). Other phthalates and their metabolites have been studied less frequently, although there is also evidence that they may be capable of inducing oxidative stress (O'Brien et al. 2001; Zhou et al. 2010).
In humans, studies of phthalates in relation to oxidative stress biomarkers have been limited. Two reports have examined the relationship among participants in the National Health and Nutrition Examination Survey, using gamma-glutamyl transferase and bilirubin in serum as markers of oxidative stress and a panel of phthalate metabolites similar to those measured in the present study. In the study of gamma-glutamyl transferase, positive associations were observed in association with MEHP only, although serum C-reactive protein, a systemic marker of inflammation that may also indicate oxidative stress, was positively associated with MiBP and MBzP (Ferguson et al. 2011). In the study of bilirubin, a potent antioxidant that may be inversely related to oxidative stress levels, DEHP and DBP metabolites as well as MCPP were found to be associated with significantly decreased bilirubin, although the strongest associations appeared to be for DEHP metabolites and MCPP (Ferguson et al. 2012).
Another cross-sectional study of urban-dwelling adults examined the relationship between phthalate metabolites and urinary MDA and 8-OHdG (n = 960) (Hong et al. 2009). The results from this analysis demonstrated significant and positive associations with DEHP metabolites as well as MBP and both oxidative stress biomarkers, although the positive associations with 8-OHdG lost significance in adjusted models. Regression coefficients were higher for MDA compared with 8-OHdG in DEHP metabolite models, but lower for MBP models. Regression coefficients were also higher for DEHP metabolites compared with MBP for both outcomes. Finally, one study in elderly subjects with measurements taken up to five times over 3 years observed a positive relationship between summed DEHP metabolites and urinary MDA levels (n = 560 subjects) (Kim et al. 2013). MBP levels were measured in the study but associations with MDA were not reported.
Our findings are somewhat consistent with these prior studies. As reported by Hong et al. (2009), we also observed positive associations between phthalates and 8-OHdG. Consistent with the studies of MDA, we observed higher levels of 8-isoprostane in association with urinary phthalate metabolites that appeared to be stronger than the associations observed with 8-OHdG. However, contrary to these studies, we observed the strongest associations for MBzP, MBP, and MiBP compared with DEHP metabolites for both outcome measures. This disparity could be attributable to differences in diet, product use, toxicant metabolism, and/or other differences between study populations. Considering our findings alone, the larger associations observed between MBzP, MBP, and MiBP with oxidative stress measures may be a result of lower temporal variability in those metabolites across pregnancy, which was observed in this study as well as others measuring phthalates at multiple time points during gestation (Adibi et al. 2008; Braun et al. 2012; Ferguson et al. 2014a). This would result in less measurement error and stronger associations, even if the true relationships between oxidative stress biomarkers and different phthalate metabolites were similar. Finally, it is possible that the associations observed are an effect of unknown confounders or other sources of error.
In conclusion, we report statistically significant increases in oxidative stress biomarkers in association with urinary phthalate metabolites during pregnancy. Our ability to detect these relationships may be attributable largely to our study design, with measurement of both urinary phthalate metabolites and oxidative stress biomarkers at up to four time points per subject across gestation. These associations with phthalate exposure may be important for pregnancy outcomes that are mediated by oxidative stress mechanisms. Additional exploration of these associations in other populations, particularly in nonpregnant women as well as men of reproductive age, children, and the elderly, may be of great importance for a range of other health outcomes that have been linked to phthalates in epidemiologic studies.
Adibi JJ, Whyatt RM, Williams PL, Calafat AM, Camann D, Herrick R, et al. 2008. Characterization of phthalate exposure among pregnant women assessed by repeat air and urine samples. Environ Health Perspect 116:467-473; doi:10.1289/ehp.10749.
Agarwal A, Aponte-Mellado A, Premkumar BJ, Shaman A, Gupta S. 2012. The effects of oxidative stress on female reproduction: a review. Reprod Biol Endocrinol 10:49; doi:10.1186/1477-7827-10-49.
Agency for Toxic Substances and Disease Registry. 2001. Toxicological Profile for Di-n-butyl Phthalate. Available: http://www.atsdr.cdc.gov/toxprofiles/ tp135.pdf [accessed 6 November 2014].
Agency for Toxic Substances and Disease Registry. 2002. Toxicological profile for Di(2-ethylhexyl) phthalate. Available: http://www.atsdr.cdc.gov/ toxprofiles/tp9.pdf [accessed 6 November 2014].
Boeniger MF, Lowry LK, Rosenberg J. 1993. Interpretation of urine results used to assess chemical exposure with emphasis on creatinine adjustments: a review. Am Ind Hyg Assoc J 54:615-627.
Braun JM, Smith KW, Williams PL, Calafat AM, Berry K, Ehrlich S, et al. 2012. Variability of urinary phthalate metabolite and bisphenol A concentrations before and during pregnancy. Environ Health Perspect 120:739-745; doi:10.1289/ehp.1104139.
Cantonwine DE, Cordero JF, Rivera-Gonzalez LO, Anzalota Del Toro LV, Ferguson KK, Mukherjee B, et al. 2013. Urinary phthalate metabolite concentrations among pregnant women in Northern Puerto Rico: distribution, temporal variability, and predictors. Environ Int 62:1-11.
Challis JR, Lockwood CJ, Myatt L, Norman JE, Strauss JF III, Petraglia F. 2009. Inflammation and pregnancy. Reprod Sci 16:206-215.
Chatterjee S, Hadi AS, Price B. 2000. Regression Analysis by Example. 3rd ed. Hoboken, NJ:John Wiley and Sons.
Erkekoglu P, Rachidi W, Yuzugullu OG, Giray B, Favier A, Ozturk M, et al. 2010. Evaluation of cytotoxicity and oxidative DNA damaging effects of di(2-ethylhexyl)-phthalate (DEHP) and mono(2-ethylhexyl)-phthalate (MEHP) on MA-10 Leydig cells and protection by selenium. Toxicol Appl Pharmacol 248:52-62.
Ferguson KK, Loch-Caruso R, Meeker JD. 2011. Urinary phthalate metabolites in relation to biomarkers of inflammation and oxidative stress: NHANES 1999-2006. Environ Res 111:718-726.
Ferguson KK, Loch-Caruso R, Meeker JD. 2012. Exploration of oxidative stress and inflammatory markers in relation to urinary phthalate metabolites: NHANES 1999-2006. Environ Sci Technol 46:477-485.
Ferguson KK, McElrath TF, Ko YA, Mukherjee B, Meeker JD. 2014a. Variability in urinary phthalate metabolite levels across pregnancy and sensitive windows of exposure for the risk of preterm birth. Environ Int 70:118-124.
Ferguson KK, McElrath TF, Meeker JD. 2014b. Environmental phthalate exposure and preterm birth. JAMA Pediatr 168:61-67.
Heazell AE, Moll SJ, Jones CJ, Baker PN, Crocker IP. 2007. Formation of syncytial knots is increased by hyperoxia, hypoxia and reactive oxygen species. Placenta 28(suppl A):S33-S40.
Hong YC, Park EY, Park MS, Ko JA, Oh SY, Kim H, et al. 2009. Community level exposure to chemicals and oxidative stress in adult population. Toxicol Lett 184:139-144.
Hornung RW, Reed L. 1990. Estimation of average concentration in the presence of nondetectable values. Appl Occup Environ Hyg 5:46-51.
Hurst CH, Waxman DJ. 2003. Activation of PPAR[alpha] and PPAR[gamma] by environmental phthalate monoesters. Toxicol Sci 74:297-308.
Il'yasova D, Scarbrough P, Spasojevic I. 2012. Urinary biomarkers of oxidative status. Clin Chim Acta 413:1446-1453.
Jiang Y, Scott AJ, Wild CJ. 2006. Secondary analysis of case-control data. Stat Med 25:1323-1339.
Kadiiska MB, Basu S, Brot N, Cooper C, Saari Csallany A, Davies MJ, et al. 2013. Biomarkers of oxidative stress study V: Ozone exposure of rats and its effect on lipids, proteins, and DNA in plasma and urine. Free Radic Biol Med 61:408-415.
Kasahara E, Sato EF, Miyoshi M, Konaka R, Hiramoto K, Sasaki J, et al. 2002. Role of oxidative stress in germ cell apoptosis induced by di(2ethylhexyl)phthalate. Biochem J 365(pt 3):849-856.
Kim JH, Park HY, Bae S, Lim YH, Hong YC. 2013. Diethylhexyl phthalates is associated with insulin resistance via oxidative stress in the elderly: a panel study. PLoS One 8:e71392; doi:10.1371/ journal.pone.0071392.
Kim SH, Kim SS, Kwon O, Sohn KH, Kwack SJ, Choi YW, et al. 2002. Effects of dibutyl phthalate and monobutyl phthalate on cytotoxicity and differentiation in cultured rat embryonic limb bud cells; protection by antioxidants. J Toxicol Environ Health A 65:461-472.
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(10):2390-2398.
Lin S, Ku HY, Su PH, Chen JW, Huang PC, Angerer J, et al. 2011. Phthalate exposure in pregnant women and their children in central Taiwan. Chemosphere 82:947-955.
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.
National Institute of Environmental Health Sciences. 2012. Biomarkers of Oxidative Stress Study. Available: http://www.niehs.nih.gov/research/ resources/databases/bosstudy/ [accessed 6 November 2014].
O'Brien ML, Cunningham ML, Spear BT, Glauert HP. 2001. Effects of peroxisome proliferators on glutathione and glutathione-related enzymes in rats and hamsters. Toxicol Appl Pharmacol 171:27-37.
Pinheiro J, Bates D, DebRoy S, Sarkar D, R Development Core Team. 2013. nlme: Linear and nonlinear mixed effects models. R package version 31-111. Available: http://cran.r-project.org/web/packages/ nlme/index.html [accessed 5 February 2015].
Potdar N, Singh R, Mistry V, Evans MD, Farmer PB, Konje JC, et al. 2009. First-trimester increase in oxidative stress and risk of small-for-gestational-age fetus. BJOG 116:637-642.
Roberts LJ, Morrow JD. 2000. Measurement of F(2)-isoprostanes as an index of oxidative stress in vivo. Free Radic Biol Med 28:505-513.
Rosado-Berrios CA, Velez C, Zayas B. 2011. Mitochondrial permeability and toxicity of diethylhexyl and monoethylhexyl phthalates on TK6 human lymphoblasts cells. Toxicol In Vitro 25:2010-2016.
Rosner B. 2011. Fundamentals of Biostatistics. 7th ed. Boston, MA:Brooks/Cole.
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.
Shono T, Taguchi T. 2014. Short-time exposure to mono-n-butyl phthalate (MBP)-induced oxidative stress associated with DNA damage and the atrophy of the testis in pubertal rats. Environ Sci Pollut Res Int 21:3187-3190.
Silva MJ, Samandar E, Preau JL Jr, Reidy JA, Needham LL, Calafat AM. 2007. Quantification of 22 phthalate metabolites in human urine. J Chromatogr B Analyt Technol Biomed Life Sci 860:106-112.
Smith KA, Shepherd J, Wakil A, Kilpatrick ES. 2011. A comparison of methods for the measurement of 8-isoPGF(2[alpha]): a marker of oxidative stress. Ann Clin Biochem 48(Pt 2):147-154.
Stein TP, Scholl TO, Schluter MD, Leskiw MJ, Chen X, Spur BW, et al. 2008. Oxidative stress early in pregnancy and pregnancy outcome. Free Radic Res 42:841-848.
Tetz LM, Cheng AA, Korte CS, Giese RW, Wang P, Harris C, et al. 2013. Mono-2-ethylhexyl phthalate induces oxidative stress responses in human placental cells in vitro. Toxicol Appl Pharmacol 268:47-54.
Vetrano AM, Laskin DL, Archer F, Syed K, Gray JP, Laskin JD, et al. 2010. Inflammatory effects of phthalates in neonatal neutrophils. Pediatr Res 68:134-139.
Woodruff TJ, Zota AR, Schwartz JM. 2011. Environmental chemicals in pregnant women in the United States: NHANES 2003-2004. Environ Health Perspect 119:878-885; doi:10.1289/ehp.1002727.
Wu LL, Chiou CC, Chang PY, Wu JT. 2004. Urinary 8-OHdG: a marker of oxidative stress to DNA and a risk factor for cancer, atherosclerosis and diabetics. Clin Chim Acta 339:1-9.
Zeman FA, Boudet C, Tack K, Floch Barneaud A, Brochot C, Pery AR, et al. 2013. Exposure assessment of phthalates in French pregnant women: results of the ELFE pilot study. Int J Hyg Environ Health 216:271-279.
Zhao Y, Ao H, Chen L, Sottas CM, Ge RS, Li L, et al. 2012. Mono-(2-ethylhexyl) phthalate affects the steroidogenesis in rat Leydig cells through provoking ROS perturbation. Toxicol In Vitro 26:950-955.
Zhou D, Wang H, Zhang J, Gao X, Zhao W, Zheng Y. 2010. Di-n-butyl phthalate (DBP) exposure induces oxidative damage in testes of adult rats. Syst Biol Reprod Med 56:413-419.
Zhou L, Beattie MC, Lin CY, Liu J, Traore K, Papadopoulos V, et al. 2013. Oxidative stress and phthalate-induced down-regulation of steroidogenesis in MA-10 Leydig cells. Reprod Toxicol 42:95-101.
Kelly K. Ferguson, (1) Thomas F. McElrath, (2) Yin-Hsiu Chen, (3) Bhramar Mukherjee, (3) and John D. Meeker (1)
(1) Department of Environmental Health Sciences, University of Michigan School of Public Health, Ann Arbor, Michigan, USA; (2) Division of Maternal-Fetal Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, USA; (3) Department of Biostatistics, University of Michigan School of Public Health, Ann Arbor, Michigan, USA
Address correspondence to J.D. Meeker, University of Michigan School of Public Health, Department of Environmental Health Sciences, 1835 SPH I, 1415 Washington Heights, Ann Arbor, MI 48109-2029 USA. Telephone: (734) 764-7184. E-mail: meekerj@ umich.edu
Supplemental Material is available online (http:// dx.doi.org/10.1289/ehp.1307996).
We thank K. Kneen, S. Clipper, G. Pace, D. Weller, and J. Bell of NSF International in Ann Arbor, Michigan, for urine phthalate analysis; E. Hurst and colleagues of Cayman Chemical in Ann Arbor, Michigan, for oxidative stress biomarker analysis; and R. Hauser of the Harvard School of Public Health for guidance on study design.
Funding was provided by the National Institute of Environmental Health Sciences, National Institutes of Health (grants R01ES018872, P42ES017198, P01ES022844, and P30ES017885).
The authors declare they have no actual or potential competing financial interests.
Received: December 2013; Accepted: 11 November 2014; Advance Publication: 14 November 2014; Final Publication: 1 March 2015.
Table 1. Demographic characteristics in weighted study population (n = 482). Characteristic Percent Race/ethnicity (n = 482) White 59 African American 16 Other 26 Education (n = 471) High school 14 Technical school 16 Junior college or 30 some college College graduate 41 Health insurance (n = 470) Private/HMO/self-pay 81 Medicaid/SSI/ 19 MassHealth Body mass index at visit 1 (n = 478) < 25 kg/[m.sup.2] 53 (underweight to normal) 25-30 kg/[m.sup.2] 27 (overweight) [greater than 20 or equal to] 30 kg/[m.sup.2] (obese to morbidly obese) Smoking during pregnancy (n = 476) Some 6 None 94 Alcohol use during pregnancy (n = 472) Some 5 None 95 Parity (n = 482) Nulliparous 45 Parous 55 Abbreviations: HMO, health maintenance organization; MassHealth, Massachusetts state health insurance provider; SSI, supplemental security income. Distributions of demographic characteristics were created from inverse probability weightings for case- control status. Table 2. Distributions of phthalate metabolites and oxidative stress biomarkers measured in urine samples collected from up to four time points during pregnancy in all samples measured from weighted population. Biomarker LOD % < LOD Geometric mean (geometric SD) MEHP ([micro]g/L) 1.0 4.7 10.6 (3.49) MEHHP ([micro]g/L) 0.1 0.0 34.2 (3.41) MEOHP ([micro]g/L) 0.1 0.1 18.3 (3.32) MECPP ([micro]g/L) 0.2 0.0 43.5 (3.40) [SIGMA]DEHP 0.39 (3.17) ([micro]mol/L) MBzP ([micro]g/L) 0.2 1.1 7.07 (3.05) MBP ([micro]g/L) 0.5 0.3 17.8 (2.40) MiBP ([micro]g/L) 0.1 0.4 7.61 (2.29) MEP ([micro]g/L) 1.0 0.1 141 (4.68) MCPP ([micro]g/L) 0.2 3.2 2.10 (3.14) 8-OHdG (ng/mL) 0.01 0.0 130 (1.66) 8-Isoprostane (pg/mL) 3.9 4.0 180 (2.64) Biomarker Percentile 25th 50th 75th 90th 95th Maximum MEHP ([micro]g/L) 4.63 9.07 21.0 56.3 106 1,555 MEHHP ([micro]g/L) 14.9 27.5 70.3 182 305 2,850 MEOHP ([micro]g/L) 8.33 15.3 37.6 93.0 152 1,128 MECPP ([micro]g/L) 17.7 34.9 99.4 231 391 3,713 [SIGMA]DEHP 0.17 0.31 0.78 2.00 3.01 21.1 ([micro]mol/L) MBzP ([micro]g/L) 3.47 6.38 13.2 29.0 55.8 465 MBP ([micro]g/L) 10.9 16.5 27.8 46.0 61.5 24,879 MiBP ([micro]g/L) 4.74 7.57 12.0 19.6 27.4 351 MEP ([micro]g/L) 47.3 121 383 1,084 2,307 48,130 MCPP ([micro]g/L) 1.03 1.68 3.50 8.60 19.6 848 8-OHdG (ng/mL) 98.4 130 173 233 288 1,339 8-Isoprostane (pg/mL) 130 210 320 460 574 2,784 Abbreviations: 8-OHdG, 8-hydroxydeoxyguanosine; LOD, limit of detection. All biomarkers were corrected for urinary specific gravity. For urinary phthalate metabolites, n = 1,693 samples, 482 subjects; for urinary oxidative stress biomarkers, n = 1,678 samples, 482 subjects. Table 3. Oxidative stress and exposure biomarker concentrations (geometric mean and geometric standard deviation) by categorical demographic characteristics in all samples measured from weighted population. Characteristic 8-OHdG (ng/mL) 8-Iso (pg/mL) Race/ethnicity White (reference) 130 (1.14) 153 (1.69) African American 133 (1.11) 277 (1.33) * Other 129 (1.15) 205 (1.43) * Education High school 146 (1.11) 289 (1.41) (reference) Technical school 133 (1.11) 217 (1.51) Junior college or 124 (1.17) * 173 (1.65) * some college College graduate 128 (1.13) 143 (1.57) * Health insurance Private 126 (1.14) 162 (1.64) insurance/HMO/ self-pay (reference) Medicaid/SSI/ 151 (1.09) * 271 (1.33) * MassHealth BMI at visit 1 < 25 kg/[m.sup.2] 128 (1.13) 155 (1.64) (reference) 25 to < 30 kg/ 132 (1.12) 181 (1.59) [m.sup.2] [greater than 136 (1.17) 270 (1.34) * or equal to] 30 kg/[m.sup.2] Tobacco use Smoked during 150 (1.11) 320 (1.47) pregnancy (reference) No smoking 129 (1.14) 173 (1.60) * during pregnancy Alcohol use Alcohol use during 138 (1.06) 201 (1.29) pregnancy (reference) No alcohol use 130 (1.14) 178 (1.63) during pregnancy Parity Nulliparous 129 (1.13) 166 (1.57) (reference) Parous 131 (1.14) 192 (1.62) * Characteristic MEHP MBP ([micro]g/L) ([micro]g/L) Race/ethnicity White (reference) 10.1 (2.20) 15.4 (1.37) African American 13.5 (2.34) * 23.6 (1.45) * Other 10.4 (2.02) 21.2 (1.63) * Education High school 9.20 (2.11) 27.8 (1.29) (reference) Technical school 9.91 (2.13) 19.6 (1.48) * Junior college or 9.57 (2.24) 16.7 (1.62) * some college College graduate 12.0 (2.11) 15.6 (1.37) * Health insurance Private 10.6 (2.19) 16.0 (1.47) insurance/HMO/ self-pay (reference) Medicaid/SSI/ 9.74 (2.05) 29.1 (1.28) * MassHealth BMI at visit 1 < 25 kg/[m.sup.2] 10.8 (2.19) 16.7 (1.37) (reference) 25 to < 30 kg/ 10.5 (2.27) 16.5 (1.54) [m.sup.2] [greater than 10.6 (2.10) 23.6 (1.56) * or equal to] 30 kg/[m.sup.2] Tobacco use Smoked during 9.87 (1.83) 26.4 (1.34) pregnancy (reference) No smoking 10.6 (2.19) 17.5 (1.47) * during pregnancy Alcohol use Alcohol use during 12.3 (2.27) 15.0 (1.19) pregnancy (reference) No alcohol use 10.5 (2.16) 18.1 (1.49) during pregnancy Parity Nulliparous 11.4 (2.23) 17.0 (1.39) (reference) Parous 10.1 (2.14) 18.6 (1.53) Abbreviations: 8-Iso, 8-isoprostane; 8-OHdG, 8-hydroxydeoxyguanosine; HMO, health maintenance organization; MassHealth, Massachusetts state health insurance provider; SSI, supplemental security income. All biomarkers were corrected for urinary specific gravity. For urinary phthalate metabolites, n = 1,693 samples, 482 subjects; for urinary oxidative stress biomarkers, n = 1,678 samples, 482 subjects. * p < 0.05 for significant difference in biomarker concentration from reference category, estimated from weighted linear mixed model with random intercepts for subject identification. Table 4. Percent difference (95% CIs) in oxidative stress biomarker in association with IQR increase in phthalate metabolite level. Metabolite IQR3 8-OHdG % difference p-Value (95% CI) MEHP 16.6 [micro]g/L 2.74 (-0.47, 6.05) 0.09 MEHHP 56.9 [micro]g/L 8.40 (4.93, 12.0) < 0.001 MEOHP 29.4 [micro]g/L 7.34 (4.01, 10.8) < 0.001 MECPP 80.5 [micro]g/L 6.53 (2.96, 10.2) < 0.001 IDEHP 0.63 [micro]mol/L 6.67 (3.23, 10.2) < 0.001 MBzP 12.5 [micro]g/L 20.7 (15.6, 26.1) < 0.001 MBP 24.8 [micro]g/L 18.1 (13.5, 22.9) < 0.001 MiBP 11.3 [micro]g/L 30.3 (24.4, 36.5) < 0.001 MEP 355 [micro]g/L 11.5 (7.32, 15.9) < 0.001 MCPP 2.98 [micro]g/L 7.23 (3.83, 10.7) < 0.001 Metabolite 8-Isoprostane % difference p-Value (95% CI) MEHP 14.1 (8.06, 20.5) < 0.001 MEHHP 15.8 (9.53, 22.4) < 0.001 MEOHP 15.9 (9.87, 22.3) < 0.001 MECPP 23.0 (16.0, 30.4) < 0.001 IDEHP 19.1 (12.7, 25.9) < 0.001 MBzP 42.7 (31.8, 54.4) < 0.001 MBP 42.0 (32.0, 52.7) < 0.001 MiBP 56.4 (43.9, 69.9) < 0.001 MEP 19.7 (11.8, 28.2) < 0.001 MCPP 20.2 (13.7, 27.1) < 0.001 8-OHdG, 8-hydroxydeoxyguanosine. Estimates were from adjusted linear mixed-effect models with random intercepts for subject identification (n = 1,604 samples, n = 464 subjects). The models were adjusted for urinary specific gravity, gestational age at sample collection, race/ethnicity, education level, health insurance provider, BMI (time-varying), time of day of urine sample collection (before vs. after 1300 hours, time-varying), and parity of infant. Models include inverse probability weights to adjust for case-control study design. "Ranges differ slightly from 25th-75th percentile differences in Table 2 because they were calculated from raw phthalate metabolite concentrations, uncorrected for urinary specific gravity.
|Printer friendly Cite/link Email Feedback|
|Author:||Ferguson, Kelly K.; McElrath, Thomas F.; Chen, Yin-Hsiu; Mukherjee, Bhramar; Meeker, John D.|
|Publication:||Environmental Health Perspectives|
|Date:||Mar 1, 2015|
|Previous Article:||Children's health in Latin America: the influence of environmental exposures.|
|Next Article:||Urinary metals and heart rate variability: a cross-sectional study of urban adults in Wuhan, China.|