Relationship of multiple metals mixture and osteoporosis in older Chinese women: An aging and longevity study.

Osteoporosis has become a major health problem in older women. Previous studies have linked individual metals exposure with osteoporosis, but combined effects remain inconclusive. We aimed to explore the individual and combined association between multiple metals mixture and osteoporosis risk in older Chinese women. A total of 2297 older women (aged ≥60) from the Hongshuihe region of Guangxi, southern China included. We measured 22 blood metal levels through inductively coupled plasma mass spectrometry. And osteoporosis was defined as a T score ≤ -2.5. The least absolute shrinkage and selection operator (LASSO) penalized regression, and Bayesian kernel machine regression (BKMR) models were performed to explore the association between blood metals and osteoporosis risk. Of 2297 older women, there were 829 osteoporosis and 1468 non-osteoporosis participants. The median age was 71 and 68 years old in the osteoporosis and the non-osteoporosis group, respectively. In the single-metal model, rubidium and vanadium were negatively associated with osteoporosis (P for trend = 0.02 and 0.002, respectively), and lead presented the reverse trend (P for trend = 0.01). The LASSO penalized regression model selected nine metals (calcium, cadmium, cobalt, lead, magnesium, rubidium, strontium, vanadium and zinc), which were included in the subsequent analysis. And the multiple-metal model presented a consistent trend with the single-metal model using the selected metals. Furthermore, we performed BKMR to explore the combined effect, and found an overall negative effect between metals mixture and osteoporosis risk when all the metals were fixed at 50th, and rubidium and vanadium were the main contributors. In addition, blood Rb and V were significantly negatively related to OP risk with other metals at different levels (25th, 50th and 75th percentiles). The study suggests metal mixture exposure and osteoporosis risk in older Chinese women, and further studies need to be conducted.

Associations of Heavy Metals with Activities of Daily Living Disability: An Epigenome-Wide View of DNA Methylation and Mediation Analysis.

BACKGROUND: Exposure to heavy metals has been reported to be associated with multiple diseases. However, direct associations and potential mechanisms of heavy metals with physical disability remain unclear. OBJECTIVES: We aimed to quantify associations of heavy metals with physical disability and further explore the potential mechanisms of DNA methylation on the genome scale. METHODS: A cross-sectional study of 4,391 older adults was conducted and activities of daily living (ADL) disability were identified using a 14-item scale questionnaire including basic and instrumental activities to assess the presence of disability (yes or no) rated on a scale of dependence. Odds ratios (ORs) and 95% confidence intervals (CI) were estimated to quantify associations between heavy metals and ADL disability prevalence using multivariate logistic regression and Bayesian kernel machine regression (BKMR) models. Whole blood-derived DNA methylation was measured using the HumanMethylationEPIC BeadChip array. An ADL disability-related epigenome-wide DNA methylation association study (EWAS) was performed among 212 sex-matched ADL disability cases and controls, and mediation analysis was further applied to explore potential mediators of DNA methylation. RESULTS: Each 1-standard deviation (SD) higher difference in log10-transformed manganese, copper, arsenic, and cadmium level was significantly associated with a 14% (95% CI: 1.05, 1.24), 16% (95% CI:1.07, 1.26), 22% (95% CI:1.13, 1.33), and 15% (95% CI:1.06, 1.26) higher odds of ADL disability, which remained significant in the multiple-metal and BKMR models. A total of 85 differential DNA methylation sites were identified to be associated with ADL disability prevalence, among which methylation level at cg220000984 and cg23012519 (annotated to IRGM and PKP3) mediated 31.0% and 31.2% of manganese-associated ADL disability prevalence, cg06723863 (annotated to ESRP2) mediated 32.4% of copper-associated ADL disability prevalence, cg24433124 (nearest to IER3) mediated 15.8% of arsenic-associated ADL disability prevalence, and cg07905190 and cg17485717 (annotated to FREM1 and TCP11L1) mediated 21.5% and 30.5% of cadmium-associated ADL disability prevalence (all p<0.05). DISCUSSION: Our findings suggested that heavy metals contributed to higher prevalence of ADL disability and that locus-specific DNA methylation are partial mediators, providing potential biomarkers for further cellular mechanism studies. https://doi.org/10.1289/EHP10602.

Associations Between Blood Pressure and Accelerated DNA Methylation Aging.

