Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Protocol
  • Published:

Profiling DNA methylome landscapes of mammalian cells with single-cell reduced-representation bisulfite sequencing

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

Rent or buy this article

Prices vary by article type

from$1.95

to$39.95

Prices may be subject to local taxes which are calculated during checkout

Figure 1: Flowchart of the experimental procedures of the scRRBS technique.
Figure 2: The size distribution of the typical scRRBS libraries.
Figure 3: Schematic of the bioinformatic analysis procedures of the scRRBS data.
Figure 4: Box plots of the mapping efficiencies and the total unique CpG sites covered in our scRRBS data.
Figure 5: The methylation status of a representative locus of sperm-specific differentially methylated regions (DMRs).
Figure 6: The average DNA methylation levels across gene bodies and the flanking intergenic regions.

Similar content being viewed by others

References

  1. Bird, A. DNA methylation patterns and epigenetic memory. Genes Dev. 16, 6–21 (2002).

    Article  CAS  Google Scholar 

  2. Smith, Z.D. et al. A unique regulatory phase of DNA methylation in the early mammalian embryo. Nature 484, 339–344 (2012).

    Article  CAS  Google Scholar 

  3. Smallwood, S.A. et al. Dynamic CpG island methylation landscape in oocytes and preimplantation embryos. Nat. Genet. 43, 811–814 (2011).

    Article  CAS  Google Scholar 

  4. Gu, T.P. et al. The role of Tet3 DNA dioxygenase in epigenetic reprogramming by oocytes. Nature 477, 606–610 (2011).

    Article  CAS  Google Scholar 

  5. Hackett, J.A. et al. Germline DNA demethylation dynamics and imprint erasure through 5-hydroxymethylcytosine. Science 339, 448–452 (2013).

    Article  CAS  Google Scholar 

  6. 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).

    Article  CAS  Google Scholar 

  7. Mayer, W., Niveleau, A., Walter, J., Fundele, R. & Haaf, T. Demethylation of the zygotic paternal genome. Nature 403, 501–502 (2000).

    Article  CAS  Google Scholar 

  8. Wossidlo, M. et al. 5-Hydroxymethylcytosine in the mammalian zygote is linked with epigenetic reprogramming. Nat. Commun. 2, 241 (2011).

    Article  Google Scholar 

  9. 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).

    Article  CAS  Google Scholar 

  10. Guo, H. et al. The DNA methylation landscape of human early embryos. Nature 511, 606–610 (2014).

    Article  CAS  Google Scholar 

  11. Smith, Z.D. et al. DNA methylation dynamics of the human preimplantation embryo. Nature 511, 611–615 (2014).

    Article  CAS  Google Scholar 

  12. 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).

    Article  Google Scholar 

  13. Seisenberger, S. et al. The dynamics of genome-wide DNA methylation reprogramming in mouse primordial germ cells. Mol. Cell 48, 849–862 (2012).

    Article  CAS  Google Scholar 

  14. Xie, W. et al. Epigenomic analysis of multilineage differentiation of human embryonic stem cells. Cell 153, 1134–1148 (2013).

    Article  CAS  Google Scholar 

  15. Reik, W. & Kelsey, G. Epigenetics: cellular memory erased in human embryos. Nature 511, 540–541 (2014).

    Article  CAS  Google Scholar 

  16. Burgess, D.J. Stem cells: epigenome reprogramming—of mice and men. Nat. Rev. Genet. 15, 571 (2014).

    Article  Google Scholar 

  17. 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).

    Article  CAS  Google Scholar 

  18. Lister, R. et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature 462, 315–322 (2009).

    Article  CAS  Google Scholar 

  19. 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).

    Article  CAS  Google Scholar 

  20. Molaro, A. et al. Sperm methylation profiles reveal features of epigenetic inheritance and evolution in primates. Cell 146, 1029–1041 (2011).

    Article  CAS  Google Scholar 

  21. Gifford, C.A. et al. Transcriptional and epigenetic dynamics during specification of human embryonic stem cells. Cell 153, 1149–1163 (2013).

    Article  CAS  Google Scholar 

  22. Ziller, M.J. et al. Charting a dynamic DNA methylation landscape of the human genome. Nature 500, 477–481 (2013).

    Article  CAS  Google Scholar 

  23. Lister, R. et al. Global epigenomic reconfiguration during mammalian brain development. Science 341, 1237905 (2013).

    Article  Google Scholar 

  24. 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).

    Article  Google Scholar 

  25. Meissner, A. et al. Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature 454, 766–770 (2008).

    Article  CAS  Google Scholar 

  26. Meissner, A. et al. Reduced representation bisulfite sequencing for comparative high-resolution DNA methylation analysis. Nucleic Acids Res. 33, 5868–5877 (2005).

    Article  CAS  Google Scholar 

  27. Gu, H. et al. Preparation of reduced representation bisulfite sequencing libraries for genome-scale DNA methylation profiling. Nat. Protoc. 6, 468–481 (2011).

    Article  CAS  Google Scholar 

  28. 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).

    Article  CAS  Google Scholar 

  29. Smith, Z.D., Gu, H., Bock, C., Gnirke, A. & Meissner, A. High-throughput bisulfite sequencing in mammalian genomes. Methods 48, 226–232 (2009).

    Article  CAS  Google Scholar 

  30. Gu, H. et al. Genome-scale DNA methylation mapping of clinical samples at single-nucleotide resolution. Nat. Methods 7, 133–136 (2010).

    Article  CAS  Google Scholar 

  31. Boyle, P. et al. Gel-free multiplexed reduced representation bisulfite sequencing for large-scale DNA methylation profiling. Genome Biol. 13, R92 (2012).

    Article  CAS  Google Scholar 

  32. Smallwood, S.A. et al. Single-cell genome-wide bisulfite sequencing for assessing epigenetic heterogeneity. Nat. Methods 11, 817–820 (2014).

    Article  CAS  Google Scholar 

  33. 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).

    Article  CAS  Google Scholar 

  34. 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).

    Article  CAS  Google Scholar 

  35. Krueger, F. & Andrews, S.R. Bismark: a flexible aligner and methylation caller for bisulfite-seq applications. Bioinformatics 27, 1571–1572 (2011).

    Article  CAS  Google Scholar 

  36. Xi, Y. & Li, W. BSMAP: whole-genome bisulfite sequence mapping program. BMC Bioinformatics 10, 232 (2009).

    Article  Google Scholar 

  37. Chen, P.Y., Cokus, S.J. & Pellegrini, M. BS Seeker: precise mapping for bisulfite sequencing. BMC Bioinformatics 11, 203 (2010).

    Article  CAS  Google Scholar 

  38. 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).

    Article  Google Scholar 

  39. Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).

    Article  Google Scholar 

  40. Tang, F. et al. RNA-seq analysis to capture the transcriptome landscape of a single cell. Nat. Protoc. 5, 516–535 (2010).

    Article  CAS  Google Scholar 

  41. Tang, F. et al. mRNA-seq whole-transcriptome analysis of a single cell. Nat. Methods 6, 377–382 (2009).

    Article  CAS  Google Scholar 

  42. 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).

    Article  CAS  Google Scholar 

  43. Evrony, G.D. et al. Single-neuron sequencing analysis of L1 retrotransposition and somatic mutation in the human brain. Cell 151, 483–496 (2012).

    Article  CAS  Google Scholar 

  44. Frumkin, D. et al. Amplification of multiple genomic loci from single cells isolated by laser micro-dissection of tissues. BMC Biotechnol. 8, 17 (2008).

    Article  Google Scholar 

  45. 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).

    Article  CAS  Google Scholar 

  46. Fan, H.C., Wang, J., Potanina, A. & Quake, S.R. Whole-genome molecular haplotyping of single cells. Nat. Biotechnol. 29, 51–57 (2011).

    Article  CAS  Google Scholar 

Download references

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

Authors

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

Correspondence to Lu Wen or Fuchou Tang.

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/nprot.2015.039

This article is cited by

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing