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.

pulseR: Versatile computational analysis of RNA turnover from metabolic labeling experiments

pulseTD: RNA life cycle dynamics analysis based on pulse model of 4sU-seq time course sequencing data

RNAflow: An Effective and Simple RNA-Seq Differential Gene Expression Pipeline Using Nextflow

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

Metabolic labeling of RNAs uncovers hidden features and dynamics of the Arabidopsis thaliana transcriptome

Halfpipe: a tool for analyzing metabolic labeling RNA-seq data to quantify RNA half-lives

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

DRUMMER-rapid detection of RNA modifications through comparative nanopore sequencing

Deep learning and direct sequencing of labeled RNA captures transcriptome dynamics

RNASeqR: an R package for automated two-group RNA-Seq analysis workflow

Impulse Model-Based Differential Expression Analysis of Time Course Sequencing Data

RNA-SeQC: RNA-seq metrics for quality control and process optimization

A Powerful and Flexible Approach to the Analysis of RNA Sequence Count Data

Count-based differential expression analysis of RNA sequencing data using R and Bioconductor

BP4RNAseq: a Babysitter Package for Retrospective and Newly Generated RNA-seq Data Analyses Using Both Alignment-Based and Alignment-Free Quantification Method.

Web-based platform for analysis of RNA folding from high throughput chemical probing data

RNA-seq analysis is easy as 1-2-3 with limma, Glimma and edgeR

A protocol for time-resolved transcriptomics through metabolic labeling and combinatorial indexing

Assessing RNA velocity stability across synthetic replicates using count splitting

Quantitative profiling of pseudouridylation dynamics in native RNAs with nanopore sequencing

L-RAPiT: A Cloud-Based Computing Pipeline for the Analysis of Long-Read RNA Sequencing Data