Miguel R. Branco

Dynamic regulation of 5-hydroxymethylcytosine in mouse ES cells and during differentiation

Methylation at the 5′ position of cytosine in DNA has important roles in genome function and is dynamically reprogrammed during early embryonic and germ cell development. The mammalian genome also contains 5-hydroxymethylcytosine (5hmC), which seems to be generated by oxidation of 5-methylcytosine (5mC) by the TET family of enzymes that are highly expressed in embryonic stem (ES) cells. Here we use antibodies against 5hmC and 5mC together with high throughput sequencing to determine genome-wide patterns of methylation and hydroxymethylation in mouse wild-type and mutant ES cells and differentiating embryoid bodies. We find that 5hmC is mostly associated with euchromatin and that whereas 5mC is under-represented at gene promoters and CpG islands, 5hmC is enriched and is associated with increased transcriptional levels. Most, if not all, 5hmC in the genome depends on pre-existing 5mC and the balance between these two modifications is different between genomic regions. Knockdown of Tet1 and Tet2 causes downregulation of a group of genes that includes pluripotency-related genes (including Esrrb, Prdm14, Dppa3, Klf2, Tcl1 and Zfp42) and a concomitant increase in methylation of their promoters, together with an increased propensity of ES cells for extraembryonic lineage differentiation. Declining levels of TETs during differentiation are associated with decreased hydroxymethylation levels at the promoters of ES cell-specific genes together with increased methylation and gene silencing. We propose that the balance between hydroxymethylation and methylation in the genome is inextricably linked with the balance between pluripotency and lineage commitment.

Quantitative Sequencing of 5-Methylcytosine and 5-Hydroxymethylcytosine at Single-Base Resolution

Distinguishing Epigenetic Marks Methylation of the cytosine base in eukaryotic DNA (5mC) is an important epigenetic mark involved in gene silencing and genome stability. Methylated cytosine can be enzymatically oxidized to 5-hydroxymethylcytosine (5hmC), which may function as a distinct epigenetic mark—possibly involved in pluripotency—and it may also be an intermediate in active DNA demethylation. To be able to detect 5hmC genome-wide and at single-base resolution, Booth et al. (p. 934, published online 26 April) developed a 5hmC sequencing chemistry that selectively oxidizes 5hmC to 5-formylcytosine and then to uracil while leaving 5mC unchanged. Using this method, mouse embryonic stem cell genomic DNA was sequenced to reveal that 5hmC is found enriched at intragenic CpG islands and long interspersed nuclear element–1 retrotransposons. A sequencing method can discriminate epigenetically modified cytosine nucleotides within embryonic stem cell DNA. 5-Methylcytosine can be converted to 5-hydroxymethylcytosine (5hmC) in mammalian DNA by the ten-eleven translocation (TET) enzymes. We introduce oxidative bisulfite sequencing (oxBS-Seq), the first method for quantitative mapping of 5hmC in genomic DNA at single-nucleotide resolution. Selective chemical oxidation of 5hmC to 5-formylcytosine (5fC) enables bisulfite conversion of 5fC to uracil. We demonstrate the utility of oxBS-Seq to map and quantify 5hmC at CpG islands (CGIs) in mouse embryonic stem (ES) cells and identify 800 5hmC-containing CGIs that have on average 3.3% hydroxymethylation. High levels of 5hmC were found in CGIs associated with transcriptional regulators and in long interspersed nuclear elements, suggesting that these regions might undergo epigenetic reprogramming in ES cells. Our results open new questions on 5hmC dynamics and sequence-specific targeting by TETs.

Uncovering the role of 5-hydroxymethylcytosine in the epigenome.

Just over 2 years ago, TET1 was found to catalyse the oxidation of 5-methylcytosine, a well-known epigenetic mark, into 5-hydroxymethylcytosine in mammalian DNA. The exciting prospect of a novel epigenetic modification that may dynamically regulate DNA methylation has led to the rapid accumulation of publications from a wide array of fields, from biochemistry to stem cell biology. Although we have only started to scratch the surface, interesting clues on the role of 5-hydroxymethylcytosine are quickly emerging.

FGF Signaling Inhibition in ESCs Drives Rapid Genome-wide Demethylation to the Epigenetic Ground State of Pluripotency

