A multivariate approach to investigate the combined biological effects of multiple exposures

Ten-Year Risk Equations for Incident Heart Failure in Established Atherosclerotic Cardiovascular Disease Populations

BACKGROUND

Ten-year risk equations for incident heart failure (HF) are available for the general population, but not for patients with established atherosclerotic cardiovascular disease (ASCVD), which is highly prevalent in HF cohorts. This study aimed to develop and validate 10-year risk equations for incident HF in patients with known ASCVD.

METHODS AND RESULTS

Ten-year risk equations for incident HF were developed using the United Kingdom Biobank cohort (recruitment 2006-2010) including participants with established ASCVD but free from HF at baseline. Model performance was validated using the Australian Baker Heart and Diabetes Institute Biobank cohort (recruitment 2000-2011) and compared with the performance of general population risk models. Incident HF occurred in 13.7% of the development cohort (n=31 446, median 63 years, 35% women, follow-up 10.7±2.7 years) and in 21.3% of the validation cohort (n=1659, median age 65 years, 25% women, follow-up 9.4±3.7 years). Predictors of HF included in the sex-specific models were age, body mass index, systolic blood pressure (treated or untreated), glucose (treated or untreated), cholesterol, smoking status, QRS duration, kidney disease, myocardial infarction, and atrial fibrillation. ASCVD-HF equations had good discrimination and calibration in development and validation cohorts, with superior performance to general population risk equations.

CONCLUSIONS

ASCVD-specific 10-year risk equations for HF outperform general population risk models in individuals with established ASCVD. The ASCVD-HF equations can be calculated from readily available clinical data and could facilitate screening and preventative treatment decisions in this high-risk group.

No-Reflow Prediction in Acute Coronary Syndrome During Percutaneous Coronary Intervention: The NORPACS Risk Score

BACKGROUND

Suboptimal coronary reperfusion (no reflow) is common in acute coronary syndrome percutaneous coronary intervention (PCI) and is associated with poor outcomes. We aimed to develop and externally validate a clinical risk score for angiographic no reflow for use following angiography and before PCI.

METHODS

We developed and externally validated a logistic regression model for prediction of no reflow among adult patients undergoing PCI for acute coronary syndrome using data from the Melbourne Interventional Group PCI registry (2005-2020; development cohort) and the British Cardiovascular Interventional Society PCI registry (2006-2020; external validation cohort).

RESULTS

A total of 30 561 patients (mean age, 64.1 years; 24% women) were included in the Melbourne Interventional Group development cohort and 440 256 patients (mean age, 64.9 years; 27% women) in the British Cardiovascular Interventional Society external validation cohort. The primary outcome (no reflow) occurred in 4.1% (1249 patients) and 9.4% (41 222 patients) of the development and validation cohorts, respectively. From 33 candidate predictor variables, 6 final variables were selected by an adaptive least absolute shrinkage and selection operator regression model for inclusion (cardiogenic shock, ST-segment-elevation myocardial infarction with symptom onset >195 minutes pre-PCI, estimated stent length ≥20 mm, vessel diameter <2.5 mm, pre-PCI Thrombolysis in Myocardial Infarction flow <3, and lesion location). Model discrimination was very good (development C statistic, 0.808; validation C statistic, 0.741) with excellent calibration. Patients with a score of ≥8 points had a 22% and 27% risk of no reflow in the development and validation cohorts, respectively.

CONCLUSIONS

The no-reflow prediction in acute coronary syndrome risk score is a simple count-based scoring system based on 6 parameters available before PCI to predict the risk of no reflow. This score could be useful in guiding preventative treatment and future trials.

Seminar: Scalable Preprocessing Tools for Exposomic Data Analysis

BACKGROUND

The exposome serves as a popular framework in which to study exposures from chemical and nonchemical stressors across the life course and the differing roles that these exposures can play in human health. As a result, data relevant to the exposome have been used as a resource in the quest to untangle complicated health trajectories and help connect the dots from exposures to adverse outcome pathways.

OBJECTIVES

The primary aim of this methods seminar is to clarify and review preprocessing techniques critical for accurate and effective external exposomic data analysis. Scalability is emphasized through an application of highly innovative combinatorial techniques coupled with more traditional statistical strategies. The Public Health Exposome is used as an archetypical model. The novelty and innovation of this seminar's focus stem from its methodical, comprehensive treatment of preprocessing and its demonstration of the positive effects preprocessing can have on downstream analytics.

