Abstract
Many analytical methods used in gut microbiome research focus on either single bacterial taxa or the whole microbiome, ignoring multibacteria relationships (microbial cliques). We present a novel analytical approach to identify microbial cliques within the gut microbiome of children at 9-11 years associated with prenatal lead (Pb) exposure. Data came from a subset of participants (n = 123) in the Programming Research in Obesity, Growth, Environment and Social Stressors cohort. Pb concentrations were measured in maternal whole blood from the second and third trimesters of pregnancy. Stool samples collected at 9-11 years old underwent metagenomic sequencing to assess the gut microbiome. Using a novel analytical approach, Microbial Co-occurrence Analysis (MiCA), we paired a machine learning algorithm with randomization-based inference to first identify microbial cliques that were predictive of prenatal Pb exposure and then estimate the association between prenatal Pb exposure and microbial clique abundance. With second-trimester Pb exposure, we identified a two-taxa microbial clique that included Bifidobacterium adolescentis and Ruminococcus callidus and a three-taxa clique that also included Prevotella clara. Increasing second-trimester Pb exposure was associated with significantly increased odds of having the two-taxa microbial clique below the median relative abundance (odds ratio (OR) = 1.03, 95% confidence interval (CI) [1.01-1.05]). Using a novel combination of machine learning and causal inference, MiCA identified a significant association between second-trimester Pb exposure and the reduced abundance of a probiotic microbial clique within the gut microbiome in late childhood.