To examine the genome-wide distribution of hydroxymethylated cytosine (5hmC) in osteoarthritic (OA) and normal chondrocytes in order to investigate the effect on OA-specific gene expression.
Cartilage was obtained from OA patients undergoing total knee arthroplasty or from control patients undergoing anterior cruciate ligament reconstruction. Genome-wide sequencing of 5hmC-enriched DNA was performed in a small cohort of normal and OA chondrocytes to identify differentially hydroxymethylated regions (DhMRs) in OA chondrocytes. Data from the genome-wide sequencing of 5hmC-enriched DNA were intersected with global OA gene expression data to define subsets of genes and pathways potentially affected by increased 5hmC levels in OA chondrocytes.
A total of 70,591 DhMRs were identified in OA chondrocytes as compared to normal chondrocytes, 44,288 (63%) of which were increased in OA chondrocytes. The majority of DhMRs (66%) were gained in gene bodies. Increased DhMRs were observed in ∼50% of genes previously implicated in OA pathology including MMP3, LRP5, GDF5, and COL11A1. Furthermore, analyses of gene expression data revealed gene body gain of 5hmC appears to be preferentially associated with activated, but not repressed, genes in OA chondrocytes.
This study provides the first genome-wide profiling of 5hmC distribution in OA chondrocytes. We had previously reported a global increase in 5hmC levels in OA chondrocytes. Gain of 5hmC in the gene body is found to be characteristic of activated genes in OA chondrocytes, highlighting the influence of 5hmC as an epigenetic mark in OA. In addition, this study identifies multiple OA-associated genes that are potentially regulated either singularly by gain of DNA hydroxymethylation or in combination with loss of DNA methylation.