A Quantile-Based g-Computation Approach to Addressing the Effects of Exposure Mixtures

Physical activity modified association of urinary metals mixture and fasting blood glucose in children: From two panel studies

There is unclear evidence available on the associations between multiple metals and fasting blood glucose (FBG) in children, and whether they could be beneficial from physical activity. We included 283 children aged 4-12 years from two panel studies with 4-consecutive morning urinary 13 essential metals and 10 non-essential metals repeated across 3 seasons. We employed multiple informant model, linear mixed-effect model, and quantile g-computation to evaluate associations of single metal and their mixture with FBG and interactions with extra-school activity. The results showed that positive relations of multiple essential metals (aluminum, chromium, copper, iron, molybdenum (Mo), nickel, selenium (Se), strontium, zinc) and non-essential metals (arsenic (As), cadmium (Cd), rubidium, titanium (Ti), thallium) with FBG were the strongest at lag 0 (the health examination day), especially in overweight & obesity children (FDR <0.05). The strongest effect presented 1-fold increment in As was related to FBG increased 1.66% (95%CI: 0.84%, 2.48%) in overweight & obesity children. Notably, modification of extra-school activity showed significant, and the effects of multiple metals on FBG were attenuated in children taking total extra-school activity ≥1 h/day, and only one type of which, low or moderate & high intensity extra-school activity reached 20 min/day (Pint <0.05). For instance, each 1-fold increased As was associated with 1.41% increased FBG in overall children taking total extra-school activity <1 h/day, while that of 0.13% in those ≥1 h/day. Meanwhile, mixture of all, essential and non-essential metals were associated with increased FBG, a trend that decreased and became nonsignificant in children having certain extra-school activity, which were dominated by Mo, Se, Ti, Cd. And such relations were substantially beneficial from extra-school activity in overweight & obesity children. Accordingly, multiple essential and non-essential metals, both individual and in mixture, were positively related to FBG in children, which might be attenuated by regular physical activity.

Combined Exposures and Mixtures Research: An Enduring NIEHS Priority

BACKGROUND

The National Institute of Environmental Health Sciences (NIEHS) continues to prioritize research to better understand the health effects resulting from exposure to mixtures of chemical and nonchemical stressors. Mixtures research activities over the last decade were informed by expert input during the development and deliberations of the 2011 NIEHS Workshop "Advancing Research on Mixtures: New Perspectives and Approaches for Predicting Adverse Human Health Effects." NIEHS mixtures research efforts since then have focused on key themes including a) prioritizing mixtures for study, b) translating mixtures data from in vitro and in vivo studies, c) developing cross-disciplinary collaborations, d) informing component-based and whole-mixture assessment approaches, e) developing sufficient similarity methods to compare across complex mixtures, f) using systems-based approaches to evaluate mixtures, and g) focusing on management and integration of mixtures-related data.

OBJECTIVES

We aimed to describe NIEHS driven research on mixtures and combined exposures over the last decade and present areas for future attention.

RESULTS

Intramural and extramural mixtures research projects have incorporated a diverse array of chemicals (e.g., polycyclic aromatic hydrocarbons, botanicals, personal care products, wildfire emissions) and nonchemical stressors (e.g., socioeconomic factors, social adversity) and have focused on many diseases (e.g., breast cancer, atherosclerosis, immune disruption). We have made significant progress in certain areas, such as developing statistical methods for evaluating multiple chemical associations in epidemiology and building translational mixtures projects that include both in vitro and in vivo models.

DISCUSSION

Moving forward, additional work is needed to improve mixtures data integration, elucidate interactions between chemical and nonchemical stressors, and resolve the geospatial and temporal nature of mixture exposures. Continued mixtures research will be critical to informing cumulative impact assessments and addressing complex challenges, such as environmental justice and climate change. https://doi.org/10.1289/EHP14340.

Metabolomic landscape of overall and common cancers in the UK Biobank: A prospective cohort study

Information about the NMR metabolomics landscape of overall, and common cancers is still limited. Based on a cohort of 83,290 participants from the UK Biobank, we used multivariate Cox regression to assess the associations between each of the 168 metabolites with the risks of overall cancer and 20 specific types of cancer. Then, we applied LASSO to identify important metabolites for overall cancer risk and obtained their associations using multivariate cox regression. We further conducted mediation analysis to evaluate the mediated role of metabolites in the effects of traditional factors on overall cancer risk. Finally, we included the 13 identified metabolites as predictors in prediction models, and compared the accuracies of our traditional models. We found that there were commonalities among the metabolic profiles of overall and specific types of cancer: the top 20 frequently identified metabolites for 20 specific types of cancer were all associated with overall cancer; most of the specific types of cancer had common identified metabolites. Meanwhile, the associations between the same metabolite with different types of cancer can vary based on the site of origin. We identified 13 metabolic biomarkers associated with overall cancer, and found that they mediated the effects of traditional factors. The accuracies of prediction models improved when we added 13 identified metabolites in models. This study is helpful to understand the metabolic mechanisms of overall and a wide range of cancers, and our results also indicate that NMR metabolites are potential biomarkers in cancer diagnosis and prevention.

Phthalate mixtures and insulin resistance: an item response theory approach to quantify exposure burden to phthalate mixtures

BACKGROUND

Molar sums are often used to quantify total phthalate exposure, but they do not capture patterns of exposure to multiple phthalates.

OBJECTIVE

To introduce an exposure burden score method for quantifying exposure to phthalate metabolites and examine the association between phthalate burden scores and Homeostatic Model Assessment for Insulin Resistance (HOMA-IR).

METHODS

We applied item response theory (IRT) to data from 3474 adults aged 20-60 years in the 2013-2018 National Health and Examination Survey (NHANES) to quantify latent phthalate exposure burden from 12 phthalate metabolites. We compared model fits of three IRT models that used different a priori groupings (general phthalate burden; low molecular weight (LMW) and high molecular weight (HMW) burdens; and LMW, HMW and DEHP burden), and used the best fitting model to estimate phthalate exposure burden scores. Regression models assessed the covariate-adjusted association between phthalate burden scores and HOMA-IR. We compared findings to those using molar sums. In secondary analyses, we examined how the IRT model could be used for data harmonization when a subset of participants are missing some phthalate metabolites, and accounted for measurement error of the phthalate burden scores in estimating associations with HOMA-IR through a resampling approach using plausible value imputation.

RESULTS

A three correlated factors model (LMW, HMW and DEHP burdens) provided the best fit. One interquartile range (IQR) increase in DEHP burden score was associated with 0.094 (95% CI: 0.022, 0.164, p = 0.010) increase in log HOMA-IR, co-adjusted for LMW and HMW burden scores. Findings were consistent when using log molar sums. Associations of phthalate burden and insulin resistance were also consistent when participants were simulated to be missing some phthalate metabolites, and when we accounted for measurement error in estimating burden scores.

CONCLUSION

Both phthalate molar sums and burden scores are sensitive to associations with insulin resistance. Phthalate burden scores may be useful for data harmonization.

Associations between short-term exposure to air pollution and thyroid function in a representative sample of the Korean population

BACKGROUND

Disruption of thyroid function can profoundly affect various organ systems. However, studies on the association between air pollution and thyroid function are relatively scarce and most studies have focused on the long-term effects of air pollution among pregnant women.

OBJECTIVES

This study aimed to explore the associations between short-term exposure to air pollution and thyroid function in the general population.

METHODS

Data from the Korea National Health and Nutrition Examination Survey (2013-2015) were analyzed (n = 5,626). Air pollution concentrations in residential addresses were estimated using Community Multiscale Air Quality models. The moving averages of air pollution over 7 days were set as exposure variables through exploratory analyses. Linear regression and quantile g-computation models were constructed to assess the effects of individual air pollutants and air pollution mixture, respectively.

RESULTS

A 10-ppb increase in NO2 (18.8-μg/m3 increase) and CO (11.5-μg/m3 increase) was associated with 2.43% [95% confidence interval (CI): 0.42, 4.48] and 0.19% (95% CI: 0.01, 0.36) higher thyroid-stimulating hormone (TSH) levels, respectively. A 10-μg/m3 increase in PM2.5 and a 10-ppb increase in O3 (19.6-μg/m3 increment) were associated with 0.87% (95% CI: 1.47, -0.27) and 0.59% (95% CI: 1.18, -0.001) lower free thyroxine (fT4) levels, respectively. A simultaneous quartile increase in PM2.5, NO2, O3, and CO levels was associated with lower fT4 but not TSH levels.

CONCLUSIONS

As the subtle changes in thyroid function can affect various organ systems, the present results may have substantial public health implications despite the relatively modest effect sizes. Because this was a cross-sectional study, it is necessary to conduct further experimental or repeated-measures studies to consolidate the current results.

The association between PM2.5 components and blood pressure changes in late pregnancy: A combined analysis of traditional and machine learning models

BACKGROUND

PM2.5 is a harmful mixture of various chemical components that pose a challenge in determining their individual and combined health effects due to multicollinearity issues with traditional linear regression models. This study aimed to develop an analytical methodology combining traditional and novel machine learning models to evaluate PM2.5's combined effects on blood pressure (BP) and identify the most toxic components.

METHODS

We measured late-pregnancy BP of 1138 women from the Heshan cohort while simultaneously analyzing 31 PM2.5 components. We utilized multiple linear regression modeling to establish the relationship between PM2.5 components and late-pregnancy BP and applied Random Forest (RF) and generalized Weighted Quantile Sum (gWQS) regression to identify the most toxic components contributing to elevated BP and to quantitatively evaluate the cumulative effect of the PM2.5 component mixtures.

RESULTS

The results revealed that 16 PM2.5 components, such as EC, OC, Ti, Fe, Mn, Cu, Cd, Mg, K, Pb, Se, Na+, K+, Cl-, NO3-, and F-, contributed to elevated systolic blood pressure (SBP), while 26 components, including two carbon components (EC, OC), fourteen metallics (Ti, Fe, Mn, Cr, Mo, Co, Cu, Zn, Cd, Na, Mg, Al, K, Pb), one metalloid (Se), and nine water-soluble ions (Na+, K+, Mg2+, Ca2+, NH4+, Cl-, NO3-, SO42-, F-), contributed to elevated diastolic blood pressure (DBP). Mn and Cr were the most toxic components for elevated SBP and DBP, respectively, as analyzed by RF and gWQS models and verified against each other. Exposure to PM2.5 component mixtures increased SBP by 1.04 mmHg (95% CI: 0.33-1.76) and DBP by 1.13 mmHg (95% CI: 0.47-1.78).

CONCLUSIONS

Our study highlights the effectiveness of combining traditional and novel models as an analytical strategy to quantify the health effects of PM2.5 constituent mixtures.

Association between multipollutant exposure and thyroid hormones in elderly people: A cross-sectional study in China

Environmental chemicals have been indicated to cause disruption of thyroid homeostasis in human populations. However, previous studies mostly focused on single group of chemicals. Herein, we investigate the independent and combined effects of multiple pollutants on thyroid homeostasis, including thyroid-stimulating hormone (TSH), total and free thyroxine (tT4 and fT4) and total and free triiodothyronine (tT3 and fT3) in elderly people. These environmental pollutants (n = 144) are from ten categories, including phenols, parabens, perfluoroalkyl substances (PFASs), polychlorinated biphenyls (PCBs), phthalate esters (PAEs), polycyclic aromatic hydrocarbons (PAHs), organochlorine pesticides (OCPs), organophosphate pesticides (OPPs), synthetic pyrethroids (SPs), herbicides, and metals. Few studies have evaluated the health risks of these 144 chemicals, especially their joint effects. In single-pollutant evaluations, multiple linear regression (MLR) models were used to estimate the independent associations between multiple exposures and thyroid biomarkers. In multi-pollutant evaluations, elastic net regression and Bayesian kernel machine regression (BKMR) models were used to estimate the combined associations. The MLR models showed that 41 chemicals were significantly related to THs levels. BKMR models revealed the most important chemical groups: metals for TSH, PAHs, SPs and PCBs for tT4, herbicides and SPs for tT3. This study will contribute to the understanding of multipollutant exposure and help prioritize specific chemical groups related to thyroid hormone disruption.

