Bioinformatics analysis of key genes, immune infiltration, and risk assessment in low bone mineral density among perimenopausal women: An observational study.

This study aimed to identify hub genes and elucidate the molecular mechanisms underlying low bone mineral density (BMD) in perimenopausal women. R software was used to normalize the dataset and screen the gene set associated with BMD in perimenopausal women from the Gene Expression Omnibus database. Cytoscape software was used to identify 7 critical genes. Gene enrichment analysis and protein interaction was employed to further analyze the core genes, and the CIBERSORT deconvolution algorithm was used to perform immune infiltration analysis of 22 immune genes in the samples. Furthermore, an analysis of the immune correlations of 7 crucial genes was conducted. Subsequently, a receiver operating characteristic curve was constructed to assess the diagnostic efficacy of these essential genes. A total of 171 differentially expressed genes were identified that were primarily implicated in the signaling pathways associated with apoptosis. Seven crucial genes (CAMP, MMP8, HMOX1, CTNNB1, ELANE, AKT1, and CEACAM8) were effectively filtered. The predominant functions of these genes were enriched in specific granules. The pivotal genes displayed robust associations with activated dendritic cells. The developed risk model showed a remarkable level of precision, as evidenced by an area under the curve of 0.8407 and C-index of 0.854. The present study successfully identified 7 crucial genes that are significantly associated with low BMD in perimenopausal women. Consequently, this research offers a solid theoretical foundation for clinical risk prediction, drug sensitivity analysis, and the development of targeted drugs specifically tailored for addressing low BMD in perimenopausal women.

Potential association of rheumatic diseases with bone mineral density and fractures: a bi-directional mendelian randomization study.

BACKGROUND: Previous studies have implicated rheumatoid arthritis as an independent risk factor for bone density loss. However, whether there is a causal relationship between rheumatic diseases and bone mineral density (BMD) and fractures is still controversial. We employed a bidirectional Mendelian analysis to explore the causal relationship between rheumatic diseases and BMD or fractures. METHODS: The rheumatic diseases instrumental variables (IVs) were obtained from a large Genome-wide association study (GWAS) meta-analysis dataset of European descent. Analyses were performed for the three rheumatic diseases: ankylosing spondylitis (AS) (n = 22,647 cases, 99,962 single nucleotide polymorphisms [SNPs]), rheumatoid arthritis (RA) (n = 58,284 cases, 13,108,512 SNPs), and systemic lupus erythematosus (SLE) (n = 14,267 cases, 7,071,163 SNPs). Two-sample Mendelian randomization (MR) analyses were carried out by using R language TwoSampleMR version 0.5.7. The inverse-variance weighted (IVW), MR-Egger, and weighted median methods were used to analyze the causal relationship between rheumatic diseases and BMD or fracture. RESULTS: The MR results revealed that there was absence of evidence for causal effect of AS on BMD or fracture. However, there is a positive causal relationship of RA with fracture of femur (95% CI = 1.0001 to 1.077, p = 0.046), and RA and fracture of forearm (95% CI = 1.015 to 1.064, p = 0.001). SLE had positive causal links for fracture of forearm (95% CI = 1.004 to 1.051, p = 0.020). Additionally, increasing in heel bone mineral density (Heel-BMD) and total bone mineral density (Total-BMD) can lead to a reduced risk of AS without heterogeneity or pleiotropic effects. The results were stable and reliable. There was absence of evidence for causal effect of fracture on RA (95% CI = 0.929 to 1.106, p = 0.759), and fracture on SLE (95% CI = 0.793 to 1.589, p = 0.516). CONCLUSIONS: RA and SLE are risk factors for fractures. On the other hand, BMD increasing can reduce risk of AS. Our results indicate that rheumatic diseases may lead to an increased risk of fractures, while increased BMD may lead to a reduced risk of rheumatic diseases. These findings provide insight into the risk of BMD and AS, identifying a potential predictor of AS risk as a reduction in BMD.

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.

Farnesyl Diphosphate Synthase Gene Associated with Loss of Bone Mass Density and Alendronate Treatment Failure in Patients with Primary Osteoporosis.

Aminobisphosphonates (NBPs) are the first-choice medication for osteoporosis (OP); NBP treatment aims at increasing bone mineral density (BMD) by inhibiting the activity of farnesyl diphosphate synthase (FDPS) enzyme in osteoclasts. Despite its efficacy, inadequate response to the drug and side effects have been reported. The A allele of the rs2297480 (A > C) SNP, found in the regulatory region of the FDPS gene, is associated with reduced gene transcription. This study evaluates the FDPS variant rs2297480 (A > C) association with OP patients' response to alendronate sodium treatment. A total of 304 OP patients and 112 controls were enrolled; patients treated with alendronate sodium for two years were classified, according to BMD variations at specific regions (lumbar spine (L1-L4), femoral neck (FN) and total hip (TH), as responders (OP-R) (n = 20) and non-responders (OP-NR) (n = 40). We observed an association of CC genotype with treatment failure (p = 0.045), followed by a BMD decrease in the regions L1-L4 (CC = -2.21% ± 2.56; p = 0.026) and TH (CC = -2.06% ± 1.84; p = 0.015) after two years of alendronate sodium treatment. Relative expression of the FDPS gene was also evaluated in OP-R and OP-NR patients. Higher expression of the FDPS gene was also observed in OP-NR group (FC = 1.84 ± 0.77; p = 0.006) when compared to OP-R. In conclusion, the influence observed of FDPS expression and the rs2897480 variant on alendronate treatment highlights the importance of a genetic approach to improve the efficacy of treatment for primary osteoporosis.

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.

Prevalence and effects of Vitamin D receptor polymorphism on bone mineral density and metabolism in patients with systemic sclerosis: a preliminary study.

