Identification of key genes and construction of CircRNA–miRNA–mRNA regulatory networks in osteoarthritis

Abstract

Background: Osteoarthritis (OA) is one of the most frequent degenerative joint diseases with high rate of disability, but its mode of action remains largely unclear. The current study was aimed at identifying key genes and molecular mechanism in OA. Gene expression datasets (GSE1919, GSE55235, GSE55457) were downloaded from Gene Expression Omnibus for integrated bioinformatics analysis. Differentially expressed genes (DEGs) in OA synovial tissues were identified using GeoDiver and GEO2R. Gene ontology enrichment analyses were undertaken via FunRich and Metascape. Also, Gene Set Enrichment Analysis was performed using miRWalk3.0. Subsequently, pathways interrelation analysis of hub genes was carried out using plug-in ClueGO v2.3.3. Additionally, circRNA- miRNA-mRNA regulatory networks were visualized using Cytoscape.

Results: A total of 508 DEGs were obtained from three GSE datasets, of these five intersection DEGs (TNFAIP3, VEGFA, GADD45B, SIK1, KLF9) were shared by three GSE datasets. Intersection DEGs were significantly enriched in LKB1 signaling events, signaling events mediated by focal adhesion kinase, and PDGFR-beta signaling pathway. Enrichment analysis for all the DEGs showed that they mainly enriched in inflammatory response, cytokine production, blood vessel development, stress response, osteoclast differentiation, and MAPK signaling pathway. A total of 39 genes were regarded as hub genes and pathways interrelation analysis indicated that hub genes mainly enriched in TNF signaling pathway, IL-17 signaling pathway, and NF-kappa B signaling pathway.

Conclusions: The current study revealed the potential key genes, pathways, and circRNA-miRNA-mRNA regulatory networks in OA, which may contribute to a more comprehensive understanding of OA pathogenesis.