Prenatal Lead Exposure is Negatively Associated with the Gut Microbiome in Childhood
====================================================================================

* Shoshannah Eggers
* Vishal Midya
* Moira Bixby
* Chris Gennings
* Libni A Torres-Olascoaga
* Ryan W. Walker
* Robert O. Wright
* Manish Arora
* Martha María Téllez-Rojo

## Abstract

**Background** Metal exposures are associated with gut microbiome (GM) composition and function, and exposures early in development may be particularly important. Considering the role of the GM in association with many adverse health outcomes, understanding the relationship between prenatal metal exposures and the GM is critically important. However, there is sparse knowledge of the association between prenatal metal exposure and GM later in childhood.

**Objectives** This analysis aims to identify associations between prenatal lead (Pb) exposure and GM composition and function in children 9-11 years old.

**Methods** Data come from the Programming Research in Obesity, Growth, Environment and Social Stressors (PROGRESS) cohort based in Mexico City, Mexico. Prenatal metal concentrations were measured in maternal whole blood drawn during the second and third trimesters of pregnancy. Stool samples collected at 9-11 years old underwent metagenomic sequencing to assess the GM. This analysis uses multiple statistical modeling approaches, including linear regression, permutational analysis of variance, weighted quantile sum regression (WQS), and individual taxa regressions, to estimate the association between maternal blood Pb during pregnancy and multiple aspects of the child GM at 9-11 years old, adjusting for relevant confounders.

**Results** Of the 123 child participants in this pilot data analysis, 74 were male and 49 were female. Mean prenatal maternal blood Pb was 33.6(SE=2.1) ug/L and 34.9(SE=2.1) ug/L at second and third trimesters, respectively. Analysis suggests a consistent negative relationship between prenatal maternal blood Pb and the GM at age 9-11, including measures of alpha and beta diversity, microbiome mixture analysis, and individual taxa. The WQS analysis showed a negative association between prenatal Pb exposure and the gut microbiome, for both second and third trimester exposures (2Tβ=-0.17,95%CI=[-0.46,0.11]; 3Tβ=-0.17,95%CI=[-0.44,0.10]). Ruminococcus gnavus, Bifidobacterium longum, Alistipes indistinctus, Bacteroides caccae, and Bifidobacterium bifidum all had weights above the importance threshold from 80% or more of the WQS repeated holdouts in association with both second and third trimester Pb exposure.

**Discussion** Pilot data analysis suggests a negative association between prenatal Pb exposure and the gut microbiome later in childhood; however, additional investigation is needed.

Keywords
*   Human Microbiome
*   Weighted Quantile Sum Regression
*   Metal Exposure
*   Environmental Epidemiology
*   Maternal Child Health

## Introduction

Lead (Pb) has been a recognized environmental hazard for centuries,(Woolley, 1984) however its etiological pathways to disease are still not entirely understood. One potential mechanistic pathway between Pb exposure and its many downstream adverse health outcomes may be via the human gut microbiome. The collection of trillions of microbes that inhabit the human gut, including bacteria, fungi, viruses, and archaea, as well as their many genetic functions, are known as the gut microbiome.(Human Microbiome Project Consortium, 2012) The normal functions of the gut microbiome include nutrient metabolism, support of the mucosal and epithelial barriers within the gut, and interactions with immune, nervous, and cardiovascular systems.(Gomaa, 2020; Wilmes et al., 2022) Xenobiotic exposures, including Pb, influence the composition and function of the gut microbiome, affecting its interaction with systemic bodily function, and may lead to altered health status.(Claus et al., 2016)

Animal studies have linked Pb exposure to alterations of the gut microbiome including shifts in individual sample diversity (α-diversity), community composition (β-diversity), and individual bacterial taxa and gene abundance (Breton et al., 2013; Wu et al., 2016; Gao et al., 2017; Xia et al., 2018). Moreover, the body of epidemiologic evidence linking Pb exposure to altered gut microbiome composition and function is growing rapidly. Studies have linked human Pb exposure to increased abundance of many specific bacteria, including Proteobacteria (Bisanz et al., 2014; Eggers et al., 2019), a common indicator of gut microbiome dysbiosis, or imbalance (Shin et al., 2015; Litvak et al., 2017). Other analyses have linked Pb exposure to bacteria that are known to affect gut barrier integrity and gut health.(Laue et al., 2020; Sitarik et al., 2020; Shen et al., 2022) Alterations in gut and blood-brain barrier integrity resulting from metal-associated shifts in the gut microbiome may allow for increased metal toxicity by allowing for greater absorption into the bloodstream. While this growing field of epidemiologic research has uncovered relationships between Pb exposure and the gut microbiome, so far the analyses have been limited to the use of 16S rRNA amplicon sequencing data, which is less accurate in assigning taxonomy and inferring gene function than metagenomic sequencing. More studies with advanced omics analysis are needed to understand the relationship between Pb exposure and the gut microbiome.

Little is known about the relationship between prenatal exposures and the gut microbiome later in childhood. In fact, there are relatively few studies that investigate the gut microbiome in general between the ages of 5 and 20, especially from populations in low to middle-income countries.(Ortega, 2022) Using data and specimens from the Programming Research in Obesity, Growth, Environment and Social Stressors (PROGRESS) cohort, we examine relationships between prenatal Pb exposure and the gut microbiome of 9-11 year old children from Mexico City, Mexico. In this study, we aim to identify critical windows of prenatal Pb exposure that are associated with the gut microbiome later in childhood. Given the existing evidence from epidemiological and animal studies, we hypothesize that elevated prenatal Pb exposure during at least one time point will be associated with altered gut microbiome composition at 9-11 years old.

## Methods

### Study Design

PROGRESS is an ongoing prospective birth cohort in Mexico City, Mexico. The study enrolled 948 women in early pregnancy who went on to a live birth through the Mexican Social Security System, followed the offspring in infancy every 6 months, and biannually thereafter. The focus of PROGRESS is on neurobehavioral development and child obesity, with emphasis on environmental exposures, like metals, in pregnancy and early life that program later life behavior and growth. Exposures and outcomes were assessed at several time points beginning in the 2nd trimester of pregnancy through a combination of surveys, physical exams, psychological and behavioral assessments. At each visit, biological specimens (including blood) were collected, processed, aliquoted and stored. Stool samples were collected from a subset of participants (n=123) when the children were between the ages of 9-11. Protocols for the main PROGRESS study, and its ancillary microbiome study were reviewed and approved by the Institutional Review Board at the Icahn School of Medicine at Mount Sinai (STUDY-12-00751A, STUDY-21-00242).

### Pb Exposure Measurement

Prenatal Pb exposure in the PROGRESS cohort has been previously described. (Heiss et al., 2020) Briefly, Pb exposure was assessed using maternal whole blood, drawn during the second and third trimesters of pregnancy, at an average of 18.3 and 31.6 weeks gestation, respectively. Pb level was measured using inductively coupled plasma mass spectrometry (ICP-MS) in the trace metals laboratory at the Icahn School of Medicine at Mount Sinai.

### Gut Microbiome Sample Collection and Processing

Participants were recruited at the PROGRESS clinic visit as part of the 9-11 year visit. Samples were collected in clinic or at home using a sampling kit provided. Once obtained, stool samples were stored in a biosafety bag in the participant’s home refrigerator at 4°C. The sample was retrieved by a driver within 24 hours, processed following the FAST(Romano et al., 2018) protocol, and stored at -70°C within 48 hours from sample deposit. Frozen samples were sent to the Microbiome Translational Center at Mount Sinai. Samples were then processed and sequenced in two batches, with n=50 and n=73 samples, respectively. Whole genome sequencing was performed using the NEBNext DNA Library Prep kit and sequenced on an Illumina HiSeq. Sequencing reads were quality trimmed with Trimmomatic (Bolger et al., 2014) and human reads removed by mapping to a reference with bowtie2. (Langmead and Salzberg, 2012) The remaining reads were processed using MetaPhlAn2 (Truong et al., 2015) and StrainPhlAn (Truong et al., 2017) to determine microbial taxonomy down to the species/strain level, and HUMAnN2 (Franzosa et al., 2018) to profile microbial gene pathways.

### Covariates

Covariates used in this analysis included child sex, child’s age at time of stool sample collection, mother’s socio-economic status (SES) during pregnancy, mother’s age at birth, mother’s body mass index (BMI) during pregnancy, and microbiome analysis batch. Mother’s height and weight were collected with a professional digital scale and stadiometer at the study visit during the second trimester of pregnancy. Weight and height were used to calculate BMI, which was modeled as a continuous covariate in regression analyses. SES during pregnancy was assessed based on the 1994 Mexican Association of Intelligence Agencies Market and Opinion (AMAI) rule 13*6. Families are classified in six levels based on 13 questions about characteristics of the household. Most families in the study were low to middle SES, thus the six categories were condensed into three: lower, middle, and higher.(Sanders et al., 2022)

### Statistical Analysis

All analysis was conducted in R version 4.0.3. Any two tailed p-value less than 0.05 was considered statistically significant.

#### Data Processing

The count data for each Taxa were converted into relative abundance. To consider any possible batch effect while measuring taxa count, (1) only those taxa with at least 5% relative abundance in both batches were further considered in primary statistical analyses; (2) All models were controlled for a batch indicator variable along with other covariates. (3) Further detailed sensitivity analysis was conducted by choosing taxa with at least 25% relative abundance in both batches. The second and third-trimester Pb concentrations were log (base = 2) transformed to meet distributional assumptions with higher confidence.

#### Alpha & Beta Diversity

We calculated Shannon alpha diversity and Bray-Curtis beta diversity.(Shannon, 1948; Bray and Curtis, 1957) To estimate the associations with second and third-trimester lead concentrations, (1) for alpha diversity, we used the Kruskal-Wallis rank sum test (without covariate adjustment) and linear regression (with covariate adjustment); and (2) for beta diversity, we used PERMANOVA with 10000 permutations and w/o covariate adjustments.(Oksanen et al., 2019)

#### Weighted Quantile Sum Regression

We used the Weighted Quantile Sum (WQSRSRH) regression (Carrico et al., 2015) with random subset (Curtin et al., 2019) and repeated holdouts (Tanner et al., 2019), as established for microbiome data (Eggers et al., 2023), to estimate the association between second and third-trimester lead exposures and the effect on the abundance of the overall mixture of microbial taxa. Since the interest lies in the directionality of the association, the WQSRSRH model was fitted with all the chosen taxa as exposures and the log transformed Pb concentration as the outcome. For ease of interpretation, the relative abundance of the taxa was converted into deciles, while a null relative abundance was kept at zero. Moreover, two WQSRSRH models were fitted at each trimester, with the overall mixture effect assumed in the negative or positive directions. The final optimal model (Liao et al., 2018) from each trimester was chosen based on the smallest Akaike information criterion (AIC) (Akaike, 1974) and the Bayesian information criterion (BIC) (Schwarz, 1978). Those microbial taxa were judged important in the final chosen models, which had weight contributions to the overall mixture index above a chance threshold (1/the number of components in the index).

Lastly, to account for any between taxa correlations and the relatively smaller sample size, each WQSRSRH model was fitted based on 200 repeated holdouts (with randomly 40% data set aside for validation) and 100 bootstrapped with-replacement sampling at each iteration.

#### Taxa-Wide Association Analysis

To estimate the effect of Pb exposure on individual relative abundance of each chosen taxa, we conducted Taxa-wide association analysis (TWAS) with generalized linear models at both trimesters. Further, the unadjusted raw p-values with respect to regression beta estimates were plotted through Volcano plots. The Bonferroni procedure was adapted for multiple comparison error correction on raw p-values. We estimated the effective number of tests for the TWAS using the eigenvalues of the relative abundance of the correlation matrix (Li and Ji, 2005; de Prado-Bert et al., 2021).

#### Gene Function Analysis

We extracted the microbial gene function pathways of the important taxa from the WQSRSRH analysis at each trimester. We selected the top 20 most frequently occurring pathways in each trimester for ease of interpretation. Through a simple Venn diagram, we further elaborated on the pathways which are common to a) both trimesters, b) only present in the second trimester, and lastly, c) only present in the third trimester.

#### Sensitivity Analysis

We also conducted a sensitivity analysis to understand whether the choice of relative abundance in both batches affected the overall mixture effect. To that end, we repeated the entire WQSRSRH analysis with only those taxa having at least 25% relative abundance in both batches. Further, we chose not to rescale the chosen taxa with at least 5% relative abundance (1) to reflect the original contributions and (2) to make the analysis robust irrespective of the chosen relative abundance cutoff.

#### Covariate Adjustment

Each model was controlled by *a priori* chosen set of covariates. Although we thought Pb exposure during childhood may be a confounder, we tested the correlation between prenatal Pb and childhood Pb at birth, 1 year, 2 years, and 4 years of age, and did not find a correlation, thus we decided to exclude these variables from the models. A few covariates had missing values (less than 5%), which were imputed by the multiple imputation chained equations as implemented in the “MICE” R package (Buuren and Groothuis-Oudshoorn, 2011).

## Results

### Study Population Characteristics

Of 123 participants in this study (Table 1), 49 were female, and 74 were male. Mean Pb concentration was 33.6 ug/L and 34.9 ug/L in the second and third trimesters of pregnancy, respectively. Mothers with low SES were more likely to be in the fourth quartile of Pb exposure for both trimesters of pregnancy.

View this table:
[Table 1.](http://medrxiv.org/content/early/2023/05/11/2023.05.10.23289802/T1)

Table 1. 
Characteristics of the analytical study population by Pb exposure in the second and third trimester of pregnancy.

### Alpha & Beta Diversity

In linear regression analysis of alpha (within individual) diversity and prenatal Pb exposure, we found slight negative associations, that were not statistically significant, for both second and third trimester Pb exposures, in unadjusted and adjusted models (Table 2). In PERMANOVA analysis of beta (between individual) diversity (Table 2), second trimester Pb exposure was associated with a non-significant R2 of less than 1% (Adjusted R2 = 0.007, p-value = 0.515). However, third trimester Pb exposure was associated with an R2 of 1.1% in the adjusted model, and was directionally trending (Adjusted R2 = 0.011, p-value = 0.066).

View this table:
[Table 2:](http://medrxiv.org/content/early/2023/05/11/2023.05.10.23289802/T2)

Table 2: 
Estimates of association between alpha and beta diversity and prenatal Pb exposure.

### Microbiome Mixture Analysis

The primary WQSRSRH analysis was run in the negative direction because we hypothesized the association between prenatal Pb exposure and the gut microbiome mixture to be negative (Figure 1). Including adjustment for covariates, second trimester Pb exposure was negatively associated with the gut microbiome mixture (β = -0.17, 95%CI = [-0.46, 0.11]), with 88% of the repeated holdout estimates below zero. Third trimester Pb exposure showed a very similar association with the gut microbiome mixture (β = -0.17, 95%CI = [-0.44, 0.10]), and had 89% of the repeated holdout estimates below zero. Within the weighted indices, taxa with a weight above 0.027 were considered important in the mixture association. Of the 20 bacterial taxa above the importance threshold for second trimester Pb exposure, 16 were also above the importance threshold in association with third trimester Pb exposure. Ruminococcus gnavus, Bifidobacterium longum, Alistipes indistinctus, Bacteroides caccae, and *Bifidobacterium bifidum* all had weights above the importance threshold from 80% or more of the repeated holdouts in association with both second and third trimester Pb exposure.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/05/11/2023.05.10.23289802/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2023/05/11/2023.05.10.23289802/F1)

Figure 1. 
WQSRSRH estimates for the association of the gut microbiome mixture with prenatal Pb exposure in the a) second and c) third trimester of pregnancy. Small grey jitter dots represent the WQSRS estimates from each repeated holdout. The violin plot around those points show the density of holdout estimates at each interval. Avserage percent weight for each taxa within the WQS index are shown for the b) second and d) third trimester Pb exposure. The green line indicates the importance threshold for weights above random chance. Taxa that are labeled in blue are above the threshold for both trimesters of exposure.

In a sensitivity analysis, we ran the same WQS model using only the bacterial taxa that were present in at least 25% (instead of 5%) of participants from both analytical batches, and found associations in the same direction, with slightly larger estimates and confidence intervals that still crossed zero (Sup. Figure 1). We conducted an additional sensitivity analysis running the WQSRSRH analysis in the positive direction and found null results. When comparing the likelihood estimates between the negative and positive WQSRSRH models, the likelihood was higher for the negative model, confirming our appropriate use of the negative model as the primary analysis.

