Abstract Single-cell technologies can measure the expression of thousands of molecular features in individual cells undergoing dynamic biological processes. While examining cells along a computationally-ordered pseudotime trajectory can reveal how changes in gene or protein expression impact cell fate, identifying such dynamic features is challenging due to the inherent noise in single-cell data. Here, we present DELVE, an unsupervised feature selection method for identifying a representative subset of molecular features which robustly recapitulate cellular trajectories. In contrast to previous work, DELVE uses a bottom-up approach to mitigate the effects of confounding sources of variation, and instead models cell states from dynamic gene or protein modules based on core regulatory complexes. Using simulations, single-cell RNA sequencing, and iterative immunofluorescence imaging data in the context of cell cycle and cellular differentiation, we demonstrate how DELVE selects features that better define cell-types and cell-type transitions. DELVE is available as an open-source python package:https://github.com/jranek/delve.
more »
« less
Dissecting transition cells from single-cell transcriptome data through multiscale stochastic dynamics
Abstract Advances in single-cell technologies allow scrutinizing of heterogeneous cell states, however, detecting cell-state transitions from snap-shot single-cell transcriptome data remains challenging. To investigate cells with transient properties or mixed identities, we present MuTrans, a method based on multiscale reduction technique to identify the underlying stochastic dynamics that prescribes cell-fate transitions. By iteratively unifying transition dynamics across multiple scales, MuTrans constructs the cell-fate dynamical manifold that depicts progression of cell-state transitions, and distinguishes stable and transition cells. In addition, MuTrans quantifies the likelihood of all possible transition trajectories between cell states using coarse-grained transition path theory. Downstream analysis identifies distinct genes that mark the transient states or drive the transitions. The method is consistent with the well-established Langevin equation and transition rate theory. Applying MuTrans to datasets collected from five different single-cell experimental platforms, we show its capability and scalability to robustly unravel complex cell fate dynamics induced by transition cells in systems such as tumor EMT, iPSC differentiation and blood cell differentiation. Overall, our method bridges data-driven and model-based approaches on cell-fate transitions at single-cell resolution.
more »
« less
- PAR ID:
- 10306177
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract Rapid growth of single-cell transcriptomic data provides unprecedented opportunities for close scrutinizing of dynamical cellular processes. Through investigating epithelial-to-mesenchymal transition (EMT), we develop an integrative tool that combines unsupervised learning of single-cell transcriptomic data and multiscale mathematical modeling to analyze transitions during cell fate decision. Our approach allows identification of individual cells making transition between all cell states, and inference of genes that drive transitions. Multiscale extractions of single-cell scale outputs naturally reveal intermediate cell states (ICS) and ICS-regulated transition trajectories, producing emergent population-scale models to be explored for design principles. Testing on the newly designed single-cell gene regulatory network model and applying to twelve published single-cell EMT datasets in cancer and embryogenesis, we uncover the roles of ICS on adaptation, noise attenuation, and transition efficiency in EMT, and reveal their trade-off relations. Overall, our unsupervised learning method is applicable to general single-cell transcriptomic datasets, and our integrative approach at single-cell resolution may be adopted for other cell fate transition systems beyond EMT.more » « less
-
ABSTRACT Inspired by Waddington's illustration of an epigenetic landscape, cell-fate transitions have been envisioned as bifurcating dynamical systems, wherein exogenous signaling dynamics couple to the enormously complex signaling and transcriptional machinery of a cell to elicit qualitative transitions in its collective state. Single-cell RNA sequencing (scRNA-seq), which measures the distributions of possible transcriptional states in large populations of differentiating cells, provides an alternate view, in which development is marked by the variations of a myriad of genes. Here, we present a mathematical formalism for rigorously evaluating, from a dynamical systems perspective, whether scRNA-seq trajectories display statistical signatures consistent with bifurcations and, as a case study, pinpoint regions of multistability along the neutrophil branch of hematopoeitic differentiation. Additionally, we leverage the geometric features of linear instability to identify the low-dimensional phase plane in gene expression space within which the multistability unfolds, highlighting novel genetic players that are crucial for neutrophil differentiation. Broadly, we show that a dynamical systems treatment of scRNA-seq data provides mechanistic insights into the high-dimensional processes of cellular differentiation, taking a step toward systematic construction of mathematical models for transcriptomic dynamics.more » « less
-
Cell state transitions are often triggered by large changes in the concentrations of transcription factors and therefore large differences in their stoichiometric ratios. Whether cells can elicit transitions using modest changes in the ratios of co-expressed factors is unclear. Here we investigate how cells in the Drosophila eye resolve state transitions by quantifying the expression dynamics of the ETS transcription factors Pnt and Yan. Eye progenitor cells maintain a relatively constant ratio of Pnt/Yan protein despite expressing both proteins with pulsatile dynamics. A rapid and sustained two-fold increase in the Pnt/Yan ratio accompanies transitions to photoreceptor fates. Genetic perturbations that modestly disrupt the Pnt/Yan ratio produce fate transition defects consistent with the hypothesis that transitions are normally driven by a two-fold shift in the ratio. A biophysical model based on cooperative Yan-DNA binding coupled with non-cooperative Pnt-DNA binding illustrates how two-fold ratio changes could generate ultrasensitive changes in target gene transcription to drive fate transitions. Thus, coupling cell state transitions to the Pnt/Yan ratio sensitizes the system to modest fold-changes, conferring robustness and ultrasensitivity to the developmental program.more » « less
-
One of the fundamental questions in developmental biology is how a cell is specified to differentiate as a specialized cell type. Traditionally, plant cell types were defined based on their function, location, morphology, and lineage. Currently, in the age of single-cell biology, researchers typically attempt to assign plant cells to cell types by clustering them based on their transcriptomes. However, because cells are dynamic entities that progress through the cell cycle and respond to signals, the transcriptome also reflects the state of the cell at a particular moment in time, raising questions about how to define a cell type. We suggest that these complexities and dynamics of cell states are of interest and further consider the roles signaling, stochasticity, cell cycle, and mechanical forces play in plant cell fate specification. Once established, cell identity must also be maintained. With the wealth of single-cell data coming out, the field is poised to elucidate both the complexity and dynamics of cell states.more » « less