Open access
Research Article
25 August 2017

Urinary Concentrations of Organophosphate Flame Retardant Metabolites and Pregnancy Outcomes among Women Undergoing in Vitro Fertilization

Publication: Environmental Health Perspectives
Volume 125, Issue 8
CID: 087018



Evidence from animal studies suggests that exposure to organophosphate flame retardants (PFRs) can disrupt endocrine function and impair embryo development. However, no epidemiologic studies have been conducted to evaluate effects on fertility and pregnancy outcomes.


We evaluated associations between urinary concentrations of PFR metabolites and outcomes of in vitro fertilization (IVF) treatment among couples recruited from an academic fertility clinic.


This analysis included 211 women enrolled in the Environment And Reproductive Health (EARTH) prospective cohort study (2005–2015) who provided one or two urine samples per IVF cycle. We measured five urinary PFR metabolites [bis(1,3-dichloro-2-propyl) phosphate (BDCIPP), diphenyl phosphate (DPHP), isopropylphenyl phenyl phosphate (ip-PPP), tert-butylphenyl phenyl phosphate (tb-PPP), and bis(1-chloro-2-propyl) phosphate (BCIPP)] using negative electrospray ionization liquid chromatography tandem mass spectrometry (LC-MS/MS). Molar concentrations of the urinary PFR metabolites were summed. We used multivariable generalized linear mixed models to evaluate the association of the PFR metabolites with IVF outcomes, accounting for multiple IVF cycles per woman.


Detection frequencies were high for BDCIPP (87%), DPHP (94%), and ip-PPP (80%), but low for tb-PPP (14%) and BCIPP (0%). We observed decreased success for several IVF outcomes across increasing quartiles of both summed and individual PFR metabolites (DPHP and ip-PPP) in our adjusted multivariable models. Significant declines in adjusted means from the lowest to highest quartile of ΣPFR were observed for the proportion of cycles resulting in successful fertilization (10% decrease), implantation (31%), clinical pregnancy (41%), and live birth (38%).


Using IVF to investigate human reproduction and pregnancy outcomes, we found that concentrations of some urinary PFR metabolites were negatively associated with proportions of successful fertilization, implantation, clinical pregnancy, and live birth.


One in six couples worldwide is affected by infertility, which is defined as the inability to get pregnant after 1 y or more of unprotected intercourse (Chandra et al. 2005), and a recent U.S. study found that pregnancy loss (miscarriage) affected approximately 28% of couples planning a pregnancy (Buck Louis et al. 2016). These figures will likely rise as the postponement of childbearing increases in developed regions of the world. Infertility has an associated health-care cost in the billions of dollars per year, not including the physical and psychological burden placed on the couple (Connolly et al. 2010). Both the high rates of fertility along with the associated high costs highlight the need to improve our understanding of risk factors that impair the ability to have a child.
One potential risk factor is environmental exposure. Several classes of endocrine-disrupting chemicals (EDCs) with widespread general population exposure, including pesticides and phthalates, have been linked to infertility and adverse pregnancy outcomes (Di Renzo et al. 2015). While there are hundreds or more EDCs, only a fraction have been evaluated for effects on infertility and pregnancy (Gore et al. 2015). Organophosphate flame retardants (PFRs) are a class of EDCs with ubiquitous exposure that have been detected in 90–100% of adult urine samples (Butt et al. 2014, 2016; Carignan et al. 2013a; Hammel et al. 2016; Hoffman et al. 2014; Meeker et al. 2013a). Over the past decade, PFRs have been used widely in the polyurethane foam of upholstered furniture as replacements for pentabromodiphenyl ether, a flame retardant mixture that was phased out of use in 2005 (Stapleton et al. 2009). These chemicals are not chemically bonded to foam and have been shown to migrate into the air and dust of indoor environments (van der Veen and de Boer 2012). The PFR triphenyl phosphate (TPHP) is used as part of flame retardant mixtures in polyurethane foam as well as in a variety of other applications, including as a plasticizer (WHO 1991).
Animal studies indicate that exposure to PFRs can disrupt endocrine function through altered thyroid action, steroidogenesis, or estrogen metabolism, and can also impair embryo development (Farhat et al. 2013; Liu et al. 2013; Wang et al. 2015). Reductions in sperm motility and increased serum total T3 levels were associated with increasing PFR exposures in a small study of men (Meeker et al. 2013b); however, no studies have investigated the effect of PFRs on pregnancy outcomes. Therefore, we explored the association between urinary concentrations of PFRs and pregnancy outcomes among women in a prospective cohort study, the Environment And Reproductive Health (EARTH) study, using assisted reproductive technologies (ART) as a model to study early developmental endpoints and pregnancy outcomes.



Study participants were women recruited into the EARTH study, which was established in June 2004 to evaluate environmental and dietary determinants of fertility from patients undergoing ART at the Massachusetts General Hospital (MGH) Fertility Center. Women between the ages of 18 and 46 were eligible to participate, and approximately 60% of those contacted by the research nurses enrolled in the study. The EARTH study was approved by the Human Studies Institutional Review Boards of the MGH and Harvard T.H. Chan School of Public Health. Participants signed an informed consent after the study procedures were explained by trained study staff and any questions were answered. Demographic information including race/ethnicity, smoking history, and education, as well as whether they had a previous pregnancy, were recorded at study entry by the participant using a questionnaire. To be included in the present analysis, women must have contributed their own oocytes and at least one urine sample for the measurement of flame retardant metabolites during an in vitro fertilization (IVF) cycle. From the 212 women (298 IVF cycles) who met these criteria, we excluded 1 woman (1 IVF cycle) with incomplete outcome data. Our final data set included 211 women with 297 IVF cycles who had complete information on the exposure and outcome variables.

Clinical Data and in Vitro Fertilization Outcomes

