Perfluoroalkyl substances (PFASs) as risk factors for breast cancer: a case–control study in Chinese population
Environmental Health volume 21, Article number: 83 (2022)
Perfluoroalkyl substances (PFASs) are a large family of synthetic chemicals, some of which are mammary toxicants and endocrine disruptors. Recent studies have implicated exposure to PFASs as a risk factor for breast cancer in Europe and America. Little is known about the role of PFASs with respect to breast cancer in the Chinese population.
Participants who were initially diagnosed with breast cancer at Tianjin Medical University Cancer Institute and Hospital between 2012 and 2016 were recruited as cases. The controls were randomly selected from the participants with available blood samples in the Chinese National Breast Cancer Screening Program (CNBCSP) cohort. Ultimately, we enrolled 373 breast cancer patients and 657 controls. Plasma PFASs were measured by an ultra-performance liquid chromatography (UPLC) system coupled to a 5500 Q-Trap triple quadrupole mass spectrometer. A logistic regression model with least absolute shrinkage and selection operator (LASSO) regularization was used to calculate odds ratios (ORs) and 95% confidence intervals (CIs) to assess the relationships between PFASs and breast cancer. The three most predictive variables in the LASSO model were selected from 17 PFASs, which was based on the optimal penalty coefficient (λ = 0.0218) identified with the minimum criterion. Additionally, Bayesian kernel machine regression (BKMR) and quantile g-computation models were applied to evaluate the associations between separate and mixed exposure to PFASs and breast cancer.
Perfluorooctanesulfonic acid (PFOS) exhibited the highest concentration in both the cases and controls. Perfluorooctanoic acid (PFOA) and perfluoro-n-decanoic acid (PFDA) were positively associated with breast cancer, and perfluoro-n-tridecanoic acid (PFTrDA) was negatively associated with breast cancer according to both the continuous-PFASs and the quartile-PFASs logistic regression models. Of note, PFOA was associated with the occurrence of estrogen receptor (ER)-, progesterone receptor (PR)-, and human epidermal growth factor receptor 2 (HER2)-positive breast cancer (ORER+ = 1.47, 95% CI: 1.19, 1.80; ORPR+ = 1.36, 95% CI: 1.09, 1.69; ORHER2 = 1.62, 95% CI: 1.19, 2.21).
Overall, we observed that PFASs were associated with breast cancer in Chinese women. Prospective cohort studies and mechanistic experiments are warranted to elucidate whether these associations are causal.
Perfluoroalkyl substances (PFASs) are a group of fluorinated chemicals that have been used in consumer products and industrial applications for their stain-, grease-, and water-repellent properties since the 1950s. Most humans carry determinable burdens of PFASs. The extensive study by the U.S. National Health and Nutrition Examination Survey (NHANES) in 2013–2014 detected these compounds in the blood of over the 98% of the participants . More than 90% of samples from the general population have detectable perfluorooctanoic acid (PFOA), perfluorooctanesulfonic acid (PFOS), and PFOS alternatives in China [2,3,4]. PFASs have attracted more concern and have been actively researched because of their extensive presence in the environment, persistent accumulation, and toxicity due to their long serum half-lives . Of note, PFOS and PFOA have been added to the Stockholm Convention on persistent organic pollutants (POPs) . Epidemiological and experimental studies have associated PFAS exposure with a series of toxicological effects that include changes in glucose and cholesterol, thyroid hormone alterations, reproductive and developmental toxicity, hepatotoxicity, and carcinogenic effects [7,8,9].
Breast cancer is the most common cancer among Chinese women. China accounts for 12.2% of new breast cancer diagnoses and 9.6% of all breast cancer deaths globally [10, 11]. In addition to genetic factors such as mutations in BRCA-1/2, risk factors for breast cancer include early menarche, late menopause, postmenopausal obesity, smoking, drinking alcohol, and high fat intake [12, 13]. Known risk factors can only explain less than one-third of the incidences, and lifestyle and environmental exposure can also alter the risk of breast cancer . Environmental chemicals may play an important role in breast cancer. PFASs have been confirmed to exert weak xenoestrogen effects . In addition, the International Agency for Research on Cancer (IARC) classified PFOA as a possible human carcinogen (Group 2B) in 2017, based on limited evidence that PFOA can cause testicular and renal cancer . Many epidemiological studies have been implemented in Europe and the United States to assess breast cancer risks associated with PFAS exposure. Nonetheless, the literature is sparse, and the findings are controversial. Furthermore, little is known about the role of PFASs in the development of breast cancer in the Chinese population. Of note, several distinct characteristics have been observed in Chinese breast cancer patients, including a younger diagnostic age (45–55 years) compared to patients in Western countries and a pattern of risk factors only partially matching the known factors in high-income countries .
Given the limited research in this area and the biological plausibility of an association, we conducted a case–control study to evaluate the associations between individual and joint PFAS exposure and breast cancer by several statistical models.
Study design and participants
This study was approved by the Ethics Board of the Tianjin Medical University Cancer Institute and Hospital. Participants who were initially diagnosed with breast cancer by pathological examinations at Tianjin Medical University Cancer Institute and Hospital between January 2012 and December 2016 were recruited as cases. The blood samples of the patients were collected within a week after breast cancer diagnosis and before treatment. The controls included in this case–control study were randomly selected from the participants in the Chinese National Breast Cancer Screening Program (CNBCSP) cohort from a time period similar to that of the cases [16, 17]. We identified 380 cases and 798 controls with available blood samples. From these participants, those who had not completed the dietary questionnaire were excluded (cases = 7, controls = 141). In total, 373 cases and 657 controls were included in this study. Moreover, the clinical data, including the morphology, metastasis, staging, and hormone receptor status of breast tumors, were collected. The flowchart summarizing the inclusion process of the case–control study is presented in Fig. 1.
According to the standard for the molecular typing of breast cancer, tumors were identified as luminal A, luminal B, human epidermal growth factor receptor 2 (HER2), and basal-like subtypes . Luminal A-like tumors were defined as estrogen receptor (ER)-positive, progesterone receptor (PR)-positive, HER2-negative, and low Ki67 expression. The remaining ER-positive tumors were classified as the luminal B subtype. HER2-driven tumors were categorized as ER-negative, PR-negative, and HER2 overexpression or amplification. Basal-like breast cancer was negatively expressed in ER, PR, and HER2 subtypes .
Sample collection and analysis of PFASs in plasma
Blood samples were collected by licenced phlebotomists into 10 mL BD® tubes using standard phlebotomy techniques. Blood samples were collected from cases after diagnosis and before treatment. In the control group, the time of blood sample collection was the same as the screening time in the CNBCSP. The samples were stored at -80 °C until analysis. The detection methodology of PFASs has been previously described  and is detailed in the Supplemental Materials. The PFASs were quantified using an ultra-performance liquid chromatography (UPLC) system coupled to a 5500 Q-Trap triple quadrupole mass spectrometer (AB Sciex, Canada). PFASs include perfluorocarboxylated acids (PFCAs) and perfluorosulfonated acids (PFSAs). In summary, the levels of 38 PFASs were measured in this study, including five linear PFSAs (PFOS, PFBS, PFHxS, PFHpS, and PFDS), 11 linear PFCAs (PFOA, PFBA, PFPeA, PFHxA, PFHpA, PFNA, PFDA, PFUdA, PFDoA, PFTrDA, and PFTeDA), 17 PFOA and PFOS isomers (P1MHpS, P3MHpS, P3MHpA, P4MHpS, P4MHpA, P5MHpS, P5MHpA, P6MHpS, P6MHpA, P55DMHxS, P55DMHxA, P44DMHxS, P44DMHxA, P45DMHxS, P45DMHxA, P35DMHxS, and P35DMHxA), two chlorinated alternatives (11CL-PF3OUdS and 9CL-PF3ONS), and three fluorotelomer sulfonates (4:2FTS, 6:2FTS, and 8:2FTS). 13C8-PFOA was used as an internal standard, and nine isotopically labelled PFASs were applied as surrogate standards (Wellington Laboratories, Guelph, Ontario, Canada). The names, acronyms, and detection frequencies of the PFASs are presented in Table S1.
Quality assurance and control
Quality assurance and control procedures included the analysis of matrix spiking samples, procedural blanks, and the monitoring of surrogate standard recoveries. The PFASs of interest were spiked into plasma and processed in five replicates along with two controls (only surrogate standards were spiked). The mean recoveries of individual PFASs ranged from 65 to 94% after subtracting the original values determined in the plasma. Two laboratory procedural blanks were processed in combination with every batch of 10 samples. All data were blank corrected prior to analysis. The recoveries of relevant surrogate standards in the samples ranged from 51 to 92%. Batch-specific blank correction was performed before data analysis. The recovery and standard deviation (SD) of each target compound are shown in Table S2.
Concentrations below the limit of detection (LOD) were imputed as the detection limit divided by the square root of two. Normality was assessed using the Shapiro–Wilk test. The continuous variables were aggregated using the median and quartile range. Frequencies and percentages were applied for categorical variables. The Mann–Whitney-Wilcoxon test and Fisher’s exact test were used to make statistical comparisons of the demographic characteristics of the various groups. Spearman’s correlation was applied to describe pairwise relationships among the 17 PFASs (Figure S1). The level of PFASs was transformed using a natural logarithm (Ln) for a more normal distribution in the continuous model. Then, scatter plots of the residual distributions were used to detect the linearity of the continuous exposure variables, which is shown in Figure S2. Unconditional logistic regression models were applied to estimate odds ratios (ORs) and 95% confidence intervals (CIs) of breast cancer with respect to Ln (PFASs) concentrations. The ORs of PFAS concentrations that were analyzed as continuous variables with one SD increase in 17 Ln (PFASs, ng/ml) and quartile variables were both reported in this study by the logistic regression model. Variables were categorized according to the distribution of the PFAS concentrations in the controls, where the lowest quartile group was used as the reference group.
The least absolute shrinkage and selection operator (LASSO) regression model was used to minimize a large number of variables in the model and, importantly, had the unique feature of penalizing the absolute value of a regression coefficient [20, 21]. The regression coefficient of the regression model was based on the value of λ (λmin or λ1se). The larger the penalization is along with the greater the shrinkage of coefficients. When the estimates of weaker factors reached zero, the unnecessary variables were automatically removed, and only the strongest variables remained in the model. In this study, the most predictive PFAS variables were selected by λ (λ1se = 0.0218, Figure S3).
The Bayesian kernel machine regression (BKMR) model was used to evaluate the individual and combined effects of exposure to PFASs on breast cancer [22,23,24]. The 17 compounds were divided into three categories: PFCAs, PFSAs, and PFSA isomers. The ∑PFCAs included PFOA, PFNA, PFDA, PFUdA, PFDoA, PFTrDA, and PFTeDA. The ∑PFSAs included PFHxS, PFHpS, and PFOS. The ∑PFSA isomers included P3MHpS, P4MHpS, P5MHpS, P6MHpS, P45DMHxS, 11CL-PF3OUdS, and 9CL-PF3ONS. The three categories of PFAS mixtures were flexibly modelled and individual PFAS-breast cancer associations were estimated using a kernel function in BKMR. This method allows for potential nonlinear exposure–response functions. In addition, the variable posterior inclusion probabilities (PIPs) estimated by BKMR for each single PFAS represented the relative importance of the PFASs to the overall mixture effect. We ln-transformed and standardized PFAS concentrations in the BKMR model.
To further improve the robustness of the joint effects of PFAS mixtures calculated by BKMR, quantile g-computation approach was used. Quantile g-computation is a recently developed approach that combines the inferential simplicity of weighted quantile sum (WQS) regression with the flexibility of g-computation [25, 26]. In this study, quantile g-computation provided an estimate of the total effect of the three categories of PFAS mixtures on breast cancer and weights for the individual PFASs, which represented their relative contribution to the overall mixture effect. The weights represented the percentage of a single PFAS’s positive or negative relationship with breast cancer, with the weights in each direction adding up to 1.0.
The identification of confounders was performed according to the following criteria. First, variables were hypothesized to be risk factors for outcomes. Second, potential confounding variables identified by professional knowledge, traditional experience, and prior research were included in the adjusted models. Third, some important covariates were included in the adjusted models regardless of how they influenced the estimates [27,28,29]. The final multivariable models included terms for age at baseline (continuous), BMI (continuous), smoking (nonsmoker vs. current or former smoker), age at menarche (< 13 vs. ≥ 13 years), age at menopause (pre/peri-menopausal, < 51, or ≥ 51 years), parity (< 2 vs. ≥ 2), breastfeeding duration (< 13 vs. ≥ 13 months), use of estrogen or estrogen replacement therapy (no vs. yes), family history of breast cancer (no vs. yes), education level (< 13 vs. ≥ 13 years of education), income (< 1000, 1000–1999, 2000–2999, or ≥ 3000 Chinese RMB per month), and food consumption (rare, occasional, or often). Consumption of red meat, pickled, fried, smoked, and barbequed foods were separated into three levels: rare (defined as less than twice per week), occasional (defined as three or four times per week), and often (defined as more than five times per week). Alcohol consumption (nondrinker vs. current or former drinker) was not adjusted due to the small number of current and former drinkers (n = 5) in the case group. The LASSO, BKMR, and quantile g-computation models were adjusted for the same covariates.
All analyses were conducted using SPSS 26.0 and R version 4.2.0 (“glmnet”, “bkmr”, and “qgcomp” packages). A two-sided test and P-value < 0.05 were considered statistically significant.
The characteristics of the study participants are presented in Table 1. This study contained 373 cases and 657 controls. The median age was 50.6 years (IQR, 44.4–59.2 years) among the breast cancer patients and 53.0 years (IQR, 50.0–57.0 years) among the controls. Additionally, cases were more likely than controls to be current or former smokers (8.0% vs. 4.9%) and have a higher percentage (6.2% vs. 3.0%) of a family history of breast cancer. In the case group, the proportion of women who had more than two children and a longer breastfeeding duration (≥ 13 months) was higher than that in the control group. The distribution of estrogen replacement therapy use in breast cancer patients (4.1%) was similar to that in controls (3.4%) in the present study. There was no significant difference between the two groups in the level of education or monthly household income per capita (p > 0.05). The proportion of daily consumption of red meat, fried food, and smoked food did not differ between the two groups (p > 0.05). The frequencies of pickled and barbequed food intake were higher in the breast cancer patients than in the controls (p < 0.05).
A total of 38 PFASs were detected in this study, and the detection rates of 17 PFASs were greater than 80%. The concentrations of the 17 PFASs (ng/mL) in the plasma samples are shown in Table 2. The proportions above the LOD of the top seven PFASs in the control samples were 100% for PFOA, 100% for 9CL-PF3ONS, 99.90% for PFHxS, 99.61% for PFDA, 99.32% for 11CL-PF3OUdS, 99.22% for PFOS, and 99.22% for PFUdA. PFOS exhibited the highest concentration in both cases and controls, followed by 9CL-PF3ONS and PFOA. Cases had higher levels of PFOA and PFDA than controls.
As shown in Figure S1, the Spearman’s rank correlations among 17 PFASs were calculated. As several of the PFASs were highly correlated, the LASSO regression model was used to reduce data dimensions and determine the most powerful PFAS variables before logistic regression model was performed. The crude and adjusted ORs of breast cancer for the 17 PFASs as continuous and quantile variables in a logistic regression model without LASSO regularization are shown in Table S3. We further evaluated the effects of the PFASs as continuous variables using LASSO regression followed by logistic regression. The LASSO logistic model minimized multicollinearity among the 17 PFASs and assessed the relationship between the selected PFASs and breast cancer. We used a LASSO model to build a PFASs classifier (Figure S3) and selected the most powerful PFAS variables (λ = 0.0218, n = 3). Finally, the logistic regression model was used to evaluate the associations between the selected variables (PFOA, PFDA, and PFTrDA) and breast cancer. The results showed that PFOA (OR = 3.32, 95% CI: 2.32, 4.75) and PFDA (OR = 2.22, 95% CI: 1.55, 3.17) were positively associated with breast cancer (Table 3). We also found that an increase in the SD in Ln (PFTrDA) was negatively associated with breast cancer (OR = 0.03, 95% CI: 0.02, 0.06; Table 3). The results from logistic regression with LASSO-selected PFASs showed that both the positive and negative associations were enhanced (ORPFOA, 3.32 vs. 1.57; ORPFDA, 2.22 vs. 1.27; ORPFTrDA, 0.03 vs. 0.10) compared to the results from single-pollutant logistic regression. To further improve the robustness of the results of the selected PFASs calculated by the LASSO regression model, BKMR was used to estimate individual PFAS-breast cancer associations (Fig. 2). Among the 17 PFASs, the substances that had the highest conditional PIPs were PFOA (1.00), PFDA (1.00), PFUdA (1.00), and PFTrDA (1.00) (Table S4). Regarding the association between individual PFAS and breast cancer, while setting all other exposure variables at their median values, positive and significant associations were observed between PFOA and PFDA and breast cancer (Fig. 2A, B). Of note, a negative association existed between PFTrDA and breast cancer (Fig. 2A, B). Compared with the ORs shown in Table 3, the results calculated by BKMR were consistent with those calculated by LASSO regression, revealing the robustness of the associations between the PFASs and breast cancer.
In addition to exploring the associations between individual PFAS exposure and breast cancer, the relationships between the three categories of PFASs including ∑PFCAs (PFOA, PFNA, PFDA, PFUdA, PFDoA, PFTrDA, and PFTeDA), ∑PFSAs (PFHxS, PFHpS, and PFOS), and ∑PFSA isomers (P3MHpS, P4MHpS, P5MHpS, P6MHpS, P45DMHxS, 11CL-PF3OUdS, and 9CL-PF3ONS), and breast cancer were also investigated in this study based on the BKMR and quantile g-computation models. As shown in Figure S4 and Figure S5, the results calculated by the two models were consistent. According to the results calculated by the quantile g-computation model, ∑PFCAs (OR = 0.86, 95% CI: 0.83, 0.88) and ∑PFSAs (OR = 0.92, 95% CI: 0.89, 0.95) were negatively associated with breast cancer. Exposure to PFSA isomers had no significant association with breast cancer (OR = 0.97, 95% CI: 0.94, 1.01).
Based on the variables selected by the LASSO model (PFOA and PFDA), the logistic regression model was further used to evaluate the associations between the plasma level of one SD increase in Ln (PFOA, ng/ml) and Ln (PFDA, ng/ml) and breast cancer stratified by subtypes or receptor expression status (Table 4). Interestingly, we observed that PFOA exposure was more likely to be associated with ER-positive (OR = 1.47, 95% CI: 1.19, 1.80) and PR-positive (OR = 1.36, 95% CI: 1.09, 1.69) breast cancer. Additionally, exposure to PFOA also showed positive associations with the luminal A (OR = 1.36, 95% CI: 1.08, 1.71) and HER2 (OR = 1.62, 95% CI: 1.19, 2.21) subtypes of breast cancer. We found that PFDA was positively related to ER-positive (OR = 1.19, 95% CI: 1.01, 1.41) breast cancer only in the crude model.
The case–control study demonstrated that plasma concentrations of PFOA and PFDA were positively associated with breast cancer, and PFTrDA was negatively associated with breast cancer in Chinese women. PFOA was more likely to be associated with the ER-, PR-, and HER2-positive breast cancer. The associations between joint and individual exposure to PFASs and breast cancer were examined by LASSO regression, BKMR, and quantile g-computation models, which proved the robustness of the conclusions.
The study population was exposed to a variety of perfluorinated compounds. PFOS, 9CL-PF3ONS, and PFOA were the top three perfluorinated compounds with the highest serum concentrations. The median serum concentrations (data not shown) measured in the study population were 12.18 ng/mL for PFOS, 3.35 ng/mL for PFOA, 0.94 ng/mL for PFDA, and 0.38 ng/ml for PFTrDA; these fit within the range of exposure levels reported in the C8 Health Project, a cross-sectional study conducted with 408 women in Shenyang, China, in which mean PFAS concentrations of 19.56 ng/ml for PFOS, 6.41 ng/ml for PFOA, 0.94 ng/ml for PFDA, and 0.68 ng/ml for PFTrDA were found .
This case–control study supported and extended the epidemiological findings on PFAS-breast cancer associations. PFOA and PFDA were positively associated with breast cancer based on this study. Two prospective studies [26, 31] both reported a positive association between PFOA and breast cancer, which was consistent with our findings. Moreover, a case–control study from the Greater Manila Area  found positive associations between breast cancer and the 3rd [OR (95% CI) = 4.09 (1.03, 21.08)] and 4th quantile groups [OR (95% CI) = 9.26 (2.54, 45.10)] of PFDA concentrations compared to the reference group. Furthermore, negative associations between PFOS and PFTrDA and breast cancer were observed both in our research and in a cross-sectional study conducted in Japan (ORPFOS = 0.15, 95% CI: 0.06, 0.39; ORPFTrDA = 0.14, 95% CI: 0.05, 0.35) . In contrast to our results, no significant association of PFOA with breast cancer was observed in case–control study with an Inuit population  and the California Teachers Study . The plasma concentration of PFOA in the Inuit population study  fit within the range of exposure levels of our population, and the inconsistencies may be explained by different genetic backgrounds and susceptibility to environmental contaminants. It should be noted that the main difference between our study population and the California population  was the time of blood sampling, which was performed after the cases received treatment. Therefore, the effect of cancer treatment on internal PFAS levels could not be eliminated. The discrepancies found in these epidemiological studies may reflect the different exposure profiles and susceptibility of the study population, as well as the characteristics of the study design. Notably, the positive association between PFOS and breast cancer observed in the Inuit population study may be partly explained by the high PFOS exposure. There are suspicions that several other POPs, which are present in Inuit people at considerably greater levels, may also have an impact on breast cancer. The risk of breast cancer in Inuit women may be increased by exposure to other carcinogenic or heterogeneous substances, such as PCBs .
Another interesting finding of this study is that PFOA is more likely to be associated with ER-, PR-, and HER2-positive breast cancer. In the French E3N cohort, low levels of PFOS and PFOA were related to receptor-negative tumors, which highlights the need to consider PFAS exposure as a potential risk factor for breast cancer . No association was found between the PFOA serum concentration and ER-positive breast cancer in a case–control study of Taiwanese women either in the young (age ≤ 50 years, OR = 1.41, 95% CI: 0.77, 2.56) or older groups (age > 50 years, OR = 0.70, 95% CI: 0.35, 1.42) . The inconsistency may be explained by the lower PFOA concentrations (1.77 ng/mL) reported in the Taiwanese study compared to those reported in our study (3.35 ng/mL). Future epidemiological studies are warranted to validate the current conclusions.
Furthermore, the negative associations between PFCAs and PFSAs and breast cancer were observed in our study. The mechanisms of a possible negative association are not well understood. Nevertheless, some suggestive evidence has been reported. This study indicated that PFTrDA may contribute most to the negative association between PFCAs and breast cancer based on the results calculated by quantile g-computation models. We speculated that the negative relationship may be partly explained by the change in thyroid function induced by PFTrDA exposure. T4 is a proliferative factor in vitro for breast cancer cells . A previous study showed that the elimination of the clinical action of endogenous T4 in breast cancer patients has favorably affected the course of cancer . Of note, Ji et al. revealed that the concentration of PFTrDA was negatively correlated with total T4 in Korean females . Our study also revealed that PFHxS contributed most to the negative association between PFSAs and breast cancer. Alterations in lipid metabolism induced by PFHxS may be involved in the negative relationship. LDL could stimulate the migration of MDA-MB-231 cells . The overexpression of the STARD3 gene, which encodes a cholesterol-binding membrane protein, could decrease cell adhesiveness, promote breast cancer metastases, and be associated with a poor prognosis for breast cancer patients . Nevertheless, previous studies have reported significant negative associations of PFHxS with TC, and non-HDL [43, 44]. Moreover, differences between the results of single and mixed PFASs exposure may be partly related to statistical models. Individual components may act differently in the presence of additional chemicals, making the overall evaluation of mixtures in the environment difficult . To be prudent, the results for PFAS must be considered with caution.
Stratified analyses of parity and breastfeeding duration were conducted in our study. We found that multi-parity women had a lower risk of breast cancer when exposed to PFOA or PFTrDA than nulliparous or single-birth women (Table S5). Women who had breastfed for a longer period of time (≥ 13 months) had a decreased risk of breast cancer than those who breastfed for a shorter amount of time due to PFOA or PFTrDA exposure (Table S6). The findings partly supported the hypothesis that parity and breastfeeding can transfer PFASs  and are associated with lower serum PFAS concentrations [47, 48], which may reduce the risk of breast cancer. Nonetheless, we observed that the multi-parity participants or women who had breastfed for more than 13 months were more sensitive to PFDA exposure than those with fewer children or shorter periods of breastfeeding.
The potential mechanisms underlying PFAS-breast cancer relationships involve four aspects: (1) PFASs in vitro can affect hormone receptors, disrupt endocrine homeostasis, and increase the risk of adverse health effects [49,50,51]; (2) PFASs can activate the peroxisome proliferator activated receptor alpha (PPARα) and induce oxidative stress [52,53,54]; (3) PFASs can induce epigenetic alterations, such as DNA methylation and histone modifications, which are implicated in tumorigenesis ; and (4) PFASs can cause reproductive toxicity, which can delay mammary gland development and increase susceptibility to carcinogens in offspring [56, 57].
Our study has some specific strengths. First, the associations between individual and joint exposure to PFASs and breast cancer were examined by several statistical models, and the consistency reflected the robustness of the conclusions. Second, this research benefited from a large number of cases that could be further grouped by hormone receptors status, and ultimately identified the breast cancer subtypes that are more sensitive to PFAS exposure, which can provide evidence for mechanistic studies. Several limitations should be noted in this study. First, the temporal relationships between PFAS exposure and breast cancer could not be ascertained based on the cross-sectional design of the research. Second, one major route of exposure to PFASs is through the ingestion of contaminated food or water . The dietary variables included in this study may be insufficient, and data for some food categories, such as seafood, takeout food, and packaged food, which may be associated with PFASs [59, 60], have not been collected. Third, only one measurement of PFASs concentrations at the time of diagnosis among the study population was obtained, which may not accurately reflect the individuals' whole exposure scenario. However, it can be assumed that the general population may be exposed to a constant level of PFASs due to the slowly changing environmental PFASs concentrations, which were caused by the PFASs’ resistance to hydrolysis, photolysis, and microbial degradation because of the high-energy carbon–fluorine bond , and a long half-life in the human body . To corroborate our findings, longitudinal follow-up studies from adolescence onwards are needed, and further exploration is necessary to determine whether any mediators cause or influence the observed relationships between PFAS levels and breast cancer.
In conclusion, we found that the plasma concentrations of PFOA and PFDA were positively associated with breast cancer, and PFTrDA was negatively associated with breast cancer in Chinese women. PFOA was more likely to be associated with ER-, PR-, and HER2-positive breast cancer. These associations warrant attention in prospective cohort studies and well-designed toxicological experiments to elucidate whether these relationships are causal.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
National Health and Nutrition Examination Survey
Persistent organic pollutants
The International Agency for Research on Cancer
Chinese National Breast Cancer Screening Program
Human epidermal growth factor receptor 2
Ultra-performance liquid chromatography
Limit of detection
Least absolute shrinkage and selection operator
Bayesian kernel machine regression
Posterior inclusion probabilities
Weighted quantile sum
Peroxisome proliferator activated receptor
Calafat AM, Kato K, Hubbard K, Jia T, Botelho JC, Wong LY. Legacy and alternative per- and polyfluoroalkyl substances in the U.S. general population: paired serum-urine data from the 2013–2014 national health and nutrition examination survey. Environ Int. 2019;131:105048.
Liu X, Chen D, Wang B, Xu F, Pang Y, Zhang L, et al. Does low maternal exposure to Per- and Polyfluoroalkyl substances elevate the risk of spontaneous preterm birth? a nested case-control study in China. Environ Sci Technol. 2020;54:8259–68.
Pan Y, Cui Q, Wang J, Sheng N, Jing J, Yao B, et al. Profiles of emerging and legacy Per-/Polyfluoroalkyl substances in matched serum and semen samples: new implications for human semen quality. Environ Health Perspect. 2019;127: 127005.
Xu H, Zhou Q, Zhang J, Chen X, Zhao H, Lu H, et al. Exposure to elevated per- and polyfluoroalkyl substances in early pregnancy is related to increased risk of gestational diabetes mellitus: a nested case-control study in Shanghai. China Environ Int. 2020;143: 105952.
Olsen GW, Burris JM, Ehresman DJ, Froehlich JW, Seacat AM, Butenhoff JL, et al. Half-life of serum elimination of perfluorooctanesulfonate, perfluorohexanesulfonate, and perfluorooctanoate in retired fluorochemical production workers. Environ Health Perspect. 2007;115:1298–305.
Miralles-Marco A, Harrad S. Perfluorooctane sulfonate: a review of human exposure, biomonitoring and the environmental forensics utility of its chirality and isomer distribution. Environ Int. 2015;77:148–59.
Rodriguez-Jorquera IA, Colli-Dula RC, Kroll K, Jayasinghe BS, Parachu Marco MV, Silva-Sanchez C, et al. Blood transcriptomics analysis of fish exposed to perfluoro alkyls substances: assessment of a non-lethal sampling technique for advancing aquatic toxicology research. Environ Sci Technol. 2019;53:1441–52.
Seo SH, Son MH, Choi SD, Lee DH, Chang YS. Influence of exposure to perfluoroalkyl substances (PFASs) on the Korean general population: 10-year trend and health effects. Environ Int. 2018;113:149–61.
Lau C, Anitole K, Hodes C, Lai D, Pfahles-Hutchens A, Seed J. Perfluoroalkyl acids: a review of monitoring and toxicological findings. Toxicol Sci. 2007;99:366–94.
Fan L, Strasser-Weippl K, Li JJ, St Louis J, Finkelstein DM, Yu KD, et al. Breast cancer in China. Lancet Oncol. 2014;15:e279-289.
Li T, Mello-Thoms C, Brennan PC. Descriptive epidemiology of breast cancer in China: incidence, mortality, survival and prevalence. Breast Cancer Res Treat. 2016;159:395–406.
Britt KL, Cuzick J, Phillips KA. Key steps for effective breast cancer prevention. Nat Rev Cancer. 2020;20:417–36.
Veronesi U, Boyle P, Goldhirsch A, Orecchia R, Viale G. Breast cancer. Lancet. 2005;365:1727–41.
Ghisari M, Eiberg H, Long M, Bonefeld-Jorgensen EC. Polymorphisms in phase I and phase II genes and breast cancer risk and relations to persistent organic pollutant exposure: a case-control study in Inuit women. Environ Health. 2014;13:19.
IARC. IARC Monographs on the Evaluation of Carcinogenic Risks to Humans. Some Chemicals Used as Solvents and in Polymer Manufacture. Lyon, France, 2017. http://monographs.iarc.fr/ENG/Monographs/vol110/index.php. 20 June 2022.
Huang Y, Tong Z, Chen K, Wang Y, Liu P, Gu L, et al. Interpretation of breast cancer screening guideline for Chinese women. Cancer Biol Med. 2019;16:825–35.
Huang Y, Dai H, Song F, Li H, Yan Y, Yang Z, et al. Preliminary effectiveness of breast cancer screening among 1.22 million Chinese females and different cancer patterns between urban and rural women. Sci Rep. 2016;6:39459.
Goldhirsch A, Winer EP, Coates AS, Gelber RD, Piccart-Gebhart M, Thurlimann B, et al. Personalizing the treatment of women with early breast cancer: highlights of the St Gallen International Expert Consensus on the primary therapy of early breast cancer 2013. Ann Oncol. 2013;24:2206–23.
Wolff AC, Hammond ME, Hicks DG, Dowsett M, McShane LM, Allison KH, et al. Recommendations for human epidermal growth factor receptor 2 testing in breast cancer: American society of clinical oncology/college of American pathologists clinical practice guideline update. J Clin Oncol. 2013;31:3997–4013.
Aung MT, Yu Y, Ferguson KK, Cantonwine DE, Zeng L, McElrath TF, et al. Cross-sectional estimation of endogenous biomarker associations with prenatal phenols, phthalates, metals, and polycyclic aromatic hydrocarbons in single-pollutant and mixtures analysis approaches. Environ Health Perspect. 2021;129:37007.
Xu C, Su X, Xu Y, Ma S, Duan W, Mo X. Exploring the associations of serum concentrations of PCBs, PCDDs, and PCDFs with walking speed in the U.S. general population: Beyond standard linear models. Environ Res. 2019;178:108666.
Bobb JF, Claus Henn B, Valeri L, Coull BA. Statistical software for analyzing the health effects of multiple concurrent exposures via Bayesian kernel machine regression. Environ Health. 2018;17:67.
Bobb JF, Valeri L, Claus Henn B, Christiani DC, Wright RO, Mazumdar M, et al. Bayesian kernel machine regression for estimating the health effects of multi-pollutant mixtures. Biostatistics. 2015;16:493–508.
Preston EV, Hivert MF, Fleisch AF, Calafat AM, Sagiv SK, Perng W, et al. Early-pregnancy plasma per- and polyfluoroalkyl substance (PFAS) concentrations and hypertensive disorders of pregnancy in the project viva cohort. Environ Int. 2022;165: 107335.
Keil AP, Buckley JP, O’Brien KM, Ferguson KK, Zhao S, White AJ. A quantile-based g-computation approach to addressing the effects of exposure mixtures. Environ Health Perspect. 2020;128:47004.
Feng Y, Bai Y, Lu Y, Chen M, Fu M, Guan X, et al. Plasma perfluoroalkyl substance exposure and incidence risk of breast cancer: a case-cohort study in the Dongfeng-Tongji cohort. Environ Pollut. 2022;306: 119345.
Liang Z, Yang Y, Qian Z, Ruan Z, Chang J, Vaughn MG, et al. Ambient PM2.5 and birth outcomes: estimating the association and attributable risk using a birth cohort study in nine Chinese cities. Environ Int. 2019;126:329–35.
Lin H, Guo Y, Zheng Y, Zhao X, Cao Z, Rigdon SE, et al. Exposure to ambient PM25 associated with overall and domain-specific disability among adults in six low- and middle-income countries. Environ Int. 2017;104:69–75.
Pan SC, Huang CC, Chin WS, Chen BY, Chan CC, Guo YL. Association between air pollution exposure and diabetic retinopathy among diabetics. Environ Res. 2020;181: 108960.
Zeng XW, Lodge CJ, Dharmage SC, Bloom MS, Yu Y, Yang M, et al. Isomers of per- and polyfluoroalkyl substances and uric acid in adults: isomers of C8 health project in China. Environ Int. 2019;133: 105160.
Mancini FR, Cano-Sancho G, Gambaretti J, Marchand P, Boutron-Ruault MC, Severi G, et al. Perfluorinated alkylated substances serum concentration and breast cancer risk: evidence from a nested case-control study in the French E3N cohort. Int J Cancer. 2020;146:917–28.
Velarde MC, Chan AFO, Sajo M, Zakharevich I, Melamed J, Uy GLB, et al. Elevated levels of perfluoroalkyl substances in breast cancer patients within the greater Manila area. Chemosphere. 2022;286: 131545.
Itoh H, Harada KH, Kasuga Y, Yokoyama S, Onuma H, Nishimura H, et al. Serum perfluoroalkyl substances and breast cancer risk in Japanese women: a case-control study. Sci Total Environ. 2021;800: 149316.
Bonefeld-Jorgensen EC, Long M, Bossi R, Ayotte P, Asmund G, Kruger T, et al. Perfluorinated compounds are related to breast cancer risk in Greenlandic Inuit: a case control study. Environ Health. 2011;10:88.
Hurley S, Goldberg D, Wang M, Park JS, Petreas M, Bernstein L, et al. Breast cancer risk and serum levels of per- and poly-fluoroalkyl substances: a case-control study nested in the California teachers study. Environ Health. 2018;17:83.
Wielsoe M, Kern P, Bonefeld-Jorgensen EC. Serum levels of environmental pollutants is a risk factor for breast cancer in Inuit: a case control study. Environ Health. 2017;16:56.
Tsai MS, Chang SH, Kuo WH, Kuo CH, Li SY, Wang MY, et al. A case-control study of perfluoroalkyl substances and the risk of breast cancer in Taiwanese women. Environ Int. 2020;142: 105850.
Tang HY, Lin HY, Zhang S, Davis FB, Davis PJ. Thyroid hormone causes mitogen-activated protein kinase-dependent phosphorylation of the nuclear estrogen receptor. Endocrinology. 2004;145:3265–72.
Hercbergs A, Mousa SA, Leinung M, Lin HY, Davis PJ. Thyroid hormone in the clinic and breast cancer. Horm Cancer. 2018;9:139–43.
Ji K, Kim S, Kho Y, Paek D, Sakong J, Ha J, et al. Serum concentrations of major perfluorinated compounds among the general population in Korea: dietary sources and potential impact on thyroid hormones. Environ Int. 2012;45:78–85.
Antalis CJ, Uchida A, Buhman KK, Siddiqui RA. Migration of MDA-MB-231 breast cancer cells depends on the availability of exogenous lipids and cholesterol esterification. Clin Exp Metastasis. 2011;28:733–41.
Vassilev B, Sihto H, Li S, Holtta-Vuori M, Ilola J, Lundin J, et al. Elevated levels of StAR-related lipid transfer protein 3 alter cholesterol balance and adhesiveness of breast cancer cells: potential mechanisms contributing to progression of HER2-positive breast cancers. Am J Pathol. 2015;185:987–1000.
Nelson JW, Hatch EE, Webster TF. Exposure to polyfluoroalkyl chemicals and cholesterol, body weight, and insulin resistance in the general U.S. population. Environ Health Perspect. 2010;118:197–202.
Bijland S, Rensen PC, Pieterman EJ, Maas AC, van der Hoorn JW, van Erk MJ, et al. Perfluoroalkyl sulfonates cause alkyl chain length-dependent hepatic steatosis and hypolipidemia mainly by impairing lipoprotein production in APOE*3-Leiden CETP mice. Toxicol Sci. 2011;123:290–303.
Rager JE, Clark J, Eaves LA, Avula V, Niehoff NM, Kim YH, et al. Mixtures modeling identifies chemical inducers versus repressors of toxicity associated with wildfire smoke. Sci Total Environ. 2021;775: 145759.
Zheng P, Liu Y, An Q, Yang X, Yin S, Ma LQ, et al. Prenatal and postnatal exposure to emerging and legacy per-/polyfluoroalkyl substances: Levels and transfer in maternal serum, cord serum, and breast milk. Sci Total Environ. 2022;812: 152446.
Liu J, Gao X, Wang Y, Leng J, Li J, Zhao Y, et al. Profiling of emerging and legacy per-/polyfluoroalkyl substances in serum among pregnant women in China. Environ Pollut. 2021;271: 116376.
Pirard C, Dufour P, Charlier C. Background contamination of perfluoralkyl substances in a Belgian general population. Toxicol Lett. 2020;333:13–21.
Kjeldsen LS, Bonefeld-Jorgensen EC. Perfluorinated compounds affect the function of sex hormone receptors. Environ Sci Pollut Res Int. 2013;20:8031–44.
Du G, Hu J, Huang H, Qin Y, Han X, Wu D, et al. Perfluorooctane sulfonate (PFOS) affects hormone receptor activity, steroidogenesis, and expression of endocrine-related genes in vitro and in vivo. Environ Toxicol Chem. 2013;32:353–60.
Liu H, Pan Y, Jin S, Sun X, Jiang Y, Wang Y, et al. Associations between six common per- and polyfluoroalkyl substances and estrogens in neonates of China. J Hazard Mater. 2021;407: 124378.
Rigden M, Pelletier G, Poon R, Zhu J, Auray-Blais C, Gagnon R, et al. Assessment of urinary metabolite excretion after rat acute exposure to perfluorooctanoic acid and other peroxisomal proliferators. Arch Environ Contam Toxicol. 2015;68:148–58.
Shi X, Zhou B. The role of Nrf2 and MAPK pathways in PFOS-induced oxidative stress in zebrafish embryos. Toxicol Sci. 2010;115:391–400.
Lu Y, Gao K, Li X, Tang Z, Xiang L, Zhao H, et al. Mass spectrometry-based metabolomics reveals occupational exposure to Per- and Polyfluoroalkyl substances relates to oxidative stress, fatty acid beta-oxidation disorder, and kidney injury in a manufactory in China. Environ Sci Technol. 2019;53:9800–9.
Pierozan P, Cattani D, Karlsson O. Perfluorooctane sulfonate (PFOS) and perfluorooctanoic acid (PFOA) induce epigenetic alterations and promote human breast cell carcinogenesis in vitro. Arch Toxicol. 2020;94:3893–906.
White SS, Kato K, Jia LT, Basden BJ, Calafat AM, Hines EP, et al. Effects of perfluorooctanoic acid on mouse mammary gland development and differentiation resulting from cross-foster and restricted gestational exposures. Reprod Toxicol. 2009;27:289–98.
White SS, Stanko JP, Kato K, Calafat AM, Hines EP, Fenton SE. Gestational and chronic low-dose PFOA exposures and mammary gland growth and differentiation in three generations of CD-1 mice. Environ Health Perspect. 2011;119:1070–6.
Trudel D, Horowitz L, Wormuth M, Scheringer M, Cousins IT, Hungerbuhler K. Estimating consumer exposure to PFOS and PFOA. Risk Anal. 2008;28:251–69.
Haug LS, Salihovic S, Jogsten IE, Thomsen C, van Bavel B, Lindstrom G, et al. Levels in food and beverages and daily intake of perfluorinated compounds in Norway. Chemosphere. 2010;80:1137–43.
Sunderland EM, Hu XC, Dassuncao C, Tokranov AK, Wagner CC, Allen JG. A review of the pathways of human exposure to poly- and perfluoroalkyl substances (PFASs) and present understanding of health effects. J Expo Sci Environ Epidemiol. 2019;29:131–47.
Liang L, Pan Y, Bin L, Liu Y, Huang W, Li R, et al. Immunotoxicity mechanisms of perfluorinated compounds PFOA and PFOS. Chemosphere. 2022;291: 132892.
We express our appreciation to all women enrolled in our study for their participation and to the researchers who have contributed to the success of this research, including Jing Pang, Jiayu He, Honglu Zhang, and Hongyue Sun.
This work was supported by grants from the National Natural Science Foundation of China (grant numbers 81973002 and 81703192).
Ethics approval and consent to participate
All study procedures were approved by the Ethics Board of the Tianjin Medical University Cancer Institute and Hospital.
Consent for publication
All the authors have reviewed and approved the manuscript 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.
About this article
Cite this article
Li, X., Song, F., Liu, X. et al. Perfluoroalkyl substances (PFASs) as risk factors for breast cancer: a case–control study in Chinese population. Environ Health 21, 83 (2022). https://doi.org/10.1186/s12940-022-00895-3