Serum glycolipids mediate the relationship of urinary bisphenols with NAFLD: analysis of a population-based, cross-sectional study
Environmental Health volume 21, Article number: 124 (2022)
Bisphenol A (BPA) and its substitutes bisphenol S (BPS) and bisphenol F (BPF) are endocrine-disrupting chemicals widely used in consumer products, which have been proposed to induce various human diseases. In western countries, one of the most common liver diseases is non-alcoholic fatty liver disease (NAFLD). However, studies on the associations of the three bisphenols with NAFLD in human beings are scarce.
We included 960 participants aged ≥ 20 years from the NHANES 2013–16 who had available data on levels of urinary BPA, BPS and BPF. The hepatic steatosis index (HSI) > 36 was used to predict NAFLD. Logistic regression analysis and mediation effect analysis were used to evaluate the associations among bisphenols, glycolipid-related markers and NAFLD.
A total of 540 individuals (56.3%) were diagnosed with NAFLD, who had higher concentrations of BPA and BPS but not BPF than those without NAFLD. An increasing trend in NAFLD risks and HSI levels was observed among BPA and BPS tertiles (p for trend < 0.05). After adjustment for confounders, elevated levels of BPA or BPS but not BPF were significantly associated with NAFLD. The odds ratio for NAFLD was 1.581 (95% confidence intervals [CI]: 1.1–2.274, p = 0.013) comparing the highest with the lowest tertile of BPA and 1.799 (95%CI: 1.2462.597, p = 0.002) for BPS. Mediation effect analysis indicated that serum high-density lipoprotein cholesterol and glucose had a mediating effect on the relationships between bisphenols and NAFLD.
The present study showed that high exposure levels of BPA and BPS increased NAFLD incidence, which might be mediated through regulating glycolipids metabolism. Further studies on the role of bisphenols in NAFLD are warranted.
Non-alcoholic fatty liver disease (NAFLD), the most common chronic liver disease, has a global prevalence of approximately 25% among adults and is increasingly recognized as a hepatic manifestation part of metabolic syndrome . Furthermore, the economic burden of NAFLD is projected to increase during the coming decades . Growing evidence indicates that NAFLD frequently coexists with other abnormal conditions and might play a synergistic role in liver injury . Importantly, NAFLD patients commonly tend to have glycolipid metabolism disorders which are shared cardiometabolic risk factors for cardiovascular disease (CVD) . Additionally, NAFLD is obviously related to an increased risk of long-term morbidity and mortality mainly attributable to CVD . Although NAFLD is more prevalent in patients with obesity, insulin resistance, and diabetes, its cause remains largely unknown . In this regard, the importance of chemicals present in consumer products that are widely distributed in the environment is a relatively unexplored factor that could influence NAFLD [6, 7].
Bisphenol A (BPA) is an ingredient of polycarbonate plastics commonly found in food and beverage containers and used as an additive for other plastics . BPA exhibits multiple endocrine-disrupting effects on people of all ages . Considering the substantial potential human health risks based on its endocrine-disrupting feature and toxicological hazards [10, 11], many regulatory agencies have banned the use of BPA . Subsequently, the number of BPA-free products on the market have increased rapidly in the past two decades, and BPA analogues, including bisphenol F (BPF) and bisphenol S (BPS), have gradually been used in some BPA-free products [13, 14]. Consequently, the measured levels of BPS and BPF exposure in human samples markedly rise from 2000 to 2014 . However, the toxicological information of BPF and BPS is limited. Despite an increasing amount of evidence from in vitro and in vivo studies indicated that BPS and BPF, as BPA alternatives, might have similar endocrine-disrupting features as BPA, clinical studies and in vitro and in vivo studies on potential mechanisms for the association of BPS and BPF exposures with human health effects are rare .
Therefore, this study was performed to explore the relationships of BPA, BPS and BPF with NAFLD risks in a population-based study from the National Health and Nutrition Examination Survey (NHANES) and to evaluate the underlying role of glycolipid-related markers in these associations.
Materials and methods
Study design and participants
NHANES is a nationally cross-sectional survey of the non-institutionalized US residents performed by the National Center for Health Statistics (NCHS) at the US Centers for Disease Control and Prevention (CDC). A stratified, multistage probability cluster design was performed in NHANES and detailed information on demographics, socioeconomic data, personal lifestyle, individual medical and health conditions, and laboratory indicators were collected. The NHANES data are publicly released on a 2-year cycle.
The study sample was collected from 2013–14 and 2015–16 continuous cycles of the NHANES given that urinary of BPS and BPF levels were examined in these cycles. The flowchart of this study was shown in Fig. 1. The population in the final analysis consisted of participants aged 20 years and older with complete information, which were used to construct the parameter for calculating the hepatic steatosis index (HSI) to predict NAFLD prevalence, along with data on urinary BPA, BPS, and BPF. Furthermore, subjects who had viral hepatitis (hepatitis B and C) or with excessive alcohol consumption (> 3 drinks/day for males and > 2 drinks/day for females) and pregnant women were excluded. Additionally, we further excluded individuals with missing demographic and biochemical information on key covariates including poverty income ratio, education status, urine creatinine, glycosylated hemoglobin A1c (HbA1c), total cholesterol (TC), triglycerides (TG) and low-density lipoprotein cholesterol (LDL-C). A total of 960 participants were included in the final analysis.
The Institutional Review Board of the NCHS approved this program for NHANES and strictly implemented it in accordance with the Declaration of Helsinki . All participants signed Informed consent forms.
Measurement of urinary BPA, BPS and BPF concentrations
One-third of random subsamples aged six years or over were selected for determination of urinary bisphenols in NHANES. Urinary levels of BPA, BPS and BPF were measured by online solid-phase extraction coupled with high performance liquid chromatography and tandem mass spectrometry at the Laboratory Science Division of the National Center for Environmental Health of the CDC. This method was described in previous studies [15, 18]. According to the description of the NHANES website, when urinary levels of the three targets were below the lower limits of detection (LLOD), the value obtained by dividing the LLOD by the square root of 2 was routinely substituted (see https://wwwn.cdc.gov/Nchs/Nhanes/2013-2014/EPHPP_H.htm). The LLODs were 0.2 ng/mL for urinary BPA, 0.1 ng/mL for urinary BPS, and 0.2 ng/mL for urinary BPF. In total, 95% of BPA measurements, 55% of BPF measurements, and 90% of BPS measurements were above the LLODs in the current study. Additionally, urinary creatinine was included in all the analysis models of all three targets to correct for urinary dilution in the current study, as done in the previous study . All participants were divided into three subgroups according to tertiles of BPA and BPS, respectively. Because the detection rate of urinary BPF was relatively low, participants for whom the concentration of urinary BPF was not detected categorized as a reference (tertile 1). The remainder of the participants were divided into two groups (tertile 2 and 3) according to the median of urinary BPF value.
Detection of serum glycolipid-related parameters levels
Blood samples from every participant were collected in the non-fasting state. Serum TG and TC concentrations were detected using specific enzymatic assays. Immunoassays were applied to detect serum high-density lipoprotein cholesterol (HDL-C) concentrations . The values of serum LDL-C were calculated by the Friedewald equation was (LDL-C = TC – HDL-C – TG/5) .
Fasting glucose levels were determined by a hexokinase-method (Roche Diagnostics, Indianapolis, IN) . HbA1c levels were quantitatively measured using a Tosoh G8 Glycohemoglobin Analyzer (Tosoh Medics, Inc. San Francisco, CA) .
The gold standard for diagnosing NAFLD is liver biopsy, but it is not practical in the general population studies due to the invasiveness and high costs of the procedure. HSI is a common non-invasive screening tool to predict NAFLD presence and is calculated based on the following specific formula = 8 × alanine aminotransferase (ALT, IU/L)/aspartate aminotransferase (AST, IU/L) + body mass index (BMI, kg/m2) + 2 (if female) + 2 (if type 2 diabetes). Furthermore, participants were defined as having NAFLD when HSI > 36 according to a previous publication .
Other key covariates
Several potential confounders identified in previous literature [24, 25] were also included in this study as covariates to get closer to reality, including age, sex, ethnicity, family poverty income ratio (< 1, ≥ 1), education status, BMI (continuous), smoking status, drinking status, diabetes and hypertension status. Smoking status was classified into three subgroups according to never smoking before questionnaire (never), quit smoking at least one day (former) and now smoking or have smoked at least one hundred cigarettes (current). Individuals who had two drinks per day for men and one drink per day for women were considered to be drinking. The diagnosis of diabetes and hypertension of participants was based on the clinical criteria or depended on their medical history (self-reported or treated with antihypertensive/antidiabetic drugs) from questionnaires.
Continuous indices were presented as the means (standard deviation) or medians (interquartile range) and categorical indices were expressed as numbers (percentage). The distribution pattern of variables was tested by Kolmogorov–Smirnov analysis. Because TG, BPA, BPS, and BPF levels exhibited skewed distributions, log-transformed levels of the four indicators as continuous variables were applied in our analysis. Student’s t test, analysis of variance, the Mann–Whitney U test, or the λ2 test were used to evaluate significant differences in variables between groups where appropriate. Pearson correlation analysis was used to examine the correlations of BPA, BPS and BPF with glycolipid-related parameters and HSI. Logistic regression analysis was applied to examine the odds ratios (ORs) with 95% confidence intervals (CIs) for NAFLD according to tertiles and log-transformed continuous concentrations of urinary BPA, BPS and BPF. Urine creatinine was included in Model 1. Furthermore, ethnicity, education status, drinking, hypertension, diabetes, TG, HDL-C, glucose, HbA1c and urine creatinine were included in model 2. Additionally, model 3 was also adjusted for BPA and/or BPS (log-transformed concentrations) based on model 2 in different conditions. Considering the possible effects of bisphenols on metabolism in participants of different sexes, we performed stratified analyses by sex.
Mediation effect analysis was utilized to evaluate an intermediate variable as a mediator in the pathway between a risk factor and an outcome, estimating the extent to which the effect of the risk factor occurred through the mediator. Mediation effect analysis was performed according to three models : path c for explaining the outcome (Y, NAFLD/HSI) of a risk factor (X, bisphenols) using a logistic regression model (Y is NAFLD) or linear regression model (Y is HSI), path a for explaining the mediator (M, glycolipid-related indicators) for the risk factor using a linear regression model, and path b for estimating the association of the mediator and the risk factor with the outcome using a logistic regression model (Y is NAFLD) or linear regression model (Y is HSI). In addition, mediation effect analysis was performed with adjustments for urine creatinine, ethnicity, education status, hypertension and diabetes (or without diabetes in the mediation effect analysis models of glucose and HbA1c) in the current study. The total effects of bisphenols consisted of the direct effect (DE) and the indirect effect (IE) on NAFLD or HSI. When propositions of a “statistically significant association between bisphenols (X) and mediators (M)” and a “statistically significant association of mediators with NAFLD or HSI (Y)” were true, the intermediate effect occurred as previous description .
To examine the robustness of these findings, we performed sensitivity analyses by excluding individuals with levels below the LLOD of urinary BPF because the detection rate of BPF was relatively low. SPSS (version 25.0) and R (version 4.1.3) were used for all statistical analyses in the present study. Here, p < 0.05 was suggested statistically significant difference.
In total, 960 adults aged 20 years and older were finally eligible to include in the present study, of whom 444 (46.2%) were male. The mean age of the whole population was 52.9 ± 17.3 years. Among them, 540 individuals, over 50% of the total number, were defined as NAFLD according to HSI > 36. The mean HSI values were 38.5 ± 8.2, 43.9 ± 6.65 and 31.5 ± 3.07 in the total population, the NAFLD and non-NAFLD groups, respectively. The baseline characteristics classified by NAFLD prevalence were presented in Table 1. The interquartile ranges of urinary levels were from 0.5 to 2.2 ng/mL with a median of 1.1 ng/mL for BPA, from 0.2 to 1.2 ng/mL with a median of 0.5 ng/mL for BPS, and from 0.14 to 0.8 ng/mL with a median of 0.25 ng/mL for BPF. Participants with NAFLD were more prone to be non-Hispanic black (25.7%), received less than a high school education (25%), and suffering from hypertension (54.8%) and diabetes (34.1%), and were prevalent in higher BMI, TG, glucose and HbA1c levels and lower of HDL-C levels compared to those without NAFLD (all p < 0.05). Moreover, individuals with NAFLD had elevated urinary concentrations of BPA and BPS compared with those without NAFLD [BPA: 1.25 (0.6, 2.32) vs 0.9 (0.5, 2) ng/mL, p = 0.001; BPS: 0.6 (0.2, 1.33) vs 0.4 (0.2, 0.95) ng/mL, p < 0.001]. However, no significant differences with regard to age, sex, family poverty income ratio, smoking status, TC, LDL-C and BPF levels were noted between the NAFLD and non-NAFLD groups (all p > 0.05). The demographic characteristics were similar between individuals with NAFLD and without NAFLD in females and males with the exception of education status (shown in Supplementary Table S1). Moreover, elevated urinary concentrations of BPA and BPS but not BPF were observed in NAFLD group compared with non-NAFLD groups in both females and males (Supplementary Table S1). These findings were consistent with the analysis results for the total population. Additionally, only urinary BPA levels were higher in females than males among the three bisphenols [BPA: 1.15 (0.5, 2.6) vs 1(0.5, 2) ng/mL, p = 0.025, Supplementary Table S2).
Additionally, as shown in Table 2, we divided the whole population into three subgroups according to the tertile concentrations of urinary BPA and BPS and the levels of BPF. Participants with the highest tertile of BPA were more likely to be non-Hispanic white, have elevated levels of BMI, BPS and BPF, and have lower levels of HDL-C than those with the lowest tertile of BPA (all p < 0.05). Younger non-Hispanic blacks and subjects with a lower poverty income ratio and higher BMI, glucose, HbA1c and BPA levels tended to have elevated urinary concentrations of BPS (all p < 0.05). The following features were noted in high BPF levels group: non-Hispanic Whites, current smokers and drinkers, higher levels of education, higher BMI, BPA and BPS levels, and lower concentrations of TC and LDL-C were prevalent in the high level of urinary BPF group (all p < 0.05). Of note, ascending trends in NAFLD incident rates and HSI levels were observed in individuals who had higher urinary levels of BPA and BPS (all p < 0.05). However, NAFLD incidence and HSI levels were not significantly different among BPF tertiles.
Associations of urinary bisphenols with NAFLD
Logistic regression models were applied to analyze the associations between the three bisphenols and NAFLD risks, as summarized in Table 3. When only adjusting for urine creatine levels in model 1, elevated urinary levels of BPA and BPS, but not BPF, were significantly related to the risk of NAFLD (all p < 0.05). Because BMI was the one of the key parts of the special equation to calculate the HSI, it was not possible to adjust BMI in multivariate analysis for the NAFLD prediction model based on the HSI, as previously recommended . In model 2, some key covariates including ethnicity, education status, drinking status, hypertension, diabetes, TG, HDL-C, glucose and HbA1c, were adjusted based on model 1. After adjustment for these covariates in model 2, comparing the highest with lowest tertile of bisphenols, the OR with 95%CI for NAFLD was 1.786 (1.255–2.542, p = 0.001) for BPA, 1.97 (1.381–2.811, p < 0.001) for BPS and 1.344 (0.937–1.929, p = 0.108) for BPF. In model 3, urinary of BPA and/or BPS levels (as log-transformation) were further mutually adjusted under different conditions based on model 2. Participants in the highest tertile of BPA and BPS had 58.1% (OR: 1.581, 95%CI: 1.1–2.274) and 79.9% (OR: 1.799, 95%CI: 1.246–2.597) elevated risks for NAFLD, respectively, compared to those in the lowest tertile subgroups. As continuous variables, per log unit increases in levels of BPA and BPS were significantly related to 1.391-fold (95%CI: 1.007–1.921) and 1.476-fold (95%CI: 1.139–1.913) risks for NAFLD, respectively, after adjusting for all multiple confounding factors. No significant relationship was observed between BPF and NAFLD regardless of whether BPF was included as a categorical or continuous variable in logistic regression models.
Furthermore, the gender stratification results of associations between three bisphenols and NAFLD risks were presented in Supplementary Table S3. There was no significant interaction between sex and three bisphenols levels (all p for interaction > 0.05), although the relationship between BPS and NAFLD seemed stronger in females than males. Comparing the highest with the lowest tertile of urinary BPS, the OR for NAFLD was 2.304 (95%CI:1.394–3.807, p = 0.002 for trend) in females and 1.623 (95%CI: 0.979–2.691, p = 0.054 for trend) in males after adjusting for multiple confounders. The associations of BPF and BPA with NAFLD did not differ by sex (Supplementary Table S3). In multivariable adjusted models, no significant effect of urinary BPF levels on NAFLD and an obvious effect of urinary BPA levels on NAFLD were observed both in women and men (Supplementary Table S3), which were consistent with that obtained for the total population.
Additionally, as shown in Supplementary Table S4, the associations of urinary bisphenols levels except BPA with NAFLD did not change appreciably in sensitivity analyses when excluding individuals with urinary BPF levels below the LLOD.
Equally important, Pearson correlation analysis showed that log-transformed urinary BPA and BPS levels were significantly and positively correlated with HSI (r = 0.1 between BPA and HSI, r = 0.11 between BPS and HSI, and r = 0.07 between BPF and HSI, all p < 0.05), as presented in Fig. 2.
Roles of glycolipid-related indices in the relationships of urinary bisphenols with NAFLD
According to Pearson correlation analysis, the correlation coefficients among urinary bisphenols, serum glucolipid-related indices and HSI were presented in Fig. 2. A significantly negative correlation of BPA with HDL (r = -0.1, p < 0.05) and a significantly positive association of BPS with glucose (r = 0.07, p < 0.05) were found in the current study. Furthermore, HSI levels were remarkably and positively correlated with TG, glucose and HbA1c levels and negatively correlated with HDL-C levels in the entire population (r = 0.27 for TG, r = 0.24 for glucose, r = 0.23 for HbA1c and r = -0.31 for HDL-C, all p < 0.05, Fig. 2). Because no significant association between BPF and NAFLD/HSI was observed, mediation effect analysis was performed to identify whether glucolipid-related indicators mediated the efects of BPA and BPS on NAFLD/HSI. The role of glucolipid-related indicators levels in the associations of bisphenols with NAFLD was explored in a model with NAFLD/HSI as the dependent variable, urinary bisphenols levels as the independent variable, and glucolipid-related indicators levels as the mediator variable, and urine creatinine, ethnicity, education status, hypertension and diabetes (or without diabetes where appropriate) as covariates. Interestingly, the results presented in Figs. 3 and 4 revealed that bisphenols may participate in the development of NAFLD by mediating serum HDL-C and glucose levels. Specifically, the mediated efficacy of HDL-C accounted for 24.5% in the relationship between BPA and NAFLD development (IE = 0.025, 95%CI: 0.007–0.046, p < 0.05, Fig. 3), whereas the mediated efficacy of glucose accounted for 15.2% of the predictive ability of BPS on NAFLD (IE = 0.014, 95%CI: 0.001–0.029, p < 0.05, Fig. 4). In addition, we further analyzed the mediated effects of glucose-mediated indices on the relationships of BPA and BPS with HSI (Fig. 5). Similarly, HDL-C and glucose had significant mediation effects on the associations between bisphenols and HSI. The estimate of IE between BPA and HSI through HDL-C was 0.425 (95%CI: 0.106–0.763, p < 0.05), with a 23.1% of mediated efficacy. Moreover, the mediation efficacy of glucose accounted for 12.1% of the relation between BPS and HSI (IE = 0.174, 95%CI: 0.014–0.362, p < 0.05).
This study explored the associations of multiple bisphenols including BPA, BPS and BPF with NAFLD prevalence and first investigated the role of glycolipid-related factors in these relationships. In this US population-based study, our results revealed that adults with elevated concentrations of urinary BPA and BPS, but not BPF, tended to develop NAFLD after adjusting for relevant confounders, including sociodemographic variables, lifestyle factors, clinical and laboratory parameters. More importantly, mediation effect analysis indicated that serum glycolipid-related indices, such as HDL-C and glucose, potentially exhibit mediated effects on the associations of bisphenols with the prevalence of NAFLD. These findings further provided critical information with regard to the potential relations between bisphenols and NAFLD and insights into the underlying mechanisms of bisphenols in NAFLD development.
In the past decades, BPA had been used in various industrial applications, especially for the production of polycarbonate and plastic epoxy . Currently, BPA has become one of the highest produced volumes among chemicals worldwide . Given the numerous applications of BPA, people have expressed concerns regarding whether it is harmful to human health. The wealth of studies have proven toxicological properties and endocrine-disrupting characteristics of BPA on human diseases, such as CVD and metabolic disorders [28, 29]. When BPA was restricted in numerous consumer products since 2016, other bisphenols to replace BPA were emerging in industrial production . Subsequently, BPS and BPF are considered suitable substitutes for BPA and are widely used in industrial applications . Although BPS and BPF are thought to be endocrine-disrupting and toxic chemicals given their similar chemical structures to BPA , less information is available regarding the role of BPS and BPF in adverse human health conditions. Notably, human exposure to these three bisphenols is widespread by a variety of ways in the general population , so biomonitoring of these chemicals and studying on their relationships and mechanisms in human diseases are a real concern.
NAFLD has an estimated prevalence of 25% in adults worldwide, with a continuously increasing trend [30, 31]. Unfortunately, during the coming decades, the economic burden of NAFLD will likely increase rapidly due to the high disease prevalence and complications [31, 32]; thus, its underlying influencing factors should be identified. Emerging evidence has revealed the detrimental role of environmental toxins in liver injury, particularly in the incidence and progression of NAFLD . Exposure to bisphenols, especially BPA, as a ubiquitous endocrine-disrupting and toxic chemical, has also been reported to be associated with obesity, diabetes, and CVD [24, 29, 34]. In addition, data on the effects of BPA on NAFLD are increasing [27, 35]. Several cross-sectional studies revealed that urinary concentrations of BPA were significantly related to NAFLD in Hispanic adolescents and US adults [27, 35]. However, population-based studies exploring the relationships of BPS and BPF with NAFLD are scarce. In the present study, our findings also supported previous results of a positive association of urinary BPA concentrations and NAFLD prevalence [27, 35]. Moreover, our study first investigated the relationships of BPA alternatives, including BPS and BPF, with NAFLD risks, and a significant and positive relationship with NAFLD was observed for BPS but not BPF in this cohort of 960 adults, further providing more information on the contribution of bisphenols to the occurrence of NAFLD. In addition, we found a stronger association between BPS and NAFLD in women than men, although there was no difference in urinary BPS concentration between men and women in the current study. The estrogenic activities of BPS may explain such phenomenon but the underlying mechanisms remain to be understood.
Of note, pathogenic pathways of NAFLD are influenced by multiple metabolic factors, including abnormal glucose metabolism (due to insulin resistance) and lipid abnormalities (for example, low HDL-C and high TG) . Moreover, given its endocrine-disrupting characteristic, BPA is thought to participate in de novo fatty acid synthesis , stimulate the accumulation of TG in adipocytes and human hepatocellular carcinoma cells  and up-regulate the expression of genes involved in lipid metabolism . Similar to BPA, BPF and BPS are considered as endocrine-disrupting chemicals and have some BPA like estrogenic activities . Furthermore, it has been proposed that BPS and BPF enhance lipid accumulation and promote the differentiation of preadipocytes into adipocytes through the PPARγ signaling pathway [38,39,40]. Additionally, the exposure to BPA and BPS resulted in the differential expression of genes related to lipid metabolism and adipogenesis and also increasing 3T3-L1 adipocyte differentiation [41,42,43]. Notably, BPA can disturb glucose homeostasis and influence the pancreatic β-cell function through multiple mechanisms . BPS also plays critical roles in disturbance of glucose metabolism . Although dysglycemia and lipid abnormalities appear to bridge the link between bisphenols and NAFLD based on the above data, whether glycolipid metabolism disorders affect the relationship between bisphenols and the occurrence and development of NAFLD remains unknown. Therefore, the current study first investigated the potential role of serum glycolipid-related parameters in the association of bisphenols and NAFLD using mediation effect analyses. Interestingly, a significant mediated efficacy for HDL-C in the relation of BPA with NAFLD and for serum glucose in the relation of BPS with NAFLD were found, indicating abnormal glucolipid metabolism may mediate the associations between bisphenols NAFLD. Further studies are needed to confirm these findings and assess possible glucolipid metabolism mechanisms explaining the effects of bisphenols on NAFLD.
Due to the relatively high percentage of participants with urinary BPF levels below the LLOD, we did sensitivity analyses in those with BPF levels above the LLOD and found that the significant association of BPA with NAFLD disappeared. Here, the relatively small sample size of sensitivity analyses (N = 525) might explain this result. Regarding the negative result between BPF and NAFLD, we had to consider the influence of that the relatively high percentage of participants with urinary BPF levels below the LLOD, and the exposure levels of BPF were approximately one-fifth of BPA levels and one-half of BPS levels in the whole population. Therefore, although we did not observe a positive association between BPF and NAFLD in the current study, attention should still be given to BPF exposure levels and its harmful effects on human health. Future prospective cohort studies including participants exposed to high concentrations of urinary BPF need to be conducted.
This study analyzed a nationally representative data from the NHANES study, which can avoid the effect on the generalizability of the results and enable our findings to be generalized to a broader population. In addition, we, for the first time, observed the relationships among three target bisphenols, NAFLD and serum glucolipid-related indicators in a population.
Nevertheless, several limitations should be noted in this study. First, given the inherent features of the cross-sectional and observational study, only the relationship between bisphenols and NAFLD as well as glycolipid-related parameters, could be clarified, the causal inference could not be assessed. Second, in the current study, after exclusion of other causes of chronic liver disease, the non-invasive index (HSI) was used to define NAFLD, which can inevitably lead to misdiagnosis of the presence or absence of NAFLD. Even though the diagnostic golden standard for NAFLD is liver biopsy, it is not applicable in large population-based studies. Moreover, the HSI, which includes important liver enzyme indicators (ALT and AST), has been recognized as a reliable non-invasive diagnostic method for NAFLD in populations . The exact correlation between bisphenols exposure and the development of NAFLD should be demonstrated in the future. Third, although we included as many confounding factors as possible based on NHANES data, the potential influences of unmeasured confounders on our study cannot be ruled out. Fourth, the proportion of subjects with urinary BPF concentrations below the LLOD was relatively high. Although this did not appreciably affect the categorization of BPF exposure in our categorical data analysis, we further performed sensitivity analyses in the population with concentrations of urinary BPF above the LLOD. Fifth, NHANES staff collected spot urine samples of participants to detect levels of bisphenols, which could bring a considerable intra-individual and inter-person temporal variability of urinary bisphenols levels. However, the population from the NHANES study is sufficiently large and urine samples of every participant are randomly collected relative to eating time and bladder-emptying time, therefore, the urinary bisphenols concentration based on the analysis of one spot urine sample can fully reflect the average exposure of the population to bisphenols as noted in a previous study . Finally, information about sources of BPA, BPS and BPF is limited in NHANES, and individual confounding patterns of the three exposures could not be determined. Further studies need to confirm these associations.
In this population-based, cross-sectional study, we found that urinary BPA and BPS levels were significantly associated with NAFLD risks and that serum HDL-C and glucose had obvious mediation effects on these relationships, suggesting that biomonitoring and avoiding daily exposure to bisphenols represent possible strategies for the prevention of NAFLD.
Availability of data and materials
The data for this study were obtained from the public database of the National Center for Health Statistics. https:// www.cdc.gov/nchs/nhanes/.
Non-alcoholic fatty liver disease
National Health and Nutrition Examination Survey
Hepatic steatosis index
Low-density lipoprotein cholesterol
Glycosylated hemoglobin A1c
Lower limits of detection
High-density lipoprotein cholesterol
Body mass index
Younossi ZM, Koenig AB, Abdelatif D, Fazel Y, Henry L, Wymer M. Global epidemiology of nonalcoholic fatty liver disease-meta-analytic assessment of prevalence, incidence, and outcomes. Hepatology. 2016;64(1):73–84.
Powell EE, Wong VW, Rinella M. Non-alcoholic fatty liver disease. Lancet. 2021;397(10290):2212–24.
Powell EE, Jonsson JR, Clouston AD. Steatosis: co-factor in other liver diseases. Hepatology. 2005;42(1):5–13.
Vernon G, Baranova A, Younossi ZM. Systematic review: the epidemiology and natural history of non-alcoholic fatty liver disease and non-alcoholic steatohepatitis in adults. Aliment Pharmacol Ther. 2011;34(3):274–85.
EASL-EASD-EASO Clinical Practice Guidelines for the management of non-alcoholic fatty liver disease. J Hepatol 2016, 64(6):1388–1402.
Dallio M, Masarone M, Errico S, Gravina AG, Nicolucci C, Di Sarno R, Gionti L, Tuccillo C, Persico M, Stiuso P, et al. Role of bisphenol A as environmental factor in the promotion of non-alcoholic fatty liver disease: in vitro and clinical study. Aliment Pharmacol Ther. 2018;47(6):826–37.
Foulds CE, Treviño LS, York B, Walker CL. Endocrine-disrupting chemicals and fatty liver disease. Nat Rev Endocrinol. 2017;13(8):445–57.
Vandenberg LN, Hauser R, Marcus M, Olea N, Welshons WV. Human exposure to bisphenol A (BPA). Reprod Toxicol. 2007;24(2):139–77.
Lee J, Ahn YA, Choi K, Park J, Moon HB, Choi G, Lee JJ, Suh E, Kim HJ, Eun SH, et al. Bisphenol A in infant urine and baby-food samples among 9- to 15-month-olds. Sci Total Environ. 2019;697: 133861.
Zhang YF, Ren XM, Li YY, Yao XF, Li CH, Qin ZF, Guo LH. Bisphenol A alternatives bisphenol S and bisphenol F interfere with thyroid hormone signaling pathway in vitro and in vivo. Environ Pollut. 2018;237:1072–9.
Wei Y, Han C, Li S, Cui Y, Bao Y, Shi W. Maternal exposure to bisphenol A during pregnancy interferes ovaries development of F1 female mice. Theriogenology. 2020;142:138–48.
Bousoumah R, Leso V, Iavicoli I, Huuskonen P, Viegas S, Porras SP, Santonen T, Frery N, Robert A, Ndaw S. Biomonitoring of occupational exposure to bisphenol A, bisphenol S and bisphenol F: a systematic review. Sci Total Environ. 2021;783: 146905.
Liao C, Liu F, Kannan K. Bisphenol s, a new bisphenol analogue, in paper products and currency bills and its association with bisphenol a residues. Environ Sci Technol. 2012;46(12):6515–22.
Liao C, Kannan K. A survey of alkylphenols, bisphenols, and triclosan in personal care products from China and the United States. Arch Environ Contam Toxicol. 2014;67(1):50–9.
Ye X, Wong LY, Kramer J, Zhou X, Jia T, Calafat AM. Urinary concentrations of Bisphenol A and three other Bisphenols in convenience samples of U.S. adults during 2000–2014. Environ Sci Technol. 2015;49(19):11834–9.
Rochester JR, Bolden AL. Bisphenol S and F: a systematic review and comparison of the hormonal activity of bisphenol a substitutes. Environ Health Perspect. 2015;123(7):643–50.
Mao W, Hu Q, Chen S, Chen Y, Luo M, Zhang Z, Geng J, Wu J, Xu B, Chen M. Polyfluoroalkyl chemicals and the risk of kidney stones in US adults: a population-based study. Ecotoxicol Environ Saf. 2021;208: 111497.
Zhou X, Kramer JP, Calafat AM, Ye X. Automated on-line column-switching high performance liquid chromatography isotope dilution tandem mass spectrometry method for the quantification of bisphenol A, bisphenol F, bisphenol S, and 11 other phenols in urine. J Chromatogr B Analyt Technol Biomed Life Sci. 2014;944:152–6.
Barr DB, Wilder LC, Caudill SP, Gonzalez AJ, Needham LL, Pirkle JL. Urinary creatinine concentrations in the U.S. population: implications for urinary biologic monitoring measurements. Environ Health Perspect. 2005;113(2):192–200.
Perak AM, Ning H, Kit BK, de Ferranti SD, Van Horn LV, Wilkins JT, Lloyd-Jones DM. Trends in Levels of Lipids and Apolipoprotein B in US Youths Aged 6 to 19 Years, 1999–2016. JAMA. 2019;321(19):1895–905.
Mehta R, Reyes-Rodríguez E, Yaxmehen Bello-Chavolla O, Guerrero-Díaz AC, Vargas-Vázquez A, Cruz-Bautista I. C AA-S: Performance of LDL-C calculated with Martin’s formula compared to the Friedewald equation in familial combined hyperlipidemia. Atherosclerosis. 2018;277:204–10.
Menke A, Casagrande S, Cowie CC. Contributions of A1c, fasting plasma glucose, and 2-hour plasma glucose to prediabetes prevalence: NHANES 2011–2014. Ann Epidemiol. 2018;28(10):681-685.e682.
Lee JH, Kim D, Kim HJ, Lee CH, Yang JI, Kim W, Kim YJ, Yoon JH, Cho SH, Sung MW, et al. Hepatic steatosis index: a simple screening tool reflecting nonalcoholic fatty liver disease. Dig Liver Dis. 2010;42(7):503–8.
Liu B, Lehmler HJ, Sun Y, Xu G, Liu Y, Zong G, Sun Q, Hu FB, Wallace RB, Bao W. Bisphenol A substitutes and obesity in US adults: analysis of a population-based, cross-sectional study. Lancet Planet Health. 2017;1(3):e114–22.
Hu Z, Li Y, Yang Y, Yu W, Xie W, Song G, Qian Y, Song Z. Serum lipids mediate the relationship of multiple polyaromatic hydrocarbons on non-alcoholic fatty liver disease: A population-based study. Sci Total Environ. 2021;780: 146563.
Baron RM, Kenny DA. The moderator-mediator variable distinction in social psychological research: conceptual, strategic, and statistical considerations. J Pers Soc Psychol. 1986;51(6):1173–82.
Kim D, Yoo ER, Li AA, Cholankeril G, Tighe SP, Kim W, Harrison SA, Ahmed A. Elevated urinary bisphenol A levels are associated with non-alcoholic fatty liver disease among adults in the United States. Liver Int. 2019;39(7):1335–42.
Zhang YF, Shan C, Wang Y, Qian LL, Jia DD, Zhang YF, Hao XD, Xu HM. Cardiovascular toxicity and mechanism of bisphenol A and emerging risk of bisphenol S. Sci Total Environ. 2020;723: 137952.
Rancière F, Lyons JG, Loh VH, Botton J, Galloway T, Wang T, Shaw JE, Magliano DJ. Bisphenol A and the risk of cardiometabolic disorders: a systematic review with meta-analysis of the epidemiological evidence. Environ Health. 2015;14:46.
Younossi Z, Anstee QM, Marietti M, Hardy T, Henry L, Eslam M, George J, Bugianesi E. Global burden of NAFLD and NASH: trends, predictions, risk factors and prevention. Nat Rev Gastroenterol Hepatol. 2018;15(1):11–20.
Younossi ZM, Blissett D, Blissett R, Henry L, Stepanova M, Younossi Y, Racila A, Hunt S, Beckerman R. The economic and clinical burden of nonalcoholic fatty liver disease in the United States and Europe. Hepatology. 2016;64(5):1577–86.
Estes C, Anstee QM, Arias-Loste MT, Bantel H, Bellentani S, Caballeria J, Colombo M, Craxi A, Crespo J, Day CP, et al. Modeling NAFLD disease burden in China, France, Germany, Italy, Japan, Spain, United Kingdom, and United States for the period 2016–2030. J Hepatol. 2018;69(4):896–904.
Rajak S, Raza S, Tewari A, Sinha RA. Environmental Toxicants and NAFLD: A Neglected yet Significant Relationship. Dig Dis Sci. 2022;67(8):3497–507.
Lang IA, Galloway TS, Scarlett A, Henley WE, Depledge M, Wallace RB, Melzer D. Association of urinary bisphenol A concentration with medical disorders and laboratory abnormalities in adults. JAMA. 2008;300(11):1303–10.
Verstraete SG, Wojcicki JM, Perito ER, Rosenthal P. Bisphenol a increases risk for presumed non-alcoholic fatty liver disease in Hispanic adolescents in NHANES 2003–2010. Environ Health. 2018;17(1):12.
Marmugi A, Ducheix S, Lasserre F, Polizzi A, Paris A, Priymenko N, Bertrand-Michel J, Pineau T, Guillou H, Martin PG, et al. Low doses of bisphenol A induce gene expression related to lipid synthesis and trigger triglyceride accumulation in adult mouse liver. Hepatology. 2012;55(2):395–407.
Wada K, Sakamoto H, Nishikawa K, Sakuma S, Nakajima A, Fujimoto Y, Kamisaki Y. Life style-related diseases of the digestive system: endocrine disruptors stimulate lipid accumulation in target cells related to metabolic syndrome. J Pharmacol Sci. 2007;105(2):133–7.
Boucher JG, Ahmed S, Atlas E. Bisphenol S Induces Adipogenesis in Primary Human Preadipocytes From Female Donors. Endocrinology. 2016;157(4):1397–407.
Zheng S, Shi JC, Hu JY, Hu WX, Zhang J, Shao B. Chlorination of bisphenol F and the estrogenic and peroxisome proliferator-activated receptor gamma effects of its disinfection byproducts. Water Res. 2016;107:1–10.
Janesick A, Blumberg B. Minireview: PPARγ as the target of obesogens. J Steroid Biochem Mol Biol. 2011;127(1–2):4–8.
Boucher JG, Gagné R, Rowan-Carroll A, Boudreau A, Yauk CL, Atlas E. Bisphenol A and Bisphenol S Induce Distinct Transcriptional Profiles in Differentiating Human Primary Preadipocytes. PLoS ONE. 2016;11(9): e0163318.
Ahmed S, Atlas E. Bisphenol S- and bisphenol A-induced adipogenesis of murine preadipocytes occurs through direct peroxisome proliferator-activated receptor gamma activation. Int J Obes (Lond). 2016;40(10):1566–73.
Fang R, Yang S, Gu X, Li C, Bi N, Wang HL. Early-life exposure to bisphenol A induces dysregulation of lipid homeostasis by the upregulation of SCD1 in male mice. Environ Pollut. 2022;304: 119201.
Farrugia F, Aquilina A, Vassallo J, Pace NP. Bisphenol A and Type 2 Diabetes Mellitus: a review of Epidemiologic, functional, and early life factors. Int J Environ Res Public Health. 2021;18(2):716.
Zhao F, Jiang G, Wei P, Wang H, Ru S. Bisphenol S exposure impairs glucose homeostasis in male zebrafish (Danio rerio). Ecotoxicol Environ Saf. 2018;147:794–802.
Unalp-Arida A, Ruhl CE. Liver fibrosis scores predict liver disease mortality in the United States population. Hepatology. 2017;66(1):84–95.
Ye X, Wong LY, Bishop AM, Calafat AM. Variability of urinary concentrations of bisphenol A in spot samples, first morning voids, and 24-hour collections. Environ Health Perspect. 2011;119(7):983–8.
The authors thanked NCHS of CDC for sharing the data.
The National Natural Science Foundation of China (No.8197021705) and the Hunan Provincial Natural Science Foundation of China (No.2022JJ30069) awarded to Dr. Qin-Lin Ma, MD, PhD were partly supported this study.
Ethics approval and consent to participate
Exempt approval was obtained from the Committee on Human Research at the University of California, San Francisco.
Consent for publication
No conflict of interest with respect to this manuscript to disclose.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Baseline characteristicsof the study population according to NAFLD and non-NAFLD stratified by sex, NHANES2013–16. Table S2. Urinary concentrations of BPA, BPF and BPS according to sex. Table S3. Associations of urinary BPA, BPS, and BPF levels with NAFLD stratified by sex. Table S4. Associations of urinary BPA, BPS, and BPF levels with NAFLD in participants with urinary BPF concentrations above the LLOD (N=525).
About this article
Cite this article
Peng, J., Du, LL. & Ma, QL. Serum glycolipids mediate the relationship of urinary bisphenols with NAFLD: analysis of a population-based, cross-sectional study. Environ Health 21, 124 (2022). https://doi.org/10.1186/s12940-022-00945-w