Urinary concentrations of phthalate biomarkers and weight change among postmenopausal women: a prospective cohort study

Background

Some phthalates are endocrine disrupting chemicals used as plasticizers in consumer products, and have been associated with obesity in cross-sectional studies, yet prospective evaluations of weight change are lacking. Our objective was to evaluate associations between phthalate biomarker concentrations and weight and weight change among postmenopausal women.

Methods

We performed cross-sectional (N = 997) and longitudinal analyses (N = 660) among postmenopausal Women’s Health Initiative participants. We measured 13 phthalate metabolites and creatinine in spot urine samples provided at baseline. Participants’ weight and height measured at in-person clinic visits at baseline, year 3, and year 6 were used to calculate body mass index (BMI). We fit multivariable multinomial logistic regression models to explore cross-sectional associations between each phthalate biomarker and baseline BMI category. We evaluated longitudinal associations between each biomarker and weight change using mixed effects linear regression models.

Results

In cross-sectional analyses, urinary concentrations of some biomarkers were positively associated with obesity prevalence (e.g. sum of di (2-ethylhexyl) phthalate metabolites [ΣDEHP] 4th vs 1st quartile OR = 3.29, 95% CI 1.80–6.03 [p trend< 0.001] vs normal). In longitudinal analyses, positive trends with weight gain between baseline and year 3 were observed for mono-(2-ethyl-5-oxohexyl) phthalate, monoethyl phthalate (MEP), mono-hydroxybutyl phthalate, and mono-hydroxyisobutyl phthalate (e.g. + 2.32 kg [95% CI 0.93–3.72] for 4th vs 1st quartile of MEP; p trend < 0.001). No statistically significant associations were observed between biomarkers and weight gain over 6 years.

Conclusions

Certain phthalates may contribute to short-term weight gain among postmenopausal women.

The potential contribution of environmental factors to obesity is of increasing interest. “Obesogens” describes environmental chemicals hypothesized to promote obesity due to altered regulation of adipogenesis and lipid metabolism. Phthalates are endocrine disrupting chemicals present in many consumer products (e.g. cosmetics, food packaging, medications) and are ubiquitous in the environment. Nearly all U.S. residents have detectable concentrations of phthalate metabolites in their urine, though concentrations vary widely [1]. Limited in vitro data suggest that certain phthalates may alter pathways that promote adipogenesis [2, 3], and thus could impact development of obesity.

Scant research evaluating associations between phthalate exposure and body weight report inconsistent findings. Some cross-sectional studies report positive associations between certain phthalate metabolite concentrations and body mass index (BMI) and obesity among adult women. Specifically, one cross-sectional study using data from the 2007–2010 National Health and Nutrition Examination Survey (NHANES) reported increased prevalence of obesity associated with higher concentrations of mono-carboxyoctyl phthalate (MCOP), mono (2-ethyl-5-carboxypentyl) phthalate (MECPP), mono (2-ethyl-5-hydroxyhexyl) phthalate (MEHHP), and the sum of di-ethylhexyl phthalate metabolites (ΣDEHP) [4]. A separate cross-sectional study using 1999–2004 NHANES data observed increased obesity prevalence associated with mono (2-ethylhexyl) phthalate (MEHP) and mono-butyl phthalate (MBP), yet found borderline significant inverse associations between MECPP, MEHHP, mono (2-ethyl-5-oxohexyl) phthalate (MEOHP), ΣDEHP and BMI category [5]. Also, other studies, including one using 1999–2002 NHANES data [6] and another within the Nurses’ Health Study (NHS) and NHS2 cohorts [7], reported inverse cross-sectional associations with MEHP [6], MBP [6, 7], mono-benzyl phthalate (MBzP) [7], and mono-isobutyl phthalate (MiBP) [7].

Phthalates are rapidly metabolized in the body and excreted in urine, and urinary phthalate metabolite concentrations reflect recent exposures [8]. Therefore, the observed cross-sectional associations may reflect confounding via exposure from sources that are themselves associated with obesity, as opposed to causal associations.

