Gaussian graphical modeling of the serum exposome and metabolome reveals interactions between environmental chemicals and endogenous metabolites

Proteomics, Human Environmental Exposure, and Cardiometabolic Risk

BACKGROUND

The biological mechanisms linking environmental exposures with cardiovascular disease pathobiology are incompletely understood. We sought to identify circulating proteomic signatures of environmental exposures and examine their associations with cardiometabolic and respiratory disease in observational cohort studies.

METHODS

We tested the relations of >6500 circulating proteins with 29 environmental exposures across the built environment, green space, air pollution, temperature, and social vulnerability indicators in ≈3000 participants of the CARDIA study (Coronary Artery Risk Development in Young Adults) across 4 centers using penalized and ordinary linear regression. In >3500 participants from FHS (Framingham Heart Study) and JHS (Jackson Heart Study), we evaluated the prospective relations of proteomic signatures of the envirome with cardiovascular disease and mortality using Cox models.

RESULTS

Proteomic signatures of the envirome identified novel/established cardiovascular disease-relevant pathways including DNA damage, fibrosis, inflammation, and mitochondrial function. The proteomic signatures of the envirome were broadly related to cardiometabolic disease and respiratory phenotypes (eg, body mass index, lipids, and left ventricular mass) in CARDIA, with replication in FHS/JHS. A proteomic signature of social vulnerability was associated with a composite of cardiovascular disease/mortality (1428 events; FHS: hazard ratio, 1.16 [95% CI, 1.08-1.24]; P=1.77×10-5; JHS: hazard ratio, 1.25 [95% CI, 1.14-1.38]; P=6.38×10-6; hazard ratio expressed as per 1 SD increase in proteomic signature), robust to adjustment for known clinical risk factors.

CONCLUSIONS

Environmental exposures are related to an inflammatory-metabolic proteome, which identifies individuals with cardiometabolic disease and respiratory phenotypes and outcomes. Future work examining the dynamic impact of the environment on human cardiometabolic health is warranted.

Association of Vitamin D with Perfluorinated Alkyl Acids in Women with and without Non-Obese Polycystic Ovary Syndrome

BACKGROUND

Perfluorinated alkyl acids (PFAAs) are persistent organic pollutants affected by BMI and ethnicity, with contradictory reports of association with vitamin D deficiency.

METHODS

Twenty-nine Caucasian women with non-obese polycystic ovary syndrome (PCOS) and age- and BMI-matched Caucasian control women (n = 30) were recruited. Paired serum samples were analyzed for PFAAs (n = 13) using high-performance liquid chromatography-tandem mass spectrometry. Tandem mass spectrometry determined levels of 25(OH)D3 and the active 1,25(OH)2D3.

RESULTS

Women with and without PCOS did not differ in age, weight, insulin resistance, or systemic inflammation (C-reactive protein did not differ), but the free androgen index was increased. Four PFAAs were detected in all serum samples: perfluorooctane sulfonic acid (PFOS), perfluorooctanoic acid (PFOA), perfluorononanoic acid (PFNA), and perfluorohexane sulfonic acid (PFHxS). Serum PFOS was higher in PCOS versus controls (geometric mean [GM] 3.9 vs. 3.1 ng/mL, p < 0.05). Linear regression modeling showed that elevated PFHxS had higher odds of a lower 25(OH)D3 (OR: 2.919, 95% CI 0.82-5.75, p = 0.04). Vitamin D did not differ between cohorts and did not correlate with any PFAAs, either alone or when the groups were combined. When vitamin D was stratified into sufficiency (>20 ng/mL) and deficiency (<20 ng/mL), no correlation with any PFAAs was seen.

CONCLUSIONS

While the analyses and findings here are exploratory in light of relatively small recruitment numbers, when age, BMI, and insulin resistance are accounted for, the PFAAs do not appear to be related to 25(OH)D3 or the active 1,25(OH)2D3 in this Caucasian population, nor do they appear to be associated with vitamin D deficiency, suggesting that future studies must account for these factors in the analysis.

