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.

  • Brief Communication
  • Published:

Single-cell multimodal profiling reveals cellular epigenetic heterogeneity

Nature Methods volume 13, pages 833–836 (2016)Cite this article

Subjects

Abstract

Sample heterogeneity often masks DNA methylation signatures in subpopulations of cells. Here, we present a method to genotype single cells while simultaneously interrogating gene expression and DNA methylation at multiple loci. We used this targeted multimodal approach, implemented on an automated, high-throughput microfluidic platform, to assess primary lung adenocarcinomas and human fibroblasts undergoing reprogramming by profiling epigenetic variation among cell types identified through genotyping and transcriptional analysis.

This is a preview of subscription content, access via your institution

Access options

Buy this article

  • Purchase on SpringerLink
  • Instant access to full article PDF
Buy now

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

Figure 1: sc-GEM analysis of human cellular reprogramming.
Figure 2: Transcriptional and DNA methylation dynamics during human cellular reprogramming.
Figure 3: sc-GEM profiling of lung adenocarcinoma samples reveals population heterogeneity.

Similar content being viewed by others

Accession codes

Primary accessions

Sequence Read Archive

Referenced accessions

Gene Expression Omnibus

References

  1. Navin, N. et al. Nature 472, 90–94 (2011).

    Article  CAS  Google Scholar 

  2. Zong, C., Lu, S., Chapman, A.R. & Xie, X.S. Science 338, 1622–1626 (2012).

    Article  CAS  Google Scholar 

  3. Tang, F. et al. Cell Stem Cell 6, 468–478 (2010).

    Article  CAS  Google Scholar 

  4. Patel, A.P. et al. Science 344, 1396–1401 (2014).

    Article  CAS  Google Scholar 

  5. Guo, H. et al. Genome Res. 23, 2126–2135 (2013).

    Article  CAS  Google Scholar 

  6. Smallwood, S.A. et al. Nat. Methods 11, 817–820 (2014).

    Article  CAS  Google Scholar 

  7. Angermueller, C. et al. Nat. Methods 13, 229–232 (2016).

    Article  CAS  Google Scholar 

  8. Hou, Y. et al. Cell Res. 26, 304–319 (2016).

    Article  CAS  Google Scholar 

  9. Hu, Y. et al. Genome Biol. 17, 88 (2016).

    Article  Google Scholar 

  10. Macaulay, I.C. et al. Nat. Methods 12, 519–522 (2015).

    Article  CAS  Google Scholar 

  11. Bock, C., Farlik, M. & Sheffield, N.C. Trends Biotechnol. 34, 605–608 (2016).

    Article  CAS  Google Scholar 

  12. Lorthongpanich, C. et al. Science 341, 1110–1112 (2013).

    Article  CAS  Google Scholar 

  13. Cheow, L.F., Quake, S.R., Burkholder, W.F. & Messerschmidt, D.M. Nat. Protoc. 10, 619–631 (2015).

    Article  CAS  Google Scholar 

  14. Kim, K. et al. Nature 467, 285–290 (2010).

    Article  CAS  Google Scholar 

  15. Nishino, K. et al. PLoS Genet. 7, e1002085 (2011).

    Article  CAS  Google Scholar 

  16. Selamat, S.A. et al. Genome Res. 22, 1197–1211 (2012).

    Article  CAS  Google Scholar 

  17. Heller, G. et al. Carcinogenesis 34, 513–521 (2013).

    Article  CAS  Google Scholar 

  18. Kerr, K.M., Galler, J.S., Hagen, J.A., Laird, P.W. & Laird-Offringa, I.A. Dis. Markers 23, 5–30 (2007).

    Article  CAS  Google Scholar 

  19. Masters, G.A. et al. J. Clin. Oncol. 33, 3488–3515 (2015).

    Article  CAS  Google Scholar 

  20. Rauch, T. et al. Proc. Natl. Acad. Sci. USA 104, 5527–5532 (2007).

    Article  CAS  Google Scholar 

  21. Onder, T.T. et al. Nature 483, 598–602 (2012).

    Article  CAS  Google Scholar 

  22. Sanchez-Freire, V., Ebert, A.D., Kalisky, T., Quake, S.R. & Wu, J.C. Nat. Protoc. 7, 829–838 (2012).

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We thank T.H. Charn and S. Samydurai for technical help with developing Fluidigm C1 Scripts, C. El Farran for bioinformatics help, H.F. Wang for help with cell culture, and the GIS-Fluidigm Single-Cell Omics Center and the IMCB DNA Sequencing Facility for core facilities support. This work was supported by Career Development grant 14302FG095 (L.F.C. and E.T.C.) and an A*STAR Investigatorship grant (Y.-H.L.) from the Joint Council Office of the Agency of Science, Technology and Research. S.R.Q. is an investigator of the Howard Hughes Medical Institute.

