Metabolome-Wide Mendelian Randomization Assessing the Causal Relationship Between Blood Metabolites and Bone Mineral Density

Mounting evidence has supported osteoporosis (OP) as a metabolic disorder. Recent metabolomics studies have discovered numerous metabolites related to bone mineral density (BMD). However, the causal effects of metabolites on BMD at distinct sites remained underexplored. Leveraging genome-wide association datasets, we conducted two-sample Mendelian randomization (MR) analyses to investigate the causal relationship between 486 blood metabolites and bone mineral density at five skeletal sites including heel (H), total body (TB), lumbar spine (LS), femoral neck (FN), and ultra-distal forearm (FA). Sensitivity analyses were performed to test the presence of the heterogeneity and the pleiotropy. To exclude the influences of reverse causation, genetic correlation, and linkage disequilibrium (LD), we further performed reverse MR, linkage disequilibrium regression score (LDSC), and colocalization analyses. In the primary MR analyses, 22, 10, 3, 7, and 2 metabolite associations were established respectively for H-BMD, TB-BMD, LS-BMD, FN-BMD, and FA-BMD at the nominal significance level (IVW, P < 0.05) and passing sensitivity analyses. Among these, one metabolite, androsterone sulfate showed a strong effect on four out of five BMD phenotypes (Odds ratio [OR] for H-BMD = 1.045 [1.020, 1.071]; Odds ratio [OR] for TB-BMD = 1.061 [1.017, 1.107]; Odds ratio [OR] for LS-BMD = 1.088 [1.023, 1.159]; Odds ratio [OR] for FN-BMD = 1.114 [1.054, 1.177]). Reverse MR analysis provided no evidence for the causal effects of BMD measurements on these metabolites. Colocalization analysis have found that several metabolite associations might be driven by shared genetic variants such as mannose for TB-BMD. This study identified some metabolites causally related to BMD at distinct sites and several key metabolic pathways, which shed light on predictive biomarkers and drug targets for OP.