Skip to main content

Gut microbiota signatures and fecal metabolites in postmenopausal women with osteoporosis



Women suffer from various distress and disturbances after menopause, including osteoporosis, a risk factor associated with multiple diseases. Altered gut microbiota has been implicated in postmenopausal osteoporosis. In this study, to understand gut microbiota signatures and fecal metabolite changes in postmenopausal women with osteoporosis, 108 postmenopausal women were recruited for intestinal microbiota and fecal metabolite detection. Among these participants, 98 patients, who met the inclusion criteria, were divided into postmenopausal osteoporosis (PMO) and non-postmenopausal osteoporosis (non-PMO) groups based on bone mineral density (BMD). The compositions of gut bacteria and fungi were examined by 16 S rRNA gene sequencing and ITS sequencing, respectively. Meanwhile, fecal metabolites were analyzed using liquid chromatography coupled with mass spectrometry (LC-MS).


We found that bacterial α-diversity and β-diversity were significantly altered in PMO compared to non-PMO patients. Interestingly, fungi composition showed larger changes, and the differences in β-diversity were more significant between PMO and non-PMO patients. Metabolomics analysis revealed that fecal metabolites, such as levulinic acid, N-Acetylneuraminic acid, and the corresponding signaling pathways were also changed significantly, especially in the alpha-Linolenic acid metabolism and selenocompound metabolism. The screened differential bacteria, fungi, and metabolites closely correlated with clinical findings between these two groups, for example, the bacterial genus, Fusobacterium, the fungal genus, Devriesia, and the metabolite, L-pipecolic acid, were significantly associated with BMD.


Our findings indicated that there were remarkable changes in gut bacteria, fungi, and fecal metabolites in postmenopausal women, and such changes were notably correlated with patients’ BMD ​​and clinical findings. These correlations provide novel insights into the mechanism of PMO development, potential early diagnostic indicators, and new therapeutic approaches to improve bone health in postmenopausal women.


Osteoporosis is a metabolic bone disease that is characterized by the decreased bone mass per volume-unit. The low bone mass results in increased bone fragility and destroys microstructure, thereby reducing bone strength and increasing the risk of fractures at different sites [1, 2]. Osteoporotic fractures which are most common in elderly women (> 55 years old) and men (> 65 years old) can significantly increase bone disease-related morbidity and mortality [3]. The increased risk of osteoporosis and the following fragility fractures are the serious consequences as women grew older. About 10% of the global population suffer from osteoporosis, among which postmenopausal women over 50 years old account for 30% [4, 5]. At present, the commonly used anti-osteoporosis drugs in clinical practice have been limited by a number of factors including unobvious efficacy, long-term medication, allergic reaction, and mandibular osteonecrosis [4]. Therefore, early identification of postmenopausal women at risk for fracture and searching for safe and effective preventive intervention strategies to reduce the risk of fracture have important clinical significance.

The gut microbiome is primarily responsible for the balance and maintenance of the interaction between host and microorganisms, mainly including bacteria, fungi and viruses [6, 7]. In healthy individuals, the mutual regulation between intestinal flora and the host helps to maintain normal gastrointestinal function [7, 8]. Apart from preventing toxins from entering the peripheral circulation, a healthy gastrointestinal tract contributes to regulating the absorption of nutrients and water, and forming an intestinal barrier [9, 10]. Due to the existence of intestinal barrier, intestinal microorganisms can safely reside in the intestine. However, disruption of the intestinal ecosystem results in various digestive ailments such as ulcerative colitis and Crohn’s disease, as well as obesity, diabetes, immune system dysregulation, and osteoporosis-related metabolic diseases [9, 11]. Among gut microbiota, the proportionately low content of fungi accounts for 0.1% of total intestinal microbes. Although fungi comprise such a small percentage of intestinal flora, these microorganisms indeed affect the occurrence and development of multi-system diseases [12, 13]. However, until now no evidence had implicated a relationship between fungi and osteoporosis.

Accumulating evidence indicate[s] that the pathological process of osteoporosis is regulated by gut microbes [14]. Clinical studies have reported that the overproliferation of intestinal flora is associated with the decreased bone mineral density (BMD). Patients with enteric bacterial over-growth syndrome generally appear low bone mineral density and osteomalacia, high levels of pro-inflammatory factors such as TNF-α and IL-1, and increased activated osteoclasts [15]. Due to the important regulatory role of intestinal flora on metabolism, lack of intestinal flora at birth leads to many physiological and metabolic changes in the body, including the reduced absorption of calories, vitamins, and nutrients, and the delayed height, weight, and organ development [16]. Gut microbiota dysbiosis causes immaturity of immune, vascular, endocrine, intestinal, and nervous systems, all of which are involved in the regulation of bone mass [17, 18]. Low-dose penicillin induces a decrease in intestinal microbiota in prepubertal mice (21 days old), resulting in changes in intestinal microbiota metabolites and the abnormalities of intestinal immune [19], while antibiotic therapy-mediated intestinal microbiota depletion accelerates weight and bone growth [20]. Therefore, it follows that there is a close link between gut flora and bone loss. Importantly, there are great differences in the intestinal flora colonized in different regions and populations. The gut microbiome of postmenopausal women reveals an altered community dynamic, with co-presentations of osteoporosis and/or osteopenia diagnoses [21]. Evidence gathered on the gut microbiota-bone axis suggests Prevotella histicola was specifically able to prevent estrogen deficiency-induced bone loss [22]. A randomized controlled trial revealed that a bioavailable isoflavone and probiotic treatment can improve bone status and estrogen metabolism in postmenopausal osteopenic women [23]. Therefore, gut microbes are directly involved in the regulation of bone metabolism in postmenopausal osteoporosis. However, the specific mechanisms between gut microbiota, their fecal metabolites, and bone metabolism remains unclear.