Summary Genome-wide erasure of DNA methylation takes place in primordial germ cells (PGCs) and early embryos and is linked with pluripotency. Inhibition of Erk1/2 and Gsk3β signaling in mouse embryonic stem cells (ESCs) by small-molecule inhibitors (called 2i) has recently been shown to induce hypomethylation. We show by whole-genome bisulphite sequencing that 2i induces rapid and genome-wide demethylation on a scale and pattern similar to that in migratory PGCs and early embryos. Major satellites, intracisternal A particles (IAPs), and imprinted genes remain relatively resistant to erasure. Demethylation involves oxidation of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), impaired maintenance of 5mC and 5hmC, and repression of the de novo methyltransferases (Dnmt3a and Dnmt3b) and Dnmt3L. We identify a Prdm14- and Nanog-binding cis-acting regulatory region in Dnmt3b that is highly responsive to signaling. These insights provide a framework for understanding how signaling pathways regulate reprogramming to an epigenetic ground state of pluripotency.

In Vivo Control of CpG and Non-CpG DNA Methylation by DNA Methyltransferases

The enzymatic control of the setting and maintenance of symmetric and non-symmetric DNA methylation patterns in a particular genome context is not well understood. Here, we describe a comprehensive analysis of DNA methylation patterns generated by high resolution sequencing of hairpin-bisulfite amplicons of selected single copy genes and repetitive elements (LINE1, B1, IAP-LTR-retrotransposons, and major satellites). The analysis unambiguously identifies a substantial amount of regional incomplete methylation maintenance, i.e. hemimethylated CpG positions, with variant degrees among cell types. Moreover, non-CpG cytosine methylation is confined to ESCs and exclusively catalysed by Dnmt3a and Dnmt3b. This sequence position–, cell type–, and region-dependent non-CpG methylation is strongly linked to neighboring CpG methylation and requires the presence of Dnmt3L. The generation of a comprehensive data set of 146,000 CpG dyads was used to apply and develop parameter estimated hidden Markov models (HMM) to calculate the relative contribution of DNA methyltransferases (Dnmts) for de novo and maintenance DNA methylation. The comparative modelling included wild-type ESCs and mutant ESCs deficient for Dnmt1, Dnmt3a, Dnmt3b, or Dnmt3a/3b, respectively. The HMM analysis identifies a considerable de novo methylation activity for Dnmt1 at certain repetitive elements and single copy sequences. Dnmt3a and Dnmt3b contribute de novo function. However, both enzymes are also essential to maintain symmetrical CpG methylation at distinct repetitive and single copy sequences in ESCs.

The long noncoding RNA Kcnq1ot1 organises a lineage-specific nuclear domain for epigenetic gene silencing

Long noncoding RNAs are implicated in a number of regulatory functions in eukaryotic genomes. The paternally expressed long noncoding RNA (ncRNA) Kcnq1ot1 regulates epigenetic gene silencing in an imprinted gene cluster in cis over a distance of 400 kb in the mouse embryo, whereas the silenced region extends over 780 kb in the placenta. Gene silencing by the Kcnq1ot1 RNA involves repressive histone modifications, including H3K9me2 and H3K27me3, which are partly brought about by the G9a and Ezh2 histone methyltransferases. Here, we show that Kcnq1ot1 is transcribed by RNA polymerase II, is unspliced, is relatively stable and is localised in the nucleus. Analysis of conditional Dicer mutants reveals that the RNAi pathway is not involved in gene silencing in the Kcnq1ot1 cluster. Instead, using RNA/DNA FISH we show that the Kcnq1ot1 RNA establishes a nuclear domain within which the genes that are epigenetically inactivated in cis are frequently found, whereas nearby genes that are not regulated by Kcnq1ot1 are localised outside of the domain. The Kcnq1ot1 RNA domain is larger in the placenta than in the embryo, consistent with more genes in the cluster being silenced in the placenta. Our results show for the first time that autosomal long ncRNAs can establish nuclear domains, which might create a repressive environment for epigenetic silencing of adjacent genes. Long ncRNAs in imprinting clusters and the Xist RNA on the inactive X chromosome thus appear to regulate epigenetic gene silencing by similar mechanisms.

Genome-wide distribution of 5-formylcytosine in embryonic stem cells is associated with transcription and depends on thymine DNA glycosylase

