Exploring intermediate cell states through the lens of single cells

Data-driven model discovery and model selection for noisy biological systems

Biological systems exhibit complex dynamics that differential equations can often adeptly represent. Ordinary differential equation models are widespread; until recently their construction has required extensive prior knowledge of the system. Machine learning methods offer alternative means of model construction: differential equation models can be learnt from data via model discovery using sparse identification of nonlinear dynamics (SINDy). However, SINDy struggles with realistic levels of biological noise and is limited in its ability to incorporate prior knowledge of the system. We propose a data-driven framework for model discovery and model selection using hybrid dynamical systems: partial models containing missing terms. Neural networks are used to approximate the unknown dynamics of a system, enabling the denoising of the data while simultaneously learning the latent dynamics. Simulations from the fitted neural network are then used to infer models using sparse regression. We show, via model selection, that model discovery using hybrid dynamical systems outperforms alternative approaches. We find it possible to infer models correctly up to high levels of biological noise of different types. We demonstrate the potential to learn models from sparse, noisy data in application to a canonical cell state transition using data derived from single-cell transcriptomics. Overall, this approach provides a practical framework for model discovery in biology in cases where data are noisy and sparse, of particular utility when the underlying biological mechanisms are partially but incompletely known.

A statistical approach for systematic identification of transition cells from scRNA-seq data

Decoding cellular state transitions is crucial for understanding complex biological processes in development and disease. While recent advancements in single-cell RNA sequencing (scRNA-seq) offer insights into cellular trajectories, existing tools primarily study expressional rather than regulatory state shifts. We present CellTran, a statistical approach utilizing paired-gene expression correlations to detect transition cells from scRNA-seq data without explicitly resolving gene regulatory networks. Applying our approach to various contexts, including tissue regeneration, embryonic development, preinvasive lesions, and humoral responses post-vaccination, reveals transition cells and their distinct gene expression profiles. Our study sheds light on the underlying molecular mechanisms driving cellular state transitions, enhancing our ability to identify therapeutic targets for disease interventions.

Morphogenesis and regeneration share a conserved core transition cell state program that controls lung epithelial cell fate

Transitional cell states are at the crossroads of crucial developmental and regenerative events, yet little is known about how these states emerge and influence outcomes. The alveolar and airway epithelia arise from distal lung multipotent progenitors, which undergo cell fate transitions to form these distinct compartments. The identification and impact of cell states in the developing lung are poorly understood. Here, we identified a population of Icam1/Nkx2-1 epithelial progenitors harboring a transitional state program remarkably conserved in humans and mice during lung morphogenesis and regeneration. Lineage-tracing and functional analyses reveal their role as progenitors to both airways and alveolar cells and the requirement of this transitional program to make distal lung progenitors competent to undergo airway cell fate specification. The identification of a common progenitor cell state in vastly distinct processes suggests a unified program reiteratively regulating outcomes in development and regeneration.

Applying single-cell and single-nucleus genomics to studies of cellular heterogeneity and cell fate transitions in the nervous system

Single-cell and single-nucleus genomic approaches can provide unbiased and multimodal insights. Here, we discuss what constitutes a molecular cell atlas and how to leverage single-cell omics data to generate hypotheses and gain insights into cell transitions in development and disease of the nervous system. We share points of reflection on what to consider during study design and implementation as well as limitations and pitfalls.

Interpreting single-cell and spatial omics data using deep neural network training dynamics

Single-cell and spatial omics datasets can be organized and interpreted by annotating single cells to distinct types, states, locations or phenotypes. However, cell annotations are inherently ambiguous, as discrete labels with subjective interpretations are assigned to heterogeneous cell populations on the basis of noisy, sparse and high-dimensional data. Here we developed Annotatability, a framework for identifying annotation mismatches and characterizing biological data structure by monitoring the dynamics and difficulty of training a deep neural network over such annotated data. Following this, we developed a signal-aware graph embedding method that enables downstream analysis of biological signals. This embedding captures cellular communities associated with target signals. Using Annotatability, we address key challenges in the interpretation of genomic data, demonstrated over eight single-cell RNA sequencing and spatial omics datasets, including identifying erroneous annotations and intermediate cell states, delineating developmental or disease trajectories, and capturing cellular heterogeneity. These results underscore the broad applicability of annotation-trainability analysis via Annotatability for unraveling cellular diversity and interpreting collective cell behaviors in health and disease.

The spatial and cellular portrait of transposable element expression during gastric cancer

Gastric Cancer (GC) is a lethal malignancy, with urgent need for the discovery of novel biomarkers for its early detection. I previously showed that Transposable Elements (TEs) become activated in early GC (EGC), suggesting a role in gene expression. Here, I follow-up on that evidence using single-cell data from gastritis to EGC, and show that TEs are expressed and follow the disease progression, with 2,430 of them being cell populations markers. Pseudotemporal trajectory modeling revealed 111 TEs associated with the origination of cancer cells. Analysis of spatial data from GC also confirms TE expression, with 204 TEs being spatially enriched in the tumor regions and the tumor microenvironment, hinting at a role of TEs in tumorigenesis. Finally, a network of TE-mediated gene regulation was modeled, indicating that ~ 2,000 genes could be modulated by TEs, with ~ 500 of them already implicated in cancer. These results suggest that TEs might play a functional role in GC progression, and highlights them as potential biomarker for its early detection.

The ability to sense the environment is heterogeneously distributed in cell populations

Channel capacity of signaling networks quantifies their fidelity in sensing extracellular inputs. Low estimates of channel capacities for several mammalian signaling networks suggest that cells can barely detect the presence/absence of environmental signals. However, given the extensive heterogeneity and temporal stability of cell state variables, we hypothesize that the sensing ability itself may depend on the state of the cells. In this work, we present an information-theoretic framework to quantify the distribution of sensing abilities from single-cell data. Using data on two mammalian pathways, we show that sensing abilities are widely distributed in the population and most cells achieve better resolution of inputs compared to an 'average cell'. We verify these predictions using live-cell imaging data on the IGFR/FoxO pathway. Importantly, we identify cell state variables that correlate with cells' sensing abilities. This information-theoretic framework will significantly improve our understanding of how cells sense in their environment.

