Exploring the Key Genes and Pathways of Osteoarthritis in Knee Cartilage in a Rat Model Using Gene Expression Profiling

PURPOSE: To compare differentially expressed genes (DEGs) mediating osteoarthritis (OA) in knee cartilage and in normal knee cartilage in a rat model of OA and to identify their impact on molecular pathways associated with OA. MATERIALS AND METHODS: A gene expression profile was downloaded from the Gene Expression Omnibus database. Analysis of DEGs was carried out using GEO2R. Enrichment analyses were performed on the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway using the Search Tool for the Retrieval of Interacting Genes database (http://www.string-db.org/). Subsequently, the regulatory interaction network of OA-associated genes was visualized using Cytoscape software (version 3.4.0; www.cytoscape.org). RESULTS: In the gene expression profile GSE103416, a total of 99 DEGs were identified. Among them, 76 DEGs (76.77%) were overexpressed, and the remaining 23 DEGs (23.23%) were underexpressed. GO and pathway enrichment analyses of target genes were performed. Using gene-gene interaction network analysis, relevant core genes, including MET, UBB, GNAI3, and GNA13, were shown to hold a potential relationship with the development of OA in cartilage. Using quantitative real-time PCR, the Gna13/cGMP-PKG signaling pathway was identified as a potential research target for therapy and for further understanding the development of OA. CONCLUSION: The results of the present study provide a comprehensive understanding of the roles of DEGs in knee cartilage in relation to the development of OA.