One prior prospective analysis, among 977 women aged 32–79 from the Nurses’ Health Study (NHS) and NHS2 [7], reported positive associations with weight gain for MBzP (+ 0.42 kg/year for 4th vs 1st quartiles) and the sum of butyl phthalate metabolites (MBP and MiBP; + 0.34 kg/year for 4th vs 1st quartiles) over a 10 year follow-up period. Weight change was not associated with concentrations of ΣDEHP or mono-ethyl phthalate (MEP) [7].

Whether phthalates affect weight gain remains an unanswered, yet critically important, question. We prospectively evaluated associations between 13 phthalate metabolites (or their sums) and weight change among 997 postmenopausal women enrolled in the Women’s Health Initiative (WHI).

Subjects and methods

Study population

We included 1257 postmenopausal women selected for a nested case-control study of phthalates and breast cancer risk within the WHI. The design of the WHI has been reported previously [9]. Briefly, from October 1, 1993 to December 21, 1998 a total of 161,808 women aged 50–79 years were enrolled in the WHI. WHI participants who were enrolled at three bone density substudy sites (Birmingham, AL; Pittsburgh, PA; Tucson/Phoenix, AZ) provided first morning void urine samples at baseline. A nested case-control study of breast cancer within the WHI quantified urinary concentrations of phthalate metabolite on 419 incident breast cancer cases and 838 matched controls selected from among these bone density substudy participants. Breast cancer cases were selected as all cases of invasive breast carcinoma that occurred among these participants after the year 3 follow-up clinic visit through 2013; controls were matched on enrollment date, length of follow-up, age at enrollment, and WHI study arm with a 1:2 ratio. This analysis includes 997 participants (337 cases, 660 controls) with complete data available (Fig. 1). The longitudinal analysis included only participants selected as controls (N = 660) in the parent study, given that weight gain is common following breast cancer treatment [10].

Selection of study population for the analysis of phthalate biomarkers and weight change in postmenopausal women

Full size image

All participants provided written informed consent upon enrollment into the WHI. The WHI was approved by institutional review boards (IRB) at each clinical center. Additionally IRB approval for the present study was obtained from the University of Massachusetts Amherst. The involvement of the Centers for Disease Control and Prevention (CDC) laboratory in the analysis of samples did not constitute engagement in human subjects research.

Quantification of urinary phthalate metabolites

WHI followed a standard collection, processing, and storage protocol at the three clinical centers that collected urine samples. First morning void urine samples were collected at home and processed within 30 min after participants arrived at the clinic. WHI recommended, but did not require, the use of phthalate-free polypropylene urine collection containers; one site used the recommended containers while the composition of the containers used at the other two clinical centers is unknown. However, all sites used polypropylene centrifuge tubes and cryovials for long-term storage. Additionally, we measured concentrations of metabolites as opposed to the parent phthalate, which should reflect endogenous exposure as opposed to contamination. Urine samples were centrifuged for 5 min at 1330×g and 1.8 mL aliquots were frozen and shipped, packed in dry ice, via overnight FedEx to McKesson Bioservices where they were stored at − 70 °C.

