Analysis and validation of biomarkers of immune cell-related genes in postmenopausal osteoporosis: An observational study.

Postmenopausal osteoporosis (PMOP) is a common metabolic inflammatory disease. In conditions of estrogen deficiency, chronic activation of the immune system leads to a hypo-inflammatory phenotype and alterations in its cytokine and immune cell profile, although immune cells play an important role in the pathology of osteoporosis, studies on this have been rare. Therefore, it is important to investigate the role of immune cell-related genes in PMOP. PMOP-related datasets were downloaded from the Gene Expression Omnibus database. Immune cells scores between high bone mineral density (BMD) and low BMD samples were assessed based on the single sample gene set enrichment analysis method. Subsequently, weighted gene co-expression network analysis was performed to identify modules highly associated with immune cells and obtain module genes. Differential analysis between high BMD and low BMD was also performed to obtain differentially expressed genes. Module genes are intersected with differentially expressed genes to obtain candidate genes, and functional enrichment analysis was performed. Machine learning methods were used to filter out the signature genes. The receiver operating characteristic (ROC) curves of the signature genes and the nomogram were plotted to determine whether the signature genes can be used as a molecular marker. Gene set enrichment analysis was also performed to explore the potential mechanism of the signature genes. Finally, RNA expression of signature genes was validated in blood samples from PMOP patients and normal control by real-time quantitative polymerase chain reaction. Our study of PMOP patients identified differences in immune cells (activated dendritic cell, CD56 bright natural killer cell, Central memory CD4 T cell, Effector memory CD4 T cell, Mast cell, Natural killer T cell, T follicular helper cell, Type 1 T-helper cell, and Type 17 T-helper cell) between high and low BMD patients. We obtained a total of 73 candidate genes based on modular genes and differential genes, and obtained 5 signature genes by least absolute shrinkage and selection operator and random forest model screening. ROC, principal component analysis, and t-distributed stochastic neighbor embedding down scaling analysis revealed that the 5 signature genes had good discriminatory ability between high and low BMD samples. A logistic regression model was constructed based on 5 signature genes, and both ROC and column line plots indicated that the model accuracy and applicability were good. Five signature genes were found to be associated with proteasome, mitochondria, and lysosome by gene set enrichment analysis. The real-time quantitative polymerase chain reaction results showed that the expression of the signature genes was significantly different between the 2 groups. HIST1H2AG, PYGM, NCKAP1, POMP, and LYPLA1 might play key roles in PMOP and be served as the biomarkers of PMOP.

Genetic influence on within-person longitudinal change in anthropometric traits in the UK Biobank.

The causes of temporal fluctuations in adult traits are poorly understood. Here, we investigate the genetic determinants of within-person trait variability of 8 repeatedly measured anthropometric traits in 50,117 individuals from the UK Biobank. We found that within-person (non-directional) variability had a SNP-based heritability of 2-5% for height, sitting height, body mass index (BMI) and weight (P ≤ 2.4 × 10-3). We also analysed longitudinal trait change and show a loss of both average height and weight beyond about 70 years of age. A variant tracking the Alzheimer's risk APOE- E 4 allele (rs429358) was significantly associated with weight loss ( ß = -0.047 kg per yr, s.e. 0.007, P = 2.2 × 10-11), and using 2-sample Mendelian Randomisation we detected a relationship consistent with causality between decreased lumbar spine bone mineral density and height loss (bxy = 0.011, s.e. 0.003, P = 3.5 × 10-4). Finally, population-level variance quantitative trait loci (vQTL) were consistent with within-person variability for several traits, indicating an overlap between trait variability assessed at the population or individual level. Our findings help elucidate the genetic influence on trait-change within an individual and highlight disease risks associated with these changes.

Interaction between MARK3 (rs11623869), PLCB4 (rs6086746) and GEMIN2 (rs2277458) variants with bone mineral density and serum 25-hidroxivitamin D levels in Mexican Mestizo women.