When examining microbial gene function pathways of the important taxa associated with each of the second and third trimester Pb exposure separately, of the top 20 most abundant gene pathways for each trimester, approximately 1/2 of pathways from each trimester were found to be unique (Figure 2). Overall, common pathways were more likely to do with nucleic acid biosynthesis, and functions essential to all bacteria, while pathways associated with only one trimester Pb exposure or the other, were more likely to be involved in amino acid biosynthesis, and more specialized metabolic functions.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/05/11/2023.05.10.23289802/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2023/05/11/2023.05.10.23289802/F2)

Figure 2. 
Venn diagram of the top 20 most abundant microbial gene functions from the highly weighted taxa in the WQSRSRH analysis.

### Individual Taxa Analysis

In analysis of each individual bacterial taxa with prenatal Pb exposure (Figure 3), we found six taxa (*Alistipes putredinis, Ruminococcus ghavus, Bacteroides caccae, B. intestinalis, Coprococcus catus*, and *A. indistinctus*) to be negatively associated, and one (*B. coprocola*) to be positively associated with second trimester exposure. With third trimester exposure, we found three taxa (*Bifidobacterium bifidum, B. longum, A. indistinctus*) to be negatively associated, and three taxa (*B. coprocola, Eubacterium eligens, B. finegoldii*) to be positively associated with Pb.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/05/11/2023.05.10.23289802/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2023/05/11/2023.05.10.23289802/F3)

Figure 3. 
Volcano plot of estimates and p-values from the taxa wide association analysis (TWAS) of bacterial abundance in association with a) second and b) third trimester Pb exposure. Blue dots are taxa associated in the negative direction and black dots are associated in the positive direction.

## Discussion

In this analysis of pilot microbiome data from the PROGRESS cohort, we examined associations between prenatal Pb exposure and several different components of gut microbiome composition and function. We found consistent negative associations between Pb exposure during the second and third trimester of pregnancy and several assessments of gut microbiome composition and function at 9-11 years old. Associations between prenatal Pb and the gut microbiome later in childhood tended to be strongest for second trimester exposure, providing support for the second trimester of pregnancy in particular as a critical window of exposure.

The focus of this analysis was on the association between prenatal Pb exposure and the gut microbiome later in childhood. There are many other potential exposures that may influence the composition and function of the gut microbiome between the prenatal exposure and the gut microbiome at 9-11 years, including continued exposure to Pb, diet, antibiotic use, animal exposure, child BMI and many others. However, none of these factors were included in this analysis because they are not confounders, as they are either 1) not associated with prenatal Pb exposure, or 2) potentially on the causal pathway between prenatal Pb exposure and the gut microbiome later in life.

The underlying mechanism of association between prenatal Pb exposure and the gut microbiome later in childhood could work in multiple ways. Because the gut microbiome is shaped in part by the host immune system, (Petersen et al., 2019; Wong et al., 2022) prenatal Pb exposure may alter the trajectory of immune system development, which can influence gut microbial composition through the life-course. Another possible mechanism is via the maternal gut microbiome. The infant gut microbiome is primarily colonized by bacteria that are transferred from the mother at birth. (Ferretti et al., 2018) Pb exposure in adults also influences gut microbiome composition,(Eggers et al., 2019) thus Pb induced changes in the maternal gut microbiome during pregnancy may be transferred to children at birth, and carried on into childhood, and even later in life.(Gschwendtner et al., 2019) Another potential mechanism may be that maternal Pb exposure is transferred to children through breastmilk, an established Pb exposure pathway (Klein et al., 2017). In that case, postnatal Pb exposure through breastmilk would be along the causal pathway between prenatal Pb exposure and the gut microbiome. Furthermore, all of these mechanisms may be working together. Further mechanistic studies are needed to better understand this complex relationship.

In previously published analyses of prenatal Pb exposure and the gut microbiome in studies with relatively small sample sizes, other researchers have seen similar results. In analysis of a birth cohort based in Detroit, MI (n=146), Sitarik et al, identified associations between prenatal Pb exposure, measured in baby teeth, and decreased abundance of several species of *Bacteroides* within the gut microbiome at 1 month and 6 months of age.(Sitarik et al., 2020) They identified *B. caccae* as being negatively associated with second trimester Pb exposure specifically. We also identified a decreased abundance of *B. caccae* in association with Pb exposure in the second trimester using TWAS, and *B. caccae* was also heavily weighted in WQSRSRH for both trimesters of prenatal Pb exposure. *B. caccae* are common fiber degraders in the human gut microbiome, with mixed health effects. (Wexler, 2007; Chen et al., 2021; Yang et al., 2022; Zhang et al., 2022) This consistent association of prenatal Pb exposure and decreased *B. caccae* abundance both early and late in childhood, from different populations in different countries, is strong evidence of association in this relatively new field.

In another epidemiologic study using data from a Canadian cohort (n=70), Shen et al, found that prenatal Pb exposure, measured in maternal blood, was associated with increased abundance of *Fusobacteriota* in the gut microbiome at 6-7 years of age.(Shen et al., 2022) They did not find a significant association between prenatal Pb and alpha or beta diversity. While the associated gut microbes identified in this analysis and the study by Shen et al, were not the same, it is important to note that links between prenatal Pb and the gut microbiome many years later in childhood were identified, even in these two small data sets. This suggests that there are likely true underlying links between prenatal Pb exposure and the gut microbiome in childhood. Differences in bacterial taxa associated may be due to differences in the study participants, or differences in sample and data processing procedures.

