Although the gut microbiota has been reported to influence osteoporosis risk, the individual species involved, and underlying mechanisms, remain largely unknown. We performed integrative analyses in a Chinese cohort of peri-/post-menopausal women with metagenomics/targeted metabolomics/whole-genome sequencing to identify novel microbiome-related biomarkers for bone health. Bacteroides vulgatus was found to be negatively associated with bone mineral density (BMD), which was validated in US white people. Serum valeric acid (VA), a microbiota derived metabolite, was positively associated with BMD and causally downregulated by B. vulgatus. Ovariectomized mice fed B. vulgatus demonstrated increased bone resorption and poorer bone micro-structure, while those fed VA demonstrated reduced bone resorption and better bone micro-structure. VA suppressed RELA protein production (pro-inflammatory), and enhanced IL10 mRNA expression (anti-inflammatory), leading to suppressed maturation of osteoclast-like cells and enhanced maturation of osteoblasts in vitro. The findings suggest that B. vulgatus and VA may represent promising targets for osteoporosis prevention/treatment.

Osteoporosis is a highly prevalent chronic aging-related disease that frequently is only detected after fracture. We hypothesized that aminobutyric acids could serve as biomarkers for osteoporosis. We developed a quick, accurate, and sensitive screening method for aminobutyric acid isomers and enantiomers yielding correlations with bone mineral density (BMD) and osteoporotic fracture. In serum, γ-aminobutyric acid (GABA) and (R)-3-aminoisobutyric acid (D-BAIBA) have positive associations with physical activity in young lean women.D-BAIBA positively associated with hip BMD in older individuals without osteoporosis/osteopenia. Lower levels of GABA were observed in 60–80 year old women with osteoporotic fractures. Single nucleotide polymorphisms in seven genes related to these metabolites associated with BMD and osteoporosis. In peripheral blood monocytes, dihydropyrimidine dehydrogenase, an enzyme essential toD-BAIBA generation, exhibited positive association with physical activity and hip BMD. Along with their signaling roles, BAIBA and GABA might serve as biomarkers for diagnosis and treatments of osteoporosis.

Individuals’ peak bone mineral density (BMD) achieved and maintained at ages 20–40 years is the most powerful predictor of low bone mass and osteoporotic fractures later in life. The aim of this study was to identify metabolomic factors associated with peak BMD variation in US Caucasian women. A total of 136 women aged 20–40 years, including 65 subjects with low and 71 with high hip BMD, were enrolled. The serum metabolites were assessed using a liquid chromatography-mass spectrometry (LC-MS) method. The partial least-squares discriminant analysis (PLS-DA) method and logistic regression models were used, respectively, to examine the associations of metabolomic profiles and individual metabolites with BMD. The low and high BMD groups could be differentiated by the detected serum metabolites using PLS-DA (P permutation = 0.008). A total of 14 metabolites, including seven amino acids and amino acid derivatives, five lipids (including three bile acids), and two organic acids, were significantly associated with the risk for low BMD. Most of these metabolites are novel in that they have never been linked with BMD in humans earlier. The prediction model including the newly identified metabolites significantly improved the classification of the groups with low and high BMD. The area under the receiver operating characteristic curve without and with metabolites were 0.88 (95% CI: 0.83–0.94) and 0.97 (95% CI: 0.94–0.99), respectively (P for the difference = 0.0004). Metabolomic profiling may improve the risk prediction of osteoporosis among Caucasian women. Our findings also suggest the potential importance of the metabolism of amino acids and bile acids in bone health.

At present, there have only been a few DNA sequencing-based studies to explore the genetic determinants of bone mineral density (BMD). We carried out the largest whole genome sequencing analysis to date for femoral neck and spine BMD (n = 4981), with one of the highest average sequencing depths implemented thus far at 22×, in a multiethnic sample (58% Caucasian and 42% African American) from the Louisiana Osteoporosis Study (LOS). The LOS samples were combined with summary statistics from the GEFOS consortium and several independent samples of various ethnicities to perform GWAS meta-analysis (n = 44 506). We identified 31 and 30 genomic risk loci for femoral neck and spine BMD, respectively. The findings substantiate many previously reported susceptibility loci (e.g. WNT16 and ESR1) and reveal several others that are either novel or have not been widely replicated in GWAS for BMD, including two for femoral neck (IGF2 and ZNF423) and one for spine (SIPA1). Although we were not able to uncover ethnicity specific differences in the genetic determinants of BMD, we did identify several loci which demonstrated sex-specific associations, including two for women (PDE4D and PIGN) and three for men (TRAF3IP2, NFIB and LYSMD4). Gene-based rare variant association testing detected MAML2, a regulator of the Notch signaling pathway, which has not previously been suggested, for association with spine BMD. The findings provide novel insights into the pathophysiological mechanisms of osteoporosis.

Background: Osteoporosis is a major global health issue, weakening bones and increasing fracture risk. Dual-energy X-ray absorptiometry (DXA) is the standard for measuring bone mineral density (BMD) and diagnosing osteoporosis, but its costliness and complexity impede widespread screening adoption. Predictive modeling using genetic and clinical data offers a cost-effective alternative for assessing osteoporosis and fracture risk. This study aims to develop BMD prediction models using data from the UK Biobank (UKBB) and test their performance across different ethnic and geographical populations. Methods and findings: We developed BMD prediction models for the femoral neck (FNK) and lumbar spine (SPN) using both genetic variants and clinical factors (such as sex, age, height, and weight), within 17,964 British white individuals from UKBB. Models based on regression with least absolute shrinkage and selection operator (LASSO), selected based on the coefficient of determination (R²) from a model selection subset of 5,973 individuals from British white population. These models were tested on 5 UKBB test sets and 12 independent cohorts of diverse ancestries, totaling over 15,000 individuals. Furthermore, we assessed the correlation of predicted BMDs with fragility fractures risk in 10 years in a case-control set of 287,183 European white participants without DXA-BMDs in the UKBB. With single-nucleotide polymorphism (SNP) inclusion thresholds at 5×10⁻⁶ and 5×10⁻⁷, the prediction models for FNK-BMD and SPN-BMD achieved the highest R² of 27.70% with a 95% confidence interval (CI) of [27.56%, 27.84%] and 48.28% (95% CI [48.23%, 48.34%]), respectively. Adding genetic factors improved predictions slightly, explaining an additional 2.3% variation for FNK-BMD and 3% for SPN-BMD over clinical factors alone. Survival analysis revealed that the predicted FNK-BMD and SPN-BMD were significantly associated with fragility fracture risk in the European white population (P < 0.001). The hazard ratios (HRs) of the predicted FNK-BMD and SPN-BMD were 0.83 (95% CI [0.79, 0.88], corresponding to a 1.44% difference in 10-year absolute risk) and 0.72 (95% CI [0.68, 0.76], corresponding to a 1.64% difference in 10-year absolute risk), respectively, indicating that for every increase of one standard deviation in BMD, the fracture risk will decrease by 17% and 28%, respectively. However, the model’s performance declined in other ethnic groups and independent cohorts. The limitations of this study include differences in clinical factors distribution and the use of only SNPs as genetic factors. Conclusions: In this study, we observed that combining genetic and clinical factors improves BMD prediction compared to clinical factors alone. Adjusting inclusion thresholds for genetic variants (e.g., 5×10⁻⁶ or 5×10⁻⁷) rather than solely considering genome-wide association study (GWAS)-significant variants can enhance the model’s explanatory power. The study highlights the need for training models on diverse populations to improve predictive performance across various ethnic and geographical groups.