Reduction of haemoglobin is related to metal mixtures exposure in Chinese preschoolers: Joint effect models

Heavy metal exposure is a known risk factor for hematologic disorders in children, yet the impact of co-exposure to multiple metals remains underexplored. This cross-sectional study investigates the relationship between urinary levels of 23 metals and haemoglobin (Hb) in 1460 Chinese preschoolers. The concentrations of the 23 urinary metals were quantified using an inductively coupled plasma mass spectrometer, while Hb levels were assessed through finger prick blood samples. To evaluate the co-exposure effects, we employed three approaches: Generalized linear regression model, joint effect models including Quantile g-Computation and Bayesian Kernel Machine Regression (BKMR). From the generalized linear regression and Quantile g-computation, urinary uranium, thallium, aluminium, iron and tungsten were correlated negatively with Hb, while urinary barium was correlated positively (all P < 0.05). Moreover, significant negative associations between metal mixtures exposure with Hb were identified in both Quantile g-computation [β (95% CI): -0.083 (-0.132, -0.033), P = 0.0012] and BKMR [90th percentile vs. 50th percentile β (95% CI): -0.238 (-0.368, -0.107), P < 0.001] with aluminium emerging as the primary contributor to this joint effect (weight in Quantile g-computation = 0.399, PIPs in BKMR = 0.896). These findings provide a potential explanation for environmental exposure to metals and Hb-related disease in preschoolers.

Association between PM2.5 constituents and cardiometabolic risk factors: Exploring individual and combined effects, and mediating inflammation

BACKGROUND

The individual and combined effects of PM2.5 constituents on cardiometabolic risk factors are sparsely investigated. Besides, the key cardiometabolic risk factor that PM2.5 constituents targeted and the biological mechanisms remain unclear.

METHOD

A multistage, stratified cluster sampling survey was conducted in two typically air-polluted Chinese cities. The PM2.5 and its constituents including sulfate, nitrate, ammonium, organic matter, and black carbon were predicted using a machine learning model. Twenty biomarkers in three category were simultaneously adopted as cardiometabolic risk factors. We explored the individual and mixture association of long-term PM2.5 constituents with these markers using generalized additive model and quantile-based g-computation, respectively. To minimize potential confounding effects, we accounted for covariates including demographic, lifestyle, meteorological, temporal trends, and disease-related information. We further used ROC curve and mediation analysis to identify the key subclinical indicators and explore whether inflammatory mediators mediate such association, respectively.

RESULT

PM2.5 constituents was positively correlated with HOMA-B, TC, TG, LDL-C and LCI, and negatively correlated with PP and RC. Further, PM2.5 constituent mixture was positive associated with DBP, MAP, HbA1c, HOMA-B, AC, CRI-1 and CRI-2, and negative associated with PP and HDL-C. The ROC analysis further reveals that multiple cardiometabolic risk factors can collectively discriminate exposure to PM2.5 constituents (AUC>0.9), among which PP and CRI-2 as individual indicators exhibit better identifiable performance for nitrate and ammonium (AUC>0.75). We also found that multiple blood lipid indicators may be affected by PM2.5 and its constituents, possibly mediated through complement C3 or hsCRP.

CONCLUSION

Our study suggested associations of individual and combined PM2.5 constituents exposure with cardiometabolic risk factors. PP and CRI-2 were the targeted markers of long-term exposure to nitrate and ammonium. Inflammation may serve as a mediating factor between PM2.5 constituents and dyslipidemia, which enhance current understanding of potential pathways for PM2.5-induced preclinical cardiovascular responses.

Power Analysis of Exposure Mixture Studies via Monte Carlo Simulations

Estimating sample size and statistical power is an essential part of a good epidemiological study design. Closed-form formulas exist for simple hypothesis tests but not for advanced statistical methods designed for exposure mixture studies. Estimating power with Monte Carlo simulations is flexible and applicable to these methods. However, it is not straightforward to code a simulation for non-experienced programmers and is often hard for a researcher to manually specify multivariate associations among exposure mixtures to set up a simulation. To simplify this process, we present the R package mpower for power analysis of observational studies of environmental exposure mixtures involving recently-developed mixtures analysis methods. The components within mpower are also versatile enough to accommodate any mixtures methods that will developed in the future. The package allows users to simulate realistic exposure data and mixed-typed covariates based on public data set such as the National Health and Nutrition Examination Survey or other existing data set from prior studies. Users can generate power curves to assess the trade-offs between sample size, effect size, and power of a design. This paper presents tutorials and examples of power analysis using mpower.

Estimated daily intake and cumulative risk assessment of organophosphate esters and associations with DNA damage among e-waste workers in Hong Kong

Organophosphate esters (OPEs) are extensively used as additives in various products, including electronic equipment, which becomes e-waste when obsolete. Nevertheless, no study has evaluated OPEs exposure levels and the related health risks among e-waste workers in Hong Kong. Therefore, 201 first-spot morning urine samples were collected from 101 e-waste workers and 100 office workers to compare eight urinary OPE metabolites (mOPEs) levels in these groups. The concentrations of six mOPEs were similar in e-waste workers and office workers, except for significantly higher levels of diphenyl phosphate (DPHP) in e-waste workers and bis(1-chloro-2propyl) phosphate (BCIPP) in office workers. Spearman correlation analysis showed that most non-chlorinated mOPEs were correlated with each other in e-waste workers (i.e., nine out of ten pairs, including di-p-cresyl phosphate (DpCP) and di-o-cresyl phosphate (DoCP), DpCP and bis(2-butoxyethyl) phosphate (BBOEP), DpCP and DPHP, DpCP and dibutyl phosphate (DBP), DoCP and BBOEP, DoCP and DPHP, DoCP and DBP, BBOEP and DPHP, DPHP and DBP), indicating that handling e-waste could be the exposure source of specific OPEs. The median values of estimated daily intake (EDI) and hazard quotient (HQ) suggested that the health risks from OPEs exposures were under the recommended thresholds. However, linear regression models, Quantile g-computation, and Bayesian kernel machine regression found that urinary mOPEs elevated 8-hydroxy-2-deoxyguanosine (8-OhdG) levels individually or as a mixture, in which DPHP contributed prominently. In conclusion, although e-waste might not elevate the internal OPEs levels among the participating Hong Kong e-waste workers, attention should be paid to the potential DNA damage stimulated by OPEs under the currently recommended thresholds.

Prenatal per- and polyfluoroalkyl substances (PFAS) and maternal oxidative stress: Evidence from the LIFECODES study

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are synthetic chemicals linked to adverse pregnancy outcomes. Although their underlying biological mechanisms are not fully understood, evidence suggests PFAS may disrupt endocrine functions and contribute to oxidative stress (OS) and inflammation.

OBJECTIVE

We examined associations between early pregnancy PFAS exposure and OS biomarkers, exploring potential effect modifications by fetal sex and maternal race.

METHODS

We used data from 469 LIFECODES participants with measured plasma PFAS (median 10 weeks gestation) and repeated measures (median 10, 18, 26, and 35 weeks gestation) of urinary OS biomarkers [8-iso-prostaglandin-F2α (8-isoprostane) and 8-hydroxydeoxyguanosine (8-OHdG)]. Protein damage biomarkers (chlorotyrosine, dityrosine, and nitrotyrosine) were additionally measured in plasma from a subset (N = 167) during the third visit. Associations between each PFAS and OS biomarkers were examined using linear mixed-effects models and multivariable linear regressions, adjusting for potential confounders, including maternal age, race, education level, pre-pregnancy BMI, insurance status, and parity. Effect modifications were evaluated by including an interaction term between each PFAS and fetal sex or maternal race in the models.

RESULTS

We observed significant positive associations between PFOS and 8-isoprostane, with a 9.68% increase in 8-isoprostane levels (95% CI: 0.10%, 20.18%) per interquartile range increase in PFOS. In contrast, PFUA was negatively associated [9.32% (95% CI: -17.68%, -0.11%)], while there were suggestive positive associations for MPAH and PFOA with 8-isoprostane. The associations of several PFAS with 8-OHdG varied by fetal sex, showing generally positive trends in women who delivered females, but negative or null in those who delivered males. No significant effect modification by maternal race was observed.

CONCLUSIONS

This study provides evidence linking PFAS exposure to OS during pregnancy, with potential sex-specific effects of certain PFAS on 8-OHdG. Further research should explore additional OS/inflammatory biomarkers and assess the modifying effects of dietary and behavioral patterns across diverse populations.

Prenatal exposure to persistent and non-persistent chemical mixtures and associations with adverse birth outcomes in the Atlanta African American Maternal-Child Cohort

BACKGROUND

African Americans (AAs) experience higher rates of preterm birth and fetal growth restriction relative to other pregnant populations. Differential in utero exposure to environmental chemicals may partially explain these health disparities, as AAs are disproportionately exposed to environmental hazards.

OBJECTIVE

We examined the individual and mixture effects of non-persistent chemicals and persistent organic pollutants (POPs) on gestational age at birth and birthweight for gestational age z-scores within a prospective cohort of pregnant AAs.

METHODS

First-trimester serum and urine samples obtained from participants within the Atlanta African American Maternal-Child cohort were analyzed for 43 environmental chemicals, including per-and polyfluoroalkyl substances (PFAS), polybrominated diphenyl ethers (PBDEs), organochlorine pesticides, pyrethroid insecticides, phthalates, bisphenol A, nicotine, and the primary metabolite of delta-9-tetrahydrocannabinol. Linear regression was used to estimate individual associations between chemicals and gestational age and birthweight z-scores (N ranging from 107 to 523). Mixture associations were estimated using quantile g-computation, principal component (PC) analyses, and hierarchical Bayesian kernel machine regression among complete cases (N = 86).

RESULTS

Using quantile g-computation, increasing all chemical exposures by one quantile was modestly associated with a reduction in gestational age (mean change per quartile increase = -0.47, 95% CI = -1.56, 0.61) and birthweight z-scores (mean change per quartile increase = -0.49, 95% CI = -1.14, 0.15). All PCs were associated with a reduction in birthweight z-scores; associations were greatest in magnitude for the two PCs reflecting exposure to combined tobacco, insecticides, PBDEs, and phthalates. In single pollutant models, we observed inconsistent and largely non-significant associations.

SIGNIFANCE

We conducted multiple targeted exposure assessment methods to quantify levels of environmental chemicals and leveraged mixture methods to quantify their joint effects on gestational age and birthweight z-scores. Our findings suggest that prenatal exposure to multiple classes of persistent and non-persistent chemicals is associated with reduced gestational age and birthweight z-scores in AAs.

IMPACT

African Americans (AAs) experience higher rates of preterm birth and fetal growth restriction relative to other pregnant populations. Differential in utero exposure to environmental chemicals may partially explain these health disparities, as AAs are disproportionately exposed to environmental hazards. In the present study, we analyzed serum and urine samples for levels of 43 environmental chemicals. We used quantile g-computation, principal component analysis, and BKMR to assess associations between chemical exposure mixtures and adverse birth outcomes. Our findings suggest that prenatal exposure to multiple classes of chemicals is associated with reduced birthweight z-scores, a proxy for fetal growth, in AAs.

Metal mixture exposures and serum lipid levels in childhood: the Rhea mother-child cohort in Greece

BACKGROUND

Growing evidence suggests that cardiovascular disease develops over the lifetime, often beginning in childhood. Metal exposures have been associated with cardiovascular disease and important risk factors, including dyslipidemia, but prior studies have largely focused on adult populations and single metal exposures.