DISCUSSION

State-of-the-art technologies are described for data harmonization and to mitigate noise, which can stymie downstream interpretation, and to select key exposomic features, without which analytics may lose focus. A main task is the reduction of multicollinearity, a particularly formidable problem that frequently arises from repeated measurements of similar events taken at various times and from multiple sources. Empirical results highlight the effectiveness of a carefully planned preprocessing workflow as demonstrated in the context of more highly concentrated variable lists, improved correlational distributions, and enhanced downstream analytics for latent relationship discovery. The nascent field of exposome science can be characterized by the need to analyze and interpret a complex confluence of highly inhomogeneous spatial and temporal data, which may present formidable challenges to even the most powerful analytical tools. A systematic approach to preprocessing can therefore provide an essential first step in the application of modern computer and data science methods. https://doi.org/10.1289/EHP12901.

A weighted quantile sum regression with penalized weights and two indices

Background

New statistical methodologies were developed in the last decade to face the challenges of estimating the effects of exposure to multiple chemicals. Weighted Quantile Sum (WQS) regression is a recent statistical method that allows estimating a mixture effect associated with a specific health effect and identifying the components that characterize the mixture effect.

Objectives

In this study, we propose an extension of WQS regression that estimates two mixture effects of chemicals on a health outcome in the same model through the inclusion of two indices, one in the positive direction and one in the negative direction, with the introduction of a penalization term.

Methods

To evaluate the performance of this new model we performed both a simulation study and a real case study where we assessed the effects of nutrients on obesity among adults using the National Health and Nutrition Examination Survey (NHANES) data.

Results

The method showed good performance in estimating both the regression parameter and the weights associated with the single elements when the penalized term was set equal to the magnitude of the Akaike information criterion of the unpenalized WQS regression. The two indices further helped to give a better estimate of the parameters [Positive direction Median Error (PME): 0.022; Negative direction Median Error (NME): -0.044] compared to the standard WQS without the penalization term (PME: -0.227; NME: 0.215). In the case study, WQS with two indices was able to find a significant effect of nutrients on obesity in both directions identifying sodium and magnesium as the main actors in the positive and negative association, respectively.

Discussion

Through this work, we introduced an extension of WQS regression that improved the accuracy of the parameter estimates when considering a mixture of elements that can have both a protective and a harmful effect on the outcome; and the advantage of adding a penalization term when estimating the weights.

Multi-omics signatures of the human early life exposome

