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.

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

Generalizing RNA Velocity to Transient Cell States Through Dynamical Modeling

Mapping transcriptomic vector fields of single cells

SDEvelo: a deep generative approach for transcriptional dynamics with cell-specific latent time and multivariate stochastic modeling

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

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

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

Spatiotemporally resolved transcriptomics reveals the subcellular RNA kinetic landscape

On the Mathematics of RNA Velocity I: Theoretical Analysis

Dissection and integration of bursty transcriptional dynamics for complex systems

RNA Velocity—current Challenges and Future Perspectives

Improving the study of RNA dynamics through advances in RNA-seq with metabolic labeling and nucleotide-recoding chemistry

On the Mathematics of RNA Velocity II: Algorithmic Aspects

Assessing RNA velocity stability across synthetic replicates using count splitting

Inferring Single-Cell RNA Kinetics from Various Biological Priors

Modeling single cell trajectory using forward-backward stochastic differential equations

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

Global transcription regulation revealed from dynamical correlations in time-resolved single-cell RNA-sequencing

Stochastic Simulation And Statistical Inference Platform For Visualization And Estimation Of Transcriptional Kinetics

SLAM-Drop-seq reveals mRNA kinetic rates throughout the cell cycle

TFvelo: gene regulation inspired RNA velocity estimation