Aberrantly hydroxymethylated differentially expressed genes and the associated protein pathways in osteoarthritis

Background The elderly population is at risk of osteoarthritis (OA), a common, multifactorial, degenerative joint disease. Environmental, genetic, and epigenetic (such as DNA hydroxymethylation) factors may be involved in the etiology, development, and pathogenesis of OA. Here, comprehensive bioinformatic analyses were used to identify aberrantly hydroxymethylated differentially expressed genes and pathways in osteoarthritis to determine the underlying molecular mechanisms of osteoarthritis and susceptibility-related genes for osteoarthritis inheritance. Methods Gene expression microarray data, mRNA expression profile data, and a whole genome 5hmC dataset were obtained from the Gene Expression Omnibus repository. Differentially expressed genes with abnormal hydroxymethylation were identified by MATCH function. Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the genes differentially expressed in OA were performed using Metascape and the KOBAS online tool, respectively. The protein–protein interaction network was built using STRING and visualized in Cytoscape, and the modular analysis of the network was performed using the Molecular Complex Detection app. Results In total, 104 hyperhydroxymethylated highly expressed genes and 14 hypohydroxymethylated genes with low expression were identified. Gene ontology analyses indicated that the biological functions of hyperhydroxymethylated highly expressed genes included skeletal system development, ossification, and bone development; KEGG pathway analysis showed enrichment in protein digestion and absorption, extracellular matrix–receptor interaction, and focal adhesion. The top 10 hub genes in the protein–protein interaction network were COL1A1, COL1A2, COL2A1, COL3A1, COL5A1, COL5A2, COL6A1, COL8A1, COL11A1, and COL24A1. All the aforementioned results are consistent with changes observed in OA. Conclusion After comprehensive bioinformatics analysis, we found aberrantly hydroxymethylated differentially expressed genes and pathways in OA. The top 10 hub genes may be useful hydroxymethylation analysis biomarkers to provide more accurate OA diagnoses and target genes for treatment of OA.