At study entry, the participant’s date of birth was collected, and her weight and height were measured by trained study staff. Body mass index (BMI) was calculated as weight (in kilograms) per height (in meters) squared. Clinical information was collected or abstracted from the patient’s electronic medical record by trained study staff at enrollment into the study and after each IVF cycle. Follicle stimulating hormone (FSH) was measured in a blood sample drawn on the third day of the menstrual cycle and analyzed for with an automated electrochemiluminescence immunoassay at the MGH Core Laboratory as previously described (Mok-Lin et al. 2010). Cause of infertility was diagnosed by a physician at the MGH Fertility Center according to the Society for Assisted Reproductive Technology (SART) using standard infertility definitions (SART 2014; Mok-Lin et al. 2010). IVF treatment protocols were assigned by a physician at the MGH Fertility Center based on clinical indications and factors such as age and infertility diagnosis. Treatment protocols include: a) luteal phase gonadotropin-releasing hormone (GnRH) agonist (low-, regular-, or high-dose leuprolide acetate; Lupron); b) follicular phase GnRH agonist/flare stimulation; or c) GnRH antagonist. Lupron dose was reduced at, or shortly after, the start of ovarian stimulation with FSH/human menopausal gonadotropin (hMG) in the luteal phase GnRH agonist protocol. Serum peak estradiol (pmol/L) was measured on the day of ovulation trigger with human chorionic gonadotropin (hCG) [estradiol (E2) trigger levels] using an automated electrochemiluminescence immunoassay at the MGH Core Laboratory. Oocytes were counted and classified by embryologists after egg retrieval as germinal vesicle, metaphase I, metaphase II, or degenerated. Fertilization was obtained by conventional IVF or intracytoplasmatic sperm injection (ICSI). Fertilization was confirmed 17–20 h after insemination by the presence of a fertilized oocyte with two pronuclei. Embryos were monitored for cell number and morphological quality [1 (best) to 5 (worst)] on days 2 and 3, and considered of best quality if they had a score of 1 or 2. Transfer occurred rarely on day 2, with the majority on day 3 or 5 and comprised of cleavage stage and blastocyst stage embryos, respectively. Most women had 1–3 embryos transferred (range=05), with the number dependent on the woman’s age, cycle number. and the day of transfer. Implantation was defined as a serum β-hCGlevel>6mIU/mL approximately 17 d (range=1520d) after egg retrieval, clinical pregnancy as the presence of an intrauterine pregnancy confirmed by ultrasound at approximately 6-wk gestation, and live birth as the birth of a neonate on or after 24-wk gestation.

Organophosphate Flame Retardant Assessment in Urine Samples

Urine samples were collected from participants enrolled in the EARTH study between May 2005 and January 2015. Each sample was collected into a sterile polypropylene cup, and up to two urine samples were collected during each IVF cycle. Urine samples were collected at a geometric mean (GM) [95% confidence interval (CI)] of 7.6 (7.3, 8.0) d apart with 178 (141, 215) d between the first and second IVF cycle. Following collection of each sample, specific gravity (SG) was measured using a Protometer (hand held) Model 100B refractometer (National Instrument Company, Inc.); the urine sample was divided into aliquots and frozen at 80°C. Samples were shipped on dry ice overnight to H.M. Stapleton’s lab at Duke University (Durham, NC) for the quantification of the PFR metabolites (Figure 1).
Figure 1. Organophosphate flame retardant parent compound and primary urinary metabolite.
Extraction and analysis methods for bis(1-chloro-2-propyl) phosphate (BCIPP), bis(1,3-dichloro-2-propyl) phosphate (BDCIPP), diphenyl phosphate (DPHP), isopropylphenyl phenyl phosphate (ip-PPP), and tert-butylphenyl phenyl phosphate (tb-PPP) followed methods previously developed by H.M. Stapleton’s laboratory (Butt et al. 2014). Briefly, urine samples were thawed, and a 2.5- to 5-mL aliquot was transferred to a clean glass test tube, where it was spiked with mass-labeled internal standards (d10-BDCIPP=80ng, d10-DPHP=60ng). After acidifying to pH <6.5 with formic acid, samples were diluted 1:1 with water and concentrated and cleaned using solid-phase extraction techniques (SPE). The SPE eluent was blown to dryness under a gentle nitrogen stream, reconstituted in 500μL of 1:1 H2O:MeOH and spiked with the recovery standard (13C2-DPHP=81.5ng). Extracts were analyzed by negative electrospray ionization liquid chromatography tandem mass spectrometry (LC-MS/MS) as previously described (Butt et al. 2014). Chromatography was achieved under gradient conditions using a Luna C18(2) column (50×2.0mm, 2.5μm particle size, Phenomenex) preceded by a SecurityGuard Polar-RP (4×2.0mm) guard cartridge (Phenomenex). The mobile phases were methanol and water (modified with 0.8mM ammonium acetate), flow rate was 300μL/min, the injection volume was 5μL, and the column oven was 45°C. Data were acquired under multiple reaction monitoring conditions using optimized parameters. Analyte responses were normalized to internal standard responses. BCIPP and BDCIPP were normalized using d10-BDCIPP, while DPHP, ip-PPP, and tb-PPP were normalized using d10-DPHP. Urinary SG ranged from 1.002 to 1.028 with a mean of 1.016.
Samples were analyzed by LC-MS/MS in 10 separate batches, and unique method detection limits (MDLs) were calculated for each analysis batch. In the urine samples, the mean recovery of the mass-labeled standards was 119% (standarderror=0.75%) for d10-DPHP and 152% (2.2%) for d10-BDCIPP. One laboratory blank (5mL MilliPore water only) sample was extracted with every batch (n=95). An in-house standard reference material (SRM) was prepared from pooled urine that was collected during previous studies. SRM samples were periodically analyzed during the extraction batches (n=18) and were generally within 10% for DPHP, 15% for BDCIPP, and 20% for ip-PPP. Two of the individual subsamples were analyzed in duplicate to assess method precision and were generally within 15% for DPHP, 25% for ip-PPP, and 35% for BDCIPP. Very low levels of DPHP (mean=0.58ng) and ip-PPP (mean=0.21ng) were commonly detected in the laboratory blanks, and analyte values were blank corrected using the mean laboratory blank values. MDLs were calculated as three times the standard deviation of laboratory blanks normalized to the volume of water extracted (5mL). MDLs ranged (n=10) from 68180pg/mL for BCIPP, 31300pg/mL for BDCIPP, 25130pg/mL for DPHP, and 23120pg/mL for ip-PPP, and 10150pg/mL for tb-PPP, respectively.

