Genome-wide methylation analysis of retrocopy-associated CpG islands and their genomic environment.

Gene duplication by retrotransposition, i.e., the reverse transcription of an mRNA and integration of the cDNA into the genome, is an important mechanism in evolution. Based on whole-genome bisulfite sequencing of monocyte DNA, we have investigated the methylation state of all CpG islands (CGIs) associated with a retrocopy (n = 1,319), their genomic environment, as well as the CGIs associated with the ancestral genes. Approximately 10% of retrocopies are associated with a CGI. Whereas almost all CGIs of the human genome are unmethylated, 68% of the CGIs associated with a retrocopy are methylated. In retrocopies resulting from multiple retrotranspositions of the same ancestral gene, the methylation state of the CGI often differs. There is a strong positive correlation between the methylation state of the CGI/retrocopy and their genomic environment, suggesting that the methylation state of the integration site determined the methylation state of the CGI/retrocopy, or that methylation of the retrocopy by a host defense mechanism has spread into the adjacent regions. Only a minor fraction of CGI/retrocopies (n = 195) has intermediate methylation levels. Among these, the previously reported CGI/retrocopy in intron 2 of the RB1 gene (PPP1R26P1) as well as the CGI associated with the retrocopy RPS2P32 identified in this study carry a maternal methylation imprint. In conclusion, these findings shed light on the evolutionary dynamics and constraints of DNA methylation.

Altering TET dioxygenase levels within physiological range affects DNA methylation dynamics of HEK293 cells.

The TET family of dioxygenases (TET1/2/3) can convert 5-methylcytosine (5 mC) into 5-hydroxymethylcytosine (5 hmC) and has been shown to be involved in active and passive DNA demethylation. Here, we demonstrate that altering TET dioxygenase levels within physiological range can affect DNA methylation dynamics of HEK293 cells. Overexpression of TET1 increased global 5 hmC levels and was accompanied by mild DNA demethylation of promoters, gene bodies and CpG islands. Conversely, the simultaneous knockdown of TET1, TET2, and TET3 led to decreased global 5 hmC levels and mild DNA hypermethylation of above-mentioned regions. The methylation changes observed in the overexpression and knockdown studies were mostly non-reciprocal and occurred with different preference depending on endogenous methylation and gene expression levels. Single-nucleotide 5 hmC profiling performed on a genome-wide scale revealed that TET1 overexpression induced 5 mC oxidation without a distribution bias among genetic elements and structures. Detailed analysis showed that this oxidation was related to endogenous 5 hmC levels. In addition, our results support the notion that the effects of TET1 overexpression on gene expression are generally unrelated to its catalytic activity.