The ability to sense the environment is heterogeneously distributed in cell populations

Channel capacity of signaling networks quantifies their fidelity in sensing extracellular inputs. Low estimates of channel capacities for several mammalian signaling networks suggest that cells can barely detect the presence/absence of environmental signals. However, given the extensive heterogeneity and temporal stability of cell state variables, we hypothesize that the sensing ability itself may depend on the state of the cells. In this work, we present an information theoretic framework to quantify the distribution of sensing abilities from single cell data. Using data on two mammalian pathways, we show that sensing abilities are widely distributed in the population and most cells achieve better resolution of inputs compared to an " average cell ". We verify these predictions using live cell imaging data on the IGFR/FoxO pathway. Importantly, we identify cell state variables that correlate with cells' sensing abilities. This information theoretic framework will significantly improve our understanding of how cells sense in their environment.

Transcriptional networks of transient cell states during human prefrontal cortex development

The human brain is divided into various anatomical regions that control and coordinate unique functions. The prefrontal cortex (PFC) is a large brain region that comprises a range of neuronal and non-neuronal cell types, sharing extensive interconnections with subcortical areas, and plays a critical role in cognition and memory. A timely appearance of distinct cell types through embryonic development is crucial for an anatomically perfect and functional brain. Direct tracing of cell fate development in the human brain is not possible, but single-cell transcriptome sequencing (scRNA-seq) datasets provide the opportunity to dissect cellular heterogeneity and its molecular regulators. Here, using scRNA-seq data of human PFC from fetal stages, we elucidate distinct transient cell states during PFC development and their underlying gene regulatory circuitry. We further identified that distinct intermediate cell states consist of specific gene regulatory modules essential to reach terminal fate using discrete developmental paths. Moreover, using in silico gene knock-out and over-expression analysis, we validated crucial gene regulatory components during the lineage specification of oligodendrocyte progenitor cells. Our study illustrates unique intermediate states and specific gene interaction networks that warrant further investigation for their functional contribution to typical brain development and discusses how this knowledge can be harvested for therapeutic intervention in challenging neurodevelopmental disorders.

Functional Resilience of Mutually Repressing Motifs Embedded in Larger Networks

Elucidating the design principles of regulatory networks driving cellular decision-making has important implications for understanding cell differentiation and guiding the design of synthetic circuits. Mutually repressing feedback loops between 'master regulators' of cell fates can exhibit multistable dynamics enabling "single-positive" phenotypes: (high A, low B) and (low A, high B) for a toggle switch, and (high A, low B, low C), (low A, high B, low C) and (low A, low B, high C) for a toggle triad. However, the dynamics of these two motifs have been interrogated in isolation in silico, but in vitro and in vivo, they often operate while embedded in larger regulatory networks. Here, we embed these motifs in complex larger networks of varying sizes and connectivity to identify hallmarks under which these motifs maintain their canonical dynamical behavior. We show that an increased number of incoming edges onto a motif leads to a decay in their canonical stand-alone behaviors. We also show that this decay can be exacerbated by adding self-inhibition but not self-activation loops on the 'master regulators'. These observations offer insights into the design principles of biological networks containing these motifs and can help devise optimal strategies for the integration of these motifs into larger synthetic networks.

spliceJAC: transition genes and state-specific gene regulation from single-cell transcriptome data

Extracting dynamical information from single-cell transcriptomics is a novel task with the promise to advance our understanding of cell state transition and interactions between genes. Yet, theory-oriented, bottom-up approaches that consider differences among cell states are largely lacking. Here, we present spliceJAC, a method to quantify the multivariate mRNA splicing from single-cell RNA sequencing (scRNA-seq). spliceJAC utilizes the unspliced and spliced mRNA count matrices to constructs cell state-specific gene-gene regulatory interactions and applies stability analysis to predict putative driver genes critical to the transitions between cell states. By applying spliceJAC to biological systems including pancreas endothelium development and epithelial-mesenchymal transition (EMT) in A549 lung cancer cells, we predict genes that serve specific signaling roles in different cell states, recover important differentially expressed genes in agreement with pre-existing analysis, and predict new transition genes that are either exclusive or shared between different cell state transitions.

Inferring structural and dynamical properties of gene networks from data with deep learning

The reconstruction of gene regulatory networks (GRNs) from data is vital in systems biology. Although different approaches have been proposed to infer causality from data, some challenges remain, such as how to accurately infer the direction and type of interactions, how to deal with complex network involving multiple feedbacks, as well as how to infer causality between variables from real-world data, especially single cell data. Here, we tackle these problems by deep neural networks (DNNs). The underlying regulatory network for different systems (gene regulations, ecology, diseases, development) can be successfully reconstructed from trained DNN models. We show that DNN is superior to existing approaches including Boolean network, Random Forest and partial cross mapping for network inference. Further, by interrogating the ensemble DNN model trained from single cell data from dynamical system perspective, we are able to unravel complex cell fate dynamics during preimplantation development. We also propose a data-driven approach to quantify the energy landscape for gene regulatory systems, by combining DNN with the partial self-consistent mean field approximation (PSCA) approach. We anticipate the proposed method can be applied to other fields to decipher the underlying dynamical mechanisms of systems from data.

DURIAN: an integrative deconvolution and imputation method for robust signaling analysis of single-cell transcriptomics data

Single-cell RNA sequencing trades read-depth for dimensionality, often leading to loss of critical signaling gene information that is typically present in bulk data sets. We introduce DURIAN (Deconvolution and mUltitask-Regression-based ImputAtioN), an integrative method for recovery of gene expression in single-cell data. Through systematic benchmarking, we demonstrate the accuracy, robustness and empirical convergence of DURIAN using both synthetic and published data sets. We show that use of DURIAN improves single-cell clustering, low-dimensional embedding, and recovery of intercellular signaling networks. Our study resolves several inconsistent results of cell-cell communication analysis using single-cell or bulk data independently. The method has broad application in biomarker discovery and cell signaling analysis using single-cell transcriptomics data sets.