In this study, the differential gut bacteria and fungi of intestinal microbiota and fecal metabolites were analyzed between postmenopausal osteoporosis (PMO) and non-postmenopausal osteoporosis (non-PMO) women using 16 S rRNA gene sequencing and ITS sequencing. The present study provides potential early diagnostic indicators to discriminate PMO and offers new strategies for treating osteoporosis.


The difference of intestinal microorganisms between PMO and non-PMO patients

To uncover the difference of intestinal microorganisms between PMO and non-PMO patients, we examined the abundance of gut bacteria and fungi in these two groups. Using 16 S rRNA sequencing, 82 exclusive bacteria genera were observed in non-PMO population, while 48 bacteria genera were screened in PMO cohort. Additionally, both groups shared 1,952 genera of bacteria (Fig. 1A). In terms of fungi, 84 fungi were unique to the non-PMO population, while 90 fungi were unique to the PMO population, and 221 fungi were shared between these two groups (Fig. 1B).

Fig. 1
figure 1

Venn diagram of gut microbiota composition in PMO and non-PMO patients. (A) Venn diagram demonstrating the differences in the distribution of gut bacteria between PMO and non-PMO populations. (B) Venn diagram demonstrating the differences in the distribution of gut fungi between PMO and non-PMO populations. PMO, postmenopausal osteoporosis; non-PMO, non-postmenopausal osteoporosis. Brown: the number of differential bacteria (n = 2000) or fungi (n = 2034) in PMO patients. Blue: the number of differential bacteria (n = 311) or fungi (n = 305) in non-PMO patients

Based on the in-depth analysis of gut bacteria, α-diversity index, Chao1 index, Abundance-based Coverage Estimator (ACE) and Shannon index were significantly decreased in PMO population compared to that in non-PMO patients (Fig. 2A). In comparison to non-PMO, although Simpson index and J index showed a decreasing trend in PMO population, there was no significant difference between them. However, there was no significant difference in α-diversity index, Shannon index, and Simpson index of fungi between non-PMO population and PMO population (Fig. 2B). β-diversity, which is an important index to evaluate the distribution of bacterial and fungal genera, showed a significant decrease in gut bacteria of PMO patients and a significant increase in gut fungi of PMO patients by contrast to non-PMO patients (Fig. 2C and D). Therefore, intestinal microorganisms changed greatly in PMO patients.

Fig. 2
figure 2

Analysis of intestinal microbial diversity in PMO and non-PMO patients. Box plot demonstrating α-diversity index, including Chao1 index, Abundance-based Coverage Estimator and Shannon index, Simpson index and J index in intestinal bacteria (A) and intestinal fungi (B) of PMO and non-PMO patients, repectively. Brown: PMO; Blue: non-PMO. PCoA plot was used to analyze the differential beta-diversity index in intestinal bacteria (C) and intestinal fungi (D) of PMO and non-PMO patients, repectively. *P < 0.05, **P < 0.01

Bacterial genus with differential abundance between PMO and non-PMO patients

The search for signature species between different populations contributes to exploring the distinct diagnostic factors and providing potential candidate targets for subsequent mechanism studies. The LDA Effect Size index was set as 2 between these two groups. Data indicated that Veillonella, Parabacteroides, and Harryflintia were mainly enriched in PMO population, while Veillonella, Prevotella, and Enterobacterium mainly appeared in non-PMO patients. Importantly, the abundance of these enriched bacteria differed significantly between the two groups (Fig. 3A and B). In terms of fungi, Pichia, Auricularia, and Myrothecium were observed in non-PMO patients, while Eurotium, Penicillium, and Chlorophyllum were mainly concentrated in PMO population (Fig. 3C and D). Thus, there were significant differences in genus levels between PMO and non-PMO patients.

Fig. 3
figure 3

The differential bacterial community analysis between PMO and non-PMO patients. (A and C) Significantly enriched bacterial taxa (A) and fungi (C) in the different groups as determined by LEfSe analysis (LDA sore > 2). (B and D) Boxplot of the top 10 genera that showed the difference of the microbiota composition in bacteria (B) and fungi (D) level between PMO and Non-PMO groups. Brown: PMO; Blue: non-PMO. *P < 0.05, **P < 0.01, ***P < 0.001

Analysis of intestinal microbial metabolic pathways in PMO and non-PMO patients

Metabolic pathway analysis was performed using the KEGG database. Based on KEGG database, PICRUSt analysis was applied to predict the functional profiling of microbial communities according to 16 S sequencing data. We confirmed that the enrichment of flora-related Adipocytokine signaling pathway, Amoebiasis pathway, and Ethylbenzene degradation pathway was significantly increased in PMO population compared with non-PMO population. However, the enrichment of pathways in bladder cancer, prion diseases, and bacterial invasion of epithelial cells were significantly decreased in PMO compared to non-PMO patients (Fig. 4A). COG database was further used to analyze the top 20 potential metabolic pathways. The results indicated that the pathway related to transposase and inactivated derivatives was enriched in non-PMO population, and the pathway associated with Holliday junction resolvasome and endonuclease subunit were enriched in PMO patients (Fig. 4B). Therefore, metabolic pathways were also changed significantly in PMO patients.

Fig. 4
figure 4

Metabolic pathways analysis between PMO and non-PMO patients. (A) KEGG was used to analyze the different metabolic pathways in intestinal bacteria. PICRUSt was used to calculate microbial abundance. Box plot demonstrating differential metabolomic pathways in intestinal bacteria of PMO and non-PMO women. (B) PICRUSt was used to calculate microbial abundance. COG was employed to analyze intestinal bacterial differential metabolic pathways. Box plot demonstrating differential metabolomic pathways (ID) of PMO and non-PMO women. Brown: PMO; Blue: non-PMO.

