Life-Course Genome-wide Association Study Meta-analysis of Total Body BMD and Assessment of Age-Specific Effects

Assessing the causal relationship between plasma proteins and osteoporosis: novel insights into pathological mechanisms and therapeutic implications

Identifying dysregulated plasma proteins in osteoporosis (OP) progression offers insights into prevention and treatment. This study found 8 such proteins associated with OP, suggesting them as therapy targets. This discovery may cut drug development costs and improve personalized treatments.

PURPOSE

This study aims to identify potential therapeutic targets for OP using summary data-based Mendelian randomization (SMR) and colocalization analysis methods. Furthermore, we seek to explore the biological significance and pharmacological value of these drug targets.

METHODS

To identify potential therapeutic targets for OP, we conducted SMR and colocalization analysis. Plasma protein (pQTL, exposure) data were sourced from the study by Ferkingstad et al. (n = 35,559). Summary statistics for bone mineral density (BMD, outcome) were obtained from the GWAS Catalog (n = 56,284). Additionally, we utilized enrichment analysis, protein-protein interaction (PPI) network analysis, drug prediction, and molecular docking to further analyze the biological significance and pharmacological value of these drug targets.

RESULTS

In the SMR analysis, while 20 proteins showed significance, only 8 potential drug targets (GCKR, ERBB3, CFHR1, GPN1, SDF2, VTN, BET1L, and SERPING1) received support from colocalization (PP.H4 > 0.8). These proteins are closely associated with immune function in terms of biological significance. Molecular docking also demonstrated favorable binding of drugs to proteins, consistent with existing structural data, further substantiating the pharmacological value of these targets.

CONCLUSIONS

The study identified 8 potential drug targets for OP. These prospective targets are believed to have a higher chance of success in clinical trials, thus aiding in prioritizing OP drug development and reducing development costs.

Causal relationship between circulating levels of cytokines and bone mineral density: A mendelian randomization study

BACKGROUND

Numerous studies have shown that various cytokines are important factors affecting bone mineral density (BMD), but the causality between the two remains uncertain.

METHODS

Genetic variants associated with 41 circulating cytokines from a genome-wide association study (GWAS) in 8,293 Finns were used as instrumental variables (IVs) for a two-sample Mendelian randomization (MR) analysis. Inverse variance weighting (IVW) was employed as the primary method to investigate whether the 41 cytokines were causally associated with BMD at five different sites [total body bone mineral density (TB-BMD), heel bone mineral density (HE-BMD), forearm bone mineral density (FA-BMD), femoral neck bone mineral density (FN-BMD), and lumbar spine bone mineral density (LS-BMD)]. Weighted median and MR-Egger were chosen to further confirm the robustness of the results. We performed MR pleiotropy residual sum and outlier test (MR-PRESSO), MR-Egger regression, and Cochran's Q test to detect pleiotropy and sensitivity testing.

RESULTS

After Bonferroni correction, two circulating cytokines had a strong causality with BMD at corresponding sites. Genetically predicted circulating hepatocyte growth factor (HGF) levels and HE-BMD were negatively correlated [β (95 % CI) -0.035(-0.055, -0.016), P=0.00038]. Circulating macrophage inflammatory protein-1α (MIP-1α) levels and TB-BMD were negatively correlated [β(95 %CI): -0.058(-0.092, -0.024), P=0.00074]. Weighted median and MR-Egger results were in line with the IVW results. We also found suggestive causal relationship (IVW P<0.05) between seven circulating cytokines and BMD at corresponding sites. No significant pleiotropy or heterogeneity was observed in our study.

CONCLUSION

Our MR analyses indicated a causal effect between two circulating cytokines and BMD at corresponding sites (HGF and HE-BMD, MIP-1α and TB-BMD), along with suggestive evidence of a potential causality between seven cytokines and BMD at the corresponding sites. These findings would provide insights into the prevention and treatment of osteoporosis, especially immunoporosis.

Enhanced fracture risk prediction: a novel multi-trait genetic approach integrating polygenic scores of fracture-related traits

The novel metaPGS, integrating multiple fracture-related genetic traits, surpasses traditional polygenic scores in predicting fracture risk. Demonstrating a robust association with incident fractures, this metaPGS offers significant potential for enhancing clinical fracture risk assessment and tailoring prevention strategies.

INTRODUCTION

Current polygenic scores (PGS) have limited predictive power for fracture risk. To improve genetic prediction, we developed and evaluated a novel metaPGS combining genetic information from multiple fracture-related traits.

METHODS

We derived individual PGS from genome-wide association studies of 16 fracture-related traits and employed an elastic-net logistic regression model to examine the association between the 16 PGSs and fractures. An optimal metaPGS was constructed by combining 11 significant individual PGSs selected by the elastic regularized regression model. We evaluated the predictive power of the metaPGS alone and in combination with clinical risk factors recommended by guidelines. The discrimination ability of metaPGS was assessed using the concordance index. Reclassification was assessed using net reclassification improvement (NRI) and integrated discrimination improvement (IDI).

RESULTS

The metaPGS had a significant association with incident fractures (HR 1.21, 95% CI 1.18-1.25 per standard deviation of metaPGS), which was stronger than previously developed bone mineral density (BMD)-related individual PGSs. Models with PGS_FNBMD, PGS_TBBMD, and metaPGS had slightly higher but statistically non-significant c-index than the base model (0.640, 0.644, 0.644 vs. 0.638). However, the reclassification analysis showed that compared to the base model, the model with metaPGS improves the reclassification of fracture.

CONCLUSIONS

The metaPGS is a promising approach for stratifying fracture risk in the European population, improving fracture risk prediction by combining genetic information from multiple fracture-related traits.

Obesity Variants in the GIPR Gene Are not Associated With Risk of Fracture or Bone Mineral Density

CONTEXT

It is not clear if antagonizing the GIP (glucose-dependent insulinotropic polypeptide) receptor (GIPR) for treatment of obesity is likely to increase the risk of fractures, or to lower bone mineral density (BMD) beyond what is expected with rapid weight loss.

OBJECTIVE

The objective of this study was to investigate the risk of fracture and BMD of sequence variants in GIPR that reduce the activity of the GIP receptor and have been associated with reduced body mass index (BMI).

METHODS

We analyzed the association of 3 missense variants in GIPR, a common variant, rs1800437 (p.Glu354Gln), and 2 rare variants, rs139215588 (p.Arg190Gln) and rs143430880 (p.Glu288Gly), as well as a burden of predicted loss-of-function (LoF) variants with risk of fracture and with BMD in a large meta-analysis of up to 1.2 million participants. We analyzed associations with fractures at different skeletal sites in the general population: any fractures, hip fractures, vertebral fractures and forearm fractures, and specifically nonvertebral and osteoporotic fractures in postmenopausal women. We also evaluated associations with BMD at the lumbar spine, femoral neck, and total body measured with dual-energy x-ray absorptiometry (DXA), and with BMD estimated from heel ultrasound (eBMD).

RESULTS

None of the 3 missense variants in GIPR was significantly associated with increased risk of fractures or with lower BMD. Burden of LoF variants in GIPR was not associated with fractures or with BMD measured with clinically validated DXA, but was associated with eBMD.

CONCLUSION

Missense variants in GIPR, or burden of LoF variants in the gene, are not associated with risk of fractures or with lower BMD.

The Causal Relationship between PCSK9 Inhibitors and Osteoporosis Based on Drug-Targeted Mendelian Combined Mediation Analysis

PCSK9 inhibitors have been shown to lower serum low density lipoprotein cholesterol (LDL-C) levels and are considered integral in the treatment of cardiovascular diseases. However, the potential association between PCSK9 inhibitors and osteoporosis is unclear now. In this study, drug-targeted mendelian randomization (MR) was utilized in conjunction with mediation analysis including bone mineral density (BMD), total 25-hydroxyvitamin D (T25(OH)D) levels and calcium supplementation to investigate the causal relationship between PCSK9 inhibitors and osteoporosis. The LDL-C level was chosen as the exposure variable in a sample size of 173,082 individuals. We conducted a MR analysis on the relationship between PCSK9 inhibitors and osteoporosis, elucidating the mediators involved. Utilizing the inverse variance weighted (IVW) method, we found the risk of osteoporosis was reduced by 0.6% in those who used PCSK9 inhibitors compared with non-users (OR: 0.994, 95%CI: 0.991-0.998, P < 0.001). In people aged 30-45 years, the risk of low BMD was 1.176 times higher among PCSK9 inhibitor users compared to non-users (OR: 1.176, 95%CI: 1.017-1.336, P = 0.045). Conversely, people aged 45-60 years who used PCSK9 inhibitors had a 14.9% lower risk of low BMD compared to non-users (OR: 0.851, 95%CI: 0.732-0.968, P = 0.007). Mediation analysis revealed that 43.33% of the impact of PCSK9 inhibitors on osteoporosis was mediated through BMD levels, with the remaining 56.67% being a direct effect. Effects of PCSK9 inhibitors on BMD levels varied in different ages. In addition, the risk of high serum T25(OH)D levels were 1.091 times among PCSK9 inhibitor users compared to non-users (OR: 1.091, 95%CI: 1.065-1.112, P < 0.001), providing valuable insights for clinicians.

Genetic Prediction of Osteoporosis by Anti-Müllerian Hormone Levels and Reproductive Factors in Women: A Mendelian Randomization Study

Previous observational studies have suggested that anti-Müllerian hormone (AMH) and reproductive factors are linked to reduced bone mineral density (BMD) and an increased risk of osteoporosis (OP) in women. However, related studies are limited, and these traditional observational studies may be subject to residual confounders and reverse causation, while also lacking a more comprehensive observation of various reproductive factors. Univariate and multivariate two-sample Mendelian randomization analyses were conducted to determine the causal associations of AMH levels and six reproductive factors with BMD and OP, using the random-effects inverse-variance weighted method. Heterogeneity was assessed using Cochran's Q-statistic, and sensitivity analyses were performed to identify causal correlations. Age at menarche (AAM) was negatively associated with total body BMD (TB-BMD) in females aged 45-60 and over 60 years, as well as with heel bone mineral density (eBMD). Conversely, age at natural menopause (ANM) was positively associated with TB-BMD in the same age ranges and with eBMD. ANM was only causally associated with self-reported OP and showed no significant correlation with definitively diagnosed OP. Neither AMH level nor other reproductive factors were significantly associated with a genetic predisposition to BMD at any age and OP. Later AAM and earlier ANM are significantly genetically causally associated with decreased BMD but not with OP. AMH levels, length of menstrual cycle, age at first birth, age at last birth, and number of live births, in terms of genetic backgrounds, are not causally related to BMD or OP.

Mendelian randomization study on association between grip strength and BMD in different age groups

INTRODUCTION