Thirteen phthalate metabolites were measured in baseline urine samples at the CDC: MEP, MBP, mono-hydroxybutyl phthalate (MHBP), MiBP, mono-hydroxyisobutyl phthalate (MHiBP), MBzP, mono (3-carboxypropyl) phthalate (MCPP), MEHP, MEHHP, MEOHP, MECPP, MCOP, and mono-carboxynonyl phthalate (MCNP). Concentrations of phthalate metabolites were quantified after enzymatic hydrolysis of the conjugated metabolites followed by on-line solid phase extraction coupled to high performance liquid chromatography-electrospray ionization-isotope dilution tandem mass spectrometry. Complete details of the analytical method are published online at https://wwwn.cdc.gov/nchs/data/nhanes/2013-2014/labmethods/PHTHTE_H_MET_Phthalates.pdf. The limits of detection (LODs) were in the low ng/mL range. Study samples were randomly distributed through the analytical batches, with cases and matched controls analyzed together. A blinded 10% quality control sample was included, and estimated CVs were as follows: MBP 5.4%, MBzP 6.1%, MCNP 4.7%, MCOP 6.3%, MCPP 5.8%, MECPP 4.3%, MEHHP 5.4%, MEHP 19.5%, MEOHP 6.0%, MEP 3.1%, MHBP 9.0%, MHiBP 21.9%, MiBP 10.3%; the higher average CVs for MEHP and MHiBP reflect small differences in absolute levels of replicates having very low concentrations. All laboratory staff were masked to the identity, disease status, and demographic and risk factor characteristics of the samples. Creatinine was also measured by using an enzymatic assay at CDC on a Roche Modular P Chemistry Analyzer (Indianapolis, IN). The LOD for creatinine was 10 mg/L and the CV of the blinded quality control sample was 2.5%.

We analyzed concentrations of each phthalate metabolite individually. For phthalates with multiple measured metabolites, we also grouped the data by parent phthalate by dividing each metabolite of a single parent by its molecular weight and then summing across metabolites [11, 12]. For example, we calculated the molar sum of DEHP metabolites (ΣDEHP) by dividing each metabolite concentration by its molar mass and then summing the individual concentrations (μmol/L): [MEHHP × (1/294.35)] + MEHP × (1/278.34)] + [MECPP × (1/308.33)] + [MEOHP × (1/292.33)]. The sum of dibutyl phthalate metabolites (ΣDBP) was calculated as the molar sum of MBP and MHBP, and the sum of di-isobutyl phthalate metabolites (ΣDiBP) was calculated as the molar sum of MiBP and MHiBP.

Measurement of weight and BMI calculation

Height and weight were measured at the baseline, year 3, and year 6 clinic visits and used to calculate BMI as weight (kg)/height² (m²) grouped as: underweight/normal weight (< 25.0 kg/m²), overweight (25.0– < 30.0 kg/m²), and obese (≥30.0 kg/m²).

Assessment of covariates

Extensive data on demographic, reproductive, medical history, and behavioral characteristics were collected in the WHI using self-administrated questionnaires at baseline. We considered the following variables as covariates: age (continuous), race/ethnicity (Caucasian, African American, Hispanic/Latino, other), education level (<high school, high school/some college, college degree and higher), income (<$10,000, $10,000-$19,999, $20,000-$34,999, $35,000-$49,999, ≥ $50000), health insurance (no insurance, military insurance, Medicare, Medicaid, private insurance), smoking status (never smoker, past smoker, current smoker), alcohol use (non-drinker, past drinkers, current drinkers), Healthy Eating Index-2005 (HEI-2005, [13]) score (continuous), dietary energy intake (kcal per day; continuous), total recreational physical activity (categorized in quartiles of Metabolic Equivalent values per week (METs/wk.; < 1.25 METs/wk., [1.25- < 6.38 METs/wk., 6.38- < 16.5 METs/wk., ≥16.5 METs/wk), oral contraceptive use (ever, never), any hormone therapy use (never, past, current), ever had diabetes (no, yes), ever had cardiovascular disease (no, yes), hypertension (never hypertensive, untreated hypertensive, treated hypertensive), and dyslipidemia (no, yes).

Statistical analyses

Phthalate biomarker concentrations were natural log transformed to improve normality. Baseline characteristics were summarized according to the BMI categories and differences assessed using analysis of variance (ANOVA) or chi square tests, as appropriate. Geometric means were calculated for each creatinine-standardized phthalate biomarker (i.e. individual metabolite or sum of metabolites of a common parent phthalate) with stratification on baseline BMI group, and differences across groups were assessed with ANOVA.

In cross-sectional analyses we included both cases and controls, given that cases were all diagnosed following the year 3 clinic visit and thus were considered “healthy” at baseline. We categorized phthalate metabolite concentrations into quartiles using the distribution among the controls. Linear regression and multinomial logistic regression analyses were used to model the relationship of each individual phthalate biomarker and baseline weight and BMI category, respectively. All models were adjusted for age and urinary creatinine concentration. We built the regression models by 1) fitting univariable linear and multinominal logistic regression models for each variable with weight and BMI, respectively, 2) including all variables with p < 0.25 in the univariable model in a preliminary multivariable model along with the phthalate biomarker, and 3) evaluating the significance of each covariate using backward selection and retaining all covariates with a p value < 0.10 or of known biological importance. A common set of covariates was included in the multivariable models to facilitate comparisons across phthalate biomarkers. Trends in the weight β coefficient and the odds ratio (OR) of overweight and obesity with increasing categories of phthalate biomarker were evaluated by testing the significance of a continuous variable including the median concentration of each biomarker quartile in the regression model. We included 997 participants with complete data on covariates, exposure, and outcomes in our analysis.