Data Analysis

Demographic and clinical characteristics were reported using mean±interquartilerange or percentages. Unquantified concentrations <MDL were substituted with a value equal to the MDL/2 (Hornung and Reed 1990). To account for urinary dilution, we adjusted to SG as described by Pearson et al. (2009): CSG=C*[(SGm1)/(SGi1)] where CSG=SG-adjustedurinarymetaboliteconcentration, C=urinarymetaboliteconcentration, SGm=meanSG for the population, and SGi=SG for an individual sample. Cycle-specific urinary metabolite concentrations were calculated using the GM of the two urinary metabolite concentrations from each IVF cycle. Cycle-specific concentrations were divided into quartiles for use in regression models. All analyses used the SG-adjusted urinary metabolite concentrations.
To evaluate the associations between the urinary metabolites and IVF outcomes, we fit multivariable generalized linear mixed models with random intercepts to account for multiple IVF cycles in the same woman. These models allow for the use of multiple outcome observations per individual while accounting for within-person correlations in outcomes. These models are also appropriate and can provide unbiased estimates in the presence of an unbalanced design (e.g., different number of cycles contributed per woman) when imbalance in the number of IVF cycles is not completely random (e.g., women with more cycles are having more difficulty getting pregnant), and the lack of balance can be accurately predicted by all measured covariates in the adjusted model. A normal distribution and identity link function was specified for peak estradiol and endometrial lining thickness. A Poisson distribution and log link function were specified for the number of mature and total oocytes as well as best quality embryo. A binomial distribution and logit link function were spcified for fertilization and the proportion of mature to total oocytes. Finally, a binary distribution with a logit link function were specified for the clinical outcomes (implantation, clinical pregnancy, and live birth). Tests for trend were conducted across quartiles using the median log-transformed urinary metabolite concentration in each quartile as a linear variable in the regression models. To allow for better interpretation of the results, all results are presented as population marginal means adjusted for covariates. Percent decrease was calculated as the difference in marginal means from Q1 to Q4 divided by the marginal mean from Q1.
We evaluated confounding using prior knowledge and descriptive statistics from our cohort. The following covariates were considered for inclusion in the final model: maternal age (continuous), race/ethnic group (black/Asian/other, white/Caucasian), BMI (continuous), smoking history (ever, never), education (high school, college graduate, graduate degree), year of IVF treatment cycle (continuous), and primary infertility diagnosis (female factor, male factor, and unexplained). Variables were included in the final model if they were associated with any of the individual PFR metabolites (BDCIPP, DPHP, ip-PPP) in our population, were suspected to be associated with exposure based on previous research, or were strong predictors of the outcome. In our primary analysis of the early developmental outcomes, we excluded seven women (16 IVF cycles) with unsuccessful oocyte retrieval, and applied this restriction in a sensitivity analysis for the clinical outcomes. We tested whether the associations of urinary metabolite concentrations with fertilization were modified by ICSI by entering a product of the metabolite quartiles and a binary variable representing the presence or absence of ICSI into the models. Given our limited sample size for detecting interactions, a suggestion of interaction was considered if the p-value for this interaction term was <0.10. Finally, we tested for a potential cohort effect by controlling for maternal year of birth, rather than maternal age, in our adjusted model.
We evaluated the relationship between urinary PFR metabolites using Spearman’s correlation and estimated the variability of urinary metabolite concentrations within a cycle and for all cycles by calculating the intraclass correlation coefficient (ICC). We conducted all statistical analyses using SAS (version 9.2; SAS Institute Inc.), and all tests other than interaction with ICSI considered two-sided significance levels less than 0.05 as statistically significant.


Study Population

Our analysis included 211 women who were on average 35.0 y of age; 87% were Caucasian, 75% had never smoked, and 34% had a prior pregnancy (Table 1). The primary SART diagnosis was approximately equally distributed between female factor (36%), male factor (29%), and unexplained (35%).
Table 1 Demographic and reproductive characteristics of 211 women [297 in vitro fertilization (IVF) cycles] in the Environment and Reproductive Health (EARTH) study.
CharacteristicsTotal cohort median (IQR) or n (%)
Demographic characteristics 
 Age, years35.0 (32.0, 37.0)
 Race/ethnic group, n (%) 
  Black/Asian/other28 (13)
  White/Caucasian183 (87)
 Body mass index, kg/m222.8 (20.8, 25.8)
 Ever smoker, n (%)52 (25)
 Education, n (%)a 
  High school/some college12 (6)
  College graduate71 (37)
  Graduate degree109 (57)
Reproductive/cycle characteristics 
 History of ever been pregnant71 (34)
 Day 3 FSH levels, IU/L6.7 (5.9, 8.1)
 Initial infertility diagnosis, n (%) 
  Female factor75 (36)
  Male factor62 (29)
  Unexplained74 (35)
 Previous intrauterine insemination, n (%)a85 (41)
 Previous in vitro fertilization, n (%)a43 (21)
 Treatment protocol, n (%) 
  Antagonist37 (13)
  Flareb54 (18)
  Luteal phase agonistc206 (69)
 E2 trigger levels, pmol/L1,965 (1,517, 2,570)
 ICSI cycles, n (%)167 (56)
 Cycle canceled prior to transfer, n (%)16 (5)
 Embryo transfer day, n (%) 
  Number of embryos transferred16 (5)
  Day 216 (5)
  Day 3163 (55)
  Day 5102 (35)
 Number of embryos transferred, n (%) 
  0 embryos16 (5)
  1 embryo56 (19)
  2 embryos162 (55)
  3 embryos63 (21)
Note: E2, estradiol; FSH, follicle stimulating hormone; ICSI, intracytoplasmatic sperm injection; IQR, interquartile range.
Has missing data (19 missing education; 3 missing previous IUI and IVF).
Follicular phase GnRH agonist/flare protocol.
Luteal phase GnRH agonist protocol.

Urinary Organophosphate Flame Retardant Metabolites

