A primary role of TET proteins in establishment and maintenance of De Novo bivalency at CpG islands

Ten Eleven Translocation (TET) protein-catalyzed 5mC oxidation not only creates novel DNA modifications, such as 5hmC, but also initiates active or passive DNA demethylation. TETs’ role in the crosstalk with specific histone modifications, however, is largely elusive. Here, we show that TET2-mediated DNA demethylation plays a primary role in the de novo establishment and maintenance of H3K4me3/H3K27me3 bivalent domains underlying methylated DNA CpG islands (CGIs). Overexpression of wild type (WT), but not catalytic inactive mutant (Mut), TET2 in low-TET-expressing cells results in an increase in the level of 5hmC with accompanying DNA demethylation at a subset of CGIs. Most importantly, this alteration is sufficient in making de novo bivalent domains at these loci. Genome-wide analysis reveals that these de novo synthesized bivalent domains are largely associated with a subset of essential developmental gene promoters, which are located within CGIs and are previously silenced due to DNA methylation. On the other hand, deletion of Tet1 and Tet2 in mouse embryonic stem (ES) cells results in an apparent loss of H3K27me3 at bivalent domains, which are associated with a particular set of key developmental gene promoters. Collectively, this study demonstrates the critical role of TET proteins in regulating the crosstalk between two key epigenetic mechanisms, DNA methylation and histone methylation (H3K4me3 and H3K27me3), particularly at CGIs associated with developmental genes.     

Genome-wide investigation of the dynamic changes of epigenome modifications after global DNA methylation editing

Chromatin properties are regulated by complex networks of epigenome modifications. Currently, it is unclear how these modifications interact and if they control downstream effects such as gene expression. We employed promiscuous chromatin binding of a zinc finger fused catalytic domain of DNMT3A to introduce DNA methylation in HEK293 cells at many CpG islands (CGIs) and systematically investigated the dynamics of the introduced DNA methylation and the consequent changes of the epigenome network. We observed efficient methylation at thousands of CGIs, but it was unstable at about 90% of them, highlighting the power of genome-wide molecular processes that protect CGIs against DNA methylation. Partially stable methylation was observed at about 1000 CGIs, which showed enrichment in H3K27me3. Globally, the introduced DNA methylation strongly correlated with a decrease in gene expression indicating a direct effect. Similarly, global but transient reductions in H3K4me3 and H3K27ac were observed after DNA methylation but no changes were found for H3K9me3 and H3K36me3. Our data provide a global and time-resolved view on the network of epigenome modifications, their connections with DNA methylation and the responses triggered by artificial DNA methylation revealing a direct repressive effect of DNA methylation in CGIs on H3K4me3, histone acetylation, and gene expression.     

TET1 is a maintenance DNA demethylase that prevents methylation spreading in differentiated cells

TET1 is a 5-methylcytosine dioxygenase and its DNA demethylating activity has been implicated in pluripotency and reprogramming. However, the precise role of TET1 in DNA methylation regulation outside of developmental reprogramming is still unclear. Here, we show that overexpression of the TET1 catalytic domain but not full length TET1 (TET1-FL) induces massive global DNA demethylation in differentiated cells. Genome-wide mapping reveals that 5-hydroxymethylcytosine production by TET1-FL is inhibited as DNA methylation increases, which can be explained by the preferential binding of TET1-FL to unmethylated CpG islands (CGIs) through its CXXC domain. TET1-FL specifically accumulates 5-hydroxymethylcytosine at the edges of hypomethylated CGIs, while knockdown of endogenous TET1 induces methylation spreading from methylated edges into hypomethylated CGIs. We also found that gene expression changes after TET1-FL overexpression are relatively small and independent of its dioxygenase function. Thus, our results identify TET1 as a maintenance DNA demethylase that does not purposely decrease methylation levels, but specifically prevents aberrant methylation spreading into CGIs in differentiated cells.     

Prolonged Treatment with DNMT Inhibitors Induces Distinct Effects in Promoters and Gene-Bodies

Treatment with the demethylating drugs 5-azacytidine (AZA) and decitabine (DAC) is now recognised as an effective therapy for patients with Myelodysplastic Syndromes (MDS), a range of disorders arising in clones of hematopoietic progenitor cells. A variety of cell models have been used to study the effect of these drugs on the methylation of promoter regions of tumour suppressor genes, with recent efforts focusing on the ability of these drugs to inhibit DNA methylation at low doses. However, it is still not clear how nano-molar drug treatment exerts its effects on the methylome. In this study, we have characterised changes in DNA methylation caused by prolonged low-dose treatment in a leukemic cell model (SKM-1), and present a genome-wide analysis of the effects of AZA and DAC. At nano-molar dosages, a one-month continuous treatment halved the total number of hypermethylated probes in leukemic cells and our analysis identified 803 candidate regions with significant demethylation after treatment. Demethylated regions were enriched in promoter sequences whereas gene-body CGIs were more resistant to the demethylation process. CGI methylation in promoters was strongly correlated with gene expression but this correlation was lost after treatment. Our results indicate that CGI demethylation occurs preferentially at promoters, but that it is not generally sufficient to modify expression patterns, and emphasises the roles of other means of maintaining cell state.     