We modeled the prospective weight change rate over 3 and 6 years by the quartiles of urinary phthalate biomarker concentrations using mixed-effect models with: a random coefficient, a fixed effect for weight and year of follow up, and including product terms between phthalate biomarkers and year of follow up (i.e. year 3 and year 6). A parsimonious multivariable model was built using the process described above. Analyses were repeated with stratification on baseline BMI to evaluate possible effect modification, and we plotted predicted weight change over time by BMI category for models including an interaction with BMI and a model without this term. We obtained p-values for linear trends by including an interaction term between each year of follow up and the median concentration of each biomarker quartile in the mixed-models as a continuous variable. We considered a P-value < 0.05 as statistically significant. All analyses were conducted using SAS version 9.4 (SAS Institute, Cary, North Carolina) and Stata version 15.0 (Stata Corp, College Station, TX).

Table 1 shows baseline characteristics of the study population by BMI category. Compared to underweight/normal weight women, obese women were older, and more likely to be Black/African American, earn <$20,000/year, and have lower educational attainment. Obese women also were less likely to drink alcohol, had lower diet quality and higher dietary energy intake, had lower physical activity, and had more hypertension and diabetes than underweight/normal weight women.

Table 2 displays the creatinine-corrected geometric mean concentrations and 95% confidence intervals for each phthalate metabolite by BMI category. MECPP and MEHHP concentrations were significantly higher among obese and overweight women compared to underweight/normal weight women, although the difference in means was small. MiBP concentrations were observed to be slightly higher among obese and overweight women compared to underweight/normal weight women. The geometric means of the other measured metabolites did not significantly differ by BMI category.

Table 3 displays cross-sectional associations between urinary phthalate biomarker concentrations and BMI category. In multivariable adjusted multinomial logistic regression models, we observed statistically significant positive trends in association between quartiles of MiBP, MCNP, MCOP, MCPP, ΣDEHP, and the DEHP metabolites MECPP, MEHHP, and MEOHP. The association between ΣDEHP and BMI category was particularly strong, with those in the 4th quartile of ΣDEHP having two-fold increased odds of overweight (OR 2.72, 95% CI 1.57–4.72) and three-fold increased odds of obesity (OR 3.29, 95% CI 1.80–6.03) compared to those in the 1st quartile. Obesity was positively associated with MBzP (OR 2.73, 95% CI 1.48–5.04 for 4th vs 1st quartile; p trend = 0.01), though the association of overweight was attenuated (OR 1.58, 95% CI 0.94–2.66; p trend = 0.18). MEP and MHBP concentrations were inversely associated with obesity, although the association with MHBP was not statistically significant (p trend = 0.04 and p trend = 0.10, respectively). Analyses with weight as the outcome showed similar trends, as expected (Additional file 1: Table S1). We observed similar associations when restricting to participants selected as controls (Additional file 1: Tables S2 and Table S3). We reran analyses excluding the four variables with the highest amounts of missing data as covariates in the model (physical activity, hypertension, high cholesterol, and cardiovascular disease history); results were generally similar for all phthalate biomarkers in this larger sample population (N = 1187) (Additional file 1: Table S4).