Background Individuals of the same chronological age may exhibit diverse susceptibilities to death. However, few studies have investigated the associations between blood pressure and the accelerated aging. Methods and Results A cross-sectional study was conducted in 288 adults aged ≥50 years. We assessed the DNA methylation-based measures of biological age using CpG sites on the Illumina HumanMethylationEPIC BeadChip. Epigenetic age acceleration metrics were derived by regressing residuals (ΔAge) and ratios (aging rate) of DNA methylation age on chronological age. Dose-response relationships between blood pressure and epigenetic age acceleration were quantified using multiple linear regression and restricted cubic regression models. We found that each 10-mm Hg increase in systolic blood pressure was associated with 0.608 (95% CI, 0.231-0.984) years increase in ΔAge and 0.007 (95% CI, 0.002-0.012) increase in aging rate; meanwhile, for pulse pressure, the increase was 1.12 (95% CI, 0.625-1.61) years for ΔAge and 0.013 (95% CI, 0.007-0.020) for aging rate. Subgroup analysis showed that the significant associations of systolic blood pressure and pulse pressure with epigenetic age acceleration appeared to be limited to women, although interactions between blood pressure and sex were not significant (P values for interaction >0.05). The combination of women and hypertension was associated with a much higher increase in ΔAge (ß [95% CI], 4.05 [1.07-7.02]) and aging rate (ß [95% CI], 0.047 [0.008-0.087]), compared with male participants without hypertension. Conclusions Our findings suggested that high systolic blood pressure and pulse pressure were associated with the epigenetic age acceleration, providing important clues for relationships between blood pressure and epigenetic aging.

Association of both prenatal and early childhood multiple metals exposure with neurodevelopment in infant: A prospective cohort study.

BACKGROUND: Impaired neurodevelopment of children has become a growing public concern; however, the associations between metals exposure and neurocognitive function have remained largely unknown. OBJECTIVES: We systematically evaluated the associations of multiple metals exposure during pregnancy and childhood on the neurodevelopment of children aged 2-3 years. METHODS: We measured 22 metals in the serum and urine among703 mother-child pairs from the Guangxi Birth Cohort Study. The neurocognitive development of children was assessed by the Gesell Development Diagnosis Scale (GDDS; Chinese version). Multiple linear regression models were used to evaluate the relationship between the metals (selected by elastic net regression) and the outcomes. The Bayesian kernel machine regression (BKMR) was used to evaluate the possible joint effect between the multiple metal mixture and the outcomes. RESULTS: Prenatal aluminum (Al) exposure was negatively associated with the fine motor developmental quotient (DQ) (ß = -1.545, 95%CI: 2.231, -0.859), adaption DQ (ß = -1.182, 95%CI: 1.632, -0.732), language DQ (ß = -1.284, 95% CI: 1.758, -0.809), and social DQ (ß = -1.729, 95% CI: 2.406, -1.052) in the multi-metal model. Prenatal cadmium (Cd) exposure was negatively associated with gross motor DQ (ß = -2.524, 95% CI: 4.060, -0.988), while postpartum Cd exposure was negatively associated with language DQ (ß = -1.678, 95% CI: 3.227, -0.129). In stratified analyses, infants of different sexes had different sensitivities to metal exposure, and neurobehavioral development was more significantly affected by metal exposure in the first and second trimester. BKMR analysis revealed a negative joint effect of the Al, Cd, and selenium (Se) on the language DQ score; postpartum Cd exposure played a major role in this relationship. CONCLUSION: Prenatal exposure to Al, Ba, Cd, molybdenum (Mo), lead (Pb), antimony (Sb), and strontium (Sr), and postpartum exposure to cobalt (Co), Cd, stannum (Sn), iron (Fe), nickel (Ni), and Se are associated with neurological development of infants. The first and second trimester might be the most sensitive period when metal exposure affects neurodevelopment.

Causal Effect of Genetically Determined Blood Copper Concentrations on Multiple Diseases: A Mendelian Randomization and Phenome-Wide Association Study.