Of all the taxa in this analysis, *Alistipes indistinctus* was the only one identified by both WQSRSRH and TWAS to have a negative association with Pb exposure in the second and third trimesters. *A. indistinctus* is a common member of the human gut microbiome, and relatively newly identified.(Nagai et al., 2010) Little is known about the health implications of *A. indistinctus*. However, one study has identified it as protective against liver fibrosis.(Shao et al., 2018), and another identified it as a keystone species for restoring a healthy gut microbiome in patients with non-alcoholic fatty liver disease.(Wu et al., 2022) Overall, the *Alistipes* genus has been shown to have both beneficial and detrimental health effects in humans.(Parker et al., 2020) Continued investigation of the health implications of *A. indistinctus* are necessary to understand the links between prenatal Pb exposure and downstream health status via the microbiome.

In the analysis of microbial gene pathways from the most highly weighted taxa in the WQSRSRH analysis, approximately 1/3 of the gene pathways were shared between the taxa associated with second and third trimester Pb exposure. With a few exceptions, the pathways that were common among these bacteria are used for nucleic acid biosynthesis and other essential functions for bacterial life.(Tsuchiya et al., 2018) The pathways that were not shared were more likely to be used in amino acid biosynthesis, fermentation, and other metabolic pathways. Because these more varied pathways provide a wider range of functional capabilities, they are more likely to influence host health, although their direct influence is not yet fully understood.

While this study added to the growing field of evidence around the negative relationship between metal exposure and the human gut microbiome, there are some limitations to consider. Because this was a pilot study, sample size was limited, which limited our power to detect associations. However, because results were primarily in the negative direction across multiple analytical approaches, and over 85% of the beta estimates from the WQS repeated holdouts were negative, we have increased confidence in reporting a negative result. Samples from this pilot study were also analyzed in two batches. Our strategy for reducing batch effects limited the breadth of microbiome data we could include in the analysis; however, similar prevalence thresholds are frequently used as data reduction steps in microbiome data analysis. Another statistical limitation was the use of prenatal Pb exposure as an outcome in the WQSRSRH models rather than a predictor. This model structure was necessary due to WQS format. The implications of this limitation are minimal however, as this analysis was used to determine association, not causation. Lastly, the use of maternal blood Pb during pregnancy to measure prenatal Pb exposure is not ideal, as it is not a direct measure of fetal Pb exposure.

In future analyses based off this work, we hope to examine additional and more nuanced relationships between prenatal environmental metal exposures and the gut microbiome in childhood. We plan to do additional sample collection in this cohort that will expand the sample size for future analyses. Furthermore, we are developing additional methodological approaches to examine relationships between chemical exposures and the gut microbiome, and their combined effect on downstream health.

In conclusion, this pilot study found a consistent negative association between prenatal Pb exposure and the gut microbiome in late childhood. These results support a growing body of evidence that human Pb exposure may alter gut microbial composition and function, leading to downstream health implications. More studies with larger sample sizes are needed to better understand this relationship.

## Supporting information

Supplemental Information [[supplements/289802_file05.docx]](pending:yes)

## Data Availability

The data used in this study can be made accessible to researchers upon appropriate request with the following restrictions to ensure the privacy of human subjects. Note that access to the data is limited due to a data-sharing agreement approved by the IRB at Mount Sinai. Researchers interested in accessing PROGRESS data must send their resume/CV and CITI training certificates to the IRB chair, Ilene Wilets (Ilene.wilets{at}mssm.edu). They must also submit a data analysis plan to the Principal Investigators for PROGRESS; Robert O. Wright (robert.wright{at}mssm.edu), Martha Tellez-Rojo (mmtellez{at}insp.mx), and Andrea Baccarelli (andrea.baccarelli{at}columbia.edu). Once this process is completed, the PROGRESS data analyst, Nia McRae (nia.x.mcrae{at}mssm.edu), will send a de-identified dataset via Box, a secure data-sharing platform. 

## 3 Conflict of Interest

MA is an employee and equity holder of Linus Biotechnology Inc., a start-up company of Mount Sinai Health System. The company develops tools for the detection of autism spectrum disorder and related conditions. The following authors report no competing interests: SE, VM, CG, MB, LT, RWW, ROW, MT.

## 4 Author Contributions

SE, VM, CG and MA contributed to conception and design of the study. LT, MT, and ROW contributed to data acquisition. VM and MB performed statistical analysis. SE, VM, CG, and RWW contributed to data interpretation. SE wrote the first draft of the manuscript. VM wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.

## 5 Funding

This work was supported by the National Institute of Environmental Health Sciences (K99ES032884, P30ES023515, R01ES013744).

## 17 Data Availability Statement

The data that was used in this study can be made accessible to researchers upon appropriate request with the following restrictions to ensure the privacy of human subjects. Note that access to the data is limited due to a data sharing agreement approved by the IRB at Mount Sinai. Researchers that are interested in accessing PROGRESS data must send their resume/CV as well as CITI training certificates to the IRB chair, Ilene Wilets (Ilene.wilets{at}mssm.edu). They must also submit a data analysis plan to the Principal Investigators for PROGRESS; Robert O. Wright (robert.wright{at}mssm.edu), Martha Tellez-Rojo (mmtellez{at}insp.mx), and Andrea Baccarelli (andrea.baccarelli{at}columbia.edu). Once this process is completed, the PROGRESS data analyst, Nia McRae (nia.x.mcrae{at}mssm.edu) will send a de-identified dataset via Box, a secure data sharing platform.

## 6 Acknowledgments

The authors would like to acknowledge the entire PROGRESS study team, as well as the participants. We would also like to thank Dr. Jeremiah Faith and the Microbiome Translational Center at the Icahn School of Medicine at Mount Sinai.

*   Received May 10, 2023.
*   Revision received May 10, 2023.
*   Accepted May 11, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  Akaike, H. (1974). A new look at the statistical model identification. IEEE Trans. Autom. Control 19, 716–723.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TAC.1974.1100705&link_type=DOI) 