Table 4 presents the estimated additional weight change (i.e. beyond the average annual weight change in the study population) associated with phthalate biomarker concentrations after three and 6 years of follow up. We observed positive associations between some phthalate biomarker concentrations and weight change over 3 years (β for 4th vs 1st quartile): MEP (β = 2.32, 95% CI 0.93–3.72; p trend = 0.001), MBP (β = 1.24, 95% CI -0.22-2.50; p trend = 0.05), MHiBP (β = 1.98, 95% CI 0.62–3.33; p trend = 0.02), and MEOHP (β = 1.44, 95% CI 0.07–2.80; p trend = 0.05). These coefficients are interpretable as additional weight gain in the 4th versus 1st quartile of exposure; for example, women in the 4th quartile of MEP gained an additional 2.32 kg over 3 years. Significant associations were observed only in the 4th quartile of each biomarker compared to the 1st quartile. Phthalate biomarker concentrations at baseline were not associated with weight change over 6 years of follow-up (all p trend > 0.10). We did not observe effect modification by baseline BMI for any of the phthalate biomarkers evaluated.

In this sample of postmenopausal women from the WHI, we report that concentrations of certain phthalate biomarkers are associated with increased BMI and weight cross-sectionally and also with weight gain over 3 years. Notably, in cross-sectional analyses, women in the highest quartile of ΣDEHP were two times more likely to be overweight (OR 2.72, 95% CI 1.57–4.72) and three times more likely to be obese (OR 3.29, 95% CI 1.80–6.03) compared to those in the lowest quartile. However, in prospective analyses ΣDEHP was not associated with weight gain, and only one DEHP metabolite, MEOHP, was associated: women in the highest quartile of MEOHP concentration gained an additional 1.44 kg over 3 years compared to women with the lowest MEOHP concentrations. Although we observed positive cross-sectional associations between MiBP, MBzP, MCNP, MCOP, MCPP, MECPP, MEHHP, and MEHP and BMI, concentrations of these phthalate biomarkers were not associated with weight change. This suggests that higher concentrations of these phthalate biomarkers might be related to increased exposure to the parent phthalates among the obese, as opposed to a causal relationship with weight gain. Interestingly, concentrations of MEP were inversely associated with obesity in cross-sectional analyses, yet displayed significant positive associations with weight gain over 3 years (4th vs 1st quartile: + 2.32 kg, 95% CI 0.93–3.72; p trend = 0.001). Diethyl phthalate (DEP), the parent compound of MEP, is primarily found in personal care products and cosmetics [14, 15]. One potential explanation of our findings is that use of such products varies by BMI category, thus confounding the cross-sectional associations between MEP and BMI, yet exposure to DEP via these products may promote weight gain. Future studies will be useful in elucidating the potential effect of DEP exposure on body weight.

The few prior cross-sectional evaluations of urinary phthalate metabolites and BMI in adult women produced inconsistent results. Overall, our findings are generally consistent with those of Buser et al [4] (using 2007–2010 NHANES data) but are in contrast with the other studies. Specifically, we and Buser et al [4] reported positive associations with MECPP, MEHHP, MEOHP, and DEHP, while Yaghjyan et al [5] (using 1999–2004 NHANES data) reported inverse associations. Yaghjyan et al [5] also reported a positive association with MEHP, which agrees with our findings; however, an early study by Hatch et al [6] (using 1999–2002 NHANES data) reported an inverse association between MEHP and BMI and Buser et al [4] observed no association. Our findings support the positive association between MCOP and BMI reported by Buser et al [4], although they reported no association with MCNP while we observed a positive association. The inverse association we observed with MEP contrasts with the non-significant positive trends reported by Hatch et al [6] and Buser et al [4], although no association was observed by Yaghjyan et al [5] or Song et al [7] (using a sub-sample of NHS/NHS2 participants). We report no significant association between MBP and BMI, in agreement with Buser et al [4], while an inverse association was observed by Hatch et al [6] and Song et al [7], and Yaghjyan et al [5] observed a positive association. Finally, the positive association we observed between MBzP and BMI is in contrast to the inverse association reported by Song et al [7] and the null findings of Hatch et al [6], Buser et al [4], and Yaghjyan et al [5].