This study aimed to use the Mendelian randomization study method to verify the causal relationship between grip strength and bone mineral density (BMD) in different ages and different parts of the body.

MATERIALS AND METHODS

The analysis was based on pooled data from genome-wide association studies (GWAS). Hand grip strength (right) was used as the exposure variable and total body bone mineral density (BMD) of different age groups was used as the outcome variable. Single-nucleotide polymorphisms highly correlated with exposure variables were used as instrumental variables. The inverse variance weighted (IVW) method was used as the primary analysis method, and the Mendelian randomization Egger (MR-Egger) regression and weighted median methods were used as supplementary evidence for the IVW results. Horizontal pleiotropy and heterogeneity tests were conducted to ensure the stability of the results.

RESULTS

Analyzing the GWAS data on osteoporosis as the outcome variable, the IVW analysis showed that osteoporosis risk was associated with decreased grip strength in the 45-60 age group and the risk of declining lumbar spine BMD was associated with decreased grip strength. However, there was no significant correlation between the risk of osteoporosis in other age groups and changes in grip strength.

CONCLUSION

A causal relationship exists between decreased grip strength and osteoporosis risk in people aged 45-60 years. The risk of BMD declining in the lumbar spine was associated with reduced grip strength.

Causal relationships between body mass index, low-density lipoprotein and bone mineral density: Univariable and multivariable Mendelian randomization

SUMMARY

Utilizing the Mendelian randomization technique, this research clarifies the putative causal relationship between body mass index (BMI) andbone mineral density (BMD), and the mediating role of low-density lipoprotein (LDL). The implications of these findings present promising opportunities for enhancing our understanding of complex bone-related characteristics and disorders, offering potential directions for treatment and intervention.

OBJECTIVE

The objective of this study is to examine the correlation between BMI and BMD, while exploring the intermediary role of LDL in mediating the causal impact of BMI on BMD outcomes via Mendelian randomization.

METHODS

In this study, we employed genome-wide association study (GWAS) data on BMI, LDL, and BMD to conduct a comparative analysis using both univariate and multivariate Mendelian randomization.

RESULTS

Our study employed a two-sample Mendelian randomization design. Considering BMI as the exposure and BMD as the outcome, our results suggest that BMI may function as a potential protective factor for BMD (β = 0.05, 95% CI 1.01 to 1.09, P = 0.01). However, when treating LDL as the exposure and BMD as the outcome, our findings indicate LDL as a risk factor for BMD (β = -0.04, 95% CI 0.92 to 0.99, P = 0.04). In our multivariate Mendelian randomization (MVMR) model, the combined influence of BMI and LDL was used as the exposure for BMD outcomes. The analysis pointed towards a substantial protective effect of LDL on BMD (β = 0.08, 95% CI 0.85 to 0.97, P = 0.006). In the analysis of mediation effects, LDL was found to mediate the relationship between BMI and BMD, and the effect was calculated at (β = 0.05, 95% CI 1.052 to 1.048, P = 0.04).

CONCLUSION

Our findings suggest that BMI may be considered a protective factor for BMD, while LDL may act as a risk factor. Moreover, LDL appears to play a mediatory role in the causal influence of BMI on BMD.

The impact of bone mineral density on the risk of falling: evidence from genetic correlation and Mendelian randomization analysis

BACKGROUND

Falls are the most common consequence of low bone mineral density (BMD). However, due to limitations inherent in observational studies, the causal relationship between the two remains unestablished.

METHODS

This study utilized Mendelian Randomization (MR) analysis to explore the causal relationship between BMD and the risk of falling, incorporating linkage disequilibrium score (LDSC) regression for genetic correlation assessment. The primary method was inverse-variance weighted (IVW), supplemented with sensitivity analyses and the causal analysis using summary effect estimates (CAUSE) to address heterogeneity and pleiotropy biases.

RESULTS

LDSC analysis indicated significant genetic correlations between BMD at various sites and falling risk (rg range: -0.82 to 0.76, all P < 0.05). IVW analysis, with False Discovery Rate (FDR) correction, showed a protective causal effect of total body BMD (OR = 0.85, 95% CI 0.82-0.88, P = 7.63 × 10-17, PFDR = 1.91 × 10-16), femoral neck BMD (OR = 0.81, 95% CI 0.75-0.88, P = 3.33 × 10-7, PFDR = 5.55 × 10-7), lumbar spine BMD (OR = 0.85, 95% CI 0.79-0.91, P = 9.56 × 10-7, PFDR = 1.20 × 10-6), and heel BMD (OR = 0.82, 95% CI 0.79-0.81, P = 1.69 × 10-39, PFDR = 8.45 × 10-39) on falling risk. No causal relationship was found for forearm BMD (OR = 1.02, 95% CI 0.94-1.11, P = 0.64, PFDR = 0.64). Replication datasets and CAUSE analysis provided causal evidence consistent with the main findings.

CONCLUSION

The study established a causal relationship between BMD at four different sites and the risk of falling, highlighting potential areas for targeted prevention strategies.

Elucidating causal relationships of diet-derived circulating antioxidants and the risk of osteoporosis: A Mendelian randomization study

Background

Osteoporosis (OP) is typically diagnosed by evaluating bone mineral density (BMD), and it frequently results in fractures. Here, we investigated the causal relationships between diet-derived circulating antioxidants and the risk of OP using Mendelian randomization (MR).

Methods

Published studies were used to identify instrumental variables related to absolute levels of circulating antioxidants like lycopene, retinol, ascorbate, and β-carotene, as well as antioxidant metabolites such as ascorbate, retinol, α-tocopherol, and γ-tocopherol. Outcome variables included BMD (in femoral neck, lumbar spine, forearm, heel, total body, total body (age over 60), total body (age 45-60), total body (age 30-45), total body (age 15-30), and total body (age 0-15)), fractures (in arm, spine, leg, heel, and osteoporotic fractures), and OP. Inverse variance weighted or Wald ratio was chosen as the main method for MR analysis based on the number of single nucleotide polymorphisms (SNPs). Furthermore, we performed sensitivity analyses to confirm the reliability of the findings.

Results

We found a causal relationship between absolute retinol levels and heel BMD (p = 7.6E-05). The results of fixed effects IVW showed a protective effect of absolute retinol levels against heel BMD, with per 0.1 ln-transformed retinol being associated with a 28% increase in heel BMD (OR: 1.28, 95% CI: 1.13-1.44). In addition, a sex-specific effect of the absolute circulating retinol levels on the heel BMD has been observed in men. No other significant causal relationship was found.

Conclusion

There is a positive causal relationship between absolute retinol levels and heel BMD. The implications of our results should be taken into account in future studies and in the creation of public health policies and OP prevention tactics.

Bone health, cardiovascular disease, and imaging outcomes in UK Biobank: a causal analysis

This study examined the association of estimated heel bone mineral density (eBMD, derived from quantitative ultrasound) with: (1) prevalent and incident cardiovascular diseases (CVDs: ischemic heart disease (IHD), myocardial infarction (MI), heart failure (HF), non-ischemic cardiomyopathy (NICM), arrhythmia), (2) mortality (all-cause, CVD, IHD), and (3) cardiovascular magnetic resonance (CMR) measures of left ventricular and atrial structure and function and aortic distensibility, in the UK Biobank. Clinical outcomes were ascertained using health record linkage over 12.3 yr of prospective follow-up. Two-sample Mendelian randomization (MR) was conducted to assess causal associations between BMD and CMR metrics using genetic instrumental variables identified from published genome-wide association studies. The analysis included 485 257 participants (55% women, mean age 56.5 ± 8.1 yr). Higher heel eBMD was associated with lower odds of all prevalent CVDs considered. The greatest magnitude of effect was seen in association with HF and NICM, where 1-SD increase in eBMD was associated with 15% lower odds of HF and 16% lower odds of NICM. Association between eBMD and incident IHD and MI was non-significant; the strongest relationship was with incident HF (SHR: 0.90 [95% CI, 0.89-0.92]). Higher eBMD was associated with a decreased risk in all-cause, CVD, and IHD mortality, in the fully adjusted model. Higher eBMD was associated with greater aortic distensibility; associations with other CMR metrics were null. Higher heel eBMD is linked to reduced risk of a range of prevalent and incident CVD and mortality outcomes. Although observational analyses suggest associations between higher eBMD and greater aortic compliance, MR analysis did not support a causal relationship between genetically predicted BMD and CMR phenotypes. These findings support the notion that bone-cardiovascular associations reflect shared risk factors/mechanisms rather than direct causal pathways.

Genetic Evaluation for Monogenic Disorders of Low Bone Mass and Increased Bone Fragility: What Clinicians Need to Know

PURPOSE OF REVIEW

The purpose of this review is to outline the principles of clinical genetic testing and to provide practical guidance to clinicians in navigating genetic testing for patients with suspected monogenic forms of osteoporosis.

RECENT FINDINGS

Heritability assessments and genome-wide association studies have clearly shown the significant contributions of genetic variations to the pathogenesis of osteoporosis. Currently, over 50 monogenic disorders that present primarily with low bone mass and increased risk of fractures have been described. The widespread availability of clinical genetic testing offers a valuable opportunity to correctly diagnose individuals with monogenic forms of osteoporosis, thus instituting appropriate surveillance and treatment. Clinical genetic testing may identify the appropriate diagnosis in a subset of patients with low bone mass, multiple or unusual fractures, and severe or early-onset osteoporosis, and thus clinicians should be aware of how to incorporate such testing into their clinical practices.

Osteoporosis and coronary heart disease: a bi-directional Mendelian randomization study

Background

Osteoporosis (OP) and cardiovascular disease (CVD) are major global public health issues, especially exacerbated by the challenges of an aging population. As these problems intensify, the associated burden on global health is expected to increase significantly. Despite extensive epidemiological investigations into the potential association between OP and CVD, establishing a clear causal relationship remains elusive.

Methods