Introduction: Understanding the genetic factors contributing to variations in bone mineral density (BMD) and vitamin D could provide valuable insights into the pathogenesis of osteoporosis. This study aimed to evaluate the association of single nucleotide variants in MARK3 (rs11623869), PLCB4 (rs6086746), and GEMIN2 (rs2277458) with BMD in Mexican women. Methods: The gene-gene interaction was evaluated in these variants in serum 25(OH)D levels and BMD. A genetic risk score (GRS) was created on the basis of the three genetic variants. Genotyping was performed using predesigned TaqMan assays. Results: A significant association was found between the rs6086746-A variant and BMD at the total hip, femoral neck, and lumbar spine, in women aged 45 years or older. However, no association was observed between the variants rs11623869 and rs2277458. The rs11623869 × rs2277458 interaction was associated with total hip (p=0.002) and femoral neck BMD (p=0.013). Similarly, for vitamin D levels, we observed an interaction between the variants rs6086746 × rs2277458 (p=0.021). GRS revealed a significant association with total hip BMD (p trend=0.003) and femoral neck BMD (p trend=0.006), as well as increased vitamin D levels (p trend=0.0003). These findings provide evidence of the individual and joint effect of the MARK3, PLCB4, and GEMIN2 variants on BMD and serum vitamin D levels in Mexican women. Discussion: This knowledge could help to elucidate the interaction mechanism between BMD-related genetic variants and 25OHD, contributing to the determination of the pathogenesis of osteoporosis and its potential implications during early interventions.

Causal relationship between intervertebral disc degeneration and osteoporosis: a bidirectional two-sample Mendelian randomization study.

Introduction: The relationship between intervertebral disc degeneration (IVDD) and osteoporosis (OP), diagnosed primarily using bone mineral density (BMD), remains unclear so far. The present study, therefore, aimed to investigate the potential relationship between osteoporosis and intervertebral disc degeneration using Mendelian randomization and genome-wide association analyses. Specifically, the impact of bone mineral density on the development of intervertebral disc degeneration was evaluated. Materials and methods: The genome-wide association studies (GWAS) summary data of OP/BMDs and IVDD were collected from the FinnGen consortium, the GEFOS consortium, and MRC-IEU. The relationship between IVDD and OP was then explored using TSMR. The inverse-variance weighted (IVW) method was adopted as the primary effect estimate, and the reliability and stability of the results were validated using various methods, including MR-Egger, weighted median, simple mode, weighted mode, and MR-PRESSO. Results: No significant causal relationship was observed between OP and IVDD (IVW, P > 0.05) or between femoral neck BMD (FA-BMD) and IVDD when OP and FA-BMD were used as exposures. However, increased levels of total body BMD (TB-BMD) and lumbar spine BMD (LS-BMD) were revealed as significant risk factors for IVDD (TB-BMD: IVW, OR = 1.201, 95% CI: 1.123-1.284, P = 8.72 × 10-8; LS-BMD: IVW, OR = 1.179, 95% CI: 1.083-1.284, P = 1.43 × 10-4). Interestingly, both heel BMD (eBMD) and femur neck BMD (FN-BMD) exhibited potential causal relationships (eBMD: IVW, OR = 1.068, 95% CI: 1.008-1.131, P = 0.0248; FN-BMD, IVW, OR = 1.161, 95% CI: 1.041-1.295, P = 0.0074) with the risk of IVDD. The reverse MR analysis revealed no statistically causal impact of IVDD on OP and the level of BMD (P > 0.05). Conclusion: OP and the level of FA-BMD were revealed to have no causal relationship with IVDD. The increased levels of TB-BMD and LS-BMD could promote the occurrence of IVDD. Both eBMD and FN-BMD have potential causal relationships with the risk of IVDD. No significant relationship exists between IVDD and the risk of OP. Further research is warranted to comprehensively comprehend the molecular mechanisms underlying the impact of OP and BMD on IVDD and vice versa.

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.

Genetically determined type 1 diabetes mellitus and risk of osteoporosis.

