MiR-210 promotes bone formation in ovariectomized rats by regulating osteogenic/adipogenic differentiation of bone marrow mesenchymal stem cells through downregulation of EPHA2

PURPOSE

In osteoporosis, the balance between osteogenic and adipogenic differentiation of mesenchymal stem cells (MSCs) is disrupted. The osteogenic differentiation of bone marrow MSCs (BMSCs) is important for improving osteoporosis. The aim of this study was to explore the role and molecular mechanism of miR-210 in the balance of osteogenic/adipogenic differentiation of BMSCs in postmenopausal osteoporosis.

METHODS

Postmenopausal osteoporosis rat models were constructed by ovariectomy (OVX). BMSCs were isolated from the femur in rats of Sham and OVX groups. MiR-210 was overexpressed and suppressed by miR-210 mimics and inhibitor, respectively. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to detect the relative mRNA expression of miR-210, ephrin type-A receptor 2 (EPHA2), alkaline phosphatase (ALP), osterix (OSX), osteocalcin (Bglap), Runt-related transcription factor 2 (Runx2), peroxisome proliferator activated receptor gamma, and fatty acid binding protein 4 (FABP4) in each group of rat femoral tissues or BMSCs. Western blot was applied to detect the protein expression level of EPHA2 in rat femoral tissues and cells. Alizarin red S staining and oil red O staining were performed to assess the osteogenic and adipogenic differentiation of BMSCs, respectively. In addition, the targeting relationship between miR-210 and EPHA2 was verified by a dual luciferase gene reporter assay.

RESULTS

The expression of miR-210 was significantly reduced in femoral tissues and BMSCs of OVX rats, and its low expression was associated with reduced bone formation. The osteogenic differentiation was enhanced in OVX rats treated with miR-210 mimic. Overexpression of miR-210 in transfected BMSCs was also found to significantly promote osteogenic differentiation and even inhibit adipogenic differentiation in BMSCs, while knockdown of miR-210 did the opposite. Further mechanistic studies showed that miR-210 could target and inhibit the expression of EPHA2 in BMSCs, thus promoting osteogenic differentiation and inhibiting adipogenic differentiation of BMSCs.

CONCLUSION

MiR-210 promotes osteogenic differentiation and inhibits adipogenic differentiation of BMSCs by down-regulating EPHA2 expression. As it plays an important role in the osteogenic/adipogenic differentiation of osteoporosis, miR-210 can serve as a potential miRNA biomarker for osteoporosis.

Long non-coding RNA Homeobox D gene cluster antisense growth-associated long noncoding RNA/microRNA-182-5p/Homeobox protein A10 alleviates postmenopausal osteoporosis via accelerating osteoblast differentiation of bone marrow mesenchymal stem cells

BACKGROUND

Studies have illuminated that long non-coding RNA (lncRNA) influences bone cell differentiation and formation. Nevertheless, whether lncRNA Homeobox D gene cluster antisense growth-associated long noncoding RNA (HAGLR) was implicated in postmenopausal osteoporosis (PMOP) was yet uncertain.

PURPOSE

The research was to explore HAGLR's role in the osteogenic differentiation (OD) process of bone marrow mesenchymal stem cells (BMSCs).

METHODS

BMSCs were isolated from mouse bone marrow tissues and identified by electron microscope and flow cytometry. HAGLR, microRNA (miR)-182-5p, and homeobox protein A10 (Hoxa10) levels in BMSCs were detected. Mouse BMSC OD process was induced, and calcium deposition and alkaline phosphatase content were analyzed, as well as expressions of runt-related transcription factor 2, osteopontin, and osteocalcin, and cell apoptosis. Bilateral ovaries were resected from mice to construct the ovariectomized model and bone mineral density, maximum bending stress, maximum load, and elastic modulus of the femur were tested, and the femur was histopathologically evaluated. Chondrocyte apoptosis in the articular cartilage of mice was analyzed. Analysis of the interaction of HAGLR, miR-182-5p with Hoxa10 was conducted.

RESULTS

HAGLR and Hoxa10 were down-regulated and miR-182-5p was elevated in PMOP patients. During the BMSC OD process, HAGLR and Hoxa10 levels were suppressed, while miR-182-5p was elevated. Promotion of HAGLR or suppression of miR-182-5p accelerated OD of BMSCs. Inhibition of miR-182-5p reversed the inhibitory effect of HAGLR on BMSC OD. In in vivo experiments, up-regulating HAGLR alleviated PMOP, while silencing Hoxa10 reversed the effects of upregulating HAGLR. HAGLR performed as a sponge for miR-182-5p, while miR-182-5p targeted Hoxa10.

CONCLUSION

In general, HAGLR boosted the OD process of BMSCs and relieved PMOP via the miR-182-5p/Hoxa10 axis. These data preliminarily reveal the key role of HAGLR in PMOP, and the research results have a certain reference for the treatment of PMOP.

Identification of ROCK1 as a novel biomarker for postmenopausal osteoporosis and pan-cancer analysis

BACKGROUND

Postmenopausal osteoporosis (PMOP) is a prevalent bone disorder with significant global impact. The elevated risk of osteoporotic fracture in elderly women poses a substantial burden on individuals and society. Unfortunately, the current lack of dependable diagnostic markers and precise therapeutic targets for PMOP remains a major challenge.

METHODS

PMOP-related datasets GSE7429, GSE56814, GSE56815, and GSE147287, were downloaded from the GEO database. The DEGs were identified by "limma" packages. WGCNA and Machine Learning were used to choose key module genes highly related to PMOP. GSEA, DO, GO, and KEGG enrichment analysis was performed on all DEGs and the selected key hub genes. The PPI network was constructed through the GeneMANIA database. ROC curves and AUC values validated the diagnostic values of the hub genes in both training and validation datasets. xCell immune infiltration and single-cell analysis identified the hub genes' function on immune reaction in PMOP. Pan-cancer analysis revealed the role of the hub genes in cancers.

RESULTS