Exposure to environmental pollutants and genetic variants related to oxidative stress and xenobiotic metabolism-Association with prostate cancer

This study assessed whether genetic variants coding for certain enzymes involved in xenobiotic detoxification, antioxidant defences and DNA repair, along with exposure to environmental chemicals, were associated with an increased prostate cancer (PCa) risk. The study population consisted of 300 men (150 PCa cases and 150 controls) which underwent prostate biopsy as their serum prostate specific antigen (PSA) levels were greater than 4 ng/ml. Genetic variants in GSTM1, GSTP1, SOD2, CAT, GPX1, XRCC1 were determined and data for chemical exposures was obtained through a structured questionnaire and by biomonitoring in a subsample of cases and controls. High serum PSA levels were associated with a greater risk of PCa, while physical exercise appears to exert a protective effect against its development. In addition, elevated urinary levels of certain organic pollutants, such as benzo(a)pyrene (BaP), bisphenol A (BPA), and ethyl-paraben (EPB), were associated with an increased risk of PCa.

Assessment of phthalate exposure at a fire site in Korean firefighters

To determine phthalate exposure in 32 firefighters, the concentrations of urinary phthalate metabolites, immediately (exposure day) and three weeks (control day) after fire suppression, were compared. Mono-(2-ethyl-5-carboxypentyl) phthalate, mono-(2-ethyl-5-hydroxyhexyl) phthalate, mono-(2-ethyl-5-oxohexyl) phthalate, mono-n-butyl phthalate (MBP), mono-n-benzyl phthalate (MBzP), and total phthalates (∑phthalates) levels, and creatinine-adjusted levels of MBP, MBzP, and ∑phthalates were significantly higher on exposure day than on control day. Phthalate concentration was significantly higher in firefighters who performed the fire extinguishing tasks (geometric mean [GM], 149.9 μg/L) than in those who performed other tasks (GM 70.8 μg/L) (p = .012). The GM concentration of firefighters who were active ≤ 50 m from the fire was 119.0 μg/L, and 37.6 μg/L for those who were > 50 m away (p = .012). The GM concentration was significantly different (p = .039) in firefighters with subjective symptoms after fire suppression (151.9 μg/L) compared to those without symptoms (81.6 μg/L). This study showed that firefighters were exposed to phthalate.

Targeted metabolomics reveals plasma biomarkers and metabolic alterations of the aging process in healthy young and older adults

With the exponential growth in the older population in the coming years, many studies have aimed to further investigate potential biomarkers associated with the aging process and its incumbent morbidities. Age is the largest risk factor for chronic disease, likely due to younger individuals possessing more competent adaptive metabolic networks that result in overall health and homeostasis. With aging, physiological alterations occur throughout the metabolic system that contribute to functional decline. In this cross-sectional analysis, a targeted metabolomic approach was applied to investigate the plasma metabolome of young (21-40y; n = 75) and older adults (65y + ; n = 76). A corrected general linear model (GLM) was generated, with covariates of gender, BMI, and chronic condition score (CCS), to compare the metabolome of the two populations. Among the 109 targeted metabolites, those associated with impaired fatty acid metabolism in the older population were found to be most significant: palmitic acid (p < 0.001), 3-hexenedioic acid (p < 0.001), stearic acid (p = 0.005), and decanoylcarnitine (p = 0.036). Derivatives of amino acid metabolism, 1-methlyhistidine (p = 0.035) and methylhistamine (p = 0.027), were found to be increased in the younger population and several novel metabolites were identified, such as cadaverine (p = 0.034) and 4-ethylbenzoic acid (p = 0.029). Principal component analysis was conducted and highlighted a shift in the metabolome for both groups. Receiver operating characteristic analyses of partial least squares-discriminant analysis models showed the candidate markers to be more powerful indicators of age than chronic disease. Pathway and enrichment analyses uncovered several pathways and enzymes predicted to underlie the aging process, and an integrated hypothesis describing functional characteristics of the aging process was synthesized. Compared to older participants, the young group displayed greater abundance of metabolites related to lipid and nucleotide synthesis; older participants displayed decreased fatty acid oxidation and reduced tryptophan metabolism, relative to the young group. As a result, we offer a better understanding of the aging metabolome and potentially reveal new biomarkers and predicted mechanisms for future study.

