Identification of genes and pathways associated with subchondral bone in osteoarthritis via bioinformatic analysis

Osteoarthritis (OA) is a high prevalent musculoskeletal problem, which can cause severe pain, constitute a huge social and economic burden, and seriously damage the quality of life. This study was intended to identify genetic characteristics of subchondral bone in patients with OA and to elucidate the potential molecular mechanisms involved. Data of gene expression profiles (GSE51588), which contained 40 OA samples and 10 normal samples, was obtained from the Gene Expression Omnibus (GEO). The raw data were integrated to obtain differentially expressed genes (DEGs) and were further analyzed with bioinformatic analysis. The protein–protein interaction (PPI) networks were built and analyzed via Search Tool for the Retrieval of Interacting Genes (STRING). The significant modules and hub genes were identified via Cytoscape. Moreover, Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analysis were performed. Totally 235 DEGs were differentially expressed in the subchondral bone from OA patients compared with those of normal individuals, of which 78 were upregulated and 157 were downregulated. Eight hub genes were identified, including DEFA4, ARG1, LTF, RETN, PGLYRP1, OLFM4, ORM1, and BPI. The enrichment analyses of the DEGs and significant modules indicated that DEGs were mainly involved in inflammatory response, extracellular space, RAGE receptor binding, and amoebiasis pathway. The present study provides a novel and in-depth understanding of pathogenesis of the OA subchondral bone at molecular level. DEFA4, ARG1, LTF, RETN, PGLYRP1, OLFM4, ORM1, and BPI may be the new candidate targets for diagnosis and therapies on patients with OA in the future.