Abstract:Abstract Nascent RNA sequencing protocols, such as PRO-seq and NET-seq, are now widely used in the study of eukaryotic transcription, and these experimental techniques have given rise to a variety of statistical and machine-learning methods for data analysis. These computational methods, however, are generally designed to address specialized signal-processing or prediction tasks, rather than directly describing the dynamics of RNA polymerases as they move along the DNA template. Here, I introduce a general probabilistic model that describes the kinetics of transcription initiation, elongation, pause release, and termination, as well as the generation of sequencing read counts. I show that this generative model enables estimation of separate pause-release rates, termination rates, and the initiation/elongation rate ratio up to a proportionality constant. Furthermore, if applied to time-course data in a nonequilibrium setting, the model can be used to estimate elongation rates. This model leads naturally to likelihood ratio tests for differences between genes, conditions, or species in various rates of interest. If read counts are assumed to be Poisson-distributed, convenient, closed-form solutions are available for both parameter estimates and likelihood-ratio-test statistics. Straightforward extensions of the model accommodate uncertainty in the pause site and steric hindrance of initiation by paused polymerases. Additional extensions address Bayesian inference under the Poisson model and a generalized linear model that can be used to discover genomic features associated with rates of elongation. Finally, I address technicalities concerning estimation of library size, normalization and sequencing replicates. Altogether, this modeling framework enables a unified treatment of many common tasks in the analysis of nascent RNA sequencing data.

Measuring Nascent Transcripts by Nascent-seq

Deciphering transcriptional dynamics in vivo by counting nascent RNA molecules

Comparative analysis of nascent RNA sequencing methods and their applications in studies of cotranscriptional splicing dynamics

Enhancing Transcriptome Mapping with Rapid PRO-seq Profiling of Nascent RNA

Transient transcriptome sequencing: computational pipeline to quantify genome-wide RNA kinetic parameters and transcriptional enhancer activity

Capturing nascent transcripts in single cells

Protocol for Quantitative Analysis of RNA 3′-End Processing Induced by Disassociated Subunits Using Chromatin-Associated RNA-seq Data

A Unified Probabilistic Modeling Framework for Eukaryotic Transcription Based on Nascent RNA Sequencing Data

Nascent RNA kinetics with complex promoter architecture: Analytic results and parameter inference

NanoTrans: an integrated computational framework for comprehensive transcriptome analysis with Nanopore direct RNA sequencing

TransRef Enables Accurate Transcriptome Assembly by Redefining Accurate Neo-Splicing Graphs.

Dissecting the nascent human transcriptome by analysing the RNA content of transcription factories

Landscape of Transcription Termination in Arabidopsis Revealed by Single-Molecule Nascent RNA Sequencing.

Single-cell nascent RNA sequencing unveils coordinated global transcription

Deep learning and direct sequencing of labeled RNA captures transcriptome dynamics

Quantification of mRNA Expression Using Single-Molecule Nanopore Sensing

Quantitative mapping of the mammalian epitranscriptome

Enhancing transcriptome expression quantification through accurate assignment of long RNA sequencing reads with TranSigner

Illuminating Enhancer Transcription at Nucleotide Resolution with Native Elongating Transcript Sequencing (NET-Seq)

RNA-eXpress annotates novel transcript features in RNA-seq data

Ab initio construction of a eukaryotic transcriptome by massively parallel mRNA sequencing