Causal impact of DNA methylation on refracture in elderly individuals with osteoporosis - a prospective cohort study.

BACKGROUND: Osteoporotic vertebral compression fractures (OVCF) in the elderly increase refracture risk post-surgery, leading to higher mortality rates. Genome-wide association studies (GWAS) have identified susceptibility genes for osteoporosis, but the phenotypic variance explained by these genes has been limited, indicating the need to explore additional causal factors. Epigenetic modifications, such as DNA methylation, may influence osteoporosis and refracture risk. However, prospective cohorts for assessing epigenetic alterations in Chinese elderly patients are lacking. Here, we propose to conduct a prospective cohort study to investigate the causal network of DNA polymorphisms, DNA methylation, and environmental factors on the development of osteoporosis and the risk of refracture. METHODS: We will collect vertebral and peripheral blood from 500 elderly OVCF patients undergoing surgery, extract DNA, and generate whole genome genotype data and DNA methylation data. Observation indicators will be collected and combined with one-year follow-up data. A healthy control group will be selected from a natural population cohort. Epigenome-wide association studies (EWAS) of osteoporosis and bone mineral density will be conducted. Differential methylation analysis will compare candidate gene methylation patterns in patients with and without refracture. Multi-omics prediction models using genetic variants and DNA methylation sites will be built to predict OVCF risk. DISCUSSION: This study will be the first large-scale population-based study of osteoporosis and bone mineral density phenotypes based on genome-wide data, multi-time point methylation data, and phenotype data. By analyzing methylation changes related to osteoporosis and bone mineral density in OVCF patients, the study will explore the feasibility of DNA methylation in evaluating postoperative osteoporosis intervention effects. The findings may identify new molecular markers for effective anti-osteoporosis treatment and inform individualized prevention and treatment strategies. TRIAL REGISTRATION: chictr.org.cn ChiCTR2200065316, 02/11/2022.

Muscle-origin creatinine-cystatin C ratio is an osteoporosis marker in individuals with normal renal function: evidence from observational and Mendelian randomization analysis.

Background: Creatinine-cystatin C ratio (CCR) has been demonstrated as an objective marker of sarcopenia in clinical conditions but has not been evaluated as an osteoporosis marker in individuals with normal renal function. Methods: We selected 271,831 participants with normal renal function from UK Biobank cohort. Multivariable linear/logistic regression and Cox proportional hazards model were used to investigate the phenotypic relationship between CCR and osteoporosis in total subjects and gender-stratified subjects. Based on the genome-wide association study (GWAS) data, linkage disequilibrium regression (LDSC) and Mendelian randomization (MR) analysis were performed to reveal the shared genetic correlations and infer the causal effects, respectively. Results: Amongst total subjects and gender-stratified subjects, serum CCR was positively associated with eBMD after adjusting for potential risk factors (all P<0.05). The multivariable logistic regression model showed that the decrease in CCR was associated with a higher risk of osteoporosis/fracture in all models (all P<0.05). In the multivariable Cox regression analysis with adjustment for potential confounders, reduced CCR is associated with the incidence of osteoporosis and fracture in both total subjects and gender-stratified subjects (all P<0.05). A significant non-linear dose-response was observed between CCR and osteoporosis/fracture risk (P non-linearity < 0.05). LDSC found no significant shared genetic effects by them, but PLACO identified 42 pleiotropic SNPs shared by CCR and fracture (P<5×10-8). MR analyses indicated the causal effect from CCR to osteoporosis/fracture. Conclusions: Reduced CCR predicted increased risks of osteoporosis/fracture, and significant causal effects support their associations. These findings indicated that the muscle-origin serum CCR was a potential biomarker to assess the risks of osteoporosis and fracture.

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.

OPN, BSP, and Bone Quality-Structural, Biochemical, and Biomechanical Assessment in OPN-/-, BSP-/-, and DKO Mice.

Osteopontin (OPN) and Bone Sialoprotein (BSP), abundantly expressed by osteoblasts and osteoclasts, appear to have important, partly overlapping functions in bone. In gene-knockout (KO, -/-) models of either protein and their double (D)KO in the same CD1/129sv genetic background, we analyzed the morphology, matrix characteristics, and biomechanical properties of femur bone in 2 and 4 month old, male and female mice. OPN-/- mice display inconsistent, perhaps localized hypermineralization, while the BSP-/- are hypomineralized throughout ages and sexes, and the low mineralization of young DKO mice recovers with age. The higher contribution of primary bone remnants in OPN-/- shafts suggests a slow turnover, while their lower percentage in BSP-/- indicates rapid remodeling, despite FTIR-based evidence in this genotype of a high maturity of the mineralized matrix. In 3-point bending assays, OPN-/- bones consistently display higher Maximal Load, Work to Max. Load and in young mice Ultimate Stress, an intrinsic characteristic of the matrix. Young male and old female BSP-/- also display high Work to Max. Load along with low Ultimate Stress. Principal Component Analysis confirms the major role of morphological traits in mechanical competence, and evidences a grouping of the WT phenotype with the OPN-/- and of BSP-/- with DKO, driven by both structural and matrix parameters, suggesting that the presence or absence of BSP has the most profound effects on skeletal properties. Single or double gene KO of OPN and BSP thus have multiple distinct effects on skeletal phenotypes, confirming their importance in bone biology and their interplay in its regulation.

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.