Capybara: A computational tool to measure cell identity and fate transitions

Measuring cell identity in development, disease, and reprogramming is challenging as cell types and states are in continual transition. Here, we present Capybara, a computational tool to classify discrete cell identity and intermediate "hybrid" cell states, supporting a metric to quantify cell fate transition dynamics. We validate hybrid cells using experimental lineage tracing data to demonstrate the multi-lineage potential of these intermediate cell states. We apply Capybara to diagnose shortcomings in several cell engineering protocols, identifying hybrid states in cardiac reprogramming and off-target identities in motor neuron programming, which we alleviate by adding exogenous signaling factors. Further, we establish a putative in vivo correlate for induced endoderm progenitors. Together, these results showcase the utility of Capybara to dissect cell identity and fate transitions, prioritizing interventions to enhance the efficiency and fidelity of stem cell engineering.

Triku: a feature selection method based on nearest neighbors for single-cell data

BACKGROUND

Feature selection is a relevant step in the analysis of single-cell RNA sequencing datasets. Most of the current feature selection methods are based on general univariate descriptors of the data such as the dispersion or the percentage of zeros. Despite the use of correction methods, the generality of these feature selection methods biases the genes selected towards highly expressed genes, instead of the genes defining the cell populations of the dataset.

RESULTS

Triku is a feature selection method that favors genes defining the main cell populations. It does so by selecting genes expressed by groups of cells that are close in the k-nearest neighbor graph. The expression of these genes is higher than the expected expression if the k-cells were chosen at random. Triku efficiently recovers cell populations present in artificial and biological benchmarking datasets, based on adjusted Rand index, normalized mutual information, supervised classification, and silhouette coefficient measurements. Additionally, gene sets selected by triku are more likely to be related to relevant Gene Ontology terms and contain fewer ribosomal and mitochondrial genes.

CONCLUSION

Triku is developed in Python 3 and is available at https://github.com/alexmascension/triku.

Epithelial-to-mesenchymal transition proceeds through directional destabilization of multidimensional attractor

How a cell changes from one stable phenotype to another one is a fundamental problem in developmental and cell biology. Mathematically, a stable phenotype corresponds to a stable attractor in a generally multi-dimensional state space, which needs to be destabilized so the cell relaxes to a new attractor. Two basic mechanisms for destabilizing a stable fixed point, pitchfork and saddle-node bifurcations, have been extensively studied theoretically; however, direct experimental investigation at the single-cell level remains scarce. Here, we performed live cell imaging studies and analyses in the framework of dynamical systems theories on epithelial-to-mesenchymal transition (EMT). While some mechanistic details remain controversial, EMT is a cell phenotypic transition (CPT) process central to development and pathology. Through time-lapse imaging we recorded single cell trajectories of human A549/Vim-RFP cells undergoing EMT induced by different concentrations of exogenous TGF-β in a multi-dimensional cell feature space. The trajectories clustered into two distinct groups, indicating that the transition dynamics proceeds through parallel paths. We then reconstructed the reaction coordinates and the corresponding quasi-potentials from the trajectories. The potentials revealed a plausible mechanism for the emergence of the two paths where the original stable epithelial attractor collides with two saddle points sequentially with increased TGF-β concentration, and relaxes to a new one. Functionally, the directional saddle-node bifurcation ensures a CPT proceeds towards a specific cell type, as a mechanistic realization of the canalization idea proposed by Waddington.

Cells with the same genetic code can take on many different formss, or phenotypes, which have distinct roles and appearances. Sometimes cells switch from one phenotype to another as part of healthy growth or during disease. One such change is the epithelial-to-mesenchymal transition (EMT), which is involved in fetal development, wound healing and the spread of cancer cells.

During EMT, closely connected epithelial cells detach from one another and change into mesenchymal cells that are able to migrate. Cells undergo a number of changes during this transition; however, the path they take to reach their new form is not entirely clear. For instance, do all cells follow the same route, or are there multiple ways that cells can shift from one state to the next?

To address this question, Wang et al. studied individual lung cancer cells that had been treated with a protein that drives EMT. The cells were then imaged at regular intervals over the course of two to three days to see how they changed in response to different concentrations of protein. Using a mathematical analysis designed to study chemical reactions, Wang et al. showed that the cells transform into the mesenchymal phenotype through two main routes.

This result suggests that attempts to prevent EMT, in cancer treatment for instance, would require blocking both paths taken by the cells. This information could be useful for biomedical researchers trying to regulate the EMT process. The quantitative approach of this study could also help physicists and mathematicians study other types of transition that occur in biology.

High content Image Analysis to study phenotypic heterogeneity in endothelial cell monolayers

Endothelial cells (EC) are heterogeneous across and within tissues, reflecting distinct, specialised functions. EC heterogeneity has been proposed to underpin EC plasticity independently from vessel microenvironments. However, heterogeneity driven by contact-dependent or short-range cell-cell crosstalk cannot be evaluated with single cell transcriptomic approaches as spatial and contextual information is lost. Nonetheless, quantification of EC heterogeneity and understanding of its molecular drivers is key to developing novel therapeutics for cancer, cardiovascular diseases and for revascularisation in regenerative medicine. Here, we developed an EC profiling tool (ECPT) to examine individual cells within intact monolayers. We used ECPT to characterise different phenotypes in arterial, venous and microvascular EC populations. In line with other studies, we measured heterogeneity in terms of cell cycle, proliferation, and junction organisation. ECPT uncovered a previously under-appreciated single-cell heterogeneity in NOTCH activation. We correlated cell proliferation with different NOTCH activation states at the single cell and population levels. The positional and relational information extracted with our novel approach is key to elucidating the molecular mechanisms underpinning EC heterogeneity.

How a cell decides its own fate: a single-cell view of molecular mechanisms and dynamics of cell-type specification