Metabolomics analysis of fecal metabolites between PMO and non-PMO patients

To further evaluate the association between metabolites and bone mineral density, we employed metabolomics to detect the distribution of metabolites in these two groups. As shown in Fig. 5A, the principal component analysis (PCA) showed a significant difference in fecal metabolite compositions between non-PMO and PMO patients (Fig. 5A). In brief, the abundances of Levulinic acid, N-Acetylneuraminic acid, Pimelic acid, Adenine, L-Lysine, Hydroxyisocaproic acid, and Azelaic acid notably raised in PMO population, while the abundance of Quinate, Kynurenic acid, Phit-val, Glycochenodeoxycholate, 1-Palmitoyl-sn-glycero-3-phosphocholine, and 1-Stearoyl-2-hydroxy-sn-glycero-3-phosphocho declined significantly in PMO patients, by contrast to that in non-PMO population (Fig. 5B). Subsequently, cluster analysis revealed that the main enrichment pathways in PMO cohort included alpha-Linolenic acid metabolism, Selenocompound metabolism, beta-Alanine metabolism and Lysine degradation, Arginine and proline metabolism, D-Glutamine and D-glutamate metabolism and Nitrogen metabolism (Fig. 5C). Possibly, the 16 S rRNA gene sequencing results showed that some components of the fecal microbiota significantly differed between the PMO and non-PMO groups; thus, we inferred that alterations in the fecal microbiota may lead to alterations in fecal metabolites.

Fig. 5
figure 5

Metabolomics analysis of fecal metabolites between PMO and non-PMO patients. (A) Unsupervised principal component analysis (PCA) score plot demonstrating clear separation between fecal metabolomic profiles of PMO and non-PMO women based on all shared metabolites. PCA was calculated using MetaboAnalyst metabolomics analysis suite and ANOVA followed by Tukey’s multiple comparisons test using PRISM. (B) Scatter diagram demonstrating differential abundance of negative ion patterns (left) and positive ion patterns (right) of metabolomics between PMO and non-PMO patients. (C) Box plot demonstrating differential metabolite-related signaling pathways (top25) between PMO and non-PMO patients. PMO, postmenopausal osteoporosis; non-PMO, non-postmenopausal osteoporosis

Association analysis of bacteria, fungi, and metabolites with clinical indicators in PMO and non-PMO patients

Next, we analyzed the correlation between intestinal bacteria/fungi/metabolites and clinical indications using Pearson correlation analysis in PMO and non-PMO patients. Data revealed that gut bacterial Fusobacterium, Parabacteroides, Anaerotruncus, Defluviitaleaceae, Acetanaerobacterium, and Leptotrichia were closely related to BMD (Fig. 6A). Gut fungal Devriesia, Montagnulaceae, and Nectriaceae were significantly correlated with BMD (Fig. 6B). Metabolites related to BMD included L-Pipecolic acid, alpha-Linolenoyl ethanolamide, D-Alanyl-D-alanine(D-Ala-D-Ala), N-Acetylmannosamine, and Serine-Valine (Fig. 6C and D). Thus, the altered intestinal bacteria/fungi/metabolites in PMO and non-PMO patients were related to BMD.

Fig. 6
figure 6

Clinical correlation analysis of bacteria, fungi, and metabolite in PMO-and non-PMO patients. (A) Clustering analysis demonstrating the correlation between intestinal bacterial and clinical indicators in PMO and non-PMO patients. (B) Clustering analysis demonstrating the correlation between the abundance of intestinal fungi and clinical indicators in PMO and non-PMO patients. (C) Clustering analysis demonstrating the correlation between metabolites and clinical indicators by mass spectrometry in the negative ion mode. (D) Clustering analysis demonstrating the correlation between metabolites and clinical indicators by mass spectrometry using positive ion-mode detection. PMO, postmenopausal osteoporosis; non-PMO, non-postmenopausal osteoporosis

Modeling analysis of PMO based on gut bacteria and fungi

Subsequently, we established an identifying model for distinguishing non-PMO and PMO patients according to the differential bacterial and fungi communities. As shown in Fig. 7A, we developed a model using characteristic bacterial Prevotella as a classification factor according to the random forest model (Fig. 7A). The results showed that the identification efficiency between PMO and non-PMO population was as high as 0.9008 (Fig. 7B). Consistently, the discriminant efficiency reached 0.8151 when fungi were used as the classification factor in the model. To verify the true validity of our model, we further collected fecal samples from a prospective cohort of patients and performed fecal flora analysis. The discriminant efficiency of the bacteria-based classification model in this cohort was 0.8962, while the discriminant efficiency of the fungi-based classification model was 0.7923 (Fig. 7C and D).

Fig. 7
figure 7

Intestinal microbiota-based model for distinguishing non-PMO patients from PMO patients. (A) Bacterial factors for intestinal bacteria-based PMO diagnosis model. (B) Line chart demonstrating the discriminant efficacy of the intestinal bacteria-based PMO diagnosis model in the exploration and validation cohorts. (C) Bacterial factors for intestinal fungi-based PMO diagnosis model. (D) Line chart demonstrating the differential efficacy of the intestinal fungi-based PMO diagnosis model in the exploration and validation cohorts. PMO, postmenopausal osteoporosis; non-PMO, non-postmenopausal osteoporosis