Urinary Paraben Concentrations and Associations with the Periconceptional Urinary Metabolome: Untargeted and Targeted Metabolomics Analyses of Participants from the Early Pregnancy Study

BACKGROUND

Parabens, found in everyday items from personal care products to foods, are chemicals with endocrine-disrupting activity, which has been shown to influence reproductive function.

OBJECTIVES

This study investigated whether urinary concentrations of methylparaben, propylparaben, or butylparaben were associated with the urinary metabolome during the periconceptional period, a critical window for female reproductive function. Changes to the periconceptional urinary metabolome could provide insights into the mechanisms by which parabens could impact fertility.

METHODS

Urinary paraben concentrations were measured in paired pre- and postconception urine samples from 42 participants in the Early Pregnancy Study, a prospective cohort of 221 women attempting to conceive. We performed untargeted and targeted metabolomics analyses using ultrahigh-performance liquid chromatography quadrupole time-of-flight mass spectrometry. We used principal component analysis, orthogonal partial least-squares discriminant analysis, and permutation testing, coupled with univariate statistical analyses, to find metabolites associated with paraben concentration at the two time points. Potential confounders were identified with a directed acyclic graph and used to adjust results with multivariable linear regression. Metabolites were identified using fragmentation data.

RESULTS

Seven metabolites were associated with paraben concentration (variable importance to projection score > 1 , false discovery rate-corrected q -value < 0.1 ). We identified four diet-related metabolites to the Metabolomics Standards Initiative (MSI) certainty of identification level 2, including metabolites from smoke flavoring, grapes, and olive oil. One metabolite was identified to the class level only (MSI level 3). Two metabolites were unidentified (MSI level 4). After adjustment, three metabolites remained associated with methylparaben and propylparaben, two of which were diet-related. No metabolomic markers of endocrine disruption were associated with paraben concentrations.

DISCUSSION

This study identified novel relationships between urinary paraben concentrations and diet-related metabolites but not with metabolites on endocrine-disrupting pathways, as hypothesized. It demonstrates the feasibility of integrating untargeted metabolomics data with environmental exposure information and epidemiological adjustment for confounders. The findings underscore a potentially important connection between diet and paraben exposure, with applications to nutritional epidemiology and dietary exposure assessment. https://doi.org/10.1289/EHP12125.

PFAS Exposures and the Human Metabolome: A Systematic Review of Epidemiological Studies

Purpose of Review

There is a growing interest in understanding the health effects of exposure to per- and polyfluoroalkyl substances (PFAS) through the study of the human metabolome. In this systematic review, we aimed to identify consistent findings between PFAS and metabolomic signatures. We conducted a search matching specific keywords that was independently reviewed by two authors on two databases (EMBASE and PubMed) from their inception through July 19, 2022 following PRISMA guidelines.

Recent Findings

We identified a total of 28 eligible observational studies that evaluated the associations between 31 different PFAS exposures and metabolomics in humans. The most common exposure evaluated was legacy long-chain PFAS. Population sample sizes ranged from 40 to 1,105 participants at different stages across the lifespan. A total of 19 studies used a non-targeted metabolomics approach, 7 used targeted approaches, and 2 included both. The majority of studies were cross-sectional (n = 25), including four with prospective analyses of PFAS measured prior to metabolomics.

Summary

Most frequently reported associations across studies were observed between PFAS and amino acids, fatty acids, glycerophospholipids, glycerolipids, phosphosphingolipids, bile acids, ceramides, purines, and acylcarnitines. Corresponding metabolic pathways were also altered, including lipid, amino acid, carbohydrate, nucleotide, energy metabolism, glycan biosynthesis and metabolism, and metabolism of cofactors and vitamins. We found consistent evidence across studies indicating PFAS-induced alterations in lipid and amino acid metabolites, which may be involved in energy and cell membrane disruption.