On its path from a fertilized egg to one of the many cell types in a multicellular organism, a cell turns the blank canvas of its early embryonic state into a molecular profile fine-tuned to achieve a vital organismal function. This remarkable transformation emerges from the interplay between dynamically changing external signals, the cell's internal, variable state, and tremendously complex molecular machinery; we are only beginning to understand. Recently developed single-cell omics techniques have started to provide an unprecedented, comprehensive view of the molecular changes during cell-type specification and promise to reveal the underlying gene regulatory mechanism. The exponentially increasing amount of quantitative molecular data being created at the moment is slated to inform predictive, mathematical models. Such models can suggest novel ways to manipulate cell types experimentally, which has important biomedical applications. This review is meant to give the reader a starting point to participate in this exciting phase of molecular developmental biology. We first introduce some of the principal molecular players involved in cell-type specification and discuss the important organizing ability of biomolecular condensates, which has been discovered recently. We then review some of the most important single-cell omics methods and relevant findings they produced. We devote special attention to the dynamics of the molecular changes and discuss methods to measure them, most importantly lineage tracing. Finally, we introduce a conceptual framework that connects all molecular agents in a mathematical model and helps us make sense of the experimental data.

Landscape and kinetic path quantify critical transitions in epithelial-mesenchymal transition

Epithelial-mesenchymal transition (EMT), a basic developmental process that might promote cancer metastasis, has been studied from various perspectives. Recently, the early warning theory has been used to anticipate critical transitions in EMT from mathematical modeling. However, the underlying mechanisms of EMT involving complex molecular networks remain to be clarified. Especially, how to quantify the global stability and stochastic transition dynamics of EMT and what the underlying mechanism for early warning theory in EMT is remain to be fully clarified. To address these issues, we constructed a comprehensive gene regulatory network model for EMT and quantified the corresponding potential landscape. The landscape for EMT displays multiple stable attractors, which correspond to E, M, and some other intermediate states. Based on the path-integral approach, we identified the most probable transition paths of EMT, which are supported by experimental data. Correspondingly, the results of transition actions demonstrated that intermediate states can accelerate EMT, consistent with recent studies. By integrating the landscape and path with early warning concept, we identified the potential barrier height from the landscape as a global and more accurate measure for early warning signals to predict critical transitions in EMT. The landscape results also provide an intuitive and quantitative explanation for the early warning theory. Overall, the landscape and path results advance our mechanistic understanding of dynamical transitions and roles of intermediate states in EMT, and the potential barrier height provides a new, to our knowledge, measure for critical transitions and quantitative explanations for the early warning theory.

Noise distorts the epigenetic landscape and shapes cell-fate decisions

The Waddington epigenetic landscape has become an iconic representation of the cellular differentiation process. Recent single-cell transcriptomic data provide new opportunities for quantifying this originally conceptual tool, offering insight into the gene regulatory networks underlying cellular development. While many methods for constructing the landscape have been proposed, by far the most commonly employed approach is based on computing the landscape as the negative logarithm of the steady-state probability distribution. Here, we use simple models to highlight the complexities and limitations that arise when reconstructing the potential landscape in the presence of stochastic fluctuations. We consider how the landscape changes in accordance with different stochastic systems and show that it is the subtle interplay between the deterministic and stochastic components of the system that ultimately shapes the landscape. We further discuss how the presence of noise has important implications for the identifiability of the regulatory dynamics from experimental data. A record of this paper's transparent peer review process is included in the supplemental information.

Dissecting transition cells from single-cell transcriptome data through multiscale stochastic dynamics

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.

VeloViz: RNA velocity-informed embeddings for visualizing cellular trajectories

MOTIVATION

Single-cell transcriptomics profiling technologies enable genome-wide gene expression measurements in individual cells but can currently only provide a static snapshot of cellular transcriptional states. RNA velocity analysis can help infer cell state changes using such single-cell transcriptomics data. To interpret these cell state changes inferred from RNA velocity analysis as part of underlying cellular trajectories, current approaches rely on visualization with principal components, t-distributed stochastic neighbor embedding and other 2D embeddings derived from the observed single-cell transcriptional states. However, these 2D embeddings can yield different representations of the underlying cellular trajectories, hindering the interpretation of cell state changes.

RESULTS

We developed VeloViz to create RNA velocity-informed 2D and 3D embeddings from single-cell transcriptomics data. Using both real and simulated data, we demonstrate that VeloViz embeddings are able to capture underlying cellular trajectories across diverse trajectory topologies, even when intermediate cell states may be missing. By considering the predicted future transcriptional states from RNA velocity analysis, VeloViz can help visualize a more reliable representation of underlying cellular trajectories.

AVAILABILITY AND IMPLEMENTATION

