Mixture Effects of Metals, PCBs, Dioxins, and Furans on Liver Function

Quantifying the mixture effects on humans exposed remains challenging because mixture components are correlated and may act bidirectionally by exhibiting nonlinear dose-response relationships, which may contribute to subclinical organ dysfunction. The liver is a vital organ in the body with broad functions, making it vulnerable to injury as it is the first organ exposed to circulating toxicants, which can precipitate hepatic damage. Our study’s objective was to evaluate the combined and component-specific associations of a multi-chemical exposure mixture of heavy metals, polychlorinated biphenyls (PCBs), polychlorinated dibenzo-p-dioxins (dioxins), and polychlorinated dibenzofurans (furans), with liver biomarkers, and to compare concentration-based results with the toxic equivalent (TEQ) potency of the weighted results for dioxin-like compounds. In an unweighted analytic sample of U.S. adults from NHANES 2003–2004 with 947 complete cases, we examined heavy metals (cadmium, lead, and mercury), PCBs (12 congeners), dioxins (7 congeners), and furans (10 congeners) in relation to eight liver biomarkers (albumin, ALP, ALT, AST, GGT, LDH, total bilirubin, and total protein). We applied multi-exposure linear regression, weighted quantile sum (WQS) regression, quantile g-computation (qgcomp), and Bayesian kernel machine regression (BKMR), with parallel TEQ-based models using WHO 2005 TEFs for dioxin-like PCBs, dioxins, and furans. Across mixture methods, the mixture structure was chemically sparse, with a limited set of recurring contributors. Total bilirubin showed the most consistent positive mixture association across qgcomp and BKMR and persisted under TEQ weighting, with prominent PCB- and dioxin-like contributions (notably PCB81/PCB TEQs and dioxin-related components). Albumin demonstrated inverse mixture patterns in BKMR and TEQ-BKMR, with dioxin-like components (notably Dioxin3 and Dioxin3_TEQ) repeatedly emerging as key drivers. For ALT, ALP, AST, GGT, LDH, and total protein, overall mixture effects were frequently attenuated or null in qgcomp despite structured component weights, indicating bidirectional sub-mixtures and internal counterbalancing. BKMR PIPs similarly concentrated on a small number of dominant predictors (e.g., lead for ALP, mercury for ALT, PCB28 for AST, and cadmium and PCB189 for LDH), while interaction summaries provided limited evidence of stable non-additivity. Using multiple complementary mixture methods, we identified outcome-specific mixture patterns suggesting hepatobiliary vulnerability. TEQ concordance supports toxicological relevance of the dioxin-like axis, while metals and non–dioxin-like mechanisms likely contribute additional pathways.