Detection frequencies were high for BDCIPP (87%), DPHP (94%), and ip-PPP (80%), but low for tb-PPP (14%) and BCIPP (0%) (Table 2). Concentrations of DPHP were 19% higher than BDCIPP on average and 252% higher than ip-PPP. Due to the pmol/L low detection frequencies for tb-PPP and BCIPP, central tendency values were not calculated and statistical analysis was not performed for these urinary metabolites. BDCIPP was weakly but significantly correlated with DPHP (rs=0.21, p=0.002) and ip-PPP (rs=0.25, p=0.0003); DPHP and ip-PPP were also weakly correlated (rs=0.19, p=0.006). ICCs indicated moderate within-person variability across cycles (BDCIPP=0.36; DPHP=0.34; ip-PPP=0.37), with somewhat higher reproducibility within a cycle (BDCIPP=0.48; DPHP=0.53; ip-PPP=0.44).
Table 2 Distribution of urinary organophosphate flame retardant metabolites (ug/L) measured among 211 women with 297 cycles and 563 urine samples in the Environment And Reproductive Health (EARTH) study.
SampleN>MDLGM (95% CI)Min10th pctl25th pctl50th pctl75th pctl90th pctl95th pctlMax
Specific gravity adjusteda          
 BDCIPP487 (87)0.68 (0.62, 0.75)<MDL<MDL0.320.691.403.305.0863.4
 DPHP527 (94)0.81 (0.75, 0.89)<MDL<MDL0.430.751.283.105.25616
 ip-PPP450 (80)0.23 (0.21, 0.25)<MDL<MDL0.120.240.440.771.1871.4
 tb-PPP81 (14)NA<MDL<MDL<MDL<MDL<MDL0.270.38351
Wet weight (not specific gravity adjusted)          
 BDCIPP487 (87)0.66 (0.59, 0.73)<MDL<MDL0.290.691.683.655.7663.4
 DPHP527 (94)0.78 (0.71, 0.87)<MDL<MDL0.400.791.493.505.92657
 ip-PPP450 (80)0.22 (0.20, 0.25)<MDL<MDL0.090.240.490.941.3779.5
 tb-PPP81 (14)NA<MDL<MDL<MDL<MDL<MDL0.110.25374
Note: All values below MDL were assigned a value equal to the MDL divided by 2. BCIPP, bis(1-chloro-2-propyl) phosphate; BDCIPP, bis(1,3-dichloro-2-propyl)phosphate; CI, confidence interval; DPHP, diphenyl phosphate; GM, geometric mean; ip-PPP, isopropylphenyl phenyl phosphate; Max, maximum; <MDL, methoddetectionlimit; Min, minimum; NA, not applicable; pctl, percentile; tb-PPP, tert-butylphenyl phenyl phosphate.
Adjusted to specific gravity, range (1.001–1.10).
Demographic variables associated with urinary PFR concentrations included age (BDCIPP and DPHP), BMI (BDCIPP), and year (ip-PPP). For every unit increase in age, there was a 4% (p=0.05) and 3% (p=0.04) decline in mean urinary BDCIPP and DPHP concentrations. For every unit increase in BMI, there was a 4% (p=0.03) increase in urinary BDCIPP, and every year, there was a 7% (p=0.002) decline in urinary ip-PPP concentrations. There were no other significant associations with the demographic variables, and no associations of urinary PFR quartiles with reproductive/cycle characteristics were observed.

In Vitro Fertilization Outcomes

We observed significant declines in adjusted means from the lowest to highest quartile of ΣPFR for the proportion of cycles resulting in successful fertilization (10% decrease, p-trend=0.04), implantation (31%, p-trend=0.02), clinical pregnancy (41%, p-trend=0.004) and live birth (38%, p-trend=0.05), adjusted for age, BMI, race/ethnicity, year of cycle, and infertility diagnosis (Figure 2, Tables 3 and 4). There were also significant declines in adjusted means from the highest to lowest quartile of DPHP for the proportion of cycles resulting in successful implantation (28% decrease, p-trend=0.02) and clinical pregnancy (36%, p-trend=0.01) as well as for ip-PPP for the proportion of cycles resulting in successful fertilization (16%, p-trend=0.0006), implantation (27%, p-trend=0.05), and live birth (34%, p-trend=0.05). Unadjusted results were similar (Table S1). Results for the clinical outcomes were similar when restricted to IVF cycles with successful oocyte retrieval, restricted to the first cycle in the EARTH study, restricted to nulliparous women, or controlling for maternal year of birth (Table S2).
Figure 2. Adjusted mean [95% confidence interval (CI)] proportion of cycles resulting in implantation, live birth, and clinical pregnancy by quartile of urinary organophosphate flame retardant (PFR) metabolite concentrations among 211 women undergoing 297 in vitro fertilization (IVF) cycles. Gray shading indicates change in means from the first and fourth quartile. Adjusted models control for maternal age (continuous), body mass index (BMI) (continuous), race/ethnicity (black/Asian/other, white/Caucasian), year of IVF treatment cycle (continuous), and primary Society for Assisted Reproductive Technology (SART) infertility diagnosis at study entry (female, male, unknown), with continuous variables at their mean level and categorical variables weighted by their frequency in the study population. *Significantly different from the lowest quartile (Q1) at the α=0.05 level.
Table 3 Adjusted mean [95% confidence interval (CI)] proportion of cycles resulting in implantation, live birth, and clinical pregnancy by quartile of urinary organophosphate flame retardant metabolite concentrations among 211 women undergoing 297 in vitro fertilization (IVF) cycles.
 Q10.71 (0.59, 0.81)0.72 (0.59, 0.81)0.69 (0.56, 0.79)0.74 (0.62, 0.83)
 Q20.62 (0.49, 0.73)0.58 (0.46, 0.70)0.64 (0.52, 0.75)0.50 (0.37, 0.62)
 Q30.57 (0.45, 0.69)0.53 (0.41, 0.65)0.56 (0.44, 0.68)0.61 (0.48, 0.72)
 Q40.49 (0.37, 0.62)*0.56 (0.43, 0.67)0.50 (0.38, 0.62)*0.54 (0.41, 0.66)*
 Q10.63 (0.50, 0.74)0.61 (0.49, 0.73)0.61 (0.49, 0.72)0.63 (0.51, 0.74)
 Q20.54 (0.42, 0.66)0.50 (0.38, 0.62)0.55 (0.42, 0.66)0.40 (0.29, 0.53)
 Q30.50 (0.38, 0.62)0.47 (0.35, 0.59)0.48 (0.36, 0.60)0.51 (0.39, 0.63)
 Q40.37 (0.26, 0.50)*0.45 (0.33, 0.57)0.39 (0.28, 0.52)*0.48 (0.36, 0.61)
