We employed our previously-described single-cell gene expression analysis CALISTA (Clustering And Lineage Inference in Single-Cell Transcriptional Analysis) to evaluate transcriptional uncertainty at the single-cell level using a stochastic mechanistic model of gene expression. We reconstructed a transcriptional uncertainty landscape during cell differentiation by visualizing single-cell transcriptional uncertainty surface over a two dimensional representation of the single-cell gene expression data. The reconstruction of transcriptional uncertainty landscapes for ten publicly available single-cell gene expression datasets from cell differentiation processes with linear, single or multi-branching cell lineage, reveals universal features in the cell differentiation trajectory that include: (i) a peak in single-cell uncertainty during transition states, and in systems with bifurcating differentiation trajectories, each branching point represents a state of high transcriptional uncertainty; (ii) a positive correlation of transcriptional uncertainty with transcriptional burst size and frequency; (iii) an increase in RNA velocity preceeding the increase in the cell transcriptional uncertainty. Finally, we provided biological interpretations of the universal rise-then-fall profile of the transcriptional uncertainty landscape, including a link with the Waddington's epigenetic landscape, that is generalizable to every cell differentiation system.
MAINCell differentiation is the process through which unspecialized stem cells become more specialized. Because of its importance in development, cellular repair, and organismal homeostasis, the molecular mechanisms of cell differentiation has been the subject of intense scrutiny. Since roughly 50 years ago -along with the promulgation of the central dogma of molecular biology by Francis Crick and the characterization of the lactose operon by François Jacob and Jacques Monod -the existence of a genetic program has become a prevailing explanation for the cell differentiation process. Although the details were originally not defined, at least not formally, such a genetic program purports a constellation of master genes (i.e., transcription factors) that orchestrate the transcription of downstream target genes in a precise spatiotemporal fashion, resulting in long-lasting alterations in the gene expression patterns 1-3 . A notable experimental evidence substantiating this view is the overexpression of myoD inducing a myogenic phenotype in seemingly naive cells 4 . Over the past few decades, the repertoire of such master genes across numerous stem cell systems, such as Nanog, Oct4, Sox2, BATF and MyoD, begin to coalesce 5-7 .Recent advances in single-cell technologies has revealed new aspects of the cell differentiation that are incompatible with the idea of ordered and programmed (i.e., deterministic) gene expression. More specifically, single-cell data paint a stochastic differentiation process that increases cell-to-cell variability of gene expression.