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.

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

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

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

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

The integrative analysis of competitive endogenous RNA regulatory networks in osteoporosis

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

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

circ_0134120: a new frontier in understanding postmenopausal osteoporosis pathogenesis

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

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

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

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

CircRNA_28313/miR-195a/CSF1 axis modulates osteoclast differentiation to affect OVX-induced bone absorption in mice

The roles of circRNA–miRNA–mRNA networks in the development and treatment of osteoporosis

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

Identification of circRNA Expression Profiles in BMSCs from Glucocorticoid-Induced Osteoporosis Model

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

Identification of the canonical and noncanonical role of miR-143/145 in estrogen-deficient bone loss

Exosomal miR-186 derived from BMSCs promote osteogenesis through hippo signaling pathway in postmenopausal osteoporosis

Exosomal circFAM63Bsuppresses bone regeneration of postmenopausal osteoporosis via regulating miR‐578/HMGA2 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