Skip to main content

Association between exposure to brominated flame retardants (BFRs) and blood lipid profiles in American adults: a cross-sectional study

Abstract

Background

Exposure to brominated flame retardants (BFRs) has been linked to alterations in human metabolism and disease processes. However, the relationship between BFR exposure and blood lipid levels remains unclear. This study aimed to investigate the potential association between BFR exposure and blood lipid profiles in American adults.

Methods

A cross-sectional study was conducted using data from the National Health and Nutrition Examination Survey (NHANES) 2005–2016. Serum concentrations of twelve BFRs, PBB153 and eleven polybrominated diphenyl ethers (PBDEs), were quantified using isotope dilution gas chromatography/high-resolution mass spectrometry (GC/HRMS). Blood lipid levels, including total cholesterol (TC) and high-density lipoprotein cholesterol (HDL-C) were measured enzymatically. The Friedewald equation was used to determine low-density lipoprotein cholesterol (LDL-C): [LDL-C] = [TC] − [HDL-C] − [TG/5]. Remnant cholesterol (RC) was calculated using the formula: [RC] = [TC] − [HDL-C] − [LDL-C]. Multivariable regression analyses were applied to examine the associations between individual BFRs and TC, HDL-C, LDL-C, and RC. The overall associations of BFR mixtures with blood lipids were evaluated using quantile g-computation (QGC) analyses and weighted quantile sum (WQS) regression. In order to identify potential gender-specific differences, stratified mixture analyses were performed by gender.

Results

A total of 3,154 eligible participants were included. Nine BFRs with a detection rate greater than 70% were included in the analysis. Individually, PBB153, PBDE209, PBDE153, and PBDE28 were positively associated with TC and RC after adjusted all covariates. Furthermore, PBB153, PBDE209, and PBDE153 were positively associated with LDL-C. No association was found between individual BFR and HDL-C. WQS and QGC analyses confirmed that BFR mixtures were positively associated with TC, LDL-C, and RC.

Conclusion

This study demonstrates that BFR exposure is associated with increased levels of TC, LDL-C, and RC, indicating an elevated risk of dyslipidemia and cardiovascular diseases.

Introduction

Dyslipidemia has a substantial financial impact worldwide and is a risk factor for atherosclerotic cardiovascular disease (ASCVD), the world's leading cause of death​ [1, 2]. High-density lipoprotein (HDL), the smallest lipoprotein in circulation, consists of proteins and lipids organized into distinct subspecies. It contains a hydrophobic core of neutral lipids surrounded by a monolayer of amphipathic lipids, and serves various physiological functions​ [3, 4]. Numerous studies conducted over the decades have consistently shown an inverse relationship between high-density lipoprotein cholesterol (HDL-C) levels and the risk of coronary heart disease (CHD) [5,6,7]. Additionally, low HDL-C levels linked to a range of non-cardiovascular conditions, including cancer, infectious diseases, and autoimmune disorders [8,9,10]. Conversely, it is well-established that reducing low-density lipoprotein cholesterol (LDL-C) levels lowers the incidence of CHD, with elevated LDL-C recognized as an independent risk factor for CHD [11].

Remnant cholesterol (RC) is an emerging lipoprotein marker that has shown promise in predicting the risk of ischemic stroke, ASCVD, and mortality [12,13,14]. RC is composed of triglyceride-rich lipoproteins, including chylomicron remnants, intermediate-density lipoprotein cholesterol (IDL-C), and very low-density lipoprotein cholesterol (VLDL-C). It can be computed by deducting LDL-C and HDL-C from total cholesterol (TC) using lipid profile data​ [15, 16]. A nationwide cohort study revealed that individuals with RC levels of 27.7 mg/dL or higher had increased hazard ratios (HRs) for all-cause mortality (HR = 1.03), ischemic stroke mortality (HR = 1.22), ischemic heart disease mortality (HR = 1.19), and cardiovascular disease mortality (HR = 1.17) compared to those with RC levels below 17.9 mg/dL [17].

Over recent decades, the widespread use of chemical flame retardants has played a key role in reducing fire-related incidents​ [18]. There are many different types of flame retardants, including the inorganic flame retardants, halogenated organic, phosphorus-containing, and nitrogen-containing [19]. Brominated flame retardants (BFRs) has occupied the main position due to their superior performance and low cost [20]. While traditional BFRs like polybrominated diphenyl ethers (PBDEs) and polybrominated biphenyls (PBBs) have enhanced fire safety, their environmental persistence and bioaccumulative potential have raised substantial health concerns [21,22,23]. Exposure to BFRs—whether through inhalation, dermal contact, or ingestion—has been linked to various adverse health outcomes [24]. Previous studies have suggested associations between BFR exposure and several negative health effects, including neurotoxicity, metabolic disorders, periodontitis, chronic kidney disease, hypertension, and osteoporosis​ [25,26,27,28,29]. Additionally, prenatal exposure to BFRs has been associated with neurodevelopmental and metabolic disturbances​ [30, 31].

Although studies on the health effects of BFR exposure are increasing, the relationship with lipid levels remains unclear. Some studies have found that BFR exposure may be associated with decreased lipid levels. For example, PBDE209 and other PBDE congeners were inversely associated with concentrations of TC and LDL-C in studies from areas producing BFRs in China [32]. However, some studies have reported positive correlation results, such as PBDE99, PBDE47, and PBDE28 being positively associated with TC levels in pregnant US women [33, 34]. In these studies, the study population was limited to specific populations, such as residents of areas where BFRs are produced or pregnant women. The sample sizes were also limited, which may have resulted in insufficient generalizability of the results.

