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.

Accurate prediction of gene expression by integration of DNA sequence statistics with detailed modeling of transcription regulation

Effective gene expression prediction from sequence by integrating long-range interactions

Understanding Distal Transcriptional Regulation from Sequence Motif, Network Inference and Interactome Perspectives

Understanding Transcriptional Regulation by Integrative Analysis of Transcription Factor Binding Data

Integrating genome sequence and structural data for statistical learning to predict transcription factor binding sites

Efficiency and versatility of distal multisite transcription regulation

Posttranscriptional regulatory networks: from expression profiling to integrative analysis of mRNA and microRNA data.

Fundamentals for predicting transcriptional regulations from DNA sequence patterns

Transcription factor binding site prediction with multivariate gene expression data

Learning the Regulatory Code of Gene Expression

Deciphering eukaryotic gene-regulatory logic with 100 million random promoters

Systems Biophysics of Gene Expression

Genome-wide in silico prediction of gene expression

Quantifying the relationship between co-expression, co-regulation and gene function

A unified-field theory of genome organization and gene regulation

EPInformer: A Scalable Deep Learning Framework for Gene Expression Prediction by Integrating Promoter-enhancer Sequences with Multimodal Epigenomic Data

Dynamics of gene expression and the regulatory inference problem

GENet: A Graph-Based Model Leveraging Histone Marks and Transcription Factors for Enhanced Gene Expression Prediction

Using sequence-specific chemical and structural properties of DNA to predict transcription factor binding sites

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

Integrative prediction of gene expression with chromatin accessibility and conformation data