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

Blumberg, Amit, Zhao, Yixin, Huang, Yi-Fei, Dukler, Noah, Rice, Edward J, Chivu, Alexandra G, Krumholz, Katie, Danko, Charles G, Siepel, Adam (February 2021) Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data. BMC Biology, 19 (1). p. 30. ISSN 1741-7007

[thumbnail of 2021.Blumberg.RNAStability.pdf] PDF
2021.Blumberg.RNAStability.pdf

Download (6MB)

Abstract

BACKGROUND: The concentrations of distinct types of RNA in cells result from a dynamic equilibrium between RNA synthesis and decay. Despite the critical importance of RNA decay rates, current approaches for measuring them are generally labor-intensive, limited in sensitivity, and/or disruptive to normal cellular processes. Here, we introduce a simple method for estimating relative RNA half-lives that is based on two standard and widely available high-throughput assays: Precision Run-On sequencing (PRO-seq) and RNA sequencing (RNA-seq). RESULTS: Our method treats PRO-seq as a measure of transcription rate and RNA-seq as a measure of RNA concentration, and estimates the rate of RNA decay required for a steady-state equilibrium. We show that this approach can be used to assay relative RNA half-lives genome-wide, with good accuracy and sensitivity for both coding and noncoding transcription units. Using a structural equation model (SEM), we test several features of transcription units, nearby DNA sequences, and nearby epigenomic marks for associations with RNA stability after controlling for their effects on transcription. We find that RNA splicing-related features are positively correlated with RNA stability, whereas features related to miRNA binding and DNA methylation are negatively correlated with RNA stability. Furthermore, we find that a measure based on U1 binding and polyadenylation sites distinguishes between unstable noncoding and stable coding transcripts but is not predictive of relative stability within the mRNA or lincRNA classes. We also identify several histone modifications that are associated with RNA stability. CONCLUSION: We introduce an approach for estimating the relative half-lives of individual RNAs. Together, our estimation method and systematic analysis shed light on the pervasive impacts of RNA stability on cellular RNA concentrations.

Item Type: Paper
Subjects: Investigative techniques and equipment
bioinformatics > genomics and proteomics > genetics & nucleic acid processing > DNA, RNA structure, function, modification > RNA expression
Investigative techniques and equipment > assays
Investigative techniques and equipment > PRO-seq
Investigative techniques and equipment > assays > RNA-seq
CSHL Authors:
Communities: CSHL labs > Siepel lab
SWORD Depositor: CSHL Elements
Depositing User: CSHL Elements
Date: 15 February 2021
Date Deposited: 04 May 2021 15:01
Last Modified: 23 Jan 2024 18:31
PMCID: PMC7885420
Related URLs:
URI: https://repository.cshl.edu/id/eprint/39980

Actions (login required)

Administrator's edit/view item Administrator's edit/view item