Environmental exposures during early life play a critical role in life-course health, yet the molecular phenotypes underlying environmental effects on health are poorly understood. In the Human Early Life Exposome (HELIX) project, a multi-centre cohort of 1301 mother-child pairs, we associate individual exposomes consisting of >100 chemical, outdoor, social and lifestyle exposures assessed in pregnancy and childhood, with multi-omics profiles (methylome, transcriptome, proteins and metabolites) in childhood. We identify 1170 associations, 249 in pregnancy and 921 in childhood, which reveal potential biological responses and sources of exposure. Pregnancy exposures, including maternal smoking, cadmium and molybdenum, are predominantly associated with child DNA methylation changes. In contrast, childhood exposures are associated with features across all omics layers, most frequently the serum metabolome, revealing signatures for diet, toxic chemical compounds, essential trace elements, and weather conditions, among others. Our comprehensive and unique resource of all associations ( https://helixomics.isglobal.org/ ) will serve to guide future investigation into the biological imprints of the early life exposome.

A review on the application of the exposome paradigm to unveil the environmental determinants of age-related diseases

Age-related diseases account for almost half of all diseases among adults worldwide, and their incidence is substantially affected by the exposome, which is the sum of all exogenous and endogenous environmental exposures and the human body’s response to these exposures throughout the entire lifespan. Herein, we perform a comprehensive review of the epidemiological literature to determine the key elements of the exposome that affect the development of age-related diseases and the roles of aging hallmarks in this process. We find that most exposure assessments in previous aging studies have used a reductionist approach, whereby the effect of only a single environmental factor or a specific class of environmental factors on the development of age-related diseases has been examined. As such, there is a lack of a holistic and unbiased understanding of the effect of multiple environmental factors on the development of age-related diseases. To address this, we propose several research strategies based on an exposomic framework that could advance our understanding—in particular, from a mechanistic perspective—of how environmental factors affect the development of age-related diseases. We discuss the statistical methods and other methods that have been used in exposome-wide association studies, with a particular focus on multiomics technologies. We also address future challenges and opportunities in the realm of multidisciplinary approaches and genome–exposome epidemiology. Furthermore, we provide perspectives on precise public health services for vulnerable populations, public communications, the integration of risk exposure information, and the bench-to-bedside translation of research on age-related diseases.

Tracing the oomycete pathogen Saprolegnia parasitica in aquaculture and the environment

Saprolegnia parasitica causes saprolegniosis, a disease responsible for significant economic losses in aquaculture and declines of fish populations in the wild, but the knowledge of its distribution and prevalence in the environment is limited. We developed a fast, sensitive and specific S. parasitica droplet digital PCR (ddPCR) assay and demonstrated its applicability for the detection and quantification of the pathogen in environmental samples: swab DNA collected from the host (trout skin, surface of eggs) and environmental DNA extracted from water. The developed assay was used to assess how abiotic (i.e. physico-chemical parameters of the water) and biotic (health status of the host) factors influence the S. parasitica load in the environment. The pathogen load in water samples was positively correlated with some site-specific abiotic parameters such as electrical conductivity (EC) and calcium, while fluorides were negatively correlated, suggesting that physico-chemical parameters are important for determining S. parasitica load in natural waters. Furthermore, skin swabs of injured trout had significantly higher pathogen load than swabs collected from healthy fish, confirming that S. parasitica is a widespread opportunistic pathogen. Our results provide new insights into various environmental factors that influence the distribution and abundance of S. parasitica.

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.

Variations in the Sporulation Efficiency of Pathogenic Freshwater Oomycetes in Relation to the Physico-Chemical Properties of Natural Waters

Oomycete pathogens in freshwaters, such as Saprolegnia parasitica and Aphanomyces astaci, are responsible for fish/crayfish population declines in the wild and disease outbreaks in aquaculture. Although the formation of infectious zoospores in the laboratory can be triggered by washing their mycelium with natural water samples, the physico-chemical properties of the water that might promote sporulation are still unexplored. We washed the mycelia of A. astaci and S. parasitica with a range of natural water samples and observed differences in sporulation efficiency. The results of Partial Least Squares Regression (PLS-R) multivariate analysis showed that SAC (spectral absorption coefficient measured at 254 nm), DOC (dissolved organic carbon), ammonium-N and fluoride had the strongest positive effect on sporulation of S. parasitica, while sporulation of A. astaci was not significantly correlated with any of the analyzed parameters. In agreement with this, the addition of environmentally relevant concentrations of humic acid, an important contributor to SAC and DOC, to the water induced sporulation of S. parasitica but not of A. astaci. Overall, our results point to the differences in ecological requirements of these pathogens, but also present a starting point for optimizing laboratory protocols for the induction of sporulation.

Associations between persistent organic pollutants and endometriosis: A multiblock approach integrating metabolic and cytokine profiling

Humans are exposed daily to complex mixtures of chemical pollutants through their environment and diet, some of which have the potential to disrupt the bodies' natural endocrine functions and contribute to reproductive diseases like endometriosis. Increasing epidemiological and experimental evidence supports the association between endometriosis and certain persistent organic pollutants (POPs) like dioxins; however, little is known about the underlying linking mechanisms. The main objective of this study is to proof the methodological applicability and discovery potential of integrating ultra-trace mass spectrometry (MS) profiling of POP biomarkers and endogenous biomarker profiling (MS metabolomics and cytokines) in a case-control study for the etiological research of endometriosis. The approach is applied in a pilot clinical-based study conducted in France where women with and without surgically confirmed endometriosis were recruited. Serum samples were analysed with high-resolution MS for about 30 polychlorinated biphenyls (PCBs), organochlorinated pesticides and perfluoroalkyl substances (PFAS). About 600 serum metabolites and lipids were identified with targeted metabolomics using tandem MS with the Biocrates MxP® Quant 500 Kit. A panel of 4 pro-inflammatory cytokines were analysed using ELISA-based 4-PLEX analyser. Statistical analysis included a battery of variable selection approaches, multivariate logistic regression for single-chemical associations, Bayesian kernel machine regressions (BKMR) to identify mixture effects of POPs and a multiblock approach to identify shared biomarker signatures among high risk clusters. The results showed the positive associations between some POPs and endometriosis risk, including the pesticide trans-nonachlor Odds Ratio (95% Confidence Interval) 3.38 (2.06-5.98), p < 0.0001 and PCB 114 OR (95% CI) 1.83 (1.17-2.93), p = 0.009. The BKMR approach showed a tendency of a positive cumulative effect of the mixture, however trans-nonachlor exhibited significant associations within the mixture and interacted with other PCBs, strengthening the effects at highest concentrations. Finally, the multiblock analysis, relating the various blocks of data, revealed a latent cluster of women with higher risk of endometrioma presenting higher concentrations of trans-nonachlor, PCB 114 and dioxin-like toxic equivalents from PCBs, together with an increased inflammatory profile (i.e. elevated interleukin-8 and monocyte chemoattractant protein-1). It was also highlighted a specific metabolic pattern characterized by dysregulation of bile acid homeostasis and lipase activity. Further research will be required with larger sample size to confirm these findings and gain insight on the underlying mechanisms between POPs and endometriosis.

Immune Response in Crayfish Is Species-Specific and Exhibits Changes along Invasion Range of a Successful Invader

Immunity is an important component of invasion success since it enables invaders' adaptation to conditions of the novel environment as they expand their range. Immune response of invaders may vary along the invasion range due to encountered parasites/microbial communities, conditions of the local environment, and ecological processes that arise during the range expansion. Here, we analyzed changes in the immune response along the invasion range of one of the most successful aquatic invaders, the signal crayfish, in the recently invaded Korana River, Croatia. We used several standard immune parameters (encapsulation response, hemocyte count, phenoloxidaze activity, and total prophenoloxidaze) to: i) compare immune response of the signal crayfish along its invasion range, and between species (comparison with co-occurring native narrow-clawed crayfish), and ii) analyze effects of specific predictors (water temperature, crayfish abundance, and body condition) on crayfish immune response changes. Immune response displayed species-specificity, differed significantly along the signal crayfish invasion range, and was mostly affected by water temperature and population abundance. Specific immune parameters showed density-dependent variation corresponding to increased investment in them during range expansion. Obtained results offer baseline insights for elucidating the role of immunocompetence in the invasion success of an invertebrate freshwater invader.

Prenatal Exposure to Multiple Air Pollutants, Mediating Molecular Mechanisms, and Shifts in Birthweight

Mechanisms underlying adverse birth and later in life health effects from exposure to air pollution during the prenatal period have not been not fully elucidated, especially in the context of mixtures. We assessed the effects of prenatal exposure to mixtures of air pollutants of particulate matter (PM), PM2.5, PM10, nitrogen oxides, NO2, NOx, ultrafine particles (UFP), and oxidative potential (OP) of PM2.5 on infant birthweight in four European birth cohorts and the mechanistic underpinnings through cross-omics of metabolites and inflammatory proteins. The association between mixtures of air pollutants and birthweight z-scores (standardized for gestational age) was assessed for three different mixture models, using Bayesian machine kernel regression (BKMR). We determined the direct effect for PM2.5, PM10, NO2, and mediation by cross-omic signatures (identified using sparse partial least-squares regression) using causal mediation BKMR models. There was a negative association with birthweight z-scores and exposure to mixtures of air pollutants, where up to -0.21 or approximately a 96 g decrease in birthweight, comparing the 75th percentile to the median level of exposure to the air pollutant mixture could occur. Shifts in birthweight z-scores from prenatal exposure to PM2.5, PM10, and NO2 were mediated by molecular mechanisms, represented by cross-omics scores. Interleukin-17 and epidermal growth factor were identified as important inflammatory responses underlyingair pollution-associated shifts in birthweight. Our results signify that by identifying mechanisms through which mixtures of air pollutants operate, the causality of air pollution-associated shifts in birthweight is better supported, substantiating the need for reducing exposure in vulnerable populations.

Bayesian Profile Regression to Deal With Multiple Highly Correlated Exposures and a Censored Survival Outcome. First Application in Ionizing Radiation Epidemiology

As multifactorial and chronic diseases, cancers are among these pathologies for which the exposome concept is essential to gain more insight into the associated etiology and, ultimately, lead to better primary prevention strategies for public health. Indeed, cancers result from the combined influence of many genetic, environmental and behavioral stressors that may occur simultaneously and interact. It is thus important to properly account for multifactorial exposure patterns when estimating specific cancer risks at individual or population level. Nevertheless, the risk factors, especially environmental, are still too often considered in isolation in epidemiological studies. Moreover, major statistical difficulties occur when exposures to several factors are highly correlated due, for instance, to common sources shared by several pollutants. Suitable statistical methods must then be used to deal with these multicollinearity issues. In this work, we focused on the specific problem of estimating a disease risk from highly correlated environmental exposure covariates and a censored survival outcome. We extended Bayesian profile regression mixture (PRM) models to this context by assuming an instantaneous excess hazard ratio disease sub-model. The proposed hierarchical model incorporates an underlying truncated Dirichlet process mixture as an attribution sub-model. A specific adaptive Metropolis-Within-Gibbs algorithm-including label switching moves-was implemented to infer the model. This allows simultaneously clustering individuals with similar risks and similar exposure characteristics and estimating the associated risk for each group. Our Bayesian PRM model was applied to the estimation of the risk of death by lung cancer in a cohort of French uranium miners who were chronically and occupationally exposed to multiple and correlated sources of ionizing radiation. Several groups of uranium miners with high risk and low risk of death by lung cancer were identified and characterized by specific exposure profiles. Interestingly, our case study illustrates a limit of MCMC algorithms to fit full Bayesian PRM models even if the updating schemes for the cluster labels incorporate label-switching moves. Then, although this paper shows that Bayesian PRM models are promising tools for exposome research, it also opens new avenues for methodological research in this class of probabilistic models.

Associations between persistent organic pollutants and endometriosis: A multipollutant assessment using machine learning algorithms

Endometriosis is a gynaecological disease characterised by the presence of endometriotic tissue outside of the uterus impacting a significant fraction of women of childbearing age. Evidence from epidemiological studies suggests a relationship between risk of endometriosis and exposure to some organochlorine persistent organic pollutants (POPs). However, these chemicals are numerous and occur in complex and highly correlated mixtures, and to date, most studies have not accounted for this simultaneous exposure. Linear and logistic regression models are constrained to adjusting for multiple exposures when variables are highly intercorrelated, resulting in unstable coefficients and arbitrary findings. Advanced machine learning models, of emerging use in epidemiology, today appear as a promising option to address these limitations. In this study, different machine learning techniques were compared on a dataset from a case-control study conducted in France to explore associations between mixtures of POPs and deep endometriosis. The battery of models encompassed regularised logistic regression, artificial neural network, support vector machine, adaptive boosting, and partial least-squares discriminant analysis with some additional sparsity constraints. These techniques were applied to identify the biomarkers of internal exposure in adipose tissue most associated with endometriosis and to compare model classification performance. The five tested models revealed a consistent selection of most associated POPs with deep endometriosis, including octachlorodibenzofuran, cis-heptachlor epoxide, polychlorinated biphenyl 77 or trans-nonachlor, among others. The high classification performance of all five models confirmed that machine learning may be a promising complementary approach in modelling highly correlated exposure biomarkers and their associations with health outcomes. Regularised logistic regression provided a good compromise between the interpretability of traditional statistical approaches and the classification capacity of machine learning approaches. Applying a battery of complementary algorithms may be a strategic approach to decipher complex exposome-health associations when the underlying structure is unknown.

Applying the exposome concept in birth cohort research: a review of statistical approaches

The exposome represents the totality of life course environmental exposures (including lifestyle and other non-genetic factors), from the prenatal period onwards. This holistic concept of exposure provides a new framework to advance the understanding of complex and multifactorial diseases. Prospective pregnancy and birth cohort studies provide a unique opportunity for exposome research as they are able to capture, from prenatal life onwards, both the external (including lifestyle, chemical, social and wider community-level exposures) and the internal (including inflammation, metabolism, epigenetics, and gut microbiota) domains of the exposome. In this paper, we describe the steps required for applying an exposome approach, describe the main strengths and limitations of different statistical approaches and discuss their challenges, with the aim to provide guidance for methodological choices in the analysis of exposome data in birth cohort studies. An exposome approach implies selecting, pre-processing, describing and analyzing a large set of exposures. Several statistical methods are currently available to assess exposome-health associations, which differ in terms of research question that can be answered, of balance between sensitivity and false discovery proportion, and between computational complexity and simplicity (parsimony). Assessing the association between many exposures and health still raises many exposure assessment issues and statistical challenges. The exposome favors a holistic approach of environmental influences on health, which is likely to allow a more complete understanding of disease etiology.

How polluted is your fat? What the study of adipose tissue can contribute to environmental epidemiology

The study of the potential contribution of low-dose exposure to environmental chemicals on the development of chronic conditions in human populations is often hampered by methodological issues, including exposure misclassification and the inability to assess biological effects in target organs. White adipose tissue (WAT) presents the unique feature of being both an advantageous matrix for assessing long-term exposure to mixtures of persistent organic pollutants and an interesting tissue to investigate early preclinical effects. Moreover, other lipophilic non-persistent chemicals and heavy metals have been recently quantified in fat, suggesting that human WAT contains chemical mixtures more complex than initially thought. However, WAT has been scarcely used in environmental epidemiology due to collection difficulties. In this essay we discuss the potential of using human WAT as a source of both exposure and effect biomarkers, with the aim of advancing the epidemiological research of obesity-related diseases, including metabolic syndrome and cancer. Overall, we discuss the implications of investigating WAT in a multidisciplinary framework combining toxicological and epidemiological knowledge in order to improve the inference of causal relationships in observational settings. We finalise by suggesting feasible designs and scenarios in which WAT samples may be reasonably collected.

Analyzing Metabolomics Data for Environmental Health and Exposome Research

The exposome is the cumulative measure of environmental influences and associated biological responses across the life span, with critical relevance for understanding how exposures can impact human health. Metabolomics analysis of biological samples offers unique advantages for examining the exposome. Simultaneous analysis of external exposures, biological responses, and host susceptibility at a systems level can help establish links between external exposures and health outcomes. As metabolomics technologies continue to evolve for the study of the exposome, metabolomics ultimately will help provide valuable insights for exposure risk assessment, and disease prevention and management. Here, we discuss recent advances in metabolomics, and describe data processing protocols that can enable analysis of the exposome. This chapter focuses on using liquid chromatography-mass spectrometry (LC-MS)-based untargeted metabolomics for analysis of the exposome, including (1) preprocessing of untargeted metabolomics data, (2) identification of exposure chemicals and their metabolites, and (3) methods to establish associations between exposures and diseases.

Urine Metabolic Signatures of Multiple Environmental Pollutants in Pregnant Women: An Exposome Approach

Exposure to environmental pollutants, particularly during pregnancy, can have adverse consequences on child development but little is known about the effects of pollutant mixtures on endogenous metabolism in pregnant women. We aimed to identify urinary metabolic signatures associated with low level exposure to multiple environmental pollutants in pregnant women from the INMA (INfancia y Medio Ambiente) birth cohort (Spain, N = 750). 35 chemical exposures were quantified in first trimester blood samples (organochlorine pesticides, PCBs, PFAS), in cord blood (mercury), and twice in urine at 12 and 32 weeks of pregnancy (metals, phthalates, bisphenol A). ¹H nuclear magnetic resonance (NMR) metabolic profiles of urine were acquired in the same samples as pollutants. We explored associations between exposures and metabolism through an exposome-metabolome wide association scan and multivariate O2PLS modeling. Novel and reproducible associations were found across two periods of pregnancy for three nonpersistent pollutants and across two subcohorts for four of the persistent pollutants. We found novel metabolic signatures associated with arsenic exposure: TMAO and dimethylamine possibly related to gut microbial methylamine metabolism and homarine related to fish intake. Tobacco smoke exposure was related to coffee metabolism and PCBs with 3-hydroxyvaleric acid, usually released under ketoacidosis. These findings will have implications for further understanding of maternal-fetal health, and health across the life-course.