Metabolomics Insights into Osteoporosis Through Association With Bone Mineral Density

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).

[1]  Kai Yu,et al.  Power calculation for the general two‐sample Mendelian randomization analysis , 2020, Genetic epidemiology.

[2]  Hui Shen,et al.  A road map for understanding molecular and genetic determinants of osteoporosis , 2019, Nature Reviews Endocrinology.

[3]  C. Cheung,et al.  Serum Metabolome of Coffee Consumption and its Association with Bone Mineral Density: The Hong Kong Osteoporosis Study. , 2019, The Journal of clinical endocrinology and metabolism.

[4]  Ji-Gang Zhang,et al.  Metabolomic profiles associated with bone mineral density in US Caucasian women , 2018, Nutrition & Metabolism.

[5]  M. Tomita,et al.  Metabolomics-based profiles predictive of low bone mass in menopausal women , 2018, Bone reports.

[6]  David S. Wishart,et al.  MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis , 2018, Nucleic Acids Res..

[7]  J. Totman,et al.  Association of Plasma Lipids and Polar Metabolites with Low Bone Mineral Density in Singaporean-Chinese Menopausal Women: A Pilot Study , 2018, International journal of environmental research and public health.

[8]  B. Neale,et al.  Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases , 2018, Nature Genetics.

[9]  C. Cheung,et al.  Cohort Profile: The Hong Kong Osteoporosis Study and the follow-up study. , 2018, International journal of epidemiology.

[10]  T. Spector,et al.  Metabolomic Pathways to Osteoporosis in Middle‐Aged Women: A Genome‐Metabolome‐Wide Mendelian Randomization Study , 2018, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[11]  Qi Zhang,et al.  Association between the metabolome and bone mineral density in pre- and post-menopausal Chinese women using GC-MS. , 2016, Molecular bioSystems.

[12]  B. Gualano,et al.  Creatine supplementation in the aging population: effects on skeletal muscle, bone and brain , 2016, Amino Acids.

[13]  G. Davey Smith,et al.  Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator , 2016, Genetic epidemiology.

[14]  M. Pirinen,et al.  Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA , 2016, Nature Communications.

[15]  G. Davey Smith,et al.  Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression , 2015, International journal of epidemiology.

[16]  S. Vollset,et al.  Plasma dimethylglycine, nicotine exposure and risk of low bone mineral density and hip fracture: the Hordaland Health Study , 2015, Osteoporosis International.

[17]  Jeng-Min Chiou,et al.  Association Between the Metabolome and Low Bone Mineral Density in Taiwanese Women Determined by 1H NMR Spectroscopy , 2014, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[18]  A. Butterworth,et al.  Mendelian Randomization Analysis With Multiple Genetic Variants Using Summarized Data , 2013, Genetic epidemiology.

[19]  D. Kiel,et al.  Plasma phosphatidylcholine concentrations of polyunsaturated fatty acids are differentially associated with hip bone mineral density and hip fracture in older adults: The framingham osteoporosis study , 2012, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[20]  Daniel L. Koller,et al.  Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture , 2012, Nature Genetics.

[21]  Jianguo Xia,et al.  Web-based inference of biological patterns, functions and pathways from metabolomic data using MetaboAnalyst , 2011, Nature Protocols.

[22]  S. Carr,et al.  Lipid profiling identifies a triacylglycerol signature of insulin resistance and improves diabetes prediction in humans. , 2011, The Journal of clinical investigation.

[23]  V. Mootha,et al.  Metabolite profiles and the risk of developing diabetes , 2011, Nature Medicine.

[24]  M. Thorén,et al.  Associations between amino acids and bone mineral density in men with idiopathic osteoporosis. , 2010, Bone.

[25]  B. Dawson-Hughes A revised clinician's guide to the prevention and treatment of osteoporosis. , 2008, The Journal of clinical endocrinology and metabolism.

[26]  M. Pencina,et al.  Evaluating the added predictive ability of a new marker: From area under the ROC curve to reclassification and beyond , 2008, Statistics in medicine.

[27]  O. Johnell,et al.  FRAX™ and the assessment of fracture probability in men and women from the UK , 2008, Osteoporosis International.

[28]  M. Lever,et al.  Homocysteine, glycine betaine, and N,N-dimethylglycine in patients attending a lipid clinic. , 2005, Metabolism: clinical and experimental.

[29]  S Gary Firestein,et al.  Kelley's Textbook of Rheumatology , 2004 .

[30]  J. Kanis,et al.  Diagnosis of osteoporosis and assessment of fracture risk , 2002, The Lancet.

[31]  Susan R. Johnson,et al.  Osteoporosis prevention, diagnosis, and therapy. , 2001, JAMA.

[32]  D. Kiel,et al.  Risk Factors for Longitudinal Bone Loss in Elderly Men and Women: The Framingham Osteoporosis Study , 2000, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[33]  R. Tibshirani Regression Shrinkage and Selection via the Lasso , 1996 .

[34]  O. Johnell,et al.  Risk factors for hip fracture in european women: The MEDOS study , 1995 .

[35]  M. Bonneu,et al.  Assessment of the risk of post-menopausal osteoporosis using clinical factors. , 1992, Clinical endocrinology.

[36]  M. Bonneu,et al.  Assessment of the risk of post‐menopausal osteoporosis using clinical factors , 1992 .

[37]  E. DeLong,et al.  Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. , 1988, Biometrics.

[38]  W. Kannel,et al.  An investigation of coronary heart disease in families. The Framingham offspring study. , 1979, American journal of epidemiology.

[39]  T. Dawber,et al.  Epidemiological approaches to heart disease: the Framingham Study. , 1951, American journal of public health and the nation's health.