Development of a Novel Output Value for Quantitative Assessment in Methylated DNA Immunoprecipitation-CpG Island Microarray Analysis

In DNA methylation microarray analysis, quantitative assessment of intermediate methylation levels in samples with various global methylation levels is still difficult. Here, specifically for methylated DNA immunoprecipitation-CpG island (CGI) microarray analysis, we developed a new output value. The signal log ratio reflected the global methylation levels, but had only moderate linear correlation (r = 0.72) with the fraction of DNA molecules immunoprecipitated. By multiplying the signal log ratio using a coefficient obtained from the probability value that took account of signals in neighbouring probes, its linearity was markedly improved (r = 0.94). The new output value, Me value, reflected the global methylation level, had a strong correlation also with the fraction of methylated CpG sites obtained by bisulphite sequencing (r = 0.88), and had an accuracy of 71.8 and 83.8% in detecting completely methylated and unmethylated CGIs. Analysis of gastric cancer cell lines using the Me value showed that methylation of CGIs in promoters and gene bodies was associated with low and high, respectively, gene expression. The degree of demethylation of promoter CGIs after 5-aza-2''-deoxycytidine treatment had no association with that of induction of gene expression. The Me value was considered to be useful for analysis of intermediate methylation levels of CGIs.     

Methylome sequencing for fibrolamellar hepatocellular carcinoma depicts distinctive features

With the goal of studying epigenetic alterations in fibrolamellar hepatocellular carcinoma (FLC) and establish an associated DNA methylation signature, we analyzed LINE-1 methylation in a cohort of FLC and performed next-generation sequencing of DNA methylation in a training set of pure-FLCs and non-cirrhotic hepatocellular carcinomas (nc-HCC). DNA methylation was correlated with gene expression. Furthermore, we established and validated an epigenetic signature differentiating pure-FLC from other HCCs. LINE-1 methylation correlated with shorter recurrence-free survival and overall survival in resected pure-FLC patients. Unsupervised clustering using CG sites located in islands distinguished pure-FLC from nc-HCC. Major DNA methylation changes occurred outside promoters, mainly in gene bodies and intergenic regions located in the vicinity of liver developmental genes (i.e., SMARCA4 and RXRA). Partially methylated domains were more prone to DNA methylation changes. Furthermore, we identified several putative tumor suppressor genes (e.g., DLEU7) and oncogenes (e.g., DUSP4). While ∼70% of identified gene promoters gaining methylation were marked by bivalent histone marks (H3K4me3/H3K27me3) in embryonic stem cells, ∼70% of those losing methylation were marked by H3K4me3. Finally, we established a pure FLC DNA methylation signature and validated it in an independent dataset. Our analysis reveals a distinct epigenetic signature of pure FLC as compared to nc-HCC, with DNA methylation changes occurring in the vicinity of liver developmental genes. These data suggest new options for targeting FLC based on cancer epigenome aberrations.     

Effects of genome architecture and epigenetic factors on susceptibility of promoter CpG islands to aberrant DNA methylation induction

Aberrant DNA methylation is induced at specific promoter CpG islands (CGIs) in contrast with mutations. The specificity is influenced by genome architecture and epigenetic factors, but their relationship is still unknown. In this study, we isolated promoter CGIs susceptible and resistant to aberrant methylation induction during prostate and breast carcinogenesis. The effect of genome architecture was more evident for promoter CGIs susceptible in both of the two tissues than for promoter CGIs susceptible only in one tissue. Multivariate analysis of promoter CGIs with tissue-nonspecific susceptibility showed that genome architecture, namely a remote location from SINE (OR = 5.98; 95% CI = 2.33-15.34) and from LINE (OR = 2.08; 95% CI = 1.03-4.21), was associated with increased susceptibility, independent of epigenetic factors such as the presence of RNA polymerase II (OR = 0.09; 95% CI = 0.02-0.48) and H3K27me3 (OR = 3.28; 95% CI = 1.17-9.21). These results showed that methylation susceptibility of promoter CGIs is determined both by genome architecture and epigenetic factors, independently.     

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 (5mC) into 5-hydroxymethylcytosine (5hmC) 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 5hmC 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 5hmC 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 5hmC profiling performed on a genome-wide scale revealed that TET1 overexpression induced 5mC oxidation without a distribution bias among genetic elements and structures. Detailed analysis showed that this oxidation was related to endogenous 5hmC levels. In addition, our results support the notion that the effects of TET1 overexpression on gene expression are generally unrelated to its catalytic activity.