Inferring population dynamics from single-cell RNA-sequencing time series data

Recent single-cell RNA-sequencing studies have suggested that cells follow continuous transcriptomic trajectories in an asynchronous fashion during development. However, observations of cell flux along trajectories are confounded with population size effects in snapshot experiments and are therefore hard to interpret. In particular, changes in proliferation and death rates can be mistaken for cell flux. Here we present pseudodynamics, a mathematical framework that reconciles population dynamics with the concepts underlying developmental trajectories inferred from time-series single-cell data. Pseudodynamics models population distribution shifts across trajectories to quantify selection pressure, population expansion, and developmental potentials. Applying this model to time-resolved single-cell RNA-sequencing of T-cell and pancreatic beta cell maturation, we characterize proliferation and apoptosis rates and identify key developmental checkpoints, data inaccessible to existing approaches.A new computational method allows key developmental checkpoints and important parameters of population dynamics to be inferred from single-cell RNA-sequencing time series data.

[1]  A. Klochendler,et al.  G0-G1 Transition and the Restriction Point in Pancreatic β-Cells In Vivo , 2014, Diabetes.

[2]  Tommi S. Jaakkola,et al.  Learning Population-Level Diffusions with Generative RNNs , 2016, ICML.

[3]  O. Ornatsky,et al.  Mass cytometry: technique for real time single cell multitarget immunoassay based on inductively coupled plasma time-of-flight mass spectrometry. , 2009, Analytical chemistry.

[4]  Hannah A. Pliner,et al.  Reversed graph embedding resolves complex single-cell trajectories , 2017, Nature Methods.

[5]  S. Bonner-Weir,et al.  Apoptosis participates in the remodeling of the endocrine pancreas in the neonatal rat. , 1997, Endocrinology.

[6]  James A. Gagnon,et al.  Whole-organism lineage tracing by combinatorial and cumulative genome editing , 2016, Science.

[7]  Caleb Weinreb,et al.  Fundamental limits on dynamic inference from single-cell snapshots , 2017, Proceedings of the National Academy of Sciences.

[8]  Tomoharu Kajiyama,et al.  Quantitative analysis of gene expression in a single cell by qPCR , 2009, Nature Methods.

[9]  F Allgöwer,et al.  On the relationship between cell cycle analysis with ergodic principles and age-structured cell population models. , 2017, Journal of theoretical biology.

[10]  J. Junker,et al.  Simultaneous lineage tracing and cell-type identification using CRISPR/Cas9-induced genetic scars , 2018, Nature Biotechnology.

[11]  Erik Sundström,et al.  RNA velocity of single cells , 2018, Nature.

[12]  A. Cook,et al.  Proliferation and lineage potential in fetal thymic epithelial progenitor cells , 2010 .

[13]  Ludo Waltman,et al.  A smart local moving algorithm for large-scale modularity-based community detection , 2013, The European Physical Journal B.

[14]  D. M. Smith,et al.  Single-Cell Transcriptome Profiling of Human Pancreatic Islets in Health and Type 2 Diabetes , 2016, Cell metabolism.

[15]  Olle Korsgren,et al.  Identification of proliferative and mature β-cells in the islets of Langerhans , 2016, Nature.

[16]  Thomas Höfer,et al.  Disparate Individual Fates Compose Robust CD8+ T Cell Immunity , 2013, Science.

[17]  Liu Yang,et al.  Deciphering Pancreatic Islet β Cell and α Cell Maturation Pathways and Characteristic Features at the Single-Cell Level. , 2017, Cell metabolism.

[18]  P. Pereira,et al.  Cutting Edge: Thymic NK Cells Develop Independently from T Cell Precursors , 2010, The Journal of Immunology.

[19]  Fabian J Theis,et al.  Lessons Learned from Quantitative Dynamical Modeling in Systems Biology , 2013, PloS one.

[20]  Fabian J Theis,et al.  PAGA: graph abstraction reconciles clustering with trajectory inference through a topology preserving map of single cells , 2019, Genome Biology.

[21]  Fabian J Theis,et al.  Diffusion pseudotime robustly reconstructs lineage branching , 2016, Nature Methods.

[22]  Fabian J. Theis,et al.  Diffusion maps for high-dimensional single-cell analysis of differentiation data , 2015, Bioinform..

[23]  Yarden Katz,et al.  A single-cell survey of the small intestinal epithelium , 2017, Nature.

[24]  M. Kirschner,et al.  Dynamics extracted from fixed cells reveal feedback linking cell growth to cell cycle , 2013, Nature.

[25]  Fabian J. Theis,et al.  Parameter estimation for dynamical systems with discrete events and logical operations , 2016, Bioinform..

[26]  Fabian J Theis,et al.  SCANPY: large-scale single-cell gene expression data analysis , 2018, Genome Biology.

[27]  Yvan Saeys,et al.  A comparison of single-cell trajectory inference methods: towards more accurate and robust tools , 2018, bioRxiv.

[28]  H. Cho,et al.  Modelling acute myeloid leukaemia in a continuum of differentiation states , 2018, Letters in biomathematics.

[29]  R. Germain.  T-cell development and the CD4–CD8 lineage decision , 2002, Nature Reviews Immunology.

[30]  A. Hindmarsh,et al.  CVODE, a stiff/nonstiff ODE solver in C , 1996 .

[31]  E. Wolf,et al.  Postnatal Development of Numbers and Mean Sizes of Pancreatic Islets and Beta-Cells in Healthy Mice and GIPRdn Transgenic Diabetic Mice , 2011, PloS one.

[32]  Allon M. Klein,et al.  Droplet Barcoding for Single-Cell Transcriptomics Applied to Embryonic Stem Cells , 2015, Cell.

[33]  Ping Xu,et al.  A Single‐Cell Transcriptomic Atlas of Thymus Organogenesis Resolves Cell Types and Developmental Maturation , 2018, Immunity.

[34]  Cole Trapnell,et al.  The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells , 2014, Nature Biotechnology.

[35]  A. Cumano,et al.  A thymic pathway of mouse natural killer cell development characterized by expression of GATA-3 and CD127 , 2006, Nature Immunology.

[36]  A. Regev,et al.  Spatial reconstruction of single-cell gene expression , 2015, Nature Biotechnology.

[37]  Evan Z. Macosko,et al.  Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets , 2015, Cell.

[38]  S. Jacobsen,et al.  Emergence of NK-cell progenitors and functionally competent NK-cell lineage subsets in the early mouse embryo. , 2012, Blood.

[39]  P. Rigollet,et al.  Optimal-Transport Analysis of Single-Cell Gene Expression Identifies Developmental Trajectories in Reprogramming , 2019, Cell.

[40]  Anneliese O. Speak,et al.  T cell fate and clonality inference from single cell transcriptomes , 2016, Nature Methods.

[41]  Jan Hasenauer,et al.  PESTO: Parameter EStimation TOolbox , 2017, Bioinform..

[42]  Ellen V. Rothenberg,et al.  Developmental gene networks: a triathlon on the course to T cell identity , 2014, Nature Reviews Immunology.