2.  Bisanz, J. E., Enos, M. K., Mwanga, J. R., Changalucha, J., Burton, J. P., Gloor, G. B., et al. (2014). Randomized open-label pilot study of the influence of probiotics and the gut microbiome on toxic metal levels in Tanzanian pregnant women and school children. mBio 5, e01580–01514. doi: 10.1128/mBio.01580-14.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/mBio.01580-14&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25293764&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 

3.  Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinforma. Oxf. Engl. 30, 2114–2120. doi: 10.1093/bioinformatics/btu170.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu170&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24695404&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000340049100004&link_type=ISI) 

4.  Bray, J. R., and Curtis, J. T. (1957). An Ordination of the Upland Forest Communities of Southern Wisconsin. Ecol. Monogr. 27, 325–349. doi: 10.2307/1942268.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/1942268&link_type=DOI) 

5.  Breton, J., Massart, S., Vandamme, P., De Brandt, E., Pot, B., and Foligné, B. (2013). Ecotoxicology inside the gut: impact of heavy metals on the mouse microbiome. BMC Pharmacol. Toxicol. 14, 62. doi: 10.1186/2050-6511-14-62.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/2050-6511-14-62&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24325943&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 

6.  Buuren, S. van, and Groothuis-Oudshoorn, K. (2011). mice: Multivariate Imputation by Chained Equations in R. J. Stat. Softw. 45, 1–67. doi: 10.18637/jss.v045.i03.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v045.i03&link_type=DOI) 

7.  Carrico, C., Gennings, C., Wheeler, D. C., and Factor-Litvak, P. (2015). Characterization of Weighted Quantile Sum Regression for Highly Correlated Data in a Risk Analysis Setting. J. Agric. Biol. Environ. Stat. 20, 100–120. doi: 10.1007/s13253-014-0180-3.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s13253-014-0180-3&link_type=DOI) 

8.  Chen, Y., Lin, H., Cole, M., Morris, A., Martinson, J., Mckay, H., et al. (2021). Signature changes in gut microbiome are associated with increased susceptibility to HIV-1 infection in MSM. Microbiome 9, 237. doi: 10.1186/s40168-021-01168-w.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s40168-021-01168-w&link_type=DOI) 

9.  Claus, S. P., Guillou, H., and Ellero-Simatos, S. (2016). The gut microbiota: a major player in the toxicity of environmental pollutants? NPJ Biofilms Microbiomes 2, 16003. doi: 10.1038/npjbiofilms.2016.3.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/npjbiofilms.2016.3&link_type=DOI) 

10. Curtin, P., Kellogg, J., Cech, N., and Gennings, C. (2019). A random subset implementation of weighted quantile sum (WQSRS) regression for analysis of high-dimensional mixtures. Commun. Stat. - Simul. Comput. , 1–16. doi: 10.1080/03610918.2019.1577971.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/03610918.2019.1577971&link_type=DOI) 

11. de Prado-Bert, P., Ruiz-Arenas, C., Vives-Usano, M., Andrusaityte, S., Cadiou, S., Carracedo, Á., et al. (2021). The early-life exposome and epigenetic age acceleration in children. Environ. Int. 155, 106683.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2021.106683&link_type=DOI) 

12. Eggers, S., Bixby, M., Renzetti, S., Curtin, P., and Gennings, C. (2023). Human Microbiome Mixture Analysis Using Weighted Quantile Sum Regression. Int. J. Environ. Res. Public. Health 20, 94. doi: 10.3390/ijerph20010094.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijerph20010094&link_type=DOI) 

13. Eggers, S., Safdar, N., Sethi, A. K., Suen, G., Peppard, P. E., Kates, A. E., et al. (2019). Urinary lead concentration and composition of the adult gut microbiota in a cross-sectional populationbased sample. Environ. Int. 133, 105122. doi: 10.1016/j.envint.2019.105122.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2019.105122&link_type=DOI) 

14. Ferretti, P., Pasolli, E., Tett, A., Asnicar, F., Gorfer, V., Fedi, S., et al. (2018). Mother-to-Infant Microbial Transmission from Different Body Sites Shapes the Developing Infant Gut Microbiome. Cell Host Microbe 24, 133-145.e5. doi: 10.1016/j.chom.2018.06.005.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chom.2018.06.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 

15. Franzosa, E. A., McIver, L. J., Rahnavard, G., Thompson, L. R., Schirmer, M., Weingart, G., et al. (2018). Species-level functional profiling of metagenomes and metatranscriptomes. Nat. Methods 15, 962–968. doi: 10.1038/s41592-018-0176-y.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41592-018-0176-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 

16. Gao, B., Chi, L., Mahbub, R., Bian, X., Tu, P., Ru, H., et al. (2017). Multi-Omics Reveals that Lead Exposure Disturbs Gut Microbiome Development, Key Metabolites, and Metabolic Pathways. Chem. Res. Toxicol. 30, 996–1005. doi: 10.1021/acs.chemrestox.6b00401.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/acs.chemrestox.6b00401&link_type=DOI) 

17. Gomaa, E. Z. (2020). Human gut microbiota/microbiome in health and diseases: a review. Antonie Van Leeuwenhoek 113, 2019–2040. doi: 10.1007/s10482-020-01474-7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10482-020-01474-7&link_type=DOI) 

18. Gschwendtner, S., Kang, H., Thiering, E., Kublik, S., Fösel, B., Schulz, H., et al. (2019). Early life determinants induce sustainable changes in the gut microbiome of six-year-old children. Sci. Rep. 9, 12675. doi: 10.1038/s41598-019-49160-7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-019-49160-7&link_type=DOI) 

