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.

  • Letter
  • Published:

The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells

Abstract

Defining the transcriptional dynamics of a temporal process such as cell differentiation is challenging owing to the high variability in gene expression between individual cells. Time-series gene expression analyses of bulk cells have difficulty distinguishing early and late phases of a transcriptional cascade or identifying rare subpopulations of cells, and single-cell proteomic methods rely on a priori knowledge of key distinguishing markers1. Here we describe Monocle, an unsupervised algorithm that increases the temporal resolution of transcriptome dynamics using single-cell RNA-Seq data collected at multiple time points. Applied to the differentiation of primary human myoblasts, Monocle revealed switch-like changes in expression of key regulatory factors, sequential waves of gene regulation, and expression of regulators that were not known to act in differentiation. We validated some of these predicted regulators in a loss-of function screen. Monocle can in principle be used to recover single-cell gene expression kinetics from a wide array of cellular processes, including differentiation, proliferation and oncogenic transformation.

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

Access options

Buy this article

Purchase on Springer Link

Instant access to full article PDF

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

Figure 1: Monocle orders single-cell RNA-Seq data of differentiating myoblasts in pseudotime.
Figure 2: Monocle orders individual cells by progress through differentiation.
Figure 3: Pseudotime ordering of cells reveals genes activated or repressed early in differentiation, along with potential upstream regulators.
Figure 4: Loss-of-function screen on selected transcription factors.

Similar content being viewed by others

Accession codes

Primary accessions

Gene Expression Omnibus