A total of 1278 DEGs were identified between PMOP patients and the healthy controls. The purple module and cyan module were selected as the key modules and 112 common genes were selected after combining the DEGs and module genes. Five Machine Learning algorithms screened three hub genes (KCNJ2, HIPK1, and ROCK1), and a PPI network was constructed for the hub genes. ROC curves validate the diagnostic values of ROCK1 in both the training (AUC = 0.73) and validation datasets of PMOP (AUC = 0.81). GSEA was performed for the low-ROCK1 patients, and the top enriched field included protein binding and immune reaction. DCs and NKT cells were highly expressed in PMOP. Pan-cancer analysis showed a correlation between low ROCK1 expression and SKCM as well as renal tumors (KIRP, KICH, and KIRC).

CONCLUSIONS

ROCK1 was significantly associated with the pathogenesis and immune infiltration of PMOP, and influenced cancer development, progression, and prognosis, which provided a potential therapy target for PMOP and tumors. However, further laboratory and clinical evidence is required before the clinical application of ROCK1 as a therapeutic target.

Hmga1-overexpressing lentivirus protects against osteoporosis by activating the Wnt/β-catenin pathway in the osteogenic differentiation of BMSCs

Postmenopausal osteoporosis is associated with bone formation inhibition mediated by the impaired osteogenic differentiation potential of bone marrow mesenchymal stem cells (BMSCs). However, identifying and confirming the essential genes in the osteogenic differentiation of BMSCs and osteoporosis remain challenging. The study aimed at revealing the key gene that regulated osteogenic differentiation of BMSCs and led to osteoporosis, thus exploring its therapeutic effect in osteoporosis. In the present study, six essential genes related to the osteogenic differentiation of BMSCs and osteoporosis were identified, namely, fibrillin 2 (Fbn2), leucine-rich repeat-containing 17 (Lrrc17), heat shock protein b7 (Hspb7), high mobility group AT-hook 1 (Hmga1), nexilin F-actin-binding protein (Nexn), and endothelial cell-specific molecule 1 (Esm1). Furthermore, the in vivo and in vitro experiments showed that Hmga1 expression was increased during the osteogenic differentiation of rat BMSCs, while Hmga1 expression was decreased in the bone tissue of ovariectomized (OVX) rats. Moreover, the expression of osteogenic differentiation-related genes, the activity of alkaline phosphatase (ALP), and the number of mineralized nodules were increased after Hmga1 overexpression, which was partially reversed by a Wnt signaling inhibitor (DKK1). In addition, after injecting Hmga1-overexpressing lentivirus into the bone marrow cavity of OVX rats, the bone loss, and osteogenic differentiation inhibition of BMSCs in OVX rats were partially reversed, while osteoclast differentiation promotion of BMSCs in OVX rats was unaffected. Taken together, the present study confirms that Hmga1 prevents OVX-induced bone loss by the Wnt signaling pathway and reveals that Hmga1 is a potential gene therapeutic target for postmenopausal osteoporosis.

Anti-osteoporosis mechanism of resistance exercise in ovariectomized rats based on transcriptome analysis: a pilot study

Postmenopausal osteoporosis is the main cause of fractures in women. Resistance exercise has a positive effect on bone mineral density in postmenopausal osteoporosis patients, but its mechanism is unclear. The purpose of this study was to explore the mechanism of resistance exercise in improving ovariectomized osteoporotic rats based on the transcriptome sequencing technique. Eighteen female Sprague-Dawley rats were randomly divided into the sham-operated group, the non-exercise group, and the resistance exercise group. The rat model of postmenopausal osteoporosis was established by bilateral ovariectomy. Ten weeks after the operation, the resistance exercise group received 2 weeks of adaptive training, and 12 weeks of resistance exercise began in the 13th week. The rats were trained 5 days per week, in 4 sets of 3 repetitions per day. After the intervention, all rats were sacrificed, and the body weight, bone mineral density, trabecular bone microarchitecture, and bone biomechanics were examined. At the same time, RNA-seq and enrichment analysis of gene ontology and Kyoto Encyclopedia of Genes and Genomes were performed on the left tibias, followed by Elisa and RT-qPCR verification. It had been found that resistance exercise can effectively counteract the weight gain of ovariectomized osteoporotic rats, and has a good effect on bone mineral density and trabecular bone microarchitecture. Enrichment analysis showed that regulation of gene expression and osteoclast differentiation is the most closely related biological process and signaling pathway shared by RE/Ovx and NE/Ovx groups. Our results revealed that resistance exercise can play a role in inhibiting osteoclast activation and preventing the enhancement of osteoclast bone resorption function in ovariectomized osteoporotic rats by inhibiting Fos/Fosb-regulated TRAP activation and relieving Calcr inhibition, which has important application value in preventing bone loss caused by estrogen deficiency.

Transcriptomic analysis the mechanisms of anti-osteoporosis of desert-living Cistanche herb in ovariectomized rats of postmenopausal osteoporosis

Desert-living Cistanche herb (DC), as a traditional Chinese medicine for tonifying kidney yang, is often used to treat postmenopausal osteoporosis (PMOP). Total phenylethanoid glycosides are instruction ingredients for discrimination and assay according to the China pharmacopoeia for DC. This research aimed to reveal the anti-osteoporosis mechanism of total phenylethanoid glycosides of DC (PGC) by transcriptomic analysis of ovariectomized rats. Serum levels of BGP were evaluated by ELISA, the bone weight was measured, and transmission electron microscopy was used to examine the ultrastructure of osteoblasts in rats. In addition, micro-CT was used to detect the bone volume (Tb.BS/BV), bone mineral density (Tb.BMD), and bone mineral content (Tb.BMC) in trabecular bone, and the ratio of cortical bone area to total area (Ct.ar/Tt.ar), and the level of bone mineral content (Ct.BMC) in cortical bone. Differential expressed genes (DEGs) after PGC treatment were analyzed by transcriptomics. Then, a bioinformatics analysis of DEGs was carried out through GO enrichment, KEGG enrichment, and selection of the nucleus gene through the protein-protein interaction network. Through qRT-PCR analysis, the DEGs were verified. The analysis results indicated that PGC increased the secretion of osteogenic markers, and ultrastructural characterization of osteoblasts and bone morphology were improved in ovariectomized rats. A total of 269 genes were differentially expressed, including 201 genes that were downregulated and 68 genes that were upregulated between the model group and the PGC group. Bioinformation analysis results prompt the conclusion that PGC could promote the bone metabolism by muscle cell development, myofibril assembly, etc. In addition, our study also found that PGC has a good effect on osteoporosis complicated with cardiomyopathy, and it also provided evidence for the correlation between sarcopenia and osteoporosis.

Identification of biomarkers associated with oxidative stress-related genes in postmenopausal osteoporosis

Osteoporosis (OP) is a prevalent metabolic disease, with aging and menopause being the major risk factors. Studies have shown that nearly one-third of postmenopausal women suffer from osteoporosis. However, there is a scarcity of research on antioxidant systems for the prevention and treatment of postmenopausal osteoporosis (PM-OP). To address this gap, we performed differential analysis using Limma to identify differentially expressed genes (DEGs) in PM-OP samples. We employed weighted correlation network analysis (WGCNA) to identify oxidative stress (OS)-related genes (OSRGs) highly correlated with PM-OP. The intersection of key modular genes and DEGs yielded differentially expressed OSRGs (DE-OSRGs) specific to PM-OP. We conducted GO and KEGG functional enrichment analyses on these genes. Additionally, we constructed a protein-protein interaction (PPI) network and utilized support vector machine recursive feature elimination (SVM-RFE) and random forest (RF) algorithms to identify signature genes. The diagnostic value of the signature genes was evaluated and validated using ROC curves. GSEA enrichment analysis was employed to explore the potential mechanisms associated with the signature genes. Finally, we constructed a regulatory network involving TF-miRNA-mRNA interactions for the signature genes and verified the biological roles of FOXO3 and DDIT3 in PM-OP and healthy groups using quantitative real-time PCR (qRT-PCR). Our analysis revealed 20 DE-OSRGs specific to PM-OP, obtained by intersecting modular and differential genes. The PPI network identified central genes (DDIT3, MAPK8, CDK2, SIRT1, and FOXO3) with more than 3 nodes. Through integration with machine learning algorithms, we identified DDIT3 and FOXO3 as signature genes. The ROC curve analysis indicated that the AUC value was greater than 0.7, suggesting the potential diagnostic value of these signature genes. Furthermore, GSEA results revealed their involvement in pathways related to the regulation of neutrophil activation, oxidative phosphorylation, MAPK signaling, mitochondrial matrix, and phagocytosis. Lastly, we constructed a regulatory network comprising 27 nodes (22 TFs, 3 miRNAs, and 2 mRNAs) and 28 edges. Additionally, qRT-PCR confirmed the significant up-regulation of FOXO3 and DDIT3 expressions in the PM-OP group compared to the healthy control group. In summary, this study employed bioinformatics analysis to identify OS-related biomarkers (DDIT3 and FOXO3) in PM-OP, providing new biological targets for clinical treatment.

Evaluation of COLIA1 gene rs1107946 polymorphism in relation to bone mineral density and fracture risk in Slovak postmenopausal women

OBJECTIVES

The aim of the study was the evaluation of the rs1107946 polymorphism of the COLIA1 gene impact on bone mineral density and fracture risk in Slovak postmenopausal women.

METHODS

One hundred and twenty-seven postmenopausal Slovak women with a diagnosis of osteopenia/osteoporosis were genotyped for rs1107946 polymorphism of the COLIA1 gene. Clinical and anthropometric data were obtained. DNA isolation was performed using a standard protocol. Genetic analyses of the rs1107946 polymorphism of the COLIA1 gene were performed by the TaqMan SNP genotyping assays.

RESULTS

The study confirmed a statistically significant relationship using an association analysis between the rs1107946 polymorphism of the COLIA1 gene genotypes and body weight of the Slovak postmenopausal women with osteopenia/osteoporosis (p = 0.03). The study revealed a significant association of the risk T allele of the rs1107946 polymorphism of the COLIA1 gene with osteoporotic fractures (p = 0.038). The odds ratio confirmed 2.060 times higher risk of osteoporotic fractures in Slovak postmenopausal women with the presence of risk T allele of the rs1107946 COLIA1 gene polymorphism (OR = 2.060; 95% CI: 1.024-4.144).

CONCLUSION

The results of this study revealed an association of T allele of the rs1107946 COLIA1 gene polymorphism with osteoporotic fractures in Slovak postmenopausal women with osteopenia/osteoporosis and suggest that the rs1107946 polymorphism of the COLIA1 gene may be a molecular biomarker usable in the management of osteoporosis.

Evaluation of Relationship Between SOD1 50-bp Deletion Gene Polymorphism, Cu, Zn Level, and Viscosity in Postmenopausal Osteoporosis Patients with Vertebral Fractures

Oxidative stress plays a role in the pathogenesis of bone loss, causing low bone mineral density (BMD) and associated osteoporotic fractures. In our study, we aimed to investigate the relationship of SOD1 50-bp insertion(Ins)/deletion(Del) polymorphism that is involved in oxidative stress metabolism, Cu and Zn element concentrations, and plasma viscosity level, with postmenopausal osteoporosis and related vertebral fractures. The study included 167 voluntary individuals. The 50-bp Ins/Del polymorphism of SOD1 was determined by allele-specific PCR. Plasma Cu and Zn levels were measured by atomic absorption spectrophotometry (AAS). The plasma viscosity was determined using the Harkness Capillary Viscometer device. In our study, the distribution of SOD1 promoter 50-bp Ins/Del polymorphism did not indicate a significant difference between the groups and in postmenopausal osteoporosis patients with and without fractures (p > 0.05). The Ins/Ins genotype was found to be common in individuals in both groups. The Cu and Zn levels of the study group were found to be between the normal reference values (p > 0.05). It was determined that plasma viscosity increased significantly in the group of osteoporotic patients and in patients with postmenopausal osteoporosis with fractures (p < 0.01). In addition, plasma viscosity was found to significantly increase in patients with Ins/Ins genotype and fractures (p < 0.01). Postmenopausal osteoporosis and associated vertebral fracture were found not to be directly related to SOD1 50-bp polymorphism and Cu and Zn element levels. Plasma viscosity levels were found to increase due to the increase in oxidative stress products. Further studies are needed to clarify the roles and relationships of SOD genes and trace elements in the development of postmenopausal osteoporosis and vertebral fracture.

