Identification of genes associated with osteoarthritis by microarray analysis

The aim of the present study was to investigate the mechanisms of osteoarthritis (OA). Raw microarray data (GSE51588) were downloaded from Gene Expression Omnibus, including samples from OA (n=20) and non‑OA (n=5) knee lateral and medial tibial plateaus. Differentially expressed genes (DEGs) were identified using Student's t‑test. Functional and pathway enrichment analyses were performed for the upregulated and downregulated DEGs. A protein‑protein interaction network (PPI) was constructed according to the Search Tool for the Retrieval of Interacting Genes/Proteins database, and module analysis of the PPI network was performed using CFinder. The protein domain enrichment analysis for genes in modules was performed using the INTERPRO database. A total of 869 upregulated and 508 downregulated DEGs were identified. The enriched pathways of downregulated and upregulated DEGs were predominantly associated with the cell cycle (BUB1, BUB1B, CCNA2, CCNB1 and CCNE1), and extracellular matrix (ECM)‑receptor interaction (CD36, COL11A2, COL1A1, COL2A1 and COL3A1). Functional enrichment analysis of the DEGs demonstrated that FGF19, KIF11 and KIF2C were involved in the response to stress and that ACAN, ADAMTS10 and BGN were associated with proteinaceous ECM. The top protein domain was IPR001752: Kinesin motor region involving three genes (KIF2C, KIF11 and KIF20A). The identified DEGs, including KIF2C, KIF11 and KIF20A, may be significant in the pathogenesis of OA.