OBJECTIVE

To investigate the individual and joint impacts of multiple metal exposures on lipid levels during childhood.

METHODS

This cross-sectional study included 291 4-year-old children from the Rhea Cohort Study in Heraklion, Greece. Seven metals (manganese, cobalt, selenium, molybdenum, cadmium, mercury, and lead) were measured in whole blood using inductively coupled plasma mass spectrometry. Serum lipid levels included total cholesterol, triglycerides, high-density lipoprotein (HDL) cholesterol, and low-density lipoprotein (LDL) cholesterol. To determine the joint and individual impacts of child metal exposures (log2-transformed) on lipid levels, Bayesian kernel machine regression (BKMR) was employed as the primary multi-pollutant approach. Potential effect modification by child sex and childhood environmental tobacco smoke exposure was also evaluated.

RESULTS

BKMR identified a positive association between the metal mixture and both total and LDL cholesterol. Of the seven metals examined, selenium (median 90.6 [IQR = 83.6, 96.5] µg/L) was assigned the highest posterior inclusion probability for both total and LDL cholesterol. A difference in LDL cholesterol of 8.22 mg/dL (95% CI = 1.85, 14.59) was observed when blood selenium was set to its 75th versus 25th percentile, holding all other metals at their median values. In stratified analyses, the positive association between selenium and LDL cholesterol was only observed among boys or among children exposed to environmental tobacco smoke during childhood.

IMPACT STATEMENT

Growing evidence indicates that cardiovascular events in adulthood are the consequence of the lifelong atherosclerotic process that begins in childhood. Therefore, public health interventions targeting childhood cardiovascular risk factors may have a particularly profound impact on reducing the burden of cardiovascular disease. Although growing evidence supports that both essential and nonessential metals contribute to cardiovascular disease and risk factors, such as dyslipidemia, prior studies have mainly focused on single metal exposures in adult populations. To address this research gap, the current study investigated the joint impacts of multiple metal exposures on lipid concentrations in early childhood.

Modeling Historic Arsenic Exposures and Spatial Risk for Bladder Cancer

Arsenic is a bladder carcinogen though less is known regarding the specific temporal relationship between exposure and bladder cancer diagnosis. In this study, we modeled time-varying mixtures of arsenic exposures at many historic temporal windows to evaluate their association with bladder cancer risk in the New England Bladder Cancer Study. We used arsenic exposure estimates up to 60 years prior to study entry and compared the goodness of fit of models using these mixtures to those using summary measures of arsenic exposures. We used the Bayesian index low rank kriging multiple membership model (LRK-MMM) to estimate the associations of these mixtures with bladder cancer and estimate cumulative spatial risk for bladder cancer using participants' residential histories. We found consistent evidence that modeling arsenic exposures as a time-varying mixture provided better fit to the data than using a single arsenic exposure summary measure. We estimated several positive though not significant associations of the time-varying arsenic mixtures with bladder cancer having odds ratios (ORs) of 1.03-1.14 and identified many significant and positive associations for an interaction among those who consumed water from a private dug well (ORs 1.28-1.60). Arsenic exposures 40-50 years before study entry received elevated importance weights in these mixtures. Additionally, we found two small areas of elevated cumulative spatial risk for bladder cancer in southern New Hampshire and in south central Maine. These results emphasize the importance of considering time-varying mixtures of exposures for diseases with long latencies such as bladder cancer.

The effect of bisphenols on sex and thyroid hormone concentrations in cord blood among newborns

This study aims to investigate the associations of bisphenols with sex and thyroid hormones in cord blood among newborns. Four bisphenols, three hormones related to gonadal function, and four parameters related to thyroid function were measured in umbilical cord blood in 378 mother-newborn pairs. Multivariable linear regression, quantile-based g-computation (QGC), and Bayesian kernel machine regression were used. In the multivariable linear regression, bisphenol A (BPA) was associated with increased testosterone (TT) (regression coefficient, β = 0.049, 95% confidence interval, CI: 0.013,0.085; p = 0.007) and free tri-iodothyronine (FT3) levels (β = 0.019, 95% CI: 0.003, 0.035; p = 0.023), and decreased thyroid peroxidase antibody (TPOAb) (β = -0.053, 95% CI: 0.098, -0.008; p = 0.021). Consistently associations were observed in males, except TT, which was observed in females, and bisphenol AF (BPAF) was associated with decreased follicle-stimulating hormone (FSH) in females. These associations were also observed in a mixture of bisphenols. Moreover, we observed maternal prepregnancy body mass index (BMI) and delivery mode disparity in the relationship between bisphenols and sex and thyroid hormones. This study suggests that bisphenols may exert effects on sex and thyroid hormones in newborns, the effect may vary with sex differences, maternal prepregnancy BMI, and delivery mode.

Per- and polyfluoroalkyl substances and bone mineral content in early adolescence: Modification by diet and physical activity

BACKGROUND

Per- and polyfluoroalkyl substance (PFAS) exposures may negatively impact bone mineral accrual, but little is known about potential mitigators of this relation. We assessed whether associations of PFAS and their mixture with bone mineral content (BMC) in adolescence were modified by diet and physical activity.

METHODS

We included 197 adolescents enrolled in a prospective pregnancy and birth cohort in Cincinnati, Ohio (2003-2006). At age 12 years, we collected serum for PFAS measurements and used dual-energy x-ray absorptiometry to measure BMC. We calculated dietary calcium intake and Health Eating Index (HEI) scores from repeated 24-h dietary recalls, physical activity scores using the Physical Activity Questionnaire for Older Children (PAQ-C), and average moderate to vigorous physical activity (MVPA) based on accelerometry. We estimated covariate-adjusted differences in BMC z-scores per interquartile range (IQR) increase of individual PFAS concentrations using linear regression and per simultaneous IQR increase in all four PFAS using g-computation. We evaluated effect measure modification (EMM) using interaction terms between each modifier and PFAS.

RESULTS

Higher serum perfluorooctanoic acid, perfluorooctanesulfonic acid, and perfluorononanoic acid concentrations and the PFAS mixture were associated with lower BMC z-scores. An IQR increase in all PFAS was associated with a 0.27 (-0.54, 0.01) lower distal radius BMC z-score. Associations with lower BMC were generally stronger among adolescents classified as < median for calcium intake, HEI scores, or MVPA compared to those ≥ median. The difference in distal radius BMC z-score per IQR increase in all PFAS was -0.38 (-0.72, -0.04) for those with

CONCLUSION

Healthy, calcium-rich diets and higher intensity physical activity may mitigate the adverse impact of PFAS on adolescent bone health.

Collapse

Key Words

• Adolescent health
• Bone mineral content
• Calcium
• Diet
• Perfluoroalkyl substances
• Physical activity

Collapse

MESH Headings

• Humans
• Female
• Fluorocarbons/blood
• Male
• Bone Density/drug effects
• Child
• Adolescent
• Diet
• Environmental Pollutants/blood
• Prospective Studies
• Ohio
• Alkanesulfonic Acids/blood
• Exercise
• Motor Activity/drug effects

Collapse

Grants

• R01 ES030078 NIEHS NIH HHS
• R01 ES033252 NIEHS NIH HHS
• R01 ES027224 NIEHS NIH HHS
• L40 ES032257 NIEHS NIH HHS
• P01 ES011261 NIEHS NIH HHS
• R01 ES032836 NIEHS NIH HHS
• R01 ES028277 NIEHS NIH HHS
• T32 ES007141 NIEHS NIH HHS
• R01 ES025214 NIEHS NIH HHS
• R01 ES020349 NIEHS NIH HHS

Collapse

Affiliation(s)

Jessie P Buckley Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA.
Junyi Zhou Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, USA
Katherine M Marquess Department of Environmental Health and Engineering, Johns Hopkins Bloomberg School of Public Health, Baltimore, USA
Bruce P Lanphear Faculty of Health Sciences, Simon Fraser University, Vancouver, Canada
Kim M Cecil Department of Radiology, Cincinnati Children's Hospital Medical Center, University of Cincinnati College of Medicine, Cincinnati, USA; Department of Pediatrics, Cincinnati Children's Hospital Medical Center, University of Cincinnati College of Medicine, Cincinnati, USA
Aimin Chen Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania Perelman School of Medicine, Philadelphia, USA
Clara G Sears Christina Lee Brown Envirome Institute, Department of Medicine, Division of Environmental Medicine, University of Louisville, KY, USA
Yingying Xu Department of Pediatrics, Cincinnati Children's Hospital Medical Center, University of Cincinnati College of Medicine, Cincinnati, USA
Kimberly Yolton Department of Pediatrics, Cincinnati Children's Hospital Medical Center, University of Cincinnati College of Medicine, Cincinnati, USA
Heidi J Kalkwarf Department of Pediatrics, Cincinnati Children's Hospital Medical Center, University of Cincinnati College of Medicine, Cincinnati, USA
Joseph M Braun Department of Epidemiology, Brown University, Providence, USA
Jordan R Kuiper Department of Environmental and Occupational Health, The George Washington University Milken Institute School of Public Health, Washington, D.C., USA

Collapse

Epigenetic footprints: Investigating placental DNA methylation in the context of prenatal exposure to phenols and phthalates

BACKGROUND

Endocrine disrupting compounds (EDCs) such as phthalates and phenols can affect placental functioning and fetal health, potentially via epigenetic modifications. We investigated the associations between pregnancy exposure to synthetic phenols and phthalates estimated from repeated urine sampling and genome wide placental DNA methylation.

METHODS

The study is based on 387 women with placental DNA methylation assessed with Infinium MethylationEPIC arrays and with 7 phenols, 13 phthalates, and two non-phthalate plasticizer metabolites measured in pools of urine samples collected twice during pregnancy. We conducted an exploratory analysis on individual CpGs (EWAS) and differentially methylated regions (DMRs) as well as a candidate analysis focusing on 20 previously identified CpGs. Sex-stratified analyses were also performed.

RESULTS

In the exploratory analysis, when both sexes were studied together no association was observed in the EWAS. In the sex-stratified analysis, 114 individual CpGs (68 in males, 46 in females) were differentially methylated, encompassing 74 genes (36 for males and 38 for females). We additionally identified 28 DMRs in the entire cohort, 40 for females and 42 for males. Associations were mostly positive (for DMRs: 93% positive associations in the entire cohort, 60% in the sex-stratified analysis), with the exception of several associations for bisphenols and DINCH metabolites that were negative. Biomarkers associated with most DMRs were parabens, DEHP, and DiNP metabolite concentrations. Some DMRs encompassed imprinted genes including APC (associated with parabens and DiNP metabolites), GNAS (bisphenols), ZIM2;PEG3;MIMT1 (parabens, monoethyl phthalate), and SGCE;PEG10 (parabens, DINCH metabolites). Terms related to adiposity, lipid and glucose metabolism, and cardiovascular function were among the enriched phenotypes associated with differentially methylated CpGs. The candidate analysis identified one CpG mapping to imprinted LGALS8 gene, negatively associated with ethylparaben.

CONCLUSIONS

By combining improved exposure assessment and extensive placental epigenome coverage, we identified several novel genes associated with the exposure, possibly in a sex-specific manner.

Social Vulnerability and Biological Aging in New York City: An Electronic Health Records-Based Study

