Identification of key mRNAs and microRNAs in the pathogenesis and progression of osteoarthritis using microarray...

Osteoarthritis (OA) is a common type of disease affecting the joints that results from the breakdown of joint cartilage and the underlying bone; currently, its pathogenesis is still unclear. The aim of the present study was to identify key mRNAs and miRNAs involved in the pathogenesis and progression of OA using microarray analysis. The gene expression profile of GSE27492 was downloaded from the Gene Expressed Omnibus database, and included 49 arthritic mouse ankle samples collected at 6 time points (0, 1, 3, 7, 12 and 18 days) following the induction of arthritis via serum transfer. Differentially expressed genes (DEGs) were identified in ankle samples taken on days 1, 3, 7, 12 and 18 following serum transfer compared with day 0 samples, and overlapping DEGs in day 3, 7, 12 and 18 samples were identified. The Database for Annotation, Visualization and Integrated Discovery online tool was used to perform functional and pathway enrichment analyses of the overlapping DEGs. The miRWalk database was used to identify potential micro (mi) RNAs regulating the selected overlapping DEGs, and regulatory miRNA‑target mRNA pairs were obtained. The Cytoscape platform was used to establish and visualize the miRNA‑mRNA regulatory network. The present results revealed that 35, 103, 62 and 75 DEGs were identified in day 3, 7, 12 and 18 samples, respectively. A total of 17 overlapping DEGs were identified among the 4 sample sets, and revealed to be enriched in 14 gene ontology terms and 3 Kyoto Encyclopedia of Genes and Genomes pathways. miRWalk analysis identified 242 potential miRNA‑mRNA regulatory pairs and 211 nodes were revealed to be involved in the miRNA‑mRNA regulatory network. The present study identified potential genes, including C‑type lectin domain family 4 member D, chemokine (C‑X‑C motif) ligand 1 and C‑C motif chemokine ligand, and pathways, including chemokine signaling pathways, cytokine‑cytokine receptor interactions and nucleotide‑binding oligomerization domain‑like receptor signaling pathways, which may be involved in the pathogenesis and progression of OA. These findings may help elucidate the molecular mechanisms underlying OA pathophysiology, and may be useful for the development of novel therapeutic targets for the treatment of patients with OA.