Abstract
Brominated flame retardants (BFRs), commonly found in consumer products, have been identified as potential hazards to liver function. While the individual effects of specific BFRs are somewhat understood, there is limited evidence on how mixtures of these chemicals, especially when influenced by demographic factors, interact to affect liver function. This study utilized data from 10,828 participants aged 12 and above from the National Health and Nutrition Examination Survey (2005-2016) to investigate the associations between BFRs (both individually and in combinations) and biomarkers of liver injury. The study focused on how age, gender, and body mass index (BMI) modify modulate these effects. Multivariate linear regression, restricted cubic spline function, weighted quantile sum (WQS) regression, and quantile g-computation (qgcomp) models were used to analyze the linear, non-linear, and joint associations between BFR levels and liver function parameters. We found positive associations between the mixed BFRs index and AST, ALT, GGT, ALP, and TBIL levels and a negative association with ALB levels. PBDE28, PBDE47, and PBB153 consistently contributed to the top weight in both the WQS and qgcomp models. Most critically, the study demonstrated that the relationship between co-exposure to BFRs and liver function parameters was modified by age, gender, and BMI. Therefore, our study highlights the importance of considering demographic diversity in assessing the risk of BFR-induced liver damage and supports the implementation of tailored preventive and intervention strategies.