Chronological age is not an accurate predictor of morbidity and mortality risk, as individuals' aging processes are diverse. Phenotypic age acceleration (PhenoAgeAccel) is a validated biological age measure incorporating chronological age and biomarkers from blood samples commonly used in clinical practice that can better reflect aging-related morbidity and mortality risk. The heterogeneity of age-related decline is not random, as environmental exposures can promote or impede healthy aging. Social Vulnerability Index (SVI) is a composite index accounting for different facets of the social, economic, and demographic environment grouped into four themes: socioeconomic status, household composition and disability, minority status and language, and housing and transportation. We aim to assess the concurrent and combined associations of the four SVI themes on PhenoAgeAccel and the differential effects on disadvantaged groups. We use electronic health records data from 31,913 patients from the Mount Sinai Health System (116,952 person-years) and calculate PhenoAge for years with available laboratory results (2011-2022). PhenoAge is calculated as a weighted linear combination of lab results and PhenoAgeAccel is the differential between PhenoAge and chronological age. A decile increase in the mixture of SVI dimensions was associated with an increase of 0.23 years (95% CI: 0.21, 0.25) in PhenoAgeAccel. The socioeconomic status dimension was the main driver of the association, accounting for 61% of the weight. Interaction models revealed a more substantial detrimental association for women and racial and ethnic minorities with differences in leading SVI themes. These findings suggest that neighborhood-level social vulnerability increases the biological age of its residents, increasing morbidity and mortality risks. Socioeconomic status has the larger detrimental role amongst the different facets of social environment.

New findings on the risk of hypertension from organophosphorus exposure under different glycemic statuses: The key role of lipids?

BACKGROUND AND OBJECTIVE

Considering the widespread use of organophosphorus pesticides (OPs) and the global prevalence of hypertension (HTN), as well as studies indicating that different glycemic statuses may respond differently to the biological effects of OPs. Therefore, this study, based on the Henan rural cohort, aims to investigate the association between OPs exposure and HTN, and further explores whether lipids mediate these associations.

METHODS

We measured the plasma levels of OPs in 2730 participants under different glycemic statuses using gas chromatography-triple quadrupole mass spectrometry (GC-MS/MS). A generalized linear model, Quantile g-computation (QGC), adaptive elastic net (AENET), and Bayesian kernel machine regression (BKMR) models were used to assess the impact of OPs exposure on HTN, with least absolute shrinkage and selection operator (LASSO) penalty regression identifying main OPs. Mediation models were used to evaluate the intermediary role of blood lipids in the OPs-HTN relationship.

RESULTS

The detection rates for all OPs were high, ranging from 76.35 % to 99.17 %. In the normal glucose tolerance (NGT) population, single exposure models indicated that malathion and phenthoate were associated with an increased incidence of HTN (P-FDR < 0.05), with corresponding odds ratios (ORs) and 95 % confidence intervals (CIs) of 1.624 (1.167,2.260) and 1.290 (1.072,1.553), respectively. QGC demonstrated a positive association between OP mixtures and HTN, with malathion and phenthoate being the primary contributors. Additionally, the AENET model's Exposure Response Score (ERS) suggested that the risk of HTN increases with higher ERS (P < 0.001). Furthermore, BKMR revealed that co-exposure to OPs increases HTN risk, with phenthoate having a significant impact. Furthermore, triglycerides (TG) mediated 6.55 % of the association between phenthoate and HTN. However, no association was observed in the impaired fasting glucose (IFG) and type 2 diabetes mellitus (T2DM) populations.

CONCLUSIONS

Our findings suggest that in the NGT population, OPs may significantly contribute to the development of HTN, proposing TG as a potential novel target for HTN prevention.

Associations of prenatal exposure to bisphenols with infant anthropometry: A prospective cohort study

BACKGROUND

Bisphenols (BPs) have been shown to exhibit developmental toxicities. Epidemiological evidence on prenatal BPs exposure and infant growth primarily confined scopes to specific BPs and birth outcomes, with few studies focusing on infant growth and reporting inconsistent findings. The joint effect of prenatal exposure to BPs mixture on infant growth was rarely studied.

OBJECTIVE

This study examined associations of prenatal exposure to individual bisphenol A (BPA) and its analogues (bisphenol F [BPF], bisphenol S [BPS], bisphenol AF [BPAF], and tetrachlorobisphenol A [TCBPA]) and their mixture with infant growth.

METHODS

Urinary concentrations of BPs in pregnant women were quantified. Weight, body mass index, skinfold thickness, and circumference measurements of infants were collected at birth, 6 and 12 months of age, rapid growth and overweight were further defined. Multiple linear regression models and Bayesian kernel machine regression models (BKMR) were used to analyze associations of exposure to individual BPs and BPs mixture with infants' anthropometric measurements, and to identify the important components among mixture. The risks for rapid growth and overweight of each BP were determined using modified Poisson regression models.

RESULTS

A general profile of higher prenatal BPs exposure (mainly BPA, BPF, and BPS) associated with higher anthropometric measurements and higher risks of overweight during infancy was found. We also observed higher risks of rapid growth in infants following prenatal BPs exposure, with risk ratios ranging from 1.46 to 1.91. The joint effect of BPs mixture and single effect of each BP from the BKMR models were consistent with findings from the linear regression models, further suggesting that associations in girls were generally driven by BPA, BPF, or BPS, while in boys mainly by BPF.

CONCLUSION

Prenatal exposure to BPs and their mixture could increase anthropometric measurements of offspring during infancy, with implications of altered growth trajectory in future.

Associations of toxic metals and their mixture with hyperuricemia in Chinese rural older adults

Previous studies have related single toxic metals (TMs) to hyperuricemia (HUA) among the general population, however, the association of the TM mixture with HUA, especially in older adults, remains poorly understood. We aimed to examine the relationships between individual TMs and their mixture and HUA in Chinese rural older adults. This study consisted of 2075 rural older adults aged 60 years or over. Blood concentrations of aluminum (Al), arsenic (As), barium (Ba), cadmium (Cd), cesium (Cs), gallium (Ga), mercury (Hg), lead (Pb), thallium (Tl), and uranium (U) were detected using inductively coupled plasma mass spectrometry. The associations of single TMs with HUA were assessed using logistic regression and restricted cubic spline (RCS) models, and the association of TM mixture with HUA was explored using the elastic net with environmental risk score (ENET-ERS), quantile g-computation (QGC), and Bayesian kernel machine regression (BKMR) models, respectively. Adjusted logistic regression model showed that Cs (OR = 1.65, 95% CI 1.37-1.99) and Pb (OR = 1.46, 95% CI 1.28-1.67) were positively related to HUA, and RCS model exhibited a positive linear association of Cs and Pb with HUA. ENET-ERS and QGC models quantified a positive correlation between the TM mixture and the odds of HUA, with estimated ORs of 1.15 (95% CI 1.11-1.19) and 1.84 (95% CI 1.37-2.47), respectively, and Cs and Pb had the most weight. BKMR model demonstrated a significant linear association between the TM mixture and increased odds of HUA, with the posterior inclusion probabilities (PIPs) of both Cs and Pb being 1.00. Moreover, we observed a positive interaction between Cs and Pb on HUA. The TM mixture is associated with increased odds of HUA in rural older adults, which may mainly be driven by Cs and Pb. Subsequent studies are warranted to confirm these findings and clarify the mechanisms linking multiple TMs with HUA.

Joint B Vitamin Intake and Type 2 Diabetes Risk: The Mediating Role of Inflammation in a Prospective Shanghai Cohort

BACKGROUND AND AIMS

Type 2 diabetes (T2D) is a global and complex public health challenge, and dietary management is acknowledged as critical in its prevention. Recent studies have highlighted the involvement of micronutrients in T2D pathophysiology; our study aims to assess the association between B vitamin intake and T2D risks and the mediating role of inflammation.

METHODS

In a prospective cohort design, data on B vitamins intake, including thiamine (B1), riboflavin (B2), niacin (B3), pyridoxine (B6), folate (B9), and cobalamin (B12), was obtained using a validated food frequency questionnaire (FFQ), and blood inflammatory biomarkers were analyzed according to standard protocol in the local hospitals at baseline from 44,960 adults in the Shanghai Suburban Adult Cohort and Biobank (SSACB). Incident T2D cases were identified according to a physician's diagnosis or medication records from the electronic medical information system. We employed logistic and weighted quantile sum regression models to explore the associations of single and combined levels of B vitamins with T2D and mediation analyses to investigate the effects of inflammation.

RESULTS

Negative correlations between B vitamins and T2D were observed in the single-exposure models, except for B3. The analyses of joint exposure (B1, B2, B6, B9, and B12) also showed an inverse association (OR 0.80, 95% CI 0.71 to 0.88), with vitamin B6 accounting for 45.58% of the effects. Further mediation analysis indicated a mediating inflammatory impact, accounting for 6.72% of the relationship.

CONCLUSIONS

Dietary intake of B vitamins (B1, B2, B6, B9, B12) was associated with a reduced T2D risk partially mediated by inflammation in Shanghai residents.

Complex interplay of heavy metals and renal injury: New perspectives from longitudinal epidemiological evidence

BACKGROUND

Epidemiological studies have reported associations between heavy metals and renal function. However, longitudinal studies are required to further validate these associations and explore the interactive effects of heavy metals on renal function and their directional influence.

METHOD

This study, conducted in Northeast China from 2016 to 2021, included a four-time repeated measures design involving 384 participants (1536 observations). Urinary concentrations of chromium (Cr), cadmium (Cd), manganese (Mn), and lead (Pb) were measured, along with renal biomarkers including urinary microalbumin (umAlb), urinary albumin-to-creatinine ratio (UACR), N-acetyl-β-D-glucosaminidase (NAG), and β2-microglobulin (β2-MG) levels. Estimated glomerular filtration rate (eGFR) was calculated. A Linear Mixed Effects Model (LME) examined the association between individual metal exposure and renal biomarkers. Subsequently, Quantile g-computation and Bayesian Kernel Machine Regression (BKMR) models assessed the overall effects of heavy metal mixtures. Marginal Effect models examined the directional impact of metal interactions in the BKMR on renal function.

RESULT

Results indicate significant impacts of individual and combined exposures of Cr, Cd, Pb, and Mn on renal biomarkers. Metal interactions in the BKMR model were observed, with synergistic effects of Cd-Cr on NAG, umAlb, UACR; Cd-Pb on NAG, UACR; Pb-Cr on umAlb, UACR, eGFR-MDRD, eGFR-EPI; and an antagonistic effect of Mn-Pb-Cr on UACR.

CONCLUSION

Both individual and combined exposures to heavy metals are associated with renal biomarkers, with significant synergistic interactions leading to renal damage. Our findings elucidate potential interactions among these metals, offering valuable insights into the mechanisms linking multiple metal exposures to renal injury.

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

BACKGROUND AND AIMS

Exposure to brominated flame retardants (BFRs) has been widely confirmed to impair the normal functioning of the human body system. However, there is a paucity of study on the effects of serum BFRs on bone mineral density (BMD). This study aims to investigate the relationship between exposure to BFRs and BMD in a nationally representative sample of U.S. adults.

METHODS

3079 participants aged between 20 and 80 years with complete data were included in the study. Serum levels of BFRs were measured using automated liquid-liquid extraction and subsequent sample clean-up. The BMD of all participants were assessed by DXA examinations. Generalize linear model, Restricted cubic spline (RCS), subgroup, weighted quantile sum (WQS) and bayesian kernel machine regression (BKMR) were used to estimate the association between serum BFRs and BMD.

RESULTS

Multivariate linear regression analyses revealed that, after adjusting for covariates, PBB153 was significantly associated with TF-BMD (β = 0.0177, 95%CI: 0.0103-0.0252), FN-BMD (β = 0.009, 95%CI: 0.0036-0.0145), TS-BMD (β = 0.0081, 95%CI: 0.0013-0.015) and L1-BMD (β = 0.0144, 95%CI: 0.0075-0.0213). However, the associations lose their statistical significance after further adjustment for sex. BFRs exhibited S-shaped or line-plateau dose-response curves with BMD. In subgroup analyses, BFRs were significantly associated with BMD in participants who were younger than 55 years, female, overweight (BMI >25 kg/m2), and less alcohol consumption. In WQS and BKMR analyses, the effects of BFRs mixtures on BMD differed by sex, and PBDE153, PBDE209 and PBB153 had the highest weights in the WQS regression model.