The intestinal tract is known as the second brain of human body, in which intestinal flora is a functional neuron that coordinates the operation of the whole-body system. Particularly, intestinal flora can participate in the regulation of postmenopausal osteoporosis [24, 25]. There is a close relationship between gut microbiome and bone turnover markers in postmenopausal women [26]. Numerous studies have explored the association between bone and intestinal flora in germ-free mice lacking intestinal flora, animal models treated with antibiotics or probiotics, and humans [27]. Although there are somewhat contradictions, the current consensus is that gut flora acts as a major regulator of BMD by influencing the immune system [28]. Supplementation with probiotics or a dietary fiber diet can regulate the distribution of intestinal flora, thereby reducing bone loss caused by estrogen loss [29, 30].

Gut microbiota changes have a significant impact on bone loss. There is growing evidence that connecting microbiome and menopause holds promise for new interventions to alleviate menopausal symptoms and for healthy ageing [31]. Lower gut microbiome diversity and a shift toward greater similarity to the male gut microbiome were confirmed in PMO patients [32]. Using 16 S rRNA gene sequencing, researchers have identified a close relationship between gut microbiota composition and osteoporosis/fracture risk in Japanese postmenopausal women [33]. However, there was no significant difference in bacterial α-diversity between the two groups. For example, there was no significant difference in bacterial α-diversity between the two groups. This is similar to an emerging report that confirmed the enrichment Lactobacillus sp. in non-PMO and the increased abundances of Peptoniphilus sp., propionic acid bacteria, and members of the Galicola genus in PMO [34]. Gut fungi accounts for a small proportion of intestinal microbes. However, it participates in the occurrence and development of multiple diseases [12, 13]. For instance, it has demonstrated that antifungals treatment can reduce liver damage in a fecal microbiome-humanized mouse model of Western diet-induced steatohepatitis [35]. Thus, fungi may play the important role in a variety of pathological processes. Our data were consistent with these findings. In the present study, species richness, diversity, and characteristic fungi appeared significantly different in PMO compared to that in non-PMO patients, indicating that there was a close correlation between fungi and osteoporosis. This finding further confirmed that fungi was an important regulatory factor among the overall biological flora in PMO patients. However, the detailed interaction between them is worth further study. In addition to that, our data proved a significant association between Fusobacterium/Parabacteroides with BMD of PMO patients. Consistently, previous studies have confirmed that the abundance of these bacteria was also altered in osteoporosis patients [36, 37]. In term of fungi, Devriesia, Montagnulaceae and Nectriaceae were associated with BMD indexes in this study. However, there is no evidence to confirm the association between these fungi species and osteoporosis. Thus, our data might provide a new target for exploring the mechanism underlying the process of PMO. Considering the positive association between Devriesia and BMD, and the negative association between Montagnulaceae and BMD, the two fungi species may play the beneficial and harmful roles during osteoporosis, respectively.

Several mechanisms underlie abnormal gut flora-mediated bone loss. Clinically, it has been found that the damaged strength in women’s spinal curvature strength is associated with a significant reduction in the gut microbiome, B cells, and T cells, which indicates a correlation between immune cell number and bone tissue characteristics [38]. Additionally, depletion of intestinal flora can lead to hyperimmune states [39]. Gut flora exerts lasting effects on the immune system either through direct contact or its metabolites. The increased bacterial abundance in the gut can raise the antigen load, leading to increased inflammatory cytokines, impaired osteoblast function, and bone loss in sickle cell disease (SCD) mice by compromising intestinal barrier of the immune system [40]. Previous study has indicated that menopause leads to increased gut permeability and inflammation, and greater gut permeability is related with more inflammation and lower bone mineral density across the menopause transition [41]. Thus, the association between gut flora and bone loss in PMO patients has the involvement of inflammation and immune regulation. In our study, significant changes in fecal metabolite abundance were observed in non-PMO and PMO populations, which may affect various metabolic pathways of the host, including alpha-Linolenic acid metabolism, etc. These metabolites and metabolic pathways have been shown to be closely related to bone metabolism and development [42, 43]. However, we also found that some novel metabolites, such as levulinic acid and pimelic acid, were significantly up-regulated in PMO population compared to non-PMO patients. Although no significant correlation was observed between these metabolites and BMD, these metabolites may still potentially influence the development of postmenopausal osteoporosis in women. Intestinal fungi are capable of developing symbiotic interactions with the hosts [12]. In brief, fungal wall components can be recognized by receptors of host cells, subsequently triggering antifungal signal transduction cascades and ultimately regulating innate and adaptive immune responses [44, 45]. Considering that intestinal bacteria-mediated the activation of immune system is closely related to the occurrence of osteoporosis, intestinal fungi may also regulate bone loss by regulating immune system.

The incidence of postmenopausal osteoporosis is about 57%. It is of great significance to effectively predict the occurrence of osteoporosis at the early stage [1, 46]. A recent study has reported 20 factors that closely affect osteoporosis. After screening with multiple feature selection methods, the differential efficiency of the random forest model is 0.921, which may be used as a practical method for the early diagnosis of postmenopausal osteoporosis [47]. Other studies have also proposed other potential predictors [48, 49]. In this study, we focused on the early diagnostic value of intestinal flora in PMO. Combined with the stochastic forest model, we confirmed that intestinal microecology including bacteria and fungi could act as a better identifier for predicting PMO. Importantly, although the efficacy of the validation cohort was reduced compared to the exploration cohort, the AUC was still higher than 0.8 for both bacteria and fungi.


Collectively, our results suggested that postmenopausal women have significant changes in gut bacteria, fungi, and metabolites that are significantly correlated with BMD values and clinical features of PMO patients. This correlation provides potential directions for exploring the mechanism of PMO development and provides potential early diagnostic indicators for PMO. This study may provide novel interventions to improve the level of bone health in postmenopausal women.

Materials and methods

Study cohort