The reasons for these differences across studies are unclear, but may relate to differences in population exposure to phthalates over time. This is unlikely to fully explain observed differences, however, given that our results most closely align with the Buser et al [4] study, which used samples donated more than a decade following ours. Additionally, the DEHP metabolites measured in NHANES varied over time, with only MEHP measured in 1999–2000, MEOHP and MEHHP first measured in 2001–2002, and MECPP first measured in 2003–2004. Also, differences in the demographic characteristics of our highly selected population compared to those of NHANES participants, which are randomly selected to be representative of the general U.S. population, and the older age of our participants as compared to NHS/NHS2 participants may have affected results. Another important factor contributing to the inconsistency across studies may be the relatively high within-person variability in concentrations of several urinary phthalate metabolites [16], which may result in substantial non-differential misclassification. In sum, our results and those of Buser et al [4] support positive cross-sectional associations between some phthalate metabolites, especially the DEHP metabolites, and BMI category.

Importantly, the observed positive prospective associations between some phthalate biomarkers and weight change over 3 years did not persist over 6 years. Since the half-lives of phthalate metabolites are short (6–12 h) and urinary concentrations exhibit within-person variability over time [16], a single measurement of phthalate metabolites taken at baseline is unlikely to fully characterize an individual’s long-term exposure. Thus, a single measurement of a given phthalate exposure biomarker may predict weight gain over short time periods, but not over longer time periods.

Only one prior study prospectively evaluated associations between phthalate exposure and weight change. In a sample of NHS and NHS2 participants, Song and colleagues [7] reported significantly higher annual weight gain associated with MBzP (β = 0.42; 95% CI: 0.26–0.57, p-trend< 0.001), phthalic acid (a non-specific biomarker of phthalate exposure) (β = 0.33; 95% CI: 0.15–0.50, p-trend = 0.001), the sum of butyl phthalate metabolites (MiBP and MBP) (β = 0.34; 95% CI: 0.18–0.50, p-trend < 0.001) and the total sum of all the phthalate metabolites measured (β = 0.17; 95% CI: 0.02–0.33, p-trend = 0.05). Our observed 1.24 kg increased weight gain over 3 years among associated with the highest quartile of MBP is consistent with their reported 0.34 kg annual weight gain, which would be equivalent to a gain of 1.02 kg over 3 years. We did not observe statistically significantly increased weight gain associated with MBzP, however. We also observed increased weight gain associated with higher concentrations of MEOHP, MEP, and MHiBP; MEP was not associated with weight gain in the Song et al [7] study, and MEOHP and MHiBP were not separately evaluated. Our study population was, on average, older and of higher baseline weight than the NHS/NHS2 study population, which could explain some differences between our findings, given that patterns of weight gain and loss vary among older age groups. Also, Song et al [7] utilized self-reported weight over a 10-year period as their outcome, while we utilized measured weight over 3- and 6-year intervals. Future work will benefit from utilizing repeated measurements of phthalate biomarkers over extended time periods.

It is noteworthy that we observed some inconsistent associations between the cross-sectional and the prospective analyses. Specifically, we observed that MEP was inversely associated with BMI in the cross-sectional analysis, yet positively associated with weight gain in the prospective analysis. Also, MCNP, MCOP, and MECPP, were positively associated with BMI in the cross-sectional analysis yet were not associated with prospective weight gain. These differences may be explained in part by differences in the study population used for the cross-sectional and prospective analyses; we excluded women later diagnosed with breast cancer from the prospective analysis, which resulted in a significantly reduced sample size (997 vs 660) and loss of statistical power. Cases and controls also may differ in ways that could have affected results, although we observed similar cross-sectional results among controls. Another potential explanation is that phthalate biomarker concentrations may reflect exposure to a source of phthalates (e.g. processed foods) that is also associated with obesity, and thus the cross-sectional associations we observed could be due to reverse causality. Overweight/obese woman may have demographic characteristics (e.g. non-White, lower socioeconomic status) and/or engage in health behaviors (e.g. alcohol use, poor diet quality, consumption of canned/processed foods) that are associated with higher phthalate exposure. Additionally, overweight and obese women are more likely to have comorbidities that result in more frequent contact with medical equipment and medications, which also may cause increased exposure to certain phthalates. While we adjusted for many factors including race/ethnicity, dietary quality and energy intake, and comorbidities, it is possible that residual confounding affected our results.