Instrumental variables were selected from summary statistics of the IEU GWAS database. Five different components of BMD (heel BMD, LS BMD, FA BMD, FN BMD, and TB BMD) were used as OP phenotypes. CHD, MI, and stroke were selected to represent CVD. Multiple analysis methods were used to evaluate the causal relationship between CVD and OP comprehensively. In addition, sensitivity analyses(Cochran's Q test, MR-Egger intercept test, and "leave one out" analysis) were performed to verify the reliability of the results.

Results

The MR showed a significant causal relationship between CHD on heel BMD and TB BMD; in the reverse analysis, there was no evidence that OP has a significant causal effect on CVD. The reliability of the results was confirmed through sensitivity analysis.

Conclusion

The study results revealed that CHD was causally associated with Heel BMD and TB BMD, while in the reverse MR analysis, the causal relationship between OP and CVD was not supported. This result posits CHD as a potential etiological factor for OP and prompts that routine bone density assessment at traditional sites (forearm, femoral neck, lumbar spine) using DAX may inadequately discern underlying osteoporosis issues in CHD patients. The recommendation is to synergistically incorporate heel ultrasound or DAX for total body bone density examinations, ensuring clinical diagnostics are both precise and reliable. Moreover, these findings provide valuable insights for public health, contributing to the development of pertinent prevention and treatment strategies.

Association between BMD and coronary artery calcification: an observational and Mendelian randomization study

Observational studies have reported inconsistent associations between bone mineral density (BMD) and coronary artery calcification (CAC). We examined the observational association of BMD with CAC in 2 large population-based studies and evaluated the evidence for a potential causal relation between BMD and CAC using polygenic risk scores (PRS), 1- and 2-sample Mendelian randomization (MR) approaches. Our study populations comprised 1414 individuals (mean age 69.9 yr, 52.0% women) from the Rotterdam Study and 2233 individuals (mean age 56.5 yr, 50.9% women) from the Framingham Heart Study with complete information on CAC and BMD measurements at the total body (TB-), lumbar spine (LS-), and femoral neck (FN-). We used linear regression models to evaluate the observational association between BMD and CAC. Subsequently, we compared the mean CAC across PRSBMD quintile groups at different skeletal sites. In addition, we used the 2-stage least squares regression and the inverse variance weighted (IVW) model as primary methods for 1- and 2-sample MR to test evidence for a potentially causal association. We did not observe robust associations between measured BMD levels and CAC. These results were consistent with a uniform random distribution of mean CAC across PRSBMD quintile groups (P-value > .05). Moreover, neither 1- nor 2-sample MR supported the possible causal association between BMD and CAC. Our results do not support the contention that lower BMD is (causally) associated with an increased CAC risk. These findings suggest that previously reported epidemiological associations of BMD with CAC are likely explained by unmeasured confounders or shared etiology, rather than by causal pathways underlying both osteoporosis and vascular calcification processes.

The Impact of Non-alcohol Fatty Liver Disease on Bone Mineral Density is Mediated by Sclerostin by Mendelian Randomization Study

Non-alcoholic fatty liver disease (NAFLD) has been found to be associated with osteoporosis (OP) in observational studies. However, the precise causal relationship between NAFLD and OP remains unclear. Here, we used Mendelian randomization (MR) to explore the causal relationship. We selected NAFLD-related single-nucleotide polymorphisms from a genome-wide meta-analysis (8434 cases and 434,770 controls) as instrumental variants. We used inverse variance weighted analysis for the primary MR analysis. Furthermore, we used similar methodologies in parallel investigations of other chronic liver diseases (CLDs). We performed sensitivity analyses to ensure the reliability of the results. We observed a causality between NAFLD and forearm bone mineral density (FABMD) (beta-estimate [β]: - 0.212; p-value: 0.034). We also found that sclerostin can act as a mediator to influence the NAFLD and FABMD pathways to form a mediated MR network (mediated proportion = 8.8%). We also identified indications of causal relationships between other CLDs and OP. However, we were unable to establish any associated mediators. Notably, our analyses did not yield any evidence of pleiotropy. Our findings have implications in the development of preventive and interventional measures aimed at managing low bone mineral density in patients with NAFLD.

Environmental perchlorate, thiocyanate, and nitrate exposures and bone mineral density: a national cross-sectional study in the US adults

Associations of perchlorate, thiocyanate, and nitrate exposures with bone mineral density (BMD) in adults have not previously been studied. This study aimed to estimate the associations of individual and concurrent exposure of the three chemicals with adult BMD. Based on National Health and Nutrition Examination Survey (NHANES, 2011-2018), 1618 non-pregnant adults (age ≥ 20 years and 47.0% female) were included in this study. Survey-weighted linear regression models were used to estimate individual urinary perchlorate, thiocyanate, and nitrate concentrations with lumbar spine BMD and total BMD in adults. Then, weighted quantile sum (WQS) regression and Bayesian kernel machine regression (BKMR) models were conducted to evaluate associations of co-occurrence of the three chemicals with adult BMD. In all participants, nitrate exposure was inversely associated with lumbar spine BMD (β =  - 0.054, 95%CI: - 0.097, - 0.010). In stratification analyses, significant inverse associations were observed in female and participants older than 40 years old. In WQS regressions, significant negative associations of the weighted sum of the three chemicals with total and lumbar spine BMD (β =  - 0.014, 95%CI: - 0.021, - 0.007; β =  - 0.011, 95%CI: - 0.019, - 0.004, respectively) were found, and the dominant contributor was nitrate. In the BKMR models, non-linear dose-response associations of nitrate exposure with lumbar spine and total BMD were observed. These findings suggested that environmental perchlorate, thiocyanate, and nitrate exposure may reduce adult BMD and nitrate is the main contributor.

Systemic Inflammatory Regulators Associated with Osteoporosis: A Bidirectional Mendelian Randomization Study

To elucidate the precise upstream and downstream regulatory mechanisms of inflammatory factors in osteoporosis (OP) progression and to establish a causal relationship between inflammatory factors and OP. We conducted bidirectional Mendelian randomization (MR) analyses using data for 41 cytokines obtained from three independent cohorts comprising 8293 Finnish individuals. Estimated bone mineral density (eBMD) data were derived from 426,824 UK Biobank White British individuals (55% female) and fracture data from 416,795 UK Biobank participants of European ancestry. The inverse variance-weighted method was the primary MR analysis approach. We employed other methods as complementary approaches for mutual corroboration. To test for pleiotropy and heterogeneity, we used the MR-Egger regression, MR-pleiotropy residual sum and outlier global test, and the Cochrane Q test. Macrophage inflammatory protein (MIP)-1α and interleukin (IL)-12p70 expression associated negatively and causally with eBMD (β = -0.017 [MIP-1α], β = -0.011 [IL-12p70]). Conversely, tumor necrosis factor-related apoptosis-inducing ligand was associated with a decreased risk of fractures (Odds Ratio: 0.980). Additionally, OP influenced the expression of multiple inflammatory factors, including growth-regulated oncogene-α, interferon-gamma, IL-6, beta nerve growth factor, and IL-2. Finally, we discovered complex bidirectional causal relationships between IL-8, IL-10, and OP. Specific inflammatory factors may contribute to OP development or may be causally affected by OP. We identified a bidirectional causal relationship between certain inflammatory factors and OP. These findings provide new perspectives for early prediction and targeted treatment of OP. Larger cohort studies are necessary in the future to further validate these findings.

Bone health index in the assessment of bone health: The Generation R Study

Bone Health Index (BHI) has been proposed as a useful instrument for assessing bone health in children. However, its relationship with fracture risk remains unknown. We aimed to investigate whether BHI is associated with bone mineral density (BMD) and prevalent fracture odds in children from the Generation R Study. We also implemented genome-wide association study (GWAS) and polygenic score (PGS) approaches to improve our understanding of BHI and its potential. In total, 4150 children (49.4 % boys; aged 9.8 years) with genotyped data and bone assessments were included in this study. BMD was measured across the total body (less head following ISCD guidelines) using a GE-Lunar iDXA densitometer; and BHI was determined from the hand DXA scans using BoneXpert®. Fractures were self-reported collected with home questionnaires. The association of BHI with BMD and fractures was evaluated using linear models corrected for age, sex, ethnicity, height, and weight. We observed a positive correlation between BHI and BMD (ρ = 0.32, p-value<0.0001). Further, every SD decrease in BHI was associated with an 11 % increased risk of prevalent fractures (OR:1.11, 95 % CI 1.00-1.24, p-value = 0.05). Our BHI GWAS identified variants (lead SNP rs1404264-A, p-value = 2.61 × 10-14) mapping to the ING3/CPED1/WNT16 locus. Children in the extreme tails of the BMD PGS presented a difference in BHI values of -0.10 standard deviations (95% CI -0.14 to -0.07; p-value<0.0001). On top of the demonstrated epidemiological association of BHI with both BMD and fracture risk, our results reveal a partially shared biological background between BHI and BMD. These findings highlight the potential value of using BHI to screen children at risk of fracture.

Diabetes and osteoporosis: a two-sample mendelian randomization study

BACKGROUND

The effects on bone mineral density (BMD)/fracture between type 1 (T1D) and type 2 (T2D) diabetes are unknown. Therefore, we aimed to investigate the causal relationship between the two types of diabetes and BMD/fracture using a Mendelian randomization (MR) design.

METHODS

A two-sample MR study was conducted to examine the causal relationship between diabetes and BMD/fracture, with three phenotypes (T1D, T2D, and glycosylated hemoglobin [HbA1c]) of diabetes as exposures and five phenotypes (femoral neck BMD [FN-BMD], lumbar spine BMD [LS-BMD], heel-BMD, total body BMD [TB-BMD], and fracture) as outcomes, combining MR-Egger, weighted median, simple mode, and inverse variance weighted (IVW) sensitivity assessments. Additionally, horizontal pleiotropy was evaluated and corrected using the residual sum and outlier approaches.

RESULTS

The IVW method showed that genetically predicted T1D was negatively associated with TB-BMD (β = -0.018, 95% CI: -0.030, -0.006), while T2D was positively associated with FN-BMD (β = 0.033, 95% CI: 0.003, 0.062), heel-BMD (β = 0.018, 95% CI: 0.006, 0.031), and TB-BMD (β = 0.050, 95% CI: 0.022, 0.079). Further, HbA1c was not associated with the five outcomes (β ranged from - 0.012 to 0.075).

CONCLUSIONS

Our results showed that T1D and T2D have different effects on BMD at the genetic level. BMD decreased in patients with T1D and increased in those with T2D. These findings highlight the complex interplay between diabetes and bone health, suggesting potential age-specific effects and genetic influences. To better understand the mechanisms of bone metabolism in patients with diabetes, further longitudinal studies are required to explain BMD changes in different types of diabetes.

Causal association of genetically determined caffeine intake from tea or coffee with bone health: a two-sample Mendelian randomization study

BACKGROUND

Relationship of caffeine intake and consumption of caffeinated beverages, such as tea and coffee, with bone health remains controversial. This study aimed to evaluate whether genetically determined caffeine intake from tea or coffee has causal effects on overall total body bone mineral density (TB-BMD) and fracture. We also assessed the association with TB-BMD in five age strata.

METHODS

Using two-sample Mendelian randomization approach, summary statistics were retrieved from genome-wide association studies (GWAS)/GWAS meta-analyses of caffeine intake from tea (n = 395 866)/coffee (n = 373 522), TB-BMD (n = 66 628), and fracture (n = 426 795). Inverse variance weighted method was adopted as the main univariable analysis. Multivariable analysis was conducted to evaluate whether the causal effect is independent.

RESULTS

In univariable analysis, genetically determined caffeine intake from tea had positive association with overall TB-BMD (per SD increase in genetically determined caffeine intake, beta of TB-BMD [in SD]: 0.166; 95% confidence interval (CI): 0.006-0.326) and inverse association with fracture (OR = 0.79; 95% CI: 0.654-0.954). Genetically determined caffeine intake from coffee was also positively associated with overall TB-BMD (beta = 0.231; 95% CI: 0.093-0.369). The association remained significant after adjustment for smoking in multivariable analysis. Genetically determined caffeine intake from tea or coffee was both positively associated with TB-BMD in the age strata of 45-60 years, but we lacked evidence of association in other strata.

CONCLUSIONS

Genetically, caffeine intake from tea or coffee may be beneficial to bone health. Due to the ascertainment method of caffeine intake from tea, our study also implied genetically higher tea consumption may improve TB-BMD and lower fracture risk.

The Relationship Between Systemic Lupus Erythematosus and Osteoporosis Based on Different Ethnic Groups: a Two-Sample Mendelian Randomization Analysis

The previous observational studies could not overcome the effects of confounding variables and reverse causality. We aimed to determine whether there is a causal relationship between systemic lupus erythematosus and osteoporosis in East Asian and European populations, respectively, by two-sample Mendelian Randomization analysis. We obtained and downloaded data from publicly available genome-wide association study databases and analyses for East Asian and European populations, including systemic lupus erythematosus (SLE), osteoporosis (OP), multisite bone mineral density (BMD), and OP with fracture. After screening for instrumental single-nucleotide polymorphisms (SNPs) significantly correlated to SLE, the inverse-variance weighted (IVW) method was used for calculating the ratio and 95% confidence interval, besides utilizing MR-Egger, weighted median, and weighted mode to assess the robustness of the primary outcome. Moreover, multiple analyses, including MR-PRESSO, MR-Egger intercept, Cochran's Q test, as well as "leave-one-out" sensitivity, were used for evaluating horizontal pleiotropy, heterogeneity, and stability. Finally, we exchanged exposure and outcome and performed a reverse MR analysis. IVW (OR = 1.05, 95% CI = 1.01-1.09, P = 0.009) indicated a significant positive correlation between genetically predicted SLE and OP in East Asians. Furthermore, neither heterogeneity nor horizontal pleiotropy was observed. In Europe, there was no significant genetically predicted causal relation between SLE and OP. Bi-directional MR analysis showed no reverse causality between SLE and OP. In the East Asian population, genetically predicted SLE may have had a positive causal relationship with OP. In Europe, there is insufficient evidence for a potential causal relation between SLE and OP or BMD and fracture, and the correlations currently observed may be attributed to a variety of confounder variables.

Association between fatty liver index and controlled attenuation parameters as markers of metabolic dysfunction-associated fatty liver disease and bone mineral density: observational and two-sample Mendelian randomization studies

Previously observational studies did not draw a clear conclusion on the association between fatty liver diseases and bone mineral density (BMD). Our large-scale studies revealed that MAFLD and hepatic steatosis had no causal effect on BMD, while some metabolic factors were correlated with BMD. The findings have important implications for the relationship between fatty liver diseases and BMD, and may help direct the clinical management of MAFLD patients who experience osteoporosis and osteopenia.

PURPOSE

Liver and bone are active endocrine organs with several metabolic functions. However, the link between metabolic dysfunction-associated fatty liver disease (MAFLD) and bone mineral density (BMD) is contradictory.

METHODS

Using the UK Biobank and National Health and Nutrition Examination Survey (NHANES) dataset, we investigated the association between MAFLD, steatosis, and BMD in the observational analysis. We performed genome-wide association analysis to identify single-nucleotide polymorphisms associated with MAFLD. Large-scale two-sample Mendelian randomization (TSMR) analyses examined the potential causal relationship between MAFLD, hepatic steatosis, or major comorbid metabolic factors, and BMD.

RESULTS

After adjusting for demographic factors and body mass index, logistic regression analysis demonstrated a significant association between MAFLD and reduced heel BMD. However, this association disappeared after adjusting for additional metabolic factors. MAFLD was not associated with total body, femur neck, and lumbar BMD in the NHANES dataset. Magnetic resonance imaging-measured steatosis did not show significant associations with reduced total body, femur neck, and lumbar BMD in multivariate analysis. TSMR analyses indicated that MAFLD and hepatic steatosis were not associated with BMD. Among all MAFLD-related comorbid factors, overweight and type 2 diabetes showed a causal relationship with increased BMD, while waist circumference and hyperlipidemia had the opposite effect.

CONCLUSION

No causal effect of MAFLD and hepatic steatosis on BMD was observed in this study, while some metabolic factors were correlated with BMD. This has important implications for understanding the relationship between fatty liver disease and BMD, which may help direct the clinical management of MAFLD patients with osteoporosis.

Causal Relationship of Obstructive Sleep Apnea with Bone Mineral Density and the Role of BMI

Background

Observational studies have yielded conflicting evidence concerning the relationships between obstructive sleep apnea (OSA) and bone mineral density (BMD). As the exact causal inferences remain inconclusive, we conducted a two-sample Mendelian randomization (MR) to identify the causal associations between OSA and BMD.

Methods

Single-nucleotide polymorphisms associated with OSA were extracted from the FinnGen study. Summary statistics for 10 BMD measured at different age or skeletal sites were obtained from the publicly available IEU GWAS database. Inverse-variance weighted (IVW) method was chosen as the primary analysis, combined with several sensitivity analyses to evaluate the robustness of results. The study design included two-sample MR and network MR.

Results

Our primary MR analysis revealed that genetically predicted OSA was positively linked to increased forearm BMD (β = 0.24, 95% confidence interval [CI]: 0.06-0.41, p = 0.009) and heel BMD (β=0.10, 95% CI = 0.02-0.18, p = 0.018), while no significant causal relationships were observed between OSA and total body BMD, lumbar spine BMD, or femoral neck BMD (all p > 0.05). Network MR suggests that OSA might act as a mediating factor in the effect of BMI on forearm BMD and heel BMD, with a mediated portion estimated at 73% and 84%, respectively.

Conclusion

Our findings provide support for a causal relationship between genetically predicted OSA and increased forearm BMD and heel BMD. Furthermore, our results suggest that OSA may play a role in mediating the influence of BMI on BMD.

Large-scale circulating proteome association study (CPAS) meta-analysis identifies circulating proteins and pathways predicting incident hip fractures

Hip fractures are associated with significant disability, high cost, and mortality. However, the exact biological mechanisms underlying susceptibility to hip fractures remain incompletely understood. In an exploratory search of the underlying biology as reflected through the circulating proteome, we performed a comprehensive Circulating Proteome Association Study (CPAS) meta-analysis for incident hip fractures. Analyses included 6430 subjects from two prospective cohort studies (Cardiovascular Health Study and Trøndelag Health Study) with circulating proteomics data (aptamer-based 5 K SomaScan version 4.0 assay; 4979 aptamers). Associations between circulating protein levels and incident hip fractures were estimated for each cohort using age and sex-adjusted Cox regression models. Participants experienced 643 incident hip fractures. Compared with the individual studies, inverse-variance weighted meta-analyses yielded more statistically significant associations, identifying 23 aptamers associated with incident hip fractures (conservative Bonferroni correction 0.05/4979, P < 1.0 × 10-5). The aptamers most strongly associated with hip fracture risk corresponded to two proteins of the growth hormone/insulin growth factor system (GHR and IGFBP2), as well as GDF15 and EGFR. High levels of several inflammation-related proteins (CD14, CXCL12, MMP12, ITIH3) were also associated with increased hip fracture risk. Ingenuity pathway analysis identified reduced LXR/RXR activation and increased acute phase response signaling to be overrepresented among those proteins associated with increased hip fracture risk. These analyses identified several circulating proteins and pathways consistently associated with incident hip fractures. These findings underscore the usefulness of the meta-analytic approach for comprehensive CPAS in a similar manner as has previously been observed for large-scale human genetic studies. Future studies should investigate the underlying biology of these potential novel drug targets.

Association of coffee intake with bone mineral density: a Mendelian randomization study

Background

In observational studies, the relationship between coffee intake and bone mineral density (BMD) is contradictory. However, residual confounding tends to bias the results of these studies. Therefore, we used a two-sample Mendelian randomization (MR) approach to further investigate the potential causal relationship between the two.

Methods

Genetic instrumental variables (IVs) associated with coffee intake were derived from genome-wide association studies (GWAS) of the Food Frequency Questionnaire (FFQ) in 428,860 British individuals and matched using phenotypes in PhenoScanner. Summarized data on BMD were obtained from 537,750 participants, including total body BMD (TB-BMD), TB-BMD in five age brackets ≥60, 45-60, 30-45, 15-30, and 0-15 years, and BMD in four body sites: the lumbar spine, the femoral neck, the heel, and the ultradistal forearm. We used inverse variance weighting (IVW) methods as the primary analytical method for causal inference. In addition, several sensitivity analyses (MR-Egger, Weighted median, MR-PRESSO, Cochran's Q test, and Leave-one-out test) were used to test the robustness of the results.

Results

After Bonferroni correction, Coffee intake has a potential positive correlation with total body BMD (effect estimate [Beta]: 0.198, 95% confidence interval [Cl]: 0.05-0.35, P=0.008). In subgroup analyses, coffee intake was potentially positively associated with TB-BMD (45-60, 30-45 years) (Beta: 0.408, 95% Cl: 0.12-0.69, P=0.005; Beta: 0.486, 95% Cl: 0.12-0.85, P=0.010). In addition, a significant positive correlation with heel BMD was also observed (Beta: 0.173, 95% Cl: 0.08-0.27, P=0.002). The results of the sensitivity analysis were generally consistent.

Conclusion

The results of the present study provide genetic evidence for the idea that coffee intake is beneficial for bone density. Further studies are needed to reveal the biological mechanisms and offer solid support for clinical guidelines on osteoporosis prevention.

Integrative high-throughput enhancer surveying and functional verification divulges a YY2-condensed regulatory axis conferring risk for osteoporosis

The precise roles of chromatin organization at osteoporosis risk loci remain largely elusive. Here, we combined chromatin interaction conformation (Hi-C) profiling and self-transcribing active regulatory region sequencing (STARR-seq) to qualify enhancer activities of prioritized osteoporosis-associated single-nucleotide polymorphisms (SNPs). We identified 319 SNPs with biased allelic enhancer activity effect (baaSNPs) that linked to hundreds of candidate target genes through chromatin interactions across 146 loci. Functional characterizations revealed active epigenetic enrichment for baaSNPs and prevailing osteoporosis-relevant regulatory roles for their chromatin interaction genes. Further motif enrichment and network mapping prioritized several putative, key transcription factors (TFs) controlling osteoporosis binding to baaSNPs. Specifically, we selected one top-ranked TF and deciphered that an intronic baaSNP (rs11202530) could allele-preferentially bind to YY2 to augment PAPSS2 expression through chromatin interactions and promote osteoblast differentiation. Our results underline the roles of TF-mediated enhancer-promoter contacts for osteoporosis, which may help to better understand the intricate molecular regulatory mechanisms underlying osteoporosis risk loci.

Intervertebral Disk Degeneration and Bone Mineral Density: A Bidirectional Mendelian Randomization Study

This study aimed to investigate the causal relationship between bone mineral density (BMD) and intervertebral disk degeneration (IVDD) using a two-sample bidirectional Mendelian randomization analysis. Summary-level data from the Genome-Wide Association Study (GWAS) were used. Instrumental variables (IVs) for IVDD were selected from the large-scale Genome-Wide Association Study (GWAS) (20,001 cases and 164,682 controls). Bone mineral density (BMD) at five different sites (heel (n = 426,824), total body (TB) (n = 56,284), forearm (FA) (n = 8143), femoral neck (FN) (n = 32,735), and lumbar spine (LS) (n = 28,498)) was used as a phenotype for OP. Bidirectional causality between IVDD and BMD was assessed using inverse variance weighting (IVW) and other methods. Related sensitivity analyses were performed. Myopia was also analyzed as a negative control result to ensure the validity of IVs. Heel bone mineral density (heel BMD), total body bone mineral density (TB-BMD), femoral neck bone mineral density (FN-BMD), and lumbar spine bone mineral density (LS-BMD) have a direct causal relationship on intervertebral disk degeneration (IVDD) [heel BMD-related analysis: beta = 0.06, p = 0.03; TB-BMD-related analysis: beta = 0.18, p = 8.72E-08; FN-BMD-related analysis: beta = 0.15, p = 4.89E-03; LS-BMD-related analysis: beta = 0.16, p = 1.43E-04]. There was no evidence of a significant causal effect of IVDD on BMD. In conclusion, our study found a significant positive causal effect of lower BMD on IVDD, and we identified significant causal effects of heel, TB-, FN-, and LS-BMD on IVDD, but there was no evidence of a significant causal effect of IVDD on BMD.

Genome-wide association study identifies genetic regulation of oestrone concentrations and association with endometrial cancer risk in postmenopausal women

BACKGROUND

Oestrone, predominantly made in fat, is the main circulating oestrogen and important for target tissue oestradiol production in women after menopause. The present study was undertaken to determine the genetic regulation of blood oestrone, measured with precision, in postmenopausal women and to explore associations between the identified genetic loci and endometrial cancer in a large, independent cohort.

METHODS

A genome-wide association study (GWAS) was undertaken in women aged at least 70 years to identify genetic associations with blood oestrone concentrations measured by liquid chromatography and tandem mass spectrometry. The GWAS included participants from the Sex Hormones in Older Women (SHOW) study, a sub-study of the longitudinal ASPREE (ASPirin in Reducing Events in the Elderly) randomised trial. Of the 6358 women providing a biobank sample at enrolment, 4951 unrelated women of European ancestry, not taking sex hormones, anti-oestrogens, anti-androgens or systemic glucocorticoids were included in the GWAS. Single nucleotide polymorphisms (SNPs) from loci identified below the genome-wide significance threshold were then tested in an independent cohort (the UK Biobank) for association with endometrial cancer risk, using logistic regression and adjusting for age, body mass index (BMI) and the top 10 genetic principal components.

FINDINGS

The median age of the 4951 women included in the GWAS was 75.9 years (range 70-94.8 years). The GWAS identified four independent SNPs associated with oestrone concentrations (p < 5 × 10-8). Among them, the effect (minor) alleles rs34670419-T, rs2846729-T and rs2414098-T were associated with lower oestrone concentrations. Carrying these effect alleles was associated with lower oestrone concentrations in a dose-dependent manner. The effect allele rs56400819-A was associated with higher oestrone concentrations. When applied to UK Biobank, carrier status for rs2414098-T associated with the CYP19A1 gene which encodes the aromatase enzyme required for oestrogen synthesis was significantly associated with lower endometrial cancer risk (adjusted odd ratio [aOR] 0.87 [95% CI 0.82-0.93]; p = 6.69 × 10-5 for women across all ages and aOR 0.89 [95% CI 0.83-0.96]; p = 0.003 for postmenopausal women). None of the models that included age, body mass index (BMI), the top 10 genetic principal components, parity and diabetes mellitus explained more than 7.6% of the variation in risk.

INTERPRETATION

We have shown genetic regulation of oestrone concentrations in postmenopausal women, and that SNPs associated with oestrone were also associated with endometrial cancer risk, independent of BMI, parity and diabetes mellitus. Although the apparent contribution was modest, the biological influence of oestrone concentrations may be greater through conversion to oestradiol in endometrial tissue.

FUNDING

The ASPREE trial was supported by the National Institute on Aging and the National Cancer Institute at the National Institutes of Health (Grant U01AG029824); the National Health and Medical Research Council (NHMRC) of Australia (Grant 34047, 1127060); Monash University (Australia); and the Victorian Cancer Agency (Australia). The ASPREE Healthy Ageing Biobank was funded by the CSIRO (Flagship Grant), the National Cancer Institute (Grant U01 AG029824) and Monash University. This analysis of sex hormones was funded by an NHMRC of Australia Project Grant (No. 1105305). SRD holds an NHMRC Investigator Grant (2016627). PL is supported by a National Heart Foundation Future Leader Fellowship (102604).

Meta-regression of genome-wide association studies to estimate age-varying genetic effects

Fixed-effect meta-analysis has been used to summarize genetic effects on a phenotype across multiple Genome-Wide Association Studies (GWAS) assuming a common underlying genetic effect. Genetic effects may vary with age (or other characteristics), and not allowing for this in a GWAS might lead to bias. Meta-regression models between study heterogeneity and allows effect modification of the genetic effects to be explored. The aim of this study was to explore the use of meta-analysis and meta-regression for estimating age-varying genetic effects on phenotypes. With simulations we compared the performance of meta-regression to fixed-effect and random -effects meta-analyses in estimating (i) main genetic effects and (ii) age-varying genetic effects (SNP by age interactions) from multiple GWAS studies under a range of scenarios. We applied meta-regression on publicly available summary data to estimate the main and age-varying genetic effects of the FTO SNP rs9939609 on Body Mass Index (BMI). Fixed-effect and random-effects meta-analyses accurately estimated genetic effects when these did not change with age. Meta-regression accurately estimated both main genetic effects and age-varying genetic effects. When the number of studies or the age-diversity between studies was low, meta-regression had limited power. In the applied example, each additional minor allele (A) of rs9939609 was inversely associated with BMI at ages 0 to 3, and positively associated at ages 5.5 to 13. Our findings challenge the assumption that genetic effects are consistent across all ages and provide a method for exploring this. GWAS consortia should be encouraged to use meta-regression to explore age-varying genetic effects.

Tea intake and total body bone mineral density of all ages: a Mendelian randomization analysis

Background

There is increasing evidence indicating that tea intake affects bone mineral density levels; however, the causality between tea intake and bone mineral density is inconclusive. This study aimed to assess the causal relationship between tea intake and total body bone mineral density (TB-BMD) through two-sample Mendelian randomization (MR) analysis.

Methods

We conducted a two-sample MR approach to estimate the potential causal effects of tea intake on TB-BMD at all ages in a European population. The analyses were performed using summary statistics obtained for single-nucleotide polymorphisms (SNPs), identified from a genome-wide association meta-analysis of tea intake (N = up to 447,485 individuals) and from the GEnetic Factors for OSteoporosis (GEFOS) Consortium's genome-wide association meta-analysis (N = up to 56,284 individuals), with baseline data collected in 2018 and populations derived from the European ancestry. The association between each SNP and TB-BMD was weighted by its association with tea intake, and estimates were combined mainly using an inverse-variance weighted meta-analysis. In addition, we explored the potential causal effects between green tea intake, herbal tea intake, and TB-BMD.

Results

The MR analysis revealed that genetically determined tea intake exerts a causal impact on TB-BMD, with an odds ratio (OR) of 1.204 (95% CI: 1.062-1.366, p = 0.004), especially in the age group of 45-60 years (OR = 1.360, 95% CI: 1.088-1.700, p = 0.007). No horizontal pleiotropy and heterogeneity were observed. However, there was no causal effect of tea intake on TB-BMD in the age groups of 0-15, 15-30, 30-45, and over 60 years. In the subgroup analysis, when green tea intake was regarded as the exposure factor, no salient associations were found between green tea consumption and TB-BMD (IVW p = 0.368). Similarly, there was also no causal association between herbal tea intake and TB-BMD (IVW p = 0.264).

Conclusion

The findings of this study support the evidence that tea consumption increases bone density and reduces the risk of osteoporosis in the age group of 45-60 years within the European population.

Association of type 2 Diabetes Mellitus and bone mineral density: a two-sample Mendelian randomization study

BACKGROUND

Observational studies have suggested that type 2 Diabetes Mellitus (DM2) is a potentially modifiable risk factor for lower BMD, but the causal relationship is unclear. This study aimed to examine whether the association of DM2 with lower BMD levels was causal by using Mendelian randomization (MR) analyses.

METHODS

We collected genome-wide association study data for DM2 and BMD of total body and different skeletal sites from the IEU database. Subsequently, we performed a two-sample Mendelian randomization analysis using the Two Sample MR package.

RESULTS

We identified a positive association between DM2 risk (61,714 DM2 cases and 596,424 controls) and total BMD, and other skeletal sites BMD, such as femoral neck BMD, ultra-distal forearm BMD and heel BMD. However, non-significant trends were observed for the effects of DM2 on lumbar-spine BMD.

CONCLUSION

In two-sample MR analyses, there was positive causal relationship between DM2 and BMD in both overall samples. In summary, while observational analyses consistently indicate a strong association between DM2 and low BMD, our MR analysis introduces a nuanced perspective. Contrary to the robust association observed in observational studies, our MR analysis suggests a significant link between DM2 and elevated BMD.

Genetic evidence of the causal relationship between chronic liver diseases and musculoskeletal disorders

BACKGROUND

Chronic liver diseases constitute a major global public health burden, posing a substantial threat to patients' daily lives and even survival due to the potential development of musculoskeletal disorders. Although the relationship between chronic liver diseases and musculoskeletal disorders has received extensive attention, their causal relationship has not been comprehensively and systematically investigated.

METHODS

This study aimed to assess the causal relationships between viral hepatitis, primary biliary cholangitis, primary sclerosing cholangitis (PSC), liver cirrhosis, and hepatocellular carcinoma (HCC) with osteoporosis, osteoarthritis, and sarcopenia through bidirectional Mendelian randomization (MR) research. The traits related to osteoporosis and osteoarthritis included both overall and site-specific phenotypes, and the traits linked to sarcopenia involved indicators of muscle mass and function. Random-effect inverse-variance weighted (IVW), weighted median, MR-Egger, and Causal Analysis Using the Summary Effect Estimates were used to evaluate causal effects, with IVW being the main analysis method. To enhance robustness, sensitivity analyses were performed using Cochran's Q test, MR-Egger intercept, MR-PRESSO global test, funnel plots, leave-one-out analyses, and latent causal variable model.

RESULTS

The forward MR analysis indicated that PSC can reduce forearm bone mineral density (beta = - 0.0454, 95% CI - 0.0798 to - 0.0110; P = 0.0098) and increase the risk of overall osteoarthritis (OR = 1.012, 95% CI 1.002-1.022; P = 0.0247), while HCC can decrease grip strength (beta = - 0.0053, 95% CI - 0.008 to - 0.0025; P = 0.0002). The reverse MR analysis did not find significant causal effects of musculoskeletal disorders on chronic liver diseases. Additionally, no heterogeneity or pleiotropy was detected.

CONCLUSIONS

These findings corroborate the causal effects of PSC on osteoporosis and osteoarthritis, as well as the causal impact of HCC on sarcopenia. Thus, the implementation of comprehensive preventive measures is imperative for PSC and HCC patients to mitigate the risk of musculoskeletal disorders, ultimately improving their quality of life.

The osteoporosis susceptibility SNP rs188303909 at 2q14.2 regulates EN1 expression by modulating DNA methylation and E2F6 binding

EN1 encodes a homeodomain-containing transcription factor and is a determinant of bone density and fracture. Previous powerful genome-wide association studies (GWASs) have identified multiple single-nucleotide polymorphisms (SNPs) near EN1 at 2q14.2 locus for osteoporosis, but the causal SNPs and functional mechanisms underlying these associations are poorly understood. The target genes regulated by the transcription factor EN1 are also unclear. In this study, we identified rs188303909, a functional CpG-SNP, as a causal SNP for osteoporosis at 2q14.2 through the integration of functional and epigenomic analyses. Functional experiments demonstrated that unmethylated rs188303909 acted as a strong allele-specific distal enhancer to regulate EN1 expression by modifying the binding of transcription factor E2F6, but rs188303909 methylation attenuated the active effect of E2F6 on EN1 expression. Importantly, transcription factor EN1 could differentially bind osteoporosis GWAS lead SNPs rs4869739-T and rs4355801-G to upregulate CCDC170 and COLEC10 expression, thus promoting bone formation. Our study provided a mechanistic insight into expression regulation of the osteoporosis susceptibility gene EN1, which could be a potential therapeutic target for osteoporosis precision medicine. KEY MESSAGES: CpG-SNP rs188303909 is a causal SNP at the osteoporosis susceptibility locus 2q14.2. Rs188303909 distally regulates EN1 expression by modulating DNA methylation and E2F6 binding. EN1 upregulates CCDC170 and COLEC10 expression through osteoporosis GWAS lead SNPs rs4869739 and rs4355801.

Lipid Metabolism, Methylation Aberrant, and Osteoporosis: A Multi-omics Study Based on Mendelian Randomization

BACKGROUND

Observational studies have shown a causal association between dyslipidemia and osteoporosis, but the genetic causation and complete mechanism of which are uncertain. The disadvantage of previous observational studies is that they are susceptible to confounding factors and bias, that makes it difficult to infer a causal link between those two diseases. Abnormal epigenetic modifications, represented by DNA methylation, are important causes of many diseases. However, there are no studies showing a bridging role for methylation modifications in blood lipid metabolism and osteoporosis.

METHODS

SNPs for lipid profile (Blood VLDL cholesterol (VLDL-C), blood LDL cholesterol (LDL-C), blood HDL cholesterol (HDL-C), blood triglycerides (TG), diagnosed pure hypercholesterolaemia, blood apolipoprotein B (Apo B), blood apolipoprotein A1(Apo A1)), and bone mineral density (BMD) in different body parts (Heel BMD, lumbar BMD, whole-body BMD, femoral neck BMD) were obtained from large meta-analyses of genome-wide association studies as instrumental variables for two-sample Mendelian randomization. Assessment of the genetic effects of lipid profile-associated methylation sites and bone mineral density was carried out using the summary-data-based Mendelian randomization (SMR) method.

RESULTS

Two-sample Mendelian randomization showed that there was a negative causal association between hypercholesterolaemia and heel BMD (p = 0.0103, OR = 0.4590), and total body BMD (p = 0.0002, OR = 0.2826). LDL-C had a negative causal association with heel BMD (p = 8.68E-05, OR = 0.9586). VLDL-C had a negative causal association with heel BMD (p = 0.035, OR = 0.9484), lumbar BMD (p = 0.0316, OR = 0.9356), and total body BMD (p = 0.0035, OR = 0.9484). HDL-C had a negative causal association with heel BMD (p = 1.25E-05, OR = 0.9548), lumbar BMD (p = 0.0129, OR = 0.9358), and total body BMD (p = 0.0399, OR = 0.9644). Apo B had a negative causal association with heel BMD (p = 0.0001, OR = 0.9647). Apo A1 had a negative causal association with heel BMD (p = 0.0132, OR = 0.9746) and lumbar BMD (p = 0.0058, OR = 0.9261). The p-values of all positive results corrected by the FDR method remained significant and sensitivity analysis showed that there was no horizontal pleiotropy in the results despite the heterogeneity in some results. SMR identified 3 methylation sites associated with lipid profiles in the presence of genetic effects on BMD: cg15707428(GREB1), cg16000331(SREBF2), cg14364472(NOTCH1).

CONCLUSION

Our study provides insights into the potential causal links and co-pathogenesis between dyslipidemia and osteoporosis. The genetic effects of dyslipidaemia on osteoporosis may be related to certain aberrant methylation genetic modifications.

Causal effects of osteoporosis on structural changes in specific brain regions: a Mendelian randomization study

Observational studies have reported that osteoporosis is associated with cortical changes in the brain. However, the inherent limitations of observational studies pose challenges in eliminating confounding factors and establishing causal relationships. And previous observational studies have not reported changes in specific brain regions. By employing Mendelian randomization, we have been able to infer a causal relationship between osteoporosis and a reduction in the surficial area (SA) of the brain cortical. This effect is partially mediated by vascular calcification. We found that osteoporosis significantly decreased the SA of global brain cortical (β = -1587.62 mm2, 95%CI: -2645.94 mm2 to -529.32 mm2, P = 0.003) as well as the paracentral gyrus without global weighted (β = - 19.42 mm2, 95%CI: -28.90 mm2 to -9.95 mm2, P = 5.85 × 10-5). Furthermore, we estimated that 42.25% and 47.21% of the aforementioned effects are mediated through vascular calcification, respectively. Osteoporosis leads to a reduction in the SA of the brain cortical, suggesting the presence of the bone-brain axis. Vascular calcification plays a role in mediating this process to a certain extent. These findings establish a theoretical foundation for further investigations into the intricate interplay between bone, blood vessels, and the brain.

Exploring causality between bone mineral density and frailty: A bidirectional Mendelian randomization study

OBJECTIVE

The bidirectional correlation between low bone mineral density (BMD) and frailty, despite its extensive documentation, still lacks a conclusive understanding. The objective of this Mendelian randomization (MR) study is to investigate the bidirectional causal relationship between BMD and frailty.

METHODS

We utilized summary statistics data for BMD at different skeletal sites-including heel BMD (e-BMD, N = 40,613), forearm BMD (FA-BMD, N = 8,143), femoral neck BMD (FN-BMD, N = 32,735), and lumbar spine BMD (LS-BMD, N = 28,489), alongside frailty index (FI, N = 175,226) data in participants of European ancestry. MR analysis in our study was conducted using well-established analytical methods, including inverse variance weighted (IVW), weighted median (WM), and MR-Egger approaches.

RESULTS

We observed negative causal estimates between genetically predicted e-BMD (IVW β = - 0.020, 95% confidence interval (CI) = - 0.038, - 0.002, P = 0.029) and FA-BMD (IVW β = -0.035, 95% CI = -0.066, -0.004, P = 0.028) with FI. However, the results did not reach statistical significance after applying the Bonferroni correction, with a significance threshold set at P < 0.0125 (0.05/4). There was no causal effect of FN-BMD (IVW β = - 0.024, 95% CI = -0.052, 0.004, P = 0.088) and LS-BMD (IVW β = - 0.005, 95% CI = -0.034, 0.024, P = 0.749) on FI. In the reverse Mendelian randomization (MR) analysis, we observed no causal effect of FI on BMD at various skeletal sites.

CONCLUSION

Our study provides support for the hypothesis that low BMD may be a potential causal risk factor for frailty, but further research is needed to confirm this relationship. However, our findings did not confirm reverse causality.

Tumor necrosis factor mediates the impact of PM2.5 on bone mineral density: Inflammatory proteome Mendelian randomization and colocalization analyses

To assess the causal effect of particulate matter 2.5 (PM2.5) on human bone mineral density (BMD) and to explore the possible mechanism and proportion mediated by inflammation-related protein. The genetic correlation between PM2.5 and BMD was assessed using the Linkage Disequilibrium Score (LDSC), and the causal effect between PM2.5 and BMD was assessed by two-sample Mendelian randomization (TSMR). A 2-step Mendelian randomization (MR) approach was employed to evaluate the potential role of inflammation-associated protein as the mediator in the causal association between PM2.5 and BMD. The multivariate Mendelian randomization (MVMR) study was designed to perform mediation analyses, exclude possible confounders and calculate the proportion of mediation. Subsequently, we used Bayesian colocalization analysis to consolidate the MR results. Finally, using drug-target MR design, we evaluated the potential repurposing of tumor necrosis factor (TNF) inhibitors for the treatment of osteoporosis (OP). The results of the analyses show that BMD is negatively influenced by PM2.5 (Inverse variance weighted [IVW] beta [β] = -0.288, 95% confidence interval [CI]: -0.534 - -0.042, P < 0.05). PM2.5 has a positive causal association with TNF (IVW β = 1.564, 95% CI: 0.155 - 2.973, P < 0.05) and a negative causal association with protachykinin-1 (TAC-1) (IVW β = -1.654, 95% CI: -3.063 - -0.244, P < 0.05). TNF has a negative causal association with BMD (Wald ratio β = -0.082, 95% CI: -0.165 - 0.000, P < 0.05) and TAC-1 has a positive causal association with BMD (IVW β = 0.042, 95% CI: 0.007 - 0.077, P < 0.05). After adjusting TNF and TAC-1, PM2.5 has no causal association with BMD (IVW β = -0.200, 95% CI: -0.579 - 0.179, P > 0.05). After adjusting PM2.5 and TAC-1, there was still a negative causal association between TNF and BMD (IVW β = -0.089, 95% CI: -0.166 - -0.012, P < 0.05). In the final drug-target MR study, the protective effect of TNF/TNF receptor 1 (TNFR1) inhibition on BMD was observed. For every 10% decrease of circulating C-reactive protein (CRP) achieved by TNF/TNF receptor 1 (TNFR1) blockade, β was 0.540 (95% CI: 0.040-1.040) for BMD. We found a negative causal association between PM2.5 and BMD and that causal association was mediated by TNF. The results of drug-target MR do support TNFR1 as a promising target for OP prevention among the general population.

Causal relationship between glycemic traits and bone mineral density in different age groups and skeletal sites: a Mendelian randomization analysis

INTRODUCTION

Previous research has confirmed that patients with type 2 diabetes mellitus tend to have higher bone mineral density (BMD), but it is unknown whether this pattern holds true for individuals without diabetes. This Mendelian randomization (MR) study aims to investigate the potential causal relationship between various glycemic trait (including fasting glucose, fasting insulin, 2-h postprandial glucose, and glycated hemoglobin) and BMD in non-diabetic individuals. The investigation focuses on different age groups (15-30, 30-45, 45-60, and 60 + years) and various skeletal sites (forearm, lumbar spine, and hip).

MATERIALS AND METHODS

We utilized genome-wide association study data from large population-based cohorts to identify robust instrumental variables for each glycemic traits parameter. Our primary analysis employed the inverse-variance weighted method, with sensitivity analyses conducted using MR-Egger, weighted median, MR-PRESSO, and multivariable MR methods to assess the robustness and potential horizontal pleiotropy of the study results.

RESULTS

Fasting insulin showed a negative modulating relationship on both lumbar spine and forearm. However, these associations were only nominally significant. No significant causal association was observed between blood glucose traits and BMD across the different age groups. The direction of fasting insulin's causal effects on BMD showed inconsistency between genders, with potentially decreased BMD in women with high fasting insulin levels and an increasing trend in BMD in men.

CONCLUSIONS

In the non-diabetic population, currently available evidence does not support a causal relationship between glycemic traits and BMD. However, further investigation is warranted considering the observed gender differences.

Osteoporosis and osteoarthritis: a bi-directional Mendelian randomization study

OBJECTIVE

To investigate the causal relationship between low bone mineral density (BMD) and osteoarthritis (OA) using Mendelian randomization (MR) design.

METHODS

Two-sample bi-directional MR analyses were performed using summary-level information on OA traits from UK Biobank and arcOGEN. Sensitivity analyses including MR-Egger, simple median, weighted median, MR pleiotropy residual sum, and outlier approaches were utilized in conjunction with inverse variance weighting (IVW). Gene ontology (GO) enrichment analyses and expression quantitative trait locus (eQTL) colocalization analyses were used to investigate the potential mechanism and shared genes between osteoporosis (OP) and OA.

RESULTS

The IVW method revealed that genetically predicted low femoral neck BMD was significantly linked with hip (β = 0.105, 95% CI: 0.023-0.188) and knee OA (β = 0.117, 95% CI: 0.049-0.184), but not with other site-specific OA. Genetically predicted low lumber spine BMD was significantly associated with OA at any sites (β = 0.048, 95% CI: 0.011-0.085), knee OA (β = 0.101, 95% CI: 0.045-0.156), and hip OA (β = 0.150, 95% CI: 0.077-0.224). Only hip OA was significantly linked with genetically predicted reduced total bone BMD (β = 0.092, 95% CI: 0.010-0.174). In the reverse MR analyses, no evidence for a causal effect of OA on BMD was found. GO enrichment analysis and eQTL analysis illustrated that DDN and SMAD-3 were the most prominent co-located genes.

CONCLUSIONS

These findings suggested that OP may be causally linked to an increased risk of OA, indicating that measures to raise BMD may be effective in preventing OA. More research is required to determine the underlying processes via which OP causes OA.

Causal Effect of Selenium Levels on Osteoporosis: A Mendelian Randomization Study

Prior research has demonstrated equivocal associations between selenium (Se) concentrations and osteoporosis (OP), yielding inconclusive findings. The purpose of the current study was to examine the potential correlation between Se levels and the risk of OP by using the Mendelian randomization (MR) study design. The genetic variants related to Se levels were obtained from a meta-analysis of a Genome-Wide Association Study (GWAS) conducted on toenail Se levels (n = 4162) and blood Se levels (n = 5477). The data summary for OP and bone mineral density (BMD) was obtained by utilizing the GWAS database. To examine the association between Se levels and BMD and OP, we employed three statistical methods: inverse variance weighted, weighted median, and MR-Egger. The reliability of the analysis was verified by sensitivity testing. All three methods of MR analysis revealed that Se levels had no effect on OP risk. In addition, the sensitivity analysis revealed no heterogeneity or pleiotropy, and the significance of the overall effect remained unaffected by single-nucleotide polymorphisms (SNPs), as determined by the leave-one-out analysis, indicating that our findings are relatively reliable. The results of our study indicate that there is no causal association between Se levels and the risk of OP. However, additional investigation is necessary to ascertain whether there is a potential association between these variables.

Zebrafish as a Model for Osteoporosis: Functional Validations of Genome-Wide Association Studies

PURPOSE OF REVIEW

GWAS, as a largely correlational analysis, requires in vitro or in vivo validation. Zebrafish (Danio rerio) have many advantages for studying the genetics of human diseases. Since gene editing in zebrafish has been highly valuable for studying embryonic skeletal developmental processes that are prenatally or perinatally lethal in mammalian models, we are reviewing pros and cons of this model.

RECENT FINDINGS

The true power for the use of zebrafish is the ease by which the genome can be edited, especially using the CRISPR/Cas9 system. Gene editing, followed by phenotyping, for complex traits such as BMD, is beneficial, but the major physiological differences between the fish and mammals must be considered. Like mammals, zebrafish do have main bone cells; thus, both in vivo stem cell analyses and in vivo imaging are doable. Yet, the "long" bones of fish are peculiar, and their bone cavities do not contain bone marrow. Partial duplication of the zebrafish genome should be taken into account. Overall, small fish toolkit can provide unmatched opportunities for genetic modifications and morphological investigation as a follow-up to human-first discovery.

Genetic and Gene Expression Resources for Osteoporosis and Bone Biology Research

PURPOSE OF REVIEW

The integration of data from multiple genomic assays from humans and non-human model organisms is an effective approach to identify genes involved in skeletal fragility and fracture risk due to osteoporosis and other conditions. This review summarizes genome-wide genetic variation and gene expression data resources relevant to the discovery of genes contributing to skeletal fragility and fracture risk.

RECENT FINDINGS

Genome-wide association studies (GWAS) of osteoporosis-related traits are summarized, in addition to gene expression in bone tissues in humans and non-human organisms, with a focus on rodent models related to skeletal fragility and fracture risk. Gene discovery approaches using these genomic data resources are described. We also describe the Musculoskeletal Knowledge Portal (MSKKP) that integrates much of the available genomic data relevant to fracture risk. The available genomic resources provide a wealth of knowledge and can be analyzed to identify genes related to fracture risk. Genomic resources that would fill particular scientific gaps are discussed.

Causal relationship between breakfast skipping and bone mineral density: a two-sample Mendelian randomized study

Objective

To explore the causal association between breakfast skipping and bone mineral density (BMD) through two-sample Mendelian randomisation (MR) analysis.

Methods

A two-sample MR approach was adopted to explore the causal relationship of breakfast skipping with BMDs (across three skeletal sites and five age groups). Publicly available genome-wide association study summary data were used for MR analysis. We used five methods to estimate the causal associations between breakfast skipping and BMDs: inverse-variance weighting (IVW), MR-Egger, weighted median, simple mode, and weighted mode. IVW was used for the main analysis and the remaining four methods were used as supplementary analyses. The heterogeneity of the MR results was determined using IVW and MR-Egger methods. The pleiotropy of the MR results was determined using MR-Egger intercept. Furthermore, a leave-one-out test was performed to determine whether the MR results were affected by a single nucleotide polymorphism.

Results

With the IVW method, we did not find any causal relationship between breakfast skipping and forearm, femoral neck, and lumbar spine BMD. Subsequently, when we included BMD data stratified by five different age groups in the analysis, the results showed that there was no apparent causal effect between breakfast skipping and age-stratified BMD. This finding was supported by all four supplementary methods (P > 0.05 for all methods). No heterogeneity or horizontal pleiotropy was detected in any of the analyses (P > 0.05). The leave-one-out tests conducted in the analyses did not identify any single nucleotide polymorphism that could have influenced the MR results, indicating the reliability of our findings.

Conclusion

No causal effect was found between breakfast skipping and BMD (across three skeletal sites and five age groups).

An atlas of genetic determinants of forearm fracture

Osteoporotic fracture is among the most common and costly of diseases. While reasonably heritable, its genetic determinants have remained elusive. Forearm fractures are the most common clinically recognized osteoporotic fractures with a relatively high heritability. To establish an atlas of the genetic determinants of forearm fractures, we performed genome-wide association analyses including 100,026 forearm fracture cases. We identified 43 loci, including 26 new fracture loci. Although most fracture loci associated with bone mineral density, we also identified loci that primarily regulate bone quality parameters. Functional studies of one such locus, at TAC4, revealed that Tac4-/- mice have reduced mechanical bone strength. The strongest forearm fracture signal, at WNT16, displayed remarkable bone-site-specificity with no association with hip fractures. Tall stature and low body mass index were identified as new causal risk factors for fractures. The insights from this atlas may improve fracture prediction and enable therapeutic development to prevent fractures.

The role of genetic and epigenetic factors in determining the risk of spinal fragility fractures: new insights in the management of spinal osteoporosis

Osteoporosis predisposes patients to spinal fragility fractures. Imaging plays a key role in the diagnosis and prognostication of these osteoporotic vertebral fractures (OVF). However, the current imaging knowledge base for OVF is lacking sufficient standardisation to enable effective risk prognostication. OVF have been shown to be more prevalent in Caucasian patient cohorts in comparison to the Eastern Asian population. These population-based differences in risk for developing OVF suggest that there could be genetic and epigenetic factors that drive the pathogenesis of osteoporosis, low bone mineral density (BMD) and OVF. Several genetic loci have been associated with a higher vertebral fracture risk, although at varying degrees of significance. The present challenge is clarifying whether these associations are specific to vertebral fractures or osteoporosis more generally. Furthermore, these factors could be exploited for diagnostic interpretation as biomarkers [including novel long non-coding (lnc)RNAs, micro (mi)RNAs and circular (circ)RNAs]. The extent of methylation of genes, alongside post-translational histone modifications, have shown to affect several interlinked pathways that converge on the regulation of bone deposition and resorption, partially through their influence on osteoblast and osteoclast differentiation. Lastly, in addition to biomarkers, several exciting new imaging modalities could add to the established dual-energy X-ray absorptiometry (DXA) method used for BMD assessment. New technologies, and novel sequences within existing imaging modalities, may be able to quantify the quality of bone in addition to the BMD and bone structure; these are making progress through various stages of development from the pre-clinical sphere through to deployment in the clinical setting. In this mini review, we explore the literature to clarify the genetic and epigenetic factors associated with spinal fragility fractures and delineate the causal genes, pathways and interactions which could drive different risk profiles. We also outline the cutting-edge imaging modalities which could transform diagnostic protocols for OVF.

A genome-wide genomic score added to standard recommended stratification tools does not improve the identification of patients with very low bone mineral density

The role of integrating genomic scores (GSs) needs to be assessed. Adding a GS to recommended stratification tools does not improve the prediction of very low bone mineral density. However, we noticed that the GS performed equally or above individual risk factors in discrimination.

PURPOSE

We aimed to investigate whether adding a genomic score (GS) to recommended stratification tools improves the discrimination of participants with very low bone mineral density (BMD).

METHODS

BMD was measured in three thoracic vertebrae using CT. All participants provided information on standard osteoporosis risk factors. GSs and FRAX scores were calculated. Participants were grouped according to mean BMD into very low (<80 mg/cm3), low (80-120 mg/cm3), and normal (>120 mg/cm3) and according to the Bone Health and Osteoporosis Foundation recommendations for BMD testing into an "indication for BMD testing" and "no indication for BMD testing" group. Different models were assessed using the area under the receiver operating characteristics curves (AUC) and reclassification analyses.

RESULTS

In the total cohort (n=1421), the AUC for the GS was 0.57 (95% CI 0.52-0.61) corresponding to AUCs for osteoporosis risk factors. In participants without indication for BMD testing, the AUC was 0.60 (95% CI 0.52-0.69) above or equal to AUCs for osteoporosis risk factors. Adding the GS to a clinical risk factor (CRF) model resulted in AUCs not statistically significant from the CRF model. Using probability cutoff values of 6, 12, and 24%, we found no improved reclassification or risk discrimination using the CRF-GS model compared to the CRF model.

CONCLUSION

Our results suggest adding a GS to a CRF model does not improve prediction. However, we noticed that the GS performed equally or above individual risk factors in discrimination. Clinical risk factors combined showed superior discrimination to individual risk factors and the GS, underlining the value of combined CRFs in routine clinics as a stratification tool.

Genetic Variants in WNT16 and PKD2L1 Locus Affect Heel Ultrasound Bone Stiffness: Analyses from the General Population and Patients Evaluated for Osteoporosis

Osteoporosis, a complex chronic disease with increasing prevalence, is characterised by reduced bone mineral density (BMD) and increased fracture risk. The high heritability of BMD suggests substantial impact of the individual genetic disposition on bone phenotypes and the development of osteoporosis. In the past years, genome-wide association studies (GWAS) identified hundreds of genetic variants associated with BMD or osteoporosis. Here, we analysed 1103 single nucleotide polymorphisms (SNPs), previously identified as associated with estimated BMD (eBMD) in the UK Biobank. We assessed whether these SNPs are related to heel stiffness index obtained by quantitative ultrasound in 5665 adult participants of the Study of Health in Pomerania (SHIP). We confirmed 45 significant associations after correction for multiple testing. Next, we analysed six selected SNPs in 631 patients evaluated for osteoporosis [rs2707518 (CPED1/WNT16), rs3779381 (WNT16), rs115242848 (LOC101927709/EN1), rs10239787 (JAZF1), rs603424 (PKD2L1) and rs6968704 (JAZF1)]. Differences in minor allele frequencies (MAF) of rs2707518 and rs3779381 between SHIP participants (higher MAF) and patients evaluated for osteoporosis (lower MAF) indicated a protective effect of the minor allele on bone integrity. In contrast, differences in MAF of rs603424 indicated a harmful effect. Co-localisation analyses indicated that the rs603424 effect may be mediated via stearoyl-CoA desaturase (SCD) expression, an enzyme highly expressed in adipose tissue with a crucial role in lipogenesis. Taken together, our results support the role of the WNT16 pathway in the regulation of bone properties and indicate a novel causal role of SCD expression in adipose tissue on bone integrity.

Evaluation of the Effect of Sodium-Glucose Cotransporter 2 Inhibition on Fracture Risk: Evidence From Mendelian Randomization and Genetic Association Study

This study aims to evaluate the causal effect of sodium-glucose cotransporter 2 (SGLT2) inhibition on bone mineral density (BMD), osteoporosis, and fracture risk using genetics. Two-sample Mendelian randomization (MR) analyses were performed utilizing two sets of genetic variants as instruments (six and two single-nucleotide polymorphisms [SNPs]) associated with SLC5A2 gene expression and glycated hemoglobin A1c levels. Summary statistics of BMD from the Genetic Factors for Osteoporosis consortium (BMD for total body, n = 66,628; femoral neck, n = 32,735; lumbar spine, n = 28,498; forearm, n = 8143) and osteoporosis (6303 cases, 325,717 controls) and 13 types of fracture (≤17,690 cases, ≤328,382 controls) data from the FinnGen study were obtained. One-sample MR and genetic association analyses were conducted in UK Biobank using the individual-level data of heel BMD (n = 256,286) and incident osteoporosis (13,677 cases, 430,262 controls) and fracture (25,806 cases, 407,081 controls). Using six SNPs as the instrument, genetically proxied SGLT2 inhibition showed little evidence of association with BMD of total body, femoral neck, lumbar spine, and forearm (all p ≥ 0.077). Similar results were observed using two SNPs as instruments. Little evidence was found for the SGLT2 inhibition effect on osteoporosis (all p ≥ 0.112) or any 11 major types of fracture (all p ≥ 0.094), except for a nominal significance for fracture of lower leg (p = 0.049) and shoulder and upper arm (p = 0.029). One-sample MR and genetic association analysis showed that both the weighted genetic risk scores constructed from the six and two SNPs were not causally associated with heel BMD, osteoporosis, and fracture (all p ≥ 0.387). Therefore, this study does not support an effect of genetically proxied SGLT2 inhibition on fracture risk. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Total body bone mineral density and various spinal disorders: a Mendelian randomization study

Introduction

Observational studies have yielded inconsistent findings regarding the correlation between bone mineral density (BMD) and various spinal disorders. To explore the relationship between total-body BMD and various spinal disorders further, we conducted a Mendelian randomization analysis to assess this association.

Methods

Two-sample bidirectional Mendelian randomization (MR) analysis was employed to investigate the association between total-body BMD and various spinal disorders. The inverse-variance weighted (IVW) method was used as the primary effect estimate, and additional methods, including weighted median, MR-Egger, simple mode, and weighted mode, were used to assess the reliability of the results. To examine the robustness of the data further, we conducted a sensitivity analysis using alternative bone-density databases, validating the outcome data.

Results

MR revealed a significant positive association between total-body BMD and the prevalence of spondylosis and spinal stenosis. When total-body BMD was considered as the exposure factor, the analysis demonstrated an increased risk of spinal stenosis (IVW odds ratio [OR] 1.23; 95% confidence interval [CI], 1.14-1.32; P < 0.001) and spondylosis (IVW: OR 1.24; 95%CI, 1.16-1.33; P < 0.001). Similarly, when focusing solely on heel BMD as the exposure factor, we found a positive correlation with the development of both spinal stenosis (IVW OR 1.13, 95%CI, 1.05-1.21; P < 0.001) and spondylosis (IVW OR 1.10, 95%CI, 1.03-1.18; P = 0.0048). However, no significant associations were found between total-body BMD and other spinal disorders, including spinal instability, spondylolisthesis/spondylolysis, and scoliosis (P > 0.05).

Conclusion

This study verified an association of total-body BMD with spinal stenosis and with spondylosis. Our results imply that when an increasing trend in BMD is detected during patient examinations and if the patient complains of numbness and pain, the potential occurrence of conditions such as spondylosis or spinal stenosis should be investigated and treated appropriately.

The associations of alcoholic liver disease and nonalcoholic fatty liver disease with bone mineral density and the mediation of serum 25-Hydroxyvitamin D: A bidirectional and two-step Mendelian randomization

BACKGROUND

Reduced bone mineral density (BMD) and osteoporosis are common in chronic liver diseases. However, the causal effect of alcoholic liver disease (ALD) and non-alcoholic fatty liver disease (NAFLD) on BMD remains uncertain.

OBJECTIVES

This study uses a two-sample Mendelian randomization (MR) design to evaluate the genetically predicted effect of ALD and NAFLD on BMDs using summary data from publically available genome-wide association studies (GWASs).

METHODS

The GWAS summary statistics of ALD (1416 cases and 213,592 controls) and NAFLD (894 cases and 217,898 controls) were obtained from the FinnGen consortium. BMDs of four sites (total body, n = 56,284; femoral neck, n = 32,735; lumbar spine, n = 28,498; forearm, n = 8143) were from the GEnetic Factors for OSteoporosis Consortium. Data for alcohol consumption (n = 112,117) and smoking (n = 33,299) and serum 25-Hydroxyvitamin D (25-OHD) level (n = 417,580) were from UK-biobank. We first performed univariate MR analysis with the Inverse Variance Weighted (IVW) method as the primary analysis to investigate the genetically predicted effect of ALD or NAFLD on BMD. Then, multivariate MR and mediation analysis were performed to identify whether the effect was mediated by alcohol consumption, smoking, or serum 25-OHD level.

RESULTS

The MR results suggested a robust genetically predicted effect of ALD on reduced BMD in the femoral neck (FN-BMD) (IVW beta = -0.0288; 95% CI: -0.0488, -0.00871; P = 0.00494) but not the other three sites. Serum 25-OHD level exhibited a significant mediating effect on the association between ALD and reduced FN-BMD albeit the proportion of mediation was mild (2.21%). No significant effects of NAFLD, alcohol consumption, or smoking on BMD in four sites, or reverse effect of BMD on ALD or NAFLD were detected.

CONCLUSION

Our findings confirm the genetically predicted effect of ALD on reduced FN-BMD, and highlight the importance of periodic BMD and serum 25-OHD monitoring and vitamin D supplementation as needed in patients with ALD. Future research is required to validate our results and investigate the probable underlying mechanisms.