Abstract:Experimental approaches for measuring single-cell gene expression can observe each cell at only one time point, requiring computational approaches for reconstructing the dynamics of gene expression during cell fate transitions. RNA velocity is a promising computational approach for this problem, but existing inference methods fail to capture key aspects of real data, limiting their utility. To address these limitations, we developed VeloVAE, a Bayesian model for RNA velocity inference. VeloVAE uses variational Bayesian inference to estimate the posterior distribution of latent time, latent cell state, and kinetic rate parameters for each cell. Our approach addresses key limitations of previous methods by inferring a global time and cell state value for each cell; explicitly modeling the emergence of multiple cell types; incorporating prior information such as time point labels; using scalable minibatch optimization; and quantifying parameter uncertainty. We show that VeloVAE significantly outperforms previous approaches in terms of data fit and accuracy of inferred differentiation directions. VeloVAE can also capture qualitative features of expression dynamics neglected by previous methods, including late induction, early repression, transcriptional boosts, and bifurcations. These improvements allow VeloVAE to accurately model gene expression dynamics in complex biological systems, including hematopoiesis, induced pluripotent stem cell reprogramming, neurogenesis, and organogenesis. We find that the latent time automatically inferred using all cells can even outperform pseudotime inferred using manually chosen cell subsets and root cells. We further use the inferred parameters to construct cell type transition graphs and fit branching differential equation models that describe the effects of cell type bifurcations on kinetic rate parameters.

Multi-omic Single-Cell Velocity Models Epigenome–transcriptome Interactions and Improves Cell Fate Prediction

Single-cell multi-omic velocity infers dynamic and decoupled gene regulation

Generalizing RNA Velocity to Transient Cell States Through Dynamical Modeling

Deep Generative Modeling of Transcriptional Dynamics for RNA Velocity Analysis in Single Cells

mmVelo: A deep generative model for estimating cell state-dependent dynamics across multiple modalities

PhyloVelo enhances transcriptomic velocity field mapping using monotonically expressed genes

TSvelo: Comprehensive RNA velocity inference by jointly modeling Transcription and Splicing

DeepVelo: deep learning extends RNA velocity to multi-lineage systems with cell-specific kinetics

RNA Velocity—current Challenges and Future Perspectives

Bayesian Inference of RNA Velocity from Multi-Lineage Single-Cell Data

TFvelo: gene regulation inspired RNA velocity estimation

GraphVelo allows inference of multi-modal single cell velocities and molecular mechanisms

NeuroVelo: interpretable learning of cellular dynamics from single-cell transcriptomic data

Challenges and Progress in RNA Velocity: Comparative Analysis Across Multiple Biological Contexts

RegVelo: gene-regulatory-informed dynamics of single cells

Mapping Cell Fate Transition in Space and Time

CRAK-Velo: Chromatin Accessibility Kinetics integration improves RNA Velocity estimation

scKINETICS: inference of regulatory velocity with single-cell transcriptomics data

Statistical inference with a manifold-constrained RNA velocity model uncovers cell cycle speed modulations

Multi-Experiment Nonlinear Mixed Effect Modeling of Single-Cell Translation Kinetics after Transfection

Mapping transcriptomic vector fields of single cells