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.

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

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

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

Exosomal circFAM63Bsuppresses bone regeneration of postmenopausal osteoporosis via regulating miR‐578/HMGA2 axis

Circ_0134944 inhibits osteogenesis through miR-127-5p/PDX1/SPHK1 pathway

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

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

circ_0134120: a new frontier in understanding postmenopausal osteoporosis pathogenesis

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

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

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

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

CircHmbox1 Targeting miRNA-1247-5p Is Involved in the Regulation of Bone Metabolism by TNF-α in Postmenopausal Osteoporosis

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

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

The integrative analysis of competitive endogenous RNA regulatory networks in osteoporosis

Circulating miR-338 Cluster activities on osteoblast differentiation: Potential Diagnostic and Therapeutic Targets for Postmenopausal Osteoporosis

MicroRNA-410 Participates in the Pathological Process of Postmenopausal Osteoporosis by Downregulating Bone Morphogenetic Protein-2

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

Silencing of miR-483-5p alleviates postmenopausal osteoporosis by targeting SATB2 and PI3K/AKT pathway

miR-18a-5p promotes osteogenic differentiation of BMSC by inhibiting Notch2