Live birth    
 Q10.48 (0.35, 0.61)0.49 (0.36, 0.61)0.47 (0.34, 0.59)0.53 (0.40, 0.65)
 Q20.42 (0.30, 0.55)0.36 (0.24, 0.48)0.42 (0.30, 0.55)0.35 (0.24, 0.47)
 Q30.38 (0.27, 0.51)0.36 (0.25, 0.49)0.37 (0.26, 0.50)0.35 (0.24, 0.48)
 Q40.30 (0.20, 0.43)*0.36 (0.25, 0.49)0.31 (0.21, 0.44)0.35 (0.24, 0.48)
Note: Adjusted models control for maternal age (continuous), body mass index (continuous), race/ethnicity (black/Asian/other, white/Caucasian), year of IVF treatment cycle (continuous), and primary Society for Assisted Reproductive Technology (SART) infertility diagnosis at study entry (female, male, unknown). Adjusted means are presented for the mean maternal age (35.2), body mass index (23.9), race/ethnicity (white), year of IVF treatment cycle (2010), and primary SART infertility diagnosis at study entry (female=1 male=0, unexplained=0). BDCIPP, bis(1,3-dichloro-2-propyl) phosphate; DPHP, diphenyl phosphate; ip-PPP, isopropylphenyl phenyl phosphate; PFR, organophosphate flame retardant.
Significantly different from Q1 at the p=0.05 level.
Table 4 Adjusted means [95% confidence interval (CI)] for early developmental outcomes by quartile of urinary organophosphate flame retardant metabolite concentrations among 205 women undergoing 281 in vitro fertilization (IVF) cycles that had oocyte retrieval.
E2 trigger levels, pmol/L    
 Q11,939 (1,766, 2,112)1,995 (1,821, 2,169)1,969 (1,794, 2,145)2,028 (1,855, 2,200)
 Q22,163 (1,992, 2,333)2,061 (1,889, 2,233)2,120 (1,946, 2,293)2,097 (1,927, 2,267)
 Q32,007 (1,838, 2,176)2,015 (1,842, 2,188)2,143 (1,970, 2,317)2,153 (1,985, 2,322)
 Q42,123 (1,947, 2,299)2,155 (1980, 2,331)1,994 (1,817, 2,171)1,950 (1,773, 2,128)
Total oocytes, count    
 Q110.6 (9.6, 11.8)9.9 (8.9, 11.0)11.0 (9.9, 12.2)10.7 (9.6, 11.9)
 Q211.1 (10.0, 12.3)11.4 (10.3, 12.7)11.1 (10.0, 12.3)10.9 (9.8, 12.1)
 Q310.8 (9.8, 12.0)11.1 (10.0, 12.3)11.3 (10.2, 12.6)11.9 (10.7, 13.1)
 Q411.5 (10.3, 12.8)11.7 (10.6, 13.0)*10.7 (9.6, 11.9)10.7 (9.5, 11.9)
Total MII oocytes, count    
 Q19.2 (8.2, 10.3)8.5 (7.6, 9.5)9.3 (8.3, 10.4)8.9 (7.9, 9.9)
 Q29.4 (8.4, 10.5)9.4 (8.4, 10.5)9.5 (8.5, 10.6)9.3 (8.4, 10.4)
 Q38.8 (7.9, 9.8)9.3 (8.3, 10.3)9.4 (8.4, 10.5)10.1 (9.1, 11.2)
 Q49.6 (8.6, 10.7)9.8 (8.8, 10.9)8.8 (7.8, 9.9)8.6 (7.7, 9.7)
Fertilization, proportion    
 Q10.77 (0.72, 0.81)0.75 (0.70, 0.79)0.75 (0.70, 0.79)0.77 (0.72, 0.81)
 Q20.72 (0.67, 0.76)0.73 (0.68, 0.77)0.68 (0.63, 0.73)0.75 (0.71, 0.79)
 Q30.71 (0.66, 0.76)0.71 (0.66, 0.76)0.78 (0.73, 0.82)0.71 (0.66, 0.75)
 Q40.70 (0.65, 0.74)*0.71 (0.66, 0.76)0.68 (0.63, 0.73)*0.65 (0.60, 0.70)*
Best quality embryos, count    
 Q11.6 (1.3, 2.1)1.5 (1.1, 1.9)1.6 (1.2, 2.1)1.5 (1.1, 1.9)
 Q21.5 (1.2, 2.0)1.3 (1.0, 1.7)1.4 (1.1, 1.9)1.8 (1.4, 2.3)
 Q31.4 (1.0, 1.8)1.6 (1.3, 2.1)1.6 (1.3, 2.1)1.5 (1.2, 2.0)
 Q41.5 (1.1, 1.9)1.5 (1.2, 2.0)1.4 (1.0, 1.8)1.2 (0.9, 1.6)
Endometrial wall thickness, mm    
 Q110.3 (9.79, 10.8)10.2 (9.7, 10.7)10.2 (9.7, 10.8)10.1 (9.6, 10.6)
 Q210.3 (9.79, 10.8)10.2 (9.7, 10.7)10.2 (9.7, 10.7)10.1 (9.6, 10.6)
 Q39.80 (9.31, 10.3)9.9 (9.4, 10.4)10.1 (9.6, 10.6)9.8 (9.3, 10.3)
 Q49.94 (9.42, 10.5)10.0 (9.5, 10.5)9.8 (9.3, 10.3)10.3 (9.8, 10.8)