19. Heiss, J. A., Téllez-Rojo, M. M., Estrada-Gutiérrez, G., Schnaas, L., Amarasiriwardena, C., Baccarelli, A. A., et al. (2020). Prenatal lead exposure and cord blood DNA methylation in PROGRESS: an epigenome-wide association study. Environ. Epigenetics 6, dvaa014. doi: 10.1093/eep/dvaa014.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/eep/dvaa014&link_type=DOI) 

20. Human Microbiome Project Consortium (2012). Structure, function and diversity of the healthy human microbiome. Nature 486, 207–214. doi: 10.1038/nature11234.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature11234&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22699609&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305189000025&link_type=ISI) 

21. Klein, L. D., Breakey, A. A., Scelza, B., Valeggia, C., Jasienska, G., and Hinde, K. (2017). Concentrations of trace elements in human milk: Comparisons among women in Argentina, Namibia, Poland, and the United States. PLoS ONE 12, e0183367. doi: 10.1371/journal.pone.0183367.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0183367&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28817665&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 

22. Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.1923&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22388286&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302218500017&link_type=ISI) 

23. Laue, H. E., Moroishi, Y., Jackson, B. P., Palys, T. J., Madan, J. C., and Karagas, M. R. (2020). Nutrient-toxic element mixtures and the early postnatal gut microbiome in a United States longitudinal birth cohort. Environ. Int. 138, 105613. doi: 10.1016/j.envint.2020.105613.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2020.105613&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 

24. Li, J., and Ji, L. (2005). Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix. Heredity 95, 221–227.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/sj.hdy.6800717&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16077740&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231438500007&link_type=ISI) 

25. Liao, J. G., Cavanaugh, J. E., and McMurry, T. L. (2018). Extending AIC to best subset regression. Comput. Stat. 33, 787–806. doi: 10.1007/s00180-018-0797-8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00180-018-0797-8&link_type=DOI) 

26. Litvak, Y., Byndloss, M. X., Tsolis, R. M., and Bäumler, A. J. (2017). Dysbiotic Proteobacteria expansion: a microbial signature of epithelial dysfunction. Curr. Opin. Microbiol. 39, 1–6. doi: 10.1016/j.mib.2017.07.003.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mib.2017.07.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28783509&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 

27. Nagai, F., Morotomi, M., Watanabe, Y., Sakon, H., and Tanaka, R. (2010). Alistipes indistinctus sp. nov. and Odoribacter laneus sp. nov., common members of the human intestinal microbiota isolated from faeces. Int. J. Syst. Evol. Microbiol. 60, 1296–1302. doi: 10.1099/ijs.0.014571-0.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1099/ijs.0.014571-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19667375&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000279369100007&link_type=ISI) 

28. Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2019). vegan: Community Ecology Package. Available at: [https://CRAN.R-project.org/package=vegan](https://CRAN.R-project.org/package=vegan) x[Accessed August 12, 2019].
    
    

29. Ortega, R. P. (2022). Microbiome data dominated by wealthy countries. Science 375, 709. doi: 10.1126/science.ada1336.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/science.ada1336&link_type=DOI) 

30. Parker, B. J., Wearsch, P. A., Veloo, A. C. M., and Rodriguez-Palacios, A. (2020). The Genus Alistipes: Gut Bacteria With Emerging Implications to Inflammation, Cancer, and Mental Health. Front. Immunol. 11. Available at: [https://www.frontiersin.org/articles/10.3389/fimmu.2020.00906](https://www.frontiersin.org/articles/10.3389/fimmu.2020.00906) x[Accessed February 23, 2023].
    
    

31. Petersen, C., Bell, R., Klag, K. A., Lee, S.-H., Soto, R., Ghazaryan, A., et al. (2019). T cell–mediated regulation of the microbiota protects against obesity. Science 365, eaat9351. doi: 10.1126/science.aat9351.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiIzNjUvNjQ1MS9lYWF0OTM1MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA1LzExLzIwMjMuMDUuMTAuMjMyODk4MDIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

32. Romano, K. A., Dill-McFarland, K. A., Kasahara, K., Kerby, R. L., Vivas, E. I., Amador-Noguez, D., et al. (2018). Fecal Aliquot Straw Technique (FAST) allows for easy and reproducible subsampling: assessing interpersonal variation in trimethylamine-N-oxide (TMAO) accumulation. Microbiome 6, 91. doi: 10.1186/s40168-018-0458-8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s40168-018-0458-8&link_type=DOI) 

33. Sanders, A. P., Gennings, C., Tamayo-Ortiz, M., Mistry, S., Pantic, I., Martinez, M., et al. (2022). Prenatal and early childhood critical windows for the association of nephrotoxic metal and metalloid mixtures with kidney function. Environ. Int. 166, 107361. doi: 10.1016/j.envint.2022.107361.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2022.107361&link_type=DOI) 

34. Schwarz, G. (1978). Estimating the dimension of a model. Ann. Stat., 461–464.
    
    

35. Shannon, C. E. (1948). A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1145/584091.584093&link_type=DOI) 

36. Shao, L., Ling, Z., Chen, D., Liu, Y., Yang, F., and Li, L. (2018). Disorganized gut microbiome contributed to liver cirrhosis progression: a meta-omics-based study. Front. Microbiol. 9, 3166.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fmicb.2018.03166&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 

37. Shen, Y., Laue, H. E., Shrubsole, M. J., Wu, H., Bloomquist, T. R., Larouche, A., et al. (2022). Associations of Childhood and Perinatal Blood Metals with Children’s Gut Microbiomes in a Canadian Gestation Cohort. Environ. Health Perspect. 130, 17007. doi: 10.1289/EHP9674.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1289/EHP9674&link_type=DOI) 

