Identification and comprehensive analysis of circRNA–miRNA–mRNA regulatory networks in osteoarthritis

Abstract

Osteoarthritis (OA) is a common orthopedic degenerative disease, leading to high disability in activities of daily living. There remains an urgent need to identify the underlying mechanisms and identify new therapeutic targets in OA diagnosis and treatment. Circular RNAs (circRNAs) play a role in the development of multiple diseases. Many studies have reported that circRNAs regulate microRNAs (miRNAs) through an endogenous competitive mechanism. However, it remains unclear if an interplay between circRNAs, miRNAs, and target genes plays a deeper regulatory role in OA. Four datasets were downloaded from the GEO database, and differentially expressed circRNAs (DECs), differentially expressed miRNAs (DEMs), and differentially expressed genes (DEGs) were identified. Functional annotation and pathway enrichment analysis of DEGs and DECs were carried out to determine the main associated mechanism in OA. A protein–protein network (PPI) was constructed to analyze the function of, and to screen out, hub DEGs in OA. Based on the artificial intelligence prediction of protein crystal structures of two hub DEGs, TOP2A and PLK1, digitoxin and oxytetracycline were found to have the strongest affinity, respectively, with molecular docking. Subsequently, overlapping DEMs and miRNAs targeted by DECs obtained target DEMs (DETMs). Intersection of DEGs and genes targeted by DEMs obtained target DEGs (DETGs). Thus, a circRNA–miRNA–mRNA regulatory network was constructed from 16 circRNAs, 32 miRNAs, and 97 mRNAs. Three hub DECs have the largest number of regulated miRNAs and were verified through in vitro experiments. In addition, the expression level of 16 DECs was validated by RT-PCR. In conclusion, we constructed a circRNA–miRNA–mRNA regulatory network in OA and three new hub DECs, hsa_circ_0027914, hsa_circ_0101125, and hsa_circ_0102564, were identified as novel biomarkers for OA.