- 1Department of Spine Surgery, Honghui Hospital, Xi’an Jiaotong University, Xi’an, China
- 2Institute of Medical Research, Northwestern Polytechnical University, Xi’an, China
- 3Department of Osteoporosis, Honghui Hospital, Xi’an Jiaotong University, Xi’an, China
- 4School of Electronics and Information, Northwestern Polytechnical University, Xi’an, China
Osteoporosis (OP) is a metabolic bone disorder characterized by low bone mass and deterioration of micro-architectural bone tissue. The most common type of OP is postmenopausal osteoporosis (PMOP), with fragility fractures becoming a global burden for women. Recently, the gut microbiota has been connected to bone metabolism. The aim of this study was to characterize the gut microbiota signatures in PMOP patients and controls. Fecal samples from 21 PMOP patients and 37 controls were collected and analyzed using amplicon sequencing of the V3-V4 regions of the 16S rRNA gene. The bone mineral density (BMD) measurement and laboratory biochemical test were performed on all participants. Two feature selection algorithms, maximal information coefficient (MIC) and XGBoost, were employed to identify the PMOP-related microbial features. Results showed that the composition of gut microbiota changed in PMOP patients, and microbial abundances were more correlated with total hip BMD/T-score than lumbar spine BMD/T-score. Using the MIC and XGBoost methods, we identified a set of PMOP-related microbes; a logistic regression model revealed that two microbial markers (Fusobacteria and Lactobacillaceae) had significant abilities in disease classification between the PMOP and control groups. Taken together, the findings of this study provide new insights into the etiology of OP/PMOP, as well as modulating gut microbiota as a therapeutic target in the diseases. We also highlight the application of feature selection approaches in biological data mining and data analysis, which may improve the research in medical and life sciences.
Introduction
Osteoporosis (OP) is the most common metabolic bone disease and is characterized by low bone mass and microarchitectural deterioration of bone tissue. Postmenopausal osteoporosis (PMOP), which results from estrogen deficiency, is the most common type of osteoporosis (Eastell et al., 2016). The human gut microbiota has been regarded as a key mediator of osteoporosis and osteogenesis (Seely et al., 2021). By regulating host metabolism, immune function and hormone secretion, the gut microbiota can affect bone metabolism (Li et al., 2019). For example, the gut microbiota can regulate bone mass and improve osteoporosis by inhibiting osteoclast proliferation and differentiation, inducing apoptosis, reducing bone resorption, or promoting osteoblast proliferation and maturation (Ding et al., 2020). By influencing multiple factors such as short-chain fatty acids (SCFAs), estrogen, immune factors and vitamin D, gut microbiome play important roles on calcium balance (Wang J et al., 2022). Supplementation of antibiotic-treated mice with microbial metabolism products SCFAs could induce insulin-like growth factor 1 (IGF-1) levels, and promote bone formation and growth (Yan et al., 2016). Thus, researchers have proposed that the gut microbiota is an overlooked factor that plays significant roles in osteoporosis (Hao et al., 2019). This evidence invokes the perspective of the microbiota-gut-bone axis and supports the contention that the gut microbiota may be a novel therapeutic target in the treatment of osteoporosis and fracture prevention (Li et al., 2021; He and Chen, 2022).
Recently, a growing number of studies in OP and PMOP research have linked bone loss to changes in gut microbiota. Studies have indicated that primary osteoporosis is related to changes in the gut microbiome, particularly the enriched genera Dialister and Faecalibacterium (Xu et al., 2020). Cinnamic acid suppresses bone loss via induction of osteoblast differentiation with the alteration of gut microbiota (Hong et al., 2022). Studies have also found that external interventions such as dietary changes can alter the gut microbiota composition and improve bone health (Yan et al., 2022). Overall, multiple factors, including hereditary, dietary and physical factors could alter the gut microbiota composition and further regulate bone metabolism (Yan et al., 2022).
Our previous studies have preliminarily investigated changes in the diversity and composition of gut microbiota in elderly Chinese osteoporosis patients (Wang et al., 2017; Wang Y et al., 2022). We found that gut microbiota changes were more highly correlated with bone mineral density (BMD) alterations in elderly females than in elderly males (Wang Y et al., 2022), which implies that the gut microbiota is a potential target for the management of PMOP patients. Here, we conducted a further analysis of the gut microbiota in postmenopausal women. It is worth mentioning that feature selection-based machine learning approaches (detailed information regarding feature selection approaches are described in the Methods) were employed to identify PMOP-related microbial features, which may provide additional clues for biological data analysis.
Materials and methods
Study population
A total of 58 postmenopausal female participants, including 37 healthy controls and 21 PMOP patients were recruited for this study. Specifically, female participants with a T-score > −2.5 for the BMD measurement were considered as healthy controls, whereas female participants with a T-score < −2.5 for the BMD measurement were diagnosed with PMOP. We excluded participants or patients with other malignancies, including chronic heart disease, liver disease, kidney disease, diabetes, and diseases related to secondary OP (hyperthyroidism, steroid abuse, Cushing’s syndrome, hyperparathyroidism, etc.,). In this study, none of the participants had a history of fractures. The participant characteristics are shown in Table 1. Our study was approved by the Institutional Review Board of Honghui Hospital, Xi’an Jiaotong University and was a part of the project “Diversity Analysis for Intestinal Flora in Patients with Primary Osteoporosis” registered at 1 as #ChiCTR-1,800,019,048#. Written informed consent was obtained from each participant.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1113174/fmicb-14-1113174-HTML/image_m/fmicb-14-1113174-t001.jpg)
Table 1. Characteristics of the 58 participants and the statistical results between PMOP and control groups.
Biochemical tests
Fasting venous blood samples were collected from the 58 participants for laboratory biochemical analysis, including serum calcium (Ca), phosphorus (P), 25-hydroxyvitamin D (vitamin D), alkaline phosphatase (ALP), C-terminal telopeptide of type I collagen (b-CTX) and procollagen type 1 N-peptide (P1NP). A 5 ml vacuum blood collection tube was used to collect fasting venous blood, which was left at room temperature for 45 (±10) min, and the serum was collected for detection after centrifugation (not less than 1.5 mL serum). Specifically, the serum contents of ALP were detected by using a Roche Cobas c501 automatic biochemical analyzer and the corresponding reagents (Roche, United States); serum contents of Ca, P, Vitamin D, b-CTX and P1NP were detected by using a Roche Cobas e601 electrochemiluminescence instrument and the corresponding reagents (Roche, United States).
BMD measurement
All 58 participants underwent BMD measurement with a HOLOGIC Discovery Dual Energy X-ray Absorptiometry (DXA) scanner (HOLOGIC, United States). Subjects meeting any one of the following criteria were not allowed to undergo BMD measurement: (1) patients who had taken orally administered drugs within the previous 2–6 d that could affect image development, (2) patients who were physically weak and unable to lie on their backs or lie flat for 5 min, (3) patients with metallic implants or severe deformities of the spine, or (4) patients who had radioisotope examinations within the previous 3 days. For all participants, both lumbar spine and total hip BMD measurements were performed. We collected information pertaining to the participant BMD values (g/cm2) and T-scores of the lumbar spine and total hip. The T-score reference ranges were calculated with data from a healthy Asian population provided by the bone densitometry equipment manufacturer. Based on the diagnostic criteria, postmenopausal females with T-scores ≤ −2.5 at any site were diagnosed as having PMOP, −2.5 < T-scores < −1.0 were diagnosed as having osteopenia (ON), and T-scores ≥ −1.0 indicated normal bone mass. In this study, we divided the subjects into two major groups according to the diagnostic criteria: the PMOP group (T-scores ≤ −2.5 at any site) and the control group (T-scores > −2.5 at any site).
Fecal sample collection and bacterial DNA extraction
We collected 58 fecal samples from the participants for gut microbiota analyses. Before the collection of fecal samples, the subjects were required to avoid the following medicines for 1 month: antibiotics, probiotic, steroids, nonsteroidal anti-inflammatory drugs, microecological adjustment preparations, immunosuppressants, and any Chinese herbal medicine preparations. Participants were required to collect a 3–5 mL fresh stool sample between 6:00 and 8:00 in the morning on the day of fecal sample collection using a disposable sterile stool collection tube. The collected sample collection tube was placed in an ice box, and the researcher transported the sample to the laboratory within 2 h after collection and stored it at −80°C until further use. The QIAamp Fast DNA Stool Mini Kit (Qiagen, Germany) was used to extract the microbial genome from fecal samples following the manufacturer’s instructions. DNA concentration and purity were monitored on 1% agarose gels. Based on the concentration, the DNA was diluted to 1 ng/μL using sterile water.
16S rRNA polymerase chain reaction (PCR)
The TransGen AP221-02 Kit (TransGen, China) was used to amplify the V3-V4 region of the microbial 16S RNA gene. The polymerase chain reaction (PCR) primers of the 16S rRNA V3-V4 region were as follows: 338F 5’-ACTCCTACGGGAGGCAGCAG-3′ and 806R 5’GGACTACHVGGGTWTCTAAT-3′(Mori et al., 2014). The same volume of 1× loading buffer (containing SYBR green) was mixed with the PCR products, and electrophoresis was performed on a 2% agarose gel for detection. The PCR products were mixed in equidensity ratios. Next, the PCR products were purified with a GeneJET Gel Extraction Kit (Thermo Scientific, United States).
Library preparation and Illumina sequencing
The sequencing libraries were generated using the NEBNext Ultra DNA Library Prep Kit for Illumina (New England Biolabs, United States) following the manufacturer’s recommendations. The library quality was assessed on a Qubit@ 2.0 Fluorometer (Thermo Scientific, United States). Finally, the library was sequenced on the Illumina NovaSeq platform (Illumina, United States). The 16S rRNA amplicon sequencing project for 58 participants is available on the NCBI-BioProject site under accession PRJNA565497, and the raw data have been deposited in the NCBI-SRA database (accessions SRX6849522-SRX6849524, SRX6849535, SRX6849539-SRX6849588, SRX6849591, SRX6849602, SRX6849613, and SRX6849624).
Data analysis for 16S rRNA amplicon sequencing
Sequence analysis was performed with UPARSE software (UPARSE v7.0.1001, http://drive5.com/uparse/) (Edgar, 2013). Sequences with ≥97% similarity were assigned to the same operational taxonomic units (OTUs). Representative sequences for each OTU were screened for further annotation. For each representative sequence, the Silva Database (Version 138, https://www.arb-silva.de/) (Quast et al., 2013) was used to annotate taxonomic information based on the Mothur Bayesian classifier algorithm. To study the phylogenetic relationships of different OTUs and the difference in the dominant species among samples (groups), multiple sequence alignment was conducted using MUSCLE software (Version 3.8.31, http://www.drive5.com/muscle/). OTU abundance information was normalized using a standard sequence number corresponding to the sample with the fewest sequences. Alpha diversity was used to analyze the complexity of species diversity for a sample through several indices, including observed species, Shannon, Simpson, Chao1 and Ace. All of these indices were calculated with QIIME (Version 1.7.0). Specifically, the Shannon and Simpson indices were used to identify community diversity; the Chao1 and Ace indices were selected to identify community richness. For beta diversity analysis, we calculated both weighted and unweighted UniFrac distances using QIIME for principal coordinate analysis (PCoA).
Using XGBoost and SHAP methods to identify microbial biomarkers for PMOP
To obtain gut microbes (as features) that have the most significant effect on PMOP, the methods of maximal information coefficient (MIC) (Reshef et al., 2011; Kinney and Atwal, 2014) and XGBoost (Chen et al., 2015; Chen and Guestrin, 2016) explained with Shapley Additive exPlanations (SHAP) (Lundberg and Lee, 2017; Lundberg et al., 2020; Parsa et al., 2020) were employed for feature selection, and the intersection of the two sets of results was selected as the final microbial feature.
Maximal information coefficient
Reshef, etc., proposed a novel correlation measurement algorithm designated MIC based on mutual information (MI). Compared with MI, MIC is more equitable for high-dimensional datasets and has a wider range of applications in feature selection. For a dataset , the number of samples is , and are positive integers greater than 1. The MIC can be computed by:
where , and is the MI of in the grid of . The correlations between features and classifications can be measured by MIC, and those features with larger MIC values were considered potential signatures for the identification of PMOP. The MIC code is available on https://github.com/minepy/minepy. In our study, we performed the MATLAB version to obtain the associations between any two variables, and the parameter of alpha for function “mine” was set to 0.6.
XGBoost explained with SHAP
XGBoost is an efficient implementation of gradient-boosted decision trees that uses a sequence of decision trees to improve the model and build a stronger learner (Parsa et al., 2020). The objective function of XGBoost is as follows:
is composed of two parts: a convex loss function and a regularization term . The loss function measures the difference between the prediction and the target , and the regular function penalizes the complexity of the model (Chen and Guestrin, 2016). The parameters of and are set manually; is the vector formed by the values of all leaf nodes of the decision tree, and is the number of leaf nodes.
Generally, models using XGBoost are considered black boxes due to their complexity, and the interpretability of the XBGoost model is critical to the prediction accuracy for many applications such as bioinformatics and medicine. In 2017, Lundberg et al. proposed Shapley Additive exPlanations (SHAP) (Lundberg and Lee, 2017; Lundberg et al., 2020), a unified framework for interpreting predictions, for interpreting the output of those black boxes, and the importance of each feature can be computed using the TreeExplainer of SHAP with a high speed (Lundberg et al., 2020). In our study, we used TreeExplainer to obtain the most valuable factors with the greatest impact on PMOP. The versions of Shap and XGBoost based on Python 3.9.16 are 0.41.0 and 1.6.2, respectively. The code using shap and XGBoost to explain the output is available on.2
Statistical analysis
The statistical tools R 3.6.0 and GraphPad Prism 5.01 were used for statistical analysis and figure generation in this study. Data on body measurements, biochemical tests, BMD measurements, and alpha diversity indices are presented as the mean ± standard deviation (SD). Differential analyses between the PMOP and control groups of body measurements, biochemical tests, BMD measurements, and alpha diversity indices were performed with the Wilcoxon test in R, with p < 0.05 indicating a significant difference. For gut microbiota composition, we collected the top 10 dominant microbial taxa at the phylum, class, order, family, genus, and species levels, with p < 0.05 in the Wilcoxon test indicating a significant difference between the PMOP and control groups and 0.05 < p < 0.1 indicating a tendency toward differences. Correlations between gut microbial composition and BMD value/T-score as well as between the abundance of functional group and BMD value/T-score were performed using the Spearman correlation test in R, with p < 0.05 indicating significant correlations and 0.05 < p < 0.1 indicating a tendency for correlations. MedCalc 20.104 and GraphPad Prism 5.01 statistical software were used for logistic regression and receiver operating characteristic (ROC) analysis based on the gut microbiota.
Results
Characteristics of the study population
The baseline demographic features, laboratory biochemical tests and bone BMD measurements of the participants are summarized in Table 1. Based on the diagnostic and exclusion criteria, a total of 58 postmenopausal female participants were included in this study. The average age and menopausal age of 58 participants were 57.47 ± 5.32 and 47.59 ± 2.36, respectively. Generally, there were no significant differences in baseline demographic features or biochemical tests between the PMOP and control groups (p > 0.05 in the Wilcoxon test). However, BMD values and T-scores for both the lumbar spine and total hip of PMOP subjects were significantly lower than those of the control subjects (p < 0.0001). Detailed information regarding demographic features, laboratory biochemical tests and BMD measurements are described in the Methods.
Gut microbial diversity and composition analysis between the PMOP and control groups
To identify the gut microbiota features of the participants and microbial composition changes in patients with PMOP, we sequenced and analyzed the 16S rRNA V3-V4 region of fecal samples from the participants using the Illumina platform. As described in the Methods, we deposited the raw data of the sequencing files into the NCBI-SRA database. The species accumulation curve of 58 samples tended to plateau with increasing sample size, which indicates that the amount of measured data was reasonable and reached the sequencing depth, as shown in Supplementary Figure S1. Generally, no significant differences in the richness and diversity of the gut microbiota were found between PMOP and control subjects, although the control group exhibited slightly higher values of observed species and higher indices of Chao1 and Ace than those of the PMOP group (p > 0.1, Supplementary Table S1). To assess the overall structural changes in the gut microbiota between the PMOP and control groups, PCoA score plots were constructed based on both weighted and unweighted UniFrac distances. As shown in Supplementary Figure S2, the permutational multivariate analysis of variance (PerMANOVA) test based on weighted (p = 0.329) and unweighted UniFrac distance measures (p = 0.025) showed that the microbial structure differed between the PMOP and control subjects.
We then identified the gut microbial compositions in all the samples. The microbial community analysis indicated that Firmicutes, Bacteroidetes, Actinobacteria and Proteobacteria accounted for the majority of the phyla with relative abundances of 71.16, 14.91, 10.02, and 3.59%, respectively, in the total study population. There were no significant differences in the abundances of the four major phyla between the PMOP and control groups (p > 0.1), although the PMOP subjects exhibited relatively higher Fusobacteria abundance than that of the control subjects (p = 0.0395). Figure 1 and Supplementary Table S2 illustrate and summarize the relative microbial compositions in PMOP and control subjects at the phylum, class, order, family, genus and species levels. Notably, the microbial compositions of Bacilli and Erysipelotrichia were higher in the PMOP group than those in the control group (p < 0.1 at the class level), and the relative abundances of Lactobacillales/Lactobacillaceae/Lactobacillus/Lactobacillus salivarius were also increased in PMOP subjects at the order/family/genus/species levels, respectively (p < 0.1). In contrast, the compositions of the bacteria Ruminococcaceae and Bacteroides eggerthii were lower in the PMOP group than those in the control group at the family and species levels (p < 0.05), respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1113174/fmicb-14-1113174-HTML/image_m/fmicb-14-1113174-g001.jpg)
Figure 1. Relative abundance of gut microbiota at different taxonomy levels. The statistical analysis was performed using Wilcoxon test in R, “*” represents p < 0.05, “^” represents 0.05 < p < 0.1. Detailed values about the boxplot were summarized in Supplementary Table S2.
Gut microbial compositions were more correlated with the BMD value/T-score of the total Hip
We subsequently assessed the correlation between gut microbiota composition and BMD value or T-score to further explore BMD-related microbiota in postmenopausal females. It is worth mentioning that there were more correlations between microbial composition and BMD value/T-score observed in the total hip than in the lumbar spine, as demonstrated through Spearman correlation coefficients. In the lumbar spine, only the abundances of two species, Lactobacillus salivarius and Bacteroides eggerthii, were correlated with BMD value/T-score of the participants (p < 0.05). In the total hip, the relative abundances of microbial taxa, including Clostridia/Clostridiales, Lachnospiraceae, Blautia, Streptococcus, Ruminococcus_sp_5_1_39BFAA and Bacteroides eggerthii, were positively correlated with BMD value/T-score; whereas Fusobacteria, Erysipelotrichia/ Erysipelotrichales, Megamonas, Bifidobacterium pseudocatenulatum and Lactobacillus salivarius compositions were negatively correlated with the BMD value/T-score of the participants, detailed results of the Spearman correlation analysis were shown in Figure 2 and Supplementary Table S3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1113174/fmicb-14-1113174-HTML/image_m/fmicb-14-1113174-g002.jpg)
Figure 2. Correlations between microbial composition and BMD value/T-score at both lumbar spine and total hip of the participants. The statistical analysis was performed using Spearman correlation test in R, “*” represents p < 0.05, “^” represents 0.05 < p < 0.1; different colors represent the Spearman correlation coefficients. Detailed values about the correlation analysis were summarized in Supplementary Table S3.
Feature selection methods identified microbial biomarkers for PMOP
To determine whether the gut microbiota can be regarded as an identifying biomarker for distinguishing PMOP from control subjects, two major feature selection methods, the MIC and XGBoost explained with SHAP, were employed, and ROC curves were utilized to test the performance of the classification (details are provided in the Methods). In the feature selection procedure, we analyzed the importance of the 60 major microbial features (as biomarkers in classifying the PMOP and control groups) at all six taxonomic levels and selected the top 20 microbial features based on both the MIC and XGBoost approaches. The maximal information coefficient is a measure of the strength of the linear or nonlinear association between microbes and PMOP, and we identified the top 20 PMOP-related microbes with MIC, as shown in Figure 3A. SHAp values of the top 20 microbial features for every sample are illustrated in Figure 3B. The plot below sorts features by the sum of the SHAP value magnitudes over all samples and uses SHAP values to show the distribution of the impacts that each feature has on the model output. We identified nine microbial features that overlapped between the two feature selection approaches, including the phyla Tenericutes and Fusobacteria; classes Clostridia, Gammaproteobacteria, Negativicutes, and Erysipelotrichia; family Lactobacillaceae, and species Bacteroides caccae and Bacteroides eggerthii. Logistic regression results of the nine microbes indicated that two microbial markers, Fusobacteria and Lactobacillaceae, had significant abilities in disease classification between the PMOP and control groups; the logistic regression equation was Y = −2.920 + 4412*(p_Fusobacteria) + 34.75*(f_Lactobacillaceae) (Table 2). We could distinguish PMOP patients from control subjects based on the combination of p_Fusobacteria and f_Lactobacillaceae, as indicated by the area under the ROC curve (AUC), which had a value of 0.7709 (p < 0.0001, Std. Error: 0.0687, 95% CI: 0.642 to 0.871, Supplementary Figure S3). Overall, the PMOP-associated microbial features captured by the feature selection method offered further evidence of the dysbiotic gut microbiome and highlighted its potential application for the detection of PMOP.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1113174/fmicb-14-1113174-HTML/image_m/fmicb-14-1113174-g003.jpg)
Figure 3. Identified the top 20 microbial features as the classification biomarkers by using feature selection methods. (A) Top 20 microbial features selected by MIC method. (B) Top 20 microbial features selected by XGBoost explained with SHAP. Each point represents a sample, the color represents the feature value (red high, blue low).
Discussion
In the current study, we demonstrated that PMOP patients had different gut microbiota profiles compared with those of control subjects and showed that the compositions of microbiota were more correlated with BMD values of the total hip than those of the lumbar spine in participants. By conducting feature selection, we identified PMOP-related microbes, which could be used as biomarkers for PMOP diagnosis or management.
Our results confirmed that gut microbiota composition were altered in PMOP patients compared with those in normal patients. The abundance of Fusobacteria increased in the PMOP patients (PMOP vs. control, 0.08% vs. 0.03%, p = 0.0395). Consistent with other studies, the relative abundance of Fusobacteria was negatively associated with the BMD value (Greenbaum et al., 2022). Fusobacteria has been shown to promote M1 macrophage production via AKT2 signaling (Liu et al., 2019), which induces inflammation and prompts the development of osteoporosis (Yang and Yang, 2019). In another randomized trial of postmenopausal Japanese women, the total hip BMD of participants increased after receiving the probiotic Bacillus subtilis (C-3102) compared with that of the placebo group (Cronin et al., 2022). Consistent decreases in Fusobacteria species were noted at 12 and 24 weeks after using C-3102, which led the researchers to speculate that the decrease in Fusobacteria might favorably influence BMD by decreasing the production of cytokines that regulate bone resorption (Cronin et al., 2022). In addition, we also observed higher levels of Lactobacillales/Lactobacillaceae/Lactobacillus/Lactobacillus salivarius in PMOP patients than those in the control group, which is contradictory to other results and requires further confirmation. In most related studies, supplementation with probiotics increased BMD values (Yu et al., 2021). For example, supplementation with Lactobacillus reuteri ATCC PTA 6475 reduced bone loss in older women with low BMD (Li et al., 2022). Lactobacillus acidophilus inhibited bone loss and increased bone heterogeneity by modulating Treg-Th17 cell function in osteoporotic mice (Dar et al., 2018). Lactobacillus rhamnosus attenuated bone loss and maintained bone health in ovariectomized mice (Sapra et al., 2021). The above evidence supports the beneficial effects of Lactobacillales/Lactobacillaceae/Lactobacillus on increasing BMD in osteoporosis patients and mouse models. It is worth mentioning that even within the same phylum, each taxon such as family, genus or species has different functions (Ozaki et al., 2021). Generally, Lactobacillales/Lactobacillaceae/Lactobacillus bacteria contribute to weight loss and thus benefit obesity treatment, and studies have shown that the beneficial effects are strain dependent (Crovesy et al., 2017; Alvarez-Arrano and Martin-Pelaez, 2021). Lactobacillus can reduce body weight and alleviate fat accumulation in mice fed a high-fat diet (Wang et al., 2020). For osteoporosis patients, increasing BMI was identified as a protective factor for bone loss (Wang Y et al., 2022; Yin et al., 2022). The results from a large cohort survey showed that the mean rate of fragility fracture was significantly increased in the underweight group compared with the obese and normal weight groups in postmenopausal women (Kim et al., 2021). In an age-adjusted analysis, men who lost weight during follow-up had a significantly greater rate of BMD loss (Sheu et al., 2009). Thus, we may initially obtain a negative correlation between Lactobacillales/Lactobacillaceae/Lactobacillus and body weight or BMI as well as a positive correlation between BMD and body weight or BMI, within a reasonable range. Although this does not represent a negative correlation between Lactobacillales/Lactobacillaceae/Lactobacillus and BMD value, as discussed earlier, most studies have shown that Lactobacillales or Lactobacillus is beneficial for increasing bone mass, and therefore we may speculate that different species of Lactobacillales/Lactobacillaceae/Lactobacillus have different functions, and some strains may be unfavorable for increasing bone mass. Subsequent validation with larger populations and mechanistic experiments should be performed to better explore this hypothesis.
In contrast to the elevated microbial taxonomies, the abundances of Ruminococcaceae (PMOP vs. control, 25.64% vs. 34.99%, p = 0.0332) and Bacteroides eggerthii (PMOP vs. control, nearly 00.00% vs. 0.62%, p = 0.0007) were lower in PMOP patients than those in control subjects. The bacteria of Ruminococcaceae include some butyric acid-producing microbes such as Faecalibacterium and Butyricicoccus (Ozaki et al., 2021). Butyrate-producing bacteria are promising probiotic candidates for gastrointestinal disorders such as inflammatory bowel disease (IBD) (Roux et al., 2015). A recent study demonstrated that Bacteroides eggerthii was decreased in OP patients and negatively correlated with bone resorption markers but positively correlated with bone formation markers and 25-OH-D3 (Qin et al., 2021), which is consistent with our findings in the current study.
We performed BMD measurements at both the lumbar spine and total hip for all the participants in this study. The Spearman correlation analysis revealed that gut microbial compositions were more related to BMD value/T-score at the total hip than those at the lumbar spine. Researchers have suggested that clinicians and densitometrists should expect that more than 40% of women screened for DXA will have T-score discordance between the spine and hip, which is an important point worthy of attention (Yoon and Kim, 2021). In clinical trials, changes in BMD have been shown to provide a reliable estimate of treatment-related antifracture effects, and evidence supports the suitability of total hip BMD as a meaningful outcome for the clinical management of patients with OP (Leslie et al., 2016; Banefelt et al., 2022). In another earlier cohort study, the total hip was the best site for overall fracture assessment, and the spine was the most useful site for the prediction of spinal fractures alone (Leslie et al., 2007). The above evidence and our findings indicate that the BMD value/T-score and BMD-related factors (such as BMD-related microbiota) are discordant among different test sites. This phenomenon should be regarded as a real and prevalent finding for physicians and researchers to develop better strategies for patients (Moayyeri et al., 2005; Leslie et al., 2020).
In the procedure of identifying PMOP-related microbial biomarkers, we highlight the application of feature selection approaches for biological data mining and analysis. Recently, the number of microbiome-related studies has notably increased the availability of data on human microbiome composition and function, which has prompted the application of machine learning in human microbiome studies, including but not limited to feature selection, biomarker identification, disease prediction and treatment (Goren et al., 2021; Marcos-Zambrano et al., 2021; Hinton and Mucha, 2022). Recently, human gut microbiota biomarkers were selected via different feature selection methods for patients with IBD (Bakir-Gungor et al., 2022). In addition to gut microbes (Cai et al., 2015; Ai et al., 2019), advanced feature selection methods also worked well for selecting the key gut microbiome genes (Chen et al., 2016). In this study, we selected several gut microbes that were most strongly related to PMOP via the feature selection methods MIC and XGBoost and further demonstrated that the combination of Fusobacteria and Lactobacillaceae performed well in distinguishing between the PMOP and control groups. We highlight the power and utility of feature selection methods in data mining and processing during biological dataset analysis for identifying diagnostic or predictive markers.
Our study has several limitations. Firstly, the GM analyses associated with this study is limited to its small sample size, which directly influences the research findings. Besides, the research was not profound enough, particularly in relying solely on fecal sample 16S rRNA sequencing to study the intestinal flora. Moreover, multiple factors such as the dietary habit, physical activity, smoking habit influent the GM structure and composition, it’s necessary to consider the confounders of individuals in larger cohort in the future study, to make the research more rigorous.
Conclusion
OP is a chronic, long-term pathologic process that is associated with heredity and environmental factors. PMOP is the most common type of OP and has become a major public health burden for women around the world. Evidence has demonstrated the linkages between OP/PMOP and gut microbiota. Our results showed that gut microbial compositions were altered in the PMOP group compared with the control group. However, the richness and diversity of the gut microbiota between the two groups were not that different. We performed BMD measurements at both the lumbar spine and total hip for all the participants. The Spearman correlation analysis revealed that gut microbial compositions were more strongly related to BMD value/T-score at the total hip than the lumbar spine. Furthermore, we employed two major feature selection methods, MIC and XGBoost, to identify PMOP-related microbial biomarkers, and the combination of specific microbes exhibited potential diagnostic value for distinguishing between PMOP and normal samples. These findings may provide new insights for revealing potential etiologies for OP/PMOP, understanding the role of gut microbiota and modulating gut microbiota as a therapeutic target in the diseases. In addition, we highlight the application of feature selection approaches in biological data mining and data analysis, which may help improve the research techniques for the medical and life sciences.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Ethics statement
The studies involving human participants were reviewed and approved by the Institutional Review Board of Honghui Hospital, Xi’an Jiaotong University (no. 201801019); and was a part of the project “Diversity Analysis for Intestinal Flora in Patients with Primary Osteoporosis” registered at www.chictr.org.cn as #ChiCTR-1,800,019,048#. The patients/participants provided their written informed consent to participate in this study.
Author contributions
DH and JW: project administration and funding acquisition. YZ: cases inclusion and exclusion, blood sample collection, and preparation. YW: data processing and analysis, manuscript writing, and supervision. QL: manuscript revision. All authors contributed to the article and approved the submitted version.
Funding
This work was by supported by National Natural Science Foundation of China (no. 81702067); Shaanxi Provincial Key Research and Development Program (np. 2020GXLH-Y-027; no. 2021SF030); Fundamental Research Funds for the Central Universities (G2020KY0516).
Acknowledgments
The authors would like to thank all the patients who participated in this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1113174/full#supplementary-material
Footnotes
References
Ai, D., Pan, H., Han, R., Li, X., Liu, G., and Xia, L. C. (2019). Using decision tree aggregation with random Forest model to identify gut microbes associated with colorectal cancer. Genes (Basel) 10:112. doi: 10.3390/genes10020112
Alvarez-Arrano, V., and Martin-Pelaez, S. (2021). Effects of probiotics and Synbiotics on weight loss in subjects with overweight or obesity: A systematic review. Nutrients 13:3627. doi: 10.3390/nu13103627
Bakir-Gungor, B., Hacilar, H., Jabeer, A., Nalbantoglu, O. U., Aran, O., and Yousef, M. (2022). Inflammatory bowel disease biomarkers of human gut microbiota selected via different feature selection methods. PeerJ 10:e13205. doi: 10.7717/peerj.13205
Banefelt, J., Timoshanko, J., Soreskog, E., Ortsater, G., Moayyeri, A., Akesson, K. E., et al. (2022). Total hip bone mineral Density as an indicator of fracture risk in bisphosphonate-treated patients in a real-world setting. J. Bone Miner. Res. 37, 52–58. doi: 10.1002/jbmr.4448
Cai, L., Wu, H., Li, D., Zhou, K., and Zou, F. (2015). Type 2 diabetes biomarkers of human gut microbiota selected via iterative sure independent screening method. PLoS One 10:e0140827. doi: 10.1371/journal.pone.0140827
Chen, T., and Guestrin, C. (2016). “Xgboost: A scalable tree boosting system” in Proceedings of the 22nd acm sigkdd International Conference on Knowledge Discovery and Data mining. eds. T. Chen and C. Guestrin, 785–794.
Chen, T., He, T., Benesty, M., Khotilovich, V., Tang, Y., Cho, H., et al. (2015). Xgboost: Extreme Gradient Boosting. R package version 0.4–2 1, 1–4.
Chen, L., Zhang, Y. H., Huang, T., and Cai, Y. D. (2016). Gene expression profiling gut microbiota in different races of humans. Sci. Rep. 6:23075. doi: 10.1038/srep23075
Cronin, O., Lanham-New, S. A., Corfe, B. M., Gregson, C. L., Darling, A. L., Ahmadi, K. R., et al. (2022). Role of the microbiome in regulating bone metabolism and susceptibility to osteoporosis. Calcif. Tissue Int. 110, 273–284. doi: 10.1007/s00223-021-00924-2
Crovesy, L., Ostrowski, M., Ferreira, D., Rosado, E. L., and Soares-Mota, M. (2017). Effect of lactobacillus on body weight and body fat in overweight subjects: a systematic review of randomized controlled clinical trials. Int. J. Obes. 41, 1607–1614. doi: 10.1038/ijo.2017.161
Dar, H. Y., Shukla, P., Mishra, P. K., Anupam, R., Mondal, R. K., Tomar, G. B., et al. (2018). Lactobacillus acidophilus inhibits bone loss and increases bone heterogeneity in osteoporotic mice via modulating Treg-Th17 cell balance. Bone Rep. 8, 46–56. doi: 10.1016/j.bonr.2018.02.001
Ding, K., Hua, F., and Ding, W. (2020). Gut microbiome and osteoporosis. Aging Dis. 11, 438–447. doi: 10.14336/AD.2019.0523
Eastell, R., O'Neill, T. W., Hofbauer, L. C., Langdahl, B., Reid, I. R., Gold, D. T., et al. (2016). Postmenopausal osteoporosis. Nat. Rev. Dis. Primers. 2:16069. doi: 10.1038/nrdp.2016.69
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Goren, E., Wang, C., He, Z., Sheflin, A. M., Chiniquy, D., Prenni, J. E., et al. (2021). Feature selection and causal analysis for microbiome studies in the presence of confounding using standardization. BMC Bioinformatics 22:362. doi: 10.1186/s12859-021-04232-2
Greenbaum, J., Lin, X., Su, K. J., Gong, R., Shen, H., Shen, J., et al. (2022). Integration of the human gut microbiome and serum metabolome reveals novel biological factors involved in the regulation of bone mineral Density. Front. Cell. Infect. Microbiol. 12:853499. doi: 10.3389/fcimb.2022.853499
Hao, M. L., Wang, G. Y., Zuo, X. Q., Qu, C. J., Yao, B. C., and Wang, D. L. (2019). Gut microbiota: an overlooked factor that plays a significant role in osteoporosis. J. Int. Med. Res. 47, 4095–4103. doi: 10.1177/0300060519860027
He, Y., and Chen, Y. (2022). The potential mechanism of the microbiota-gut-bone axis in osteoporosis: a review. Osteoporos. Int. 33, 2495–2506. doi: 10.1007/s00198-022-06557-x
Hinton, A. L., and Mucha, P. J. (2022). A simultaneous feature selection and compositional association test for detecting sparse associations in high-dimensional metagenomic data. Front. Microbiol. 13:837396. doi: 10.3389/fmicb.2022.837396
Hong, S., Cha, K. H., Park, J. H., Jung, D. S., Choi, J. H., Yoo, G., et al. (2022). Cinnamic acid suppresses bone loss via induction of osteoblast differentiation with alteration of gut microbiota. J. Nutr. Biochem. 101:108900. doi: 10.1016/j.jnutbio.2021.108900
Kim, J., Lee, S., Kim, S. S., Lee, J. P., Kim, J. S., Jung, J. G., et al. (2021). Association between body mass index and fragility fracture in postmenopausal women: a cross-sectional study using Korean National Health and nutrition examination survey 2008-2009 (KNHANES IV). BMC Womens Health 21:60. doi: 10.1186/s12905-021-01209-4
Kinney, J. B., and Atwal, G. S. (2014). Equitability, mutual information, and the maximal information coefficient. Proc. Natl. Acad. Sci. 111, 3354–3359. doi: 10.1073/pnas.1309933111
Leslie, W. D., Lix, L. M., Tsang, J. F., and Caetano, P. A., Manitoba Bone Density Program (2007). Single-site vs multisite bone density measurement for fracture prediction. Arch. Intern. Med. 167, 1641–1647. doi: 10.1001/archinte.167.15.1641
Leslie, W. D., Majumdar, S. R., Morin, S. N., and Lix, L. M. (2016). Change in bone mineral Density is an indicator of treatment-related Antifracture effect in routine clinical practice: a registry-based cohort study. Ann. Intern. Med. 165, 465–472. doi: 10.7326/M15-2937
Leslie, W. D., Morin, S. N., Lix, L. M., and Binkley, N. (2020). Impact of spine-hip discordance on fracture risk assessment and treatment qualification in Canada: the Manitoba BMD registry. Arch. Osteoporos. 15:85. doi: 10.1007/s11657-020-00763-y
Li, R., Boer, C. G., Oei, L., and Medina-Gomez, C. (2021). The gut microbiome: a New frontier in musculoskeletal research. Curr. Osteoporos. Rep. 19, 347–357. doi: 10.1007/s11914-021-00675-x
Li, P., Ji, B., Luo, H., Sundh, D., Lorentzon, M., and Nielsen, J. (2022). One-year supplementation with lactobacillus reuteri ATCC PTA 6475 counteracts a degradation of gut microbiota in older women with low bone mineral density. NPJ Biofilms Microbiomes 8:84. doi: 10.1038/s41522-022-00348-2
Li, L., Rao, S., Cheng, Y., Zhuo, X., Deng, C., Xu, N., et al. (2019). Microbial osteoporosis: the interplay between the gut microbiota and bones via host metabolism and immunity. Microbiology 8:e00810. doi: 10.1002/mbo3.810
Liu, L., Liang, L., Liang, H., Wang, M., Lu, B., Xue, M., et al. (2019). Fusobacterium nucleatum aggravates the progression of colitis by regulating M1 macrophage polarization via AKT2 pathway. Front. Immunol. 10:1324. doi: 10.3389/fimmu.2019.01324
Lundberg, S. M., Erion, G., Chen, H., DeGrave, A., Prutkin, J. M., Nair, B., et al. (2020). From local explanations to global understanding with explainable AI for trees. Nat. Mach. Intell. 2, 56–67. doi: 10.1038/s42256-019-0138-9
Lundberg, S. M., and Lee, S.-I. (2017). A unified approach to interpreting model predictions. Adv. Neural Inf. Proces. Syst. 30, 4765–4774. doi: 10.48550/arXiv.1705.07874
Marcos-Zambrano, L. J., Karaduzovic-Hadziabdic, K., Loncar Turukalo, T., Przymus, P., Trajkovik, V., Aasmets, O., et al. (2021). Applications of machine learning in human microbiome studies: a review on feature selection, biomarker identification, disease prediction and treatment. Front. Microbiol. 12:634511. doi: 10.3389/fmicb.2021.634511
Moayyeri, A., Soltani, A., Tabari, N. K., Sadatsafavi, M., Hossein-Neghad, A., and Larijani, B. (2005). Discordance in diagnosis of osteoporosis using spine and hip bone densitometry. BMC Endocr. Disord. 5:3. doi: 10.1186/1472-6823-5-3
Mori, H., Maruyama, F., Kato, H., Toyoda, A., Dozono, A., Ohtsubo, Y., et al. (2014). Design and experimental application of a novel non-degenerate universal primer set that amplifies prokaryotic 16S rRNA genes with a low possibility to amplify eukaryotic rRNA genes. DNA Res. 21, 217–227. doi: 10.1093/dnares/dst052
Ozaki, D., Kubota, R., Maeno, T., Abdelhakim, M., and Hitosugi, N. (2021). Association between gut microbiota, bone metabolism, and fracture risk in postmenopausal Japanese women. Osteoporos. Int. 32, 145–156. doi: 10.1007/s00198-020-05728-y
Parsa, A. B., Movahedi, A., Taghipour, H., Derrible, S., and Mohammadian, A. (2020). Toward safer highways, application of XGBoost and SHAP for real-time accident detection and feature analysis. Accid. Anal. Prev. 136:105405. doi: 10.1016/j.aap.2019.105405
Qin, Q., Yan, S., Yang, Y., Chen, J., Yan, H., Li, T., et al. (2021). The relationship between osteoporosis and intestinal microbes in the Henan Province of China. Front. Cell Dev. Biol. 9:752990. doi: 10.3389/fcell.2021.752990
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Reshef, D. N., Reshef, Y. A., Finucane, H. K., Grossman, S. R., McVean, G., Turnbaugh, P. J., et al. (2011). Detecting novel associations in large data sets. Science 334, 1518–1524. doi: 10.1126/science.1205438
Roux, A. L., El Sayed, F., Duffiet, P., Bauer, T., Heym, B., Gaillard, J. L., et al. (2015). Ruminococcus gnavus total hip arthroplasty infection in a 62-year-old man with ulcerative colitis. J. Clin. Microbiol. 53, 1428–1430. doi: 10.1128/JCM.03040-14
Sapra, L., Dar, H. Y., Bhardwaj, A., Pandey, A., Kumari, S., Azam, Z., et al. (2021). Lactobacillus rhamnosus attenuates bone loss and maintains bone health by skewing Treg-Th17 cell balance in Ovx mice. Sci. Rep. 11:1807. doi: 10.1038/s41598-020-80536-2
Seely, K. D., Kotelko, C. A., Douglas, H., Bealer, B., and Brooks, A. E. (2021). The human gut microbiota: a key mediator of osteoporosis and osteogenesis. Int. J. Mol. Sci. 22:9452. doi: 10.3390/ijms22179452
Sheu, Y., Cauley, J. A., Wheeler, V. W., Patrick, A. L., Bunker, C. H., Kammerer, C. M., et al. (2009). Natural history and correlates of hip BMD loss with aging in men of African ancestry: the Tobago bone health study. J. Bone Miner. Res. 24, 1290–1298. doi: 10.1359/jbmr.090221
Wang, Y., Gao, X., Lv, J., Zeng, Y., Li, Q., Wang, L., et al. (2022). Gut microbiome signature are correlated with bone mineral Density alterations in the Chinese elders. Front. Cell. Infect. Microbiol. 12:827575. doi: 10.3389/fcimb.2022.827575
Wang, J., Wang, Y., Gao, W., Wang, B., Zhao, H., Zeng, Y., et al. (2017). Diversity analysis of gut microbiota in osteoporosis and osteopenia patients. PeerJ 5:e3450. doi: 10.7717/peerj.3450
Wang, J., Wu, S., Zhang, Y., Yang, J., and Hu, Z. (2022). Gut microbiota and calcium balance. Front. Microbiol. 13:1033933. doi: 10.3389/fmicb.2022.1033933
Wang, M., Zhang, B., Hu, J., Nie, S., Xiong, T., and Xie, M. (2020). Intervention of five strains of lactobacillus on obesity in mice induced by high-fat diet. J. Funct. Foods 72:104078. doi: 10.1016/j.jff.2020.104078
Xu, Z., Xie, Z., Sun, J., Huang, S., Chen, Y., Li, C., et al. (2020). Gut microbiome reveals specific Dysbiosis in primary osteoporosis. Front. Cell. Infect. Microbiol. 10:160. doi: 10.3389/fcimb.2020.00160
Yan, Q., Cai, L., and Guo, W. (2022). New advances in improving bone health based on specific gut microbiota. Front. Cell. Infect. Microbiol. 12:821429. doi: 10.3389/fcimb.2022.821429
Yan, J., Herzog, J. W., Tsang, K., Brennan, C. A., Bower, M. A., Garrett, W. S., et al. (2016). Gut microbiota induce IGF-1 and promote bone formation and growth. Proc. Natl. Acad. Sci. U. S. A. 113, E7554–E7563. doi: 10.1073/pnas.1607235113
Yang, D. H., and Yang, M. Y. (2019). The role of macrophage in the pathogenesis of osteoporosis. Int. J. Mol. Sci. 20:2093. doi: 10.3390/ijms20092093
Yin, Y., Lu, X., Li, Z., Liu, S., Shao, L., Cao, L., et al. (2022). Risk factors for worsening of bone loss in patients newly diagnosed with inflammatory bowel disease. Gastroenterol. Res. Pract. 2022, 1498293–1498210. doi: 10.1155/2022/1498293
Yoon, B. H., and Kim, D. Y. (2021). Discordance between hip and spine bone mineral Density: a point of care. J. Bone Metab. 28, 249–251. doi: 10.11005/jbm.2021.28.4.249
Keywords: postmenopausal osteoporosis, gut microbiota, bone mineral density, feature selection, microbial biomarker
Citation: Huang D, Wang J, Zeng Y, Li Q and Wang Y (2023) Identifying microbial signatures for patients with postmenopausal osteoporosis using gut microbiota analyses and feature selection approaches. Front. Microbiol. 14:1113174. doi: 10.3389/fmicb.2023.1113174
Edited by:
Katri Korpela, University of Helsinki, FinlandCopyright © 2023 Huang, Wang, Zeng, Li and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qingmei Li, liqingmei09@163.com; Yangyang Wang, wangyang2154@mail.nwpu.edu.cn
†These authors have contributed equally to this work