Laboratory evidence suggests several mechanisms by which phthalates could act as obesogens. Some phthalates may induce adipogenesis by activating metabolic sensors, such as the peroxisome proliferator-activated receptors (PPARs) [2]. PPARs are transcription factors that play an essential role in energy metabolism. PPARγ, expressed in the endothelial, vascular and smooth muscle cells, controls fatty acids storage in the adipose tissue by promoting adipocyte differentiation [17]. Evidence from experimental studies has shown that MBzP, mono-sec-butyl phthalate, and MEHP [2, 3], directly activate PPARγ and stimulate a dose dependent increase in adipocyte differentiation [3], which in turn promotes adipogenesis. Furthermore, MEP has been shown to increase activity of human PPARγ [18]. Phthalates also are known to disrupt the thyroid hormone system, which also could affect weight gain. Triiodothyronine (T3), a thyroid hormone, regulates adipogenesis by controlling genes that are involved in lipolysis and lipogenesis [19, 20]. In a recent study DBP and DEHP were shown to inhibit the regular activity of T3, with DBP acting as a stronger antagonist than DEHP [21]. Such findings are consistent with our reported positive associations between some DBP and DEHP metabolites (MBP, MEOHP) and weight gain.

Our study has several important strengths. We utilized objectively measured data on height and weight, which reduces misclassification in the outcome measures. We also considered a broad range of variables as potential confounders using the extensive WHI data resource. We considered a broad panel of 13 phthalate metabolites, which allowed for a more complete evaluation of phthalate exposure in our population than in prior studies.

Our findings should be interpreted in light of additional relevant limitations. We used a single measure of urinary phthalate biomarkers to estimate the study participants’ long-term exposure to phthalates. Phthalates are rapidly metabolized and excreted within hours, and a single measurement may not accurately reflect the long-term exposure levels of our participants. Furthermore, the within person variability of the phthalate metabolites could lead to non-differential misclassification of the exposure and attenuate the true associations; therefore, we may have failed to detect some true associations and others may be even stronger than we have reported. Additionally, urinary creatinine is an imperfect correction for urine dilution and could be associated with BMI and/or factors predictive of BMI [22]. Type I error also is possible given the number of statistical tests performed. Also, under-reporting of energy intake by overweight/obese participants has been reported in WHI [23, 24], which could affect the HEI-2005 scores and cause residual confounding, though we expect this to have minimal effects on our estimates because diet is not the major source of phthalate exposure. Some participants were excluded from the full analyses because of missing data on covariates. While multiple imputation is useful when missing data are present, we were concerned that the resulting increased confidence interval width would further obscure any potential associations given the measurement error that already affected our analyses. In sensitivity analyses of the cross-sectional study which included 1187 participants instead of the 997 reported in the Results section, results were generally similar for all phthalate biomarkers, thus we anticipate that any potential bias from our complete case analysis was minimal. Lastly, we utilized a highly selected population of postmenopausal women, with limited racial/ethnic diversity. Future work in populations including men, younger women, and increased minority representation will be useful to fully understand whether phthalates can affect weight gain in other populations.

We provide evidence that certain phthalates may contribute to weight gain among postmenopausal women at 3 years of follow-up, but not at 6 years of follow-up. Additional research to characterize the potential role of phthalates in development of obesity is important given the continued obesity epidemic and the ubiquity of these chemicals.

Program Office: Jacques Rossouw, Shari Ludlam, Joan McGowan, Leslie Ford, and Nancy Geller; National Heart, Lung, and Blood Institute, Bethesda, MD.

