5'-Hydroxymethylcytosine Precedes Loss of CpG Methylation in Enhancers and Genes Undergoing Activation in Cardiomyocyte Maturation

Cardiomyocytes undergo major changes in DNA methylation during maturation and transition to a non-proliferative state after birth. 5'-hydroxylation of methylated cytosines (5hmC) is not only involved in DNA loss of CpG methylation but is also thought to be an epigenetic mark with unique distribution and functions. Here, we sought to get insight into the dynamics of 5'-hydroxymethylcytosine in newborn and adult cardiomyocytes.

Simultaneous analysis of 5hmC-Seq, MethylC-Seq, RNA-Seq and ChIP-Seq data at two defined time points of cardiomyocyte maturation demonstrates that 5hmC is positively associated with gene expression and decorates sites of subsequent DNA loss of CpG methylation.

Cardiomyocyte nuclei from newborn and adult C57BL/6 mice were purified by flow cytometric sorting. 5hmC-containing DNA was captured by selective chemical labeling, followed by deep sequencing. Sequencing reads of library replicates were mapped independently (n = 3 for newborn, n = 2 for adult mice) and merged for further analysis steps. 5hmC coverage was normalized to read length and the total number of mapped reads (RPKM). MethylC-Seq, ChIP-Seq and RNA-Seq data sets of newborn and adult cardiomyocytes served to elucidate specific features of 5hmC at gene bodies and around low methylated regions (LMRs) representing regulatory genomic regions with enhancer function.

163,544 and 315,220 5hmC peaks were identified in P1 and adult cardiomyocytes, respectively. Of these peaks, 66,641 were common between P1 and adult cardiomyocytes with more than 50% reciprocal overlap. P1 and adult 5hmC peaks were overrepresented in genic features such as exons, introns, 3'- and 5'-untranslated regions (UTRs), promotors and transcription end sites (TES). During cardiomyocyte maturation, 5hmC was found to be enriched at sites of subsequent DNA loss of CpG methylation such as gene bodies of upregulated genes (i.e. Atp2a2, Tnni3, Mb, Pdk4). Additionally, centers of postnatally established enhancers were premarked by 5hmC before DNA loss of CpG methylation. Conclusions: Simultaneous analysis of 5hmC-Seq, MethylC-Seq, RNA-Seq and ChIP-Seq data at two defined time points of cardiomyocyte maturation demonstrates that 5hmC is positively associated with gene expression and decorates sites of subsequent DNA loss of CpG methylation.