Investigating the causal relationship between thyroid dysfunction diseases and osteoporosis: a two-sample Mendelian randomization analysis.

The prevalence of thyroid dysfunction diseases (TDFDs) and osteoporosis (OP) is high. Previous studies have indicated a potential association between TDFDs and OP, yet the causal direction remains unclear. This study aimed to investigate the potential causal relationship between TDFDs and the risk of developing OP and related fractures. We obtained pooled data from genome-wide association studies (GWASs) conducted on TDFDs and OP in European populations and identified single-nucleotide polymorphisms (SNPs) with genome-wide significance levels associated with exposure to TDFDs as instrumental variables. Inverse variance weighted (IVW) was employed as the primary method for Mendelian randomization (MR) analysis, supplemented by MR‒Egger, weighted median, simple mode and weighted mode methods. Sensitivity analyses were conducted to evaluate the robustness of the findings. The IVW method demonstrated an increased risk of OP in patients with TDFDs, including hyperthyroidism and hypothyroidism (TDFDs: OR = 1.11; 95% CI 1.09, 1.13; hypothyroidism: OR = 1.14; 95% CI 1.10, 1.17; hyperthyroidism: OR = 1.09; 95% CI 1.06, 1.12). These findings were supported by supplementary analysis, which revealed a positive correlation between TDFDs and the risk of OP. Multiple sensitivity analyses confirmed the absence of horizontal pleiotropy in the study, thus indicating the robustness of our results. The causal relationship between TDFDs and increased risk of OP implies the need for early bone mineral density (BMD) screening and proactive prevention and treatment strategies for individuals with TDFDs.

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.

GREM2 inactivation increases trabecular bone mass in mice.

Osteoporosis is a common skeletal disease affecting millions of individuals world-wide, with an increased risk of fracture, and a decreased quality of life. Despite its well-known consequences, the etiology of osteoporosis and optimal treatment methods are not fully understood. Human genetic studies have identified genetic variants within the FMN2/GREM2 locus to be associated with trabecular volumetric bone mineral density (vBMD) and vertebral and forearm fractures, but not with cortical bone parameters. GREM2 is a bone morphogenetic protein (BMP) antagonist. In this study, we employed Grem2-deficient mice to investigate whether GREM2 serves as the plausible causal gene for the fracture signal at the FMN2/GREM2 locus. We observed that Grem2 is moderately expressed in bone tissue and particularly in osteoblasts. Complete Grem2 gene deletion impacted mouse survival and body growth. Partial Grem2 inactivation in Grem2+/- female mice led to increased trabecular BMD of femur and increased trabecular bone mass in tibia due to increased trabecular thickness, with an unchanged cortical thickness, as compared with wildtype littermates. Furthermore, Grem2 inactivation stimulated osteoblast differentiation, as evidenced by higher alkaline phosphatase (Alp), osteocalcin (Bglap), and osterix (Sp7) mRNA expression after BMP-2 stimulation in calvarial osteoblasts and osteoblasts from the long bones of Grem2-/- mice compared to wildtype littermates. These findings suggest that GREM2 is a possible target for novel osteoporotic treatments, to increase trabecular bone mass and prevent osteoporotic fractures.

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