Exposures to copper have become a health concern. We aim to explore the broad clinical effects of blood copper concentrations. A total of 376,346 Caucasian subjects were enrolled. We performed a Mendelian randomization and phenome-wide association study (MR-PheWAS) to evaluate the causal association between copper and a wide range of outcomes in UK Biobank, and we constructed a protein-protein interaction network. We found association between blood copper concentrations and five diseases in the overall population and nine diseases in male. MR analysis implicated a causal role of blood copper in five diseases (overall population), including prostate cancer (OR = 0.87, 95% CI 0.77-0.98), malignant and unknown neoplasms of the brain and nervous system (OR = 0.58, 95% CI 0.38-0.89), and hypertension (OR = 0.94, 95% CI 0.90-0.98), essential hypertension (OR = 0.94, 95% CI 0.90-0.98) and cancer of brain and nervous system (OR = 0.63, 95% CI 0.41-0.98). For male, except for dysphagia being newly associated with blood copper (OR = 1.39, 95% CI 1.18-1.63), other MR results were consistent with the overall population. In addition, the PPI network showed possible relationship between blood copper and four outcomes, namely brain cancer, prostate cancer, hypertension, and dysphagia. Blood copper may have causal association with prostate cancer, malignant and unknown neoplasms of the brain and nervous system, hypertension, and dysphagia. Considering that copper is modifiable, exploring whether regulation of copper levels can be used to optimize health outcomes might have public health importance. Supplementary Information: The online version contains supplementary material available at 10.1007/s43657-022-00052-3.

Sex-specific associations of exposure to metal mixtures with telomere length change: Results from an 8-year longitudinal study.

Studies on the relationships between exposure to metal mixtures and telomere length (TL) are limited, particularly longitudinal studies. Few studies are available on the potential sex-specific associations between metal exposures and TL change. We examined blood metal concentrations and TL at baseline (August 2012) and follow-up (June 2020) among 316 participants in a ferro-manganese refinery. The least absolute shrinkage and selection operator (LASSO) followed by the generalized linear model (GLM) was applied to evaluate the associations between multiple-metal exposures and TL change (TL in 2012 minus TL in 2020). Bayesian kernel machine regression (BKMR) was applied to cope with metal mixtures and evaluate their joint effects on TL change. Among men, three statistical methods consistently showed rubidium was negatively associated with TL change (ß [95% CI] = -2.755 [-5.119, -0.391] in the GLM) and dominated the negative overall effects of 10 metal mixtures (magnesium, manganese, iron, cobalt, copper, zinc, selenium, rubidium, cadmium, and lead) on TL change (posterior inclusion probabilities = 0.816). Among women, the GLM (ß [95% CI] = 4.463 [0.943, 7.983]) and LASSO (ß = 4.289) showed rubidium was positively associated with TL change. Interestingly, no significant association was observed between exposure to metal mixtures and TL change in overall participants (P > 0.05). Furthermore, stratified analysis showed significant relationships between rubidium and TL change in men (ß = -2.744), women (ß = 3.624), and current smokers (ß = -3.266) (both P interaction <0.05). In summary, our findings underlined the steady and negative association between rubidium and TL change among men with potential sex-dependent heterogeneities. Further experimental studies are required to expound the underlying mechanisms.

MR-PheWAS for the causal effects of serum magnesium on multiple disease outcomes in Caucasian descent.

Magnesium is integral to many physiological processes, whereas variations in its levels, even within the normal range, can have critical implications for health. To explore the broad clinical effects of varying serum magnesium levels, we performed a two-sample Mendelian randomization and phenome-wide association study (MR-PheWAS) in the UK Biobank cohort. In total, MR-PheWAS analysis implicated a causal role of serum magnesium levels in five disease groups and six disease outcomes. In addition, our study indicated the gender-specific effects of nine disease groups/outcomes in MR estimated effects. The protein-protein interaction network demonstrated an interaction between the serum magnesium-associated gene DCDC1 and the cataract- associated gene PAX6. The present study verified several previously reported disease outcomes and identified novel potential disease outcomes for serum magnesium levels. The DCDC1 gene and the PAX6 gene may be the new targets for promoting the treatments of cataracts using magnesium intervention.

Combined exposure to multiple metals and cognitive function in older adults.

Single toxic metal exposure has been reported to be associated with impaired cognitive function, but less is known about the effects of combined exposure to multiple metals. The aim of the study was to investigate the potential associations and interactions of multiple metals with cognitive function in older adults using multi-pollutants approach. A cross-sectional study was conducted in a total of 2879 participants aged ≥ 60 years old. We systematically measured levels of 22 blood metals and used the Mini-Mental State Examination (MMSE) to assess the cognitive function. The least absolute shrinkage and selection operator (LASSO) penalized regression was applied to identify independently main metals. Adjusted estimates of cognitive function with selected metals were investigated by generalized linear regression in the multi-metal model. We found that calcium, titanium, vanadium, copper, zinc, arsenic, selenium, rubidium, molybdenum, cadmium, barium, and lead were independently identified based on LASSO penalized regression. The multi-metal model showed a higher MMSE of 0.384 (95% CI: 0.122-0.646) for a 1-SD increment in log-transformed rubidium and a lower MMSE of 0.460 (95% CI: - 0.706 to - 0.214) for a 1-SD increment in log-transformed cadmium (P < 0.05). The significantly negative associations between cadmium and cognitive function were attenuated to null accompanying with increasing concentrations of rubidium (P interaction = 0.256). Our findings suggested that blood rubidium and cadmium were mainly associated with cognitive function when accounting for co-exposure to other metals and higher level of rubidium appeared to attenuate the toxic effects of cadmium on cognitive function in older adults.

