Abstract:The time-resolved scRNA-seq (tscRNA-seq) provides the possibility to infer physically meaningful kinetic parameters, e.g., the transcription, splicing or RNA degradation rate constants with correct magnitudes, and RNA velocities by incorporating temporal information. Previous approaches utilizing the deterministic dynamics and steady-state assumption on gene expression states are insufficient to achieve favorable results for the data involving transient process. We present a dynamical approach, Storm (Stochastic models of RNA metabolic-labeling), to overcome these limitations by solving stochastic differential equations of gene expression dynamics. The derivation reveals that the new mRNA sequencing data obeys different types of cell-specific Poisson distributions when jointly considering both biological and cell-specific technical noise. Storm deals with measured counts data directly and extends the RNA velocity methodology based on metabolic labeling scRNA-seq data to transient stochastic systems. Furthermore, we relax the constant parameter assumption over genes/cells to obtain gene-cell-specific transcription/splicing rates and gene-specific degradation rates, thus revealing time-dependent and cell-state-specific transcriptional regulations. Storm will facilitate the study of the statistical properties of tscRNA-seq data, eventually advancing our understanding of the dynamic transcription regulation during development and disease. Intricate regulation of RNA biogenesis, such as RNA transcription, splicing and degradation, plays a critical role in most biological processes. Previous approaches have leveraged a deterministic model of spliced and unspliced RNAs to estimate kinetic parameters and quantify RNA velocity, the rate of changes in gene expression states in single cells. However, accurately estimating meaningful kinetic parameters and RNA velocity is hindered by biased capture of unspliced RNA and absence of temporal information in conventional scRNA-seq. Significant advances have been made in better measuring RNA kinetics with metabolic labeling enabled scRNA-seq, however, computational tools to analyze them lag far behind. Prior work of Dynamo provides one of the early solutions to properly model RNA metabolic labeling data, but its method still largely uses a deterministic model that only utilizes part of extreme cells and is unable to analyze datasets with significant transient dynamics. To address these challenges, we developed Storm that explicitly models transient stochastic RNA dynamics. Importantly, Storm models RNA kinetics with stochastic differential equations that explicitly account for biological and cell-specific noises. Storm is generally applicable to many metabolic labeling scRNA-seq datasets and we demonstrate excellent performance of Storm in fitting the data to capture the transient dynamics under various noise models.

Assessing RNA velocity stability across synthetic replicates using count splitting

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

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

Generalizing RNA Velocity to Transient Cell States Through Dynamical Modeling

Pumping the brakes on RNA velocity by understanding and interpreting RNA velocity estimates

Quantifying uncertainty in RNA velocity

RNA Velocity—current Challenges and Future Perspectives

On the Mathematics of RNA Velocity II: Algorithmic Aspects

On the Mathematics of RNA Velocity I: Theoretical Analysis

Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data

Improving the RNA velocity approach using long-read single cell sequencing

Improving the RNA Velocity Approach with Single-Cell RNA Lifecycle (nascent, Mature and Degrading RNAs) Sequencing Technologies.

RNA velocity of single cells

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

Accurate RNA velocity estimation based on multibatch network reveals complex lineage in batch scRNA-seq data

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

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

Storm: Incorporating transient stochastic dynamics to infer the RNA velocity with metabolic labeling information

Predicting the stability of profiling signals of small RNAs

Estimating RNA dynamics using one time point for one sample in a single-pulse metabolic labeling experiment