BackgroundMethylation of cytosine in DNA (5mC) is an important epigenetic mark that is involved in the regulation of genome function. During early embryonic development in mammals, the methylation landscape is dynamically reprogrammed in part through active demethylation. Recent advances have identified key players involved in active demethylation pathways, including oxidation of 5mC to 5-hydroxymethylcytosine (5hmC) and 5-formylcytosine (5fC) by the TET enzymes, and excision of 5fC by the base excision repair enzyme thymine DNA glycosylase (TDG). Here, we provide the first genome-wide map of 5fC in mouse embryonic stem (ES) cells and evaluate potential roles for 5fC in differentiation.ResultsOur method exploits the unique reactivity of 5fC for pulldown and high-throughput sequencing. Genome-wide mapping revealed 5fC enrichment in CpG islands (CGIs) of promoters and exons. CGI promoters in which 5fC was relatively more enriched than 5mC or 5hmC corresponded to transcriptionally active genes. Accordingly, 5fC-rich promoters had elevated H3K4me3 levels, associated with active transcription, and were frequently bound by RNA polymerase II. TDG down-regulation led to 5fC accumulation in CGIs in ES cells, which correlates with increased methylation in these genomic regions during differentiation of ES cells in wild-type and TDG knockout contexts.ConclusionsCollectively, our data suggest that 5fC plays a role in epigenetic reprogramming within specific genomic regions, which is controlled in part by TDG-mediated excision. Notably, 5fC excision in ES cells is necessary for the correct establishment of CGI methylation patterns during differentiation and hence for appropriate patterns of gene expression during development.

Oxidative bisulfite sequencing of 5-methylcytosine and 5-hydroxymethylcytosine

To uncover the function of and interplay between the mammalian cytosine modifications 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), new techniques and advances in current technology are needed. To this end, we have developed oxidative bisulfite sequencing (oxBS-seq), which can quantitatively locate 5mC and 5hmC marks at single-base resolution in genomic DNA. In bisulfite sequencing (BS-seq), both 5mC and 5hmC are read as cytosines and thus cannot be discriminated; however, in oxBS-seq, specific oxidation of 5hmC to 5-formylcytosine (5fC) and conversion of the newly formed 5fC to uracil (under bisulfite conditions) means that 5hmC can be discriminated from 5mC. A positive readout of actual 5mC is gained from a single oxBS-seq run, and 5hmC levels are inferred by comparison with a BS-seq run. Here we describe an optimized second-generation protocol that can be completed in 2 d.

Promoter DNA methylation couples genome-defence mechanisms to epigenetic reprogramming in the mouse germline

Mouse primordial germ cells (PGCs) erase global DNA methylation (5mC) as part of the comprehensive epigenetic reprogramming that occurs during PGC development. 5mC plays an important role in maintaining stable gene silencing and repression of transposable elements (TE) but it is not clear how the extensive loss of DNA methylation impacts on gene expression and TE repression in developing PGCs. Using a novel epigenetic disruption and recovery screen and genetic analyses, we identified a core set of germline-specific genes that are dependent exclusively on promoter DNA methylation for initiation and maintenance of developmental silencing. These gene promoters appear to possess a specialised chromatin environment that does not acquire any of the repressive H3K27me3, H3K9me2, H3K9me3 or H4K20me3 histone modifications when silenced by DNA methylation. Intriguingly, this methylation-dependent subset is highly enriched in genes with roles in suppressing TE activity in germ cells. We show that the mechanism for developmental regulation of the germline genome-defence genes involves DNMT3B-dependent de novo DNA methylation. These genes are then activated by lineage-specific promoter demethylation during distinct global epigenetic reprogramming events in migratory (∼E8.5) and post-migratory (E10.5-11.5) PGCs. We propose that genes involved in genome defence are developmentally regulated primarily by promoter DNA methylation as a sensory mechanism that is coupled to the potential for TE activation during global 5mC erasure, thereby acting as a failsafe to ensure TE suppression and maintain genomic integrity in the germline.

Complex multi-enhancer contacts captured by genome architecture mapping

The organization of the genome in the nucleus and the interactions of genes with their regulatory elements are key features of transcriptional control and their disruption can cause disease. Here we report a genome-wide method, genome architecture mapping (GAM), for measuring chromatin contacts and other features of three-dimensional chromatin topology on the basis of sequencing DNA from a large collection of thin nuclear sections. We apply GAM to mouse embryonic stem cells and identify enrichment for specific interactions between active genes and enhancers across very large genomic distances using a mathematical model termed SLICE (statistical inference of co-segregation). GAM also reveals an abundance of three-way contacts across the genome, especially between regions that are highly transcribed or contain super-enhancers, providing a level of insight into genome architecture that, owing to the technical limitations of current technologies, has previously remained unattainable. Furthermore, GAM highlights a role for gene-expression-specific contacts in organizing the genome in mammalian nuclei.