Mediated relationships between multiple metals exposure and fasting blood glucose by reproductive hormones in Chinese men.

Previous studies have reported metals exposure contribute to the change of fasting blood glucose (FBG) level. However, the roles of reproductive hormones in their associations have not been fully elucidated. The aim of the study is to investigate the associations of multiple serum metals with reproductive hormones, and to further explore potential roles of reproductive hormones in relationships between metals exposure and FBG level. A total of 1911 Chinese Han men were analyzed by a cross-sectional study. We measured serum levels of 22 metals by inductively coupled plasma mass spectrometer (ICP-MS). FBG, total testosterone (TT), estradiol (E2), follicle stimulating hormone (FSH), and sex hormone-binding globulin (SHBG) levels were determined. Least absolute shrinkage and selection operator (LASSO) regression models were conducted to select important metals, and restricted cubic spline models were then used to estimate dose-response relationships between selected metals and reproductive hormones. We also conducted mediation analyses to evaluate whether reproductive hormones played mediating roles in the associations between metals and FBG. We found significant inverse dose-dependent trends of copper, tin and zinc with E2; zinc with SHBG; copper and nickel with TT, while significant positive dose-dependent trend of iron with E2, respectively. Moreover, approximately inverted U-shaped associations existed between lead and SHBG, iron and TT. In addition, E2, SHBG and TT were negatively associated with FBG level. In mediation analyses, the association of copper with FBG was mediated by E2 and TT, with a mediation ratio of 10.4% and 22.1%, respectively. Furthermore, E2 and SHBG mediated the relationship of zinc with FBG, with a mediation ratio of 7.8% and 14.5%, respectively. E2 mediated 11.5% of positive relationship between tin with FBG. Our study suggested that the associations of metals exposure with FBG may be mediated by reproductive hormones.

Sex-specific associations of plasma metals and metal mixtures with glucose metabolism: An occupational population-based study in China.

Studies with multi-pollutant approach on the relationships between multiple metals and fasting plasma glucose (FPG) are limited. Few studies are available on the potential sex-specific associations between metal exposures and glucose metabolism. We explored the associations between 22 plasma metals and FPG level among the 769 participants from the manganese-exposed workers healthy cohort in China. We applied a sparse partial least squares (sPLS) regression followed by ordinary least-squares regression to evaluate multi-pollutant association. Bayesian kernel machine regression (BKMR) model was used to deal with metal mixtures and evaluate their joint effects on FPG level. In the sPLS model, negative associations on FPG levels were observed for plasma iron (belta = -0.066), cobalt (belta = -0.075), barium (belta = -0.109), and positive associations for strontium (belta = 0.082), and selenium (belta = 0.057) in men, which overlapped with the results among the overall participants. Among women, plasma copper (belta = 0.112) and antimony (belta = 0.137) were positively associated with elevated FPG level. Plasma magnesium was negatively associated with FPG level in both sexes (belta = -0.071 in men and belta = -0.144 in women). The results of overlapped for plasma magnesium was selected as the significant contributor to decreasing FPG level in the multi-pollutant, single-metal, and multi-metal models. BKMR model showed a significantly negative over-all effect of six metal mixtures (magnesium, iron, cobalt, selenium, strontium and barium) on FPG level among the overall participants from all the metals fixed at 50th percentile. In summary, our findings underline the probable role of metals in glucose homeostasis with potential sex-dependent heterogeneities, and suggest more researches are needed to explore the sex-specific associations of metal exposures with risk of diabetes.

The associations of multiple metals mixture with accelerated DNA methylation aging.