Environmental pollutants usually exist as mixtures, and their combined toxicity may have synergistic or antagonistic effects [35, 36]. Interactions between environmental pollutants can have significant effects on health, so assessing the mixing effects of pollutants is essential for a comprehensive assessment of their health risks. To comprehensively examine the association between single and mixed BFR exposure and lipid levels, a cross-sectional study was conducted using data from the large, representative population of the National Health and Nutrition Examination Survey (NHANES) and systematically adjusted for multiple confounding variables, including demographic characteristics and health status, to improve the robustness of the results.

Methods

Study population

This study extracted data spanning six NHANES cycles. A total of 34,180 adult participants (≥ 20 years) were initially considered from the 2005 to 2016 NHANES cohorts. The following are the exclusion criteria on participants: (1) missing data on serum BFRs (N = 24,296); (2) missing data on blood lipids (N = 5,230); (3) participants taking statins within 30 days prior to the examination (including rosuvastatin, fluvastatin, atorvastatin, pitavastatin, pravastatin, lovastatin, and simvastatin, N = 826), pregnant women (N = 75), participants with cancer (N = 282), participants missing data on body mass index (BMI) (N = 39), ratio of family income to poverty (PIR) (N = 271), education level (N = 2), smoking status (N = 3), diabetes (N = 2). In conclusion, the study included 3,154 participants who met the eligibility criteria (Fig. 1).

Fig. 1
figure 1

Flowchart of inclusion criteria for NHANES study participants (2005–2016)

Serum brominated flame retardants

The serum BFRs were quantified using isotope dilution gas chromatography/high-resolution mass spectrometry (GC/HRMS) following the protocol outlined by Sjodin et al. [37]. Serum samples underwent automated liquid-liquid extraction and chromatographic purification to remove contaminants. Isotope-labeled internal standards were added to ensure accurate quantification. PBDEs were separated by gas chromatography and quantified using high-resolution mass spectrometry in selective ion monitoring mode. Please refer to the Supplementary file for detailed information. (Description of laboratory methodology) Although the analytical method was consistent across NHANES cycles, slight variations in the limit of detection (LOD) occurred due to updates in instrumentation and calibration. Serum BFR concentrations below the lower limit of detection (LLOD) were imputed as the LLOD value divided by the square root of 2, following NHANES guidelines. To address these differences, the highest LOD for each BFR across all cycles was used as the reference threshold. Details of the proportion of samples below the LOD across NHANES cycles are provided in supplement file. (Table S1).

In order to guarantee the precision of the analytical outcomes, only BFRs with a total detection rate higher than 70% were included in this study. PBDE183, PBDE66, and PBDE17 were excluded from further consideration due to their relatively low detection rates of 46.0%, 14.0%, and 5.5% respectively. Ultimately, eight PBDEs including Decabromodiphenyl ether (PBDE209, LOD: 71.4%), 2, 2´, 4, 4´, 5, 6´-Hexabromodiphenyl ether (PBDE154, LOD: 76.1%), 2, 2´, 4, 4´, 5, 5´-Hexabromodiphenyl ether (PBDE153, LOD: 100%), 2, 2´, 4, 4´, 6-Pentabromodiphenyl ether (PBDE100, LOD: 100%), 2, 2´, 4, 4´, 5-Pentabromodiphenyl ether (PBDE99, LOD: 100%), 2, 2´, 3, 4, 4´-Tentabromodiphenyl ether (PBDE85, LOD: 74.9%), 2, 2´, 4, 4´-Tetrabromodiphenyl ether (PBDE47, LOD: 100%), and 2, 4, 4´-Tribromodiphenyl ether (PBDE28, LOD: 97.3%) and 2, 2´, 4, 4´, 5, 5´-Hexabromobiphenyl (PBB-153, LOD: 98.4%) were included in the further analyses.

Blood lipid profiles

HDL-C (mg/dL) was measured using a method in which non-HDL cholesterol forms water-soluble complexes with magnesium/dextran sulfate, followed by the enzymatic conversion of HDL cholesterol esters into HDL cholesterol, which was then quantified photometrically. TC (mg/dL) was assessed using an enzymatic assay that converts esterified cholesterol into free cholesterol, producing hydrogen peroxide, which reacts with 4-aminophenazone to form a measurable product at 505 nm. The Friedewald equation was used to determine LDL-C (mg/dL): [LDL-C] = [TC] − [HDL-C] − [TG/5]​ [38]. According to guidelines, RC (mg/dL) was derived using the formula: [RC] = [TC] − [HDL-C] − [LDL-C] [39].

Covariates

On the basis of previous studies, a range of covariates were included in the analyses [32, 40]. Demographic covariates included gender (male, female), age (< 50 years, > 50 years), race (Mexican American, other Hispanic, non-Hispanic White, non-Hispanic Black, and other race), education level (under high school, high school, and college graduate), and PIR (< 1, 1–3, and > 3). Other covariates including BMI (< 25, 25–30, and > 30 kg/m2), smoking status (yes, no), diabetes (yes, no, and borderline), and research cycle (2005–2006, 2007–2008, 2009–2010, 2011–2012, 2013–2014, and 2015–2016).

Statistical analyses

Continuous and categorical variables were described by mean ± standard deviation (SD) and number (%), respectively. The chi-square test and t-test were employed for the purpose of evaluating the continuous and categorical variables of participants according to gender. The continuous variables were transformed using natural logarithms in order to fit a normal distribution. MEC weights were applied in the analyses to conform to the complex, multistage sampling design [41].