Author information

Author notes

  1. Elise T Courtois, Yuliana Tan, Ramya Viswanathan and Qiaorui Xing: These authors contributed equally to this work.

Authors and Affiliations

  1. Institute of Molecular and Cell Biology, Singapore

    Lih Feng Cheow, Ramya Viswanathan, Qiaorui Xing, Yuin-Han Loh & William F Burkholder

  2. Genome Institute of Singapore, Singapore

    Elise T Courtois, Yuliana Tan, Daniel S W Tan & Paul Robson

  3. School of Biological Sciences, Nanyang Technological University, Singapore

    Qiaorui Xing

  4. Bioinformatics Institute, Singapore

    Rui Zhen Tan

  5. National Cancer Centre Singapore, Singapore

    Daniel S W Tan

  6. Department of Biological Sciences, National University of Singapore, Singapore

    Paul Robson & Yuin-Han Loh

  7. The Jackson Laboratory for Genomic Medicine, Farmington, Connecticut, USA

    Paul Robson

  8. Departments of Bioengineering and Applied Physics, Stanford University, Stanford, California, USA

    Stephen R Quake

  9. Howard Hughes Medical Institute, Stanford, California, USA

    Stephen R Quake

Authors
  1. Lih Feng Cheow
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Elise T Courtois
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yuliana Tan
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Ramya Viswanathan
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Qiaorui Xing
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Rui Zhen Tan
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Daniel S W Tan
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Paul Robson
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Yuin-Han Loh
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Stephen R Quake
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. William F Burkholder
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

L.F.C. conceived the study. E.T.C., Y.-H.L., S.R.Q. and W.F.B. assisted in experimental design. L.F.C., E.T.C., Y.T., R.V., Q.X. and R.Z.T. performed the experiments and analyzed the data. D.S.W.T., P.R., Y.-H.L., S.R.Q. and W.F.B. supervised the study. L.F.C. wrote the manuscript. E.T.C., Y.T., R.V. and Q.X. contributed equally to the study. All authors discussed the results and commented on the manuscript.

Corresponding authors

Correspondence to Lih Feng Cheow, Stephen R Quake or William F Burkholder.

Ethics declarations

Competing interests

The authors declare no competing financial interests.

Integrated supplementary information

Supplementary Figure 1 Workflow for multimodal analysis using sc-GEM on a programmable microfluidic device (Fluidigm).

1) Cells are captured and lysed, 2) RNA from lysed single cell is reverse-transcribed to obtain complementary DNA (cDNA), 3) genomic DNA (gDNA) is digested with a methylation-sensitive restriction enzyme, 4) multiplexed PCR using primers specific for target cDNA (for gene expression assays) and target gDNA (for DNA methylation and genotype assays) is performed on-chip. The products from the chip are then collected, and a portion of each sample is reamplified with genotyping primers and tagged with sample-specific barcoded sequencing adapters. Multiple amplicons are pooled together and sequenced on a next generation sequencing platform to obtain the targeted genotype of the corresponding single cells. Another portion of each product collected from the chip is amplified with nested primers in a microfluidics qPCR chip (the Fluidigm Dynamic Array IFC) to obtain gene expression data and the DNA methylation status of the target loci. 5mC: Methylated CpG sites, blue strands: cellular RNA, red strands: cellular DNA, green strands: complementary DNA, orange strands: primers.

Supplementary Figure 2 Quality metrics of sc-GEM data obtained from unreprogrammed and reprogrammed human fibroblast.

(a-f) Sensitivity of sc-GEM gene expression detection was measured on bulk RNA samples, using a 1:1 mixture of RNA from fibroblast (BJ) and human ES cells. Titration curves on 10-fold dilution of bulk RNA showed that sc-GEM is able to amplify multiple target mRNA over a wide dynamic range. A parallel analysis on single cells undergoing reprogramming indicated that the average amount of target mRNA per cell is within sc-GEM’s linear range of analysis across multiple genes, including somatic genes (a, TFAP2A; b, LUM), pluripotent genes (c, NANOG; d, SOX2; e, ZFP42) and housekeeping gene (f, ACTB), as visualized by comparing the distribution of qPCR threshold cycle number values obtained on single cells to the dynamic range of sc-GEM on bulk RNA standards (horizontal dotted line: mean Ct value +/- 1 standard deviation in single cells [35 fibroblast cells in A-B and 58 human iPS and ES cells in c-f]). (g) Sensitivity of sc-GEM single copy DNA detection was measured at 22 loci that did not contain a CCGG site in 90 single cells. 1972 out of the 1980 single cell PCR reactions produced a specific product, giving a very high estimated successful amplification rate of 99.6%. None of these loci amplified in three microfluidic chambers in the same chip that did not capture any cells, indicating that human contamination is not a major issue in sc-GEM.

