Abstract:Background: Over the past decade, experimental procedures such as metabolic labeling for determining RNA turnover rates at the transcriptome-wide scale have been widely adopted and are now turning to single cell measurements. Several computational methods to estimate RNA synthesis, processing and degradation rates from such experiments have been suggested, but they all require several RNA sequencing samples. Here we present a method that can estimate those three rates from a single sample. Methods: Our method relies on the analytical solution to the Zeisel model of RNA dynamics. It was validated on metabolic labeling experiments performed on mouse embryonic stem cells. Resulting degradation rates were compared both to previously published rates on the same system and to a state-of-the-art method applied to the same data. Results: Our method is computationally efficient and outputs rates that correlate well with previously published data sets. Using it on a single sample, we were able to reproduce the observation that dynamic biological processes tend to involve genes with higher metabolic rates, while stable processes involve genes with lower rates. This supports the hypothesis that cells control not only the mRNA steady-state abundance, but also its responsiveness, i.e., how fast steady state is reached. Moreover, degradation rates obtained with our method compare favourably with the other tested method. Conclusions: In addition to saving experimental work and computational time, estimating rates for a single sample has several advantages. It does not require an error-prone normalization across samples and enables the use of replicates to estimate uncertainty and assess sample quality. Finally the method and theoretical results described here are general enough to be useful in other contexts such as nucleotide conversion methods and single cell metabolic labeling experiments.

Quantifying mRNA synthesis and decay rates using small RNAs

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

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

A spatiotemporally resolved atlas of mRNA decay in the embryo reveals differential regulation of mRNA stability across stages and cell types

Quantifying the noise in bursty gene expression under regulation by small RNAs

Ribosomes modulate transcriptome abundance via generalized frameshift and out-of-frame mRNA decay

Nuclear mRNA decay: regulatory networks that control gene expression

Ribosome dynamics and mRNA turnover, a complex relationship under constant cellular scrutiny

Improved RNA stability estimation through Bayesian modeling reveals most transcripts have sub-minute half-lives

Global SLAM-seq for accurate mRNA decay determination and identification of NMD targets

Phosphorylation-dependent tuning of mRNA deadenylation rates

The Ccr4-Not complex monitors the translating ribosome for codon optimality

Multi-species stochastic models: small RNA regulation of gene expression

Non-Invasive Analysis of Recombinant mRNA Stability in Escherichia coli by a Combination of Transcriptional Inducer Wash-Out and qRT-PCR

Subcellular mRNA kinetic modeling reveals nuclear retention as rate-limiting

A spatiotemporally resolved atlas of mRNA decay in the C. elegans embryo reveals differential regulation of mRNA stability across stages and cell types

Arabidopsis mRNA decay landscape arises from specialized RNA decay substrates, decapping-mediated feedback, and redundancy

Genome‐wide Study of Mrna Degradation and Transcript Elongation in E Scherichia Coli

Stochastic modeling of regulation of gene expression by multiple small RNAs

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

Determination of in vivo target search kinetics of regulatory noncoding RNA