Identification of Transcripts with Shared Roles in the Pathogenesis of Postmenopausal Osteoporosis and Cardiovascular Disease.

Epidemiological evidence suggests existing comorbidity between postmenopausal osteoporosis (OP) and cardiovascular disease (CVD), but identification of possible shared genes is lacking. The skeletal global transcriptomes were analyzed in trans-iliac bone biopsies (n = 84) from clinically well-characterized postmenopausal women (50 to 86 years) without clinical CVD using microchips and RNA sequencing. One thousand transcripts highly correlated with areal bone mineral density (aBMD) were further analyzed using bioinformatics, and common genes overlapping with CVD and associated biological mechanisms, pathways and functions were identified. Fifty genes (45 mRNAs, 5 miRNAs) were discovered with established roles in oxidative stress, inflammatory response, endothelial function, fibrosis, dyslipidemia and osteoblastogenesis/calcification. These pleiotropic genes with possible CVD comorbidity functions were also present in transcriptomes of microvascular endothelial cells and cardiomyocytes and were differentially expressed between healthy and osteoporotic women with fragility fractures. The results were supported by a genetic pleiotropy-informed conditional False Discovery Rate approach identifying any overlap in single nucleotide polymorphisms (SNPs) within several genes encoding aBMD- and CVD-associated transcripts. The study provides transcriptional and genomic evidence for genes of importance for both BMD regulation and CVD risk in a large collection of postmenopausal bone biopsies. Most of the transcripts identified in the CVD risk categories have no previously recognized roles in OP pathogenesis and provide novel avenues for exploring the mechanistic basis for the biological association between CVD and OP.

Falls, fracture and frailty risk in multiple sclerosis: a Mendelian Randomization study to identify shared genetics.

INTRODUCTION: Patients with multiple sclerosis (MS) commonly present musculoskeletal disorders characterized by lower bone mineral density (BMD) and muscle weakness. However, the underlying etiology remains unclear. Our objective is to identify shared pleiotropic genetic effects and estimate the causal relationship between MS and musculoskeletal disorders. MATERIALS AND METHODS: We conducted linkage disequilibrium score regression (LDSR), colocalization, and Mendelian randomization (MR) analyses using summary statistics from recent large-scale genome-wide association studies (GWAS), encompassing MS, falls, fractures, and frailty. Additional MR analyses explored the causal relationship with musculoskeletal risk factors, such as BMD, lean mass, grip strength, and vitamin D. RESULTS: We observed a moderate genetic correlation between MS and falls (RG = 0.10, P-value = 0.01) but not between MS with fracture or frailty in the LDSR analyses. MR revealed MS had no causal association with fracture and frailty but a moderate association with falls (OR: 1.004, FDR q-value = 0.018). We further performed colocalization analyses using nine SNPs that exhibited significant associations with both MS and falls in MR. Two SNPs (rs7731626 on ANKRD55 and rs701006 on OS9 gene) showed higher posterior probability of colocalization (PP.H4 = 0.927), suggesting potential pleiotropic effects between MS and falls. The nine genes are associated with central nervous system development and inflammation signaling pathways. CONCLUSION: We found potential pleiotropic genetic effects between MS and falls. However, our analysis did not reveal a causal relationship between MS and increased risks of falls, fractures, or frailty. This suggests that the musculoskeletal disorders frequently reported in MS patients in clinical studies are more likely attributed to secondary factors associated with disease progression and treatment, rather than being directly caused by MS itself.

Evolutionary and functional analyses of LRP5 in archaic and extant modern humans.

BACKGROUND: The human lineage has undergone a postcranial skeleton gracilization (i.e. lower bone mass and strength relative to body size) compared to other primates and archaic populations such as the Neanderthals. This gracilization has been traditionally explained by differences in the mechanical load that our ancestors exercised. However, there is growing evidence that gracilization could also be genetically influenced. RESULTS: We have analyzed the LRP5 gene, which is known to be associated with high bone mineral density conditions, from an evolutionary and functional point of view. Taking advantage of the published genomes of archaic Homo populations, our results suggest that this gene has a complex evolutionary history both between archaic and living humans and within living human populations. In particular, we identified the presence of different selective pressures in archaics and extant modern humans, as well as evidence of positive selection in the African and South East Asian populations from the 1000 Genomes Project. Furthermore, we observed a very limited evidence of archaic introgression in this gene (only at three haplotypes of East Asian ancestry out of the 1000 Genomes), compatible with a general erasing of the fingerprint of archaic introgression due to functional differences in archaics compared to extant modern humans. In agreement with this hypothesis, we observed private mutations in the archaic genomes that we experimentally validated as putatively increasing bone mineral density. In particular, four of five archaic missense mutations affecting the first ß-propeller of LRP5 displayed enhanced Wnt pathway activation, of which two also displayed reduced negative regulation. CONCLUSIONS: In summary, these data suggest a genetic component contributing to the understanding of skeletal differences between extant modern humans and archaic Homo populations.

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.

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

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.

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.

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.

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.

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.

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.

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.