Systematic Analysis of mRNAs and ncRNAs in BMSCs of Senile Osteoporosis Patients

Yiyun Geng; Jinfu Chen; Chongfei Chang; Yifen Zhang; Li Duan; Weimin Zhu; Lisha Mou; Jianyi Xiong; Daping Wang

doi:10.3389/fgene.2021.776984

Systematic Analysis of mRNAs and ncRNAs in BMSCs of Senile Osteoporosis Patients

Front Genet. 2021 Dec 20:12:776984. doi: 10.3389/fgene.2021.776984. eCollection 2021.

Authors

Yiyun Geng^{1

2}, Jinfu Chen¹, Chongfei Chang¹, Yifen Zhang¹, Li Duan^{1

3}, Weimin Zhu¹, Lisha Mou¹, Jianyi Xiong^{1

3}, Daping Wang^{1

3}

Affiliations

¹ Shenzhen Key Laboratory of Tissue Engineering, Shenzhen Second People's Hospital (The First Hospital Affiliated to Shenzhen University), Shenzhen, China.
² School of Biotechnology and Food Engineering, Changshu Institute of Technology, Suzhou, China.
³ Guangdong Provincial Research Center for Artificial Intelligence and Digital Orthopedic Technology, Shenzhen, China.

Abstract

Senile osteoporosis (SOP) is a worldwide age-related disease characterized by the loss of bone mass and decrease in bone strength. Bone mesenchymal stem cells (BMSCs) play an important role in the pathology of senile osteoporosis. Abnormal expression and regulation of non-coding RNA (ncRNA) are involved in a variety of human diseases. In the present study, we aimed to identify differentially expressed mRNAs and ncRNAs in senile osteoporosis patient-derived BMSCs via high-throughput transcriptome sequencing in combination with bioinformatics analysis. As a result, 415 mRNAs, 30 lncRNAs, 6 circRNAs and 27 miRNAs were found to be significantly changed in the senile osteoporosis group. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were applied to analyze the function of differentially expressed mRNAs and ncRNAs. The circRNA-miRNA-mRNA regulatory network was constructed using the cytoHubba plugin based on the Cytoscape software. Interestingly, circRNA008876-miR-150-5p-mRNA was the sole predicted circRNA-miRNA-mRNA network. The differential expression profile of this ceRNA network was further verified by qRT-PCR. The biological function of this network was validated by overexpression and knockdown experiments. In conclusion, circRNA008876-miR-150-5p-mRNA could be an important ceRNA network involved in senile osteoporosis, which provides potential biomarkers and therapeutic targets for senile osteoporosis.

Keywords: ceRNA network; human bone marrow mesenchymal stem cells; non-coding RNA; senile osteoporosis; whole transcriptome sequencing.