Patients with systemic sclerosis (SSc) have a disproportionately high prevalence of reduced bone mineral density (BMD). Polymorphisms of the vitamin D receptor (VDR) gene have been associated with osteoporosis in patients with autoimmune diseases. The aim of this study was to investigate the prevalence and possible effects of VDR polymorphism on BMD and bone metabolism in patients with SSc. In patients with SSc measurement of BMD was performed using dual-energy X-ray absorptiometry. VDR polymorphisms (FokI, BsmI) were genotyped using restriction fragment length polymorphism analysis. Markers of bone metabolism (calcium, osteocalcin, ß-crosslaps) were determined. Primary endpoint was the prevalence of VDR gene polymorphisms and the association with reduced BMD. Secondary endpoints included associations between bone metabolism and VDR gene polymorphism. 79 Caucasian patients with SSc were included. Overall, 83.5% had reduced BMD (51.9% osteopenia, 31.6% osteoporosis). The prevalence of VDR gene polymorphism (73% BsmI, 77% FokI) was comparable to studies in healthy and rheumatic populations. The homozygous presence of FokI polymorphism, but not BsmI, was significantly associated with reduced axial BMD. Fokl polymorphism was significantly associated with reduced CTX levels, although changes remained within the reference limits. VDR polymorphisms can frequently be found in patients with SSc in comparable prevalence to healthy and rheumatic populations. The homozygous presence of FokI polymorphism, but not BsmI, was significantly associated with reduced axial BMD. This could be a possible contributor for the high prevalence of reduced BMD in 83.5% of patients with SSc in this study.Trial registration. DRKS00032768, date: 05.10.2023, retrospectively registered.

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.

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.

Apolipoprotein ɛ4 Is Associated With Increased Risk of Fall- and Fracture-Related Hospitalization: The Perth Longitudinal Study of Ageing Women.

Apolipoprotein ɛ4 (APOE ɛ4) may be a genetic risk factor for reduced bone mineral density (BMD) and muscle function, which could have implications for fall and fracture risk. We examined the association between APOE ɛ4 status and long-term fall- and fracture-related hospitalization risk in older women. A total of 1 276 community-dwelling women from the Perth Longitudinal Study of Aging Women (mean age ±â€…SD = 75.2 ±â€…2.7 years) were included. At baseline, women underwent APOE genotyping and detailed phenotyping for covariates including prevalent falls and fractures, as well as health and lifestyle factors. The association between APOE ɛ4 and fall-, any fracture-, and hip fracture-related hospitalizations, obtained over 14.5 years from linked health records, was examined using multivariable-adjusted Cox-proportional hazard models. Over 14.5 years, 507 (39.7%) women experienced a fall-related hospitalization and 360 (28.2%) women experienced a fracture-related hospitalization, including 143 (11.2%) attributed to a hip fracture. In multivariable-adjusted models, compared to noncarriers, APOE ɛ4 carriers (n = 297, 23.3%) had greater risk for a fall- (hazard ratio [HR] 1.48, 95% CI: 1.22-1.81), fracture- (HR 1.28, 95% CI: 1.01-1.63), or hip fracture-related hospitalization (HR 1.83, 95% CI: 1.29-2.61). The estimates remained similar when specific fall and fracture risk factors (fear of falling, plasma 25-hydroxyvitamin D, grip strength, timed up-and-go, hip BMD, vitamin K status, prevalent diabetes, HbA1c, cholesterol, and abbreviated mental test score) were added to the multivariable model. In conclusion, APOE ɛ4 is a potential risk factor for fall- and fracture-related hospitalization in community-dwelling older women. Screening for APOE ɛ4 could provide clinicians an opportunity to direct higher-risk individuals to appropriate intervention strategies.

Outdoor walking, genetic predisposition, and the risk of incident osteoporosis among older adults: A prospective large population-based cohort study.

This large-scale prospective study showed that a significant association between longer duration of daily outdoor walking and reduced osteoporosis risk was found among older adults, particularly among those with a low genetic predisposition to osteoporosis, which highlighted the importance of outdoor walking as a simple, cost-effective adjunct for preventing osteoporosis. PURPOSE: The available cross-sectional data and small-scale studies indicate that outdoor walking benefits bone metabolism. Nevertheless, there is a scarcity of comprehensive prospective research investigating the enduring correlation between outdoor walking and osteoporosis. This study aims to conduct a prospective analysis of the correlation between outdoor walking and osteoporosis while also examining potential variations influenced by genetic susceptibility to osteoporosis. METHODS: 24,700 older adults without osteoporosis at baseline were enrolled. These individuals were followed up until December 31, 2021, during which data on outdoor walking was gathered. The genetic risk score for osteoporosis was comprised of 14 single-nucleotide polymorphisms. RESULTS: 4,586 cases of osteoporosis were identified throughout a median follow-up period of 37.3 months. Those who walked outside for > 30 but ≤ 60 min per day had a hazard ratio (HR) of 0.83 (95% confidence interval (CI): 0.72-0.95) for incident osteoporosis, whereas those who walked outside for > 60 min per day had an HR of 0.60 (95% CI: 0.39-0.92). We found that osteoporosis risk exhibited a declining trend in individuals with low genetic risk. Individuals walking outside for > 60 min per day tended to have the lowest overall osteoporosis risk among those with high genetic risk. CONCLUSIONS: A significant negative correlation exists between an extended period of daily outdoor walking and osteoporosis incidence risk. This correlation is particularly pronounced among individuals with low genetic risk. The results above underscore the significance of outdoor walking as a simple and economical adjunct to public health programs to prevent osteoporosis.

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.

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.

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.

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.

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.

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.

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.

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.