Aging is a leading cause of mortality for the elderly and DNA methylation age is reported to be predictive of biological aging. However, few studies have investigated the associations between multiple metals exposure and accelerated aging in the elderly. We performed a pilot study of 288 elderly participants aged 50-115 years and measured genome-wide DNA methylation and 22 blood metals concentrations. Measures of DNA methylation age were estimated using CpGs from Illumina HumanMethylation EPIC BeadChip. Linear mixed regression and Bayesian kernel machine regression (BKMR) models were used to estimate the individual and overall associations between multiple metals and accelerated methylation aging. Single metal models revealed that each 1-standard deviance (SD) increase in log-transformed vanadium, cobalt, nickel, zinc, arsenic, and barium was associated with a -2.256, -1.318, 1.004, -1.926, 1.910 and -1.356 changes in ΔAge, respectively; meanwhile, for aging rate, the change was -0.019, -0.013, 0.010, -0.018, 0.023, and -0.012, respectively (all P < 0.05). The BKMR models showed reverse U-shaped associations of the overall metals mixture with ΔAge and aging rate. Downward trends of ΔAge and aging rate were observed for increasing quantiles of essential metals mixture, but upward trends were observed for non-essential metals mixture. Further individual analysis of the BKMR revealed that the 95% confidence interval of ΔAge and aging rate associated with vanadium, zinc, and arsenic did not cross 0, when other metals concentrations set at 25th, 50th, and 75th percentile. Our findings suggest reverse U-shaped associations of the overall metals mixture with accelerated methylation aging for the first time, and vanadium, zinc, and arsenic may be major contributors driving the associations.

Associations of prenatal exposure to multiple metals with testicular volume and anogenital distance in infant boys: A longitudinal cohort study.

BACKGROUND: Human are widely exposed to multiple metals, some of which have suspected reproductive toxicity, but no human studies have investigated the developmental effects of prenatal metal exposure. OBJECTIVES: We aimed to evaluate the associations between prenatal multiple metal exposure and reproductive development in boys at 2-3 years using multi-pollutant approach. METHODS: This prospective study used data of 564 mother-child pairs recruited from the Guangxi Birth Cohort Study. Twenty serum metal concentrations were measured. Least absolute shrinkage and selection operator (LASSO) penalized regression was used to identify independent associations between prenatal multiple metal exposure and testicular volume (TV), and anogenital distance (AGD). Adjusted estimates were then obtained using multiple linear regression analysis, and the regression tree method was used to explore the interactions. RESULTS: Boys in the highest quartile of prenatal lead exposure had a 0.064 mL (95% CI: -0.124, -0.004) smaller ln-transformed TV, 0.060 cm (95% CI: -0.110, -0.011) shorter ln-transformed anopenile distance (AGDap), and 0.115 cm (95% CI: -0.190, -0.039) shorter ln-transformed anoscrotal distance (AGDas) than boys in the lowest quartile (all Ptrend < 0.05). Chromium was inversely with ln-transformed AGDap (ß = -0.078, 95% CI: -0.127, -0.030) and ln-transformed AGDas (ß = -0.113, 95% CI: -0.188, -0.038), while stibium was positivity associated with ln-transformed AGDap (ß = 0.091, 95% CI: 0.046, 0.136) and strontium was positivity associated with ln-transformed AGDas (ß = 0.120, 95% CI: 0.051, 0.189) (all Ptrend < 0.05). And the critical window of vulnerability may be the late pregnancy (the second and third trimester). Moreover, we detected interaction effects between lead, chromium and stibium on AGDap; lead, chromium and strontium on AGDas. CONCLUSIONS: The results suggest that prenatal exposure to lead, chromium, stibium and strontium may affect TV and/or AGD in infant boys. Potential mechanisms for the complex metal interactive effects during vulnerable periods are worthy of further investigation.

Occupational exposure to manganese and risk of creatine kinase and creatine kinase-MB elevation among ferromanganese refinery workers.