Firefighters and the liver: Exposure to PFAS and PAHs in relation to liver function and serum lipids (CELSPAC-FIREexpo study)

INTRODUCTION

Firefighting is one of the most hazardous occupations due to exposure to per- and polyfluoroalkyl substances (PFAS) and polycyclic aromatic hydrocarbons (PAHs). Such exposure is suspected to affect the cardiometabolic profile, e.g., liver function and serum lipids. However, only a few studies have investigated the impact of this specific exposure among firefighters.

METHODS

Men included in the CELSPAC-FIREexpo study were professional firefighters (n = 52), newly recruited firefighters in training (n = 58), and controls (n = 54). They completed exposure questionnaires and provided 1-3 samples of urine and blood during the 11-week study period to allow assessment of their exposure to PFAS (6 compounds) and PAHs (6 compounds), and to determine biomarkers of liver function (alanine aminotransferase (ALT), gamma-glutamyl transferase (GGT), aspartate aminotransferase (AST), alkaline phosphatase (ALP) and total bilirubin (BIL)) and levels of serum lipids (total cholesterol (CHOL), low-density lipoprotein cholesterol (LDL) and triglycerides (TG)). The associations between biomarkers were investigated both cross-sectionally using multiple linear regression (MLR) and Bayesian weighted quantile sum (BWQS) regression and prospectively using MLR. The models were adjusted for potential confounders and false discovery rate correction was applied to account for multiplicity.

RESULTS

A positive association between exposure to PFAS and PAH mixture and BIL (β = 28.6%, 95% CrI = 14.6-45.7%) was observed by the BWQS model. When the study population was stratified, in professional firefighters and controls the mixture showed a positive association with CHOL (β = 29.5%, CrI = 10.3-53.6%) and LDL (β = 26.7%, CrI = 8.3-48.5%). No statistically significant associations with individual compounds were detected using MLR.

CONCLUSIONS

This study investigated the associations between exposure to PFAS and PAHs and biomarkers of cardiometabolic health in the Czech men, including firefighters. The results suggest that higher exposure to a mixture of these compounds is associated with an increase in BIL and the alteration of serum lipids, which can result in an unfavourable cardiometabolic profile.

Extending Nontargeted Discovery of Environmental Chemical Exposures during Pregnancy and Their Association with Pregnancy Complications-A Cross-Sectional Study

BACKGROUND

Nontargeted analysis (NTA) methods identify novel exposures; however, few chemicals have been quantified and interrogated with pregnancy complications.

OBJECTIVES

We characterized levels of nine exogenous and endogenous chemicals in maternal and cord blood identified, selected, and confirmed in prior NTA steps, including linear and branched isomers perfluorooctane sulfonate (PFOS), perfluorohexane sulfonate (PFHxS), monoethylhexyl phthalate, 4-nitrophenol, tetraethylene glycol, tridecanedioic acid, octadecanedioic acid, and deoxycholic acid. We evaluated relationships between maternal and cord levels and between gestational diabetes mellitus (GDM) and hypertensive disorders of pregnancy in a diverse pregnancy cohort in San Francisco.

METHODS

We collected matched maternal and cord serum samples at delivery from 302 pregnant study participants from the Chemicals in Our Bodies cohort in San Francisco. Chemicals were identified via NTA and quantified using targeted approaches. We calculated distributions and Spearman correlation coefficients testing the relationship of chemicals within and between the maternal and cord blood matrices. We used adjusted logistic regression to calculate the odds of GDM and hypertensive disorders of pregnancy associated with an interquartile range increase in maternal chemical exposures.

RESULTS

