- Research article
- Open Access
Urinary phthalate exposures and risk of breast cancer: the Multiethnic Cohort study
Breast Cancer Research volume 23, Article number: 44 (2021)
The epidemiologic evidence from observational studies on breast cancer risk and phthalates, endocrine disrupting chemicals, has been inconsistent. In the only previous study based on pre-diagnostic urinary phthalates and risk of breast cancer, results were null in mostly white women.
We examined the association between pre-diagnostic urinary phthalates and breast cancer in a nested case-control study within the Multiethnic Cohort (MEC) study, presenting the first data from five major racial/ethnic groups in the USA. We measured 10 phthalate metabolites and phthalic acid, using a sensitive liquid chromatography mass spectrometry assay on 1032 women with breast cancer (48 African Americans, 77 Latinos, 155 Native Hawaiians, 478 Japanese Americans, and 274 Whites) and 1030 matched controls. Conditional logistic regression was used to examine risk with individual metabolites and ratios of primary (MEHP, mono-2-ethylhexyl-phthalate) to secondary (MEHHP, mono(2-ethyl-5-hydroxyhexyl); MEOHP, mono(2-ethyl-5-oxohexy)) metabolites of di-2-ethylhexyl phthalate (DEHP), a widely used plasticizer. In addition, we investigated risk associations with high (∑HMWP) and low molecular weight (∑LMWP) phthalates, as well as total phthalates which included high and low molecular weight phthalates with phthalic acid (∑LMHMPA) or without phthalic acid in molar ratios (∑LMHMmolar) and adjusted for creatinine and potential confounders.
Among all women, breast cancer risk was higher for those in tertile 2 and tertile 3 of primary to secondary metabolites of DEHP (MEHP/(MEHHP + MEOHP)) in comparison to those in tertile 1; the respective odds ratios were 1.32 (95% CI 1.04–1.68) and 1.26 (95% CI 0.96–1.66) (Ptrend = 0.05). Risk among Native Hawaiian women increased with exposures to eight of ten individual phthalates and total phthalates (∑LMHMPA ORT3 vs T1 = 2.66, 95% CI 1.39–5.12, Ptrend = 0.001). In analysis by hormone receptor (HR) status, exposure above the median of ∑LMWP was associated with an increased risk of HR-positive breast cancer (OR = 1.30, 95% CI 1.05–1.60) while above the median exposure to phthalic acid was associated with an increased risk of HR-negative breast cancer (ORabove vs below median = 1.59, 95% CI 1.01–2.48).
Further investigations of suggestive associations of elevated breast cancer risk with higher ratios of primary to secondary metabolites of DEHP, and differences in risk patterns by race/ethnicity and HR status are warranted.
Phthalates are industrial chemicals that are present in numerous consumer products and solvents, as additives, and plasticizers [1,2,3], and have known endocrine disrupting properties . Results from the Third National Health and Nutrition Examination Survey (NHANES III) showed ubiquitous exposure of the US population to a wide variety of phthalates . Mono-ethyl phthalate (MEP), mono-n-butyl phthalate (MBP), and mono-benzyl phthalate (MBzP), metabolites of phthalates, frequently found in personal care products, were detectable in over 97% of urine samples while mono-2-ethylhexyl-phthalate (MEHP), a metabolite of di-2-ethylhexyl phthalate (DEHP), found commonly in plasticizers, food packaging, and household products, was detectable in over 75% of the samples. These high exposures were observed across all ages, in whites and nonwhites, although there were age, sex, and racial/ethnic differences that likely reflected variations in exposure patterns  and differences in the metabolism of phthalates . Phthalates are typically metabolized by phase I hydrolysis to their respective monoesters, followed by phase II conjugation, which depending on the specific phthalates can be further metabolized via oxidation to secondary metabolites .
Phthalates have been implicated to influence developmental and reproductive processes and exert carcinogenic effects [1, 8,9,10]. However, epidemiological evidence on the association between phthalate exposure and breast cancer risk remains inconsistent. The first two case-control studies reported significantly elevated risk in relation to exposure to MEP, mono(2ethyl-5-carboxy-pentyl) phthalate (MECPP)  and MEHP , but they were small studies (75 to 233 breast cancers), and measurements were conducted in post-diagnosis samples. Urinary phthalate exposures were unrelated to risk in two larger studies (~ 400–700 breast cancers) [13, 14]; only one investigated exposures before diagnosis . These studies differed in study design, sample size, inclusion of in situ and invasive breast cancers [13, 14], and race/ethnicity composition. Positive associations were reported in studies conducted in Northern Mexico  and Alaska (mostly Eskimos, Indians, and Aleut)  whereas null results were found in studies of mainly whites from the Women’s Health Initiative (WHI)  and the Long Island Breast Cancer Study Project (LIBCSP) . Breast cancer risk was also associated with di-methyl-phthalate (DMP) exposure in a Danish prospective registry study which quantified phthalate exposure using prescription files .
To further investigate the role of pre-diagnostic urinary phthalates and breast cancer risk, we conducted a nested case-control in the Multiethnic Cohort (MEC) study  that included 274 whites and 758 nonwhites (478 Japanese Americans, 155 Native Hawaiians, 77 Latinos, and 48 African Americans) diagnosed with incident breast cancer, and individually matched control women.
Materials and methods
The MEC is a large prospective cohort that included 96,810 men and 118,441 women aged 45–75 years from five different racial/ethnic groups (African Americans, Latinos, Native Hawaiians, Japanese Americans, and whites) living in Hawaii and California (primarily from Los Angeles County) at enrollment between 1993 and 1996 . Participants completed a baseline questionnaire which assessed demographics, lifestyle, diet, and anthropometrics, and for women, menstrual and reproductive histories and hormone therapy use. Participants were followed prospectively for diagnosis of incident breast cancer (invasive and in situ) through routine yearly linkage with the California and Hawaii statewide cancer registries and for vital status through yearly linkages to the National Death Index and state death certificate files which was through 2014 for this nested case-control study. Stage at diagnosis and estrogen/progesterone receptor (ER/PR) status were obtained from cancer registries.
In 2001–2006, a prospective biorepository was established by collecting pre-diagnostic urine and blood specimens from 67,594 MEC cohort members . A short questionnaire was administered at biospecimen collection which assessed weight, use of hormone therapy and medications, and other factors. For this nested case-control study, cases (n = 1032) were diagnosed with incident breast cancer from 2001 through 2014 after urine collection; 22% were in situ and 78% were invasive cancers. The mean time between urine collection and breast cancer diagnosis was 5.5 years (standard deviation (SD) = 3.3). For each case, we selected one control, who was alive and free of breast cancer at the age of breast cancer diagnosis, and individually matched controls to cases on area (Hawaii or California), birth year (± 1 year), race and ethnicity (white, Native Hawaiian, African American, Latino, Japanese American), urine type (first morning from California, overnight and first morning from Hawaii), date (± 1 year) of urine collection, hours of fasting (8–10, > 10 h), and time of blood draw (± 2 h). As described previously, the overnight urine collection started between 5:00 pm and 9:00 pm (depending on participant) and included all urine passed during the night and the first morning urine sample covering a period of 12 h . Controls were sampled from the representative pool of subjects with existing data on obesity-related and inflammation biomarkers and genotype array data. In total, 1030 controls were identified and each control was individually matched to one case except that one white and one Native Hawaiian control were each matched to the two cases of the same race/ethnic group. There were 30 controls (19 whites, 7 Japanese Americans, 2 Native Hawaiians, 2 Latinos, 0 African Americans) and 44 cases (14 whites, 16 Japanese Americans, 6 Native Hawaiians, 4 Latinos, 4 African Americans) who had other cancers with a mean of 6.3 years (SD 5.2) and 8.4 years (SD 5.8), before the donation of urine collection, respectively. Results were unchanged without these individuals (data not shown) and they were included in the analyses.
Laboratory measurement of urinary phthalate metabolites
Phthalate metabolites (MEP, MMP, MBP, MiBP, MBzP, MCHP, MEHP, MEHHP, MEOHP, MECPP, MCMHP) and phthalic acid measurements were conducted at the University of Hawaii Cancer Center Analytical Biochemistry Shared Resource. Dr. Adrian Franke supervised the day-to-day activities and quality assurance and quality control of phthalate measurements using state-of-the-art sensitive isotope-dilution orbitrap-based high-resolution accurate-mass liquid chromatography mass spectrometry (LCMS) assay after enzymatic hydrolysis and liquid-liquid extraction [19, 20]. In brief, 0.1 mL urine was mixed with 0.01 mL of a mixture of isotopically labeled analyte that was used as internal standards followed by incubation with a glucuronidase/sulfatase mixture at 37 °C for 90 min, extraction with methyl tertiary butyl ether, and LCMS analysis. Our lower limits of quantitation were 0.5 ng/mL for each analyte. Individual phthalate all share phthalic acid as a common metabolite as all the phthalate diesters can be metabolized to phthalic acid [21, 22]. One metabolite (mono 2-carboxy-hexyl phthalate, MCMHP) was not measured reliably and was excluded from data analysis. Personnel were blinded to case-control status, and matched pairs of cases and controls were assayed in the same batch. Replicate samples of pooled urines (5%) were included in each of 37 batches for quality control measures and coefficients of variation (CV) were calculated. The CV% (SD/mean concentration × 100) within-batch was 26.7% for phthalic acid and was < 25% for eight of the ten metabolites with a median of 22.6% for the individual metabolites (range was 5.9% (MCHP) to 35.0% (MMP)). The larger CVs likely reflect several samples close to the lower limit of detection (LLOD). The mean CV of the non-blinded pool samples was 11.7% (SD 6.9%). All analytes were adjusted to urinary creatinine [19, 20] and are shown as micrograms per gram (μg/g) of creatinine. Analytes below the LLOD were assigned a value half of the LLOD. Eight of the ten metabolites and phthalic acid were detected at levels above the LLOD in over 92% of women; this was lower for MMP (80%) and MCHP (89%) (Supplementary Table 1). We also present the geometric mean concentration for each phthalate metabolite collected in the overnight urine samples (i.e., all were from Hawaii) and first morning urine samples (i.e., all but 16 were from Los Angeles County) (Supplementary Table 1).
We conducted conditional logistic regression, with the matched sets as strata (1028 pairs and 2 triplets), and modeled phthalate variables as tertiles using selected cutoff points based on the distribution among all controls. Odds ratios (ORtertiles) and 95% confidence intervals (CIs) were the primary statistics of interest, and inference was based on the Wald test. We found no evidence of a nonlinear relationship (on the log odds scale) between phthalates and risk using restricted cubic splines (data not shown). Therefore, log-transformed phthalate variables were used as trend variables to test for dose-response relationships. Models were adjusted for potential confounders that were not matching factors (e.g., established breast cancer risk factors) via indicator variables for tertiles of propensity scores for exposure to phthalates [23, 24] in order to maximize power. Specifically, an ordinal logistic regression for tertiles of each phthalate variable was performed using the following independent variables: age at urine collection, education, number of children, age at menarche, menopausal status, body mass index (BMI) at urine collection, neighborhood socioeconomic status (nSES)  at urine collection, smoking, alcohol intake, and Mediterranean diet energy adjusted total score . The nSES index is a composite measure created by principal component analysis (PCA) of US Census data that incorporated census block group data on education, occupation, unemployment, household income, poverty, rent, and house values . PCAs were conducted separately for California and Hawaii. Results across states were similar for Eigen vectors and variance explained with a single component identified. This nSES measure was categorized into quintiles based on the nSES distribution of Los Angeles County and Hawaii block groups for California and Hawaii MEC participants, respectively. The propensity for exposure was determined for each individual as the weighted average = 1 × ρ1 + 2 × ρ2 + 2 × ρ3, where ρi is the model-based probability of exposure to tertile i. Heterogeneity of the associations by race/ethnicity was assessed by a global test of the interaction terms between race and the phthalate trend variable. We repeated subgroup analyses for hormone receptor-positive (HR+, ER+, or PR+) and hormone receptor-negative (HR−, ER−, and PR−) cancer, as well as BMI and use of hormone therapy at urine collection. In addition, waist-hip ratio (WHR) was obtained in a follow-up questionnaire and was available for 333 cases prior to diagnosis and 946 control women; we explored risk association results by WHR. We also conducted a sensitivity analysis by restricting to invasive cases only (n = 798) and by lag time between time of urine collection and breast cancer diagnosis (≤ 5 years versus > 5 years), as well as excluding 187 cases that were diagnosed within 2 years of urine collection to minimize the potential effect of pre-diagnostic breast cancer on phthalate levels. Associations with P < 0.05 were considered statistically significant and with 0.05 ≤ P < 0.10 were considered suggestive. The correlations among phthalate metabolites and phthalic acid were examined using Spearman’s rho.
The short-branched phthalates, DMP and DEP, are excreted in urine as the unconjugated monoesters, MMP and MEP, respectively. The longer-branched phthalates such as DEHP are hydrolyzed first to MEHP and subsequently metabolized to MEHHP, MEOHP, MECPP, and other oxidative metabolites [7, 27, 28] (Fig. 1). Oxidation of MEHP to other secondary metabolites effectively decreases the internal body burden of MEHP to presumably less toxic secondary metabolites [29, 30]. Thus, subjects with higher MEHP% and higher ratios of MEHP to major secondary DEHP metabolites (MEHHP, MEOHP, MECPP) may be at higher risk (MEHP% was calculated by converting DEHP metabolites into nanomoles (nmol) using their respective molecular weight and dividing the molar mass of MEHP by the mass of the sum of all four metabolites and then multiplying by 100). For all women, and separately in African Americans and Latinos combined, Native Hawaiians, whites, and Japanese Americans, we examined risk in relation to the 10 metabolites and phthalic acid, MEHP%, the ratios of MEHP to secondary DEHP metabolites, and summary variables, including the sum of all major DEHP metabolites (∑DEHP), low molecular weight phthalates (∑LMWP: MBP, MiBP, MEP, MMP), high molecular weight phthalates (∑HMWP: MBzP, MCHP, ∑DEHP), and total phthalate represented by the sum of all 10 phthalate metabolites and phthalic acid (∑LMHMPA) and the sum based on molar ratios of phthalate metabolites without phthalic acid (∑LMHMmolar). We divided the concentration of each metabolite by its molecular weight to obtain the molar equivalent (micromoles/liter, μmol/L) and then summed the concentrations to get total μm/L of metabolites per creatinine unit.
The majority of participants (86%) were from Hawaii, who donated overnight (875 cases, 873 controls) or first morning urines (8 cases, 8 controls), whereas only first morning samples (149 cases, 149 controls) were collected in Los Angeles County (Table 1). Cases compared to control women were more likely to be nulliparous (15.8% vs 10.8%) and had higher BMI at urine collection (27.0 ± 5.6 vs 26.0 ± 5.6 kg/m2), but were otherwise comparable.
The individual phthalate metabolites were modestly correlated (rho’s were 0.2 to 0.4; P < 0.05) whereas the DEHP metabolites were more strongly correlated (rho’s were 0.5 to ~ 0.8). The correlations between phthalic acid and the 10 phthalates ranged from 0.15 to 0.44. (Supplementary Table 2).
The risk of breast cancer in all women combined tended to decrease with increasing exposure to MBzP (Ptrend = 0.03). In contrast, risk in all women was positively associated with high ratios of MEHP/MEOHP (Ptrend = 0.04), MEHP/(MEOHP+ MEHHP) (Ptrend = 0.052), and high MEHP% (Ptrend = 0.092); the range of ORs associated with exposure in the upper tertile was 1.18 to 1.26 compared to those in the lowest tertile of exposure. These patterns were apparent among white women (Table 2). Although there were no statistically significant racial/ethnic differences in results, some suggestive differences emerged. Among Native Hawaiian women, risks increased in association with eight of the ten individual phthalates, including MiBP (Ptrend = 0.05), MEHP (Ptrend = 0.01), MEHHP (Ptrend = 0.09), MECPP (Ptrend = 0.06), and phthalic acid (Ptrend = 0.03); ORs ranged from 1.94 to 2.73 for Native Hawaiians in the upper tertile of phthalate exposures compared to those in the lower tertile. In contrast, the individual metabolites were not associated with risk in African Americans and Latinos combined or in Japanese Americans. Among whites, below unity ORs were found for MBzP (Ptrend = 0.05) and for all three DEHP secondary metabolites, with strongest inverse trends for MEOHP (Ptrend = 0.06).
Breast cancer risk in all women combined was not associated with ∑DEHP and ∑HMWP but exposure to ∑LMWP, ∑LMHMPA, and ∑LMHMmolar was positively associated with risk in all women combined, with 18 to 23% higher risk for those in the upper tertile of exposure (Table 3). The positive associations were more prominent among Native Hawaiians who displayed statistically significant or suggestive elevated risks with four of the five summation exposures: ∑LMHMPA (Ptrend = 0.001), ∑LMHMmolar (Ptrend = 0.03), ∑LMWP (Ptrend = 0.10), and ∑HMWP (Ptrend = 0.10). Native Hawaiian women in the upper two tertiles of ∑LMHMPA exposure showed a significant 2.4- to 2.6-fold higher risk than their counterparts in the lowest tertile of exposure.
Results did not suggest statistically significant differences in comparisons between HR+ (n = 694) and HR− (n = 96) invasive breast cancers in all women combined (Table 4). Risk of HR+ breast was significantly positively associated with ∑LMWP exposure (OR above vs below median = 1.30, 95% CI 1.05–1.60), but this positive association was also observed for HR− breast cancer (OR above vs below median = 1.28, 95% CI 0.81–2.00) (Pheterogeneity = 0.95). In contrast, the risk of HR− breast cancers was increased in association with phthalic acid exposure (OR above vs below median = 1.59, 95% CI 1.01–2.46) but not HR+ breast cancer (OR above vs below median = 1.02, 95% CI 0.83–1.26) (Pheterogeneity = 0.08). Above the median exposure to ∑LMHMPA was also associated with HR− breast cancer (OR above vs below median = 1.54, 95% CI 0.98–2.41, P = 0.06); this increased risk was observed in whites, Native Hawaiians, and African Americans and Latinos combined, but it was particularly prominent among Native Hawaiian women (OR above vs below median = 4.92, 95% CI 1.33–18.2) (Supplementary Table 3).
Phthalate-breast cancer associations did not differ by years of follow-up after urine collection, or between all (invasive and in situ combined) versus invasive breast cancers only (data not shown). There were also no suggestive differences in breast cancer associations with phthalate exposures by use of hormone therapy or BMI at urine collection (data not shown). However, in the subgroup of women with data on WHR, breast cancer risks increased among those with high (≥ 0.854) WHR in association with exposure to MEHP%, ratio variables of MEHP to secondary metabolites of DEHP, phthalate acid, and ∑LMHMPA (Table 5). The ORs for ∑LMHMPA were elevated among those with high WHR (ORT3 vs T1 = 2.01, 95% CI 1.22–3.32, Ptrend = 0.01) but not among those with low (< 0.854) WHR (ORT3 vs T1 = 0.79, 95% CI 0.47–1.31) (Pheterogeneity = 0.01).
This nested case-control study adds new information on the association of phthalate exposures and breast cancer risk. This is only the second study with data on pre-diagnostic urinary exposures, and the first to examine potential racial/ethnic differences in risk associations in a single study population where similar study methods were used. Although our study still lacked adequate statistical power to detect statistically significant differences in racial/ethnic-specific comparisons, the suggestive risk pattern differences between whites and nonwhites, and between the different groups of nonwhites (Native Hawaiians, African Americans, Latinos, and Japanese Americans), highlight the challenges in studying these ubiquitous chemicals in relation to breast cancer risk. The suggestive differences in associations between HR+ and HR− breast cancer and among those with high WHR warrant confirmation in future studies.
The most consistent findings in all women combined were the inverse associations with MBzP, the increased risks associated with high MEHP%, and high ratios of DEHP metabolites, namely MEHP to MEOHP and MEHHP. The significant inverse association between MBzP exposure and risk is consistent with results observed in Mexico  and among white women in the WHI  and LIBCSP . While estrogenic effects of BBP, the precursor of MBzP, have been reported in cell culture and other experimental studies [31,32,33], these estrogenic effects may not be applicable to in vivo environments , and anti-estrogenic effects of BBP have been reported in some studies .
Our findings on the individual DEHP metabolites and ratios of primary to secondary DEHP metabolites suggest potential racial/ethnic differences in exposure and/or metabolism. Exposures to all four DEHP metabolites were inversely associated with risk in whites but were positively associated with risk in Native Hawaiians; the results in whites are supportive of the null/inverse associations reported previously in whites [13, 14], while the results in Native Hawaiians are broadly consistent with findings among Mexican  and Native Alaska women . Risk in all women appeared to be higher in association with high MEHP% and high ratios of MEHP/MEOHP and MEHP/MEHHP. As shown in Fig. 1, DEHP and other high molecular weight long-branch phthalates undergo several biotransformations, including further hydroxylation and oxidation before excretion [36,37,38]. MEHP, the hydrolysis product of DEHP, is not a major metabolite, while the oxidative metabolites (MEHHP, MEOHP, and MECPP) are the major metabolites of DEHP. As MEHP may be more toxic than the oxidative metabolites, higher MEHP% and higher ratios of MEHP to MEHHP or MEHP to MEOHP may implicate greater physiologic effects as compared with individual metabolites; this seemed to be particularly apparent among women with high WHR. While we are not aware of previous investigations of MEHP% and ratios of DEHP metabolites in relation to breast cancer risk, higher MEHP% and higher MEHP ratios may represent more toxicity in terms of male infertility , obesity based on BMI and waist circumferences , inflammation, and oxidative stress represented by gamma glutamyltransferase . Oxidative stress-stimulated inflammation may be part of the etiologic pathway for DEHP-induced carcinogenesis.
Our finding between HR+ breast cancer and ∑LMWP is supportive of suggestive effects of phthalates on ER+ breast cancer [32, 40]. In contrast, the risk associations with exposure to phthalic acid, ∑LMHMPA, and ∑LMHMmolar appeared to be stronger for HR− than HR+ breast cancers but this was based on only 96 HR− breast cancers. Although these risk associations were particularly strong among Native Hawaiians, we are cautious in our interpretation because of the very wide confidence intervals (Supplementary Table 3). Further investigations of whether HR− breast cancer may be differentially affected by phthalate exposures are needed, as fewer than 300 HR− breast cancers have been investigated in this (n = 96) and combined previous studies on urinary phthalates (n = 58 in WHI  and n = 109 in LIBCSP ).
Finally, although there was no significant evidence of statistical heterogeneity by race/ethnicity, risk association patterns between phthalate exposures and breast cancer risk were not uniform across the racial/ethnic groups. In particular, risk associations were more prominent among Native Hawaiians; in that group, risk was positively associated with eight of the ten individual phthalate metabolites, phthalic acids, and total phthalates (∑LMHMPA, ∑LMHMmolar) while the associations were less consistent in the other racial/ethnic groups. We are not aware of previous studies on exposures of phthalates among Native Hawaiians and Japanese Americans in the USA, but the divergent risk patterns may be due to differences in their use of phthalate containing products as suggested by racial/ethnic differences in profiles of individual urinary phthalates (Supplementary Table 1). Variations in genetic susceptibility to metabolize phthalates and other factors affecting metabolism may be also important and warrant investigation.
Strengths of this study include being the first prospective study to investigate urinary phthalate metabolites among five racial/ethnic groups in the same study, and providing the first of such data in Native Hawaiians and Japanese Americans, and carefully considering potential confounders (Supplementary Table 4). We used a highly sensitive assay and carefully examined risk associations with both individual metabolites and all metabolites combined as they may represent different sources and routes of exposure. However, there are several limitations. We relied on a single urine sample measurement; the within-person variability is modest for phthalate metabolites in this study as in other studies  and may have reduced the statistical power of our study. Sample sizes of African Americans and Latinas were also modest, which precluded examining their risk patterns separately. In addition, information on pre-diagnostic WHR was available on only a subset of breast cancer cases. Due to the number of tests performed, some of the findings are likely due to chance. All the samples from Los Angeles County were first morning samples, whereas almost all urine specimens from Hawaii were overnight specimens. While misclassification of exposure is unavoidable, we believe that non-differential misclassification of exposure would tend to attenuate the overall results and underestimate the true association.
Results from this large nested case-control study suggest that exposure to phthalates may be associated with breast cancer risk. However, these relationships are complex and may differ by race/ethnicity and hormone receptor status of breast cancer, perhaps reflecting different exposure patterns to phthalate containing products as well as metabolism. Our novel findings of increased risks with ratios of primary DEHP metabolite (MEHP) to secondary DEHP metabolites (MEHHP, MEOHP, MECPP) emphasize the need of further investigation of not only the individual metabolites, but also ratios of the metabolites as well as summation of these metabolites. While additional studies with pre-diagnostic samples are needed to fill our gaps in understanding the effects of these ubiquitous chemicals, continued efforts to reduce exposure burden to phthalates are needed.
Availability of data and materials
Request for additional details of the data used in the manuscript can be directed to the corresponding author.
Butyl benzyl phthalate
Body mass index kg/m2
Coefficient of variation
Hormone receptor-positive status
Hormone receptor-negative status
High molecular weight
Lower limit of detection
Long Island Breast Cancer Study Project
Low molecular weight
Mono 2-carboxy-hexyl phthalate
- NHANES III:
Third National Health and Nutrition Examination Survey
Neighborhood socioeconomic status
Surveillance, Epidemiology, and End Results Program
Women’s Health Initiative
Wang Y, Zhu H, Kannan K. A review of biomonitoring of phthalate exposures. Toxics. 2019;7(2):21. https://doi.org/10.3390/toxics7020021.
Silva MJ, Malek NA, Hodge CC, Reidy JA, Kato K, Barr DB, Needham LL, Brock JW. Improved quantitative detection of 11 urinary phthalate metabolites in humans using liquid chromatography-atmospheric pressure chemical ionization tandem mass spectrometry. J Chromatogr B Analyt Technol Biomed Life Sci. 2003;789(2):393–404. https://doi.org/10.1016/S1570-0232(03)00164-8.
Rodgers KM, Udesky JO, Rudel RA, Brody JG. Environmental chemicals and breast cancer: an updated review of epidemiological literature informed by biological mechanisms. Environ Res. 2018;160:152–82. https://doi.org/10.1016/j.envres.2017.08.045.
Woodruff TJ. Bridging epidemiology and model organisms to increase understanding of endocrine disrupting chemicals and human health effects. J Steroid Biochem Mol Biol. 2011;127(1–2):108–17. https://doi.org/10.1016/j.jsbmb.2010.11.007.
Silva MJ, Barr DB, Reidy JA, Malek NA, Hodge CC, Caudill SP, Brock JW, Needham LL, Calafat AM. Urinary levels of seven phthalate metabolites in the U.S. population from the National Health and Nutrition Examination Survey (NHANES) 1999-2000. Environ Health Perspect. 2004;112(3):331–8.
Liu X, Cao YF, Ran RX, Dong PP, Gonzalez FJ, Wu X, Huang T, Chen JX, Fu ZW, Li RS, Liu YZ, Sun HZ, Fang ZZ. New insights into the risk of phthalates: inhibition of UDP-glucuronosyltransferases. Chemosphere. 2016;144:1966–72. https://doi.org/10.1016/j.chemosphere.2015.10.076.
Frederiksen H, Skakkebaek NE, Andersson AM. Metabolism of phthalates in humans. Mol Nutr Food Res. 2007;51(7):899–911. https://doi.org/10.1002/mnfr.200600243.
Morgan M, Deoraj A, Felty Q, Roy D. Environmental estrogen-like endocrine disrupting chemicals and breast cancer. Mol Cell Endocrinol. 2017;457:89–102. https://doi.org/10.1016/j.mce.2016.10.003.
Fabjan E, Hulzebos E, Mennes W, Piersma AH. A category approach for reproductive effects of phthalates. Crit Rev Toxicol. 2006;36(9):695–726. https://doi.org/10.1080/10408440600894914.
Brody JG, Rudel RA. Environmental pollutants and breast cancer. Environ Health Perspect. 2003;111(8):1007–19. https://doi.org/10.1289/ehp.6310.
Lopez-Carrillo L, Hernandez-Ramirez RU, Calafat AM, Torres-Sanchez L, Galvan-Portillo M, Needham LL, Ruiz-Ramos R, Cebrian ME. Exposure to phthalates and breast cancer risk in northern Mexico. Environ Health Perspect. 2010;118(4):539–44. https://doi.org/10.1289/ehp.0901091.
Holmes AK, Koller KR, Kieszak SM, Sjodin A, Calafat AM, Sacco FD, Varner DW, Lanier AP, Rubin CH. Case-control study of breast cancer and exposure to synthetic environmental chemicals among Alaska native women. Int J Circumpolar Health. 2014;73(1):25760. https://doi.org/10.3402/ijch.v73.25760.
Parada H Jr, Gammon MD, Chen J, Calafat AM, Neugut AI, Santella RM, Wolff MS, Teitelbaum SL. Urinary phthalate metabolite concentrations and breast cancer incidence and survival following breast cancer: the Long Island Breast Cancer Study Project. Environ Health Perspect. 2018;126(4):047013. https://doi.org/10.1289/EHP2083.
Reeves KW, Diaz Santana M, Manson JE, Hankinson SE, Zoeller RT, Bigelow C, Sturgeon SR, Spiegelman D, Tinker L, Luo J, et al. Urinary phthalate biomarker concentrations and postmenopausal breast cancer risk. J Natl Cancer Inst. 2019;111(10):1059–67. https://doi.org/10.1093/jnci/djz002.
Ahern TP, Broe A, Lash TL, Cronin-Fenton DP, Ulrichsen SP, Christiansen PM, Cole BF, Tamimi RM, Sorensen HT, Damkier P. Phthalate exposure and breast cancer incidence: a Danish nationwide cohort study. J Clin Oncol. 2019;37(21):1800–9. https://doi.org/10.1200/JCO.18.02202.
Kolonel LN, Henderson BE, Hankin JH, Nomura AM, Wilkens LR, Pike MC, Stram DO, Monroe KR, Earle ME, Nagamine FS. A multiethnic cohort in Hawaii and Los Angeles: baseline characteristics. Am J Epidemiol. 2000;151(4):346–57. https://doi.org/10.1093/oxfordjournals.aje.a010213.
Epplein M, Franke AA, Cooney RV, Morris JS, Wilkens LR, Goodman MT, Murphy SP, Henderson BE, Kolonel LN, Le Marchand L. Association of plasma micronutrient levels and urinary isoprostane with risk of lung cancer: the multiethnic cohort study. Cancer Epidemiol Biomark Prev. 2009;18(7):1962–70. https://doi.org/10.1158/1055-9965.EPI-09-0003.
Derby KS, Cuthrell K, Caberto C, Carmella SG, Franke AA, Hecht SS, Murphy SE, Le Marchand L. Nicotine metabolism in three ethnic/racial groups with different risks of lung cancer. Cancer Epidemiol Biomark Prev. 2008;17(12):3526–35. https://doi.org/10.1158/1055-9965.EPI-08-0424.
Song Y, Hauser R, Hu FB, Franke AA, Liu S, Sun Q. Urinary concentrations of bisphenol A and phthalate metabolites and weight change: a prospective investigation in US women. Int J Obes. 2014;38(12):1532–7. https://doi.org/10.1038/ijo.2014.63.
Sun Q, Cornelis MC, Townsend MK, Tobias DK, Eliassen AH, Franke AA, Hauser R, Hu FB. Association of urinary concentrations of bisphenol A and phthalate metabolites with risk of type 2 diabetes: a prospective investigation in the Nurses’ Health Study (NHS) and NHSII cohorts. Environ Health Perspect. 2014;122(6):616–23. https://doi.org/10.1289/ehp.1307201.
Bang du Y, Lee IK, Lee BM: Toxicological characterization of phthalic acid. Toxicol Res 2011, 27(4):191–203.
Albro PW, Jordan S, Corbett JT, Schroeder JL. Determination of total phthalate in urine by gas chromatography. Anal Chem. 1984;56(2):247–50. https://doi.org/10.1021/ac00266a029.
Austin PC. An introduction to propensity score methods for reducing the effects of confounding in observational studies. Multivariate Behav Res. 2011;46(3):399–424. https://doi.org/10.1080/00273171.2011.568786.
Joffe MM, Rosenbaum PR. Invited commentary: propensity scores. Am J Epidemiol. 1999;150(4):327–33. https://doi.org/10.1093/oxfordjournals.aje.a010011.
Conroy SM, Shariff-Marco S, Yang J, Hertz A, Cockburn M, Shvetsov YB, Clarke CA, Abright CL, Haiman CA, Le Marchand L, et al. Characterizing the neighborhood obesogenic environment in the Multiethnic Cohort: a multi-level infrastructure for cancer health disparities research. Cancer Causes Control. 2018;29(1):167–83. https://doi.org/10.1007/s10552-017-0980-1.
Harmon BE, Boushey CJ, Shvetsov YB, Ettienne R, Reedy J, Wilkens LR, Le Marchand L, Henderson BE, Kolonel LN. Associations of key diet-quality indexes with mortality in the Multiethnic Cohort: the Dietary Patterns Methods Project. Am J Clin Nutr. 2015;101(3):587–97. https://doi.org/10.3945/ajcn.114.090688.
Silva MJ, Samandar E, Preau JL Jr, Needham LL, Calafat AM. Urinary oxidative metabolites of di(2-ethylhexyl) phthalate in humans. Toxicology. 2006;219(1–3):22–32. https://doi.org/10.1016/j.tox.2005.10.018.
Koch HM, Bolt HM, Preuss R, Angerer J. New metabolites of di(2-ethylhexyl)phthalate (DEHP) in human urine and serum after single oral doses of deuterium-labelled DEHP. Arch Toxicol. 2005;79(7):367–76. https://doi.org/10.1007/s00204-004-0642-4.
Hauser R. Urinary phthalate metabolites and semen quality: a review of a potential biomarker of susceptibility. Int J Androl. 2008;31(2):112–7. https://doi.org/10.1111/j.1365-2605.2007.00844.x.
Yaghjyan L, Sites S, Ruan Y, Chang SH. Associations of urinary phthalates with body mass index, waist circumference and serum lipids among females: National Health and Nutrition Examination Survey 1999-2004. Int J Obes. 2015;39(6):994–1000. https://doi.org/10.1038/ijo.2015.8.
Chen FP, Chien MH. Lower concentrations of phthalates induce proliferation in human breast cancer cells. Climacteric. 2014;17(4):377–84. https://doi.org/10.3109/13697137.2013.865720.
Chen FP, Chien MH, Chen HY, Ng YT. Effects of phthalates on normal human breast cells co-cultured with different fibroblasts. PLoS One. 2018;13(6):e0199596. https://doi.org/10.1371/journal.pone.0199596.
Hsieh TH, Tsai CF, Hsu CY, Kuo PL, Hsi E, Suen JL, Hung CH, Lee JN, Chai CY, Wang SC, Tsai EM. n-Butyl benzyl phthalate promotes breast cancer progression by inducing expression of lymphoid enhancer factor 1. PLoS One. 2012;7(8):e42750. https://doi.org/10.1371/journal.pone.0042750.
Picard K, Lhuguenot JC, Lavier-Canivenc MC, Chagnon MC. Estrogenic activity and metabolism of n-butyl benzyl phthalate in vitro: identification of the active molecule(s). Toxicol Appl Pharmacol. 2001;172(2):108–18. https://doi.org/10.1006/taap.2001.9141.
Okubo T, Suzuki T, Yokoyama Y, Kano K, Kano I. Estimation of estrogenic and anti-estrogenic activities of some phthalate diesters and monoesters by MCF-7 cell proliferation assay in vitro. Biol Pharm Bull. 2003;26(8):1219–24. https://doi.org/10.1248/bpb.26.1219.
Erkekoglu P, Rachidi W, Yuzugullu OG, Giray B, Favier A, Ozturk M, Hincal F. 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. 2010;248(1):52–62. https://doi.org/10.1016/j.taap.2010.07.016.
Meeker JD, Calafat AM, Hauser R. Di(2-ethylhexyl) phthalate metabolites may alter thyroid hormone levels in men. Environ Health Perspect. 2007;115(7):1029–34. https://doi.org/10.1289/ehp.9852.
Hauser R, Meeker JD, Singh NP, Silva MJ, Ryan L, Duty S, Calafat AM. DNA damage in human sperm is related to urinary levels of phthalate monoester and oxidative metabolites. Hum Reprod. 2007;22(3):688–95. https://doi.org/10.1093/humrep/del428.
Ferguson KK, Loch-Caruso R, Meeker JD. Urinary phthalate metabolites in relation to biomarkers of inflammation and oxidative stress: NHANES 1999-2006. Environ Res. 2011;111(5):718–26. https://doi.org/10.1016/j.envres.2011.02.002.
Kang SC, Lee BM. DNA methylation of estrogen receptor alpha gene by phthalates. J Toxicol Environ Health A. 2005;68(23–24):1995–2003. https://doi.org/10.1080/15287390491008913.
Townsend MK, Franke AA, Li X, Hu FB, Eliassen AH. Within-person reproducibility of urinary bisphenol A and phthalate metabolites over a 1 to 3 year period among women in the Nurses’ Health Studies: a prospective cohort study. Environ Health. 2013;12(1):80. https://doi.org/10.1186/1476-069X-12-80.
We thank all the study participants in the Multiethnic Cohort Study and Diana Chingos who served as patient advocate on this study. We also thank Laurie Custer and Karly Tori for technical assistance in LCMS analysis.
This work was supported by the California Breast Cancer Research Program (22UB-2309 for all authors listed) and the National Cancer Institute (U01 CA164973 for CH, LW, LM, CC), USC Norris Comprehensive Cancer Center (Core) Support (P30 CA014089 for AW and CH), Environmental Exposures, Host Factors, and Human Disease (P30 ES007048021 for AW), and University of Hawaii Cancer Center CCSG (P30 CA71789 for AF).
Ethics approval and consent to participate
All women provided consent forms when they donated biological samples. The study protocol was approved by the Institutional Review Boards of each participating study site.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Molecular weight of phthalate metabolites and phthalic acid and samples with lower limit of detection (LLOD), and geometric mean (95%CI) of phthalate metabolites and summary variables of phthalates (μg/g creatinine) by study area and by race/ethnicity. Supplementary Table 2. Pairwise Spearman’s correlations between individual phthalate biomarkers and summation variables among controls. Supplementary Table 3. Associations of risk (OR and 95% CI) a of hormone receptor positive (HR+) and HR negative (HR-) invasive breast cancers with summary phthalate exposures by race/ethnicity. Supplemental Table 4. Distribution of covariates by tertiles of total phthalate exposure (∑LMHMPA).
About this article
Cite this article
Wu, A.H., Franke, A.A., Wilkens, L.R. et al. Urinary phthalate exposures and risk of breast cancer: the Multiethnic Cohort study. Breast Cancer Res 23, 44 (2021). https://doi.org/10.1186/s13058-021-01419-6
- Breast cancer
- Urinary phthalates
- Hormone receptor status
- Waist-hip ratio