Abstract
OBJECTIVE: This study is aimed at screening the differential expression profiles of mRNA under weightlessness osteoporosis through high-throughput sequencing technology, as well as investigating the pathogenesis of weightlessness osteoporosis at the molecular level especially in bone marrow mesenchymal stem cells (BMSCs). METHODS: The mouse bone marrow mesenchymal stem cell line was divided into ground group and simulated microgravity (SMG) group. BMP-2 was used to induce osteogenic differentiation, and SMG group was placed into 2D-gyroscope to simulate weightless condition. Transcriptome sequencing was performed by Illumina technology, DEGs between ground and SMG group was conducted using the DEseq2 algorithm. Molecular functions and signaling pathways enriched by DEGs were then comprehensively analyzed via multiple bioinformatic approaches including but not limited to GO, KEGG, GSEA, and PPI analysis. RESULTS: A total of 263 DEGs were identified by comparing these 2 groups, including 186 upregulated genes and 77 downregulated genes. GO analysis showed that DEGs were enriched in osteoblasts, osteoclasts cell proliferation, differentiation, and apoptosis; KEGG analysis revealed that DEGs were significantly enriched in the TNF signaling pathway and FoxO signaling pathway; the enrichment results from Reactome database displayed that DEGs were mainly involved in the transcription of Hoxb3 gene, RUNX1 recruitment KMT2A gene, and activation of Hoxa2 chromatin signaling pathway. The four genes, IL6, CXCR4, IGF1, and PLOD2, were identified as hub genes for subsequent analysis. CONCLUSIONS: This study elucidated the significance of 10 hub genes in the development of weightlessness osteoporosis. In addition, the results of this study provide a theoretical basis and novel ideas for the subsequent research of the pathogenesis and clinical treatment of weightlessness osteoporosis.