This study was a designed case-control study. This study was approved by the Ethics Committee of Zhongshan Hospital of Xiamen University (No. 201,808). Non-PMO subjects (n = 58) and newly diagnosed PMO patients (n = 40) who were admitted to Zhongshan Hospital Xiamen University from September 2021 up to now were recruited. Clinical information was presented in Table 1. For modeling analysis, a prospective cohort of non-PMO subjects (n = 10) and PMO patients (n = 13) were collected from Xinyu People’s Hospital (Jiangxi Province, CN). The trial number was ChiCTR1900027187. Clinical information of enrolled patients (modeling analysis) was presented in Table 2. Written informed consent was obtained from PMO patients and healthy people. The inclusion criteria were as follows: postmenopausal women were aged 50–70 years and were diagnosed with osteoporosis; All cohorts participated voluntarily in the study and signed the informed consent. The exclusion criteria were as below: people with organ dysfunction, neurological diseases (such as Parkinson’s syndrome, dementia, stroke, etc.), or rheumatic immune diseases were excluded; Patients with gastrointestinal diseases such as total parenteral nutrition, inflammatory bowel disease, and gastrointestinal surgery were also excluded from this study; patients with other serious diseases, such as malignant tumors or infectious diseases were not suitable for this study; patients who used antibiotics, other microecological preparations and gastrointestinal motility drugs that could affect intestinal flora within 30 days prior to enrollment. Patients who were taking or recently using Chinese and Western medicines such as calcitonin and zoledronic acid that could affect bone metabolism; patients with secondary osteoporosis were excluded; persons with mental or legal disabilities, additional patients unsuitable for inclusion and patients who are participating in other clinical trials were not suitable for the present research. The control group was postmenopausal non-osteoporosis patients (non-PMO). According to previous literature, the incidence of osteoporosis in postmenopausal women is 57%. Based on these, the effect size D1, α-value and β-value of difference in microflora between PMO and non-PMO populations were set as 0.57, 0.05, and 0.05, respectively. The expulsion rate was set as 5%. The sample size in the two groups was defined as 38 according to SPASS 15 software. Due to the 5% expulsion rate, 40 patients in each group were enrolled in this study.

Table 1 The clinicopathological factors of non-PMO (n = 58) patients and PMO (n = 40)
Table 2 The clinicopathological factors of non-PMO (n = 10) patients and PMO (n = 13)

Clinical data

Basic information such as age, height, and weight of all subjects were recorded, and body mass index (BMI) was calculated. Blood samples were collected for serological testing in the morning after more than 6-hour fasting. Serum levels of 25 (OH) D, estradiol (E2), osteocalcin (OC), C-terminal peptide of type I collagen (CTX-I), N-terminal propeptide of type 1 procollagen (P1NP), and parathyroid hormone (PTH) were measured using the Roche Diagnostics GmbH Electrochemical Luminescence System (Roche Diagnostics GmbH, Germany). Bone mineral density (BMD) in the lumbar spine (LS: L1-4) and total hip joint (femoral neck), trochanteric and intertrochanteric areas were measured using a diurnal calibrated Hologic 4500 dual-energy X-ray absorptiometry scanner (Lunar Expert 1313, Lunar Corp, USA).

Fecal sample collection, DNA extraction, and 16 S rRNA gene sequencing

Stool samples from each volunteer were immediately stored in a -80 ℃ freezer after collection. After the stool samples were thawed and homogenized, total DNA was extracted from each sample (0.25 g) using the QIAamp Rapid DNA Stool Mini Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. The concentration and purity of extracted DNA were measured by Multiskan™ GO full-wavelength enzyme marker (Thermo Fisher Scientific, US). DNA integrity was tested by agarose-gel electrophoresis. Bacterial and fungal communities were amplified by targeting the V4 region of 16 S rRNA gene and ITS2 fragment, respectively. The forward primer sequence of 16 S was 5’-GTGCCAGCMGCCGCGGTAA-3’, and the reverse primer sequence was 5’-GGACTACNVGGGTWTCTAAT-3’. The forward primer sequence of ITS2 was 5’-GCATCGATGAAGAACGCAGC-3’, and reverse primer sequence was 5’-TCCTCCGCTTATTGATATGC-3’. Polymerase chain reaction (PCR) products were purified and evaluated using Qubit 3.0 (Thermo Fisher Scientific, US). The PCR products were purified with Qiagen Gel Extraction Kit(Qiagen, Germany). Sequencing libraries were generated usingTruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, USA) following manufacturer’s recommendations and index codes were added. The library quality was assessed on the Qubit@ 3.0 Fluorometer (Thermo Fisher Scientific, US) and Agilent Bioanalyzer 2100 system. At last, the library was sequenced on an Illumina NovaSeq platform and 250 bp paired-end reads were generated. Wilcoxon rank-sum test, that was set at a Monte-Carlo significance level α = 0.05 to calculate LDA scores, was used to detect features with different abundance levels between assigned taxa based on a normalized relative abundance matrix. All tests were performed using 999 permutations.

Bioinformatics analysis

Paired-end reads were assembled using flash software. The primer sequences and the lower readings were removed. Checks for chimeric sequences and OTU clustering are performed using clean reads. All reads were demultiplexed into a single file, clustered with 97% similarity, and then the UNITE UCHIME reference dataset (version 7) was checked for Chimera examination using UCHIME in reference mode. The representative sequence was generated; the monomers were removed; and the final OTU table was created. Representative sequences of OTU were compared on the UNITE ITS database and classified by RDP classifier.


