
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39251661
71731
10.1038/s41598-024-71731-6
Article
Associations between gut microbiota and osteoporosis or osteopenia in a cohort of Chinese Han youth
Lai Junren 1
Gong Li 1
Liu Yan 1
Zhang Xuelian 12
Liu Wenqi 1
Han Meng 1
Zhou Duoqi duoqizhou@aqnu.edu.cn

13
Shi Shuiqin shishuiqin0511@163.com

13
1 https://ror.org/0127ytz78 grid.411412.3 0000 0001 0400 4349 Anhui Province Key Laboratory of the Biodiversity Study and Ecology Conservation in Southwest Anhui College of Life Sciences, Anqing Normal University, 246133 Anqing, Anhui People’s Republic of China
2 grid.9227.e 0000000119573309 State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, 14430 Urumqi, Xinjiang People’s Republic of China
3 Present Address: School of Life Sciences, 1318 North jixian Road, 246133 Anqing, Anhui People’s Republic of China
9 9 2024
9 9 2024
2024
14 209488 4 2024
30 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Osteoporosis (OP) is a common metabolic bone disease characterized by low bone mass and microstructural deterioration of bone. Changes in the composition and structure of gut microbiota (GM) are related to changes of bone mass and bone microstructure. However, the relationship between GM and bone mineral density (BMD) is complex, and data are especially scarce for Chinese Han youth. Therefore, 62 Chinese Han youth participants were recruited. Furthermore, according to the T–score evaluation criteria of the World Health Organization (WHO), we divided the BMD levels of participants into three groups: osteoporosis\BDL, osteopenia\BDM, normal bone density\BDH, and the associations between GM community and BMD groups were conducted. According to alpha and beta diversity analysis, significant differences were found in the microbial richness and composition between groups. The dominant phyla of GM in a cohort of Chinese Han youth were Bacteroidota (50.6%) and Firmicutes (41.6%). Anaerobic microorganisms, such as g_Faecalibacterium and g_Megamonas, account for the largest proportion in the gut, which were mainly Firmicutes phylum. The dominant genera and species in the three BMD groups were g_Prevotella, g_Bacteroides, g_Faecalibacterium, g_Megamonas, s_Prevotella copri, s_unclassified_g_Faecalibacterium, s_unclassified_g_Prevotella, s_unclassified_g_Bacteroides and s_Bacteroides plebeius. g_Faecalibacterium, g_Bacteroides and g_Ruminococcus differed between the BDH and BDL groups as well as between the BDH and BDM groups. LEfSe showed three genus communities and eight species communities were enriched in the three BMD groups, respectively. The associations between microbial relative abundance and T–score was not statistically significant by Spearman and regression analysis. In conclusion, the alpha diversity indexes in the BDH group were higher than in the BDL group, and several taxa were identified that may be the targets for diagnosis and therapy of OP.

Keywords

Associations
A cohort of Chinese Han youth
Gut microbiota
Osteoporosis
Osteopenia
Subject terms

Microbiology
Physiology
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The gut microbiota (GM) is considered to be a multifunctional “organ”, involved in immunity, metabolic regulation and other life activities. Studies have proved that there is an intimate relationship between the GM and human organs. The signal pathways between GM and different organs, can be divided into the gut–brain axis, gut–kidney axis, gut–liver axis, and the gut–bone axis1. Bone is a supporting organ composed of collagen and calcium phosphate, and it is also the main organ for storing minerals such as calcium and phosphorus. Homeostasis of the bone environment is maintained by osteoblasts and osteoclasts which are responsible for the dynamic process of bone remodeling2. When the balance is broken, it will result in low bone mass, destruction of bone tissue microstructure and eventually increased bone fragility, which is known as osteoporosis (OP) 3. The prevalence of OP in the Chinese population increased by 14.94% from 2003 to 2015, moreover, the current incidence of OP in the Chinese population aged over 50 is more than twice that in 20064. In view of the increasing prevalence of OP, prevention and control measures are even more important.

Osteocytes, the most abundant cells in bone tissue, are essential for the renewal and maintenance of the bone matrix. Osteocytes produce a variety of signals that induce osteoblast formation, such as prostaglandin E 2 (PGE 2), growth factors, and glycoproteins5. Recent advances have shown that, in addition to the respective signaling pathways and enhanced expression of genes in these osteocytes, further differentiation of osteoblasts and osteoclasts depends on metabolic reprogramming6. The GM and its metabolites, short–chain Fatty Acids (SCFAs), are intimately related to bone metabolism. A study has shown that SCFAs prevent bone loss in OP mouse model by metabolic reprogramming osteoclasts to inhibit osteoclast formation at an early stage7. SCFAs may also exert their effects by binding to free fatty acid receptor 2 (FFAR2) receptors expressed mainly by osteocytes8. These findings suggest that SCFAs have great potential to regulate bone health, and that GM metabolism produces large amounts of SCFAs. In a recent mouse model study, in vivo colonization of ovariectomized mice with Lactobacillus acidophilus enhanced trabecular and cortical bone microstructure, simultaneously increasing bone mineral density (BMD) and bone heterogeneity9. In a population study, the BMD and GM composition of 1776 people were analyzed using X–ray absorptiometry and 16S rRNA, and the results showed that BMD was related to the β diversity and species composition of the GM10. The diversity, composition and SCFAs of the GM were found to be closely related to the changes of bone mass and the development of bone metabolic diseases. Especially, the study of specific intestinal bacteria is very important for the prevention and treatment of OP and osteopenia.

Since the GM is affected by various factors such as diet, genetics, living environment and living habits, and some articles have pointed out that there are differences in the composition and structure of GM among different Chinese ethnic groups11. There were significant differences in the structure and composition of GM between the Plateau Tibetan and Plain Han populations with different BMD and T–score12. We only selected a cohort of Chinese Han youth of university students in Anhui Province, China, and which is located in the Yangtze River Delta region of China, one of the most Han populous regions. It is necessary to conduct research on population of different regions and nations, as this may supplement and validate previous research findings, especially for a nation with a large population base.

Materials and methods

Research participants

Non–sports majors of Han nationality were recruited as participants in Anqing Normal University, and the participants were all from Anhui Province, China. All participants were voluntary and gave informed consent to the content of the experiment. After excluding patients with diabetes/hypertension, osteoarthritis, fractures, gastroenteritis and peptic ulcer, 62 college students (28 men and 34 women) aged 19–25 years (21.60 ± 1.38 years) were finally recruited as the study participants. The basic information about the participants is shown (Table 1). This study was performed in line with the principles of the Declaration of Helsinki. Approval was granted by the Ethics Committee of Anqing Normal University (approval number: AQNU2018018H).Table 1 Participants basic information.

Gender	Number	Age (years)	Height (cm)	Weight (Kg)	
Male	28	20.00 ± 0.94	176.04 ± 5.68	70.30 ± 10.50	
Female	34	20.33 ± 1.16	160.59 ± 5.48	53.02 ± 7.06	

Standard for measurement and classification of BMD

