Abstract: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.

Exosomal hsa_circ_0006859 is a potential biomarker for postmenopausal osteoporosis and enhances adipogenic versus osteogenic differentiation in human bone marrow mesenchymal stem cells by sponging miR-431-5p

Circ-Plod2 Destabilizes Mpo Mrna by Binding to IGF2BP2 to Promote Osteogenic Differentiation of Bone Marrow Mesenchymal Stem Cells.

Microarray analysis of circRNAs sequencing profile in exosomes derived from bone marrow mesenchymal stem cells in postmenopausal osteoporosis patients

Circular RNA_0062582 promotes osteogenic differentiation of human bone marrow mesenchymal stem cells via regulation of microRNA-145/CBFB axis

Circular RNA Hsa_circ_0006766 targets microRNA miR-4739 to regulate osteogenic differentiation of human bone marrow mesenchymal stem cells

CircRNA hsa_circ_0001421 promotes the osteoblast differentiation of human adipose mesenchymal stem cells through the miR-608/SP7 axis

Circular RNA hsa_circ_0076906 competes with OGN for miR-1305 biding site to alleviate the progression of osteoporosis

Comprehensive analysis of circRNA expression profiles in postmenopausal women differing in bone mineral density

Identification of circRNA-associated ceRNA network in BMSCs of OVX models for postmenopausal osteoporosis

Involvement of circRNA_0007059 in the regulation of postmenopausal osteoporosis by promoting the microRNA‐378/BMP‐2 axis

Circular RNA circ-3626 promotes bone formation by modulating the miR-338-3p/Runx2 axis

Circ_0062582 Promotes Osteogenic Differentiation of Human Bone Marrow Mesenchymal Stem Cells in Vitro by Elevating SMAD5 Expression Through Sponging Mir-197-3p.

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

Hsa_circ_0026827 Promotes Osteoblast Differentiation of Human Dental Pulp Stem Cells Through the Beclin1 and RUNX1 Signaling Pathways by Sponging Mir-188-3p.

CircFOXP1/FOXP1 promotes osteogenic differentiation in adipose‐derived mesenchymal stem cells and bone regeneration in osteoporosis via miR‐33a‐5p

Identification of Suitable Reference Gene and Biomarkers of Serum Mirnas for Osteoporosis

Circ_0019693 promotes osteogenic differentiation of bone marrow mesenchymal stem cell and enhances osteogenesis-coupled angiogenesis via regulating microRNA-942-5p-targeted purkinje cell protein 4 in the development of osteoporosis

Hsa_circ_0036872 has an important promotional effect in enhancing osteogenesis of dental pulp stem cells by regulating the miR-143-3p/IGF2 axis

RUNX3 derived hsa_circ_0005752 accelerates the osteogenic differentiation of adipose-derived stem cells via the miR-496/MDM2-p53 pathway

Circ_0027885 sponges miR-203-3p to regulate RUNX2 expression and alleviates osteoporosis progression

Differential Expression Profiles of Circular RNAs During Osteogenic Differentiation of Mouse Adipose-Derived Stromal Cells