References

  1. Bendall, S.C. et al. Single-cell mass cytometry of differential immune and drug responses across a human hematopoietic continuum. Science 332, 687–696 (2011).

    Article  CAS  Google Scholar 

  2. Shalek, A.K. et al. Single-cell transcriptomics reveals bimodality inexpression and splicing in immune cells. Nature 498, 236–240 (2013).

    Article  CAS  Google Scholar 

  3. Guo, G. et al. Resolution of cell fate decisions revealed by single-cell gene expression analysis from zygote to blastocyst. Dev. Cell 18, 675–685 (2010).

    Article  CAS  Google Scholar 

  4. Hashimshony, T., Wagner, F., Sher, N. & Yanai, I. CEL-Seq: single-cell RNA-Seq by multiplexed linear amplification. Cell Reports 2, 666–673 (2012).

    Article  CAS  Google Scholar 

  5. Buganim, Y. et al. Single-cell expression analyses during cellular reprogramming reveal an early stochastic and a late hierarchic phase. Cell 150, 1209–1222 (2012).

    Article  CAS  Google Scholar 

  6. Tang, F. et al. Tracing the derivation of embryonic stem cells from the inner cell mass by single-cell RNA-Seq analysis. Cell Stem Cell 6, 468–478 (2010).

    Article  CAS  Google Scholar 

  7. Ramskold, D. et al. Full-length mRNA-seq from single-cell levels of RNA and individual circulating tumor cells. Nat. Biotechnol. 30, 777–782 (2012).

    Article  Google Scholar 

  8. Simpson, E.H. The interpretation of interaction in contingency tables. J. R. Stat. Soc. Series B Stat. Methodol. 13, 238–241 (1951).

    Google Scholar 

  9. Shapiro, E., Biezuner, T. & Linnarsson, S. Single-cell sequencing-based technologies will revolutionize whole-organism science. Nat. Rev. Genet. 14, 618–630 (2013).

    Article  CAS  Google Scholar 

  10. Magwene, P.M., Lizardi, P. & Kim, J. Reconstructing the temporal ordering of biological samples using microarray data. Bioinformatics 19, 842–850 (2003).

    Article  CAS  Google Scholar 

  11. Gupta, A. & Bar-Joseph, Z. Extracting dynamics from static cancer expression data. IEEE/ACM Trans. Comput. Biol. and Bioinf. 5, 172–182 (2008).

    Article  CAS  Google Scholar 

  12. Qiu, P., Gentles, A.J. & Plevritis, S.K. Discovering biological progression underlying microarray samples. PLOS Comput. Biol. 7, e1001123 (2011).

    Article  CAS  Google Scholar 

  13. Qiu, P. et al. Extracting a cellular hierarchy from high-dimensional cytometry data with SPADE. Nat. Biotechnol. 29, 886–891 (2011).

    Article  CAS  Google Scholar 

  14. Abmayr, S.M. & Pavlath, G.K. Myoblast fusion: lessons from flies and mice. Development 139, 641–656 (2012).

    Article  CAS  Google Scholar 

  15. Tapscott, S.J. The circuitry of a master switch: Myod and the regulation of skeletal muscle gene transcription. Development 132, 2685–2695 (2005).

    Article  CAS  Google Scholar 

  16. Tomczak, K.K. Expression profiling and identification of novel genes involved in myogenic differentiation. FASEB J. 10.1096/fj.03-0568fje (2004).

  17. Hyvärinen, A. & Oja, E. Independent component analysis: algorithms and applications. Neural Netw. 411–430 (2000).

  18. Amir, E.-A.D. et al. visNe enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia. Nat. Biotechnol. 31, 545–552 (2013).

    Article  CAS  Google Scholar 

  19. Joe, A.W.B. et al. Muscle injury activates resident fibro/adipogenic progenitors that facilitate myogenesis. Nat. Cell Biol. 12, 153–163 (2010).

    Article  CAS  Google Scholar 

  20. Blais, A. An initial blueprint for myogenic differentiation. Genes Dev. 19, 553–569 (2005).

    Article  CAS  Google Scholar 

  21. Neph, S. et al. An expansive human regulatory lexicon encoded in transcription factor footprints. Nature 489, 83–90 (2012).

    Article  CAS  Google Scholar 

  22. Ernst, J. & Kellis, M. ChromHMM: automating chromatin state discovery and characterization. Nat. Methods 9, 215–216 (2012).

    Article  CAS  Google Scholar 

  23. Lun, Y., Sawadogo, M. & Perry, M. Autoactivation of Xenopus MyoD transcription and its inhibition by USF. Cell Growth Differ. 8, 275–282 (1997).

    CAS  PubMed  Google Scholar 

  24. Sansregret, L. & Nepveu, A. The multiple roles of CUX1: insights from mouse models and cell-based assays. Gene 412, 84–94 (2008).

    Article  CAS  Google Scholar 

  25. Acosta-Alvear, D. et al. XBP1 controls diverse cell type- and condition-specific transcriptional regulatory networks. Mol. Cell 27, 53–66 (2007).

    Article  CAS  Google Scholar 

  26. De Figueiredo, L.H. & de Miranda Gomes, J. Computational morphology of curves. Vis. Comput. 11, 105–112 (1994).

    Article  Google Scholar 

  27. Giesen, J. Curve reconstruction in arbitrary dimension and the traveling salesman problem. Proc. 8th Discrete Geometry and Computational Imagery Conference(DGCI) 164–176 (1999).

  28. Giesen, J. Curve reconstruction, the traveling salesman problem, and Menger's theorem on length. Discrete Comput. Geom. 24, 577–603 (2000).

    Google Scholar 

  29. Booth, K.S. & Lueker, G.S. Testing for the consecutive ones property, interval graphs, and graph planarity using PQ-tree algorithms. J. Comput. Syst. Sci. 13, 335–379 (1976).

    Article  Google Scholar 

  30. Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578 (2012).

    Article  CAS  Google Scholar 

  31. Kim, D. et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and genefusions. Genome Biol. 14, R36 (2013).

    Article  Google Scholar 

  32. Langmead, B. & Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).

    Article  CAS  Google Scholar 

  33. Harrow, J. et al. GENCODE: producing a reference annotation for ENCODE. Genome Biol. 7 (suppl. 1), S4.1–S4.9 (2006).

    Google Scholar 

  34. Trapnell, C. et al. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat. Biotechnol. 31, 46–53 (2013).

    Article  CAS  Google Scholar 

  35. Tobin, J. Estimation of relationships for limited dependent variables. Econometrica 26, 24–36 (1958).

    Article  Google Scholar 

  36. Hastie, T. & Tibshirani, R. Generalized additive models, Stat. Sci. 1, 297–318 (1986).

    Article  Google Scholar 

  37. Yee, T. W. & Wild, C.J. Generalized Additive Models. J. R. Stat. Soc. Series B Stat. Methodol. 58, 481–493 (1996).

    Google Scholar 

  38. Moffat, J. et al. A lentiviral RNAi library for human and mouse genes applied to an arrayed viral high-content screen. Cell 124, 1283–1298 (2006).

    Article  CAS  Google Scholar 

  39. Cacchiarelli, D. et al. miR-31 modulates dystrophin expression: new implications for Duchenne muscular dystrophy therapy. EMBO Rep. 12, 136–141 (2011).

    Article  CAS  Google Scholar 

  40. Zhao, S. & Fernald, R.D. Comprehensive algorithm for quantitative real-time polymerase chain reaction. J. Comput. Biol. 12, 1047–1064 (2005).

    Article  CAS  Google Scholar 

  41. Quinlan, A.R. & Hall, I.M. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842 (2010).

    Article  CAS  Google Scholar 

  42. Wu, D. & Smyth, G.K. Camera: a competitive gene set test accounting for inter-gene correlation. Nucleic Acids Res. 40, e133 (2012).

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We are grateful to S. Bordbar, C. Zhu, A. Wagers and the Broad RNAi platform for technical assistance and to M. Soumillon for discussions. C.T. was supported by a Damon Runyon Cancer Research Foundation Fellowship. D.C. was supported by a Human Frontier Science Program Fellowship. D.C. and T.S.M. were supported by the Harvard Stem Cell Institute. This work was supported by US National Institutes of Health grants 1DP2OD00667, P01GM099117 and P50HG006193-01. This work was also supported in part by the Single Cell Genomics initiative, a collaboration between the Broad Institute and Fluidigm Inc.

Author information

Authors and Affiliations

Authors

Contributions

C.T. and D.C. conceived the strategy of ordering individual cells by developmental progress. C.T. designed and wrote Monocle and performed the computational analysis. D.C., C.T., J.G., P.P., S.L. and M.M. performed the experiments. D.C., C.T. and J.L.R. designed the study. C.T., D.C., J.G., N.J.L., K.J.L., T.S.M. and J.L.R. wrote the manuscript.

Corresponding author

Correspondence to John L Rinn.

Ethics declarations

Competing interests

K.J.L. and S.L. are employees of Fluidigm, Inc. Fluidigm manufactures an instrument used in this study.

Supplementary information

Supplementary Text and Figures

Supplementary Figures 1–13 and Supplementary Tables 1 and 2 (PDF 37106 kb)

Supplementary Source Code

Source code (ZIP 10217 kb)

Rights and permissions

Reprints and permissions

About this article

Cite this article

Trapnell, C., Cacchiarelli, D., Grimsby, J. et al. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat Biotechnol 32, 381–386 (2014). https://doi.org/10.1038/nbt.2859

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/nbt.2859

This article is cited by

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