CONCLUSION

This study showed that serum BFRs negatively predicted BMD in men, but not in women or the general population. PBDE153, PBDE209, and PBB153 were significant BMD factors, especially in younger, overweight, and less alcohol consumption individuals.

Associations between organophosphate esters and bone mineral density in adults in the United States: 2011-2018 NHANES

BACKGROUND

Organophosphate esters (OPEs) are used extensively as flame retardants and plasticizers. Laboratory studies have shown that OPEs exhibit osteotoxicity by inhibiting osteoblast differentiation; however, little is known about how OPEs exposure is associated with bone health in humans.

OBJECTIVES

We conducted a cross-sectional study to investigate the association between OPEs exposure and bone mineral density (BMD) in adults in the United States using data from the 2011-2018 National Health and Nutrition Examination Survey (NHANES).

METHODS

Multivariate linear regression models were used to assess the association between concentrations of individual OPE metabolites and BMDs. We also used the Bayesian kernel machine regression (BKMR) and quantile g-computation (qgcomp) models to estimate joint associations between OPE mixture exposure and BMDs. All the analyses were stratified according to gender.

RESULTS

A total of 3546 participants (median age, 40 years [IQR, 30-50 years]; 50.11% male) were included in this study. Five urinary OPE metabolites with a detection rate of > 50% were analyzed. After adjusting for the potential confounders, OPE metabolite concentrations were associated with decreased total-body BMD and lumbar spine BMD in males, although some associations only reached significance for bis(1-chloro-2-propyl) phosphate (BCPP), dibutyl phosphate (DBUP), and bis(2-chloroethyl) phosphate (BCEP) (β = -0.013, 95% CI: -0.026, -0.001 for BCPP and total-body BMD; β = -0.022, 95% CI: -0.043, -0.0001 for DBUP and lumbar spine BMD; β=-0.018, 95% CI: -0.034, -0.002 for BCEP and lumbar spine BMD). OPE mixture exposure was also inversely associated with BMD in males, as demonstrated in the BMKR and qgcomp models.

CONCLUSIONS

This study provides preliminary evidence that urinary OPE metabolite concentrations are inversely associated with BMD. The results also suggested that males were more vulnerable than females. However, further studies are required to confirm these findings.

Associations of phthalates with accelerated aging and the mitigating role of physical activity

Phthalates are positioned as potential risk factors for health-related diseases. However, the effects of exposure to phthalates on accelerated aging and the potential modifications of physical activity remain unclear. A total of 2317 participants containing complete study-related information from the National Health and Nutrition Examination Survey 2007-2010 were included in the current study. We used two indicators, the Klemera-Doubal method biological age acceleration (BioAgeAccel) and phenotypic age acceleration (PhenoAgeAccel), to assess the accelerated aging status of the subjects. Multiple linear regression (single pollutant models), weighted quantile sum (WQS) regression, Quantile g-computation, and Bayesian kernel machine regression (BKMR) models were utilized to explore the associations between urinary phthalate metabolites and accelerated aging. Three groups of physical activity with different intensities were used to evaluate the modifying effects on the above associations. Results indicated that most phthalate metabolites were significantly associated with BioAgeAccel and PhenoAgeAccel, with effect values (β) ranging from 0.16 to 0.21 and 0.16-0.37, respectively. The WQS indices were positively associated with BioAgeAccel (0.33, 95% CI: 0.11, 0.54) and PhenoAgeAccel (0.50, 95% CI: 0.19, 0.82). Quantile g-computation indicated that phthalate mixtures were associated with accelerated aging, with effect values of 0.15 (95% CI: 0.02, 0.28) for BioAgeAccel and 0.39 (95% CI: 0.12, 0.67) for PhenoAgeAccel respectively. The BKMR models indicated a significant positive association between the concentrations of urinary phthalate mixtures with the two indicators. In addition, we found that most phthalate metabolites showed the strongest effects on accelerated aging in the no physical activity group and that the effects decreased gradually with increasing levels of physical activity (P < 0.05 for trend). Similar results were also observed in the mixed exposure models (WQS and Quantile g-computation). This study indicates that phthalates exposure is associated with accelerated aging, while physical activity may be a crucial barrier against phthalates exposure-related aging.

Trimester-specific effect of maternal co-exposure to organophosphate esters and phthalates on preschooler cognitive development: The moderating role of gestational vitamin D status

Organophosphate esters (OPEs) and phthalate acid esters (PAEs) are prevalent endocrine-disrupting chemicals (EDCs). Humans are often exposed to OPEs and PAEs simultaneously through multiple routes. Given that fetal stage is a critical period for neurodevelopment, it is necessary to know whether gestational co-exposure to OPEs and PAEs affects fetal neurodevelopment. However, accessible epidemiological studies are limited. The present study included 2, 120 pregnant women from the Ma'anshan Birth Cohort (MABC) study. The concentrations of tris (2-chloroethyl) phosphate (TCEP), 6 OPE metabolites and 7 PAE metabolites were measured in the first, second and third trimester using ultra-performance liquid chromatography-tandem mass spectrometry (LC-MS). Cognitive development of preschooler was assessed based on the Wechsler Preschool and Primary Scale of Intelligence-Fourth Edition (WPPSI-IV) of the Chinese version. Generalized estimating equations (GEEs), restricted cubic spline (RCS) and generalized additive models (GAMs) were employed to explore the associations between individual OPE exposure and preschooler cognitive development. The quantile-based g-computation (QGC) method was used to estimate the joint effect of PAEs and OPEs exposure on cognitive development. GEEs revealed significant adverse associations between diphenyl phosphate (DPHP) (β: -0.58, 95% CI: -1.14, -0.01), bis (2-butoxyethyl) phosphate(BBOEP) (β: -0.44, 95% CI: -0.85, -0.02), bis(1-chloro-2-propyl) phosphate (BCIPP) (β: -0.81, 95%CI: -1.43, -0.20) and full-scale intelligence quotient (FSIQ) in the first trimester; additionally, TCEP and bis(2-ethylhexyl) phosphate (BEHP) in the second trimester, as well as DPHP in the third trimester, were negatively associated with cognitive development. Through the QGC analyses, mixture exposure in the first trimester was negatively associated with FSIQ scores (β: -1.70, 95% CI: -3.06, -0.34), mono-butyl phthalate (MBP), BCIPP, and DPHP might be the dominant contributors after controlling for other OPEs and PAEs congeners. Additionally, the effect of OPEs and PAEs mixture on cognitive development might be driven by vitamin D deficiency.

Multiple xenoestrogen air pollutants and breast cancer risk: Statistical approaches to investigate combined exposures effect

Studies suggested that exposure to air pollutants, with endocrine disrupting (ED) properties, have a key role in breast cancer (BC) development. Although the population is exposed simultaneously to a mixture of multiple pollutants and ED pollutants may act via common biological mechanisms leading to synergic effects, epidemiological studies generally evaluate the effect of each pollutant separately. We aimed to assess the complex effect of exposure to a mixture of four xenoestrogen air pollutants (benzo-[a]-pyrene (BaP), cadmium, dioxin (2,3,7,8-Tétrachlorodibenzo-p-dioxin TCDD)), and polychlorinated biphenyl 153 (PCB153)) on the risk of BC, using three recent statistical methods, namely weighted quantile sum (WQS), quantile g-computation (QGC) and Bayesian kernel machine regression (BKMR). The study was conducted on 5222 cases and 5222 matched controls nested within the French prospective E3N cohort initiated in 1990. Annual average exposure estimates to the pollutants were assessed using a chemistry transport model, at the participants' residence address between 1990 and 2011. We found a positive association between the WQS index of the joint effect and the risk of overall BC (adjusted odds ratio (OR) = 1.10, 95% confidence intervals (CI): 1.03-1.19). Similar results were found for QGC (OR = 1.11, 95%CI: 1.03-1.19). Despite the association did not reach statistical significance in the BKMR model, we observed an increasing trend between the joint effect of the four pollutants and the risk of BC, when fixing other chemicals at their median concentrations. BaP, cadmium and PCB153 also showed positive trends in the multi-pollutant mixture, while dioxin showed a modest inverse trend. Despite we found a clear evidence of a positive association between the joint exposure to pollutants and BC risk only from WQS and QGC regression, we observed a similar suggestive trend using BKMR. This study makes a major contribution to the understanding of the joint effects of air pollution.

Urinary phthalate/DINCH metabolites associations with kisspeptin and reproductive hormones in teenagers: A cross-sectional study from the HBM4EU aligned studies

BACKGROUND

Exposure to phthalate/DINCH metabolites can induce human reproductive toxicity, however, their endocrine-disrupting mechanisms are not fully elucidated.

OBJECTIVE

To investigate the association between concentrations of phthalate/DINCH metabolites, serum kisspeptin, and reproductive hormones among European teenagers from three of the HBM4EU Aligned Studies.

METHODS

In 733 Belgian (FLEHS IV study), Slovak (PCB cohort follow-up), and Spanish (BEA study) teenagers, ten phthalate and two DINCH metabolites were measured in urine by high-performance liquid chromatography-tandem mass spectrometry. Serum kisspeptin (kiss54) protein, follicle-stimulating hormone (FSH), total testosterone (TT), estradiol (E2), and sex hormone-binding globulin (SHBG) levels were measured by immunosorbent assays. Free Androgen Index (FAI) was calculated as a proxy of free testosterone. Adjusted sex-stratified linear regression models for individual studies, mixed effect models (LME) accounting for random effects for pooled studies, and g-computation and Bayesian kernel machine regression (BKMR) models for the phthalate/DINCH mixture were performed.

RESULTS

The LME suggested that each IQR increase in ln-transformed levels of several phthalates was associated with lower kisspeptin [MnBP: %change (95%CI): -2.8 (-4.2;-0.4); MEHP: -1.4 (-3.4,0.2)] and higher FSH [∑DINP: 11.8 (-0.6;25.1)] levels in females from pooled studies. G-computation showed that the phthalates/DINCH mixture was associated with lower kisspeptin [-4.28 (-8.07;-0.34)] and higher FSH [22.13 (0.5;48.4)] also in females; BKMR showed similar although non-significant pattern. In males, higher phthalates metabolites [MEHP: -12.22 (-21.09;-1.18); oxo-MEHP: -12.73 (-22.34;-1.93)] were associated with lower TT and FAI, although higher DINCH [OH-MINCH: 16.31 (6.23;27.35), cx-MINCH: 16.80 (7.03;27.46), ∑DINCH: 17.37 (7.26;29.74)] were associated with higher TT levels. No mixture associations were found in males.

CONCLUSION

We observed sex-specific associations between urinary concentrations of phthalate/DINCH metabolites and the panel of selected effect biomarkers (kisspeptin and reproductive hormones). This suggests that exposure to phthalates would be associated with changes in kisspeptin levels, which would affect the HPG axis and thus influence reproductive health. However, further research is needed, particularly for phthalate replacements such as DINCH.

Association between long-term exposure to fine particulate matter and its chemical constituents and premature death in individuals living with HIV/AIDS

Long-term exposure to fine particulate matter (PM2.5) is associated with an increased total mortality. However, the association of PM2.5 with mortality in people living with human immunodeficiency virus/acquired immunodeficiency syndrome (HIV/AIDS, PLWHA) and the relationship between its constituents and adverse outcomes remain unknown. In this cohort study, 28,140 PLWHA were recruited from the HIV/AIDS Comprehensive Response Information Management System of the Hubei Provincial Centre for Disease Control and Prevention in China between 2001 and 2020. The annual PM2.5 chemical composition data, including sulfate (SO42-), nitrate (NO3-), ammonium (NH4+), black carbon (BC), and organic matter (OM), was extracted from the Tracking Air Pollution (TAP) dataset in China. A Cox proportional hazard model with time-varying exposure and time-to-event quantile-based generalized (g) computation was used to assess the associations between PM2.5 chemical constituents, and mortality in PLWHA. A multivariate Cox proportional hazard model estimated an excess hazard ratio (eHR) of 0.32% [95% confidence interval (CI): (0.01%, 0.64%)] for AIDS-related death (ARD), associated with 1 μg/m3 rise in PM2.5 exposure. An increase of 1 μg/m3 in NH4+ was associated with 5.13% [95% CI: (2.89%, 7.43%)] and 2.97% [95% CI: (1.52%, 4.44%)] increase in the risk of ARD and all-cause deaths (ACD), respectively. When estimated using survival-based quantile g-computation, the eHR for ARD with a joint change in a decile increase in all five components was 6.10% [95% CI: 3.77%, 8.48%)]. Long-term exposure to PM2.5 chemical composition, particularly NH4+ increased the risk of death in PLWHA. This study provides epidemiological evidence that SO42- and NH4+ increased the risk of ARD and that NH4+ increased the risk of ACD in PLWHA. Multi-constituent analyses further suggested that NH4+ may be a key component in increasing the risk of premature death in patients with HIV/AIDS. Individuals aged ≥65 with HIV/AIDS are more vulnerable to SO42-, and consequent ACD.

Risk of dementia due to Co-exposure to air pollution and neighborhood disadvantage

BACKGROUND

Co-exposure to air pollution and neighborhood disadvantage may influence cognition decline. We tested these associations in the context of dementia risk.

METHODS

We leveraged a cohort of adults ≥65 years (n = 5397) enrolled from 2011 to 2018 in the National Health and Aging Trends Study (NHATS). Particulate matter (PM) ≤ 10 μm in diameter, PM ≤ 2.5 μm in diameter, carbon monoxide, nitric oxide, and nitrogen dioxide - and neighborhood disadvantage were tested for joint associations with dementia risk. Pollutant concentrations at the 2010 census tract level were assigned using the US Environmental Protection Agency's Community Multiscale Air Quality Modeling System. Neighborhood disadvantage was defined using the tract Social Deprivation Index (SDI). Dementia was determined through self- or proxy-report or scores indicative of "probable dementia" according to NHATS screening tools. Joint effects of air pollutants and SDI were tested using quantile g-computation Cox proportional hazards models. We also stratified joint air pollution effects across SDI tertiles. Analyses adjusted for age at enrollment, sex, education, partner status, urbanicity, income, race and ethnicity, years at residence, census segregation, and census region.

RESULTS

SDI score (aHR = 1.08; 95% CI 0.96, 1.22), joint air pollution (aHR = 1.03, 95% CI 0.92, 1.16) and joint SDI with air pollution (aHR = 1.04, 95% CI 0.89, 1.22) were not associated with dementia risk. After accounting for competing risk of death, joint SDI with air pollution was not associated with dementia risk (aHR = 1.06; 95% CI 0.87, 1.29). In stratified models, joint air pollution was associated with greater risk of dementia at high (aHR = 1.19; 95% CI 0.87, 1.63), but not at medium or low SDI.

CONCLUSION

Air pollution was associated with greater dementia risk in disadvantaged areas after accounting for competing risks. Air pollution associations with dementia incidence may be attenuated when other risk factors are more prominent in disadvantaged neighborhoods.

Environmental triggers of autoimmunity: The association between bisphenol analogues and systemic lupus erythematosus

The aim of this research was to examine the correlation between the exposure to bisphenol analogues (BPs), such as bisphenol A (BPA), bisphenol F (BPF), and bisphenol S (BPS), and the risk of developing systemic lupus erythematosus (SLE). Ultra performance liquid chromatography/tandem mass spectrometry (UPLC-MS/MS) was utilized to measure the levels of BPA, BPF, and BPS in the urine of 168 female participants diagnosed with SLE and 175 female participants who were deemed healthy controls. Logistic regression models were utilized to assess the connections between levels of bisphenol and the risk of SLE. The findings indicated that levels of BPA and BPF in the urine of individuals with SLE were markedly elevated compared to those in the control group. Higher exposure to BPA and BPF exhibited positive dose-response relationships with increased SLE risk. No significant associations were identified between BPS and the risk of SLE. These findings suggest exposure to BPA and BPF may be implicated as novel environmental triggers in the development of autoimmunity such as SLE. The significantly increased levels of these bisphenol analogues detected in SLE patients versus healthy controls, along with the associations between higher exposures and elevated SLE risk, which offers crucial hints for comprehending how endocrine-disrupting substances contribute to the genesis of autoimmune illnesses. Further research using robust longitudinal assessments of bisphenol analogue exposures is warranted to corroborate these epidemiological findings. Overall, this study highlights potential environmental risk factors for SLE while calling for additional investigation into the impact of bisphenol exposures on autoimmunity development.

Environmental exposures to organophosphorus flame retardants in early pregnancy and risks of gestational diabetes mellitus: a nested case-control study

OPFRs are emerging environmental pollutants with reproductive and endocrine toxicity. This study aimed to examine the association between environmental exposure to OPFRs during early pregnancy and GDM. This nested case-control study was based on a birth cohort that was constructed at a maternal and child health hospital, including 74 cases of GDM among 512 pregnant women. The OPFRs, including TBP, TBEP, TCEP, TDCPP, TMCP, TOCP, and TPHP during 10-14 weeks of pregnancy were determined using GC-MS. The association between the OPFRs and GDM was assessed using WQS and BKMR models. The levels of OPFRs were significantly elevated in GDM patients (60) compared with the controls (90). The WQS analysis showed that mixtures of the OPFRs were significantly associated with GDM (OR 1.370, 95% CI 1.036-1.810, P = 0.027), and TBP, TPHP, and TMCP were the major contributors to the mixed exposure effect. In the BKMR model, individual exposure to TBP, TPHP, and TMCP, and the interaction of TMCP with TBP and TPHP were significantly associated with GDM. Environmental exposure to OPFRs is positively associated with GDM. These findings provide evidence for the adverse effects of OPFR exposure on the health of pregnant women.

Sexual dimorphism association of combined exposure to volatile organic compounds (VOC) with kidney damage

BACKGROUND

Epidemiological evidence emphasizes air pollutants' role in chronic kidney disease (CKD). Volatile organic compounds (VOCs) contribute to air pollution, yet research on VOCs and kidney damage, especially gender disparities, is limited.

METHODS

This study analyzed NHANES data to explore associations between urinary VOC metabolite mixtures (VOCMs) and key kidney-related parameters: estimated glomerular filtration rate (eGFR), albumin-to-creatinine ratio (ACR), chronic kidney disease (CKD), and albuminuria. Mediation analyses assessed the potential mediating roles of biological aging (BA) and serum albumin in VOCM mixtures' effects on kidney damage. Sensitivity analyses were also conducted.

RESULTS

The mixture analysis unveiled a noteworthy positive association between VOCM mixtures and the risk of developing CKD, coupled with a significant negative correlation with eGFR within the overall participant cohort. These findings remained consistent when examining the female subgroup. However, among male participants, no significant link emerged between VOCM mixtures and CKD or eGFR. Furthermore, in both the overall and female participant groups, there was an absence of a significant correlation between VOCM mixtures and either ACR or albuminuria. On the other hand, in male participants, while no significant correlation was detected with albuminuria, a significant positive correlation was observed with ACR. Pollutant analysis identified potential links between kidney damage and 1,3-butadiene, toluene, ethylbenzene, styrene, xylene, acrolein, crotonaldehyde and propylene oxide. Mediation analyses suggested that BA might partially mediate the relationship between VOCM mixtures and kidney damage.

CONCLUSION

The current findings highlight the widespread exposure to VOCs among the general U.S. adult population and indicate a potential correlation between exposure to VOC mixtures and compromised renal function parameters, with notable gender disparities. Females appear to exhibit greater sensitivity to impaired renal function resulting from VOCs exposure. Anti-aging treatments may offer some mitigation against kidney damage due to VOCs exposure.

Exposure to iodine, essential and non-essential trace element through seaweed consumption in humans

Seaweed consumption has gained popularity due to its nutritional value and potential health benefits. However, concerns regarding the bioaccumulation of several trace elements highlight the need for comprehensive studies on exposure associated with seaweed consumption. To address this gap in knowledge, we carried out a feeding intervention study of three common edible seaweeds (Nori, Kombu, and Wakame) in 11 volunteers, aiming to elucidate the extent of both beneficial and harmful trace element exposure through seaweed consumption in humans. Concentrations of total arsenic, cobalt, copper, cadmium, iodine, molybdenum, selenium, and zinc were measured in urine samples before and following seaweed consumption. Elements concentrations were also measured in the seaweeds provided for the study. Descriptive analysis for each element were conducted and we used quantile g-computation approach to assess the association between the 8-element mixture and seaweed consumption. Differences in urine element concentrations and seaweed consumption were analyzed using generalized estimating equations (GEE). Urinary concentrations of iodine and total arsenic increased after seaweed consumption. When we analyze the 8-element mixture, the largest weight was observed for iodine after Kombu consumption while for total arsenic was observed after Wakame consumption. Similar results were observed when we compared the mean differences between the elements before and after seaweed consumption through the GEE. Seaweed consumption relates with increased urinary iodine and total arsenic concentrations, particularly after Kombu and Wakame consumption.

Associations of Urinary Heavy Metal Mixtures with High Remnant Cholesterol among US Adults: Evidence from the National Health and Nutrition Examination Survey (1998-2018)

The main objective of our study is to explore the associations between combined exposure to urinary heavy metals and high remnant cholesterol (HRC), a known cardiovascular risk factor. Utilizing data from the National Health and Nutrition Examination Survey (NHANES) from 1999 to 2018, we conducted a cross-sectional analysis of 5690 participants, assessing urinary concentrations of ten heavy metals. Ten heavy metals in urine were measured by inductively coupled plasma mass spectrometry (ICP-MS). Fasting residual cholesterol ≥0.8 mmol/L was defined as HRC (using blood samples). Statistical analyses included weighted multivariable logistic regression, weighted quantile sum (WQS) regression, quantile g-computation (qgcomp), and Bayesian kernel machine regression (BKMR) to evaluate the associations of heavy metal exposure with HRC. Stratified analyses based on individual characteristics were also conducted. Multivariable logistic regression found that the four metals (OR Q4 vs. Q1: 1.33, 95% CI: 1.01-1.75 for barium (Ba); OR Q4 vs. Q1: 1.50, 95% CI: 1.16-1.94 for cadmium (Cd); OR Q4 vs. Q1: 1.52, 95% CI: 1.15-2.01 for mercury (Hg); OR Q4 vs. Q1: 1.35, 95% CI: 1.06-1.73 for lead (Pb)) were positively correlated with the elevated risk of HRC after adjusting for covariates. In addition, all three mixed models, including WQS (OR: 1.25; 95% CI: 1.07-1.46), qgcomp (OR: 1.17; 95% CI: 1.03-1.34), and BKMR, consistently showed a significant positive correlation between co-exposure to heavy metal mixtures and HRC, with Ba and Cd being the main contributors within the mixture. These associations were more pronounced in younger adults (20 to 59 years), males, and those with a higher body mass index status (≥25 kg/m2). Our findings reveal a significant relationship between exposure to the mixture of heavy metals and HRC among US adults, with Ba and Cd being the major contributors to the mixture's overall effect. Public health efforts aimed at reducing heavy metal exposure can help prevent HRC and, in turn, cardiovascular disease.

Association between pesticide exposure and thyroid function: analysis of Chinese and NHANES databases

Background

Pesticides are widely used in agricultural activities. Although pesticide use is known to cause damage to the human body, its relationship with thyroid function remains unclear. Therefore, this study aimed to investigate the association between pesticide exposure and thyroid function.

Methods