Spearman correlation analysis was used to assess the relationships between nine serum BFRs. Pairwise correlation coefficients range from −1 to 1, indicating the strength and direction of the associations. The results were visualized in a heatmap, with color intensity reflecting the magnitude of the correlations.

Weighted multivariable regression analyses were used to evaluate the association between serum BFRs and TC, HDL-C, LDL-C, and RC resulting in beta values and a 95% CI. Three models were constructed in the analyses: model 1: non-adjusted; model 2: adjusted for gender, age, race, education level, PIR, smoking status, and research cycle; and model 3: adjusted for all covariates. Restricted Cubic Spline (RCS) regression was utilized to investigate the nonlinear relationship between serum BFRs and blood lipids [42]. The median concentration of BFRs was selected as the reference point, and four knots were positioned at the 5th, 35th, 65th, and 95th percentiles of the BFR distribution. In order to identify potential gender-specific differences, stratified mixture analyses were performed by gender.

The overall association between multiple BFRs and blood lipids was assessed using two models including Quantile G-computation (QGC) analyses and weighted quantile sum (WQS) regression [43, 44]. WQS regression and QGC analyses are advanced statistical techniques employed to evaluate the combined effects of multiple environmental exposures. In the WQS model, the BFRs were divided into quartiles (q = 4), meaning the observed effect corresponds to an increment from one quartile to the next. The dataset was partitioned into a training (60%) and validation (40%) subsets. To estimate the weights of the BFRs, 1,000 bootstrap iterations were performed. In the QGC model, exposures were divided into deciles (q = 10), with the estimates representing the effect of a one-decile increase in the exposure mixture. A generalized linear model was employed to estimate the joint effect, and bootstrap resampling was performed to calculate confidence intervals for the joint effect estimates.

Statistical analyses were performed using the R software 4.3.2 and EmpowerStats. P-value less than 0.05 was defined as statistical significance.

Results

Baseline characteristics

This study encompassed 3,154 participants (48.32% male and 51.68% female) with an average age of 44.91 ± 16.29 years. The baseline characteristics of participants by gender are summarized in Table 1. Women had a higher proportion of individuals aged ≥ 50 years. Regarding BMI, women were more likely to have a BMI ≥ 30, whereas men had a higher proportion with BMI in the 25–30 range. Men also exhibited higher levels of serum BFRs, including PBB153, PBDE153, PBDE100, PBDE99, PBDE85, PBDE47 and PBDE28. Women had higher levels of TC and HDL, while men exhibited higher levels of RC (Table 1). The overall correlation between serum BFRs was positive (Fig. 2).

Table 1 Baseline characteristics of participants from NHANES 2005–2016
Fig. 2
figure 2

Spearman correlation coefficient plot showing the relationships between serum BFRs. Correlation coefficients (range: −1 to 1) are visualized by color intensity, indicating the strength of associations

Association between individual serum BFR and blood lipids

The results showed that serum PBB153 (β = 0.034), PBDE209 (β = 0.021), PBDE153 (β = 0.024), and PBDE28 (β = 0.015) were positively associated with TC after adjusting for covariates. PBB153 (β = 0.039), PBDE209 (β = 0.029), and PBDE153 (β = 0.030) were positively associated with LDL-C. PBB153 (β = 0.064), PBDE209 (β = 0.046), PBDE153 (β = 0.042), and PBDE28 (β = 0.036) were positively associated with RC. However, the association between single BFR and HDL-C was not statistically significant (Table 2). Please refer to supplement file for detailed results on multivariable regression model. (Table S2).

Table 2 Associations of individual BFRs with blood lipid levels in the general population

Restricted cubic spline (RCS) curves were used to further explore the non-linear relationships between serum BFRs and blood lipid profiles. The results obtained revealed notable threshold effects for specific serum BFRs. For PBB153, threshold effects were observed for TC at approximately 3.5, LDL-C at 3.0, and RC at 3.2, beyond which the effects plateaued. Similar threshold patterns were noted for PBDE209, with TC, LDL-C, and RC showing diminishing effects at higher concentrations. Detailed results are shown in Supplemental File. (Figure S1-4).

Association between mixed serum BFRs with blood lipids

The results of WQS regression showed that mixed serum BFRs were positively associated with TC (β = 0.035), LDL-C (β = 0.031), RC (β = 0.057), and HDL-C (β = 0.020) after adjusting for all covariates. (Table 3) In addition, WQS regression suggested that PBB153 contributed most to the overall mixture effect on TC (54.4%), HDL-C (69.2%), and LDL-C (46.7%). PBDE28 contributed most on RC (38.3%) (Fig. 3A-D). Scatter plot showed that serum BFRs were positively correlated with TC, LDL-C, RC and HDL-C overall. (Figure S5).

Table 3 Association between mixed serum BFRs and blood lipids In the general population by WQS regression and QCG analyses
Fig. 3
figure 3

Effects of individual BFRs to blood lipid levels in the general population by WQS regression. The following covariates were adjusted: gender, age, race, education level, PIR, BMI, smoking status, diabetes, and research cycle. Proportion of effect of single BFR on TC; Proportion of effect of single BFR on HDL-C; Proportion of effect of single BFR on LDL-C; Proportion of effect of single BFR on RC. QGC: quantile g-computation; WQS: weighted quantile sum

