Generally-healthy individuals with aberrant bowel movement frequencies show enrichment for microbially-derived blood metabolites associated with reduced kidney function - PubMed
- ️Mon Jan 01 2024
Generally-healthy individuals with aberrant bowel movement frequencies show enrichment for microbially-derived blood metabolites associated with reduced kidney function
Johannes P Johnson-Martínez et al. bioRxiv. 2024.
Update in
-
Johnson-Martínez JP, Diener C, Levine AE, Wilmanski T, Suskind DL, Ralevski A, Hadlock J, Magis AT, Hood L, Rappaport N, Gibbons SM. Johnson-Martínez JP, et al. Cell Rep Med. 2024 Jul 16;5(7):101646. doi: 10.1016/j.xcrm.2024.101646. Cell Rep Med. 2024. PMID: 39019013 Free PMC article.
Abstract
Bowel movement frequency (BMF) has been linked to changes in the composition of the human gut microbiome and to many chronic conditions, like metabolic disorders, neurodegenerative diseases, chronic kidney disease (CKD), and other intestinal pathologies like irritable bowel syndrome and inflammatory bowel disease. Lower BMF (constipation) can lead to compromised intestinal barrier integrity and a switch from saccharolytic to proteolytic fermentation within the microbiota, giving rise to microbially-derived toxins that may make their way into circulation and cause damage to organ systems. However, the connections between BMF, gut microbial metabolism, and the early-stage development and progression of chronic disease remain underexplored. Here, we examined the phenotypic impact of BMF variation in a cohort of generally-healthy, community dwelling adults with detailed clinical, lifestyle, and multi-omic data. We showed significant differences in microbially-derived blood plasma metabolites, gut bacterial genera, clinical chemistries, and lifestyle factors across BMF groups that have been linked to inflammation, cardiometabolic health, liver function, and CKD severity and progression. We found that the higher plasma levels of 3-indoxyl sulfate (3-IS), a microbially-derived metabolite associated with constipation, was in turn negatively associated with estimated glomerular filtration rate (eGFR), a measure of kidney function. Causal mediation analysis revealed that the effect of BMF on eGFR was significantly mediated by 3-IS. Finally, we identify self-reported diet, lifestyle, and psychological factors associated with BMF variation, which indicate several common-sense strategies for mitigating constipation and diarrhea. Overall, we suggest that aberrant BMF is an underappreciated risk factor in the development of chronic diseases, even in otherwise healthy populations.
Conflict of interest statement
DECLARATION OF INTERESTS L.H. is a former shareholder of Arivale. A.T.M. was a former employee of Arivale. Arivale is no longer a commercially operating company as of April 2019. The remaining authors report no competing interests.
Figures

Arivale participants were sampled for blood plasma and stool, in addition to filling out extensive diet, health, and lifestyle questionnaires. Clinical chemistries, untargeted metabolomics, and proteomics data were generated from blood plasma samples. Gut microbiome 16S rRNA amplicon sequencing data were generated from stool samples collected using at-home kits. BMF data were extracted from the questionnaire data as self-reported frequencies per week or day.

(A-F). POLR was used to regress BMF against the covariates (sex, age, BMI, eGFR, LDL, CRP, A1C, plus the first three principal components of genetic ancestry in the cohort, PC1, PC2, PC3). The result was that sex (p = 3.68E-16), BMI (p = 1.82E-3), age (p = 2.075E-3), and PCs1–3 (p < 0.00001, respectively) were significantly associated with BMF.

(A) Richness of amplicon sequence variants (ASVs) across BMF categories (ordinal BMF variable, Linear Regression, p = 2.85E-3). (B) Shannon Diversity across BMF categories (ordinal BMF variable, Linear Regression, p = 1.07E-3). (C) Pielou’s Evenness across BMF categories (ordinal BMF variable, Linear Regression, p = 8.5E-2).

46 significant taxa, in order of decreasing average relative abundance, with their z-scored, CLR-transformed abundances averaged within each BMF category plotted as a heatmap. Covariates included sex, age, BMI, eGFR, LDL, CRP, A1C, and PCs1–3. Asterisks denote the individual FDR-corrected significance threshold for the Wald Test p-value of the βBMF-coefficient for each BMF category, relative to the high-normal reference category. Rows without asterisks showed a significant overall model (FDR p-value <0.05), despite a lack of significance for the individual coefficients. (***): p < 0.0001, (**): 0.0001 < p < 0.01, (*): 0.01 < p < 0.05.

within each BMF category for all metabolites significantly associated with BMF. 11 significant blood plasma metabolites, with average z-scores within each BMF category plotted as a heatmap. Significant associations were identified using LIMMA, with FDR-corrected p-values of the ratio test between the main model and the null model. Here, the covariates included sex, age, BMI, eGFR, LDL, CRP, A1C, and PCs1–3. Asterisks denote metabolites with significant βBMF coefficient(s) in the linear regression model after FDR correction. (***): p < 0.0001, (**): 0.0001 < p < 0.01, (*): 0.01 < p < 0.05.

22 BMF-associated chemistries, identified using LIMMA models with FDR-corrected p-values of the ratio test between the main model and the null model, with average z-scores within each BMF category plotted as a heatmap. Here, the covariates included sex, age, BMI, eGFR, LDL, CRP, A1C, and PCs1–3. Asterisks denote FDR-corrected p-value thresholds for metabolites with significant βBMF coefficient(s) in the linear regression model. (***): p < 0.0001, (**): 0.0001 < p < 0.01, (*): 0.01 < p < 0.05.

Variables are colored by category: questions related to diet, exercise, and lifestyle (Diet/Lifestyle), and questions related to current digestive symptoms/function and health history (Health/Digestion). The BMF reference category was “high-normal” BMF (7–21 bowel movements per week). Each tick on the vertical axes represents a directional association in likelihood across the horizontal axis. The center line over the plots at x = 1.0 represents an equal likelihood of reporting an increase in number, intensity, frequency, or agreement (depending on the response variable) between the left side of the arrow on the vertical axis tick and the right side of the arrow on the vertical axis tick. A confidence interval that does not span the center line is significantly associated with the independent variable on the vertical axis tick. (*): FDR-corrected p-value < 0.05.

The average direct effect (ADE) of BMF on eGFR and the average causal mediated effect (ACME) of BMF on eGFR via 3-IS were found to be significant (N = 572; ADE −4.458, p = 0.012; ACME 1.343 p < 2E-16). The total effect and the proportion mediated terms did not pass our significance threshold of alpha=0.05.
Similar articles
-
Johnson-Martínez JP, Diener C, Levine AE, Wilmanski T, Suskind DL, Ralevski A, Hadlock J, Magis AT, Hood L, Rappaport N, Gibbons SM. Johnson-Martínez JP, et al. Cell Rep Med. 2024 Jul 16;5(7):101646. doi: 10.1016/j.xcrm.2024.101646. Cell Rep Med. 2024. PMID: 39019013 Free PMC article.
-
Bacterial metabolites and cardiovascular risk in children with chronic kidney disease.
Schlender J, Behrens F, McParland V, Müller D, Wilck N, Bartolomaeus H, Holle J. Schlender J, et al. Mol Cell Pediatr. 2021 Oct 22;8(1):17. doi: 10.1186/s40348-021-00126-8. Mol Cell Pediatr. 2021. PMID: 34677718 Free PMC article. Review.
-
Lai H, Li Y, He Y, Chen F, Mi B, Li J, Xie J, Ma G, Yang J, Xu K, Liao X, Yin Y, Liang J, Kong L, Wang X, Li Z, Shen Y, Dang S, Zhang L, Wu Q, Zeng L, Shi L, Zhang X, Tian T, Liu X. Lai H, et al. Gut Microbes. 2023 Jan-Dec;15(1):2197837. doi: 10.1080/19490976.2023.2197837. Gut Microbes. 2023. PMID: 37078654 Free PMC article. Clinical Trial.
-
The Impact of CKD on Uremic Toxins and Gut Microbiota.
Rysz J, Franczyk B, Ławiński J, Olszewski R, Ciałkowska-Rysz A, Gluba-Brzózka A. Rysz J, et al. Toxins (Basel). 2021 Mar 31;13(4):252. doi: 10.3390/toxins13040252. Toxins (Basel). 2021. PMID: 33807343 Free PMC article. Review.
-
Sun CY, Lin CJ, Pan HC, Lee CC, Lu SC, Hsieh YT, Huang SY, Huang HY. Sun CY, et al. Clin Nutr. 2019 Dec;38(6):2945-2948. doi: 10.1016/j.clnu.2018.11.029. Epub 2018 Dec 27. Clin Nutr. 2019. PMID: 30612852
References
Publication types
LinkOut - more resources
Full Text Sources
Research Materials
Miscellaneous