Dataset Information

Plasma high-resolution metabolomics identifies linoleic acid and linked metabolic pathways associated with bone mineral density.

ABSTRACT:

Background & aims

There is a considerable degree of variation in bone mineral density (BMD) within populations. Use of plasma metabolomics may provide insight into established and novel determinants of BMD variance, such as nutrition and gut microbiome composition, to inform future prevention and treatment strategies for loss of BMD. Using high-resolution metabolomics (HRM), we examined low-molecular weight plasma metabolites and nutrition-related metabolic pathways associated with BMD.

Methods

This cross-sectional study included 179 adults (mean age 49.5 ± 10.3 yr, 64% female). Fasting plasma was analyzed using ultra-high-resolution mass spectrometry with liquid chromatography. Whole body and spine BMD were assessed by dual energy X-ray absorptiometry and expressed as BMD (g/cm²) or Z-scores. Multiple linear regression, pathway enrichment, and module analyses were used to determine key plasma metabolic features associated with bone density.

Results

Of 10,210 total detected metabolic features, whole body BMD Z-score was associated with 710 metabolites, which were significantly enriched in seven metabolic pathways, including linoleic acid, fatty acid activation and biosynthesis, and glycerophospholipid metabolism. Spine BMD was associated with 970 metabolites, significantly enriched in pro-inflammatory pathways involved in prostaglandin formation and linoleic acid metabolism. In module analyses, tryptophan- and polyamine-derived metabolites formed a network that was significantly associated with spine BMD, supporting a link with the gut microbiome.

Conclusions

Plasma HRM provides comprehensive information relevant to nutrition and components of the microbiome that influence bone health. This data supports pro-inflammatory fatty acids and the gut microbiome as novel regulators of postnatal bone remodeling.

SUBMITTER: Bellissimo MP

PROVIDER: S-EPMC7714706 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Osteoporosis, a disease characterized by low bone mineral density (BMD), increases the risk for fractures. Conventional risk factors alone do not completely explain measured BMD or osteoporotic fracture risk. Metabolomics may provide additional information. We aim to identify BMD-associated metabolomic markers that are predictive of fracture risk. We assessed 209 plasma metabolites by liquid chromatography with tandem mass spectrometry (LC-MS/MS) in 1552 Framingham Offspring Study participants, and measured femoral neck (FN) and lumbar spine (LS) BMD 2 to 10 years later using dual-energy X-ray absorptiometry. We assessed osteoporotic fractures up to 27-year follow-up after metabolomic profiling. We identified 27 metabolites associated with FN-BMD or LS-BMD by LASSO regression with internal validation. Incorporating selected metabolites significantly improved the prediction and the classification of osteoporotic fracture risk beyond conventional risk factors (area under the curve [AUC] = 0.74 for the model with identified metabolites and risk factors versus AUC = 0.70 with risk factors alone, p = .001; net reclassification index = 0.07, p = .03). We replicated significant improvement in fracture prediction by incorporating selected metabolites in 634 participants from the Hong Kong Osteoporosis Study (HKOS). The glycine, serine, and threonine metabolism pathway (including four identified metabolites: creatine, dimethylglycine, glycine, and serine) was significantly enriched (false discovery rate [FDR] p value = .028). Furthermore, three causally related metabolites (glycine, phosphatidylcholine [PC], and triacylglycerol [TAG]) were negatively associated with FN-BMD, whereas PC and TAG were negatively associated with LS-BMD through Mendelian randomization analysis. In summary, metabolites associated with BMD are helpful in osteoporotic fracture risk prediction. Potential causal mechanisms explaining the three metabolites on BMD are worthy of further experimental validation. Our findings may provide novel insights into the pathogenesis of osteoporosis. © 2021 American Society for Bone and Mineral Research (ASBMR).

Project description:We previously demonstrated that skeletal mass, structure, and biomechanical properties vary considerably in heterogeneous stock (HS) rat strains. In addition, we observed strong heritability for several of these skeletal phenotypes in the HS rat model, suggesting that it represents a unique genetic resource for dissecting the complex genetics underlying bone fragility. The purpose of this study was to identify and localize genes associated with bone mineral density in HS rats. We measured bone phenotypes from 1524 adult male and female HS rats between 17 and 20 weeks of age. Phenotypes included dual-energy X-ray absorptiometry (DXA) measurements for bone mineral content and areal bone mineral density (aBMD) for femur and lumbar spine (L3-L5), and volumetric BMD measurements by CT for the midshaft and distal femur, femur neck, and fifth lumbar vertebra (L5). A total of 70,000 polymorphic single-nucleotide polymorphisms (SNPs) distributed throughout the genome were selected from genotypes obtained from the Affymetrix rat custom SNPs array for the HS rat population. These SNPs spanned the HS rat genome with a mean linkage disequilibrium coefficient between neighboring SNPs of 0.95. Haplotypes were estimated across the entire genome for each rat using a multipoint haplotype reconstruction method, which calculates the probability of descent for each genotyped locus from each of the eight founder HS strains. The haplotypes were tested for association with each bone density phenotype via a mixed model with covariate adjustment. We identified quantitative trait loci (QTLs) for BMD phenotypes on chromosomes 2, 9, 10, and 13 meeting a conservative genomewide empiric significance threshold (false discovery rate [FDR]?= 5%; p < 3 × 10(-6)). Importantly, most QTLs were localized to very small genomic regions (1-3 megabases [Mb]), allowing us to identify a narrow set of potential candidate genes including both novel genes and genes previously shown to have roles in skeletal development and homeostasis.