Twenty-five milligrams of fecal samples were mixed with 500 µL of extraction solution (methanol: acetonitrile: water = 2:2:1, internal standard mixture was labeled with isotope). Then, the samples were homogenized at 35 Hz for 4 min, followed by 5 min ultrasound in an ice water bath. The homogenization and ultrasound cycle were repeated for 3 times. After incubation at -40℃ for 1 h, the mixture was centrifuged at 4℃ at 12,000 rpm for 15 min, and the supernatant was transferred to fresh glass vials for subsequent analysis. Fecal metabolomics was measured using the Vanquish (Thermo Fisher Scientific) ultra-high performance liquid chromatograph. The target compounds were separated by Waters ACQUITY UPLC BEH Amide (2.1 mm × 100 mm, 1.7 μm) liquid chromatography column. Sample plate temperature maintained at 4℃ and sample volume was set as 2 µL. Mass spectrometry data was collected by the Orbitrap Exploris 120 mass spectrometer. Detailed parameters are as follows: sheath gas flow, 50 Arb; auxiliary gas flow, 15 Arb; capillary temperature, 320℃; full ms resolution, 60,000; MS/MS resolution,15,000; NCE mode collision energy, 10/30/60; spray voltage, 3.8 kV (positive) or -3.4 kV (negative).

Functional analysis based on bacterial taxonomy

The unobserved State Reconstruction Community Phylogenetic Survey (PICRUSt) was used to predict metagenomic functional content. The sequencing data of 16 S was used to predict the presence of genes. First, the reference set GreenGenes database was read and a closed reference OTU table was constructed using QIIME software. The generated OTU table is normalized by copy number. Metagenomes were predicted using ANOVA was used for statistical difference analysis. The results were visualized using a custom R script based on ggplot2.

Machine learning models

Random Forest (R package and caret) models were trained on data from multiple omics analyses, including 16 S and ITS, respectively. The models were used to test whether data based on bacteria or fungi could predict osteoporosis. Firstly, the importance of common genera was ranked according to their average decline in accuracy. Next, we performed stepwise feature selection using a five-fold cross-validation approach to avoid over-fitting and over-optimistic estimates. This method is used to select and predict microbial characteristics and eliminate non-information characteristics. The area under ROC curve (AUC) was calculated to evaluate the differentiation of characteristic OTU.

Statistical analysis and visualization

The estimation of α-diversity was based on a uniform and sparse OTU abundance matrix. Significant differences between α-diversity were examined using the nonparametric Kruskal-Wallis test and Benjamini Hochbery correction. The β-diversity that could estimate differences in community structure between samples was measured using the Bray-Curtis distance based on a uniform sparse OTU abundance table. Statistical differences in β-diversity measure were determined using R-package. OTU was calculated and visualized using the VennDiagram in R-package. Taxonomic abundance was measured and plotted using ggplot2. Different taxa abundance in different populations was monitored using LEfSe analysis. A genus-based index analysis was performed using R packaging to label species from genus information. Finally, a custom R script based on ggplot2 was used to visualize the results and the results were analyzed using R v3.4.1.

Data availability

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



postmenopausal osteoporosis


non-postmenopausal osteoporosis


bone mineral density