Bioinformatics identification and experimental validation of m6A-related diagnostic biomarkers in the subtype classification of blood monocytes from postmenopausal osteoporosis patients

BACKGROUND

Postmenopausal osteoporosis (PMOP) is a common bone disorder. Existing study has confirmed the role of exosome in regulating RNA N6-methyladenosine (m6A) methylation as therapies in osteoporosis. However, it still stays unclear on the roles of m6A modulators derived from serum exosome in PMOP. A comprehensive evaluation on the roles of m6A modulators in the diagnostic biomarkers and subtype identification of PMOP on the basis of GSE56815 and GSE2208 datasets was carried out to investigate the molecular mechanisms of m6A modulators in PMOP.

METHODS

We carried out a series of bioinformatics analyses including difference analysis to identify significant m6A modulators, m6A model construction of random forest, support vector machine and nomogram, m6A subtype consensus clustering, GO and KEGG enrichment analysis of differentially expressed genes (DEGs) between different m6A patterns, principal component analysis, and single sample gene set enrichment analysis (ssGSEA) for evaluation of immune cell infiltration, experimental validation of significant m6A modulators by real-time quantitative polymerase chain reaction (RT-qPCR), etc.

RESULTS

In the current study, we authenticated 7 significant m6A modulators via difference analysis between normal and PMOP patients from GSE56815 and GSE2208 datasets. In order to predict the risk of PMOP, we adopted random forest model to identify 7 diagnostic m6A modulators, including FTO, FMR1, YTHDC2, HNRNPC, RBM15, RBM15B and WTAP. Then we selected the 7 diagnostic m6A modulators to construct a nomogram model, which could provide benefit with patients according to our subsequent decision curve analysis. We classified PMOP patients into 2 m6A subtypes (clusterA and clusterB) on the basis of the significant m6A modulators via a consensus clustering approach. In addition, principal component analysis was utilized to evaluate the m6A score of each sample for quantification of the m6A subgroups. The m6A scores of patients in clusterB were higher than those of patients in clusterA. Moreover, we observed that the patients in clusterA had close correlation with immature B cell and gamma delta T cell immunity while clusterB was linked to monocyte, neutrophil, CD56dim natural killer cell, and regulatory T cell immunity, which has close connection with osteoclast differentiation. Notably, m6A modulators detected by RT-qPCR showed generally consistent expression levels with the bioinformatics results.

CONCLUSION

In general, m6A modulators exert integral function in the pathological process of PMOP. Our study of m6A patterns may provide diagnostic biomarkers and immunotherapeutic strategies for future PMOP treatment.

Can Polymorphisms in NLRP3 Inflammasome Complex Be Associated with Postmenopausal Osteoporosis Severity?

The immune system plays a critical role in bone homeostasis and, consequently, in the pathophysiology of postmenopausal osteoporosis (OP) since estrogen deficiency induces the inflammasome and increases production of pro-inflammatory cytokines, such as IL-1β and IL-18. NLRP3 inflammasome complex genes have been related with bone homeostasis in cellular and animal models. Here, we performed an association study evaluating SNVs (single-nucleotide variants) in inflammasome NLRP3 pathway genes (NLRP3, CARD8, CASP1, IL-18, and IL-1β) to assess whether variants in these genes could be related to susceptibility to primary OP in postmenopausal women.

METHODS

We genotyped 196 postmenopausal OP patients and 103 healthy controls using SNV-specific Taqman^® probes. Data and statistical analyses were performed using the SNPstats and GraphPad Prism 8 software.

RESULTS

We showed an association between NLRP3 rs35829419 CA genotype and lower bone mineral density (BMD) mean at the lumbar spine (p = 0.001); we also observed an association between IL-1β rs16944 AA genotype and higher BMD mean at the total hip (p = 0.009). The IL-1β rs16944 GG was associated with lower alkaline phosphatase levels (ALP) (p = 0.009), and the IL-18 rs1946519 AA was associated with lower vitamin D levels (p = 0.018). Additionally, OP patients presented deficient vitamin D and parathyroid hormone (PTH).

CONCLUSIONS

The NLRP3 inflammasome complex SNVs were associated with OP severity, possibly indicating these genes' participation in bone metabolism and its dysregulation.

Network pharmacology-based strategy to investigate pharmacological mechanism of Liuwei Dihuang Pill against postmenopausal osteoporosis

Postmenopausal osteoporosis (PMOP) has became 1 of most prevalent bone disorders with aging population. Liuwei Dihuang (LWDH) Pill, a classical kidney-tonifying prescription, is extensively used to treat PMOP in China. The aim of this study is to explore the pharmacological mechanisms of LWDH Pill against PMOP via network pharmacological strategy. The active ingredients of LWDH Pill were screened out from the Traditional Chinese Medicine System Pharmacology, Encyclopedia of Traditional Chinese Medicine and Bioinformatics Analysis Tool for Molecular mechANism of Traditional Chinese Medicine Databases, and their related target genes were fished in the UniProt database. Simultaneously, the GeneCards and DisGeNET databases were used to identify the target genes of PMOP. Through establishing a protein-protein interaction network, the overlapping genes between LWDH Pill and PMOP were identified to analyze their interactions and the hub target genes. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed to predict the underlying biological processes (BP) and signaling pathways, respectively. A total of 64 active ingredients and 653 related target genes were identified in LWDH Pill, and 292 target genes were closely associated with PMOP. After matching the target genes between LWDH Pill and PMOP, 84 overlapping targets were obtained and considered as therapeutically relevant. Through construction of a protein-protein interaction network, we identified 20 hub target genes including IL6, INS, tumor necrosis factor, AKT1, vascular endothelial growth factor A, IGF1, TP53, IL1B, MMP9, JUN, LEP, CTNNB1, EGF, PTGS2, PPARG, CXCL8, IL10, CCL2, FOS and ESR1. Gene Ontology enrichment analysis suggested that LWDH Pill exerted anti-PMOP effects via regulating multiple BP including cell proliferation and apoptosis, oxidative stress, inflammation and angiogenesis. Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed several pathways, such as PI3K-AKT pathway, mitogen-activated protein kinase pathway, hypoxia-inducible factors-1 pathway, tumor necrosis factor pathway, interleukin-17 (IL-17) pathway and FoxO pathway that might be involved in modulating the above BP. Through network pharmacological approach, we investigated the potential therapeutic mechanism of LWDH Pill against postmenopausal osteoporosis in a systemic perspective. These identified multi-targets and multi-pathways provide promising directions for further revealing more exact mechanisms.

MiR-221-5p/Smad3 axis in osteoclastogenesis and its function: Potential therapeutic target for osteoporosis

OBJECTIVE

To probe the role of miR-221-5p in osteoclastogenesis and the underlying mechanism.

METHODS

Serum from patients with postmenopausal osteoporosis and healthy controls was collected for determination of miR-221-5p expression. For in vitro experiment, RAW264.7 macrophages, in which the expression of miR-221-5p and/or Smad3 was altered, were induced by RANKL to differentiate into osteoclasts. For in vivo experiment, ovariectomy was performed to construct osteoporosis mouse models, followed by tail vein injection of miR-221-5p agomir. qRT-PCR and/or western blot were applied to measure the expression of miR-221-5p, Smad3, and osteoclastogenesis-related genes (NFATc1 and TRAF6). TRAP staining was utilized for assessment of osteoclast formation, MTT assay for assessment of osteoclast viability, and H&E staining for observation of histomorphological changes. The targeting relationship between miR-221-5p and Smad3 was verified by dual-luciferase reporter gene assay.

RESULTS

Compared with healthy controls, patients with postmenopausal osteoporosis had decreased miR-221-5p expression and lower lumbar vertebra bone mineral density. MiR-221-5p expression was decreased and Smad3 level was increased during osteoclastogenesis. The osteoclastogenesis was suppressed by miR-221-5p and promoted by Smad3, as evidenced by diminished number and viability of osteoclasts following overexpression of miR-221-5p or knockdown of Smad3. MiR-221-5p negatively mediated Smad3 expression. Smad3 suppression nullified the pro-osteoclastogenesis effect of miR-221-5p inhibition. Consistent results were observed in osteoporosis mouse models.

CONCLUSION

MiR-221-5p may alleviate postmenopausal osteoporosis through suppressing osteoclastogenesis via Smad3, which provides new ideas for molecule-targeted therapy of osteoporosis.

Physical Activity and Post-Transcriptional Regulation of Aging Decay: Modulation of Pathways in Postmenopausal Osteoporosis

Background and Objectives: Bones and the skeletal muscle play a key role in human physiology as regulators of metabolism in the whole organism. Bone tissue is identified as a complex and dynamic living unit that could react to physical activity. Hormones, growth factors, signaling factors, and environmental factors control osteogenesis, and it could be regulated at a post-transcriptional level. MicroRNAs (miRNAs) can interfere with mRNAs translation. Increasing data suggest that miRNAs, through different pathways, are involved in the regulation of bone marrow mesenchymal stem cells (BMSCs) differentiation and physical activity-induced bone remodeling. The purpose of this narrative review is to investigate the potential protective role played by physical activity in affecting miRNAs expression in close tissues and elaborate on the complex network of interplay that could drive various metabolic responses of the bone to physical activity. Materials and Methods: A bibliographic search of the scientific literature was carried out in scientific databases to investigate the possible effect of physical activity on age-related features detected in the musculoskeletal system. Results: Several studies suggested that the musculoskeletal system interacting at a biomolecular level could establish crosstalk between bone and muscle in an endocrine or paracrine way through myokines released by muscle at the periosteal interface or in the bloodstream, such as irisin. Mechanical stimuli have a key role in bone formation and resorption, increasing osteogenesis and downregulating adipogenesis of BMSC via regulation of expression of runt-related transcription factor 2 (Runx2) and peroxisome proliferator-activated receptor gamma (PPARγ), respectively. Conclusions: Increasing data suggest that miRNAs, through different pathways, are involved in the regulation of BMSCs differentiation and physical activity-induced bone remodeling. Modulation of miRNAs following physical exercise represents an interesting field of investigation since these non-coding RNAs may be considered defenders against degenerative diseases and as well as useful prognostic markers in skeletal and muscle-skeletal diseases, such as osteoporosis.

Mechanical-Stress-Related Epigenetic Regulation of ZIC1 Transcription Factor in the Etiology of Postmenopausal Osteoporosis

Mechanical loading exerts a profound influence on bone density and architecture, but the exact mechanism is unknown. Our study shows that expression of the neurological transcriptional factor zinc finger of the cerebellum 1 (ZIC1) is markedly increased in trabecular bone biopsies in the lumbar spine compared with the iliac crest, skeletal sites of high and low mechanical stress, respectively. Human trabecular bone transcriptome analyses revealed a strong association between ZIC1 mRNA levels and gene transcripts characteristically associated with osteoblasts, osteocytes and osteoclasts. This supposition is supported by higher ZIC1 expression in iliac bone biopsies from postmenopausal women with osteoporosis compared with age-matched control subjects, as well as strongly significant inverse correlation between ZIC1 mRNA levels and BMI-adjusted bone mineral density (BMD) (Z-score). ZIC1 promoter methylation was decreased in mechanically loaded vertebral bone compared to unloaded normal iliac bone, and its mRNA levels correlated inversely with ZIC1 promoter methylation, thus linking mechanical stress to epigenetic control of gene expression. The findings were corroborated in cultures of rat osteoblast progenitors and osteoblast-like cells. This study demonstrates for the first time how skeletal epigenetic changes that are affected by mechanical forces give rise to marked alteration in bone cell transcriptional activity and translate to human bone pathophysiology.

Identification of Serum Exosome-Derived circRNA-miRNA-TF-mRNA Regulatory Network in Postmenopausal Osteoporosis Using Bioinformatics Analysis and Validation in Peripheral Blood-Derived Mononuclear Cells

BACKGROUND

Osteoporosis is one of the most common systemic metabolic bone diseases, especially in postmenopausal women. Circular RNA (circRNA) has been implicated in various human diseases. However, the potential role of circRNAs in postmenopausal osteoporosis (PMOP) remains largely unknown. The study aims to identify potential biomarkers and further understand the mechanism of PMOP by constructing a circRNA-associated ceRNA network.

METHODS

The PMOP-related datasets GSE161361, GSE64433, and GSE56116 were downloaded from the Gene Expression Omnibus (GEO) database and were used to obtain differentially expressed genes (DEGs). Gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were applied to determine possible relevant functions of differentially expressed messenger RNAs (mRNAs). The TRRUST database was used to predict differential transcription factor (TF)-mRNA regulatory pairs. Afterwards, combined CircBank and miRTarBase, circRNA-miRNA as well as miRNA-TF pairs were constructed. Then, a circRNA-miRNA-TF-mRNA network was established. Next, the correlation of mRNAs, TFs, and PMOP was verified by the Comparative Toxicogenomics Database. And expression levels of key genes, including circRNAs, miRNAs, TFs, and mRNAs in the ceRNA network were further validated by quantitative real-time PCR (qRT-PCR). Furthermore, to screen out signaling pathways related to key mRNAs of the ceRNA network, Gene Set Enrichment Analysis (GSEA) was performed.

RESULTS

A total of 1201 DE mRNAs, 44 DE miRNAs, and 1613 DE circRNAs associated with PMOP were obtained. GO function annotation showed DE mRNAs were mainly related to inflammatory responses. KEGG analysis revealed DE mRNAs were mainly enriched in osteoclast differentiation, rheumatoid arthritis, hematopoietic cell lineage, and cytokine-cytokine receptor interaction pathways. We first identified 26 TFs and their target mRNAs. Combining DE miRNAs, miRNA-TF/mRNA pairs were obtained. Combining DE circRNAs, we constructed the ceRNA network contained 6 circRNAs, 4 miRNAs, 4 TFs, and 12 mRNAs. The expression levels of most genes detected by qRT-PCR were generally consistent with the microarray results. Combined with the qRT-PCR validation results, we eventually identified the ceRNA network that contained 4 circRNAs, 3 miRNAs, 3 TFs, and 9 mRNAs. The GSEA revealed that 9 mRNAs participate in many important signaling pathways, such as "olfactory transduction", "T cell receptor signaling pathway", and "neuroactive ligand-receptor interaction". These pathways have been reported to the occurrence and development of PMOP. To sum up, key mRNAs in the ceRNA network may participate in the development of osteoporosis by regulating related signal pathways.

CONCLUSIONS

A circRNA-associated ceRNA network containing TFs was established for PMOP. The study may help further explore the molecular mechanisms and may serve as potential biomarkers or therapeutic targets for PMOP.

Complement receptor C5aR1 on osteoblasts regulates osteoclastogenesis in experimental postmenopausal osteoporosis

In recent years, evidence has accumulated that the complement system, an integral part of innate immunity, may be involved in the regulation of bone homeostasis as well as inflammatory bone loss, for example, in rheumatoid arthritis and periodontitis. Complement may also contribute to osteoporosis development, but investigation of the mechanism is limited. Using mice with a conditional deletion of the complement anaphylatoxin receptor C5aR1, we here demonstrated that C5aR1 in osteoblasts (C5aR1 ^Runx2-Cre mice) or osteoclasts (C5aR1 ^LysM-Cre mice) did not affect physiological bone turnover or age-related bone loss in either sex, as confirmed by micro-computed tomography, histomorphometry, and biomechanical analyses of the bone and by the measurement of bone turnover markers in the blood serum. When female mice were subjected to ovariectomy (OVX), a common model for postmenopausal osteoporosis, significant bone loss was induced in C5aR1 ^fl/fl and C5aR1 ^LysM-Cre mice, as demonstrated by a significantly reduced bone volume fraction, trabecular number and thickness as well as an increased trabecular separation in the trabecular bone compartment. Confirming this, the osteoclast number and the receptor activator of nuclear factor k-B (RANK) ligand (RANKL) serum level were significantly elevated in these mouse lines. By contrast, C5aR1 ^Runx2-Cre mice were protected from bone loss after OVX and the serum RANKL concentration was not increased after OVX. These data suggested that bone cell-specific C5aR1 may be redundant in bone homeostasis regulation under physiological conditions. However, C5aR1 on osteoblasts was crucial for the induction of bone resorption under osteoporotic conditions by stimulating RANKL release, whereas C5aR1 on osteoclasts did not regulate OVX-induced bone loss. Therefore, our results implicate C5aR1 on osteoblasts as a potential target for treating postmenopausal osteoporosis.

Does the VDR gene polymorphism influence the efficacy of denosumab therapy in postmenopausal osteoporosis?

INTRODUCTION

One of the challenges of personalized medicine is a departure from traditional pharmacology toward individualized, genotype-based therapies. Postmenopausal osteoporosis is a prevalent condition requiring intensive treatment, whose effects are measurable only after a long time, and the goal is bone fracture prevention. This study aimed to determine the influence of VDR gene variation on anti-osteoporotic one-year treatment with denosumab in 63 Polish women with postmenopausal osteoporosis.

MATERIALS AND METHODS

The correlation between bone mineral density (BMD) of the lumbar vertebral column (L1-L4) and femoral neck, and genotype distributions for the ApaI, BsmI, FokI, and TaqI variants of the VDR gene was analyzed. Bone fractures during denosumab therapy were also investigated.

RESULTS

In the case of the Bsml polymorphism, female patients with BB and Bb genotypes had statistically significantly higher values of BMD and T-score/Z-score indicators, which persisted after a year of denosumab treatment. Our results indicated that the Bsml polymorphism contributes to better bone status, and, consequently, to more efficient biological therapy. The study did not reveal significant differences between changes (delta) in BMD and genotypes for the analyzed VDR gene loci. In the entire study group, one bone fracture was observed in one patient throughout the yearlong period of denosumab therapy.

CONCLUSIONS

BB and Bb genotypes of the Bsml polymorphism of the VDR gene determine higher DXA parameter values both before and after one-year denosumab therapy in postmenopausal women with osteoporosis.

No genetic causal association between iron status and osteoporosis: A two-sample Mendelian randomization

OBJECTIVE

To explore the genetic causal association between osteoporosis (OP) and iron status through Mendelian randomization (MR).

METHODS

Publicly available genome-wide association study (GWAS) summary data were used for MR analysis with four iron status-related indicators (ferritin, iron, total iron binding capacity, and transferrin saturation) as exposures and three different types of OP (OP, OP with pathological fracture, and postmenopausal OP with pathological fracture) as outcomes. The inverse-variance weighted (IVW) method was used to analyze the genetic causal association between the four indicators of iron status and OP. The heterogeneity of MR results was determined using IVW and MR-Egger methods. The pleiotropy of MR results was determined using MR-Egger regression. A leave-one-SNP-out test was performed to determine whether the MR results were affected by a single nucleotide polymorphism (SNP). The weighted median method was conducted to further validate our results.

RESULTS

Based on IVW, MR-Egger and weighted median models, we found no causal association between iron status (ferritin, iron, total iron binding capacity, or transferrin saturation) and OP (P_beta > 0.05 in all models). IVW and MR-Egger analysis of OP with pathological fracture and iron status indicators showed no potential genetic causal association (P_beta> 0.05 in the two analyses). The results of the weighted median were consistent with those of IVW (P_beta> 0.05 in all analyses). There was no potential genetic causal association between iron status and postmenopausal OP with pathological fracture based on serum iron (P_beta>0.05 in all models). No heterogeneity or horizontal pleiotropy was found in any of the analyses. None of the leave-one-out tests in the analyses found any SNP that could affect the results of MR.

CONCLUSION

Our results demonstrate that there is no genetic causal association between OP and iron status, but the effects of other factors were not excluded.

MicroRNAs in Serum Exosomes as Circulating Biomarkers for Postmenopausal Osteoporosis

Postmenopausal osteoporosis (PMOP) is the most common skeletal disease in postmenopausal women and has become a global public health issue. Emerging evidence demonstrated the important relationship between microRNAs and PMOP. However, miRNAs have not yet been reported in PMOP. Hence, the present study aimed to investigate the differences in miRNA expression profiles in PMOP with fragility fractures to identify the key circulating miRNAs in serum exosomes and to validate these molecules as potential biomarkers. Postmenopausal women with osteoporotic fracture and normal bone mass were enrolled. Serum exosomes were isolated by traditional differential ultracentrifugation from participants. Isolated exosomes were identified by electron microscopy, western blotting and nanoparticle-tracking analysis and then examined for exosomal small RNA sequencing. The expression of miRNAs was compared by sRNA deep sequencing and bioinformatics analysis. Three miRNAs (mir-324-3p, mir-766-3p and mir-1247-5p) were found to be associated with BMD of L1-L4, FN (femur neck) and TH (total hip), while mir-330-5p and mir-3124-5p were associated with BMD of FN and TH. Furthermore, mir-330-5p was found to promote the ALP activity of hBMSCs, while mir-3124-5p showed the opposite result. The results showed that serum exosomal miRNAs were differentially expressed in postmenopausal osteoporosis patients with fragility fractures. Our study provides the first evidence that exosomal miRNA profiling revealed aberrant circulating miRNA in postmenopausal osteoporosis. Mir-324-3p, mir-766-3p, mir-1247-5p, mir-330-5p and mir-3124-5p, which were associated with bone mineral density (BMD), may serve as candidate diagnostic biomarkers as well as potentially contribute to pathophysiology of PMOP.

Identification of ferroptosis-associated biomarkers for the potential diagnosis and treatment of postmenopausal osteoporosis

OBJECTIVE

Postmenopausal osteoporosis (PMOP) is one of the most commonly occurring conditions worldwide and is characterized by estrogen deficiency as well as persistent calcium loss with age. The aim of our study was to identify significant ferroptosis-associated biomarkers for PMOP.

METHODS AND MATERIALS

We obtained our training dataset from the Gene Expression Omnibus (GEO) database using GSE56815 expression profiling data. Meanwhile, we extracted ferroptosis-associated genes for further analysis. Differentially expressed ferroptosis-associated genes (DEFAGs) between OP patients and normal controls were selected using the "limma" package. We established a ferroptosis-associated gene signature using training models, specifically, random forest (RF) and support vector machine (SVM) models. It was further validated in another dataset (GSE56814) which also showed a high AUC: 0.98, indicating high diagnostic value. Using consensus clustering, the OP patient subtypes were identified. A ferroptosis associated gene (FAG)-Scoring scheme was developed by PCA. The important candidate genes associated with OP were also compared between different ferrclusters and geneclusters.

RESULTS

There were significant DEFAGs acquired, of which five (HMOX1, HAMP, LPIN1, MAP3K5, FLT3) were selected for establishing a ferroptosis-associated gene signature. Analyzed from the ROC curve, our established RF model had a higher AUC value than the SVM model (RF model AUC:1.00). Considering these results, the established RF model was chosen to be the most appropriate training model. Later, based on the expression levels of the five DEFAGs, a clinical application nomogram was established. The OP patients were divided into two subtypes (ferrcluster A, B and genecluster A, B, respectively) according to the consensus clustering method based on DEFAGs and differentially expressed genes (DEGs). Ferrcluster B and genecluster B had higher ferroptosis score than ferrcluster A and genecluster A, respectively. The expression of COL1A1 gene was significantly higher in ferrcluster B and gencluster B compared with ferrcluster A and gencluster A, respectively, while there is no statistical difference in term of VDR gene, COL1A2 genes, and PTH gene expressions between ferrcluster A and B, together with gencluster A and B.

CONCLUSIONS

On the basis of five explanatory variables (HMOX1, HAMP, LPIN1, MAP3K5 and FLT3), we developed a diagnostic ferroptosis-associated gene signature and identified two differently categorized OP subtypes that may potentially be applied for the early diagnosis and individualized treatment of PMOP. The ER gene, VDR gene, IL-6 gene, COL1A1 and COL1A2 genes, and PTH gene are important candidate gene of OP, however, more studies are still anticipated to further elucidate the relationship between these genes and ferroptosis in OP.

The cross-talk between matrix metalloproteinase-9, RANKL/OPG system and cardiovascular risk factors in ovariectomized rat model of postmenopausal osteoporosis

Epidemiology and pathogenesis of cardiovascular diseases (CVD) and osteoporosis are strikingly overlapping. This study presents matrix metalloproteinase-9 (MMP-9), as a simple molecular link more consistently associated with the pathophysiology of both osteoporosis and CVD risk factors. 40 adult female rats were randomly distributed into 4 groups [control sham-operated, untreated osteoporosis, carvedilol-treated osteoporosis and alendronate-treated osteoporosis]. After 8 weeks, blood samples were collected to estimate Lipid profile (Total cholesterol, HDL, Triglycerides), inflammatory markers (IL-6, TNF alpha, CRP and NO), and Bone turnover markers (BTM) (Alkaline phosphatase, osteocalcin and pyridinoline). The tibias were dissected to estimate MMP-9 and NF-kB gene expression, OPG, RANKL levels and for histological examination. Induction of osteoporosis resulted in a significant elevation in BTM, inflammatory markers and dyslipidemia. MMP-9 was significantly elevated and positively correlated with BTM, inflammation and dyslipidemia markers. Carvedilol and alendronate exerted a bone preservative role and attenuated dyslipidaemia and inflammation in accordance with their respective effect on MMP-9.

Identification and comparison of novel circular RNAs with associated co-expression and competing endogenous RNA networks in postmenopausal osteoporosis

BACKGROUND

Circular RNAs (circRNAs) are emerging as crucial regulators in various human diseases. So far, the expression profile and regulatory mechanism of circRNAs in postmenopausal osteoporosis (PMOP) are less studied and should be deciphered urgently. Herein, we aimed to reveal key circRNAs affecting PMOP and clarify their compounding regulatory actions.

METHODS

To reveal key circRNAs affecting PMOP and clarify their compounding regulatory actions, whole transcriptome sequencing and bioinformatics analysis were performed to identify differentially expressed circRNAs (DECs). The expression pattern and regulatory networks of DECs in peripheral blood mononuclear cells (PBMCs) were unearthed.

RESULTS

A total of 373 DECs comprising 123 intronic, 100 antisense, 70 exonic, 55 intergenic, and 25 sense-overlapping circRNAs were identified. Among these, 73 circRNAs were upregulated and 300 were downregulated. These DECs exerted pivotal functions in the pathogenesis of PMOP as demonstrated by Gene Ontology (GO) annotation and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The circRNA-miRNA-mRNA co-expression network comprising 28 DECs, 145 miRNAs, and 175 differentially expressed mRNAs predicted the possible mechanism of the pathogenesis and progression of PMOP.

CONCLUSION

The results of the present study provided a further comprehension of circRNA-associated competing endogenous RNA regulatory mechanism in PMOP. The steadily expressed and disease-specific DECs may serve as promising diagnostic and prognostic biomarkers for PMOP.

Oxidative Stress-Induced Hypermethylation of KLF5 Promoter Mediated by DNMT3B Impairs Osteogenesis by Diminishing the Interaction with β-Catenin

Aims: Emerging evidence suggests that the pathogenesis of osteoporosis, characterized by impaired osteogenesis, is shifting from estrogen centric to oxidative stress. Our previous studies have shown that the zinc-finger transcription factor krüppel-like factor 5 (KLF5) plays a key role in the degeneration of nucleus pulposus and cartilage. However, its role in osteoporosis remains unknown. We aimed to investigate the effect and mechanism of KLF5 on osteogenesis under oxidative stress. Results: First, KLF5 was required for osteogenesis and stimulated osteogenic differentiation of bone marrow mesenchymal stem cells (BMSCs). KLF5 was hypermethylated and downregulated in ovariectomy-induced osteoporosis mice and in BMSCs treated with H2O2. Interestingly, DNA methyltransferases 3B (DNMT3B) upregulation mediated the hypermethylation of KLF5 induced by oxidative stress, thereby impairing osteogenic differentiation. The inhibition of KLF5 hypermethylation using DNMT3B siRNA or 5-AZA-2-deoxycytidine (5-AZA) protected osteogenic differentiation of BMSCs from oxidative stress. Regarding the downstream mechanism, KLF5 induced β-catenin expression. More importantly, KLF5 promoted the nuclear translocation of β-catenin, which was mediated by the armadillo repeat region of β-catenin. Consistently, oxidative stress-induced KLF5 hypermethylation inhibited osteogenic differentiation by reducing the expression and nuclear translocation of β-catenin. Innovation: We describe the novel effect and mechanism of KLF5 on osteogenesis under oxidative stress, which is linked to osteoporosis for the first time. Conclusion: Our results suggested that oxidative stress-induced hypermethylation of KLF5 mediated by DNMT3B impairs osteogenesis by diminishing the interaction with β-catenin, which is likely to contribute to osteoporosis. Targeting the hypermethylation of KLF5 might be a new strategy for the treatment of osteoporosis. Antioxid. Redox Signal. 35, 1-20.