QGC analyses suggested that mixed serum BFRs were positively associated with TC (β = 0.036), LDL-C (β = 0.036), and RC (β = 0.065). However, the association between mixed BFRs and HDL-C was also not statistically significant. (Table 3) The results of the single exposure weighting in the QGC analyses were similar with the WQS regression. The results suggested that PBB153 contributed most on TC and LDL-C. PBDE28 contributed most on RC. (Figures 4A-D and 5 (A–D) illustrated the associations between joint exposure to BFRs and blood lipid levels based on QGC analyses. Taking the first decile of the BFR mixture as a reference, TC, LDL-C and RC showed a significant upward trend with increasing mixture concentration. (P < 0.05). The association with HDL-C was positive but not statistically significant (Fig. 5A-D).

Fig. 4
figure 4

Effects of individual BFRs on blood lipids in the general population by QGC analysis. The following covariates were adjusted: gender, age, race, education level, PIR, BMI, smoking status, diabetes, and research cycle. A: Effect of single BFR on TC; B: Effect of single BFR on HDL-C; C: Effect of single BFR on LDL-C; D: Effect of single BFR on RC. QGC: quantile g-computation; WQS: weighted quantile sum

Fig. 5
figure 5

Overall effect of mixed BFRs on blood lipids in the general population by QGC analyses. The following covariates were adjusted: gender, age, race, education level, PIR, BMI, smoking status, diabetes, and research cycle. A: Overall effect of BFRs on TC; B: Overall effect of BFRs on HDL-C; C: Overall effect of BFRs on LDL-C; D: Overall effect of BFRs on RC. QGC: quantile g-computation; WQS: weighted quantile sum

Stratified analyses by gender

The associations between individual BFRs and blood lipids were examined using multivariable regression models. The results indicated significant positive associations between certain BFRs (e.g., PBB153) and TC, LDL-C, and RC in both males and females after adjusted covariates. However, the magnitude of these associations varied by gender and BFR type. Detailed results are provided in supplemental file. (Table S3 and S4) The stratified mixture analyses revealed significant associations between mixed serum BFRs and blood lipid outcomes in both males and females. For males, WQS and QGC models consistently showed positive associations with TC and RC, while associations with HDL-C were not statistically significant. Mixed serum BFRs were positively associated with TC, LDL-C, and HDL-C in females (Table 4). There were some differences between the results of the stratified mixed analysis and the results of the overall participants.

Table 4 Gender-stratified associations between mixed serum BFRs and blood lipids using WQS and QGC models

Discussion

This study explored the association between exposure to BFRs and blood lipid profiles in American adults using NHANES data. The results showed that individual PBB153, PBDE209, PBDE153, and PBDE28 were positively associated with TC and RC after adjusted all covariates. In addition, PBB153, PBDE209, and PBDE153 were positively associated with LDL-C, while no significant association was found between any individual BFR and HDL-C. Further analysis using WQS regression and QGC analyses confirmed that BFR mixtures were positively associated with TC, LDL-C, and RC. The main results of the stratified analyses were consistent with those in the overall population. The findings suggest that BFR exposure may lead to increased levels of TC, LDL-C, and RC, indicating a potential link to dyslipidemia and increased risk of cardiovascular diseases.

The findings of some previous studies contrast with our results. For instance, a study of 172 participants from Laizhou Bay, China, a major flame retardant production area, found that serum levels of various PBDE congeners—including PBDE209, PBDE154, PBDE100, PBDE99, and total PBDEs—were marginally to significantly negatively correlated with TC and LDL-C [32]. Similarly, another study of 150 female participants from Shantou, China, reported that PBDE-190 levels in adipose tissue were inversely associated with TC, while total PBDEs were positively associated with TC [45]. Moreover, research has indicated that exposure to PBDE99 during the prenatal period is associated with reduced levels of TG in childhood [31]. In contrast, several studies have yielded results consistent with our findings. A study reported a positive association between PBDE99, PBDE47, and PBDE28 with TC among 388 pregnant women based on the HOME project [33, 34]. However, the current research on the relationship between BFR exposure and blood lipids remains limited, and the findings are inconclusive. The discrepancies across studies may be attributable to variations in the sources of BFRs, ethnic differences among participants, gender differences, and the control of confounding factors in the analyses. Some studies have investigated the potential causal relationship between BFR exposure and blood lipids using animal models. For example, one study reported elevated plasma TC levels in female offspring of C57BL/6N mice exposed perinatally to PBDE-71, although similar changes were not observed in male offspring [46].

One potential mechanism linking BFR exposure to altered lipid profiles is through the disruption of thyroid hormone pathways, which regulate hepatic cholesterol metabolism and lipoprotein receptor gene expression [47,48,49]. Additionally, a study demonstrated that PBDE-47 upregulated microRNA-34a-5p, leading to NAD + deficiency, impaired mitophagy, and induced mitochondrial dysfunction and oxidative damage in the liver [50]. BFR exposure may also induce oxidative stress and cause inflammation [51, 52]. All of these may disrupt lipid metabolism. Future studies are needed to clarify the mechanism of blood lipid changes induced by exposure to BFRs.

This study proposed a potential association between BFR exposure and elevated RC levels. Our results showed that PBB153 and three PBDEs included in the analyses were positively associated with RC, with a stronger effect than that observed for TC and LDL-C after adjusting for all covariates. A large-scale study conducted in the Danish population found that individuals with RC levels above 1 mg/dL, compared to those with levels below 0.5 mg/dL, had hazard ratios (HRs) for mortality of 2.2 (95% CI: 1.3–3.5) from cardiovascular disease and 2.1 (95% CI: 1.4–3.3) from other causes [53, 54]. Additionally, many studies have reported a potential causal association between RC and various cardiovascular diseases, including ischemic heart disease, myocardial infarction, and aortic stenosis through Mendelian randomization (MR) [55, 56]. Identifying risk factors for elevated RC is therefore of great significance, as controlling RC levels could play a crucial role in promoting human health.

The overall effect of BFR exposure on lipid profiles was assessed using two distinct analytical methods: WQS regression and QGC analysis. Both approaches showed that mixed serum BFRs were positively correlated with RC, LDL-C, and TC. The findings indicated that PBB153 was the most influential BFR contributing to the mixtures for TC and LDL-C. Additionally, PBDE28 was identified as the most influential BFR contributing to the mixtures for RC, while PBB153 exhibited a comparatively minimal contribution. Following a stratified analysis based on sex, the results demonstrated that, for males, the WQS and QGC models consistently exhibited positive associations between serum BFRs and TC and RC, while associations with HDL-C were not statistically significant. However, mixed serum BFRs were positively associated with TC, LDL-C, and HDL-C in females. The extant studies on this subject is limited, and further research is required to explore the mechanism causing this gender difference.

The WQS regression analysis demonstrated a positive correlation between mixed BFRs and HDL-C, yet this relationship did not attain statistical significance in QGC analysis. In QGC analyses, interaction and non-linear effects between exposures are captured, potentially leading to more significant impacts of certain exposures. In contrast, WQS regression combines exposures through linear weighting, which may overlook these complex interactions. Additionally, WQS assumes that all exposures affect the outcome in the same direction, potentially masking opposing effects, while QGC allows each exposure to have an independent direction [43, 44]. Moreover, the inability to account for NHANES sampling weights in WQS regression and QGC analyses could lead to differences in results compared to those from the multivariable regression model. Therefore, it is important to integrate findings from multiple models when evaluating the effect of multiple factors on health outcomes.

This study has several strengths. It was the first large-scale cross-sectional study to systematically explore the association between exposure to BFRs and TC, LDL-C, RC, and HDL-C. This study utilized data from the NHANES, which is representative of the general population. Furthermore, our analyses employed novel methods, including WQS and QGC models, to investigate the mixed effects of BFRs. The primary outcomes of the stratified analyses were consistent with the overall population, thereby enhancing the robustness of the study. The findings offered compelling evidence that BFRs are detrimental to cardiovascular health.

This study also has some limitations. Based on the cross-sectional study design, our results can only demonstrate an association between exposure to BFRs and blood lipid profiles, but not a causal relationship. The study population was drawn exclusively from the United States, and therefore the results cannot be generalized to other countries. Many participants were excluded because of missing information, and participants taking statins, those with a history of cancer, and pregnant women were also excluded, which may prevent the results from adequately reflecting the characteristics of the general population. Despite adjustments made for primary confounders, the possibility of residual confounding remains unconfirmed. Lifestyle includes diet, work and other factors, a potential source of both BFR exposure and blood lipid variation, were not included due to limitations in the data. Similarly, other co-occurring endocrine disruptors, such as heavy metals and phthalates, were not considered, potentially underestimating the synergistic or antagonistic effects of environmental pollutants. In addition, Due to the inherent limitation of the detection rate, BFRs with a detection rate below 70% were excluded from our analyses to minimize the impact of measurement error. However, for the included BFRs, some samples 897still exhibited concentrations the LOD, which were imputed as LOD/√2 according to NHANES guidelines. This imputation approach, while widely used, may introduce some degree of error and potentially bias the results. Furthermore, the analyses were limited to nine major PBDE congeners and PBB153 available in the NHANES. While these congeners account for a significant portion of human exposure, analyzing additional PBDE congeners could offer a more comprehensive understanding of their associations with blood lipids. Future studies with a broader range of BFRs and larger datasets are needed to validate our findings and further elucidate their health impacts.

Conclusion

This study demonstrates that BFR exposure is associated with increased levels of TC, LDL-C, and RC, indicating an elevated risk of dyslipidemia and cardiovascular diseases. These results highlight the potential public health implications of BFR exposure and underscore the need for further research into the underlying mechanisms and preventive strategies.

Data availability

The data are available from NHANES (https://www.cdc.gov/nchs/nhanes/).

References

  1. GBD 2017 Risk Factor Collaborators. Global, regional, and national comparative risk assessment of 84 behavioural, environmental and occupational, and metabolic risks or clusters of risks for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018;392(10159):1923–94.

  2. Grundy SM, Stone NJ, Bailey AL, Beam C, Birtcher KK, Blumenthal RS, Braun LT, de Ferranti S, Faiella-Tommasino J, Forman DE, et al. 2018 AHA/ACC/AACVPR/AAPA/ABC/ACPM/ADA/AGS/APhA/ASPC/NLA/PCNA guideline on the management of blood cholesterol: executive summary: a report of the American College of Cardiology/American Heart Association Task Force on clinical practice guidelines. Circulation. 2019;139(25):e1046-81.

    PubMed  Google Scholar 

  3. von Eckardstein A, Nordestgaard BG, Remaley AT, Catapano AL. High-density lipoprotein revisited: biological functions and clinical relevance. Eur Heart J. 2023;44(16):1394–407.

    Article  Google Scholar 

  4. Crea F. High-density lipoproteins, lipoprotein(a), and remnant cholesterol: new opportunities for reducing residual cardiovascular risk. Eur Heart J. 2023;44(16):1379–82.

    Article  PubMed  Google Scholar 

  5. Miller GJ, Miller NE. Plasma-high-density-lipoprotein concentration and development of ischaemic heart-disease. Lancet. 1975;1(7897):16–9.

    Article  CAS  PubMed  Google Scholar 

  6. Di Angelantonio E, Sarwar N, Perry P, Kaptoge S, Ray KK, Thompson A, Wood AM, Lewington S, Sattar N, Packard CJ, et al. Major lipids, apolipoproteins, and risk of vascular disease. JAMA. 2009;302(18):1993–2000.

    Article  PubMed  Google Scholar 

  7. Castelli WP, Garrison RJ, Wilson PW, Abbott RD, Kalousdian S, Kannel WB. Incidence of coronary heart disease and lipoprotein cholesterol levels. The Framingham Study. JAMA. 1986;256(20):2835–8.

    Article  CAS  PubMed  Google Scholar 

  8. Trinder M, Walley KR, Boyd JH, Brunham LR. Causal inference for genetically determined levels of high-density lipoprotein cholesterol and risk of infectious disease. Arterioscler Thromb Vasc Biol. 2020;40(1):267–78.

    Article  CAS  PubMed  Google Scholar 

  9. Pedersen KM, Çolak Y, Bojesen SE, Nordestgaard BG. Low high-density lipoprotein and increased risk of several cancers: 2 population-based cohort studies including 116,728 individuals. J Hematol Oncol. 2020;13(1):129.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Madsen CM, Varbo A, Nordestgaard BG. Low HDL cholesterol and high risk of autoimmune disease: two population-based cohort studies including 117341 individuals. Clin Chem. 2019;65(5):644–52.

    Article  CAS  PubMed  Google Scholar 

  11. Mudd JO, Borlaug BA, Johnston PV, Kral BG, Rouf R, Blumenthal RS, Kwiterovich PO Jr. Beyond low-density lipoprotein cholesterol: defining the role of low-density lipoprotein heterogeneity in coronary artery disease. J Am Coll Cardiol. 2007;50(18):1735–41.

    Article  CAS  PubMed  Google Scholar 

  12. Quispe R, Martin SS, Michos ED, Lamba I, Blumenthal RS, Saeed A, Lima J, Puri R, Nomura S, Tsai M, et al. Remnant cholesterol predicts cardiovascular disease beyond LDL and ApoB: a primary prevention study. Eur Heart J. 2021;42(42):4324–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Li W, Huang Z, Fang W, Wang X, Cai Z, Chen G, Wu W, Chen Z, Wu S, Chen Y. Remnant cholesterol variability and incident ischemic stroke in the general population. Stroke. 2022;53(6):1934–41.

    Article  CAS  PubMed  Google Scholar 

  14. Huang H, Guo Y, Liu Z, Zeng Y, Chen Y, Xu C. Remnant cholesterol predicts long-term mortality of patients with metabolic dysfunction-associated fatty liver disease. J Clin Endocrinol Metab. 2022;107(8):e3295–303.

    Article  PubMed  Google Scholar 

  15. Nordestgaard BG, Varbo A. Triglycerides and cardiovascular disease. Lancet. 2014;384(9943):626–35.

    Article  CAS  PubMed  Google Scholar 

  16. Wang Y, Li H, Zhang Z, Wu F, Liu J, Zhu Z, Xiang H. The association between vitamin E intake and remnant cholesterol, total cholesterol, high-density lipoprotein cholesterol, and low-density lipoprotein cholesterol in US adults: a cross-sectional study. Lipids Health Dis. 2024;23(1):325.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Tian Y, Wu Y, Qi M, Song L, Chen B, Wang C, Lu J, Yang Y, Zhang X, Cui J, et al. Associations of remnant cholesterol with cardiovascular and cancer mortality in a nationwide cohort. Sci Bull (Beijing). 2024;69(4):526–34.

    Article  CAS  PubMed  Google Scholar 

  18. Birnbaum LS, Staskal DF. Brominated flame retardants: cause for concern? Environ Health Perspect. 2004;112(1):9–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Chen L, Yu K, Huang C, Yu L, Zhu B, Lam PK, Lam JC, Zhou B. Prenatal transfer of polybrominated diphenyl ethers (PBDEs) results in developmental neurotoxicity in zebrafish larvae. Environ Sci Technol. 2012;46(17):9727–34.

    Article  CAS  PubMed  Google Scholar 

  20. Bradman A, Castorina R, Sjödin A, Fenster L, Jones RS, Harley KG, Chevrier J, Holland NT, Eskenazi B. Factors associated with serum polybrominated diphenyl ether (PBDE) levels among school-age children in the CHAMACOS cohort. Environ Sci Technol. 2012;46(13):7373–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Dvoršćak M, Jakovljević I, Jagić K, Tariba Lovaković B, Klinčić D. Polybrominated diphenyl ethers and polycyclic aromatic hydrocarbons in dust from different indoor environments in Zagreb, Croatia: levels and human exposure assessment. Indoor Air. 2022;32(11):e13145.

    Article  PubMed  Google Scholar 

  22. Hohenblum P, Steinbichl P, Raffesberg W, Weiss S, Moche W, Vallant B, Scharf S, Haluza D, Moshammer H, Kundi M, et al. Pollution gets personal! A first population-based human biomonitoring study in Austria. Int J Hyg Environ Health. 2012;215(2):176–9.

    Article  CAS  PubMed  Google Scholar 

  23. Rani M, Keshu, Meenu, Sillanpää M, Shanker U. An updated review on environmental occurrence, scientific assessment and removal of brominated flame retardants by engineered nanomaterials. J Environ Manage. 2022;321:115998.

  24. Zuiderveen EAR, Slootweg JC, de Boer J. Novel brominated flame retardants - a review of their occurrence in indoor air, dust, consumer goods and food. Chemosphere. 2020;255:126816.

    Article  CAS  PubMed  Google Scholar 

  25. Tao W, Nian W, Li L. Analysis of brominated flame retardants exposure-associated chronic kidney disease risk in the US population from the NHANES. Ecotoxicol Environ Saf. 2024;286:117159.

    Article  CAS  PubMed  Google Scholar 

  26. Pan Y, Chen Q, Yu Y, Yang H, Liu Z, Xie B, Huang Y, He B, Yan F, Chen F, et al. Association between brominated flame retardants (BFRs) and periodontitis: Results from a large population-based study. BMC Oral Health. 2024;24(1):1025.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Bai T, Li X, Zhang H, Yang W, Lv C, Du X, Xu S, Zhao A, Xi Y. The association between brominated flame retardants exposure with bone mineral density in US adults: a cross-sectional study of the national health and nutrition examination survey (NHANES) 2005–2014. Environ Res. 2024;251(Pt 1):118580.

    Article  CAS  PubMed  Google Scholar 

  28. Che Z, Jia H, Chen R, Pan K, Fan Z, Su C, Wu Z, Zhang T. Associations between exposure to brominated flame retardants and metabolic syndrome and its components in U.S. adults. Sci Total Environ. 2023;858(Pt 2):159935.

    Article  CAS  PubMed  Google Scholar 

  29. Cheng D, Chen Z, Zhou J, Cao Y, Xie X, Wu Y, Li X, Wang X, Yu J, Yang B. Association between brominated flame retardants (PBDEs and PBB153) exposure and hypertension in U.S. adults: results from NHANES 2005–2016. Environ Health. 2024;23(1):64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Tung EWY, Kawata A, Rigden M, Bowers WJ, Caldwell D, Holloway AC, Robaire B, Hales BF, Wade MG. Gestational and lactational exposure to an environmentally-relevant mixture of brominated flame retardants: effects on neurodevelopment and metabolism. Birth Defects Res. 2017;109(7):497–512.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Boutot ME, Whitcomb BW, Abdelouahab N, Baccarelli AA, Boivin A, Caku A, Gillet V, Martinez G, Pasquier JC, Zhu J et al. In utero exposure to persistent organic pollutants and childhood lipid levels. Metabolites. 2021;11(10):657. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/metabo11100657.

  32. Zhao X, Yang X, Du Y, Li R, Zhou T, Wang Y, Chen T, Wang D, Shi Z. Polybrominated diphenyl ethers in serum from residents living in a brominated flame retardant production area: Occurrence, influencing factors, and relationships with thyroid and liver function. Environ Pollut. 2021;270:116046.

    Article  CAS  PubMed  Google Scholar 

  33. Braun JM, Buckley JP, Cecil KM, Chen A, Kalkwarf HJ, Lanphear BP, Xu Y, Woeste A, Yolton K. Adolescent follow-up in the Health Outcomes and Measures of the Environment (HOME) Study: cohort profile. BMJ Open. 2020;10(5):e034838.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Vuong AM, Braun JM, Sjödin A, Calafat AM, Yolton K, Lanphear BP, Chen A. Exposure to endocrine disrupting chemicals (EDCs) and cardiometabolic indices during pregnancy: The HOME Study. Environ Int. 2021;156:106747.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Martin O, Scholze M, Ermler S, McPhie J, Bopp SK, Kienzler A, Parissis N, Kortenkamp A. Ten years of research on synergisms and antagonisms in chemical mixtures: a systematic review and quantitative reappraisal of mixture studies. Environ Int. 2021;146:106206.

    Article  CAS  PubMed  Google Scholar 

  36. Cedergreen N. Quantifying synergy: a systematic review of mixture toxicity studies within environmental toxicology. PLoS One. 2014;9(5):e96580.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Sjödin A, Jones RS, Lapeza CR, Focant JF, McGahee EE 3rd, Patterson DG Jr. Semiautomated high-throughput extraction and cleanup method for the measurement of polybrominated diphenyl ethers, polybrominated biphenyls, and polychlorinated biphenyls in human serum. Anal Chem. 2004;76(7):1921–7.

    Article  PubMed  Google Scholar 

  38. Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem. 1972;18(6):499–502.

    Article  CAS  PubMed  Google Scholar 

  39. Mach F, Baigent C, Catapano AL, Koskinas KC, Casula M, Badimon L, Chapman MJ, De Backer GG, Delgado V, Ference BA, et al. 2019 ESC/EAS Guidelines for the management of dyslipidaemias: lipid modification to reduce cardiovascular risk. Eur Heart J. 2020;41(1):111–88.

    Article  PubMed  Google Scholar 

  40. Huang SS, Ding Y, Mao HY, Jin SF, Zheng MX, Yang ZX. Response of platelet count in the US population following exposure to BDEs and BB: A cross-sectional analysis of database of NHANES (2007–2016). Ecotoxicol Environ Saf. 2024;288:117384.

    Article  CAS  PubMed  Google Scholar 

  41. Akinbami LJ, Chen TC, Davy O, Ogden CL, Fink S, Clark J, Riddles MK, Mohadjer LK. National Health and Nutrition Examination Survey, 2017-March 2020 prepandemic file: sample design, estimation, and analytic guidelines. Vital Health Stat 1. 2022;190:1–36.

    Google Scholar 

  42. Desquilbet L, Mariotti F. Dose-response analyses using restricted cubic spline functions in public health research. Stat Med. 2010;29(9):1037–57.

    Article  PubMed  Google Scholar 

  43. Carrico C, Gennings C, Wheeler DC, Factor-Litvak P. Characterization of weighted quantile sum regression for highly correlated data in a risk analysis setting. J Agric Biol Environ Stat. 2015;20(1):100–20.

    Article  PubMed  Google Scholar 

  44. 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(4):47004.

    Article  PubMed  Google Scholar 

  45. Zhang Q, Peng J, Huang A, Zheng S, Shi X, Li B, Huang W, Tan W, Wang X, Wu K. Associations between polybrominated diphenyl ethers (PBDEs) levels in adipose tissues and blood lipids in women of Shantou, China. Environ Res. 2022;214(Pt 3):114096.

    Article  CAS  PubMed  Google Scholar 

  46. Kozlova EV, Denys ME, Benedum J, Valdez MC, Enriquez D, Bishay AE, Chinthirla BD, Truong E, Krum JM, DiPatrizio NV, et al. Developmental exposure to indoor flame retardants and hypothalamic molecular signatures: Sex-dependent reprogramming of lipid homeostasis. Front Endocrinol (Lausanne). 2022;13:997304.

    Article  PubMed  Google Scholar 

  47. Khani L, Martin L, Pułaski Ł. Cellular and physiological mechanisms of halogenated and organophosphorus flame retardant toxicity. Sci Total Environ. 2023;897:165272.

    Article  CAS  PubMed  Google Scholar 

  48. Yeshoua B, Romero Castillo H, Monaghan M, van Gerwen M. A review of the association between exposure to flame retardants and thyroid function. Biomedicines. 2024;12(6):1365. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/biomedicines12061365.

  49. Duntas LH, Brenta G. A renewed focus on the association between thyroid hormones and lipid metabolism. Front Endocrinol (Lausanne). 2018;9:511.

    Article  PubMed  Google Scholar 

  50. Chen F, Feng L, Zheng YL, Lu J, Fan SH, Shan Q, Zheng GH, Wang YJ, Wu DM, Li MQ, et al. 2, 2’, 4, 4’-tetrabromodiphenyl ether (BDE-47) induces mitochondrial dysfunction and related liver injury via eliciting miR-34a-5p-mediated mitophagy impairment. Environ Pollut. 2020;258:113693.

    Article  CAS  PubMed  Google Scholar 

  51. Wang Z, Zhou Y, Xiao X, Liu A, Wang S, Preston RJS, Zaytseva YY, He G, Xiao W, Hennig B, et al. Inflammation and cardiometabolic diseases induced by persistent organic pollutants and nutritional interventions: effects of multi-organ interactions. Environ Pollut. 2023;339:122756.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Chen JC, Baumert BO, Li Y, Li Y, Pan S, Robinson S, Rubbo B, Costello E, He J, Hampson H, et al. Associations of per- and polyfluoroalkyl substances, polychlorinated biphenyls, organochlorine pesticides, and polybrominated diphenyl ethers with oxidative stress markers: a systematic review and meta-analysis. Environ Res. 2023;239(Pt 1):117308.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Wadström BN, Pedersen KM, Wulff AB, Nordestgaard BG. Elevated remnant cholesterol, plasma triglycerides, and cardiovascular and non-cardiovascular mortality. Eur Heart J. 2023;44(16):1432–45.

    Article  PubMed  Google Scholar 

  54. Varbo A, Nordestgaard BG. Nonfasting triglycerides, low-density lipoprotein cholesterol, and heart failure risk: two cohort studies of 113 554 individuals. Arterioscler Thromb Vasc Biol. 2018;38(2):464–72.

    Article  CAS  PubMed  Google Scholar 

  55. Guan B, Wang A, Xu H. Causal associations of remnant cholesterol with cardiometabolic diseases and risk factors: a mendelian randomization analysis. Cardiovasc Diabetol. 2023;22(1):207.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Kaltoft M, Langsted A, Nordestgaard BG. Triglycerides and remnant cholesterol associated with risk of aortic valve stenosis: Mendelian randomization in the Copenhagen General Population Study. Eur Heart J. 2020;41(24):2288–99.

    Article  CAS  PubMed  Google Scholar 

Download references

Funding

This work was supported by the Young Taishan Scholars Program (tsqn201909190); National Natural Science Foundation of China (82172478); Shandong Higher Education Young Science and Technology Support Program (2021KJ048); Youth Talent Assistance Program of Medical College of Qingdao University; “Medicine +” Subject Cluster Joint Exploration Project of Medical School, Qingdao University; Qingdao Key Technology Research Project (2025) (25-1-1-gjgg-61-nsh).

Author information

Authors and Affiliations

Authors

Contributions

Y.W. and Z.H.Z. contributed to the of the work and writing of the Manuscript. N.S., X.Q., H.L.and F.W. contributed to the acquisition and analysisof data. Z.Z.Z. and J.L. contributed to Production of graphs and tables. H.X. contributed to the conception of the work and Revision of the manuscript. All authors reviewed the manuscript.

Corresponding author

Correspondence to Hongfei Xiang.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, Y., Zhang, Z., Shen, N. et al. Association between exposure to brominated flame retardants (BFRs) and blood lipid profiles in American adults: a cross-sectional study. Lipids Health Dis 24, 120 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12944-025-02527-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12944-025-02527-4