The Chinese database used included 60 patients with pyrethroid poisoning and 60 participants who underwent health checkups between June 2022 and June 2023. The NHANES database included 1,315 adults enrolled from 2007 to 2012. The assessed pesticide and their metabolites included 2,4-dichlorophenoxyacetic acid (2,4-D), 4-fluoro-3-phenoxybenzoic acid (4F3PB), para-nitrophenol (PN), 3-phenoxybenzoic acid (3P), and trans-dichlorovinyl-dimethylcyclopropane carboxylic acid (TDDC). The evaluated indicators of thyroid function were measured by the blood from the included population. The relationship between pesticide exposure and thyroid function indexes was investigated using linear regression, Bayesian kernel machine regression (BKMR), restricted cubic spline (RCS), and weighted quantile sum (WQS) models.

Results

The Chinese data showed that pesticide exposure was negatively correlated with the thyroid function indicators FT4, TT4, TgAb, and TPOAb (all p < 0.05). The BKMR model analysis of the NHANES data showed that the metabolic mixture of multiple pesticides was negatively associated with FT4, TSH, and Tg, similar to the Chinese database findings. Additionally, linear regression analysis demonstrated positive correlations between 2,4-D and FT3 (p = 0.041) and 4F3PB and FT4 (p = 0.003), whereas negative associations were observed between 4F3PB and Tg (p = 0.001), 4F3PB and TgAb (p = 0.006), 3P and TgAB (p = 0.006), 3P and TPOAb (p = 0.03), PN and TSH (p = 0.003), PN and TT4 (p = 0.031), and TDDC and TPOAb (p < 0.001). RCS curves highlighted that most pesticide metabolites were negatively correlated with thyroid function indicators. Finally, WQS model analysis revealed significant differences in the weights of different pesticide metabolites on the thyroid function indexes.

Conclusion

There is a significant negative correlation between pesticide metabolites and thyroid function indicators, and the influence weights of different pesticide metabolites on thyroid function indicators are significantly different. More research is needed to further validate the association between different pesticide metabolites and thyroid disease.

Heterogeneous associations of multiplexed environmental factors and multidimensional aging metrics

Complicated associations between multiplexed environmental factors and aging are poorly understood. We manipulated aging using multidimensional metrics such as phenotypic age, brain age, and brain volumes in the UK Biobank. Weighted quantile sum regression was used to examine the relative individual contributions of multiplexed environmental factors to aging, and self-organizing maps (SOMs) were used to examine joint effects. Air pollution presented a relatively large contribution in most cases. We also found fair heterogeneities in which the same environmental factor contributed inconsistently to different aging metrics. Particulate matter contributed the most to variance in aging, while noise and green space showed considerable contribution to brain volumes. SOM identified five subpopulations with distinct environmental exposure patterns and the air pollution subpopulation had the worst aging status. This study reveals the heterogeneous associations of multiplexed environmental factors with multidimensional aging metrics and serves as a proof of concept when analyzing multifactors and multiple outcomes.

Exposure to volatile organic compounds and mortality in US adults: A population-based prospective cohort study

Volatile organic compounds (VOCs) are ubiquitous in both indoor and outdoor environments. Evidence on the associations of individual and joint VOC exposure with all-cause and cause-specific mortality is limited. Measurements of 15 urinary VOC metabolites were available to estimate exposure to 12 VOCs in the National Health and Nutritional Examination Survey (NHANES) 2005-2006 and 2011-2018. The environment risk score (ERS) was calculated using LASSO regression to reflect joint exposure to VOCs. Follow-up data on death were obtained from the NHANES Public-Use Linked Mortality File through December 31, 2019. Cox proportional hazard models and restricted cubic spline models were applied to evaluate the associations of individual and joint VOC exposures with all-cause and cause-specific mortality. Population attributable fractions were calculated to assess the death burden attributable to VOC exposure. During a median follow-up of 6.17 years, 734 (8.34 %) deaths occurred among 8799 adults. Urinary metabolites of acrolein, acrylonitrile, 1,3-butadiene, and ethylbenzene/styrene were significantly associated with all-cause, cardiovascular disease (CVD), respiratory disease (RD), and cancer mortality in a linear dose-response manner. Linear and robust dose-response relationships were also observed between ERS and all-cause and cause-specific mortality. Each 1-unit increase in ERS was associated with a 33.6 %, 39.1 %, 109.8 %, and 67.8 % increase for all-cause, CVD, RD, and cancer mortality risk, respectively. Moreover, joint exposure to VOCs contributed to 17.95 % of all-cause deaths, 13.49 % of CVD deaths, 35.65 % of RD deaths, and 33.85 % of cancer deaths. Individual and joint exposure to VOCs may enhance the risk of all-cause and cause-specific mortality. Reducing exposure to VOCs may alleviate the all-cause and cause-specific death burden.

Air pollution and survival in patients with malignant mesothelioma and asbestos-related lung cancer: a follow-up study of 1591 patients in South Korea

BACKGROUND

Despite significant advancements in treatments such as surgery, radiotherapy, and chemotherapy, the survival rate for patients with asbestos-related cancers remains low. Numerous studies have provided evidence suggesting that air pollution induces oxidative stress and inflammation, affecting acute respiratory diseases, lung cancer, and overall mortality. However, because of the high case fatality rate, there is limited knowledge regarding the effects of air pollution exposures on survival following a diagnosis of asbestos-related cancers. This study aimed to determine the effect of air pollution on the survival of patients with malignant mesothelioma and asbestos-related lung cancer.

METHODS

We followed up with 593 patients with malignant mesothelioma and 998 patients with lung cancer identified as asbestos victims between 2009 and 2022. Data on five air pollutants-sulfur dioxide, carbon monoxide, nitrogen dioxide, fine particulate matter with a diameter < 10 μm, and fine particulate matter with a diameter < 2.5 μm-were obtained from nationwide atmospheric monitoring stations. Cox proportional hazard models were used to estimate the association of cumulative air pollutant exposure with patient mortality, while adjusting for potential confounders. Quantile-based g-computation was used to assess the combined effect of the air pollutant mixture on mortality.

RESULTS

The 1-, 3-, and 5-year survival rates for both cancer types decreased with increasing exposure to all air pollutants. The estimated hazard ratios rose significantly with a 1-standard deviation increase in each pollutant exposure level. A quartile increase in the pollutant mixture was associated with a 1.99-fold increase in the risk of malignant mesothelioma-related mortality (95% confidence interval: 1.62, 2.44). For lung cancer, a quartile increase in the pollutant mixture triggered a 1.87-fold increase in the mortality risk (95% confidence interval: 1.53, 2.30).

CONCLUSION

These findings support the hypothesis that air pollution exposure after an asbestos-related cancer diagnosis can negatively affect patient survival.

Estimating effects of longitudinal and cumulative exposure to PFAS mixtures on early adolescent body composition

Few methods have been used to characterize repeatedly measured biomarkers of chemical mixtures. We applied latent profile analysis (LPA) to serum concentrations of 4 perfluoroalkyl and polyfluoroalkyl substances (PFAS) measured at 4 time points from gestation to age 12 years. We evaluated the relationships between profiles and z scores of height, body mass index, fat mass index, and lean body mass index at age 12 years (n = 218). We compared LPA findings with an alternative approach for cumulative PFAS mixtures using g-computation to estimate the effect of simultaneously increasing the area under the receiver operating characteristic curve (AUC) for all PFAS. We identified 2 profiles: a higher PFAS profile (35% of sample) and a lower PFAS profile (relative to each other), based on their average PFAS concentrations at all time points. The higher PFAS profile had generally lower z scores for all outcomes, with somewhat larger effects for males, though all 95% CIs crossed the null. For example, the higher PFAS profile was associated with a 0.50-unit lower (β = -0.50; 95% CI, -1.07 to 0.08) BMI z score among males but not among females (β = 0.04; 95% CI, -0.45 to 0.54). We observed similar patterns with AUCs. We found that a higher childhood PFAS profile and higher cumulative PFAS mixtures may be associated with altered growth in early adolescence. This article is part of a Special Collection on Environmental Epidemiology.

Relationships between minerals' intake and blood homocysteine levels based on three machine learning methods: a large cross-sectional study

BACKGROUND

Blood homocysteine (Hcy) level has become a sensitive indicator in predicting the development of cardiovascular disease. Studies have shown an association between individual mineral intake and blood Hcy levels. The effect of mixed minerals' intake on blood Hcy levels is unknown.

METHODS

Data were obtained from the baseline survey data of the Shanghai Suburban Adult Cohort and Biobank(SSACB) in 2016. A total of 38273 participants aged 20-74 years met our inclusion and exclusion criteria. Food frequency questionnaire (FFQ) was used to calculate the intake of 10 minerals (calcium, potassium, magnesium, sodium, iron, zinc, selenium, phosphorus, copper and manganese). Measuring the concentration of Hcy in the morning fasting blood sample. Traditional regression models were used to assess the relationship between individual minerals' intake and blood Hcy levels. Three machine learning models (WQS, Qg-comp, and BKMR) were used to the relationship between mixed minerals' intake and blood Hcy levels, distinguishing the individual effects of each mineral and determining their respective weights in the joint effect.

RESULTS

Traditional regression model showed that higher intake of calcium, phosphorus, potassium, magnesium, iron, zinc, copper, and manganese was associated with lower blood Hcy levels. Both Qg-comp and BKMR results consistently indicate that higher intake of mixed minerals is associated with lower blood Hcy levels. Calcium exhibits the highest weight in the joint effect in the WQS model. In Qg-comp, iron has the highest positive weight, while manganese has the highest negative weight. The BKMR results of the subsample after 10,000 iterations showed that except for sodium, all nine minerals had the high weights in the joint effect on the effect of blood Hcy levels.

CONCLUSION

Overall, higher mixed mineral's intake was associated with lower blood Hcy levels, and each mineral contributed differently to the joint effect. Future studies are available to further explore the mechanisms underlying this association, and the potential impact of mixed minerals' intake on other health indicators needs to be further investigated. These efforts will help provide additional insights to deepen our understanding of mixed minerals and their potential role in health maintenance.

Considerations and targeted approaches to identifying bad actors in exposure mixtures

Variable importance is a key statistical issue in exposure mixtures, as it allows a ranking of exposures as potential targets for intervention, and helps to identify bad actors within a mixture. In settings where mixtures have many constituents or high between-constituent correlations, estimators of importance can be subject to bias or high variance. Current approaches to assessing variable importance have major limitations, including reliance on overly strong or incorrect constraints or assumptions, excessive model extrapolation, or poor interpretability, especially regarding practical significance. We sought to overcome these limitations by applying an established doubly-robust, machine learning-based approach to estimating variable importance in a mixtures context. This method reduces model extrapolation, appropriately controls confounding, and provides both interpretability and model flexibility. We illustrate its use with an evaluation of the relationship between telomere length, a measure of biologic aging, and exposure to a mixture of polychlorinated biphenyls (PCBs), dioxins, and furans among 979 US adults from the National Health and Nutrition Examination Survey (NHANES). In contrast with standard approaches for mixtures, our approach selected PCB 180 and PCB 194 as important contributors to telomere length. We hypothesize that this difference could be due to residual confounding in standard methods that rely on variable selection. Further empirical evaluation of this method is needed, but it is a promising tool in the search for bad actors within a mixture.

Prenatal metal concentrations and physical abnormalities in the Japan Environment and Children's Study

BACKGROUND

The association between prenatal metal exposure and congenital anomalies is unclear. We aimed to examine the association between exposure to cadmium, lead, mercury, selenium, and manganese and physical abnormalities.

METHODS

Data from 89,887 pregnant women with singleton pregnancies who participated in the Japan Environment and Children's Study (JECS) were used. The correlation between maternal blood metal concentrations and physical abnormalities during the second or third trimester was investigated using logistic regression models. Physical anomalies included those observed at birth or at 1 month, primarily from ICD-10 Chapter 17, particularly congenital anomalies associated with environmental factors (e.g., hypospadias, cryptorchidism, cleft lip and palate, digestive tract atresia, congenital heart disease, and chromosomal abnormalities) and minor abnormalities.

RESULTS

After adjusting for covariates, the OR (95% CIs) of physical abnormalities for a one-unit rise in Mn concentrations in all individuals were 1.26 (1.08, 1.48). The OR (95% CIs) of physical abnormalities in the 4th quartile (≥18.7 ng/g) were 1.06 (1.01, 1.13) (p-value for the trend = 0.034) compared with those in the 1st quartile (≤12.5 ng/g).

CONCLUSION

In Japan, maternal blood Mn concentrations above threshold during pregnancy may slightly increase the incidence of physical abnormalities.

IMPACT

Physical abnormalities (including minor anomalies and congenital anomalies) are associated with prenatal manganese concentrations. They are not associated with cadmium, lead, mercury, and selenium concentrations.

Association between per- and poly-fluoroalkyl substances and nonalcoholic fatty liver disease: A nested case-control study in northwest China

Per- and poly-fluoroalkyl substances (PFAS) have been reported to have hepatotoxic effects. However, it is unclear whether they are linked to non-alcoholic fatty liver disease (NAFLD). This nested case-control study focused on the epidemiological links between PFAS and the prevalence of NAFLD. We selected 476 new cases of NAFLD and 952 age- and sex-matched controls from the Jinchang cohort population between 2014 and 2019. Serum concentrations of PFAS were measured using high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS). Only PFAS with a detection rate of ≥90 % were included for analysis, which included PFPeA, PFOA, PFNA, PFHxS, PFOS, and 9Cl-PF3ONS. The relationship between single and co-exposure to PFAS and the occurrence of NAFLD was evaluated using conditional logistic regression, Quantile g-computation (QgC), and Bayesian kernel machine regression (BKMR) model. Logistic regression indicated that PFPeA, PFOA, and 9Cl-PF3ONS were positive correlation with the incidence of NAFLD after adjusting for confounders, with odds ratios (OR) and 95 % confidence interval (CI) of 3.13 (95 % CI: 2.53, 3.86), 1.39 (95 % CI: 1.12, 1.73), and 1.41 (95 % CI: 1.20, 1.66), respectively. PFNA, PFHxS, and PFOS were nonlinearly and negatively associated with the incidence of NAFLD, with OR (95 % CI) of 0.53 (0.46, 0.62), 0.83 (0.73, 0.95), and 0.52 (0.44, 0.61), respectively. QgC showed a significant joint effect of PFAS mixture on NAFLD onset (OR: 1.52, 95 % CI: 1.24, 1.88). BKMR showed a weak positive trend between PFAS mixtures and NAFLD incidence. Positive correlations were primarily driven by PFPeA and 9Cl-PF3ONS, while negative correlations were mainly influenced by PFNA and PFOS. The BKMR model also suggested that there was an interaction between PFOS and PFNA and other four PFAS compounds. In conclusion, our findings suggest that individual and co-exposure to PFAS is associated with a risk of NAFLD onset.

Association between exposure to a mixture of organochlorine pesticides and hyperuricemia in U.S. adults: A comparison of four statistical models

The association between the exposure of organochlorine pesticides (OCPs) and serum uric acid (UA) levels remained uncertain. In this study, to investigate the combined effects of OCP mixtures on hyperuricemia, we analyzed serum OCPs and UA levels in adults from the National Health and Nutrition Examination Survey (2005-2016). Four statistical models including weighted logistic regression, weighted quantile sum (WQS), quantile g-computation (QGC), and bayesian kernel machine regression (BKMR) were used to assess the relationship between mixed chemical exposures and hyperuricemia. Subgroup analyses were conducted to explore potential modifiers. Among 6,529 participants, the prevalence of hyperuricemia was 21.15%. Logistic regression revealed a significant association between both hexachlorobenzene (HCB) and trans-nonachlor and hyperuricemia in the fifth quintile (OR: 1.54, 95% CI: 1.08-2.19; OR: 1.58, 95% CI: 1.05-2.39, respectively), utilizing the first quintile as a reference. WQS and QGC analyses showed significant overall effects of OCPs on hyperuricemia, with an OR of 1.25 (95% CI: 1.09-1.44) and 1.20 (95% CI: 1.06-1.37), respectively. BKMR indicated a positive trend between mixed OCPs and hyperuricemia, with HCB having the largest weight in all three mixture analyses. Subgroup analyses revealed that females, individuals aged 50 years and above, and those with a low income were more vulnerable to mixed OCP exposure. These results highlight the urgent need to protect vulnerable populations from OCPs and to properly evaluate the health effects of multiple exposures on hyperuricemia using mutual validation approaches.

Associations of per- and polyfluoroalkyl substances with maternal early second trimester sex-steroid hormones

BACKGROUND/AIMS

Pregnant women are exposed to persistent environmental contaminants, including per- and polyfluoroalkyl substances (PFAS) that disrupt thyroid function. However, it is unclear if PFAS alter maternal sex-steroid hormone levels, which support pregnancy health and fetal development.

METHODS

In Illinois women with relatively high socioeconomic status (n = 460), we quantified perfluorononanoic (PFNA), perfluorooctane sulfonic (PFOS), perfluorooctanoic (PFOA), methyl-perfluorooctane sulfonamide acetic acid, perfluorohexanesulphonic (PFHxS), perfluorodecanoic (PFDeA), and perfluoroundecanoic (PFUdA) acid concentrations in fasting serum samples at median 17 weeks gestation, along with plasma progesterone, testosterone, and estradiol. We evaluated covariate-adjusted associations of ln-transformed hormones with each ln-transformed PFAS individually using linear regression and with the PFAS mixture using quantile-based g-computation (QGComp).

RESULTS

Interquartile range (IQR) increases in PFOS were associated with higher progesterone (%Δ 3.0; 95%CI: -0.6, 6.6) and estradiol (%Δ: 8.1; 95%CI: 2.2, 14.4) levels. Additionally, PFHxS was positively associated with testosterone (%Δ: 10.2; 95%CI: 4.0, 16.7), whereas both PFDeA and PFUdA were inversely associated with testosterone (%Δ: -5.7; 95%CI: -10.3, -0.8, and %Δ: -4.1; 95%CI: -7.6, -0.4, respectively). The IQR-standardized PFAS mixture was not associated with progesterone (%Δ: 1.6; 95%CI: -5.8, 9.2), due equal partial positive (%Δ: 9.2; driven by PFOA) and negative (%Δ: -7.4; driven by PFOS) mixture associations. Similarly, the mixture was not associated with testosterone (%Δ: 5.3; 95%CI: -9.0, 20.1), due to similar partial positive (%Δ: 23.6; driven by PFHxS) and negative (%Δ: -17.4; driven by PFDeA) mixture associations. However, we observed a slightly stronger partial positive (%Δ: 25.6; driven by PFOS and PFUdA) than negative (%Δ: -16.3; driven by PFOA) association resulting in an overall non-significant positive trend between the mixture and estradiol (%Δ: 8.5; 95%CI: -3.7, 20.9).

CONCLUSION

PFAS mixture modeled using QGComp was not associated with maternal sex-steroid hormones due to potential opposing effects of certain PFAS. Additional prospective studies could corroborate these findings.

Individual and joint associations of per- and polyfluoroalkyl substances (PFAS) with gallstone disease in adults: A cross-sectional study

Disturbances in the enterohepatic circulation are important biological mechanisms for causing gallstones and also have important effects on the metabolism of Per- and polyfluoroalkyl substances (PFAS). Moreover, PFAS is associated with sex hormone disorder which is another important cause of gallstones. However, it remains unclear whether PFAS is associated with gallstones. In this study, we used logistic regression, restricted cubic spline (RCS), quantile g-computation (qg-comp), Bayesian kernel machine regression (BKMR), and subgroup analysis to assess the individual and joint associations of PFAS with gallstones and effect modifiers. We observed that the individual associations of perfluorodecanoic acid (PFDeA) (OR: 0.600, 95% CI: 0.444 to 0.811), perfluoroundecanoic acid (PFUA) (OR: 0.630, 95% CI: 0.453 to 0.877), n-perfluorooctane sulfonic acid (n-PFOS) (OR: 0.719, 95% CI: 0.571 to 0.906), and perfluoromethylheptane sulfonic acid isomers (Sm-PFOS) (OR: 0.768, 95% CI: 0.602 to 0.981) with gallstones were linearly negative. Qg-comp showed that the PFAS mixture (OR: 0.777, 95% CI: 0.514 to 1.175) was negatively associated with gallstones, but the difference was not statistically significant, and PFDeA had the highest negative association. Moreover, smoking modified the association of perfluorononanoic acid (PFNA) with gallstones. BKMR showed that PFDeA, PFNA, and PFUA had the highest groupPIP (groupPIP = 0.93); PFDeA (condPIP = 0.82), n-perfluorooctanoic acid (n-PFOA) (condPIP = 0.68), and n-PFOS (condPIP = 0.56) also had high condPIPs. Compared with the median level, the joint association of the PFAS mixture with gallstones showed a negative trend; when the PFAS mixture level was at the 70th percentile or higher, they were negatively associated with gallstones. Meanwhile, when other PFAS were fixed at the 25th, 50th, and 75th percentiles, PFDeA had negative associations with gallstones. Our evidence emphasizes that PFAS is negatively associated with gallstones, and more studies are needed in the future to definite the associations of PFAS with gallstones and explore the underlying biological mechanisms.

Prenatal organophosphate esters exposure and neurodevelopment trajectory in infancy: Evidence from the Shanghai Maternal-Child Pairs Cohort

BACKGROUND

Concerns remain about the neurotoxic properties of the ubiquitous organophosphate esters (OPEs), the replacement of the toxicant polybrominated diphenyl ethers.

OBJECTIVES

We examined the associations of prenatal exposure to OPEs and their mixtures with early-life neurodevelopment trajectories.

METHODS

Totally 1276 mother-child pairs were recruited from the Shanghai Maternal-Child Pairs Cohort. A high-performance liquid chromatography-triple quadrupole mass spectrometer was used to measure the levels of 7 OPEs in cord serum. Ages and Stages Questionnaires was used to examine children's neuropsychological development at 2, 6, 12, and 24 months of age. Group-based trajectory models were applied to derive the neurodevelopmental trajectories. Multiple linear regression and logistic regression model were performed to assess the relationships between OPEs exposure and neurodevelopment and trajectories. Mixtures for widely detected OPEs (n = 4) were investigated using quantile-based g-computation.

RESULTS

Tributyl phosphate (TBP), tris (2-butoxy ethyl) phosphate (TBEP), tris(1,3-dichloro-2-propyl) phosphate (TDCPP), and 2-ethylhexyl diphenyl phosphate (EHDPP), had detection rates >50 %. TDCPP had the highest median concentration (1.02 μg/L) in cord serum. EHDPP concentrations were negatively associated with scores in most domains at 12 months of age, with effect values (β) ranging from -1.89 to -0.57. EHDPP could negatively affect the total ASQ (OR = 1.07, 95 % CI: 1, 1.15) and gross-motor (OR = 1.09, 95 % CI: 1.02, 1.17) trajectory in infancy. Joint exposure to OPEs was associated with decreased scores in the total ASQ, gross-motor, fine-motor and problem-solving domain of 12-month-old infants, with β ranging from -5.93 to -1.25. In addition, the qgcomp models indicated significant positive associations between the concentrations of OPEs mixtures and risks of the persistently low group of the total ASQ, gross-motor and fine-motor development in early childhood. The impact of OPEs was more pronounced in boys.

DISCUSSION

Our findings suggested OPEs, especially EHDPP, had a persistently negative effect on neurodevelopment during the first 2 years.