38. Shin, N.-R., Whon, T. W., and Bae, J.-W. (2015). Proteobacteria: microbial signature of dysbiosis in gut microbiota. Trends Biotechnol. 33, 496–503. doi: 10.1016/j.tibtech.2015.06.011.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tibtech.2015.06.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26210164&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 

39. Sitarik, A. R., Arora, M., Austin, C., Bielak, L. F., Eggers, S., Johnson, C. C., et al. (2020). Fetal and early postnatal lead exposure measured in teeth associates with infant gut microbiota. Environ. Int. 144, 106062. doi: 10.1016/j.envint.2020.106062.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2020.106062&link_type=DOI) 

40. Tanner, E. M., Bornehag, C.-G., and Gennings, C. (2019). Repeated holdout validation for weighted quantile sum regression. MethodsX 6, 2855–2860. doi: 10.1016/j.mex.2019.11.008.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mex.2019.11.008&link_type=DOI) 

41. Truong, D. T., Franzosa, E. A., Tickle, T. L., Scholz, M., Weingart, G., Pasolli, E., et al. (2015). MetaPhlAn2 for enhanced metagenomic taxonomic profiling. Nat. Methods 12, 902–903. doi: 10.1038/nmeth.3589.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.3589&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26418763&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 

42. Truong, D. T., Tett, A., Pasolli, E., Huttenhower, C., and Segata, N. (2017). Microbial strain-level population structure and genetic diversity from metagenomes. Genome Res. 27, 626–638. doi: 10.1101/gr.216242.116.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjg6IjI3LzQvNjI2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDUvMTEvMjAyMy4wNS4xMC4yMzI4OTgwMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

43. Tsuchiya, Y., Zhyvoloup, A., Bakovic, J., Thomas, N., Yu, B. Y. K., Das, S., et al. (2018). Protein CoAlation and antioxidant function of coenzyme A in prokaryotic cells. Biochem. J. 475, 1909–1937. doi: 10.1042/BCJ20180043.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6InBwYmlvY2hlbWoiO3M6NToicmVzaWQiO3M6MTE6IjQ3NS8xMS8xOTA5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDUvMTEvMjAyMy4wNS4xMC4yMzI4OTgwMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

44. Wexler, H. M. (2007). Bacteroides: the Good, the Bad, and the Nitty-Gritty. Clin. Microbiol. Rev. 20, 593–621. doi: 10.1128/CMR.00008-07.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiY21yIjtzOjU6InJlc2lkIjtzOjg6IjIwLzQvNTkzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDUvMTEvMjAyMy4wNS4xMC4yMzI4OTgwMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

45. Wilmes, P., Martin-Gallausiaux, C., Ostaszewski, M., Aho, V. T. E., Novikova, P. V., Laczny, C. C., et al. (2022). The gut microbiome molecular complex in human health and disease. Cell Host Microbe 30, 1201–1206. doi: 10.1016/j.chom.2022.08.016.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chom.2022.08.016&link_type=DOI) 

46. Wong, C. K., Yusta, B., Koehler, J. A., Baggio, L. L., McLean, B. A., Matthews, D., et al. (2022). Divergent roles for the gut intraepithelial lymphocyte GLP-1R in control of metabolism, microbiota, and T cell-induced inflammation. Cell Metab. 34, 1514-1531.e7. doi: 10.1016/j.cmet.2022.08.003.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cmet.2022.08.003&link_type=DOI) 

47. Woolley, D. E. (1984). A perspective of lead poisoning in antiquity and the present. Neurotoxicology 5, 353–361.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=6395049&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 

48. Wu, D., Liu, L., Jiao, N., Zhang, Y., Yang, L., Tian, C., et al. (2022). Targeting keystone species helps restore the dysbiosis of butyrate-producing bacteria in nonalcoholic fatty liver disease. iMeta 1, e61. doi: 10.1002/imt2.61.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/imt2.61&link_type=DOI) 

49. Wu, J., Wen, X. W., Faulk, C., Boehnke, K., Zhang, H., Dolinoy, D. C., et al. (2016). Perinatal Lead Exposure Alters Gut Microbiota Composition and Results in Sex-specific Bodyweight Increases in Adult Mice. Toxicol. Sci. Off. J. Soc. Toxicol. 151, 324–333. doi: 10.1093/toxsci/kfw046.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/toxsci/kfw046&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26962054&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F11%2F2023.05.10.23289802.atom) 

50. Xia, J., Lu, L., Jin, C., Wang, S., Zhou, J., Ni, Y., et al. (2018). Effects of short term lead exposure on gut microbiota and hepatic metabolism in adult zebrafish. Comp. Biochem. Physiol. Toxicol. Pharmacol. CBP 209, 1–8. doi: 10.1016/j.cbpc.2018.03.007.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cbpc.2018.03.007&link_type=DOI) 

51. Yang, X., Xiu, W.-B., Wang, J.-X., Li, L.-P., He, C., and Gao, C.-P. (2022). CO2 Is Beneficial to Gut Microbiota Homeostasis during Colonoscopy: Randomized Controlled Trial. J. Clin. Med. 11, 5281. doi: 10.3390/jcm11185281.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/jcm11185281&link_type=DOI) 

52. Zhang, Y., Fan, Q., Hou, Y., Zhang, X., Yin, Z., Cai, X., et al. (2022). Bacteroides species differentially modulate depression-like behavior via gut-brain metabolic signaling. Brain. Behav. Immun. 102, 11–22. doi: 10.1016/j.bbi.2022.02.007.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbi.2022.02.007&link_type=DOI)