We detected linear PFOS, PFHxS, octadecanedioic acid, and deoxycholic acid in at least 97% of maternal samples. Correlations ranged between - 0.1 and 0.9. We observed strong correlations between cord and maternal levels of PFHxS, linear PFOS, and branched PFOS (coefficient = 0.9 , 0.8, and 0.8, respectively). An interquartile range increase in linear and branched PFOS, tridecanedioic acid, octadecanedioic acid, and deoxycholic acid was associated with increased odds ratio (OR) of GDM [OR = 1.33 (95% CI: 0.89, 2.01), 1.24 (95% CI: 0.86, 1.80), 1.26 (95% CI: 0.93, 1.73), 1.24 (95% CI: 0.86, 1.80), and 1.23 (95% CI: 0.87, 1.75), respectively]. Tridecanedioic acid was positively associated with hypertensive disorders of pregnancy [OR = 1.28 (95% CI: 0.90, 1.86)].

DISCUSSION

We identified both exogenous and endogenous chemicals seldom quantified in pregnant study participants that were also related to pregnancy complications and demonstrated the utility of NTA to identify chemical exposures of concern. https://doi.org/10.1289/EHP11546.

Childhood exposure to non-persistent endocrine disrupting chemicals and multi-omic profiles: A panel study

BACKGROUND

Individuals are exposed to environmental pollutants with endocrine disrupting activity (endocrine disruptors, EDCs) and the early stages of life are particularly susceptible to these exposures. Previous studies have focused on identifying molecular signatures associated with EDCs, but none have used repeated sampling strategy and integrated multiple omics. We aimed to identify multi-omic signatures associated with childhood exposure to non-persistent EDCs.

METHODS

We used data from the HELIX Child Panel Study, which included 156 children aged 6 to 11. Children were followed for one week, in two time periods. Twenty-two non-persistent EDCs (10 phthalate, 7 phenol, and 5 organophosphate pesticide metabolites) were measured in two weekly pools of 15 urine samples each. Multi-omic profiles (methylome, serum and urinary metabolome, proteome) were measured in blood and in a pool urine samples. We developed visit-specific Gaussian Graphical Models based on pairwise partial correlations. The visit-specific networks were then merged to identify reproducible associations. Independent biological evidence was systematically sought to confirm some of these associations and assess their potential health implications.

RESULTS

950 reproducible associations were found among which 23 were direct associations between EDCs and omics. For 9 of them, we were able to find corroborating evidence from previous literature: DEP - serotonin, OXBE - cg27466129, OXBE - dimethylamine, triclosan - leptin, triclosan - serotonin, MBzP - Neu5AC, MEHP - cg20080548, oh-MiNP - kynurenine, oxo-MiNP - 5-oxoproline. We used these associations to explore possible mechanisms between EDCs and health outcomes, and found links to health outcomes for 3 analytes: serotonin and kynurenine in relation to neuro-behavioural development, and leptin in relation to obesity and insulin resistance.

CONCLUSIONS

This multi-omics network analysis at two time points identified biologically relevant molecular signatures related to non-persistent EDC exposure in childhood, suggesting pathways related to neurological and metabolic outcomes.

Perfluoroalkyl substances and lipid concentrations in the blood: A systematic review of epidemiological studies

BACKGROUND

Perfluoroalkyl substances (PFAS) are widely used synthetic aliphatic compounds. This systematic review aims to assess PFAS associations with low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL), total cholesterol (TC) and total triglyceride (TG) concentrations in human populations.

METHOD

We systematically searched four online databases, PubMed, Scopus, Embase, and Cochrane Library for relevant peer-reviewed English language articles published until July 2021. Additional relevant articles identified were also included in the search results. We categorised populations into adults (≥18 years old) and children. Primary findings were the associations between PFAS concentrations and LDL, HDL, TC, and TG concentrations in the serum, plasma, or whole blood; secondary findings were the associations between PFAS concentrations and the odds of lipid-related health outcomes. Quantitative synthesis was done by vote counting of the effect directions between concentrations of PFAS and lipids/health outcomes, repeated on articles with sample size >1000. Sign tests were performed to assess the statistical significance of the differences between positive and negative associations. Sensitivity analysis was performed by separating out articles with populations having high concentrations of perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS). Quality was assessed with the STROBE checklist and NHBLI Study Quality Assessment Tool.

RESULTS

