For the best experience on the Abcam website please upgrade to a modern browser such as Google Chrome
If you continue without changing your cookie settings, we'll assume you’re happy with this.
Throughout DNA, chemical modifications add a layer of regulation to the expression of genes encoded within the DNA sequence. The most well-studied of these chemical modifications is 5-methylcytosine (5mC), a modification most commonly recognized as a stable, repressive regulator of gene expression. The human genome consists of approximately 1% methylated cytosine making it the most abundant and widespread DNA modification (Moore et al 2012). There a several methods available to sequence 5mC throughout the genome, all of which have pros and cons, which we will discuss later in this guide. These methods include high-resolution approaches, such as whole-genome bisulfite sequencing, and antibody-dependent DNA immunoprecipitation (DIP) or MeDIP.
5mC was initially discovered to reside within CpG islands – stretches of DNA commonly found within promoter regions enriched in CpG dinucleotides. It is within these promoter regions that 5mC acts as a stable epigenetic mark repressing gene transcription. Within the mammalian genome, methylated cytosine is initially incorporated into the DNA during early development by the de novo methyltransferase enzymes DNMT3a and DNMT3b (Okano et al 1999). These methylation marks are maintained throughout the genome by an additional methyltransferase, DNMT1, which copies DNA methylation patterns to daughter strands during DNA replication (Vertino et al 1996).
Today the notion of 5mC being an entirely stable DNA modification is less concrete. Many methylated cytosines throughout the genome, particularly within gene bodies, undergo a process known as DNA demethylation – a process that ultimately results in the removal of 5mC back to an unmodified cytosine (C). DNA demethylation can occur in one of two ways: passive DNA demethylation, where methylated cytosine is diluted from the genome due to an absence of methylation maintenance enzymes. Or active DNA demethylation, which involves the oxidation of 5mC by ten-eleven translocation (TET) enzymes into oxidized derivatives of 5mC (reviewed in Wu et al 2017).
Active DNA demethylation occurs in a cycle, starting with 5mC and finishing with an unmodified C. 5mC is initially oxidized to 5-hydroxymethlcytosine (5hmC), which is further oxidized to 5-formylcytosine (5fC), and finally, this is oxidized once more to 5-carboxylcytosine (5caC). 5fC and 5caC can be removed from DNA by thymine DNA glycosylase (TDG) in combination with base excision repair (BER) to result in an unmodified C (figure 8). 5hmC, 5fC, and 5caC have been the focus of many recent epigenetic studies. More and more are being found out about these epigenetic marks, including the potential for them to have stable epigenetic roles. Many sequencing methods have been developed to distinguish these marks throughout the genome including variations on MeDIP using 5hmC, 5fC, and 5caC antibodies, and variations on bisulfite sequencing such as TET assisted bisulfite sequencing (TAB-seq). The differences between these methods will be discussed later in this guide.
Figure 1. The cycle of DNA demethylation. Active DNA demethylation occurs by thymine DNA glycosylase (TDG) coupled with base excision repair (BER) or replication-dependent dilution of 5hmC, 5fC or 5caC. Active modification–passive dilution (AM–PD). active modification–active removal (AM–AR).
It is not possible to detect 5mC using traditional DNA amplification approaches because the mark is not maintained during sample preparation and amplification. Bisulfite conversion is one of the most widely used approaches to convert DNA methylation marks into a suitable template for amplification and downstream analysis. Bisulfite conversion uses the treatment of DNA with NaOH and sodium bisulfite in a chemical reaction that converts cytosine bases into uracil (U), while methylated cytosines are protected from the conversion (figure 9).
During downstream analysis such as PCR or sequencing, unmethylated C bases that undergo deamination in the bisulfite reaction will be interpreted as thymine (T), whereas 5-mC bases will remain unchanged and still be detected as a C by the sequencing output. This allows you to determine the locations in the genome containing methylated cytosine (Frommer et al., 1992)
Figure 2. Bisulfite conversion. Treatment of DNA with bisulfite (sulphonation) leads to the deamination of cytosine residues and converts them to uracil, while 5-methylcytosine residues remain the same
Bisulfite conversion has become the basis for several variations and applications designed for high throughput applications or the investigation of broader, whole genome-scale regions.
Here are some examples of bisulfite-based methods.
Genome-wide DNA methylation analysis
Whole-genome bisulfite sequencing (WGBS; Lister et al 2009)
Reduced representation bisulfite sequencing (RRBS; Meissner et al., 2005)
Targeted DNA methylation analysis
Methylation-specific PCR (MS-PCR; Herman et al., 1996)
Pyrosequencing (Colella et al., 2003; Tost et al., 2007)
High resolution melting (HRM) analysis (Wojdacz and Dobrovic, 2007)
Bisulfite conversion: technical considerations
Incomplete conversion. Bisulfite conversion is a very powerful method because it is relatively simple to perform, and it can deliver single-base resolution of DNA methylation status. However, the method does have some drawbacks: incomplete conversion (or on occasion, over-conversion) can occur under sub-optimal reaction conditions leading to insufficient DNA denaturation, or when the DNA strands re-anneal before completion of the reaction.
Distinguishing 5hmC. DNA degradation is often a byproduct of the harsh bisulfite conversion reaction conditions, which can make working with smaller samples challenging. Insufficient desulfonation of the reaction will leave behind residues that can inhibit DNA polymerases used in PCR. Recent evidence indicates that bisulfite conversion does not distinguish between 5mC and 5hmC Bisulfite conversion therefore lowers the overall complexity of the DNA sequence. This reduction sequence complexity can complicate primer design for downstream PCR-based interrogation or introduce challenges when attempting to uniquely map sequencing reads to a reference genome.
Another method commonly used to map the location of DNA methylation marks is DIP. DIP relies heavily on having antibodies capable of recognizing the DNA modifications of interest. However, once you have this, DIP is a straightforward and effective method. It is also considerably cheaper and easier to analyze compared to WGBS sequencing, which requires the whole genome to be sequenced. DIP only requires sequencing of the small sheared DNA regions pulled down in your IP step.
DIP has been successfully carried out for the most well-characterized DNA modifications: 5mC, 5hmC, 5fC, and 5caC (Pastor et al., 2011, Shen et al., 2013). It has been used in a range of samples, including embryonic stem (ES) cells, brain tissue, and zebrafish fish embryos. The method is similar to ChIP, but your starting material is raw genomic DNA with no chromatin required. This genomic DNA will undergo shearing to approximately 150–300bp, and then this sheared DNA can undergo heat denaturation. This step is essential as the antibody will only be able to access the modifications within denatured (open) DNA.
After DNA denaturation the sheared DNA is incubated with the antibody recognizing your modification of interest, usually overnight, and then the samples undergo an IP step to pull down all the DNA bound to the antibody and washing away any unbound DNA. We recommend using magnetic beads for this type of IP step. When you carry out DIP, it is important to treat your initial genomic DNA with RNase to remove any RNA from the samples.
Figure 3. DIP methodology. Genomic DNA is sheared, and immunoprecipitation is carried out using antibodies against your DNA modification. Pulldown DNA and input samples can then be used for qPCR, microarray, or NGS.
Genome-wide DIP analysis
DIP-sequencing (Pastor et al., 2011)
Targeted DIP analysis
DIP-PCR (Pastor et al., 2011)
DIP: technical considerations
Shear your samples appropriately. Unlike WGBS, DIP is not single-base resolution. When you are shearing your DNA samples, it is important to get these DNA fragments to a good size of between 150–300bp, to try to improve the resolution of your DIP sequencing. Having larger fragments means you will inevitably pull down more DNA flanking your DNA modification of interest and not physically bound to it. This results in broad, unspecific peaks in your sequencing analysis.
Get good antibodies. Another problem with DIP is that you need to have an antibody specific for your modification of interest. You need to make sure there is minimal cross-reactivity with similar modifications, for example, if your 5fC antibody also recognizes 5hmC, this is not ideal for mapping the location of 5fC throughout the genome. The use of antibodies for this type of sequencing also has many advantages. You are only limited by the antibodies available to you. If you wanted to investigate a modification not previously characterized in DNA, eg m6A (more commonly associated with RNA), you could do so provided that you have a specific m6A antibody.
The biggest drawback of traditional bisulfite sequencing is that it is unable to distinguish the oxidized derivatives of 5mC and will profile only 5mC itself. Fortunately, there have been many variations on bisulfite sequencing and some entirely new approaches to tackling the problem of sequencing 5hmC, 5fC, and 5caC. Here we look at some of these new methods in more detail.
Tet-assisted bisulfite sequencing (TAB-seq; Yu et al., 2012)
Oxidative bisulfite sequencing (oxBS-seq; Booth et al., 2012)
hMe-Seal. (Song et al., 2011)
Selective chemical labeling with exonuclease (SCL-exo; Sérandour et al., 2016)
M.SssI methylase-assisted bisulfite sequencing (MAB-seq; Wu et al., 2014)
5fC chemically assisted bisulfite sequencing (fCAB-seq; Song et al., 2013)
5caC chemically assisted bisulfite sequencing (caCAB-seq; Lu et al., 2013)
chemical-labeling-enabled C-to-T conversion sequencing (CLEVER-seq; Zhu et al., 2017)
It is important that you choose the best method for detecting DNA modifications that suit your needs. Consider things like whether you need single-base resolution, if you need to be able to quantify the absolute levels of the modification, and how feasible the method will be to use in your model system or sample type. Below you can find a table where we have summarized these key features for some of the available methods for sequencing 5hmC, 5fC, and 5caC.
Single base resolution?
Allows absolute quantification of the modification?
5hmC mapping only
Using 5hmC specific antibodies to enrich for 5hmC.
Pastor, W. A. et al Nature 2011
5hmC is converted to 5gmC to protect it. 5mC is converted to 5caC by TET enzymes. After Bisulphite conversion 5hmC is read as C. 5mC and 5caC are read as T.
Yu, M. et al Cell, 2012
Chemical conversion of 5hmC to 5fC using KRuO4 allows the differentiation of 5mC and 5hmC at single
Booth, M. J. et al Science, 2012
Glucosylation of 5hmC with an azide-containing glucose molecule and biotin allows for 5hmC enrichment using a biotin/streptavidin pulldown.
Song, C. X. et al. Nature Biotechnology 2011
Azide-glucose glycosylation of 5hmC followed by a biotin reaction allows endonuclease activity to stall at biotin-5gmCs.
Sérandour, A. A. et al. Genome Biology 2016
5fC and 5caC mapping
Using 5fC and 5caC specific antibodies to enrich for these marks.
Shen, L. et al. Cell 2013
M.SssI treatment of DNA converts all C into 5mC. Bisulfite conversion will then cause all C, 5mC, and 5hmC to read as C. All 5fC and 5caC will read as T’s.
Wu, H., Nature Biotechnology 2014
EtONH2 protects all 5fC in the genome from oxidation after bisulfite treatment.
Song, X. et al. Cell 2013
EDC is used to catalyze the formation of amide bonds to 5caC preventing deamination of 5caC on bisulfite conversion.
Lu, X. et al. JACS 2013
Malononitrile selectively labels 5fC creating a 5fC-M adduct which is read as a T in the sequencing.
Zhu, C. et al. Cell Stem Cell 2017
Table 6: DNA modification sequencing methods
If you have access to LC-MS/MS, then this is the best way quantify the amount of a DNA modification within total genomic DNA (Le et al 2011 and Fernandez et al., 2018). Using absolute quantification methods, LC-MS/MS gives you parallel quantification of all the DNA modifications found in total DNA from any organism and cell type (Zhang et al 2012). For absolute quantification, you are only limited which isotopic standards you have available to use as a standard to measure your sample against.
Using this technique combined with DIP (DIP-MS) allows you to determine if your DNA modification antibody is binding to your modification of interest and it will also allow you to see if it binds any other non-specific modifications. If you generate LC-MS/MS data of your DIP input and pull-down samples, you should see an enrichment of your modification of interest in the pulldown sample compared to the input. You can also then check other modifications with these same data to see if anything else came out as enriched in your samples to test for non-specific antibody binding. There is software being developed now that can even help you with this type of analysis.
LC/MS-MS: technical considerations
Technically challenging. MS equipment is costly and very specialized. The machine itself will require an enormous amount of maintenance and often requires its own technician to keep on top of things. Operating the machine is complicated and requires specialized training so it may be difficult to obtain this type of MS data on your own. Consider obtaining this data through collaborations or paid services if it is not feasible for you to purchase your own LC/MS-MS equipment.
It is also possible to carry out IHC/ICC for DNA modifications. This can be done with a few simple additions to your standard IHC/ICC protocol. The most significant difference you will need to consider is that antibodies against DNA modifications cannot access and bind to the modification if it sits within double-stranded DNA. This means that you will need to denature the DNA making it single-stranded and accessible by the antibodies.
The most common of DNA denaturation is to treat your samples with acid. This is usually 4N hydrochloric acid (HCL) applied directly to you IHC/ICC slides (Yamaguchi et al., 2013 and Kaefer et al., 2016). The best time to add this step to your protocol is before the addition of the primary antibody. Once you have permeabilized your cells or tissues with a detergent (eg PBS 0.1% Triton) you can wash and add 4N HCl to denature the DNA strands. It is important to thoroughly wash the acid off once the step is complete and neutralize the acid with an alkali (eg 100 µM NaOH in PBS). After the acid is washed and neutralized you can proceed with your usual IHC/ICC steps and add the primary antibody.
When carrying out an IHC/ICC for DNA modifications you should also be wary that your antibody may recognize very similar modifications on RNA (eg 5mC on DNA and m5C on RNA). To avoid this problem, you can treat your samples with an RNase step to remove all RNA present. Again, this step should be optimized as leaving your sample in RNase for too long can also cause damage to the DNA present.
DNA modification IHC/ICC: technical considerations
Time your acid step.
Choose the right DNA stain.
5mC and its oxidized derivatives play an important role in gene silencing and promoting gene expression after DNA demethylation. It is now known that some of these DNA modifications can act as markers to recruit proteins to specific DNA sites, altering gene expression and acting as epigenetic marks. MBD3 and methyl CpG binding protein 2 (MECP2) have both been shown to bind 5hmC in addition to 5mC. Once bound to 5hmC they play a role in DNA accessibility and activation of transcription (Yildirim et al., 2011 and Mellén et al., 2012).
A common method to screen for binders of a DNA modification is to use a pull-down technique followed by MS to screen for any proteins pulled down. This method has been successfully used to find binders of 5mC, 5hmC, and 5fC (Iurlaro et al., 2013 and Sprujit et al., 2013). For this experiment, you need to create a synthetic DNA bait containing the modification you are interested in as well as baits containing other modifications and unmodified cytosine to act as controls. This DNA bait should be linked to a biotin molecule at one end that can be used to tether the bait to streptavidin-linked magnetic beads. Protein extract from your sample of interest can then be added to the tethered bait and flushed through with various wash steps to remove any non-specifically bound proteins. After this, you can elute the remaining proteins and carry out MS analysis to find out what your specific binders are.
MBDs: Technical considerations
The number of modifications.
New DNA modifications could still be out there, just not discovered yet. It has been demonstrated that some modifications traditionally considered to be RNA modifications may also be present within DNA. One good example of this is N6-adenine methylation, known as m6A within RNA and 6mA within DNA. This modification is one of the most famous and abundant RNA modifications, but now it’s known to also reside within DNA. One of the first studies to show this was from John Gurdon’s lab in 2016 (Koziol et al., 2016). They show that 6mA is within Xenopus laevis, mice, and the human genome using an antibody against 6mA to carry out DIP-seq.
Since this study, there have been several more claims that 6mA is present within DNA in zebrafish and pig genomes (Liu et al., 2016), the mouse brain following environmental stress (Yao et al., 2017), and within the Arabidopsis thaliana genome (Liang et al., 2018). One study from 2018 took this one step further and uncovered the enzymes responsible for 6mA methylation and demethylation N6AMT1 and ALKBH1 respectively (Xiao et al., 2018). The presence of enzymes actively adding and removing the DNA modification suggests that it has a real purpose to be there and potentially its own epigenetic function.
Novel DNA modifications: technical considerations
Booth, M. J., Branco, M. R., Ficz, G., Oxley, D., Krueger, F., Reik, W., & Balasubramanian, S. (2012). Quantitative sequencing of 5-methylcytosine and 5-hydroxymethylcytosine at single-base resolution. Science, 336(6083), 934–937.
Booth, M. J., Ost, T. W. B., Beraldi, D., Bell, N. M., & Branco, M. R. (2014). Europe PMC Funders Group Oxidative bisulfite sequencing of 5-methylcytosine and 5- hydroxymethylcytosine, 8(10), 1841–1851.
Colella, S., Shen, L., Baggerly, K. A., Issa, J. P. J., & Krahe, R. (2003). Sensitive and quantitative universal PyrosequencingTM methylation analysis of CpG sites. BioTechniques, 35(1), 146–150.
Fernandez, A. F., Valledor, L., Vallejo, F., Cañal, M. J., & Fraga, M. F. (2018). Quantification of Global DNA Methylation Levels by Mass Spectrometry. In J. Tost (Ed.), DNA Methylation Protocols (pp. 49–58). New York, NY: Springer New York.
ML., G., & PA., J. (1997). Rapid quantitation of methylation differences at specific sites using methylation-sensitive single nucleotide primer extension (Ms-SNuPE). Nucleic Acids Research, 25(12 PG-2529-31), 2529–2531.
Herman, J. G., Graff, J. R., Myohanen, S., Nelkin, B. D., & Baylin, S. B. (1996). Methylation-specific PCR: A novel PCR assay for methylation status of CpG islands (DNA methylation/tumor suppressor genes/pl6/p15). Proceedings of the National Academy of Sciences of the United States of America, 93(September), 9821–9826.
Kafer, G. R., Li, X., Horii, T., Suetake, I., Tajima, S., Hatada, I., & Carlton, P. M. (2016). 5-Hydroxymethylcytosine Marks Sites of DNA Damage and Promotes Genome Stability. Cell Reports, 14(6), 1283–1292.
Koziol, M. J., Bradshaw, C. R., Allen, G. E., Costa, A. S. H., Frezza, C., & Gurdon, J. B. (2016). Identification of methylated deoxyadenosines in vertebrates reveals diversity in DNA modifications. Nature Structural and Molecular Biology, 23(1), 24–30.
Le, T., Kim, K. P., Fan, G., & Faull, K. F. (2011). A sensitive mass spectrometry method for simultaneous quantification of DNA methylation and hydroxymethylation levels in biological samples. Analytical Biochemistry, 412(2), 203–209.
Liang, Z., Shen, L., Cui, X., Bao, S., Geng, Y., Yu, G., … Yu, H. (2018). DNA N6-Adenine Methylation in Arabidopsis thaliana. Developmental Cell, 45(3), 406–416.e3.
Lister, R., Pelizzola, M., Dowen, R. H., Hawkins, R. D., Hon, G., Tonti-filippini, J., … Ecker, J. R. (2009). Human DNA methylomes at base resolution show widespread epigenomic differences. Nature, 462(7271), 315–322.
Liu, J., Zhu, Y., Luo, G. Z., Wang, X., Yue, Y., Wang, X., … He, C. (2016). Abundant DNA 6mA methylation during early embryogenesis of zebrafish and pig. Nature Communications, 7(866), 1–7.
Lu, X., Song, C., Szulwach, K., Wang, Z., Weidenbacher, P., Jin, P., & He, C. (2013). Chemical Modi fi cation-Assisted Bisul fi te Sequencing ( CAB-Seq ) for 5 ‑ Carboxylcytosine Detection in DNA. J. Am. Chem. Soc., 135(25), 9315–9317.
Iurlaro, M., Ficz, G., Oxley, D., Raiber, E. A., Bachman, M., Booth, M. J., … Reik, W. (2013). A screen for hydroxymethylcytosine and formylcytosine binding proteins suggests functions in transcription and chromatin regulation. Genome Biology, 14(10), 1–11.
Meissner, A., Gnirke, A., Bell, G. W., Ramsahoye, B., Lander, E. S., & Jaenisch, R. (2005). Reduced representation bisulfite sequencing for comparative high-resolution DNA methylation analysis. Nucleic Acids Research, 33(18), 5868–5877.
Mellén, M., Ayata, P., Dewell, S., Kriaucionis, S., & Heintz, N. (2012). MeCP2 binds to 5hmC enriched within active genes and accessible chromatin in the nervous system. Cell, 151(7), 1417–1430.
Moore, L. D., Le, T., & Fan, G. (2013). DNA methylation and its basic function. Neuropsychopharmacology, 38(1), 23–38.
Okano, M., Bell, D. W., Haber, D. A., & Li, E. (1999). Cell 1999 Okano, 99, 1–11.
Pastor, W. A., Pape, U. J., Huang, Y., Henderson, H. R., Lister, R., Ko, M., … Rao, A. (2011). Genome-wide mapping of 5-hydroxymethylcytosine in embryonic stem cells. Nature, 473(7347), 394–397.
Sérandour, A. A., Avner, S., Mahé, E. A., Madigou, T., Guibert, S., Weber, M., & Salbert, G. (2016). Single-CpG resolution mapping of 5-hydroxymethylcytosine by chemical labeling and exonuclease digestion identifies evolutionarily unconserved CpGs as TET targets. Genome Biology, 17(1), 1–12.
Shen, L., Wu, H., Diep, D., Yamaguchi, S., D’Alessio, A. C., Fung, H. L., … Zhang, Y. (2013). Genome-wide analysis reveals TET- and TDG-dependent 5-methylcytosine oxidation dynamics. Cell, 153(3), 692–706.
Song, C. X., Szulwach, K. E., Fu, Y., Dai, Q., Yi, C., Li, X., … He, C. (2011). Selective chemical labeling reveals the genome-wide distribution of 5-hydroxymethylcytosine. Nature Biotechnology, 29(1), 68–75.
Song, C. X., Szulwach, K. E., Dai, Q., Fu, Y., Mao, S. Q., Lin, L., … He, C. (2013). Genome-wide profiling of 5-formylcytosine reveals its roles in epigenetic priming. Cell, 153(3), 678–691.
Spruijt, C. G., Gnerlich, F., Smits, A. H., Pfaffeneder, T., Jansen, P. W. T. C., Bauer, C., … Vermeulen, M. (2013). Dynamic readers for 5-(Hydroxy)methylcytosine and its oxidized derivatives. Cell, 152(5), 1146–1159.
Vertino, P. M., Yen, R. W., Gao, J., & Baylin, S. B. (1996). De novo methylation of CpG island sequences in human fibroblasts overexpressing DNA (cytosine-5-)-methyltransferase. Molecular and Cellular Biology, 16(8), 4555–4565.
Wu, H., Wu, X., Shen, L., & Zhang, Y. (2014). Single-base resolution analysis of active DNA demethylation using methylase-assisted bisulfite sequencing. Nature Biotechnology, 32(12), 1231–1240.
Wu, X., & Zhang, Y. (2017). TET-mediated active DNA demethylation: Mechanism, function and beyond. Nature Reviews Genetics, 18(9), 517–534.
Xiao, C. Le, Zhu, S., He, M., De Chen, D., Zhang, Q., Chen, Y., … Yan, G. R. (2018). N6-Methyladenine DNA Modification in the Human Genome. Molecular Cell, 1–13.
Yamaguchi, S., Hong, K., Liu, R., Inoue, A., Shen, L., Zhang, K., & Zhang, Y. (2013). Dynamics of 5-methylcytosine and 5-hydroxymethylcytosine during germ cell reprogramming. Cell Research, 23(3), 329–339.
Yao, B., Cheng, Y., Wang, Z., Li, Y., Chen, L., Huang, L., … Jin, P. (2017). DNA N6-methyladenine is dynamically regulated in the mouse brain following environmental stress. Nature Communications, 8(1), 1–10.
Yildirim, O., Li, R., Hung, J. H., Chen, P. B., Dong, X., Ee, L. S., … Fazzio, T. G. (2011). Mbd3/NURD complex regulates expression of 5-hydroxymethylcytosine marked genes in embryonic stem cells. Cell, 147(7), 1498–1510.
Yu, M., Hon, G. C., Szulwach, K. E., Song, C. X., Zhang, L., Kim, A., … He, C. (2012). Base-resolution analysis of 5-hydroxymethylcytosine in the mammalian genome. Cell, 149(6), 1368–1380.
Zhang, L., Zhang, L., Zhou, K., Ye, X., Zhang, J., Xie, A., … Cai, C. (2012). Simultaneous determination of global DNA methylation and hydroxymethylation levels by hydrophilic interaction liquid chromatography-tandem mass spectrometry. Journal of Biomolecular Screening, 17(7), 877–884.
Zhu, C., Gao, Y., Guo, H., Xia, B., Song, J., Wu, X., … Kee, K. (2017). Single-Cell 5-Formylcytosine Landscapes of Mammalian Early Embryos and ESCs at Single-Base Resource Single-Cell 5-Formylcytosine Landscapes of Mammalian Early Embryos and ESCs at Single-Base Resolution. Stem Cell, 1–12.