Note: Adjusted models control for maternal age (continuous), body mass index (continuous), race/ethnicity (black/Asian/other, white/Caucasian), year of IVF treatment cycle (continuous), and primary Society for Assisted Reproductive Technology (SART) infertility diagnosis at study entry (female, male, unknown). Adjusted means are presented for the mean maternal age (35.2), body mass index (24.0), race/ethnicity (white), year of IVF treatment cycle (2010), and primary SART infertility diagnosis at study entry (female=1, male=0 unexplained=0). BDCIPP, bis(1,3-dichloro-2-propyl) phosphate; DPHP, diphenyl phosphate; ip-PPP, isopropylphenyl phenyl phosphate; PFR, organophosphate flame retardant.
Significantly different from Q1 at the p=0.05 level.
The proportion of fertilized oocytes was significantly lower in the highest quartile (Q4) of urinary ip-PPP compared to the lowest (Q1) with an adjusted difference in proportions (95% CI) of 0.10 (0.02, 0.19) using conventional fertilization that was smaller and nonsignificant when restricted to fertilization using ICSI [adjusted difference in proportions=0.04 (95%CI=0.02, 0.10)]. However, the interaction term by ICSI on the association of urinary ip-PPP with fertilization was not statistically significant (p of interaction=0.53).
There was a significant positive association between urinary BDCIPP and the number of total oocytes (18% increase, p-trend=0.04) and a nonsignificant increase in the number of mature oocytes (15% increase, p-trend=0.10) retrieved in an IVF cycle. For endometrial lining thickness, we observed a nonsignificant decline in adjusted means from the highest to lowest quartile of urinary DPHP (4% decrease, p-trend=0.15). Similar trends with oocyte count were not observed for DPHP, ip-PPP, or ΣPFR (p-trends >0.3). No other associations were observed between urinary metabolites and the other early developmental outcomes, and unadjusted results were similar (Table S3).


As far as we are aware, this is the first study to explore associations in an epidemiologic study on the effects of PFRs on female reproduction. We used the model of IVF to investigate human reproduction and early pregnancy outcomes, ranging chronologically from oocyte retrieval, oocyte fertilization, embryo quality, and implantation to clinical pregnancy and live birth. The sum of the urinary PFR metabolites was associated with reduced probability of successful fertilization, implantation, clinical pregnancy, and live birth. These findings are clinically relevant, as the adjusted proportion of live births in the highest quartile of ΣPFR was 0.30 as compared to 0.48 in the lowest quartile (adjusted difference in proportions=0.18). GM concentrations of the PFR metabolites in our study population were similar to or lower than other adult populations in the United States, Norway, and Australia; thus, these exposure levels are not abnormal or high (Butt et al. 2014; Carignan et al. 2013b; Cequier et al. 2015; Cooper et al. 2011; Dodson et al. 2014; Hoffman et al. 2015; Meeker et al. 2013b; Van den Eede et al. 2015) (Table S4).
Animal studies indicate that PFRs may adversely affect female reproduction through disruption of regulatory pathways mediated by the hypothalamus–pituitary–gonadal axis. For example, studies in zebrafish have observed decreased fecundity (hatching and survival) with independent exposures to TDCIPP and TPHP (Liu et al. 2013; Wang et al. 2015). Evidence of hormone disruption in zebrafish includes increased plasma E2, testosterone, and vitellogenin as well as increased triiodothyronine (T3) and decreased thyroxine (T4) (Liu et al. 2013; Wang et al. 2015, 2013). Studies of chicken embryos have observed delayed hatching as well as declines in plasma thyroxine and cholesterol (Farhat et al. 2014, 2013). Human studies have also observed an association between low pregnancy levels of estradiol and fetal loss (Schindler 2004), and it is well known that subclinical hypothyroidism can adversely affect fertility (Abdel Rahman et al. 2010; Bussen et al. 2000; Scoccia et al. 2012; Velkeniers et al. 2013).
Strengths of our study include the prospective study design, preconception measurement of exposure using repeated urine samples (which is necessary for short-lived biomarkers like the PFR metabolites), state-of-the-art measurement of PFR exposure biomarkers, assessment of early developmental outcomes (i.e., fertilization, implantation) that are not observable in non-IVF populations, clinical outcome data obtained from electronic medical records, and control for potential confounders. One limitation of this analysis is that we did not consider the male partner's exposure, which may be correlated with his female partner and could contribute to the observed association. Our findings are generalizable to the infertility clinic population, which is sizable (Thoma et al. 2013), and may be more broadly generalizable, assuming that women undergoing IVF have similar biological responses to PFR exposure as women not undergoing IVF. Our findings may also be more relevant and generalizable to older women, as the mean age of women in our study population was 35 y. As each outcome is dependent on the previous, a larger sample would be required to differentiate between independent effects on implantation, clinical pregnancy, and live birth.
In conclusion, using IVF as a model to investigate human reproduction and pregnancy outcomes, we found that concentrations of some urinary PFR metabolites were inversely associated with proportions of successful fertilization, implantation, clinical pregnancy, and live birth. These results highlight the potential reproductive effects of low-level exposure (i.e., background exposure levels of the general population) to PFRs and adverse IVF outcomes. Future research should focus on potential interactions between PFRs and other chemicals (i.e., mixtures analysis) that adversely affect reproductive health and also explore the potential effect of PFRs on male reproductive health.


The authors gratefully acknowledge S. Pecha for sample extraction and all members of the EARTH study team, specifically the Harvard T.H. Chan School of Public Health research nurses J.B. Ford and M.G. Keller, research staff R. Dadd and P. Morey, physicians and staff at Massachusetts General Hospital fertility center, and a special thanks to all the study participants.
The work was supported by grants ES009718, ES022955, ES000002, and T32ES007069 from the National Institute of Environmental Health Sciences (NIEHS).

Article Notes