liquid chromatography coupled with mass spectrometry.


  1. Muñoz M, Robinson K, Shibli-Rahhal A. Bone health and osteoporosis Prevention and Treatment. Clin Obstet Gynecol. 2020;63(4):770–87.

    Article  PubMed  Google Scholar 

  2. Rachner TD, Khosla S, Hofbauer LC. Osteoporosis: now and the future. Lancet (London England). 2011;377(9773):1276–87.

    Article  CAS  PubMed  Google Scholar 

  3. Compston JE, McClung MR, Leslie WD, Osteoporosis. Lancet (London England). 2019;393(10169):364–76.

    Article  CAS  PubMed  Google Scholar 

  4. Black DM, Rosen CJ. Clinical practice. Postmenopausal osteoporosis. N Engl J Med. 2016;374(3):254–62.

    Article  CAS  PubMed  Google Scholar 

  5. Yu F, Xia W. The epidemiology of osteoporosis, associated fragility fractures, and management gap in China. Archives of osteoporosis. 2019;14(1):32.

    Article  PubMed  Google Scholar 

  6. Helmink BA, Khan MAW, Hermann A, Gopalakrishnan V, Wargo JA. The microbiome, cancer, and cancer therapy. Nat Med. 2019;25(3):377–88.

    Article  CAS  PubMed  Google Scholar 

  7. De Luca F, Shoenfeld Y. The microbiome in autoimmune diseases. Clin Exp Immunol. 2019;195(1):74–85.

    Article  PubMed  Google Scholar 

  8. DeGruttola AK, Low D, Mizoguchi A, Mizoguchi E. Current understanding of Dysbiosis in Disease in Human and Animal Models. Inflamm Bowel Dis. 2016;22(5):1137–50.

    Article  PubMed  Google Scholar 

  9. Takiishi T, Fenero CIM, Câmara NO. Intestinal barrier and gut microbiota: shaping our immune responses throughout life. Tissue barriers. 2017;5(4):e1373208.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Peterson LW, Artis D. Intestinal epithelial cells: regulators of barrier function and immune homeostasis. Nat Rev Immunol. 2014;14(3):141–53.

    Article  CAS  PubMed  Google Scholar 

  11. Heintz-Buschart A, Wilmes P. Human gut microbiome: function matters. Trends Microbiol. 2018;26(7):563–74.

    Article  CAS  PubMed  Google Scholar 

  12. Zhang F, Aschenbrenner D, Yoo JY, Zuo T. The gut mycobiome in health, disease, and clinical applications in association with the gut bacterial microbiome assembly. Lancet Microbe. 2022;3(12):e969–83.

    Article  CAS  PubMed  Google Scholar 

  13. Ost KS, O’Meara TR, Stephens WZ, Chiaro T, Zhou H, Penman J, Bell R, Catanzaro JR, Song D, Singh S, Call DH, Hwang-Wong E, Hanson KE, Valentine JF, Christensen KA, O’Connell RM, Cormack B, Ibrahim AS, Palm NW, Noble SM. Round, J. L. Adaptive immunity induces mutualism between commensal eukaryotes. Nature. 2021;596(7870):114–118.

  14. Ding K, Hua F, Ding W. Gut microbiome and osteoporosis. Aging Dis. 2020;11(2):438–47.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Stotzer PO, Johansson C, Mellström D, Lindstedt G, Kilander AF. Bone mineral density in patients with small intestinal bacterial overgrowth. Hepatogastroenterology. 2003;50(53):1415–8.

    PubMed  Google Scholar 

  16. Wallace JG, Gohir W, Sloboda DM. The impact of early life gut colonization on metabolic and obesogenic outcomes: what have animal models shown us? J Dev origins health disease. 2016;7(1):15–24.

    Article  CAS  Google Scholar 

  17. Kennedy PJ, Cryan JF, Dinan TG, Clarke G. Kynurenine pathway metabolism and the microbiota-gut-brain axis. Neuropharmacology. 2017;112(Pt B):399–412.

    Article  CAS  PubMed  Google Scholar 

  18. Deng Y, Zhou M, Wang J, Yao J, Yu J, Liu W, Wu L, Wang J, Gao R. Involvement of the microbiota-gut-brain axis in chronic restraint stress: disturbances of the kynurenine metabolic pathway in both the gut and brain. Gut Microbes. 2021;13(1):1–16.

    Article  PubMed  Google Scholar 

  19. Cox LM, Yamanishi S, Sohn J, Alekseyenko AV, Leung JM, Cho I, Kim SG, Li H, Gao Z, Mahana D, Zárate Rodriguez JG, Rogers AB, Robine N, Loke P, Blaser MJ. Altering the intestinal microbiota during a critical developmental window has lasting metabolic consequences. Cell. 2014;158(4):705–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Nobel YR, Cox LM, Kirigin FF, Bokulich NA, Yamanishi S, Teitler I, Chung J, Sohn J, Barber CM, Goldfarb DS, Raju K, Abubucker S, Zhou Y, Ruiz VE, Li H, Mitreva M, Alekseyenko AV, Weinstock GM, Sodergren E, Blaser MJ. Metabolic and metagenomic outcomes from early-life pulsed antibiotic treatment. Nat Commun. 2015;6:7486.

    Article  PubMed  Google Scholar 

  21. Rettedal EA, Ilesanmi-Oyelere BL, Roy NC, Coad J, Kruger MC. The gut microbiome is altered in Postmenopausal Women with osteoporosis and Osteopenia. JBMR Plus. 2021;5(3):e10452.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Wang Z, Chen K, Wu C, Chen J, Pan H, Liu Y, Wu P, Yuan J, Huang F, Lang J, Du J, Xu J, Jin K, Chen L. 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–13.

    Article  PubMed  Google Scholar 

  23. Lambert MNT, Thybo CB, Lykkeboe S, Rasmussen LM, Frette X, Christensen LP, Jeppesen PB. Combined bioavailable isoflavones and probiotics improve bone status and estrogen metabolism in postmenopausal osteopenic women: a randomized controlled trial. Am J Clin Nutr. 2017;106(3):909–20.

    Article  CAS  PubMed  Google Scholar 

  24. He J, Xu S, Zhang B, Xiao C, Chen Z, Si F, Fu J, Lin X, Zheng G, Yu G, Chen J. Gut microbiota and metabolite alterations associated with reduced bone mineral density or bone metabolic indexes in postmenopausal osteoporosis. Aging. 2020;12(9):8583–604.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Lorenzo J. From the gut to bone: connecting the gut microbiota with Th17 T lymphocytes and postmenopausal osteoporosis. J Clin Investig. 2021;131(5).

  26. Chen L, Yan S, Yang M, Yu F, Wang J, Wang X, Xu H, Shi J, Pan L, Zeng Y, Li S, Li L, You L, Peng Y. The gut microbiome is associated with bone turnover markers in postmenopausal women. Am J Transl Res. 2021;13(11):12601–13.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Locantore P, Del Gatto V, Gelli S, Paragliola RM. Pontecorvi, A. The interplay between Immune System and Microbiota in osteoporosis. Mediat Inflamm. 2020;2020:3686749.

    Article  Google Scholar 

  28. Xu Q, Li D, Chen J, Yang J, Yan J, Xia Y, Zhang F, Wang X, Cao H. Crosstalk between the gut microbiota and postmenopausal osteoporosis: mechanisms and applications. Int Immunopharmacol. 2022;110:108998.

    Article  CAS  PubMed  Google Scholar 

  29. Huidrom S, Beg MA, Masood T. Post-menopausal osteoporosis and Probiotics. Curr Drug Targets. 2021;22(7):816–22.

    Article  CAS  PubMed  Google Scholar 

  30. Rizzoli R, Biver E. Are Probiotics the New Calcium and vitamin D for Bone Health? Curr Osteoporos Rep. 2020;18(3):273–84.

    Article  PubMed  Google Scholar 

  31. Laniewski P, Herbst-Kralovetz MM. Connecting microbiome and menopause for healthy ageing. Nat Microbiol. 2022;7(3):354–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Peters BA, Santoro N, Kaplan RC, Qi Q. Spotlight on the gut Microbiome in Menopause: current insights. Int J Womens Health. 2022;14:1059–72.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Ozaki D, Kubota R, Maeno T, Abdelhakim M, Hitosugi N. Association between gut microbiota, bone metabolism, and fracture risk in postmenopausal japanese women. Osteoporos Int. 2021;32(1):145–56.

    Article  CAS  PubMed  Google Scholar 

  34. Yang X, Chang T, Yuan Q, Wei W, Wang P, Song X, Yuan H. Changes in the composition of gut and vaginal microbiota in patients with postmenopausal osteoporosis. Front Immunol. 2022;13:930244.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Demir M, Lang S, Hartmann P, Duan Y, Martin A, Miyamoto Y, Bondareva M, Zhang X, Wang Y, Kasper P, Bang C, Roderburg C, Tacke F, Steffen HM, Goeser T, Kruglov A, Eckmann L, Starkel P, Fouts DE. Schnabl, B. The fecal mycobiome in non-alcoholic fatty liver disease. J Hepatol. 2022;76(4):788–99.

    Article  CAS  PubMed  Google Scholar 

  36. Hernandez-Vigueras S, Martinez-Garriga B, Sanchez MC, Sanz M, Estrugo-Devesa A, Vinuesa T, Lopez-Lopez J, Vinas M. Oral microbiota, Periodontal Status, and osteoporosis in postmenopausal females. J Periodontol. 2016;87(2):124–33.

    Article  PubMed  Google Scholar 

  37. Wang J, Wang Y, Gao W, Wang B, Zhao H, Zeng Y, Ji Y, Hao D. Diversity analysis of gut microbiota in osteoporosis and osteopenia patients. PeerJ. 2017;5:e3450.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Guss JD, Horsfield MW, Fontenele FF, Sandoval TN, Luna M, Apoorva F, Lima SF, Bicalho RC, Singh A, Ley RE, van der Meulen MC, Goldring SR, Hernandez CJ. Alterations to the gut Microbiome Impair Bone Strength and tissue material Properties. J bone mineral research: official J Am Soc Bone Mineral Res. 2017;32(6):1343–53.

    Article  CAS  Google Scholar 

  39. Hathaway-Schrader JD, Steinkamp HM, Chavez MB, Poulides NA, Kirkpatrick JE, Chew ME, Huang E, Alekseyenko AV, Aguirre JI, Novince CM. Antibiotic perturbation of Gut Microbiota Dysregulates Osteoimmune Cross talk in postpubertal skeletal development. Am J Pathol. 2019;189(2):370–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Tavakoli S, Xiao L. Depletion of intestinal Microbiome partially rescues bone loss in Sickle Cell Disease Male mice. Sci Rep. 2019;9(1):8659.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Shieh A, Epeldegui M, Karlamangla AS, Greendale GA. Gut permeability, inflammation, and bone density across the menopause transition. JCI Insight. 2020;5(2).

  42. Kim Y, Ilich JZ. Implications of dietary α-linolenic acid in bone health. Nutr (Burbank Los Angeles Cty Calif). 2011;27(11–12):1101–7.

    Article  CAS  Google Scholar 

  43. Song J, Jing Z, Hu W, Yu J, Cui X. α-Linolenic acid inhibits receptor activator of NF-κB ligand Induced (RANKL-Induced) Osteoclastogenesis and prevents inflammatory bone loss via downregulation of Nuclear Factor-KappaB-Inducible nitric oxide synthases (NF-κB-iNOS) signaling pathways. Med Sci monitor: Int Med J experimental Clin Res. 2017;23:5056–69.

    Article  Google Scholar 

  44. Richard ML, Sokol H. The gut mycobiota: insights into analysis, environmental interactions and role in gastrointestinal diseases. Nat Rev Gastroenterol Hepatol. 2019;16(6):331–45.

    PubMed  Google Scholar 

  45. Luisa Gil M, Murciano C, Yanez A, Gozalbo D. Role of toll-like receptors in systemic Candida albicans infections. Front Biosci (Landmark Ed). 2016;21(2):278–302.

    Article  PubMed  Google Scholar 

  46. ACOG Clinical Practice Guideline No. 2: management of postmenopausal osteoporosis: correction. Obstet Gynecol. 2022;140(1):138.

    Article  Google Scholar 

  47. Kwon Y, Lee J, Park JH, Kim YM, Kim SH, Won YJ, Kim HY. Osteoporosis pre-screening using Ensemble Machine Learning in Postmenopausal Korean Women. Healthc (Basel Switzerland). 2022;10(6).

  48. Momeni M, Asadzadeh M, Mowla K, Hanafi MG, Gharibvand MM, Sahraeizadeh A. Sensitivity and specificity assessment of DWI and ADC for the diagnosis of osteoporosis in postmenopausal patients. Radiol Med. 2020;125(1):68–74.

    Article  PubMed  Google Scholar 

  49. Migliorini F, Maffulli N, Spiezia F, Peretti GM, Tingart M, Giorgino R. Potential of biomarkers during pharmacological therapy setting for postmenopausal osteoporosis: a systematic review. J Orthop Surg Res. 2021;16(1):351.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


Not applicable.


This work was supported by National Natural Science Foundation of China (81972034, 92068104 to RX and 82002262 to NL), Natural Science Foundation of Fujian Province (2022J06003 to RX) and Project of Xiamen Cell Therapy Research (Grant No. 3502Z20214001).

Author information

Authors and Affiliations



WH, LJ, CJ, and XR planned and designed the research. WH, LJ, WZX, ZYY, CM, SBH, CBL, CN, GH, LN, CJ, and XR conducted the experiments and analyzed the data. WH and XR wrote the manuscript. All authors were involved in revising the manuscript critically. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Jian Chen or Ren Xu.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Zhongshan Hospital of Xiamen University (No. 201808). Written informed consent was obtained from the patients.

Consent for publication

All participants provided written informed consent to participate in the studies and for publication of the findings of the study. All authors give consent to publish the research.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, H., Liu, J., Wu, Z. et al. Gut microbiota signatures and fecal metabolites in postmenopausal women with osteoporosis. Gut Pathog 15, 33 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: