Osteoporosis (OP) is a disease characterized by bone mass loss, bone microstructure damage, increased bone fragility, and easy fracture. The molecular mechanism underlying OP remains unclear. In this study, we identified 217 genes associated with OP, and formed a gene set [OP-related genes gene set (OPgset)]. The highly enriched GOs and pathways showed OPgset genes were significantly involved in multiple biological processes (skeletal system development, ossification, and osteoblast differentiation), and several OP-related pathways (Wnt signaling pathway, osteoclast differentiation, steroid hormone biosynthesis, and adipocytokine signaling pathway). Besides, pathway crosstalk analysis indicated three major modules, with first module consisted of pathways mainly involved in bone development-related signaling pathways, second module in Wnt-related signaling pathway and third module in metabolic pathways. Further, we calculated degree centrality of a node and selected ten key genes/proteins, including TGFB1, IL6, WNT3A, TNF, PTH, TP53, WNT1, IGF1, IL10, and SERPINE1. We analyze the K-core and construct three k-core sub-networks of OPgset genes. In summary, we for the first time explored the molecular mechanism underlying OP via network- and pathway-based methods, results from our study will improve our understanding of the pathogenesis of OP. In addition, these methods performed in this study can be used to explore pathogenesis and genes related to a specific disease.

An imbalance between the adipogenesis and osteogenesis of bone marrow mesenchymal stem cells (BMSCs) is considered the basic pathogenesis of osteoporosis. Although actin cytoskeleton remodelling plays a crucial role in the differentiation of BMSCs, the role of actin cytoskeleton remodelling in the adipogenesis of BMSCs and postmenopausal osteoporosis (PMOP) has remained elusive. Glia maturation factor-beta (GMFB) has a unique role in remodelling the polymerization/depolymerization cycles of actin. We observed that GMFB expression was increased in bone tissue from both ovariectomized (OVX) rats and PMOP patients. GMFB knockout inhibited the accumulation of bone marrow adipocytes and increased bone mass in the OVX rat model. The inhibition of adipocyte differentiation in GMFB knockout BMSCs was mediated via actin cytoskeleton remodelling and the Ca²⁺-calcineurin-NFATc2 axis. Furthermore, we found that GMFB shRNA treatment in vivo had favourable effects on osteoporosis induced by OVX. Together, these findings suggest a pathological association of the GMFB with PMOP and highlight the potential of the GMFB as a therapeutic target for osteoporosis patients.

Osteoporosis (OP) has the characteristics of a systematically impaired bone mass, strength, and microstructure. Long non-coding RNAs (lncRNAs) are longer than 200 nt, and their functions in osteoporosis is yet not completely understood. We first harvested the bone marrow mesenchymal stem cells (BMSCs) from ovariectomy (OVX) and sham mice. Then, we systematically analyzed the differential expressions of lncRNAs and messenger RNAs (mRNAs) and constructed lncRNA–mRNA coexpression network in order to identify the function of lncRNA in osteoporosis. Totally, we screened 743 lncRNAs (461 upregulated lncRNAs and 282 downregulated lncRNAs) and 240 mRNAs (128 upregulated and 112 downregulated) with significantly differential expressions in OP compared to normal. We conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional analyses to investigate the functions and pathways of the differential expression of messenger RNAs (mRNAs), a coexpressed network of lncRNA/mRNA. Quantitative PCR (qPCR) validated that the expressions of NONMMUT096150.1, NONMMUT083450.1, and NONMMUT029743.2 were all downregulated, whereas NONMMUT026970.2, NONMMUT051734.2, NONMMUT003617.2, and NONMMUT034049.2 were all upregulated in the OVX group. NONMMUT096150.1, as a key lncRNA in OP, was identified to modulate the adipogenesis of BMSCs. Further analysis suggested that NONMMUT096150.1 might modulate the adipogenesis of BMSCs via the peroxisome proliferator-activated receptor (PPAR) signaling pathway, AMPK signaling pathway, and the lipolysis regulation in adipocyte and adipocytokine signaling pathway. Our study expands the understanding of lncRNA in the pathogenesis of OP.