Supplementary Figure 3 Correlation between methylation state of a single CpG site and the surrounding region of interest.

Bisulfite sequencing data from ENCODE showed that methylation state of a single CCGG site (dotted line indicate the CCGG site assayed in sc-GEM) is a good proxy for the regional methylation state in representative regions that undergo (a) demethylation and (b) de novo methylation during human cellular reprogramming. Red and green colors indicate methylated and unmethylated CpG sites respectively.

Supplementary Figure 4 Correlation between single cell and bulk gene expression bulk populations.

Single cell ensemble averaged gene expression results during human cellular reprogramming correlate well with bulk gene expression assays results (r = 0.79).

Supplementary Figure 5 Hierarchical clustering of single cell gene expression during human cellular reprogramming.

Hierarchical clustering separated the cells into two distinct clusters. TRA-1-60+ cells in cluster 1 expressed somatic markers (e.g. TFAP2A, TBX3, FOSL1) and early pluripotency markers (e.g. NANOG, SALL4, DNMT3B) whereas cells in cluster 2 expressed early pluripotency markers and late pluripotency markers (e.g. ZIC3, ZFP42, TFCP2L1). TRA-1-60+ cells in cluster #1 and #2 are classified as early pluripotent cells and late pluripotent cells respectively.

Supplementary Figure 6 Comparison of DNA methylation results obtained using SCRAM and bisulfite sequencing.

Single cell ensemble averaged DNA methylation results during human cellular reprogramming is comparable to results obtained using reduced representation bisulfite sequencing from ENCODE (GSE27584).

Supplementary Figure 7 Cell-to-cell correlation matrix of single cell gene expression and single cell DNA methylation.

The data shows the presence of subpopulations of late pluripotent cells within TRA-1-60+ populations at days 16 and 24 and the initiation of DNA methylation changes when cells expressed late pluripotent markers.

Supplementary Figure 8 Boxplot of SALL4 gene expression in d16 and d24 single cells.

Cells were grouped according to the methylation state of the SALL4 promoter measured by sc-GEM. SALL4 gene expression was significantly higher in cells with unmethylated promoter (p = 3.32 x 10-11, Welch Two Sample One-tailed t-test). Rectangle in boxplot spans the first to third quartile. Bold line inside the rectangle shows the median and whiskers above and below the box indicates the maximum and minimum, n is the number of cells in each group.

Supplementary Figure 9 Workflow for tumor sample processing in preparation for sc-GEM.

Lung tumor or non tumor lung samples were surgically removed, dissociated into single cells and stored at -80 °C. Samples were thawed on the day of experiment and stained with propidium iodide. Propidium iodide negative cells were FACS sorted directly onto the Fluidigm C1 chips for sc-GEM experiments.

Supplementary Figure 10 Single cell combined genotype and DNA methylation profiling reveals tumor-specific patterns of DNA methylation.

Hierarchical clustering of single cell DNA methylation profiles from three tumor samples (LUAD1, LUAD2, LUAD3) and a non-tumor lung sample (NTL2) yielded two distinct clusters. The origins of single cells are indicated on the top color bars and their genotypes are indicated on the bottom color bars. Mutation in Exon 21 is a substitution (L858R) while mutation in Exon 19 is a deletion (L747-T751 del).

Supplementary Figure 11 t-SNE analysis of single cell DNA methylation data from tumor and normal tissues.

Samples groups are colored according to a) tissue of origin or b) mutation status.

Supplementary Figure 12 Silhouette analysis of LUAD 1 hierarchical clustering based on gene expression.

Cells are optimally separated into two clusters by silhouette analysis.

Supplementary Figure 13 Comparison between mutation calling with next generation sequencing and Sanger sequencing.

Data is based on one PCR amplicon (EGFR Exon 21) in 60 cells from tumor sample of one patient (A165). A mutant call is made if the mutant read ratio with next generation sequencing exceeds 0.05. High concordance (95%) in mutation calling is observed with these two methods.

Supplementary information

Supplementary Text and Figures

Supplementary Figures 1–13 and Supplementary Table 2 (PDF 1922 kb)

Primer List

Primer sequences for sc-GEM assay (XLSX 27 kb)

Source data

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cheow, L., Courtois, E., Tan, Y. et al. Single-cell multimodal profiling reveals cellular epigenetic heterogeneity. Nat Methods 13, 833–836 (2016). https://doi.org/10.1038/nmeth.3961

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/nmeth.3961

This article is cited by

Search

Advanced search

Quick links

Nature Briefing: Translational Research

Sign up for the Nature Briefing: Translational Research newsletter — top stories in biotechnology, drug discovery and pharma.

Get what matters in translational research, free to your inbox weekly. Sign up for Nature Briefing: Translational Research