Frontiers in Endocrinology (Sep 2023)

A two-cohort study on the association between the gut microbiota and bone density, microarchitecture, and strength

  • Paul C. Okoro,
  • Eric S. Orwoll,
  • Curtis Huttenhower,
  • Curtis Huttenhower,
  • Curtis Huttenhower,
  • Curtis Huttenhower,
  • Xochitl Morgan,
  • Thomas M. Kuntz,
  • Lauren J. McIver,
  • Alyssa B. Dufour,
  • Alyssa B. Dufour,
  • Mary L. Bouxsein,
  • Mary L. Bouxsein,
  • Lisa Langsetmo,
  • Lisa Langsetmo,
  • Samaneh Farsijani,
  • Samaneh Farsijani,
  • Deborah M. Kado,
  • Deborah M. Kado,
  • Roberto Pacifici,
  • Shivani Sahni,
  • Shivani Sahni,
  • Douglas P. Kiel,
  • Douglas P. Kiel

DOI
https://doi.org/10.3389/fendo.2023.1237727
Journal volume & issue
Vol. 14

Abstract

Read online

The gut microbiome affects the inflammatory environment through effects on T-cells, which influence the production of immune mediators and inflammatory cytokines that stimulate osteoclastogenesis and bone loss in mice. However, there are few large human studies of the gut microbiome and skeletal health. We investigated the association between the human gut microbiome and high resolution peripheral quantitative computed tomography (HR-pQCT) scans of the radius and tibia in two large cohorts; Framingham Heart Study (FHS [n=1227, age range: 32 – 89]), and the Osteoporosis in Men Study (MrOS [n=836, age range: 78 – 98]). Stool samples from study participants underwent amplification and sequencing of the V4 hypervariable region of the 16S rRNA gene. The resulting 16S rRNA sequencing data were processed separately for each cohort, with the DADA2 pipeline incorporated in the16S bioBakery workflow. Resulting amplicon sequence variants were assigned taxonomies using the SILVA reference database. Controlling for multiple covariates, we tested for associations between microbial taxa abundances and HR-pQCT measures using general linear models as implemented in microbiome multivariable association with linear model (MaAslin2). Abundance of 37 microbial genera in FHS, and 4 genera in MrOS, were associated with various skeletal measures (false discovery rate [FDR] ≤ 0.1) including the association of DTU089 with bone measures, which was independently replicated in the two cohorts. A meta-analysis of the taxa-bone associations further revealed (FDR ≤ 0.25) that greater abundances of the genera; Akkermansia and DTU089, were associated with lower radius total vBMD, and tibia cortical vBMD respectively. Conversely, higher abundances of the genera; Lachnospiraceae NK4A136 group, and Faecalibacterium were associated with greater tibia cortical vBMD. We also investigated functional capabilities of microbial taxa by testing for associations between predicted (based on 16S rRNA amplicon sequence data) metabolic pathways abundance and bone phenotypes in each cohort. While there were no concordant functional associations observed in both cohorts, a meta-analysis revealed 8 pathways including the super-pathway of histidine, purine, and pyrimidine biosynthesis, associated with bone measures of the tibia cortical compartment. In conclusion, our findings suggest that there is a link between the gut microbiome and skeletal metabolism.

Keywords