BACKGROUND: Observational evidence suggests that type 1 diabetes mellitus (T1DM) is associated with the risk of osteoporosis (OP). Nevertheless, it is not apparent whether these correlations indicate a causal relationship. To elucidate the causal relationship, a two-sample Mendelian randomization (MR) analysis was performed. METHODS: T1DM data was obtained from the large genome-wide association study (GWAS), in which 6683 cases and 12,173 controls from 12 European cohorts were involved. Bone mineral density (BMD) samples at four sites were extracted from the GEnetic Factors for OSteoporosis (GEFOS) consortium, including forearm (FA) (n = 8,143), femoral neck (FN) (n = 32,735), lumbar spine (LS) (n = 28,498), and heel (eBMD) (n = 426,824). The former three samples were from mixed populations and the last one was from European. Inverse variance weighting, MR-Egger, and weighted median tests were used to test the causal relationship between T1DM and OP. A series of sensitivity analyses were then conducted to verify the robustness of the results. RESULTS: Twenty-three independent SNPs were associated with FN-BMD and LS-BMD, twenty-seven were associated with FA-BMD, and thirty-one were associated with eBMD. Inverse variance-weighted estimates indicated a causal effect of T1DM on FN-BMD (odds ratio (OR) =1.033, 95 % confidence interval (CI): 1.012-1.054, p = 0.002) and LS-BMD (OR = 1.032, 95 % CI: 1.005-1.060, p = 0.022) on OP risk. Other MR methods, including weighted median and MR-Egger, calculated consistent trends. While no significant causation was found between T1DM and the other sites (FA-BMD: OR = 1.008, 95 % CI: 0.975-1.043, p = 0.632; eBMD: OR = 0.993, 95 % CI: 0.985-1.001, p = 0.106). No significant heterogeneity (except for eBMD) or horizontal pleiotropy was found for instrumental variables, suggesting these results were reliable and robust. CONCLUSIONS: This study shows a causal relationship between T1DM and the risk of some sites of OP (FN-BMD, LS-BMD), allowing for continued research to discover the clinical and experimental mechanisms of T1DM and OP. It also contributes to the recommendation if patients with T1DM need targeted care to promote bone health and timely prevention of osteoporosis.

Metabolic and Skeletal Characterization of the KK/Ay Mouse Model-A Polygenic Mutation Model of Obese Type 2 Diabetes.