Clinical Coordinating Center: Garnet Anderson, Ross Prentice, Andrea LaCroix, and Charles Kooperberg; Fred Hutchinson Cancer Research Center, Seattle, WA.

Investigators and Academic Centers: JoAnn E. Manson, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA; Barbara V. Howard, MedStar Health Research Institute/Howard University, Washington, DC; Marcia L. Stefanick, Stanford Prevention Research Center, Stanford, CA; Rebecca Jackson, The Ohio State University, Columbus, OH; Cynthia A. Thomson, University of Arizona, Tucson/Phoenix, AZ; Jean Wactawski-Wende, University at Buffalo, Buffalo, NY; Marian Limacher, University of Florida, Gainesville/Jacksonville, FL; Jennifer Robinson, University of Iowa, Iowa City/Davenport, IA; Lewis Kuller, University of Pittsburgh, Pittsburgh, PA; Sally Shumaker, Wake Forest University School of Medicine, Winston-Salem, NC; Robert Brunner, University of Nevada, Reno, NV; Karen L. Margolis, University of Minnesota, Minneapolis, MN.

Centers for Disease Control and Prevention, Atlanta, GA: Xiaoyun Ye, Manori J. Silva, Ella Samandar, James Preau, Jr., and Tao Jia.

Funding

This work was supported by the National Institute of Environmental Health Sciences (R01ES024731, R01ES024731S1). The Women’s Health Initiative is funded by the National Heart, Lung, and Blood Institute, National Institutes of Health, U.S. Department of Health and Human Services through contracts HHSN268201600018C, HHSN268201600001C, HHSN268201600002C, HHSN268201600003C, and HHSN268201600004C.

Availability of data and materials

The data that support the findings of this study are available from the Women’s Health Initiative Coordinating Center, but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of the Women’s Health Initiative Coordinating Center.

Disclaimer

The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention (CDC). Use of trade names is for identification only and does not imply endorsement by the CDC, the Public Health Service, or the US Department of Health and Human Services.

Affiliations

1. Department of Biostatistics and Epidemiology, University of Massachusetts Amherst, 411 Arnold House, 715 North Pleasant Street, Amherst, MA, 01003, USA

   • Mary V. Díaz Santana
   • , Susan E. Hankinson
   • , Carol Bigelow
   • , Susan R. Sturgeon
   •  & Katherine W. Reeves

2. Department of Biology, University of Massachusetts Amherst, Amherst, MA, USA

   • R. Thomas Zoeller

3. Cancer Prevention Program, Fred Hutchinson Cancer Research Center, Seattle, WA, USA

   • Lesley Tinker

4. Department of Medicine, Brigham and Women’s Hospital, Harvard Medical School, and the Harvard T.H. Chan School of Public Health, Boston, MA, USA

   • Jo Ann E. Manson

5. Division of Laboratory Sciences, National Center for Environmental Health, Centers for Disease Control and Prevention, Atlanta, GA, USA

   • Antonia M. Calafat

6. Program in Public Health, Department of Family Population and Preventive Medicine, Stony Brook University, Stony Brook, USA

   • Jaymie R. Meliker

Contributions

KR, SH, CB, JM, RTZ, and LT conceptualized the study and obtained funding to conduct the research. AC performed quantification of phthalate metabolite levels in the participant urine samples. MDS, KR, and CB performed the statistical analyses. MDS and KR had primary responsibility for writing the manuscript. All authors contributed to the writing of the manuscript and read and approved the final manuscript.

Corresponding author

Correspondence to Katherine W. Reeves.

Ethics approval and consent to participate

All participants provided written informed consent upon enrollment into the WHI. The WHI was approved by institutional review boards (IRB) at each clinical center. Additionally IRB approval for the present study was obtained from the University of Massachusetts Amherst. The involvement of the Centers for Disease Control and Prevention (CDC) laboratory in the analysis of samples did not constitute engagement in human subjects research.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Reprints and Permissions