The normalized abundance mOTU table (metagenomics) and normalized peak intensities table (metabolomics/exposomics) were integrated using Data Integration Analysis for Biomarker discovery using latent components (DIABLO) within the mixOmics package in R. Details about this analysis are available in the GitLab repository . Spearman's correlation between the mOTUs and chemicals was computed in R, using the cor.test function of the base R stats package, while the Corrplot package, was used to perform the correlation plots.