Supplemental Material is available online (
The authors declare they have no actual or potential competing financial interests.
Note 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 [email protected]. Our staff will work with you to assess and meet your accessibility needs within 3 working days.

Supplementary Material

File (ehp1021.s001.acco.pdf)
File (ehp1021.smcontents.508.pdf)


Abdel Rahman AH, Aly Abbassy H, Abbassy AA. 2010. Improved in vitro fertilization outcomes after treatment of subclinical hypothyroidism in infertile women. Endocr Pract 16(5):792–797.
Buck Louis GM, Sapra KJ, Schisterman EF, Lynch CD, Maisog JM, Grantz KL, et al. 2016. Lifestyle and pregnancy loss in a contemporary cohort of women recruited before conception: The LIFE Study. Fertil Steril 106(1):180–188.
Bussen S, Steck T, Dietl J. 2000. Increased prevalence of thyroid antibodies in euthyroid women with a history of recurrent in-vitro fertilization failure. Hum Reprod 15(3):545–548.
Butt CM, Congleton J, Hoffman K, Fang M, Stapleton HM. 2014. Metabolites of organophosphate flame retardants and 2-ethylhexyl tetrabromobenzoate in urine from paired mothers and toddlers. Environ Sci Technol 48(17):10432–10438.
Butt CM, Hoffman K, Chen A, Lorenzo A, Congleton J, Stapleton HM. 2016. Regional comparison of organophosphate flame retardant (PFR) urinary metabolites and tetrabromobenzoic acid (TBBA) in mother-toddler pairs from California and New Jersey. Environ Int 94:627–634.
Carignan CC, McClean MD, Cooper EM, Watkins DJ, Fraser AJ, Heiger-Bernays W, et al. 2013. Predictors of tris(1,3-dichloro-2-propyl) phosphate metabolite in the urine of office workers. Environ Int 55:56–61.
Cequier E, Sakhi AK, Marcé RM, Becher G, Thomsen C. 2015. Human exposure pathways to organophosphate triesters - a biomonitoring study of mother-child pairs. Environ Int 75:159–165.
Chandra A, Martinez GM, Mosher WD, Abma JC, Jones J. 2005. Fertility, family planning, and reproductive health of U.S. women: data from the 2002 National Survey of Family Growth. Vital Health Stat 23(25):1–160.
Connolly MP, Hoorens S, Chambers GM. 2010. The costs and consequences of assisted reproductive technology: an economic perspective. Hum Reprod Update 16(6):603–613.
Cooper EM, Covaci A, van Nuijs AL, Webster TF, Stapleton HM. 2011. Analysis of the flame retardant metabolites bis(1,3-dichloro-2-propyl) phosphate (BDCPP) and diphenyl phosphate (DPP) in urine using liquid chromatography-tandem mass spectrometry. Anal Bioanal Chem 401(7):2123–2132.
Di Renzo GC, Conry JA, Blake J, DeFrancesco MS, DeNicola N, Martin JN Jr, et al. 2015. International Federation of Gynecology and Obstetrics opinion on reproductive health impacts of exposure to toxic environmental chemicals. Int J Gynaecol Obstet 131(3):219–225.
Dodson RE, Van den Eede N, Covaci A, Perovich LJ, Brody JG, Rudel RA. 2014. Urinary biomonitoring of phosphate flame retardants: levels in California adults and recommendations for future studies. Environ Sci Technol 48(23):13625–13633.
Farhat A, Buick JK, Williams A, Yauk CL, O'Brien JM, Crump D, et al. 2014. Tris(1,3-dichloro-2-propyl) phosphate perturbs the expression of genes involved in immune response and lipid and steroid metabolism in chicken embryos. Toxicol Appl Pharm 275(2):104–112.
Farhat A, Crump D, Chiu S, Williams KL, Letcher RJ, Gauthier LT, et al. 2013. In Ovo effects of two organophosphate flame retardants–TCPP and TDCPP–on pipping success, development, mRNA expression, and thyroid hormone levels in chicken embryos. Toxicol Sci 134(1):92–102.
Gore AC, Chappell VA, Fenton SE, Flaws JA, Nadal A, Prins GS, et al. 2015. Executive Summary to EDC-2: The Endocrine Society's Second Scientific Statement on Endocrine-Disrupting Chemicals. Endocr Rev 36(6):593–602.
Hammel SC, Hoffman K, Webster TF, Anderson KA, Stapleton HM. 2016. Measuring personal exposure to organophosphate flame retardants using silicone wristbands and hand wipes. Environ Sci Technol 50(8):4483–4491.
Hoffman K, Daniels JL, Stapleton HM. 2014. Urinary metabolites of organophosphate flame retardants and their variability in pregnant women. Environ Int 63:169–172.
Hoffman K, Garantziotis S, Birnbaum LS, Stapleton HM. 2015. Monitoring indoor exposure to organophosphate flame retardants: hand wipes and house dust. Environ Health Perspect 123(2):160–165.
Hornung RW, Reed LD. 1990. Estimation of average concentration in the presence of nondetectable values. Appl Occ Environ Hyg 5(1):46–51.
Liu X, Ji K, Jo A, Moon H-B, Choi K. 2013. Effects of TDCPP or TPP on gene transcriptions and hormones of HPG axis, and their consequences on reproduction in adult zebrafish (Danio rerio). Aquat Toxicol 134-135:104–111.
Meeker JD, Cooper EM, Stapleton HM, Hauser R. 2013a. Urinary metabolites of organophosphate flame retardants: temporal variability and correlations with house dust concentrations. Environ Health Perspect 121(5):580–585.
Meeker JD, Cooper EM, Stapleton HM, Hauser R. 2013b. Exploratory analysis of urinary metabolites of phosphorus-containing flame retardants in relation to markers of male reproductive health. Endocr Disrupt 1(1):e26306.
Mok-Lin E, Ehrlich S, Williams PL, Petrozza J, Wright DL, Calafat AM, et al. 2010. Urinary bisphenol A concentrations and ovarian response among women undergoing IVF. Int J Androl 33(2):385–393.
Pearson MA, Lu C, Schmotzer BJ, Waller LA, Riederer AM. 2009. Evaluation of physiological measures for correcting variation in urinary output: Implications for assessing environmental chemical exposure in children. J Expos Sci Environ Epidemiol 19(3):336–342.
SART (Society for Assisted Reproductive Technologies). 2014. Clinical Summary Report: All SART Member Clinics., [accessed 2 August 2016].
Schindler AE. 2004. First trimester endocrinology: consequences for diagnosis and treatment of pregnancy failure. Gynecol Endocrinol 18(1):51–57.
Scoccia B, Demir H, Kang Y, Fierro MA, Winston NJ. 2012. In vitro fertilization pregnancy rates in levothyroxine-treated women with hypothyroidism compared to women without thyroid dysfunction disorders. Thyroid 22(6):631–636.
Stapleton HM, Klosterhaus S, Eagle S, Fuh J, Meeker JD, Blum A, et al. 2009. Detection of organophosphate flame retardants in furniture foam and U.S. house dust. Environ Sci Technol 43(19):7490–7495.,
Thoma ME, McLain AC, Louis JF, King RB, Trumble AC, Sundaram R, et al. 2013. Prevalence of infertility in the United States as estimated by the current duration approach and a traditional constructed approach. Fert Steril 99(5):1324–1331.
Van den Eede N, Heffernan AL, Aylward LL, Hobson P, Neels H, Mueller JF, et al. 2015. Age as a determinant of phosphate flame retardant exposure of the Australian population and identification of novel urinary PFR metabolites. Environ Int 74:1–8.
van der Veen I, de Boer J. 2012. Phosphorus flame retardants: properties, production, environmental occurrence, toxicity and analysis. Chemosphere 88(10):1119–1153.
Velkeniers B, Van Meerhaeghe A, Poppe K, Unuane D, Tournaye H, Haentjens P. 2013. Levothyroxine treatment and pregnancy outcome in women with subclinical hypothyroidism undergoing assisted reproduction technologies: systematic review and meta-analysis of RCTs. Hum Repro Update 19(3):251–258.
Wang Q, Lam JC, Han J, Wang X, Guo Y, Lam PK, et al. 2015. Developmental exposure to the organophosphorus flame retardant tris(1,3-dichloro-2-propyl) phosphate: estrogenic activity, endocrine disruption and reproductive effects on zebrafish. Aquat Toxicol 160:163–171.
Wang Q, Liang K, Liu J, Yang L, Guo Y, Liu C, et al. 2013. Exposure of zebrafish embryos/larvae to TDCPP alters concentrations of thyroid hormones and transcriptions of genes involved in the hypothalamic-pituitary-thyroid axis. Aquat Toxicol 126:207–213.
WHO (World Health Organization). 1991. Environmental Health Criteria: 111: Triphenyl Phosphate. [accessed 31 November 2014].

Information & Authors


Published In

Environmental Health Perspectives
Volume 125Issue 816 August 2017
PubMed: 28858831


Received: 26 August 2016
Revision received: 10 April 2017
Accepted: 2 June 2017
Published online: 25 August 2017



Courtney C. Carignan
Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA
Lidia Mínguez-Alarcón
Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA
Craig M. Butt
Nicholas School of the Environment, Duke University, Durham, North Carolina, USA
Paige L. Williams
Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA
Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA
John D. Meeker
Department of Environmental Health Sciences, University of Michigan School of Public Health, Ann Arbor, Michigan, USA
Heather M. Stapleton
Nicholas School of the Environment, Duke University, Durham, North Carolina, USA
Thomas L. Toth
Department of Obstetrics and Gynecology, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts, USA
Jennifer B. Ford
Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA
Russ Hauser
Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA
Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA
Department of Obstetrics and Gynecology, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts, USA
for the EARTH Study Team


Address correspondence to C. Carignan, Department of Environmental Health, Building 1, 14th Floor, 665 Huntington Avenue, Boston, MA 02115 USA. Phone: (617) 432-4572. Email: [email protected]

Metrics & Citations


About Article Metrics


Download citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click DOWNLOAD.

Cited by

  • Suspect and non-target screening of chemicals in household cleaning products, and their toxicity assessment, Environmental Engineering Research, 10.4491/eer.2023.123, 29, 1, (230123-0), (2023).
  • Using in silico tools to predict flame retardant metabolites for more informative exposomics-based approaches, Frontiers in Toxicology, 10.3389/ftox.2023.1216802, 5, (2023).
  • Urinary levels of organophosphate flame retardants metabolites in a young population from Southern Taiwan and potential health effects, Frontiers in Endocrinology, 10.3389/fendo.2023.1173449, 14, (2023).
  • Amniogenesis in Human Amniotic Sac Embryoids after Exposures to Organophosphate Flame Retardants, Environmental Health Perspectives, 10.1289/EHP11958, 131, 4, (2023).
  • Parental Exposure to Environmentally Relevant Concentrations of Bisphenol-A Bis(diphenyl phosphate) Impairs Vascular Development in Offspring through DNA/RNA Methylation-Dependent Transmission, Environmental Science & Technology, 10.1021/acs.est.3c03579, 57, 43, (16176-16189), (2023).
  • Organophosphate Flame Retardants in Pregnant Women: Sources, Occurrence, and Potential Risks to Pregnancy Outcomes, Environmental Science & Technology, 10.1021/acs.est.2c06503, 57, 18, (7109-7128), (2023).
  • Integrating Physiologically Based Pharmacokinetic Modeling-Based Forward Dosimetry and in Vitro Bioassays to Improve the Risk Assessment of Organophosphate Esters on Human Health, Environmental Science & Technology, 10.1021/acs.est.2c04576, 57, 4, (1764-1775), (2023).
  • Emerging organophosphate ester resorcinol bis(diphenyl phosphate) exerts estrogenic effects via estrogen receptor pathways, Toxicology, 10.1016/j.tox.2023.153649, 499, (153649), (2023).
  • Chromatin modifiers: A new class of pollutants with potential epigenetic effects revealed by in vitro assays and transcriptomic analyses, Toxicology, 10.1016/j.tox.2022.153413, 484, (153413), (2023).
  • Dietary intake assessment of known and unknown organophosphate esters (OPEs) in foodstuffs via high-resolution mass spectrometry, Science of The Total Environment, 10.1016/j.scitotenv.2022.158452, 854, (158452), (2023).
  • See more

View Options

View options


View PDF

Get Access

Restore your content access

Enter your email address to restore your content access:

Note: This functionality works only for purchases done as a guest. If you already have an account, log in to access the content to which you are entitled.







Copy the content Link

Share on social media