Source code is available on GitHub (https://github.com/JEFworks-Lab/veloviz) and Bioconductor (https://bioconductor.org/packages/veloviz) with additional tutorials at https://JEF.works/veloviz/. Datasets used can be found on Zenodo (https://doi.org/10.5281/zenodo.4632471).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Modeling the effects of EMT-immune dynamics on carcinoma disease progression

During progression from carcinoma in situ to an invasive tumor, the immune system is engaged in complex sets of interactions with various tumor cells. Tumor cell plasticity alters disease trajectories via epithelial-to-mesenchymal transition (EMT). Several of the same pathways that regulate EMT are involved in tumor-immune interactions, yet little is known about the mechanisms and consequences of crosstalk between these regulatory processes. Here we introduce a multiscale evolutionary model to describe tumor-immune-EMT interactions and their impact on epithelial cancer progression from in situ to invasive disease. Through simulation of patient cohorts in silico, the model predicts that a controllable region maximizes invasion-free survival. This controllable region depends on properties of the mesenchymal tumor cell phenotype: its growth rate and its immune-evasiveness. In light of the model predictions, we analyze EMT-inflammation-associated data from The Cancer Genome Atlas, and find that association with EMT worsens invasion-free survival probabilities. This result supports the predictions of the model, and leads to the identification of genes that influence outcomes in bladder and uterine cancer, including FGF pathway members. These results suggest new means to delay disease progression, and demonstrate the importance of studying cancer-immune interactions in light of EMT.

Cross-species integration of single-cell RNA-seq resolved alveolar-epithelial transitional states in idiopathic pulmonary fibrosis

Single-cell transcriptomics analyses of the fibrotic lung uncovered two cell states critical to lung injury recovery in the alveolar epithelium-a reparative transitional cell state in the mouse and a disease-specific cell state (KRT5^-/KRT17⁺) in human idiopathic pulmonary fibrosis (IPF). The murine transitional cell state lies between the differentiation from type 2 (AT2) to type 1 pneumocyte (AT1), and the human KRT5^-/KRT17⁺ cell state may arise from the dysregulation of this differentiation process. We review major findings of single-cell transcriptomics analyses of the fibrotic lung and reanalyzed data from seven single-cell RNA sequencing studies of human and murine models of IPF, focusing on the alveolar epithelium. Our comparative and cross-species single-cell transcriptomics analyses allowed us to further delineate the differentiation trajectories from AT2 to AT1 and AT2 to the KRT5^-/KRT17⁺ cell state. We observed AT1 cells in human IPF retain the transcriptional signature of the murine transitional cell state. Using pseudotime analysis, we recapitulated the differentiation trajectories from AT2 to AT1 and from AT2 to KRT5^-/KRT17⁺ cell state in multiple human IPF studies. We further delineated transcriptional programs underlying cell-state transitions and determined the molecular phenotypes at terminal differentiation. We hypothesize that in addition to the reactivation of developmental programs (SOX4, SOX9), senescence (TP63, SOX4) and the Notch pathway (HES1) are predicted to steer intermediate progenitors to the KRT5^-/KRT17⁺ cell state. Our analyses suggest that activation of SMAD3 later in the differentiation process may explain the fibrotic molecular phenotype typical of KRT5^-/KRT17⁺ cells.

On the origins and conceptual frameworks of natural plasticity-Lessons from single-cell models in C. elegans

How flexible are cell identities? This problem has fascinated developmental biologists for several centuries and can be traced back to Abraham Trembley's pioneering manipulations of Hydra to test its regeneration abilities in the 1700s. Since the cell theory in the mid-19th century, developmental biology has been dominated by a single framework in which embryonic cells are committed to specific cell fates, progressively and irreversibly acquiring their differentiated identities. This hierarchical, unidirectional and irreversible view of cell identity has been challenged in the past decades through accumulative evidence that many cell types are more plastic than previously thought, even in intact organisms. The paradigm shift introduced by such plasticity calls into question several other key traditional concepts, such as how to define a differentiated cell or more generally cellular identity, and has brought new concepts, such as distinct cellular states. In this review, we want to contribute to this representation by attempting to clarify the conceptual and theoretical frameworks of cell plasticity and identity. In the context of these new frameworks we describe here an atlas of natural plasticity of cell identity in C. elegans, including our current understanding of the cellular and molecular mechanisms at play. The worm further provides interesting cases at the borderlines of cellular plasticity that highlight the conceptual challenges still ahead. We then discuss a set of future questions and perspectives arising from the studies of natural plasticity in the worm that are shared with other reprogramming and plasticity events across phyla.

An expanded universe of cancer targets

The characterization of cancer genomes has provided insight into somatically altered genes across tumors, transformed our understanding of cancer biology, and enabled tailoring of therapeutic strategies. However, the function of most cancer alleles remains mysterious, and many cancer features transcend their genomes. Consequently, tumor genomic characterization does not influence therapy for most patients. Approaches to understand the function and circuitry of cancer genes provide complementary approaches to elucidate both oncogene and non-oncogene dependencies. Emerging work indicates that the diversity of therapeutic targets engendered by non-oncogene dependencies is much larger than the list of recurrently mutated genes. Here we describe a framework for this expanded list of cancer targets, providing novel opportunities for clinical translation.

Phylodynamics for cell biologists

Multicellular organisms are composed of cells connected by ancestry and descent from progenitor cells. The dynamics of cell birth, death, and inheritance within an organism give rise to the fundamental processes of development, differentiation, and cancer. Technical advances in molecular biology now allow us to study cellular composition, ancestry, and evolution at the resolution of individual cells within an organism or tissue. Here, we take a phylogenetic and phylodynamic approach to single-cell biology. We explain how "tree thinking" is important to the interpretation of the growing body of cell-level data and how ecological null models can benefit statistical hypothesis testing. Experimental progress in cell biology should be accompanied by theoretical developments if we are to exploit fully the dynamical information in single-cell data.

Differential Network Analysis Reveals Regulatory Patterns in Neural Stem Cell Fate Decision

Deciphering regulatory patterns of neural stem cell (NSC) differentiation with multiple stages is essential to understand NSC differentiation mechanisms. Recent single-cell transcriptome datasets became available at individual differentiation. However, a systematic and integrative analysis of multiple datasets at multiple temporal stages of NSC differentiation is lacking. In this study, we propose a new method integrating prior information to construct three gene regulatory networks at pair-wise stages of transcriptome and apply this method to investigate five NSC differentiation paths on four different single-cell transcriptome datasets. By constructing gene regulatory networks for each path, we delineate their regulatory patterns via differential topology and network diffusion analyses. We find 12 common differentially expressed genes among the five NSC differentiation paths, with one common regulatory pattern (Gsk3b_App_Cdk5) shared by all paths. The identified regulatory pattern, partly supported by previous experimental evidence, is essential to all differentiation paths, but it plays a different role in each path when regulating other genes. Together, our integrative analysis provides both common and specific regulatory mechanisms for each of the five NSC differentiation paths.

Noise and the molecular processes underlying cell fate decision-making

Cell fate decision-making events involve the interplay of many molecular processes, ranging from signal transduction to genetic regulation, as well as a set of molecular and physiological feedback loops. Each aspect offers a rich field of investigation in its own right, but to understand the whole process, even in simple terms, we need to consider them together. Here we attempt to characterise this process by focussing on the roles of noise during cell fate decisions. We use a range of recent results to develop a view of the sequence of events by which a cell progresses from a pluripotent or multipotent to a differentiated state: chromatin organisation, transcription factor stoichiometry, and cellular signalling all change during this progression, and all shape cellular variability, which becomes maximal at the transition state.

A Spectrum of Cell States During the Epithelial-to-Mesenchymal Transition

The epithelial-to-mesenchymal transition (EMT) encompasses a complex cascade of events through which a cell transits to reduce its epithelial characteristics and become migratory. Classically, this transition has been considered complete upon loss of molecular markers characteristic of an "epithelial" state and acquisition of those associated with "mesenchymal" cells. Recently, however, evidence from both developmental and cancer EMT contexts suggest that cells undergoing EMT are often heterogeneous, concomitantly expressing both epithelial and mesenchymal markers to varying degrees; rather, cells frequently display a "partial" EMT phenotype and do not necessarily require full "mesenchymalization" to become migratory. Here, we offer a brief perspective on recent important advances in our fundamental understanding of the spectrum of cellular states that occur during partial EMT in the context of development and cancer metastasis.

Inference and multiscale model of epithelial-to-mesenchymal transition via single-cell transcriptomic data

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.

DNF: A differential network flow method to identify rewiring drivers for gene regulatory networks

Differential network analysis has become an important approach in identifying driver genes in development and disease. However, most studies capture only local features of the underlying gene-regulatory network topology. These approaches are vulnerable to noise and other changes which mask driver-gene activity. Therefore, methods are urgently needed which can separate the impact of true regulatory elements from stochastic changes and downstream effects. We propose the differential network flow (DNF) method to identify key regulators of progression in development or disease. Given the network representation of consecutive biological states, DNF quantifies the essentiality of each node by differences in the distribution of network flow, which are capable of capturing comprehensive topological differences from local to global feature domains. DNF achieves more accurate driver-gene identification than other state-of-the-art methods when applied to four human datasets from The Cancer Genome Atlas and three single-cell RNA-seq datasets of murine neural and hematopoietic differentiation. Furthermore, we predict key regulators of crosstalk between separate networks underlying both neuronal differentiation and the progression of neurodegenerative disease, among which APP is predicted as a driver gene of neural stem cell differentiation. Our method is a new approach for quantifying the essentiality of genes across networks of different biological states.

Multi-stability in cellular differentiation enabled by a network of three mutually repressing master regulators

Identifying the design principles of complex regulatory networks driving cellular decision-making remains essential to decode embryonic development as well as enhance cellular reprogramming. A well-studied network motif involved in cellular decision-making is a toggle switch-a set of two opposing transcription factors A and B, each of which is a master regulator of a specific cell fate and can inhibit the activity of the other. A toggle switch can lead to two possible states-(high A, low B) and (low A, high B)-and drives the 'either-or' choice between these two cell fates for a common progenitor cell. However, the principles of coupled toggle switches remain unclear. Here, we investigate the dynamics of three master regulators A, B and C inhibiting each other, thus forming three-coupled toggle switches to form a toggle triad. Our simulations show that this toggle triad can lead to co-existence of cells into three differentiated 'single positive' phenotypes-(high A, low B, low C), (low A, high B, low C) and (low A, low B, high C). Moreover, the hybrid or 'double positive' phenotypes-(high A, high B, low C), (low A, high B, high C) and (high A, low B, high C)-can coexist together with 'single positive' phenotypes. Including self-activation loops on A, B and C can increase the frequency of 'double positive' states. Finally, we apply our results to understand cellular decision-making in terms of differentiation of naive CD4⁺ T cells into Th1, Th2 and Th17 states, where hybrid Th1/Th2 and hybrid Th1/Th17 cells have been reported in addition to the Th1, Th2 and Th17 ones. Our results offer novel insights into the design principles of a multi-stable network topology and provide a framework for synthetic biology to design tristable systems.

DensityPath: an algorithm to visualize and reconstruct cell state-transition path on density landscape for single-cell RNA sequencing data

MOTIVATION

Visualizing and reconstructing cell developmental trajectories intrinsically embedded in high-dimensional expression profiles of single-cell RNA sequencing (scRNA-seq) snapshot data are computationally intriguing, but challenging.

RESULTS

We propose DensityPath, an algorithm allowing (i) visualization of the intrinsic structure of scRNA-seq data on an embedded 2-d space and (ii) reconstruction of an optimal cell state-transition path on the density landscape. DensityPath powerfully handles high dimensionality and heterogeneity of scRNA-seq data by (i) revealing the intrinsic structures of data, while adopting a non-linear dimension reduction algorithm, termed elastic embedding, which can preserve both local and global structures of the data; and (ii) extracting the topological features of high-density, level-set clusters from a single-cell multimodal density landscape of transcriptional heterogeneity, as the representative cell states. DensityPath reconstructs the optimal cell state-transition path by finding the geodesic minimum spanning tree of representative cell states on the density landscape, establishing a least action path with the minimum-transition-energy of cell fate decisions. We demonstrate that DensityPath can ably reconstruct complex trajectories of cell development, e.g. those with multiple bifurcating and trifurcating branches, while maintaining computational efficiency. Moreover, DensityPath has high accuracy for pseudotime calculation and branch assignment on real scRNA-seq, as well as simulated datasets. DensityPath is robust to parameter choices, as well as permutations of data.

AVAILABILITY AND IMPLEMENTATION

DensityPath software is available at https://github.com/ucasdp/DensityPath.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

scRCMF: Identification of Cell Subpopulations and Transition States From Single-Cell Transcriptomes

Single cell technologies provide an unprecedented opportunity to explore the heterogeneity in a biological process at the level of single cells. One major challenge in analyzing single cell data is to identify cell subpopulations, stable cell states, and cells in transition between states. To elucidate the transition mechanisms in cell fate dynamics, it is highly desirable to quantitatively characterize cellular states and intermediate states. Here, we present scRCMF, an unsupervised method that identifies stable cell states and transition cells by adopting a nonlinear optimization model that infers the latent substructures from a gene-cell matrix. We incorporate a random coefficient matrix-based regularization into the standard nonnegative matrix decomposition model to improve the reliability and stability of estimating latent substructures. To quantify the transition capability of each cell, we propose two new measures: single-cell transition entropy (scEntropy) and transition probability (scTP). When applied to two simulated and three published scRNA-seq datasets, scRCMF not only successfully captures multiple subpopulations and transition processes in large-scale data, but also identifies transition states and some known marker genes associated with cell state transitions and subpopulations. Furthermore, the quantity scEntropy is found to be significantly higher for transition cells than other cellular states during the global differentiation, and the scTP predicts the "fate decisions" of transition cells within the transition. The present study provides new insights into transition events during differentiation and development.

A plausible accelerating function of intermediate states in cancer metastasis

Epithelial-to-mesenchymal transition (EMT) is a fundamental cellular process and plays an essential role in development, tissue regeneration, and cancer metastasis. Interestingly, EMT is not a binary process but instead proceeds with multiple partial intermediate states. However, the functions of these intermediate states are not fully understood. Here, we focus on a general question about how the number of partial EMT states affects cell transformation. First, by fitting a hidden Markov model of EMT with experimental data, we propose a statistical mechanism for EMT in which many unobservable microstates may exist within one of the observable macrostates. Furthermore, we find that increasing the number of intermediate states can accelerate the EMT process and that adding parallel paths or transition layers may accelerate the process even further. Last, a stabilized intermediate state traps cells in one partial EMT state. This work advances our understanding of the dynamics and functions of EMT plasticity during cancer metastasis.

Integrative Transcriptomic Analysis Reveals a Multiphasic Epithelial-Mesenchymal Spectrum in Cancer and Non-tumorigenic Cells

Epithelial-mesenchymal transition (EMT), the conversion between rigid epithelial cells and motile mesenchymal cells, is a reversible cellular process involved in tumorigenesis, metastasis, and chemoresistance. Numerous studies have found that several types of tumor cells show a high degree of cell-to-cell heterogeneity in terms of their gene expression signatures and cellular phenotypes related to EMT. Recently, the prevalence and importance of partial or intermediate EMT states have been reported. It is unclear, however, whether there is a general pattern of cancer cell distribution in terms of the overall expression of epithelial-related genes and mesenchymal-related genes, and how this distribution is related to EMT process in normal cells. In this study, we performed integrative transcriptomic analysis that combines cancer cell transcriptomes, time course data of EMT in non-tumorigenic epithelial cells, and epithelial cells with perturbations of key EMT factors. Our statistical analysis shows that cancer cells are widely distributed in the EMT spectrum, and the majority of these cells can be described by an EMT path that connects the epithelial and the mesenchymal states via a hybrid expression region in which both epithelial genes and mesenchymal genes are highly expressed overall. We found that key patterns of this EMT path are observed in EMT progression in non-tumorigenic cells and that transcription factor ZEB1 plays a key role in defining this EMT path via diverse gene regulatory circuits connecting to epithelial genes. We performed Gene Set Variation Analysis to show that the cancer cells at hybrid EMT states also possess hybrid cellular phenotypes with both high migratory and high proliferative potentials. Our results reveal critical patterns of cancer cells in the EMT spectrum and their relationship to the EMT process in normal cells, and provide insights into the mechanistic basis of cancer cell heterogeneity and plasticity.

Cell lineage and communication network inference via optimization for single-cell transcriptomics

The use of single-cell transcriptomics has become a major approach to delineate cell subpopulations and the transitions between them. While various computational tools using different mathematical methods have been developed to infer clusters, marker genes, and cell lineage, none yet integrate these within a mathematical framework to perform multiple tasks coherently. Such coherence is critical for the inference of cell–cell communication, a major remaining challenge. Here, we present similarity matrix-based optimization for single-cell data analysis (SoptSC), in which unsupervised clustering, pseudotemporal ordering, lineage inference, and marker gene identification are inferred via a structured cell-to-cell similarity matrix. SoptSC then predicts cell–cell communication networks, enabling reconstruction of complex cell lineages that include feedback or feedforward interactions. Application of SoptSC to early embryonic development, epidermal regeneration, and hematopoiesis demonstrates robust identification of subpopulations, lineage relationships, and pseudotime, and prediction of pathway-specific cell communication patterns regulating processes of development and differentiation.

Calibration, Selection and Identifiability Analysis of a Mathematical Model of the in vitro Erythropoiesis in Normal and Perturbed Contexts

The in vivo erythropoiesis, which is the generation of mature red blood cells in the bone marrow of whole organisms, has been described by a variety of mathematical models in the past decades. However, the in vitro erythropoiesis, which produces red blood cells in cultures, has received much less attention from the modelling community. In this paper, we propose the first mathematical model of in vitro erythropoiesis. We start by formulating different models and select the best one at fitting experimental data of in vitro erythropoietic differentiation obtained from chicken erythroid progenitor cells. It is based on a set of linear ODE, describing 3 hypothetical populations of cells at different stages of differentiation. We then compute confidence intervals for all of its parameters estimates, and conclude that our model is fully identifiable. Finally, we use this model to compute the effect of a chemical drug called Rapamycin, which affects all states of differentiation in the culture, and relate these effects to specific parameter variations. We provide the first model for the kinetics of in vitro cellular differentiation which is proven to be identifiable. It will serve as a basis for a model which will better account for the variability which is inherent to the experimental protocol used for the model calibration.

Uncoupling Traditional Functionalities of Metastasis: The Parting of Ways with Real-Time Assays

The experimental evaluation of metastasis overly focuses on the gain of migratory and invasive properties, while disregarding the contributions of cellular plasticity, extra-cellular matrix heterogeneity, niche interactions, and tissue architecture. Traditional cell-based assays often restrict the inclusion of these processes and warrant the implementation of approaches that provide an enhanced spatiotemporal resolution of the metastatic cascade. Time lapse imaging represents such an underutilized approach in cancer biology, especially in the context of disease progression. The inclusion of time lapse microscopy and microfluidic devices in routine assays has recently discerned several nuances of the metastatic cascade. Our review emphasizes that a complete comprehension of metastasis in view of evolving ideologies necessitates (i) the use of appropriate, context-specific assays and understanding their inherent limitations; (ii) cautious derivation of inferences to avoid erroneous/overestimated clinical extrapolations; (iii) corroboration between multiple assay outputs to gauge metastatic potential; and (iv) the development of protocols with improved in situ implications. We further believe that the adoption of improved quantitative approaches in these assays can generate predictive algorithms that may expedite therapeutic strategies targeting metastasis via the development of disease relevant model systems. Such approaches could potentiate the restructuring of the cancer metastasis paradigm through an emphasis on the development of next-generation real-time assays.

Investigating epithelial-to-mesenchymal transition with integrated computational and experimental approaches

The transition between epithelial and mesenchymal (EMT) is a fundamental cellular process that plays critical roles in development, cancer metastasis, and tissue wound healing. EMT is not a binary process but involves multiple partial EMT states that give rise to a high degree of cell state plasticity. Here, we first reviewed several studies on theoretical predictions and experimental verification of these intermediate states, the role of partial EMT on kidney fibrosis development, and how quantitative signaling information controls cell commitment to partial or full EMT upon transient signals. Next, we summarized existing knowledge and open questions on the coupling between EMT and other biological processes, such as the cell cycle, epigenetic regulation, stemness, and apoptosis. Taken together, EMT is a model system that has attracted increasing interests for quantitative experimental and theoretical studies.

An enriched network motif family regulates multistep cell fate transitions with restricted reversibility

Multistep cell fate transitions with stepwise changes of transcriptional profiles are common to many developmental, regenerative and pathological processes. The multiple intermediate cell lineage states can serve as differentiation checkpoints or branching points for channeling cells to more than one lineages. However, mechanisms underlying these transitions remain elusive. Here, we explored gene regulatory circuits that can generate multiple intermediate cellular states with stepwise modulations of transcription factors. With unbiased searching in the network topology space, we found a motif family containing a large set of networks can give rise to four attractors with the stepwise regulations of transcription factors, which limit the reversibility of three consecutive steps of the lineage transition. We found that there is an enrichment of these motifs in a transcriptional network controlling the early T cell development, and a mathematical model based on this network recapitulates multistep transitions in the early T cell lineage commitment. By calculating the energy landscape and minimum action paths for the T cell model, we quantified the stochastic dynamics of the critical factors in response to the differentiation signal with fluctuations. These results are in good agreement with experimental observations and they suggest the stable characteristics of the intermediate states in the T cell differentiation. These dynamical features may help to direct the cells to correct lineages during development. Our findings provide general design principles for multistep cell linage transitions and new insights into the early T cell development. The network motifs containing a large family of topologies can be useful for analyzing diverse biological systems with multistep transitions.

The functions of cells are dynamically controlled in many biological processes including development, regeneration and disease progression. Cell fate transition, or the switch of cellular functions, often involves multiple steps. The intermediate stages of the transition provide the biological systems with the opportunities to regulate the transitions in a precise manner. These transitions are controlled by key regulatory genes of which the expression shows stepwise patterns, but how the interactions of these genes can determine the multistep processes was unclear. Here, we present a comprehensive analysis on the design principles of gene circuits that govern multistep cell fate transition. We found a large network family with common structural features that can generate systems with the ability to control three consecutive steps of the transition. We found that this type of networks is enriched in a gene circuit controlling the development of T lymphocyte, a crucial type of immune cells. We performed mathematical modeling using this gene circuit and we recapitulated the stepwise and irreversible loss of stem cell properties of the developing T lymphocytes. Our findings can be useful to analyze a wide range of gene regulatory networks controlling multistep cell fate transitions.

Intermediate cell states in epithelial-to-mesenchymal transition

The transition of epithelial cells into a mesenchymal state (epithelial-to-mesenchymal transition or EMT) is a highly dynamic process implicated in various biological processes. During EMT, cells do not necessarily exist in 'pure' epithelial or mesenchymal states. There are cells with mixed (or hybrid) features of the two, which are termed as the intermediate cell states (ICSs). While the exact functions of ICS remain elusive, together with EMT it appears to play important roles in embryogenesis, tissue development, and pathological processes such as cancer metastasis. Recent single cell experiments and advanced mathematical modeling have improved our capability in identifying ICS and provided a better understanding of ICS in development and disease. Here, we review the recent findings related to the ICS in/or EMT and highlight the challenges in the identification and functional characterization of ICS.

Transition state characteristics during cell differentiation

Models describing the process of stem-cell differentiation are plentiful, and may offer insights into the underlying mechanisms and experimentally observed behaviour. Waddington's epigenetic landscape has been providing a conceptual framework for differentiation processes since its inception. It also allows, however, for detailed mathematical and quantitative analyses, as the landscape can, at least in principle, be related to mathematical models of dynamical systems. Here we focus on a set of dynamical systems features that are intimately linked to cell differentiation, by considering exemplar dynamical models that capture important aspects of stem cell differentiation dynamics. These models allow us to map the paths that cells take through gene expression space as they move from one fate to another, e.g. from a stem-cell to a more specialized cell type. Our analysis highlights the role of the transition state (TS) that separates distinct cell fates, and how the nature of the TS changes as the underlying landscape changes-change that can be induced by e.g. cellular signaling. We demonstrate that models for stem cell differentiation may be interpreted in terms of either a static or transitory landscape. For the static case the TS represents a particular transcriptional profile that all cells approach during differentiation. Alternatively, the TS may refer to the commonly observed period of heterogeneity as cells undergo stochastic transitions.