Abstract
The heterogeneity of DNA methylation within a population of cells necessitates DNA methylome profiling at single-cell resolution. Recently, we developed a single-cell reduced-representation bisulfite sequencing (scRRBS) technique in which we modified the original RRBS method by integrating all the experimental steps before PCR amplification into a single-tube reaction. These modifications enable scRRBS to provide digitized methylation information on ∼1 million CpG sites within an individual diploid mouse or human cell at single-base resolution. Compared with the single-cell bisulfite sequencing (scBS) technique, scRRBS covers fewer CpG sites, but it provides better coverage for CpG islands (CGIs), which are likely to be the most informative elements for DNA methylation. The entire procedure takes ∼3 weeks, and it requires strong molecular biology skills.
This is a preview of subscription content, access via your institution
Access options
Subscribe to this journal
Receive 12 print issues and online access
£169.00 per year
only £14.08 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
References
Bird, A. DNA methylation patterns and epigenetic memory. Genes Dev. 16, 6–21 (2002).
Smith, Z.D. et al. A unique regulatory phase of DNA methylation in the early mammalian embryo. Nature 484, 339–344 (2012).
Smallwood, S.A. et al. Dynamic CpG island methylation landscape in oocytes and preimplantation embryos. Nat. Genet. 43, 811–814 (2011).
Gu, T.P. et al. The role of Tet3 DNA dioxygenase in epigenetic reprogramming by oocytes. Nature 477, 606–610 (2011).
Hackett, J.A. et al. Germline DNA demethylation dynamics and imprint erasure through 5-hydroxymethylcytosine. Science 339, 448–452 (2013).
Iqbal, K., Jin, S.G., Pfeifer, G.P. & Szabo, P.E. Reprogramming of the paternal genome upon fertilization involves genome-wide oxidation of 5-methylcytosine. Proc. Natl. Acad. Sci. USA 108, 3642–3647 (2011).
Mayer, W., Niveleau, A., Walter, J., Fundele, R. & Haaf, T. Demethylation of the zygotic paternal genome. Nature 403, 501–502 (2000).
Wossidlo, M. et al. 5-Hydroxymethylcytosine in the mammalian zygote is linked with epigenetic reprogramming. Nat. Commun. 2, 241 (2011).
Okada, Y., Yamagata, K., Hong, K., Wakayama, T. & Zhang, Y. A role for the elongator complex in zygotic paternal genome demethylation. Nature 463, 554–558 (2010).
Guo, H. et al. The DNA methylation landscape of human early embryos. Nature 511, 606–610 (2014).
Smith, Z.D. et al. DNA methylation dynamics of the human preimplantation embryo. Nature 511, 611–615 (2014).
Hackett, J.A. & Surani, M.A. DNA methylation dynamics during the mammalian life cycle. Philos. Trans. R. Soc. Lond. B Biol. Sci. 368, 20110328 (2013).
Seisenberger, S. et al. The dynamics of genome-wide DNA methylation reprogramming in mouse primordial germ cells. Mol. Cell 48, 849–862 (2012).
Xie, W. et al. Epigenomic analysis of multilineage differentiation of human embryonic stem cells. Cell 153, 1134–1148 (2013).
Reik, W. & Kelsey, G. Epigenetics: cellular memory erased in human embryos. Nature 511, 540–541 (2014).
Burgess, D.J. Stem cells: epigenome reprogramming—of mice and men. Nat. Rev. Genet. 15, 571 (2014).
Guo, H. et al. Single-cell methylome landscapes of mouse embryonic stem cells and early embryos analyzed using reduced representation bisulfite sequencing. Genome Res. 23, 2126–2135 (2013).
Lister, R. et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature 462, 315–322 (2009).
Shirane, K. et al. Mouse oocyte methylomes at base resolution reveal genome-wide accumulation of non-CpG methylation and role of DNA methyltransferases. PLoS Genet. 9, e1003439 (2013).
Molaro, A. et al. Sperm methylation profiles reveal features of epigenetic inheritance and evolution in primates. Cell 146, 1029–1041 (2011).
Gifford, C.A. et al. Transcriptional and epigenetic dynamics during specification of human embryonic stem cells. Cell 153, 1149–1163 (2013).
Ziller, M.J. et al. Charting a dynamic DNA methylation landscape of the human genome. Nature 500, 477–481 (2013).
Lister, R. et al. Global epigenomic reconfiguration during mammalian brain development. Science 341, 1237905 (2013).
Wen, L. et al. Whole-genome analysis of 5-hydroxymethylcytosine and 5-methylcytosine at base resolution in the human brain. Genome Biol. 15, R49 (2014).
Meissner, A. et al. Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature 454, 766–770 (2008).
Meissner, A. et al. Reduced representation bisulfite sequencing for comparative high-resolution DNA methylation analysis. Nucleic Acids Res. 33, 5868–5877 (2005).
Gu, H. et al. Preparation of reduced representation bisulfite sequencing libraries for genome-scale DNA methylation profiling. Nat. Protoc. 6, 468–481 (2011).
Smallwood, S.A. & Kelsey, G. Genome-wide analysis of DNA methylation in low cell numbers by reduced representation bisulfite sequencing. Methods Mol. Biol. 925, 187–197 (2012).
Smith, Z.D., Gu, H., Bock, C., Gnirke, A. & Meissner, A. High-throughput bisulfite sequencing in mammalian genomes. Methods 48, 226–232 (2009).
Gu, H. et al. Genome-scale DNA methylation mapping of clinical samples at single-nucleotide resolution. Nat. Methods 7, 133–136 (2010).
Boyle, P. et al. Gel-free multiplexed reduced representation bisulfite sequencing for large-scale DNA methylation profiling. Genome Biol. 13, R92 (2012).
Smallwood, S.A. et al. Single-cell genome-wide bisulfite sequencing for assessing epigenetic heterogeneity. Nat. Methods 11, 817–820 (2014).
Miura, F., Enomoto, Y., Dairiki, R. & Ito, T. Amplification-free whole-genome bisulfite sequencing by post-bisulfite adaptor tagging. Nucleic Acids Res. 40, e136 (2012).
Guo, F. et al. Active and passive demethylation of male and female pronuclear DNA in the mammalian zygote. Cell Stem Cell 15, 447–458 (2014).
Krueger, F. & Andrews, S.R. Bismark: a flexible aligner and methylation caller for bisulfite-seq applications. Bioinformatics 27, 1571–1572 (2011).
Xi, Y. & Li, W. BSMAP: whole-genome bisulfite sequence mapping program. BMC Bioinformatics 10, 232 (2009).
Chen, P.Y., Cokus, S.J. & Pellegrini, M. BS Seeker: precise mapping for bisulfite sequencing. BMC Bioinformatics 11, 203 (2010).
Langmead, B., Trapnell, C., Pop, M. & Salzberg, S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25 (2009).
Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).
Tang, F. et al. RNA-seq analysis to capture the transcriptome landscape of a single cell. Nat. Protoc. 5, 516–535 (2010).
Tang, F. et al. mRNA-seq whole-transcriptome analysis of a single cell. Nat. Methods 6, 377–382 (2009).
Yan, L. et al. Single-cell RNA-seq profiling of human preimplantation embryos and embryonic stem cells. Nat. Struct. Mol. Biol. 20, 1131–1139 (2013).
Evrony, G.D. et al. Single-neuron sequencing analysis of L1 retrotransposition and somatic mutation in the human brain. Cell 151, 483–496 (2012).
Frumkin, D. et al. Amplification of multiple genomic loci from single cells isolated by laser micro-dissection of tissues. BMC Biotechnol. 8, 17 (2008).
Ma, C. et al. A clinical microchip for evaluation of single immune cells reveals high functional heterogeneity in phenotypically similar T cells. Nat. Med. 17, 738–743 (2011).
Fan, H.C., Wang, J., Potanina, A. & Quake, S.R. Whole-genome molecular haplotyping of single cells. Nat. Biotechnol. 29, 51–57 (2011).
Acknowledgements
We thank J. Qiao and L. Yan for their great help. The project was supported by the National Science Foundation of China (31322037 and 31271543) and the National Basic Research Program of China (2012CB966704 and 2011CB966303). This work is supported by a collaborative grant from the Center for Molecular and Translational Medicine.
Author information
Authors and Affiliations
Contributions
L.W. and F.T. conceived the experiments and supervised the project. H.G., F.G., X.L., X.W. and X.F. carried out all of the experiments. P.Z. conducted the bioinformatic analyses. H.G., P.Z., L.W. and F.T. wrote the manuscript with contributions from all of the authors.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Integrated supplementary information
Supplementary Figure 1 The number of intersecting CpG sites and the union CpG sites of all possible combinations of all eight individual mESC scRRBS samples.
The blue fitting curve in panel a represents the average number of intersecting CpG sites. The red fitting curve in panel b represents the average number of union CpG sites. Only the CpG sites with more than 3-fold coverage were taken into consideration. The union CpG sites are the combined CpG sites, only if they are covered by at least one sample. The right y-axis indicates the proportion of all CpG sites in the mouse genome.
Supplementary Figure 2 Histogram of the distribution of the CpG sites with different methylation levels showing the digitized characteristics of DNA methylation measurements in our scRRBS dataset.
(a) Histograms of the distribution of CpG sites with different methylation levels in bulk hESC and human sperm cells. (b-c) Histograms of the distribution of CpG sites with different methylation levels in single human metaphase II oocytes (panel b) and single female pronuclei (panel c).
Supplementary Figure 3 Histogram of the distribution of CpG sites with different methylation levels showing the digitized characteristics of DNA methylation measurements in our scRRBS dataset.
(a) human sperm cells; (b) human male pronuclei.
Supplementary Figure 4 A heatmap view of a section of two chromosomes indicating the demethylation process of the parental genome after fertilization.
The color key from light blue to dark blue indicates the DNA methylation level from low to high, respectively. The white regions in the left panels indicate a lack of DNA methylation information. The red (panel a) or green (panel b) bars in the right panel represent the average DNA methylation level of the corresponding genomic region. The DNA methylation levels were calculated and presented based on 30 kb windows, only if these windows have more than 5 CpG sites covered.
Supplementary Figure 5 DNA methylation dynamics after fertilization.
The average DNA methylation levels of male and female pronuclei at different time points after intra-cytoplasmic sperm injection (ICSI), which indicates the demethylation pattern in paternal and maternal genomes after fertilization, respectively. The diamond and triangle in blue represent single sperm cells and single male pronuclei at different time points after ICSI, whereas ellipses and filled circles in red represent single metaphase II oocytes and single female pronuclei at different time points after ICSI, respectively. The gray triangle represents an outlier of the male pronuclei.
Supplementary information
Supplementary Data
Archive containing custom scripts used in this protocol. The overview of custom scripts, including the script names, functions and the corresponding protocol steps in which they should be used, is listed in Table 2. (ZIP 2 kb)
Supplementary Text and Figures
Supplementary Figures 1–5 (PDF 598 kb)
Rights and permissions
About this article
Cite this article
Guo, H., Zhu, P., Guo, F. et al. Profiling DNA methylome landscapes of mammalian cells with single-cell reduced-representation bisulfite sequencing. Nat Protoc 10, 645–659 (2015). https://doi.org/10.1038/nprot.2015.039
Published:
Issue Date:
DOI: https://doi.org/10.1038/nprot.2015.039
This article is cited by
-
Single-cell DNA methylation sequencing by combinatorial indexing and enzymatic DNA methylation conversion
Cell & Bioscience (2023)
-
Smart-RRBS for single-cell methylome and transcriptome analysis
Nature Protocols (2021)
-
DNA methylation studies in cattle
Journal of Applied Genetics (2021)
-
Single-cell sequencing reveals novel mechanisms of Aflatoxin B1-induced hepatotoxicity in S phase-arrested L02 cells
Cell Biology and Toxicology (2020)
-
Photothermal and gene therapy combined with immunotherapy to gastric cancer by the gold nanoshell-based system
Journal of Nanobiotechnology (2019)