The BMD at the lumbar spine (L1 to L4), total hip joint and femoral neck were measured by dual–energy X–ray absorptiometry (GE Lunar iDXA, GE Healthcare, Little Chalfont, UK) through a professional hospital. Before testing, participants were checked to ensure they had not had a barium meal, a radionuclide injection, or an injected or orally–administered contrast agent for CT or MRI within the previous 7 days. Using the World Health Organization (WHO) criteria, BMD is generally expressed as a T–score13. OP is diagnosed when the T–score of the lumbar spine, femoral neck and total hip joint is at least one T < −2.514. The specific criteria were that T < −2.5 was classified as osteoporosis\BDL (Hereinafter referred to as BDL), −2.5 < T < –1 was classified as osteopenia\BDM (BDM), and –1 < T < 1 was considered normal bone density\BDH (BDH) (Table 2).Table 2 Baseline characteristics of the participants.

BMD groups characteristics	BDH(n = 8)	BDM(n = 23)	BDL(n = 31)	P	
Malea	3(37.50)	13(56.52)	12(38.71)	0.130	
Femalea	5(62.50)	10(43.48)	19(61.29)	0.526	
Height cmb	167.07 ± 9.25	167.32 ± 9.57	167.22 ± 9.47	0.380	
Weight kgb	60.31 ± 12.37	60.72 ± 12.20	60.31 ± 12.00	0.532	
BMI/(kgm2)b	21.42 ± 2.64	21.50 ± 2.54	21.38 ± 2.47	0.350	
WHRb	0.827 ± 0.038	0.831 ± 0.039	0.830 ± 0.041	0.040*	
T-scoreb	 − 2.13 ± 0.98	 − 2.29 ± 0.88	 − 2.24 ± 0.93	0.758	
One-way ANOVA or Chi-square test were used to determine significance.

*P < 0.05; a: sex data are population (constituent ratio); b: height, weight, BMI, WHR and T-score values are mean ± SD.

Fecal collection and DNA extraction

We used disposable fecal collection boxes (GUHE Info Technology Co., Ltd, Hangzhou, China) to collect feces. The collection box consists of a sterile cotton swab and a fecal collection tube containing bacterial lysis preservation solution. Participants needed to use sterile cotton swabs to dip a small amount of feces into a fecal collection tube and kept in a portable Styrofoam box with ice packs (0–4 ℃). Fecal total DNA was extracted using the PowerMax (fecal/soil) DNA isolation kit (MoBio Laboratories, Carlsbad, CA). The quantity and quality of extracted DNA were measured using a NanoDrop ND–1000 photometer (Thermo Fisher Scientific, Waltham, MA, USA) and agarose gel electrophoresis respectively, to ensure that all samples met the sequencing requirements. Each faecal sample was frozen at −80 ℃ until sent to GUHE Info Technology Co., Ltd (Hangzhou, China) for sequencing analysis.

16S rRNA sequencing and quality control

In our study, the forward primer 515F (5′-GTGCCAGCMGCCGCGGTAA-3′) and the reverse primer 806R (5′-GGACTACHVGGGTWTCTAAT-3′) were used to amplify the V4 region of the bacterial 16S rRNA marker gene. PCR amplification was performed as follows: PCR reaction system (50 μL): Phusion High-Fidelity PCR Master Mix with HF Buffer 25 μL, 515 F and 806 R Primers (10 μmol/L) 3 μL each, DNA template 10 μL, DMSO 3 μL, DDH2O 6 μL. PCR reaction conditions: 98 ℃ 30 s, 98 ℃ 15 s, 58 ℃ 15 s, 72 ℃ 15 s, 25 cycles; 72 ℃ 1 min. The PCR amplification products were recovered and quantified, then purified using an Agencourt AMPure XP 60 mL Kit (Beckman Coulter, Brea, CA, USA) and quantified using the Qubit dsDNA HS Assay Kit (Vazyme, Nanjing, China) and finally two–terminal sequencing was performed using the Illlumina Hiseq 4000 Pair–End (Pair–End) 2 × 150 BP platform. According to Barcode sequence and primer sequence, the offline data were divided into the data of each sample. After truncating Barcode and primer sequences, the reads of each sample were spliced using VSEARCH v2.4.4 to obtain the Raw Tags. On this basis, the quality of the sequences was controlled and filtered, and the chimeric sequences were removed by UCHIME to obtain effective data, sequences with a similarity of 97% were clustered into operational taxonomic units (OTU’s) and the classification information for each sequence was annotated with the RDP classifier.

Bioinformatics analysis

Based on the clustering results, a Venn diagram was generated using the R package “Venn Diagram”. Moreover, QIIME and R Studio–based EasyAmplicon v1.14 software was used to compare the number of common and unique OTU’s among BDH, BDM and BDL groups. The R package “ggplot2” and its expansion package “ggtern” were used to generate histograms and ternary phase diagrams to analyze the composition of GM in BDH, BDM and BDL groups. The Kruskal–Wallis test was carried out with the function of “Kruskal. test” in R to analyze the taxa which differed significantly among BDH, BDM and BDL groups. The linear discriminant analysis (LDA) effect size (LEfSe) software v1.0 was used to do LEfSe analysis (LDA score threshold: 2) so as to find out the biomarkers. The co–occurrence network was analyzed by Cytoscape v3.9.1 software (https://cytoscape.org/), which reflects the interaction between genera and BMD groups. FastTree v2.1.3 (http://www.microbesonline.org/fast–tree/) was used to analyze the phylogeny and evolution of GM in a cohort of Chinese Han youth. The high–level phenotype prediction of GM in a cohort of Chinese Han youth was carried out by using BugBase (https://bugbase.cs.umn.edu/)15.

Statistical methods

The statistical tools R v3.5.3 and R Studio–based EasyAmplicon v1.14 software were used for statistical analysis and figure generation in this study. Data on height, weight, T–score, body mass index (BMI), and waist–to–hip ratio (WHR) are presented as the mean ± standard deviation (SD). The statistical differences were analyzed by the one–way ANOVA (conforming to a normal distribution), Chi–square test (categorical values) and Kruskal–Wallis (not conforming to a normal distribution). The Kruskal‒Wallis test was used for between-group comparisons of relative abundances at the phylum, genus and species levels and plotted by package “ggplot2” in R. R conducted alpha and beta diversity analysis with package “vegan”, which produced alpha diversity indexes (Chao1, ACE, Simpson and Shannon) and Beta diversity index (PCoA) at the OTU’s level. The significance of differences in the alpha diversity indexes was measured by the Kruskal–Wallis test between groups and was visualized as a box diagram using “phyloseq” package. Principal coordinate analysis (PCoA) was performed to evaluate microbial beta diversity based on Bray‒Curtis distance using R software (“phyloseq and vegan” package). Permutational multivariate analysis of variance (PERMANOVA) was performed to assess the between‒group differences in microbial communities using the “Adonis” function in the R package “vegan” with 9999 permutations. To explore the taxa differences in the three BMD groups, “Shapiro–Wilk. test” and “Levene’s. test” funtions in R were used for “normal distribution test” and “homogeneity test for variance”. According to the test results, the “Kruskal. test” function was used to compare the taxa differences in the three BMD groups, and the “pairwise–wilcox. test” function was used to calculate pairwise comparisons between groups. Post–hoc Benjamini–Hochberg Procedure correction was performed to determine pairwise differences. Spearman correlation analysis was applied to measure the correlations between microbiota abundance and T–score. Age, gender, height, weight, BMI and WHR as covariates, a regression model was constructed, and under the influence of the additive model (ADD, i.e. eliminating the influence of covariates). The results of regression were visually analyzed by GraphPad Prism v9.0.

Results

Basic information of participants in different BMD groups

There was a significant difference in WHR values. Not statistically significant differences were found in gender, height, weight, BMI, and T–score among the three groups (Table 2).

Comparisons of GM diversity between the three groups

The number of common and differential OTU’s among BDH, BDM and BDL groups was analyzed by sequence annotation for 62 samples. The total number of OTU’s obtained were 266 in the BDH group, 285 in the BDM group, 217 in the BDL group, with 197 in common (Fig. 1a).Fig. 1 Venn diagram and alpha\beta diversity analysis at OTU’s level. (a) Venn diagram. (b–e) Between-group comparisons of the alpha diversity of gut microbiota (ACE, Chao1, Shannon, Simpson). The results are presented in a box diagram, and P values were calculated by using the Kruskal–Wallis test. (f) Principal coordinate analysis (PCoA) of bacterial beta diversity based on the Bray‒Curtis distance. P values were calculated by Adonis. test with 9999 permutations.

Analysis of alpha diversity indicated that the community richness and diversity in the BDH group were significantly higher than that in the BDL group (Fig. 1b–e), while beta diversity, the PCoA plot did not show an obvious separation between groups. PERMANOVA (R2 = 0.02, P = 0.03) has showed significant differences in gut microbial composition between groups (Fig. 1f).

Gut microbial community and taxonomic composition

Taxonomic composition of the GM was clearly different in the three BMD groups. The dominant phyla in a cohort of Chinese Han youth were Bacteroidota (50.6%) and Firmicutes (41.6%). The relative abundances of phyla Bacteroidota, Firmicutes and Proteobacteria in the BDH group were 51.7, 37.9 and 9.2%, While in the BDM group they were 53.9, 37.9 and 6.5%. The relative abundance of the dominant phyla in the BDL group were Bacteroidota (46.1%), Firmicutes (49.2%) and Proteobacteria (4.0%) (Fig. 2a).Fig. 2 Taxonomic composition of the gut microbial community. (a) Phylum–level abundance plot. (b) Ternary phase diagram. Each dot represents an OTU’s, and the size, color, and position of the OTU’s represent its relative abundance, genus–level name, and grouping, respectively. (c) Genus–level abundance plot. (d) Species–level abundance plot.

Ternary phase diagram at the genus–level, g_Prevotella (31.2%), g_Bacteroides (16.9%), g_Faecalibacterium (13.2%) and g_Megamonas (4.6%) were the top 4 dominant genera in terms of average relative abundance for analysis (Fig. 2b). The relative abundance of g_Bacteroides (20.1%) in the BDH group was significantly higher than that in the BDM group (11.5%, P < 0.05) and higher than the BDL group (19.0%). The relative abundance of g_Prevotella (39.8%) in the BDM group was significantly higher than that in the BDH and BDL groups (28.1%, 25.6%, P < 0.01), and the relative abundance of g_Megamonas (7.8%) was significantly higher than that of the BDH and BDL groups (2.1%, 3.8%, P < 0.05). The relative abundance of g_Faecalibacterium (19.1%) in the BDL group was significantly higher than that of the BDH and BDM groups (11.8%, 8.6%, P < 0.01), and the relative abundance of g_Agathobacter (4.6%) was higher than that of the BDH and BDM groups (3.2%, 3.8%) (Fig. 2c).

Sample sequence annotation results at the species–level indicated that the dominant species were s_Prevotella copri (18.2%), s_unclassified_g_Faecalibacterium (12.7%), s_unclassified_g_Prevotella (9.2%), s_unclassified_g_Bacteroides (8.3%) and s_Bacteroides plebeius (7.2%). The relative abundance of s_unclassified_g_Bacteroides (11.1%) in the BDH group was significantly higher than that in the BDM group (4.1%, P < 0.01) and higher than the BDL group (10.0%). The relative abundance of s_unclassified_g_Prevotella (12.5%) in the BDM group was significantly higher than that in the BDH and BDL groups (8.2%, 6.9%, P < 0.01), while the relative abundance of s_Prevotella copri (22.3%) was also significantly higher than that of the BDH and BDL groups (16.8%, 15.5%, P < 0.01). The relative abundance of s_unclassified_g_Faecalibacterium (19.2%) in the BDL group was significantly higher than that of neither the BDH nor BDM groups (11.2%, 7.7%, P < 0.05), and the relative abundance of s_Bacteroides plebeius (7.6%) was higher than that of the BDH and BDM groups (7.2%, 6.9%) (Fig. 2d).

Identification of significant taxa differences and phenotype predication

The normal distribution and variance homogeneity test of genera and species in the three BMD groups were carried out, and the P values were all less than 0.05. Hence the data did not conform to the normal distribution and irregular variance that the Kruskal–Wallis test was used for the differences among the three groups. The pairwise-wilcox test was used to calculate pairwise comparisons between groups. The results showed that there were significant differences among the three groups in 9 genera and 1 specie: g_Faecalibacterium, g_Bacteroides, g_Prevotellaceae UCG–003, g_Lachnospiraceae UCG–010, g_Ruminococcaceae UCG–002, g_Ruminococcaceae NK4A214 group, g_Lachnospiraceae NK4A136 group, g_Lachnospira, s_Odoribacter splanchnicus and g_Ruminococcaceae UCG–005 (P < 0.05). In addition, 4 genera (g_Subdoligranulum, s_Ruminococcus bicirculans, g_Ruminococcus and g_Coprococcus) showed significant differences (P < 0.01) (Fig. 3a). The mean proportions of g_Faecalibacterium, g_Subdoligranulum and g_Ruminococcus were significantly higher in the BDL group than in the BDH group (P < 0.05), and the mean proportions of g_Faecalibacterium, g_Bacteroides, g_Ruminococcus and s_Ruminococcus bicirculans were significantly lower in the BDM group than in the BDH group (P < 0.05) (Fig. 3a).Fig. 3 Identification of significant taxa differences. (a) The Kruskal‒Wallis test was used for between-group comparisons of relative abundances at the genus-level and species-level. Only the taxa with a P value < 0.05 are shown. (b) LEfSe plot. The different colors represent different groups. The default LDA Score is 2, meaning that an LDA value greater than 2 indicates a statistical difference between groups. Note: * P < 0.05, ** P < 0.01.

According to the LEfSe analysis, g_Bacteroides, s_unclassified_g_Bacteroides, s_gut_metagenome_g__Lachnospira and s_Lachnospiraceae bacterium TF01–11 were enriched in the BDH group. In the BDM group, g_Kocuria, s_Desulfovibrio and s_unclassified_g_Kocuria were enriched. s_unclassified_g_Butyricimonas, s_unclassified_bacterium_g__norank_f__Saccharimonadaceae, s_unclassified_organism_g__UCG–002 and s_unclassified_g_Ruminococcus were enriched in the BDL group (Fig. 3b).

We used BugBase phenotype prediction, which showed that g_Faecalibacterium accounted for the highest proportion of the anaerobic phenotype, and which in the BDL group was higher than in the BDH and BDM groups (Fig. 4a). The relative abundance of g_Megamonas and g_Prevotella of the anaerobic phenotype were higher in the BDM group than in the BDH and BDL groups (Fig. 4a). The prevailing gut microbial communities within the facultative anaerobic phenotype comprised g_Haemophilus, g_unclassified_f__Enterobacteriaceae, g_Erysipelotrichaceae UCG–003, g_Lactobacillus and g_Escherichia Shigella (Fig. 4b, Table S1). We ascertained that g_Sutterella, the predominant genus existing in the aerobic phenotype, was enriched in the BDH group. Meanwhile, in the BDM group of the aerobic phenotype, the relative abundance of g_Rothia was greater than that in the BDH and BDL groups (Fig. 4c).Fig. 4 Phenotype prediction and co–occurrence network analysis of bacterial taxa at genus-level. (a) Anaerobic. (b) Facultative Anaerobic. (c) Aerobic. (d) Co–occurrence network. The co–occurrence network explores the interaction between genera and BMD groups.

Co–occurrence network analysis explores the correlation between genera and BMD groups, reflecting the range and strength of correlation between genera and BMD groups. The connecting lines between genera and BMD groups are considered as degrees, The number of degrees represents the influence of a genus in one or more BMD groups. The weighted degrees represent the strength of correlation the genera are with BMD groups. The higher the weighted degrees, the more strength of correlation the genera are with BMD groups, and the deeper the impact. The degrees of g_Megamonas, g_Faecalibacterium, g_Roseburia, g_Lachnospira, g_Prevotella, g_Bacteroides, g_Agathobacter, g_Phascolarctobacterium, g_Escherichia Shigella, and g_Dialister were three, which all had an impact in the three BMD groups. The weighted degrees of g_Prevotella, g_Bacteroides, g_Faecalibacterium, g_Megamonas, and g_Agathobacter were 45,800.62, 28,527.11, 20,643.70, 7525.19, 4998.01 respectively, with genera of Bacteroidota and Firmicutes phyla being highly intimately related to the three BMD groups. g_Alloprevotella, g_unclassified_f__Enterobacteriaceae, g_Ruminococcus, g_Sutterella and g_Haemophilus with one degree for the BDH group, while g_Eubacterium eligens group and g_Subdoligranulum with one degree for the BDL group, and which mean those genera only have an impact in specific BMD group. The BDH and BDL groups shared g_unclassified_f__Lachnospiraceae, while the BDM and BDL groups shared g_Lachnoclostridium (Fig. 4d, Table S2).

Associations of GM relative abundance and T–score at phylum–level and genus–level

No associations were found between gut microbiota abundance and T-score at phylum and genus levels by spearman correlation analysis (Table S3). After adjusting for age, BMI and gender in regression analysis, also no associations between gut microbiota abundance and T-score were found (Fig. 5, Table S3).Fig. 5 Regression analysis. (a): Phylum-level; (b): genus-level. Regression analysis reflects the associations between GM relative abundance and T-score.

Analysis of the evolutionary relationship

At the genus–level, Fast–Tree software was used to analyze the evolutionary relationship of GM in the cohort of Chinese Han youth. The results showed that, with the exception of g_Prevotella and g_Bacteroides, the GM of the cohort of Chinese Han youth belonged to the same branch. The dominant genus was g_Faecalibacterium, which had the highest homology with the other genus, followed by g_Escherichia Shigella (Fig. 6).Fig. 6 Phylogenetic tree of bacterial taxa. Each branch of the phylogenetic tree represents a genus, and the branches are colored according to the high taxonomic level of the genera. The branch length is the phylogenetic distance between the two genera, that is, the degree of genera difference. The bar chart shows the proportion of reads of genera in different BMD groups.

Discussion

The GM is closely related to bone health, and the diversity and richness of the GM are one of the important factors affecting bone health. Studies have shown that poor lifestyle, such as drinking alcohol, insufficient sleep, smoking, and irregular diet, are the critical factors leading to OP, and the proportion of contemporary youth with unhealthy lifestyles is increasing16,17. Therefore, our study analyzed the diversity, composition and structure, phenotypes and evolutionary relationship of GM community among different BMD levels in a cohort of Chinese Han youth, explored their differences and found some interesting phenomena.

Our findings performed that there were significant differences in the alpha diversity of GM between groups, which aligns with certain earlier studies18,19. PERMANOVA of beta diversity revealed that the microbial composition between groups was significantly different, and which is consistent previous study20. But there was also a similar study conducted contradict results21. The contradictory outcomes of alpha and beta diversity between different studies might be attributed to the vast heterogeneity of cohorts, also including age, sex, dietary habits, etc. On the other hand, this is an advantage of selecting a cohort of Chinese Han youth in Yangtze River Delta region as the participants, where they may have similar dietary habits, and to exclude the effects of age, dietary habits, locations.

We further analyzed the taxonomic differences. Compared with the BDH group, the relative abundance of Bacteroidota and Proteobacteria in the BDL group was decreased, while the relative abundance of Firmicutes was increased. A study of Japanese women showed that the incidence of fractures in the lower Bacteroides group was 5.6 times higher than that in the normal group22. Another study examined stool samples from 102 elderly Chinese patients with OP using 16S rRNA sequencing and showed that the abundance of Firmicutes was negatively correlated with BMD and T–score23. These findings suggested that changes in the proportion of gut dominant phyla were associated with BMD.

The g_Prevotella is one of the most abundant taxa in the oral and GM, and they are key participants in maintaining human metabolic homeostasis24. In the present study, we also annotated g_Prevotella as the dominant genus (with an average relative abundance of 31.2%). After groups annotating, it was found that the relative abundance of g_Prevotella in the BDL group (25.6%) was lower than that in the BDH group (28.1%), while the relative abundance of g_Prevotella in the BDM group (39.8%) was significantly higher than that in the BDH group (28.1%). Studies have shown that the dysbiosis of GM in OP mice induced by ovariectomized (OVX) was regulated by Prevotellaceae, which ameliorated bone loss in OP mice and has potential for preventing and treating OP25,26. Another OVX mice experiment showed that the abundance of g_Prevotella increased rapidly in the early and middle stages of OP and remained high for a certain period of time, but gradually decreased in the later stages of OP27, which is similar to our result. g_Prevotella may be a potential diagnostic biomarker.

Another dominant genus in our study is g_Faecalibacterium (with an average relative abundance of 13.2%), which is also the dominant genus in the BDL group. Literature shows that compared to the normal control group, Faecalibacterium prausnitzii is significantly enriched in patients with primary osteoporosis (PO) (0.331 vs. 0.132)28. We also found that the relative abundance of g_Megamonas in the BDM group (7.8%) was significantly higher than that in the BDH and BDL groups (2.1%, 3.8%). The relative abundance of Megamonas is positively correlated with femoral BMD29. The relative abundance of Megamonas in the gut of postmenopausal osteoporosis (PMOP) women is higher than that of the normal control group30. The reason for the difference may be a lack of estrogen.

We found that s_unclassified__g_Ruminococcus and s_unclassified__g_Butyrisimonas were enriched in the BDL group, but s_gut_metagenome__g_Lachnospira was enriched in the BDH group. Genus Butyricimonas, Lachnospira and Ruminococcus are the main members of the Clostridiales. This paradoxical result has also been reflected in previous studies. A study of the GM in an elderly Chinese population reported that g_Clostridium plays an important role in maintaining bone health, while another assessing study of analysis estimated heel bone mineral density (eBMD) genome–wide association research (GWAS) in 426,824 participants found a negative correlation between g_Clostridium abundance and BMD19,31. g_Lachnospira is closely related to heel BMD32. Probiotic metabolite Kefir fermentation peptide can regulate GM in OVX mice, leading to an increase in g_Ruminococcus33. SCFAs produced by g_Butyrisimonas participate in bone homeostasis regulation34. The results indicate that different genus of Clostridiales has different effects on regulating BMD, and their impact on BMD may be constrained by various factors such as region, population, and age.

In our study, the relative abundance of g_Kocuria in the BDM group was significantly higher than in the BDH and BDL groups, and g_Subdoligranulum was related to the BDL group. g_Kocuria infection is associated with low immunity, which has been shown to induce diseases such as peritonitis, myelodysplastic syndromes (MDS) and patients with short bowel syndrome are at increased risk of opportunistic infections such as with g_Kocuria, which may increase the likelihood of developing foot osteomyelitis35. g_Kocuria infection is most likely related to the expression of IL–6, IL–8, IFN–γ36. IFN–γ promotes osteoclast maturation in the late stage of osteoclast formation37. g_Subdoligranulum was found to stimulate TH17–cells proliferation, serum RA–relevant IgG autoantibodies, and joint swelling, while mice colonized with g_Subdoligranulum isolates developed OA, further analysis found that g_Subdoligranulum elicited CD4 T–cells response38. CD4 T–cells enter the thymus gland and differentiate to eventually form CD4 + T–cells, CD4 + T–cells can produce proinflammatory cytokines involved in inflammatory bone loss, increase osteoblast apoptosis, and promote osteoclastogenesis39. We can see that there is an intimate link among g_Kocuria and g_Subdoligranulum, the immune system, and BMD, and that g_Kocuria and g_Subdoligranulum may disrupt bone homeostasis by affecting the immune system. Interestingly, we annotated the s_Eubacterium eligens group (g_Eubacterium) was associated with the BDL group. g_Eubacterium is one of the core taxa of the human intestinal community, and its fermented glucose products typically contain large amounts of butyrate, acetate, or propionate40. Butyrate, acetate and propionate are the most common SCFAs, and SCFAs stimulate bone formation by increasing the production of bone Sialoprotein (BSP) and osteopontin (OPN) in osteoblasts, while also inhibiting osteoclast differentiation by inhibiting protein deacetylase and activating free fatty acid receptors (FFARs)41,42. Therefore, the appearance of g_Eubacterium in the gut tract of the BDL individuals is a very positive signal, indicating that g_Eubacterium has potential application value in protecting bone mass and regulating bone metabolism. High–fiber diets and Mediterranean diets have been shown to increase the relative abundance of g_Eubacterium43,44. This may provide a new idea for preventing and treating OP.

In the present study, Anaerobic bacteria were the dominant in the GM community, and a higher relative abundance of facultative anaerobic and aerobic bacteria in the BDH group compared to the BDM and BDL groups. The relative abundance of anaerobic bacteria was the highest between osteoarthritis (OA) and the control group, and the relative abundance of facultative anaerobic bacteria in OA was lower than that in the control group, but the relative abundance of aerobic bacteria in OA was higher than that in the control group45. Exercise, diet, and medication may all affect the gut environment by altering intestinal permeability, stimulating the production of immune factors46. A GWAS study investigated the composition of the GM in 181 elderly adults and found that some of the bacteria of Proteobacteria, such as g_Sutterella, f_Enterobacteriaceae and g_Citrobacter were highly enriched in the elder OP group21. Other studies have found that g_Sutterella secrete IgA protease, which can break down IgA in the intestinal mucosa, destroy the intestinal antibacterial immune system, and may cause intestinal inflammation. Intestinal inflammation increases intestinal permeability, decreases intestinal absorption of SCFAs and Ca2+, and increases the risk of osteoporosis47,48. In our study, we found that g_Sutterella was related to the BDH group, and the relative abundance of g_Sutterella was the highest in the aerobic phenotype. Some nutritional and psychiatric studies have pointed out that increased intake of dietary fatty acids leads to an increase in the abundance of g_Sutterella. g_Sutterella is a major component of the GM in more than half of children with autism, and altered intestinal permeability may be a cause of autism in children49,50. Therefore, we hypothesize that the high abundance of g_Sutterella was influenced by dietary and emotional changes, but this hypothesis needs further investigation. In addition, it is necessary to pay intimate attention to the psychological and mental state of the participants in the BDH group to avoid possible autistic tendencies.

In conclusion, our study demonstrates that the alpha diversity of GM in the BDL group are low. g_Prevotella and g_Sutterella play an important role in the three BMD groups, and these two taxa may be good candidates for monitoring and controlling OP.

However, our study has several limitations. First, the relatively small sample size might also have reduced the statistical power to detect subtle differences among the taxa, this may be one of the reasons why we did not found associations between GM and T-score, and differences in identification of dominant taxa at different perspectives. But the sample size of our study approached that of most similar studies20,51. Second, high–level BMI or WHR may lead to differences in alpha and beta diversity52,53, and some individuals of this study had high–level of BMI and there was a significant difference in WHR between groups. After regression analysis, not significant associations were found between GM abundance and T–score, and caution should be taken in interpreting these results. In addition, the study was not in depth enough, especially since it relied solely on 16S rRNA sequencing of stool samples to study the GM, a causal link between changes in GM and BMD cannot be determined. Therefore, it is necessary to consider the confounding factors of individuals in larger cohorts in future studies in order to ensure that the study is more rigorous.

Conclusion

The results of this study showed that the composition and structure of GM were changed. The richness and diversity of GM in the BDL were lower than in the BDH and BDM groups. Genus Prevotella, Eubacterium and Sutterella are potential candidate microbiota with diagnosis and therapy of OP.

Supplementary Information

Supplementary Tables.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-71731-6.

Acknowledgements

We appreciate all the volunteers who participated in this study.

Author contributions

Contributors: Author Duoqi Zhou designed the study. He is guarantor. Author Junren Lai prepared the first draft of the paper. Authors Yan Liu, Xuelian Zhang and Wenqi Liu contributed to the collected data. Author Shuiqin Shi and Li Gong were responsible for statistical analysis of the data. All authors revised the paper critically for intellectual content and approved the final version. All authors agree to be accountable for the work and to ensure that any questions relating to the accuracy and integrity of the paper are investigated and properly resolved.

Data availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Butel MJ Probiotics, gut microbiota and health Medecine et maladies infectieuses 2014 44 1 1 8 10.1016/j.medmal.2013.10.002 24290962
Butel, M. J. Probiotics, gut microbiota and health. Medecine et maladies infectieuses 44(1), 1–8 (2014).24290962 10.1016/j.medmal.2013.10.002
2. Yan Q Cai L Guo W New advances in improving bone health based on specific gut microbiota Front. Cell. Infect. Microbiol. 2022 12 821429 10.3389/fcimb.2022.821429 35860378
Yan, Q., Cai, L. & Guo, W. New advances in improving bone health based on specific gut microbiota. Front. Cell. Infect. Microbiol. 12, 821429 (2022).35860378 10.3389/fcimb.2022.821429
3. Lane JM Russell L Khan SN Osteoporosis Clin. Orthop. Relat. Res. 2000 372 139 150 10.1097/00003086-200003000-00016
Lane, J. M., Russell, L. & Khan, S. N. Osteoporosis. Clin. Orthop. Relat. Res. 372, 139–150 (2000).10.1097/00003086-200003000-00016
4. Chen P Li Z Hu Y Prevalence of osteoporosis in China: A meta-analysis and systematic review BMC Public Health 2016 16 1 1039 10.1186/s12889-016-3712-7 27716144
Chen, P., Li, Z. & Hu, Y. Prevalence of osteoporosis in China: A meta-analysis and systematic review. BMC Public Health 16(1), 1039 (2016).27716144 10.1186/s12889-016-3712-7
5. Compton JT Lee FY A review of osteocyte function and the emerging importance of sclerostin J. Bone Joint Surg. Am. 2014 96 19 1659 1668 10.2106/JBJS.M.01096 25274791
Compton, J. T. & Lee, F. Y. A review of osteocyte function and the emerging importance of sclerostin. J. Bone Joint Surg. Am. 96(19), 1659–1668 (2014).25274791 10.2106/JBJS.M.01096
6. Shi K Li F Tao Y Direct reprogramming in bone and joint degenerative diseases: Applications, obstacles and directions Current Stem Cell Res. Therapy 2023 18 6 766 778 10.2174/1574888X17666220810142943
Shi, K. et al. Direct reprogramming in bone and joint degenerative diseases: Applications, obstacles and directions. Current Stem Cell Res. Therapy 18(6), 766–778 (2023).10.2174/1574888X17666220810142943
7. Lucas S Omata Y Hofmann J Short-chain fatty acids regulate systemic bone mass and protect from pathological bone loss Nat Commun. 2018 9 1 55 10.1038/s41467-017-02490-4 29302038
Lucas, S. et al. Short-chain fatty acids regulate systemic bone mass and protect from pathological bone loss. Nat Commun. 9(1), 55 (2018).29302038 10.1038/s41467-017-02490-4
8. Montalvany-Antonucci CC Duffles LF de Arruda JAA Short-chain fatty acids and FFAR2 as suppressors of bone resorption Bone 2019 125 112 121 10.1016/j.bone.2019.05.016 31100533
Montalvany-Antonucci, C. C. et al. Short-chain fatty acids and FFAR2 as suppressors of bone resorption. Bone 125, 112–121 (2019).31100533 10.1016/j.bone.2019.05.016
9. Dar HY Shukla P Mishra PK Lactobacillus acidophilus inhibits bone loss and increases bone heterogeneity in osteoporotic mice via modulating Treg-Th17 cell balance Bone Rep. 2018 8 46 56 10.1016/j.bonr.2018.02.001 29955622
Dar, H. Y. et al. Lactobacillus acidophilus inhibits bone loss and increases bone heterogeneity in osteoporotic mice via modulating Treg-Th17 cell balance. Bone Rep. 8, 46–56 (2018).29955622 10.1016/j.bonr.2018.02.001
10. Ling CW Miao Z Xiao ML The association of gut microbiota with osteoporosis is mediated by amino acid metabolism: Multiomics in a large cohort J. Clin. Endocrinol. Metab 2021 106 10 e3852 e3864 10.1210/clinem/dgab492 34214160
Ling, C. W. et al. The association of gut microbiota with osteoporosis is mediated by amino acid metabolism: Multiomics in a large cohort. J. Clin. Endocrinol. Metab 106(10), e3852–e3864 (2021).34214160 10.1210/clinem/dgab492
11. Cheng M Ge X Zhong C Micro-coevolution of host genetics with gut microbiome in three Chinese ethnic groups J. Genetics Genom. 2021 48 11 972 983 10.1016/j.jgg.2021.09.002
Cheng, M. et al. Micro-coevolution of host genetics with gut microbiome in three Chinese ethnic groups. J. Genetics Genom. 48(11), 972–983 (2021).10.1016/j.jgg.2021.09.002
12. Fan, W. Q. Study on the differences in gut microbiota and short chain fatty acids in feces between Tibetan and Han populations with different bone densities in Hongyuan and Chengdu. 26(2), 26-30. 10.27342/d.cnki.gscdu.2021.006747 (Sichuan University, 2021).
13. Kanis JA Assessment of fracture risk and its application to screening for postmenopausal osteoporosis: synopsis of a WHO report WHO Study Group. Osteoporos. Int. 1994 4 6 368 381 10.1007/BF01622200 7696835
Kanis, J. A. Assessment of fracture risk and its application to screening for postmenopausal osteoporosis: synopsis of a WHO report WHO Study Group.. Osteoporos. Int. 4(6), 368–381 (1994).7696835 10.1007/BF01622200
14. Nayak S Edwards DL Saleh AA Systematic review and meta-analysis of the performance of clinical risk assessment instruments for screening for osteoporosis or low bone density Osteoporos. Int. 2015 26 5 1543 1554 10.1007/s00198-015-3025-1 25644147
Nayak, S. et al. Systematic review and meta-analysis of the performance of clinical risk assessment instruments for screening for osteoporosis or low bone density. Osteoporos. Int. 26(5), 1543–1554 (2015).25644147 10.1007/s00198-015-3025-1
15. Duan G Li L Deciphering the mechanism of jujube vinegar on hyperlipoidemia through gut microbiome based on 16S rRNA, BugBase analysis, and the stamp analysis of KEEG Front. Nutr. 2023 10 1160069 10.3389/fnut.2023.1160069 37275638
Duan, G. & Li, L. Deciphering the mechanism of jujube vinegar on hyperlipoidemia through gut microbiome based on 16S rRNA, BugBase analysis, and the stamp analysis of KEEG. Front. Nutr. 10, 1160069 (2023).37275638 10.3389/fnut.2023.1160069
16. Ambrosio MR Aliberti L Gagliardi I Bone health in adolescence Minerva Obstet. Gynecol. 2021 73 6 662 677 10.23736/S2724-606X.20.04713-9 34905874
Ambrosio, M. R. et al. Bone health in adolescence. Minerva Obstet. Gynecol. 73(6), 662–677 (2021).34905874 10.23736/S2724-606X.20.04713-9
17. Pitukcheewanont P Austin J Chen P Bone health in children and adolescents: risk factors for low bone density Pediatr. Endocrinol. Rev. PER 2013 10 3 318 335 23724439
Pitukcheewanont, P. et al. Bone health in children and adolescents: risk factors for low bone density. Pediatr. Endocrinol. Rev. PER 10(3), 318–335 (2013).23724439
18. Wang J Wang Y Gao W Diversity analysis of gut microbiota in osteoporosis and osteopenia patients PeerJ. 2017 5 e3450 10.7717/peerj.3450 28630804
Wang, J. et al. Diversity analysis of gut microbiota in osteoporosis and osteopenia patients. PeerJ. 5, e3450 (2017).28630804 10.7717/peerj.3450
19. Wang Y Gao X Lv J Gut microbiome signature are correlated with bone mineral density alterations in the Chinese elders Front. Cell. Infect. Microbiol. 2022 12 827575 10.3389/fcimb.2022.827575 35433497
Wang, Y. et al. Gut microbiome signature are correlated with bone mineral density alterations in the Chinese elders. Front. Cell. Infect. Microbiol. 12, 827575 (2022).35433497 10.3389/fcimb.2022.827575
20. Rettedal EA Ilesanmi-Oyelere BL Roy NC The Gut microbiome is altered in postmenopausal women with osteoporosis and osteopenia JBMR plus 2021 5 3 e10452 10.1002/jbm4.10452 33778322
Rettedal, E. A. et al. The Gut microbiome is altered in postmenopausal women with osteoporosis and osteopenia. JBMR plus 5(3), e10452 (2021).33778322 10.1002/jbm4.10452
21. Das M Cronin O Keohane DM Gut microbiota alterations associated with reduced bone mineral density in older adults Rheumatol. (Oxford, England) 2019 58 12 2295 2304 10.1093/rheumatology/kez302
Das, M. et al. Gut microbiota alterations associated with reduced bone mineral density in older adults. Rheumatol. (Oxford, England) 58(12), 2295–2304 (2019).10.1093/rheumatology/kez302
22. Ozaki D Kubota R Maeno T Association between gut microbiota, bone metabolism, and fracture risk in postmenopausal Japanese women Osteoporos. Int. 2021 32 1 145 156 10.1007/s00198-020-05728-y 33241467
Ozaki, D. et al. Association between gut microbiota, bone metabolism, and fracture risk in postmenopausal Japanese women. Osteoporos. Int. 32(1), 145–156 (2021).33241467 10.1007/s00198-020-05728-y
23. Li C Huang Q Yang R Gut microbiota composition and bone mineral loss-epidemiologic evidence from individuals in Wuhan, China Osteoporos. Int. 2019 30 5 1003 1013 10.1007/s00198-019-04855-5 30666372
Li, C. et al. Gut microbiota composition and bone mineral loss-epidemiologic evidence from individuals in Wuhan, China. Osteoporos. Int. 30(5), 1003–1013 (2019).30666372 10.1007/s00198-019-04855-5
24. Balakrishnan B Luckey D Bodhke R Prevotella histicola protects from arthritis by expansion of allobaculum and augmenting butyrate production in humanized mice Front. Immunol. 2021 12 609644 10.3389/fimmu.2021.609644 34017324
Balakrishnan, B. et al. Prevotella histicola protects from arthritis by expansion of allobaculum and augmenting butyrate production in humanized mice. Front. Immunol. 12, 609644 (2021).34017324 10.3389/fimmu.2021.609644
25. Wang Z Chen K Wu C An emerging role of Prevotella histicola on estrogen deficiency-induced bone loss through the gut microbiota-bone axis in postmenopausal women and in ovariectomized mice Am. J. Clin. Nutr. 2021 114 4 1304 1313 10.1093/ajcn/nqab194 34113963
Wang, Z. et al. An emerging role of Prevotella histicola on estrogen deficiency-induced bone loss through the gut microbiota-bone axis in postmenopausal women and in ovariectomized mice. Am. J. Clin. Nutr. 114(4), 1304–1313 (2021).34113963 10.1093/ajcn/nqab194
26. Zhang YW Cao MM Li YJ Fecal microbiota transplantation ameliorates bone loss in mice with ovariectomy-induced osteoporosis via modulating gut microbiota and metabolic function J. Orthop. Transl. 2022 37 46 60
Zhang, Y. W. et al. Fecal microbiota transplantation ameliorates bone loss in mice with ovariectomy-induced osteoporosis via modulating gut microbiota and metabolic function. J. Orthop. Transl. 37, 46–60 (2022).
27. Ma S Qin J Hao Y Structural and functional changes of gut microbiota in ovariectomized rats and their correlations with altered bone mass Aging 2020 12 11 10736 10753 10.18632/aging.103290 32484785
Ma, S. et al. Structural and functional changes of gut microbiota in ovariectomized rats and their correlations with altered bone mass. Aging 12(11), 10736–10753 (2020).32484785 10.18632/aging.103290
28. Xu Z Xie Z Sun J Gut microbiome reveals specific dysbiosis in primary osteoporosis Front. Cell. Infect. microbiol. 2020 10 160 10.3389/fcimb.2020.00160 32373553
Xu, Z. et al. Gut microbiome reveals specific dysbiosis in primary osteoporosis. Front. Cell. Infect. microbiol. 10, 160 (2020).32373553 10.3389/fcimb.2020.00160
29. Li J Liu C Xu Y Gut microbiota alterations in adolescent idiopathic scoliosis are associated with aberrant bone homeostasis Orthop. Surg. 2024 10.1111/os.14019 39238187
Li, J. et al. Gut microbiota alterations in adolescent idiopathic scoliosis are associated with aberrant bone homeostasis. Orthop. Surg.10.1111/os.14019 (2024).39238187 10.1111/os.14019
30. Yang X Chang T Yuan Q Changes in the composition of gut and vaginal microbiota in patients with postmenopausal osteoporosis Front. Immunol. 2022 13 930244 10.3389/fimmu.2022.930244 36032115
Yang, X. et al. Changes in the composition of gut and vaginal microbiota in patients with postmenopausal osteoporosis. Front. Immunol. 13, 930244 (2022).36032115 10.3389/fimmu.2022.930244
31. Ni JJ Yang XL Zhang H Assessing causal relationship from gut microbiota to heel bone mineral density Bone 2021 143 115652 10.1016/j.bone.2020.115652 32971307
Ni, J. J. et al. Assessing causal relationship from gut microbiota to heel bone mineral density. Bone 143, 115652 (2021).32971307 10.1016/j.bone.2020.115652
32. Wei J Karsenty G An overview of the metabolic functions of osteocalcin Curr. Osteoporos. Rep. 2015 13 3 180 185 10.1007/s11914-015-0267-y 25809656
Wei, J. & Karsenty, G. An overview of the metabolic functions of osteocalcin. Curr. Osteoporos. Rep. 13(3), 180–185 (2015).25809656 10.1007/s11914-015-0267-y
33. Sun P Zhang Y Wei Z Effect of Qing’e decoction on Leptin/Leptin receptor and bone metabolism in naturally aging rats Evid.Based Complementary Altern. Med. eCAM 2020 2020 2532081 10.1155/2020/2532081
Sun, P. et al. Effect of Qing’e decoction on Leptin/Leptin receptor and bone metabolism in naturally aging rats. Evid.Based Complementary Altern. Med. eCAM 2020, 2532081 (2020).10.1155/2020/2532081
34. Rosenblatt A Faintuch J Cecconello I Androgen and estrogen shifts in men before and after bariatric surgery and links to vitamins and trace elements Int. J. Vitam. Nutr. Res. 2016 86 3–4 198 241 28436761
Rosenblatt, A., Faintuch, J. & Cecconello, I. Androgen and estrogen shifts in men before and after bariatric surgery and links to vitamins and trace elements. Int. J. Vitam. Nutr. Res. 86(3–4), 198–241 (2016).28436761
35. Kandi V Palange P Vaish R Emerging bacterial infection: Identification and clinical significance of kocuria species Cureus 2016 8 8 e731 27630804
Kandi, V. et al. Emerging bacterial infection: Identification and clinical significance of kocuria species. Cureus 8(8), e731 (2016).27630804
36. Xu J Chen X Yu S Effects of early intervention with sodium butyrate on gut microbiota and the expression of inflammatory cytokines in neonatal piglets PloS ONE 2016 11 9 e0162461 10.1371/journal.pone.0162461 27611998
Xu, J. et al. Effects of early intervention with sodium butyrate on gut microbiota and the expression of inflammatory cytokines in neonatal piglets. PloS ONE 11(9), e0162461 (2016).27611998 10.1371/journal.pone.0162461
37. Zhang W Dang K Huai Y Osteoimmunology: The regulatory roles of T lymphocytes in osteoporosis Front. Endocrinol. 2020 11 465 10.3389/fendo.2020.00465
Zhang, W. et al. Osteoimmunology: The regulatory roles of T lymphocytes in osteoporosis. Front. Endocrinol. 11, 465 (2020).10.3389/fendo.2020.00465
38. Chriswell ME Lefferts AR Clay MR Clonal IgA and IgG autoantibodies from individuals at risk for rheumatoid arthritis identify an arthritogenic strain of Subdoligranulum Sci. Transl. Med. 2022 14 668 eabn5166 10.1126/scitranslmed.abn5166 36288282
Chriswell, M. E. et al. Clonal IgA and IgG autoantibodies from individuals at risk for rheumatoid arthritis identify an arthritogenic strain of Subdoligranulum. Sci. Transl. Med. 14(668), eabn5166 (2022).36288282 10.1126/scitranslmed.abn5166
39. Zhang W Gao R Rong X Immunoporosis: Role of immune system in the pathophysiology of different types of osteoporosis Front. Endocrinol. 2022 13 965258 10.3389/fendo.2022.965258
Zhang, W. et al. Immunoporosis: Role of immune system in the pathophysiology of different types of osteoporosis. Front. Endocrinol. 13, 965258 (2022).10.3389/fendo.2022.965258
40. Mukherjee A Lordan C Ross RP Gut microbes from the phylogenetically diverse genus Eubacterium and their various contributions to gut health Gut Microbes 2020 12 1 1802866 10.1080/19490976.2020.1802866 32835590
Mukherjee, A. et al. Gut microbes from the phylogenetically diverse genus Eubacterium and their various contributions to gut health. Gut Microbes 12(1), 1802866 (2020).32835590 10.1080/19490976.2020.1802866
41. Katono T Kawato T Tanabe N Sodium butyrate stimulates mineralized nodule formation and osteoprotegerin expression by human osteoblasts Arch. oral biol. 2008 53 10 903 909 10.1016/j.archoralbio.2008.02.016 18406397
Katono, T. et al. Sodium butyrate stimulates mineralized nodule formation and osteoprotegerin expression by human osteoblasts. Arch. oral biol. 53(10), 903–909 (2008).18406397 10.1016/j.archoralbio.2008.02.016
42. Rahman MM Kukita A Kukita T Two histone deacetylase inhibitors, trichostatin A and sodium butyrate, suppress differentiation into osteoclasts but not into macrophages Blood 2003 101 9 3451 3459 10.1182/blood-2002-08-2622 12511413
Rahman, M. M. et al. Two histone deacetylase inhibitors, trichostatin A and sodium butyrate, suppress differentiation into osteoclasts but not into macrophages. Blood 101(9), 3451–3459 (2003).12511413 10.1182/blood-2002-08-2622
43. Duncan SH Belenguer A Holtrop G Reduced dietary intake of carbohydrates by obese subjects results in decreased concentrations of butyrate and butyrate-producing bacteria in feces Appl. Environ. Microbial. 2007 73 4 1073 1078 10.1128/AEM.02340-06
Duncan, S. H. et al. Reduced dietary intake of carbohydrates by obese subjects results in decreased concentrations of butyrate and butyrate-producing bacteria in feces. Appl. Environ. Microbial. 73(4), 1073–1078 (2007).10.1128/AEM.02340-06
44. Ghosh TS Rampelli S Jeffery IB Mediterranean diet intervention alters the gut microbiome in older people reducing frailty and improving health status: The NU-AGE 1-year dietary intervention across five European countries Gut 2020 69 7 1218 1228 10.1136/gutjnl-2019-319654 32066625
Ghosh, T. S. et al. Mediterranean diet intervention alters the gut microbiome in older people reducing frailty and improving health status: The NU-AGE 1-year dietary intervention across five European countries. Gut 69(7), 1218–1228 (2020).32066625 10.1136/gutjnl-2019-319654
45. Liu S Li G Zhu Y Analysis of gut microbiome composition, function, and phenotype in patients with osteoarthritis Front. Microbiol. 2022 13 980591 10.3389/fmicb.2022.980591 36504782
Liu, S. et al. Analysis of gut microbiome composition, function, and phenotype in patients with osteoarthritis. Front. Microbiol. 13, 980591 (2022).36504782 10.3389/fmicb.2022.980591
46. Wilson AS Koller KR Ramaboli MC Diet and the human gut microbiome: An international review Dig. Dis. Sci. 2020 65 3 723 740 10.1007/s10620-020-06112-w 32060812
Wilson, A. S. et al. Diet and the human gut microbiome: An international review. Dig. Dis. Sci. 65(3), 723–740 (2020).32060812 10.1007/s10620-020-06112-w
47. Kaakoush NO Sutterella Species, IgA-degrading bacteria in ulcerative colitis Trends Microbial. 2020 28 7 519 522 10.1016/j.tim.2020.02.018
Kaakoush, N. O. Sutterella Species, IgA-degrading bacteria in ulcerative colitis. Trends Microbial. 28(7), 519–522 (2020).10.1016/j.tim.2020.02.018
48. Zaiss MM Jones RM Schett G The gut-bone axis: How bacterial metabolites bridge the distance J. Clin. Investing. 2019 129 8 3018 3028 10.1172/JCI128521
Zaiss, M. M. et al. The gut-bone axis: How bacterial metabolites bridge the distance. J. Clin. Investing. 129(8), 3018–3028 (2019).10.1172/JCI128521
49. Boukthir S Matoussi N Belhadj A Abnormal intestinal permeability in children with autism La Tunisie medicale 2010 88 9 685 686 20812190
Boukthir, S. et al. Abnormal intestinal permeability in children with autism. La Tunisie medicale 88(9), 685–686 (2010).20812190
50. Williams BL Hornig M Parekh T Application of novel PCR-based methods for detection, quantitation, and phylogenetic characterization of Sutterella species in intestinal biopsy samples from children with autism and gastrointestinal disturbances mBio 2012 10.1128/mBio.00261-11 23131829
Williams, B. L. et al. Application of novel PCR-based methods for detection, quantitation, and phylogenetic characterization of Sutterella species in intestinal biopsy samples from children with autism and gastrointestinal disturbances. mBio10.1128/mBio.00261-11 (2012).23131829 10.1128/mBio.00261-11
51. He J Xu S Zhang B Gut microbiota and metabolite alterations associated with reduced bone mineral density or bone metabolic indexes in postmenopausal osteoporosis Aging 2020 12 9 8583 8604 10.18632/aging.103168 32392181
He, J. et al. Gut microbiota and metabolite alterations associated with reduced bone mineral density or bone metabolic indexes in postmenopausal osteoporosis. Aging 12(9), 8583–8604 (2020).32392181 10.18632/aging.103168
52. Aljuraiban GS Alfhili MA Aldhwayan MM Metagenomic shotgun sequencing reveals specific human gut microbiota associated with insulin resistance and body fat distribution in Saudi women Biomolecules 2023 13 4 640 10.3390/biom13040640 37189387
Aljuraiban, G. S. et al. Metagenomic shotgun sequencing reveals specific human gut microbiota associated with insulin resistance and body fat distribution in Saudi women. Biomolecules 13(4), 640 (2023).37189387 10.3390/biom13040640
53. Alsuhaibani MN Aljuraiban GS Aljazairy EA Dietary polyphenols in relation to gut microbiota composition in Saudi Arabian females Metabolites 2022 13 1 6 10.3390/metabo13010006 36676929
Alsuhaibani, M. N. et al. Dietary polyphenols in relation to gut microbiota composition in Saudi Arabian females. Metabolites 13(1), 6 (2022).36676929 10.3390/metabo13010006