A total of 58 articles were included for review. There was evidence that PFAS exposure is associated with higher concentrations of LDL, HDL, and TC, particularly for PFOA-LDL, PFOA-TC, PFOS-TC, and PFNA-LDL. Associations between PFAS and TG tended to be negative, especially for perfluoroundecanoic acid (PFUnDA). Associations between PFAS concentration and the odds of secondary outcomes generally supported a positive association between PFAS and cholesterol concentrations.

CONCLUSION

We found evidence of associations between the concentrations of some PFAS-lipid pairs in human populations. Future research should be conducted on the less well-studied PFAS to explore their effects on human health and in regions where such studies are currently lacking. (300 words).

State-of-the-art methods for exposure-health studies: Results from the exposome data challenge event

The exposome recognizes that individuals are exposed simultaneously to a multitude of different environmental factors and takes a holistic approach to the discovery of etiological factors for disease. However, challenges arise when trying to quantify the health effects of complex exposure mixtures. Analytical challenges include dealing with high dimensionality, studying the combined effects of these exposures and their interactions, integrating causal pathways, and integrating high-throughput omics layers. To tackle these challenges, the Barcelona Institute for Global Health (ISGlobal) held a data challenge event open to researchers from all over the world and from all expertises. Analysts had a chance to compete and apply state-of-the-art methods on a common partially simulated exposome dataset (based on real case data from the HELIX project) with multiple correlated exposure variables (P > 100 exposure variables) arising from general and personal environments at different time points, biological molecular data (multi-omics: DNA methylation, gene expression, proteins, metabolomics) and multiple clinical phenotypes in 1301 mother-child pairs. Most of the methods presented included feature selection or feature reduction to deal with the high dimensionality of the exposome dataset. Several approaches explicitly searched for combined effects of exposures and/or their interactions using linear index models or response surface methods, including Bayesian methods. Other methods dealt with the multi-omics dataset in mediation analyses using multiple-step approaches. Here we discuss features of the statistical models used and provide the data and codes used, so that analysts have examples of implementation and can learn how to use these methods. Overall, the exposome data challenge presented a unique opportunity for researchers from different disciplines to create and share state-of-the-art analytical methods, setting a new standard for open science in the exposome and environmental health field.

Molecular Gatekeeper Discovery: Workflow for Linking Multiple Exposure Biomarkers to Metabolomics

The exposome reflects multiple exposures across the life-course that can affect health. Metabolomics can reveal the underlying molecular basis linking exposures to health conditions. Here, we explore the concept and general data analysis framework of "molecular gatekeepers"─key metabolites that link single or multiple exposure biomarkers with correlated clusters of endogenous metabolites─to inform health-relevant biological targets. We performed untargeted metabolomics on plasma from 152 adolescent girls participating in the Growing Up Healthy Study in New York City. We then performed network analysis to link metabolites to exposure biomarkers including five trace elements (Cd, Mn, Pb, Se, and Hg) and five perfluorinated chemicals (PFCs; n-PFOS, Sm-PFOS, n-PFOA, PFHxS, and PFNA). We found 144 molecular gatekeepers and annotated 22 of them. Lysophosphatidylcholine (16:0) and taurodeoxycholate were correlated with both n-PFOA and n-PFOS, suggesting a shared dysregulation from multiple xenobiotic exposures. Sphingomyelin (d18:2/14:0) was significantly associated with age at menarche; yet, no direct association was detected between any exposure biomarkers and age at menarche. Thus, molecular gatekeepers can also discover molecular linkages between exposure biomarkers and health outcomes that may otherwise be obscured by complex interactions in direct measurements.

[Chemical exposome and non-targeted approaches]

The technological advances in high-resolution mass spectrometry (HRMS), associated with the development of bioinformatics tools, allows the simultaneous detection of tens of thousands of chemical signals in biological matrices, including exogenous (i.e. xenobiotics) and endogenous molecules. These novel approaches based on HRMS, called "non-targeted" approaches, provide a unique opportunity to capture exposures to a wide range of chemicals (i.e. the internal chemical exposome) in populations, and to better understand the links between chemical exposures and the occurrence of chronic diseases.