Background: Osteoporosis (OP) and cardiovascular disease (CVD) are major global public health issues, especially exacerbated by the challenges of an aging population. As these problems intensify, the associated burden on global health is expected to increase significantly. Despite extensive epidemiological investigations into the potential association between OP and CVD, establishing a clear causal relationship remains elusive. Methods: Instrumental variables were selected from summary statistics of the IEU GWAS database. Five different components of BMD (heel BMD, LS BMD, FA BMD, FN BMD, and TB BMD) were used as OP phenotypes. CHD, MI, and stroke were selected to represent CVD. Multiple analysis methods were used to evaluate the causal relationship between CVD and OP comprehensively. In addition, sensitivity analyses(Cochran's Q test, MR-Egger intercept test, and "leave one out" analysis) were performed to verify the reliability of the results. Results: The MR showed a significant causal relationship between CHD on heel BMD and TB BMD; in the reverse analysis, there was no evidence that OP has a significant causal effect on CVD. The reliability of the results was confirmed through sensitivity analysis. Conclusion: The study results revealed that CHD was causally associated with Heel BMD and TB BMD, while in the reverse MR analysis, the causal relationship between OP and CVD was not supported. This result posits CHD as a potential etiological factor for OP and prompts that routine bone density assessment at traditional sites (forearm, femoral neck, lumbar spine) using DAX may inadequately discern underlying osteoporosis issues in CHD patients. The recommendation is to synergistically incorporate heel ultrasound or DAX for total body bone density examinations, ensuring clinical diagnostics are both precise and reliable. Moreover, these findings provide valuable insights for public health, contributing to the development of pertinent prevention and treatment strategies.

Identification of shared gene signatures and pathways for diagnosing osteoporosis with sarcopenia through integrated bioinformatics analysis and machine learning.

BACKGROUND: Prior studies have suggested a potential relationship between osteoporosis and sarcopenia, both of which can present symptoms of compromised mobility. Additionally, fractures among the elderly are often considered a common outcome of both conditions. There is a strong correlation between fractures in the elderly population, decreased muscle mass, weakened muscle strength, heightened risk of falls, and diminished bone density. This study aimed to pinpoint crucial diagnostic candidate genes for osteoporosis patients with concomitant sarcopenia. METHODS: Two osteoporosis datasets and one sarcopenia dataset were obtained from the Gene Expression Omnibus (GEO). Differential expression genes (DEGs) and module genes were identified using Limma and Weighted Gene Co-expression Network Analysis (WGCNA), followed by functional enrichment analysis, construction of protein-protein interaction (PPI) networks, and application of a machine learning algorithm (least absolute shrinkage and selection operator (LASSO) regression) to determine candidate hub genes for diagnosing osteoporosis combined with sarcopenia. Receiver operating characteristic (ROC) curves and column line plots were generated. RESULTS: The merged osteoporosis dataset comprised 2067 DEGs, with 424 module genes filtered in sarcopenia. The intersection of DEGs between osteoporosis and sarcopenia module genes consisted of 60 genes, primarily enriched in viral infection. Through construction of the PPI network, 30 node genes were filtered, and after machine learning, 7 candidate hub genes were selected for column line plot construction and diagnostic value assessment. Both the column line plots and all 7 candidate hub genes exhibited high diagnostic value (area under the curve ranging from 1.00 to 0.93). CONCLUSION: We identified 7 candidate hub genes (PDP1, ALS2CL, VLDLR, PLEKHA6, PPP1CB, MOSPD2, METTL9) and constructed column line plots for osteoporosis combined with sarcopenia. This study provides reference for potential peripheral blood diagnostic candidate genes for sarcopenia in osteoporosis patients.

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.

MicroRNAs regulate the vicious cycle of vascular calcification-osteoporosis in postmenopausal women.

Menopause is a normal physiological process accompanied by changes in various physiological states. The incidence of vascular calcification (VC) increases each year after menopause and is closely related to osteoporosis (OP). Although many studies have investigated the links between VC and OP, the interaction mechanism of the two under conditions of estrogen loss remains unclear. MicroRNAs (miRNAs), which are involved in epigenetic modification, play a critical role in estrogen-mediated mineralization. In the past several decades, miRNAs have been identified as biomarkers or therapeutic targets in diseases. Thus, we hypothesize that these small molecules can provide new diagnostic and therapeutic approaches. In this review, we summarize the close interactions between VC and OP and the role of miRNAs in their interplay.

Causal roles of educational duration in bone mineral density and risk factors for osteoporosis: a Mendelian randomization study.

BACKGROUND: Educational duration might play a vital role in preventing the occurrence and development of osteoporosis(OP). PURPOSE: To assess the causal effect of educational duration on bone mineral density(BMD) and risk factors for OP by Mendelian randomization(MR) study. METHODS: The causal relationship was analyzed using data from genome-wide association study(GWAS). Inverse variance weighting (IVW) was used as the main analysis method. Horizontal pleiotropy was identified by MR-Egger intercept test, MR pleiotropy residual sum and outlier (MR-PRESSO) test. The leave-one-out method was used as a sensitivity analysis. RESULTS: The IVW results indicated that there was a positive causal relationship between educational duration and BMD (OR = 1.012, 95%CI:1.003-1.022), physical activity(PA) (OR = 1.156, 95%CI:1.032-1.295), calcium consumption (OR = 1.004, 95%CI:1.002-1.005), and coffee intake (OR = 1.019, 95%CI:1.014-1.024). There was a negative association between whole body fat mass (OR = 0.950, 95%CI:0.939-0.961), time for vigorous PA (OR = 0.955, 95%CI:0.939-0.972), sunbath (OR = 0.987, 95%CI:0.986-0.989), salt consumption (OR = 0.965, 95%CI:0.959-0.971), fizzy drink intake (OR = 0.985, 95%CI:0.978-0.992), smoking (OR = 0.969, 95%CI:0.964-0.975), and falling risk (OR = 0.976, 95%CI:0.965-0.987). There was no significant association between educational duration and lean mass, time for light-to-moderate PA, milk intake, and alcohol intake. Horizontal pleiotropy was absent in this study. The results were robust under sensitivity analyses. CONCLUSION: A longer educational duration was causally linked with increased BMD. No causal relationship had been found between educational duration and lean mass, time for light-to-moderate PA, milk intake, and alcohol consumption as risk factors for osteoporosis.

A Causal Relationship Between Type 1 Diabetes and Risk of Osteoporosis: A Univariable and Multivariable Mendelian Randomization Study.

Objective: This Mendelian randomization (MR) analysis aims to investigate the causal relationship between type 1 diabetes (T1D) and osteoporosis (OP). Methods: Single nucleotide polymorphisms (SNPs) associated with T1D were selected from the summary statistics of the genome-wide association study (GWAS) in European ancestry as instrumental variables (IVs) for univariable MR (UVMR) to explore the causal relationship between T1D and OP. Inverse variance weighting (IVW) was the primary method used to assess possible causality between T1D and OP. MR-PRESSO and MR-Egger intercepts were used to assess the horizontal pleiotropy of the IVs, and Q tests and the "leave-one-out" method were used to test for heterogeneity of MR results. Multivariable MR (MVMR) analysis was used to account for potential confounders such as smoking, obesity, drinking, and serum 25-hydroxyvitamin D (25OHD) concentrations. Result: Inverse variance weighted estimates suggest T1D may increase risk of OP (UVMR: OR = 1.06, 95% CI: 1.02-1.10, p = 0.002) (MVMR: OR = 1.50, 95% CI: 1.07-1.90, p < 0.001). Conclusion: Our findings suggest that T1D can increase the risk of OP.

Effects of circulating inflammatory proteins on osteoporosis and fractures: evidence from genetic correlation and Mendelian randomization study.

Objective: There is a controversy in studies of circulating inflammatory proteins (CIPs) in association with osteoporosis (OP) and fractures, and it is unclear if these two conditions are causally related. This study used MR analyses to investigate the causal associations between 91 CIPs and OP and 9 types of fractures. Methods: Genetic variants data for CIPs, OP, and fractures were obtained from the publicly available genome-wide association studies (GWAS) database. We used inverse variance weighted (IVW) as the primary analysis, pleiotropy, and heterogeneity tests to analyze the validity and robustness of causality and reverse MR analysis to test for reverse causality. Results: The IVW results with Bonferroni correction indicated that CXCL11 (OR = 1.2049; 95% CI: 1.0308-1.4083; P = 0.0192) can increase the risk of OP; IL-4 (OR = 1.2877; 95% CI: 1.1003-1.5070; P = 0.0016), IL-7 (OR = 1.2572; 95% CI: 1.0401-1.5196; P = 0.0180), IL-15RA (OR = 1.1346; 95% CI: 1.0163-1.2668; P = 0.0246), IL-17C (OR = 1.1353; 95% CI: 1.0272-1.2547; P = 0.0129), CXCL10 (OR = 1.2479; 95% CI: 1.0832-1.4377; P = 0.0022), eotaxin/CCL11 (OR = 1.1552; 95% CI: 1.0525-1.2678; P = 0.0024), and FGF23 (OR = 1.9437; 95% CI: 1.1875-3.1816; P = 0.0082) can increase the risk of fractures; whereas IL-10RB (OR = 0.9006; 95% CI: 0.8335-0.9730; P = 0.0080), CCL4 (OR = 0.9101; 95% CI: 0.8385-0.9878; P = 0.0242), MCP-3/CCL7 (OR = 0.8579; 95% CI: 0.7506-0.9806; P = 0.0246), IFN-γ [shoulder and upper arm (OR = 0.7832; 95% CI: 0.6605-0.9287; P = 0.0049); rib(s), sternum and thoracic spine (OR = 0.7228; 95% CI: 0.5681-0.9197; P = 0.0083)], ß-NGF (OR = 0.8384; 95% CI: 0.7473-0.9407; P = 0.0027), and SIRT2 (OR = 0.5167; 95% CI: 0.3296-0.8100; P = 0.0040) can decrease fractures risk. Conclusion: Mendelian randomization (MR) analyses indicated the causal associations between multiple genetically predicted CIPs and the risk of OP and fractures.

Exploration of shared gene signatures and molecular mechanisms between type 2 diabetes and osteoporosis.

Type 2 diabetes mellitus (T2D) and osteoporosis (OP) are systemic metabolic diseases and often coexist. The mechanism underlying this interrelationship remains unclear. We downloaded microarray data for T2D and OP from the Gene Expression Omnibus (GEO) database. Using weighted gene co-expression network analysis (WGCNA), we identified co-expression modules linked to both T2D and OP. To further investigate the functional implications of these associated genes, we evaluated enrichment using ClueGO software. Additionally, we performed a biological process analysis of the genes unique in T2D and OP. We constructed a comprehensive miRNA-mRNA network by incorporating target genes and overlapping genes from the shared pool. Through the implementation of WGCNA, we successfully identified four modules that propose a plausible model that elucidates the disease pathway based on the associated and distinct gene profiles of T2D and OP. The miRNA-mRNA network analysis revealed co-expression of PDIA6 and SLC16A1; their expression was upregulated in patients with T2D and islet ß-cell lines. Remarkably, PDIA6 and SLC16A1 were observed to inhibit the proliferation of pancreatic ß cells and promote apoptosis in vitro, while downregulation of PDIA6 and SLC16A1 expression led to enhanced insulin secretion. This is the first study to reveal the significant roles of PDIA6 and SLC16A1 in the pathogenesis of T2D and OP, thereby identifying additional genes that hold potential as indicators or targets for therapy.

Integration of bioinformatics and machine learning approaches for the validation of pyrimidine metabolism-related genes and their implications in immunotherapy for osteoporosis.

BACKGROUND: Osteoporosis (OP), the "silent epidemic" of our century, poses a significant challenge to public health, predominantly affecting postmenopausal women and the elderly. It evolves from mild symptoms to pronounced severity, stabilizing eventually. Unique among OP's characteristics is the altered metabolic profile of affected cells, particularly in pyrimidine metabolism (PyM), a crucial pathway for nucleotide turnover and pyrimidine decomposition. While metabolic adaptation is acknowledged as a therapeutic target in various diseases, the specific role of PyM genes (PyMGs) in OP's molecular response remains to be clarified. METHODS: In pursuit of elucidating and authenticating PyMGs relevant to OP, we embarked on a comprehensive bioinformatics exploration. This entailed the integration of Weighted Gene Co-expression Network Analysis (WGCNA) with a curated list of 37 candidate PyMGs, followed by the examination of their biological functions and pathways via Gene Set Variation Analysis (GSVA). The Least Absolute Shrinkage and Selection Operator (LASSO) technique was harnessed to identify crucial hub genes. We evaluated the diagnostic prowess of five PyMGs in OP detection and explored their correlation with OP's clinical traits, further validating their expression profiles through independent datasets (GSE2208, GSE7158, GSE56815, and GSE35956). RESULTS: Our analytical rigor unveiled five PyMGs-IGKC, TMEM187, RPS11, IGLL3P, and GOLGA8N-with significant ties to OP. A deeper dive into their biological functions highlighted their roles in estrogen response modulation, cytosolic calcium ion concentration regulation, and GABAergic synaptic transmission. Remarkably, these PyMGs emerged as potent diagnostic biomarkers for OP, distinguishing affected individuals with substantial accuracy. CONCLUSIONS: This investigation brings to light five PyMGs intricately associated with OP, heralding new avenues for biomarker discovery and providing insights into its pathophysiological underpinnings. These findings not only deepen our comprehension of OP's complexity but also herald the advent of more refined diagnostic and therapeutic modalities.

CircCOX6A1 suppresses osteogenic differentiation and aggravates osteoporosis via miR-512-3p/DYRK2 axis.

BACKGROUND: Osteoporosis (OP), characterized by compromised bone integrity and increased fracture risk, poses a significant health challenge. Circular RNAs (circRNAs) have emerged as crucial regulators in various pathophysiological processes, prompting investigation into their role in osteoporosis. This study aimed to elucidate the involvement of circCOX6A1 in OP progression and understand its underlying molecular mechanisms. The primary objective was to explore the impact of circCOX6A1 on bone marrow-derived mesenchymal stem cells (BMSCs) and its potential interactions with miR-512-3p and DYRK2. METHODS: GSE161361 microarray analysis was employed to assess circCOX6A1 expression in OP patients. We utilized in vitro and in vivo models, including BMSC cultures, osteogenic differentiation assays, and an OVX-induced mouse model of OP. Molecular techniques such as quantitative RT-PCR, western blotting, and functional assays like alizarin red staining (ARS) were employed to evaluate circCOX6A1 effects on BMSC proliferation, apoptosis, and osteogenic differentiation. The interaction between circCOX6A1, miR-512-3p, and DYRK2 was investigated through dual luciferase reporter assays, RNA immunoprecipitation, and RNA pull-down assays. RESULTS: CircCOX6A1 was found to be upregulated in osteoporosis patients, and its expression inversely correlated with osteogenic differentiation of BMSCs. CircCOX6A1 knockdown enhanced osteogenic differentiation, as evidenced by increased mineralized nodule formation and upregulation of osteogenic markers. In vivo, circCOX6A1 knockdown ameliorated osteoporosis progression in OVX mice. Mechanistically, circCOX6A1 acted as a sponge for miR-512-3p, subsequently regulating DYRK2 expression. CONCLUSION: This study provides compelling evidence for the role of circCOX6A1 in osteoporosis pathogenesis. CircCOX6A1 negatively regulates BMSC osteogenic differentiation through the miR-512-3p/DYRK2 axis, suggesting its potential as a therapeutic target for mitigating OP progression.

Role of epigenetic regulatory mechanisms in age-related bone homeostasis imbalance.

Alterations to the human organism that are brought about by aging are comprehensive and detrimental. Of these, an imbalance in bone homeostasis is a major outward manifestation of aging. In older adults, the decreased osteogenic activity of bone marrow mesenchymal stem cells and the inhibition of bone marrow mesenchymal stem cell differentiation lead to decreased bone mass, increased risk of fracture, and impaired bone injury healing. In the past decades, numerous studies have reported the epigenetic alterations that occur during aging, such as decreased core histones, altered DNA methylation patterns, and abnormalities in noncoding RNAs, which ultimately lead to genomic abnormalities and affect the expression of downstream signaling osteoporosis treatment and promoter of fracture healing in older adults. The current review summarizes the impact of epigenetic regulation mechanisms on age-related bone homeostasis imbalance.

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.

RNA sequencing-based approaches to identifying disulfidptosis-related diagnostic clusters and immune landscapes in osteoporosis.

Disulfidptosis, a newly recognized cell death triggered by disulfide stress, has garnered attention for its potential role in osteoporosis (OP) pathogenesis. Although sulfide-related proteins are reported to regulate the balance of bone metabolism in OP, the precise involvement of disulfidptosis regulators remains elusive. Herein, leveraging the GSE56815 dataset, we conducted an analysis to delineate disulfidptosis-associated diagnostic clusters and immune landscapes in OP. Subsequently, vertebral bone tissues obtained from OP patients and controls were subjected to RNA sequencing (RNA-seq) for the validation of key disulfidptosis gene expression. Our analysis unveiled seven significant disulfidptosis regulators, including FLNA, ACTB, PRDX1, SLC7A11, NUBPL, OXSM, and RAC1, distinguishing OP samples from controls. Furthermore, employing a random forest model, we identified four diagnostic disulfidptosis regulators including FLNA, SLC7A11, NUBPL, and RAC1 potentially predictive of OP risk. A nomogram model integrating these four regulators was constructed and validated using the GSE35956 dataset, demonstrating promising utility in clinical decision-making, as affirmed by decision curve analysis. Subsequent consensus clustering analysis stratified OP samples into two different disulfidptosis subgroups (clusters A and B) using significant disulfidptosis regulators, with cluster B exhibiting higher disulfidptosis scores and implicating monocyte immunity, closely linked to osteoclastogenesis. Notably, RNA-seq analysis corroborated the expression patterns of two disulfidptosis modulators, PRDX1 and OXSM, consistent with bioinformatics predictions. Collectively, our study sheds light on disulfidptosis patterns, offering potential markers and immunotherapeutic avenues for future OP management.

Carboxypeptidase M modulates BMSCs osteogenesis-adipogenesis via the MAPK/ERK pathway: An integrated single-cell and bulk transcriptomic study.

The pathogenesis of osteoporosis (OP) is closely associated with the disrupted balance between osteogenesis and adipogenesis in bone marrow-derived mesenchymal stem cells (BMSCs). We analyzed published single-cell RNA sequencing (scRNA-seq) data to dissect the transcriptomic profiles of bone marrow-derived cells in OP, reviewing 56 377 cells across eight scRNA-seq datasets from femoral heads (osteoporosis or osteopenia n = 5, osteoarthritis n = 3). Seventeen genes, including carboxypeptidase M (CPM), were identified as key osteogenesis-adipogenesis regulators through comprehensive gene set enrichment, differential expression, regulon activity, and pseudotime analyses. In vitro, CPM knockdown reduced osteogenesis and promoted adipogenesis in BMSCs, while adenovirus-mediated CPM overexpression had the reverse effects. In vivo, intraosseous injection of CPM-overexpressing BMSCs mitigated bone loss in ovariectomized mice. Integrated scRNA-seq and bulk RNA sequencing analyses provided insight into the MAPK/ERK pathway's role in the CPM-mediated regulation of BMSC osteogenesis and adipogenesis; specifically, CPM overexpression enhanced MAPK/ERK signaling and osteogenesis. In contrast, the ERK1/2 inhibitor binimetinib negated the effects of CPM overexpression. Overall, our findings identify CPM as a pivotal regulator of BMSC differentiation, which provides new clues for the mechanistic study of 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.

The hidden dangers of short-term glucocorticoid use in children: A genomic analysis.

OBJECTIVE: Glucocorticoid (GC) administration has been associated with adverse drug reactions (ADRs) affecting multiple organ systems. While long-term use is widely recognized as a significant independent predictor of ADRs, it is important to note that even short-term use can lead to serious ADRs. The considerable inter-individual variability in ADRs occurrence may be influenced by genetic factors. This study, we present a case of a child who experienced significant weight gain and osteoporosis, following a brief administration of GC. METHODS: To comprehensively investigate the underlying mechanisms, we conducted a genomic analysis utilizing the whole exome sequencing (WES) technique. This analysis encompassed the examination of phase I and phase II metabolism, influx transport, efflux transport, and drug targeting. Additionally, a comprehensive analysis was conducted on a cohort of 52,119 children to determine their ABCB1 rs1045642 genotype, and an additional 37,884 children were tested for their CYP3A5 rs776746 genotype. RESULTS: The pharmacogenetic analysis unveiled the presence of a high-risk variant in ABCB1 rs1045642 and a slow metabolism variant in CYP3A5 rs776746, both of which have the potential to substantially contribute to ADRs. The findings of this study indicate that the prevalence of ABCB1 rs1045642 CT type among patients was 47.58%, with TT type accounting for 15.69 % and CC type accounting for 36.73 %. Furthermore, the distribution of CYP3A5 rs776746 CC genotype was observed in 50.54 % of individuals, while CT and TT genotypes were present in 41.15 % and 8.31 % of the population respectively. The distribution of ABCB1 and CYP3A5 genotypes among the pediatric population in China displays notable features. Specifically, for the ABCB1 rs1045642 genotype, less than 50 % of children exhibit intermediate metabotypes. Conversely, among children with the CYP3A5 rs776746 genotype, the predominant cause for enzyme activity is the slow metabolic type, accounting for up to 90 % of cases. CONCLUSIONS: Consequently, it is imperative to thoroughly evaluate the impact of allele mutation on the effectiveness and safety of glucocorticoid drugs or other medications metabolized by the ABCB1 and CYP3A5, particularly in the context of Chinese pediatric patients.