Type 2 diabetes (T2D) increases fracture incidence and fracture-related mortality rates (KK.Cg-Ay/J. The Jackson Laboratory; Available from: https://www.jax.org/strain/002468 ). While numerous mouse models for T2D exist, few effectively stimulate persistent hyperglycemia in both sexes, and even fewer are suitable for bone studies. Commonly used models like db/db and ob/ob have altered leptin pathways, confounding bone-related findings since leptin regulates bone properties (Fajardo et al. in Journal of Bone and Mineral Research 29(5): 1025-1040, 2014). The Yellow Kuo Kondo (KK/Ay) mouse, a polygenic mutation model of T2D, is able to produce a consistent diabetic state in both sexes and addresses the lack of a suitable model of T2D for bone studies. The diabetic state of KK/Ay stems from a mutation in the agouti gene, responsible for coat color in mice. This mutation induces ectopic gene expression across various tissue types, resulting in diabetic mice with yellow fur coats (Moussa and Claycombe in Obesity Research 7(5): 506-514, 1999). Male and female KK/Ay mice exhibited persistent hyperglycemia, defining them as diabetic with blood glucose (BG) levels consistently exceeding 300 mg/dL. Notably, male control mice in this study were also diabetic, presenting a significant limitation. Nevertheless, male and female KK/Ay mice showed significantly elevated BG levels, HbA1c, and serum insulin concentration when compared to the non-diabetic female control mice. Early stages of T2D are characterized by hyperglycemia and hyperinsulinemia resulting from cellular insulin resistance, whereas later stages may feature hypoinsulinemia due to ß-cell apoptosis (Banday et al. Avicenna Journal of Medicine 10(04): 174-188, 2020 and Klein et al. Cell Metabolism 34(1): 11-20, 2022). The observed hyperglycemia, hyperinsulinemia, and the absence of differences in ß-cell mass suggest that KK/Ay mice in this study are modeling the earlier stages of T2D. While compromised bone microarchitecture was observed in this study, older KK/Ay mice, representing more advanced stages of T2D, might exhibit more pronounced skeletal manifestations. Compared to the control group, the femora of KK/Ay mice had higher cortical area and cortical thickness, and improved trabecular properties which would typically be indicative of greater bone strength. However, KK/Ay mice displayed lower cortical tissue mineral density in both sexes and increased cortical porosity in females. Fracture instability toughness of the femora was lower in KK/Ay mice overall compared to controls. These findings indicate that decreased mechanical integrity noted in the femora of KK/Ay mice was likely due to overall bone quality being compromised.

SharePro: an accurate and efficient genetic colocalization method accounting for multiple causal signals.

MOTIVATION: Colocalization analysis is commonly used to assess whether two or more traits share the same genetic signals identified in genome-wide association studies (GWAS), and is important for prioritizing targets for functional follow-up of GWAS results. Existing colocalization methods can have suboptimal performance when there are multiple causal variants in one genomic locus. RESULTS: We propose SharePro to extend the COLOC framework for colocalization analysis. SharePro integrates linkage disequilibrium (LD) modeling and colocalization assessment by grouping correlated variants into effect groups. With an efficient variational inference algorithm, posterior colocalization probabilities can be accurately estimated. In simulation studies, SharePro demonstrated increased power with a well-controlled false positive rate at a low computational cost. Compared to existing methods, SharePro provided stronger and more consistent colocalization evidence for known lipid-lowering drug target proteins and their corresponding lipid traits. Through an additional challenging case of the colocalization analysis of the circulating abundance of R-spondin 3 GWAS and estimated bone mineral density GWAS, we demonstrated the utility of SharePro in identifying biologically plausible colocalized signals. AVAILABILITY AND IMPLEMENTATION: SharePro for colocalization analysis is written in Python and openly available at https://github.com/zhwm/SharePro_coloc.

Comprehensive analysis of circRNA expression profiles in postmenopausal women differing in bone mineral density.

Postmenopausal osteoporosis (PMOP) seriously endangers the bone health of older women. Although there are currently indicators to diagnose PMOP, early diagnostic biomarkers are lacking. Circular ribonucleic acid (circRNA) has a stable structure, regulates gene expression, participates in the pathological process of disease, and has the potential to become a biomarker. The purpose of this study was to investigate circRNAs that could be used to predict patients with early PMOP. Ribonucleic acid (RNA) sequencing was performed on peripheral blood leukocytes from 15 female patients to identify differential circRNAs between different groups. Using bioinformatics analysis, enrichment analysis was performed to discover relevant functions and pathways. CircRNA-micro ribonucleic acid (miRNA) interaction analysis and messenger ribonucleic acid (mRNA) prediction and network construction help us to understand the relationship between circRNA, miRNA, and mRNA. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to validate the gene expression of candidate circRNAs. We screened out 2 co-expressed differential circRNAs, namely hsa_circ_0060849 and hsa_circ_0001394. By analyzing the regulatory network, a total of 54 miRNAs and 57 osteoporosis-related mRNAs were identified, which, as potential downstream target genes of hsa_circ_0060849 and hsa_circ_0001394, may play a key role in the occurrence and development of PMOP. The occurrence and development of PMOP is regulated by circRNAs, and hsa_circ_0060849 and hsa_circ_0001394 can be used as new diagnostic markers and therapeutic targets for early PMOP.

GBA1 as a risk gene for osteoporosis in the specific populations and its role in the development of Gaucher disease.

BACKGROUND: Osteoporosis and its primary complication, fragility fractures, contribute to substantial global morbidity and mortality. Gaucher disease (GD) is caused by glucocerebrosidase (GBA1) deficiency, leading to skeletal complications. This study aimed to investigate the impact of the GBA1 gene on osteoporosis progression in GD patients and the specific populations. METHODS: We selected 8115 patients with osteoporosis (T-score ≤ - 2.5) and 55,942 healthy individuals (T-score > - 1) from a clinical database (N = 95,223). Monocytes from GD patients were evaluated in relation to endoplasmic reticulum (ER) stress, inflammasome activation, and osteoclastogenesis. An in vitro model of GD patient's cells treated with adeno-associated virus 9 (AAV9)-GBA1 to assess GBA1 enzyme activity, chitotriosidase activity, ER stress, and osteoclast differentiation. Longitudinal dual-energy X-ray absorptiometry (DXA) data tracking bone density in patients with Gaucher disease (GD) undergoing enzyme replacement therapy (ERT) over an extended period. RESULTS: The GBA1 gene variant rs11264345 was significantly associated [P < 0.002, Odds Ratio (OR) = 1.06] with an increased risk of bone disease. Upregulation of Calnexin, NOD-, LRR- and pyrin domain-containing protein 3 (NLRP3) and Apoptosis-associated speck-like protein containing a C-terminal caspase recruitment domain (ASC) was positively associated with osteoclastogenesis in patients with GD. In vitro AAV9-GBA1 treatment of GD patient cells led to enhanced GBA1 enzyme activity, reduced chitotriosidase activity, diminished ER stress, and decreased osteoclast differentiation. Long-term bone density data suggests that initiating ERT earlier in GD leads to greater improvements in bone density. CONCLUSIONS: Elevated ER stress and inflammasome activation are indicative of osteoporosis development, suggesting the need for clinical monitoring of patients with GD. Furthermore, disease-associated variant in the GBA1 gene may constitute a risk factor predisposing specific populations to osteoporosis.

Identification of bone mineral density associated genes with shared genetic architectures across multiple tissues: Functional insights for EPDR1, PKDCC, and SPTBN1.

Recent studies suggest a shared genetic architecture between muscle and bone, yet the underlying molecular mechanisms remain elusive. This study aims to identify the functionally annotated genes with shared genetic architecture between muscle and bone using the most up-to-date genome-wide association study (GWAS) summary statistics from bone mineral density (BMD) and fracture-related genetic variants. We employed an advanced statistical functional mapping method to investigate shared genetic architecture between muscle and bone, focusing on genes highly expressed in muscle tissue. Our analysis identified three genes, EPDR1, PKDCC, and SPTBN1, which are highly expressed in muscle tissue and previously unlinked to bone metabolism. About 90% and 85% of filtered Single-Nucleotide Polymorphisms were in the intronic and intergenic regions for the threshold at P≤5×10-8 and P≤5×10-100, respectively. EPDR1 was highly expressed in multiple tissues, including muscles, adrenal glands, blood vessels, and the thyroid. SPTBN1 was highly expressed in all 30 tissue types except blood, while PKDCC was highly expressed in all 30 tissue types except the brain, pancreas, and skin. Our study provides a framework for using GWAS findings to highlight functional evidence of crosstalk between multiple tissues based on shared genetic architecture between muscle and bone. Further research should focus on functional validation, multi-omics data integration, gene-environment interactions, and clinical relevance in musculoskeletal disorders.

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.

Potential causal association between leisure sedentary behaviors and osteoporosis: A two-sample Mendelian randomization analysis.

Previous observational studies have observed a correlation between sedentary behavior and osteoporosis. However, conclusions from these studies have been contradictory. To explore the potential causal relationship between sedentary behavior and osteoporosis, we conducted a Mendelian randomization analysis. A two-sample Mendelian randomization was adopted to explore the causal relationship of leisure sedentary behavior with osteoporosis. We employed 5 methods to estimate the causal associations between leisure sedentary behavior and osteoporosis. Univariable Mendelian randomization results provided evidence for the causal relationship of the time spent on computer-use with the bone mineral density estimated by heel quantitative ultrasound (eBMD) (inverse variance weighted [IVW]: ß (95% confidence interval [CI]) - 0.150 (-0.270 to -0.031), P = .013; weighted median: ß (95%CI) - 0.195 (-0.336 to -0.055), P = .006). Similar associations were observed in the driving forearm bone mineral density (FABMD) (IVW: ß (95%CI) - 0.933 (-1.860 to -0.007), P = .048) and driving lumbar spine bone mineral density (IVW: ß (95%CI) - 0.649 (-1.175 to -0.124), P = .015). However, we did not find a significant causal relationship between the time spent on watching TV and bone mineral density. Research showed that there was a causal relationship between the time spent on computer use and driving time and eBMD, FABMD, and lumbar spine bone mineral density.

Associations of residential greenness with bone mineral density and osteoporosis: the modifying effect of genetic susceptibility.

OBJECTIVES: To investigate the associations of residential greenness with bone mineral density and incident osteoporosis, and further evaluate the potential modifying effect of genetic susceptibility. METHODS: We used the Normalised Difference Vegetation Index (NDVI) at various buffer distances, including 300 m (NDVI300m), 500 m (NDVI500m), 1000 m (NDVI1000m) and 1500 m (NDVI1500m), to serve as indicators of greenness. We fitted linear regression, logistic regression and Cox proportional hazard models to assess the associations of residential greenness with estimated bone mineral density (eBMD), prevalent osteoporosis and incident osteoporosis, respectively. With the Polygenic Risk Score (PRS) for osteoporosis, we further assessed the joint effects of genetic risk and greenness on the risk of osteoporosis. We conducted causal mediation analyses to explore potential mediators. RESULTS: Each IQR increase in NDVI300m was associated with 0.0007 (95% CI 0.0002 to 0.0013) increase in eBMD, 6% lower risk of prevalent osteoporosis (OR 0.94; 95% CI 0.92 to 0.97) and 5% lower risk of incident osteoporosis (HR 0.95; 95% CI 0.93 to 0.98). The joint effects of greenness and PRS on the risk of osteoporosis displayed a clear dose-response pattern. Compared with individuals exposed to low NDVI levels and high genetic risk, those exposed to high NDVI levels and low genetic risk had a 56% (95% CI 51% to 61%) lower risk of osteoporosis. The primary mediators in the association between greenness and incident osteoporosis were identified as PM2.5 and NO2. CONCLUSIONS: Residential greenness was associated with higher bone mineral density and decreased risk of incident osteoporosis.

Femoral neck width genetic risk score is a novel independent risk factor for hip fractures.

Femoral neck width (FNW) derived from DXA scans may provide a useful adjunct to hip fracture prediction. Therefore, we investigated whether FNW is related to hip fracture risk independently of femoral neck bone mineral density (FN-BMD), using a genetic approach. FNW was derived from points automatically placed on the proximal femur using hip DXA scans from 38 150 individuals (mean age 63.8 yr, 48.0% males) in UK Biobank (UKB). Genome-wide association study (GWAS) identified 71 independent genome-wide significant FNW SNPs, comprising genes involved in cartilage differentiation, hedgehog, skeletal development, in contrast to SNPs identified by FN-BMD GWAS which primarily comprised runx1/Wnt signaling genes (MAGMA gene set analyses). FNW and FN-BMD SNPs were used to generate genetic instruments for multivariable Mendelian randomization. Greater genetically determined FNW increased risk of all hip fractures (odds ratio [OR] 1.53; 95% CI, 1.29-1.82 per SD increase) and femoral neck fractures (OR 1.58;1.30-1.92), but not trochanteric or forearm fractures. In contrast, greater genetically determined FN-BMD decreased fracture risk at all 4 sites. FNW and FN-BMD SNPs were also used to generate genetic risk scores (GRSs), which were examined in relation to incident hip fracture in UKB (excluding the FNW GWAS population; n = 338 742, 3222 cases) using a Cox proportional hazards model. FNW GRS was associated with increased risk of all incident hip fractures (HR 1.08;1.05-1.12) and femoral neck fractures (hazard ratio [HR] 1.10;1.06-1.15), but not trochanteric fractures, whereas FN-BMD GRS was associated with reduced risk of all hip fracture types. We conclude that the underlying biology regulating FNW and FN-BMD differs, and that DXA-derived FNW is causally related to hip fractures independently of FN-BMD, adding information beyond FN-BMD for hip fracture prediction. Hence, FNW derived from DXA analyses or a FNW GRS may contribute clinically useful information beyond FN-BMD for hip fracture prediction.

Femoral neck width (FNW) derived from DXA scans may provide useful information about hip fracture prediction, over and above that provided by BMD measurements. Therefore, we investigated whether FNW is related to hip fracture risk independently of BMD, using a genetic approach. FNW was derived from points automatically placed on the hip in DXA scans obtained from 38 150 individuals (mean age 63.8 yr, 48.0% males) in UK Biobank. Seventy-one distinct genetic factors were found to be associated with FNW. Individuals who were predicted by their genes to have greater FNW had a higher risk of hip but not forearm fractures. In contrast, those with greater genetically determined BMD of the femoral neck had a lower risk of both hip and forearm fractures. We conclude that the underlying biology regulating FNW and BMD of the femoral neck differs, and that FNW derived from DXA analyses may contribute clinically useful information beyond BMD for hip fracture prediction.

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.

Decreased bone mineral density, the determinant of osteoporosis, and increased coronary artery calcification are common in people at an advanced age and share some common risk factors. Some studies have reported a higher risk for coronary artery calcification in people with osteoporosis than in people without, whereas others failed to find evidence for this relationship. Recently, Mendelian randomization has emerged as an important epidemiological tool that offers a simple way to distinguish causation, minimizing the confounding present in observational studies, leveraging individual genetic data and the findings from robust genome-wide association studies. We combined data from the participants of both the Rotterdam Study and the Framingham Heart Study, and did not observe sufficient evidence for the association between bone mineral density at different skeletal sites and coronary artery calcification. Also, when using Mendelian randomization, we concluded there was no causal relation between bone deterioration and the build-up of calcium in the coronary arteries. Although more research is needed, we conclude that the associations between decreased bone mineral density and increased coronary artery calcification reported in previous studies are likely attributed to other confounders rather than a causal relationship between these traits.

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.

Association between alcohol intake and bone mineral density: results from the NHANES 2005-2020 and two-sample Mendelian randomization.

We used the data from the NHANES cross-sectional study among 14,113 participants and indicated a positive correlation between alcohol intake frequency and bone mineral density in different body sites. Mendelian randomization was conducted, and no causal relationship is significant between these two variables. The study can provide some suggestions on the daily consumption of alcohol for osteoporosis patients. PURPOSE: The effect of alcohol intake on bone mineral density (BMD) remains unclear. This study explored the association and causality between alcohol intake and BMD. METHODS: Based on the 2005-2020 National Health and Nutrition Examination Survey including 14,113 participants, we conducted co-variate-adjusted multilinear regression analyses to explore the association between alcohol intake levels and spine or femur BMD. To evaluate the causal association between alcohol intake frequency and bone mineral density, the inverse variance weighted approach of two-sample Mendelian randomization (MR) was used with genetic data from the Medical Research Council Integrative Epidemiology Unit (462,346 cases) for alcohol intake frequency and the Genetic Factors for Osteoporosis Consortium (28,496 cases) for lumbar spine and femur neck BMD (32,735 cases). RESULTS: Compared with non-drinkers, total femur BMDs but not total spine BMD increased with daily alcohol intake in males (ß = 3.63*10-2 for mild drinkers, ß = 4.21*10-2 for moderate drinkers, and ß = 4.26*10-2 for heavy drinkers). By contrast, the higher total spine BMD in females was related to higher alcohol intake levels (ß = 2.15*10-2 for mild drinkers, ß = 2.59*10-2 for moderate drinkers, and ß = 3.88*10-2 for heavy drinkers). Regarding the two-sample MR results, no causal relationship was observed between alcohol intake frequency and lumbar spine BMD (odds ratio [OR] = 1.016, P = 0.789) or femur neck BMD (OR = 1.048, P = 0.333). CONCLUSION: This study suggests a positive association between alcohol intake frequency and BMD, although the causal relationship was not significant.

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.

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.