Felix Krueger

Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications.

Summary: A combination of bisulfite treatment of DNA and high-throughput sequencing (BS-Seq) can capture a snapshot of a cells epigenomic state by revealing its genome-wide cytosine methylation at single base resolution. Bismark is a flexible tool for the time-efficient analysis of BS-Seq data which performs both read mapping and methylation calling in a single convenient step. Its output discriminates between cytosines in CpG, CHG and CHH context and enables bench scientists to visualize and interpret their methylation data soon after the sequencing run is completed. Availability and implementation: Bismark is released under the GNU GPLv3+ licence. The source code is freely available from www.bioinformatics.bbsrc.ac.uk/projects/bismark/. Contact: felix.krueger@bbsrc.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

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.

The Dynamics of Genome-wide DNA Methylation Reprogramming in Mouse Primordial Germ Cells

Summary Genome-wide DNA methylation reprogramming occurs in mouse primordial germ cells (PGCs) and preimplantation embryos, but the precise dynamics and biological outcomes are largely unknown. We have carried out whole-genome bisulfite sequencing (BS-Seq) and RNA-Seq across key stages from E6.5 epiblast to E16.5 PGCs. Global loss of methylation takes place during PGC expansion and migration with evidence for passive demethylation, but sequences that carry long-term epigenetic memory (imprints, CpG islands on the X chromosome, germline-specific genes) only become demethylated upon entry of PGCs into the gonads. The transcriptional profile of PGCs is tightly controlled despite global hypomethylation, with transient expression of the pluripotency network, suggesting that reprogramming and pluripotency are inextricably linked. Our results provide a framework for the understanding of the epigenetic ground state of pluripotency in the germline.

Dynamic CpG island methylation landscape in oocytes and preimplantation embryos

Elucidating how and to what extent CpG islands (CGIs) are methylated in germ cells is essential to understand genomic imprinting and epigenetic reprogramming. Here we present, to our knowledge, the first integrated epigenomic analysis of mammalian oocytes, identifying over a thousand CGIs methylated in mature oocytes. We show that these CGIs depend on DNMT3A and DNMT3L but are not distinct at the sequence level, including in CpG periodicity. They are preferentially located within active transcription units and are relatively depleted in H3K4me3, supporting a general transcription-dependent mechanism of methylation. Very few methylated CGIs are fully protected from post-fertilization reprogramming but, notably, the majority show incomplete demethylation in embryonic day (E) 3.5 blastocysts. Our study shows that CGI methylation in gametes is not entirely related to genomic imprinting but is a strong factor in determining methylation status in preimplantation embryos, suggesting a need to reassess mechanisms of post-fertilization demethylation.

Resetting Transcription Factor Control Circuitry toward Ground-State Pluripotency in Human

Summary Current human pluripotent stem cells lack the transcription factor circuitry that governs the ground state of mouse embryonic stem cells (ESC). Here, we report that short-term expression of two components, NANOG and KLF2, is sufficient to ignite other elements of the network and reset the human pluripotent state. Inhibition of ERK and protein kinase C sustains a transgene-independent rewired state. Reset cells self-renew continuously without ERK signaling, are phenotypically stable, and are karyotypically intact. They differentiate in vitro and form teratomas in vivo. Metabolism is reprogrammed with activation of mitochondrial respiration as in ESC. DNA methylation is dramatically reduced and transcriptome state is globally realigned across multiple cell lines. Depletion of ground-state transcription factors, TFCP2L1 or KLF4, has marginal impact on conventional human pluripotent stem cells but collapses the reset state. These findings demonstrate feasibility of installing and propagating functional control circuitry for ground-state pluripotency in human cells.

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.

DNA methylome analysis using short bisulfite sequencing data

Bisulfite conversion of genomic DNA combined with next-generation sequencing (BS-seq) is widely used to measure the methylation state of a whole genome, the methylome, at single-base resolution. However, analysis of BS-seq data still poses a considerable challenge. Here we summarize the challenges of BS-seq mapping as they apply to both base and color-space data. We also explore the effect of sequencing errors and contaminants on inferred methylation levels and recommend the most appropriate way to analyze this type of data.

Parallel single-cell sequencing links transcriptional and epigenetic heterogeneity.

We report scM&T-seq, a method for parallel single-cell genome-wide methylome and transcriptome sequencing that allows for the discovery of associations between transcriptional and epigenetic variation. Profiling of 61 mouse embryonic stem cells confirmed known links between DNA methylation and transcription. Notably, the method revealed previously unrecognized associations between heterogeneously methylated distal regulatory elements and transcription of key pluripotency genes.

Large Scale Loss of Data in Low-Diversity Illumina Sequencing Libraries Can Be Recovered by Deferred Cluster Calling

Massively parallel DNA sequencing is capable of sequencing tens of millions of DNA fragments at the same time. However, sequence bias in the initial cycles, which are used to determine the coordinates of individual clusters, causes a loss of fidelity in cluster identification on Illumina Genome Analysers. This can result in a significant reduction in the numbers of clusters that can be analysed. Such low sample diversity is an intrinsic problem of sequencing libraries that are generated by restriction enzyme digestion, such as e4C-seq or reduced-representation libraries. Similarly, this problem can also arise through the combined sequencing of barcoded, multiplexed libraries. We describe a procedure to defer the mapping of cluster coordinates until low-diversity sequences have been passed. This simple procedure can recover substantial amounts of next generation sequencing data that would otherwise be lost.