BACKGROUND: Elevated exposure to manganese (Mn) could induce cardiovascular dysfunction. However, limited research is available on the effects of occupational Mn exposure on myocardial enzymes. We aimed to evaluate the relationships between Mn exposure and myocardial enzyme elevation among Mn-exposed workers. METHODS: Data were from a follow-up investigation of a Mn-exposed workers healthy cohort in 2017. A total of 744 workers were divided into low-exposure and high-exposure groups according to Mn time-weighted average (Mn-TWA) of less than or equal to 0.15 mg/m3 or greater than 0.15 mg/m3 , respectively. Serum levels of myocardial enzymes, including creatine kinase (CK) and creatine kinase-MB (CK-MB), lactic dehydrogenase, α-hydroxybutyrate dehydrogenase, and aspartate transaminase, were assessed, as well as airborne Mn exposure levels. RESULTS: After adjustment for sex, body mass index, seniority, education, current smoking status, current drinking status, and hypertension, Mn-TWA levels were positively associated with the risk of CK elevation (odds ratio [OR] = 1.47 (95% confidence interval [CI]: 1.18-1.83) per interquartile range [IQR] increase), and risk of CK-MB elevation [OR = 1.57 (95% CI: 1.03-2.38) per IQR increase]. In a stratified analysis, Mn-TWA levels were positively correlated with CK elevation in workers of seniority greater than 19 years, male workers, current smokers, and drinkers. CONCLUSION: Our results suggest that occupational exposure to Mn is associated with increased risk of CK and CK-MB elevation. The potential mechanisms of the associations between airborne exposure to Mn and increased risk of these myocardial enzyme elevations warrant further investigation.

Associations between manganese exposure and multiple immunological parameters in manganese-exposed workers healthy cohort.

BACKGROUND: Manganese (Mn) ions play a crucial role in the immune response. The immunotoxicity of Mn is rarely reported compared with the neurotoxicity of Mn. OBJECTIVES: The purpose of this study was to investigate the associations between chronic Mn exposure and immunological parameters in occupational Mn-exposed workers. METHODS: A total of 538 workers were selected from the follow-up of manganese-exposed workers healthy cohort (MEWHC) in 2017. We divided the workers into the low-exposure group and the high-exposure group by the cutoff of the manganese-time weighted average (Mn-TWA) setting at 0.15 mg/m3. We examined serum immunological parameters by the immunoturbidimetric method and leukocyte counts and ratios in blood routine. Then we used the generalized linear model analyses and spline analyses to explore the associations between external exposure of Mn and multiple immunological parameters adjusted for variables. Based on the epidemiological analyses, we used Elisa (enzyme-linked immune sorbent assay) to detect plasma complement C3 of Mn-exposed rats. RESULTS: In male workers, the mean value of complement C3 was 1.20 ±â€¯0.16 g/L in the high-exposure group, which was significantly lower as compared to the low-exposure group (1.25 ±â€¯0.18 g/L, P = 0.023). The generalize linear models' analyses showed that complement C3 value had a significantly negative association with external exposure of Mn included adjustment for variables (ß = -0.04, P = 0.035). Moreover, in male rats, the high-exposure group also had a lower level of complement C3 compared with the low-exposure group (P < 0.001). None significant association was observed in immunological parameters among female workers and rats (all P > 0.05). CONCLUSIONS: Mn exposure from inhalable dust was associated with decreased complement C3 among occupationally Mn-exposed male individuals but not in female workers, which was further confirmed by the rat model. Further research into the possible mechanism of C3 reduction is needed in the future.

Plasma metals and serum bilirubin levels in workers from manganese-exposed workers healthy cohort (MEWHC).

Few studies specifically address the possible associations between multiple-metal exposures and liver damage among the occupational population. This study aimed to explore the cross-sectional relationships of plasma metals with liver function parameters. For 571 on-the-spot workers in the manganese-exposed workers healthy cohort (MEWHC), we determined liver function parameters: total bilirubin (TBILI), direct bilirubin (DBILI), indirect bilirubin (IBILI), alanine transaminase (ALT) and aspartate transaminase (AST). Total concentrations of 22 plasma metals were measured by ICP-MS. The LASSO (least absolute shrinkage and selection operator) penalized regression model was applied for selecting plasma metals independently associated with liver function parameters. Multiple linear regression analyses and restricted cubic spline (RCS) were utilized for identifying the exposure-response relationship of plasma metals with liver function parameters. After adjusting for covariates and selected metals, a 1-SD increase in log-10 transformed levels of iron was associated with increases in the levels of TBILI, DBILI and IBILI by 20.3%, 12.1% and 23.7%, respectively; similar increases in molybdenum for decreases in levels of TBILI, DBILI and IBILI by 6.1%, 2.6% and 8.3%, respectively. The effect of a 1-SD increase in plasma copper corresponded decreases of 3.2%, 3.4% and 5.0% in TBILI, AST and ALT levels, respectively. The spline analyses further clarified the non-linear relationships between plasma iron and bilirubin whilst negative linear relationships for plasma molybdenum and bilirubin. Plasma iron was positively whilst plasma molybdenum was negatively associated with increased serum bilirubin levels. Further studies are needed to validate these associations and uncover the underlying mechanisms.