Our study aimed to explore more mechanistic insights into the epigenetic regulation of osteoarthritis (OA).
The expression profiles (accession number: GSE64393 and GSE64394) were downloaded from the Gene Expression Omnibus database. The differentially hydroxymethylated regions (DhMRs) and differentially expressed genes (DEGs) between OA and control groups were identified. The distribution of DhMRs in the whole genome and the correlation between DhMRs and DEGs were analyzed. Functional module mining for the DEGs and DhMRs was conducted, followed by protein-protein interaction (PPI) analysis. The transcriptional factor (TF) was predicted.
Total 52,282 DhMRs were obtained, among which 31,452 ones were annotated to 9726 genes. Additionally, 1806 DEGs were selected. Hydroxymethylation mainly occurred in gene body region. Correlation analysis revealed that more than 70% of DhMRs were uncorrelated with DEGs expression. Functional module mining for the DEGs and DhMRs identified 2 functional modules, which were involved in pathways of regulation of actin cytoskeleton, and TGF-β signaling pathway. A PPI network was constructed, and ITGB3 had the highest degree. Furthermore, 7 TFs were predicted, which regulated 12 candidate genes, such as HES1-PTEN.
The onset and progression of OA may be associated with the upregulated hydroxymethylation in gene body region of PTEN. HES1 may be important TF in the pathogenesis of OA. Additionally, pathways of regulation of actin cytoskeleton, and TGF-beta signaling pathway may also play important roles in OA progression.