Summary: A combination of bisulfite treatment of DNA and high-throughput sequencing (BS-Seq) can capture a snapshot of a cell's 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.ukSupplementary information: Supplementary data are available at Bioinformatics online.

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.

SummaryGenome-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.

SummaryCurrent 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.

We report a single-cell bisulfite sequencing method (scBS-Seq) capable of accurately measuring DNA methylation at up to 48.4% of CpGs. We observed that ESCs grown in serum or 2i both display epigenetic heterogeneity, with “2i-like” cells present in serum cultures. In silico integration of 12 individual mouse oocyte datasets largely recapitulates the whole DNA methylome, making scBS-Seq a versatile tool to explore DNA methylation in rare cells and heterogeneous populations.

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.

Elucidating how and to what extent CpG islands (CGIs) are methylated in germ cells is essential to understand genomic imprinting and epigenetic reprogramming1-3. Here, we present 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 DNMT3L4-5, but are not distinct at the sequence level, including in CpG periodicity6. 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-fertilisation reprogramming but, surprisingly, the majority exhibits incomplete demethylation in E3.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.

We report scM&T-seq, a method for parallel single-cell genome-wide methylome and transcriptome sequencing, allowing 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 reveals novel associations between heterogeneously methylated distal regulatory elements and transcription of key pluripotency genes.