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

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 disrupti...

Full description

Bibliographic Details
Main Authors: Amit Blumberg, Yixin Zhao, Yi-Fei Huang, Noah Dukler, Edward J. Rice, Alexandra G. Chivu, Katie Krumholz, Charles G. Danko, Adam Siepel
Format: Article
Language:English
Published: BMC 2021-02-01
Series:BMC Biology
Subjects:
Online Access:https://doi.org/10.1186/s12915-021-00949-x
_version_ 1818411874439397376
author Amit Blumberg
Yixin Zhao
Yi-Fei Huang
Noah Dukler
Edward J. Rice
Alexandra G. Chivu
Katie Krumholz
Charles G. Danko
Adam Siepel
author_facet Amit Blumberg
Yixin Zhao
Yi-Fei Huang
Noah Dukler
Edward J. Rice
Alexandra G. Chivu
Katie Krumholz
Charles G. Danko
Adam Siepel
author_sort Amit Blumberg
collection DOAJ
description 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.
first_indexed 2024-12-14T10:38:21Z
format Article
id doaj.art-402f5011b83a418da6494e4ad564b208
institution Directory Open Access Journal
issn 1741-7007
language English
last_indexed 2024-12-14T10:38:21Z
publishDate 2021-02-01
publisher BMC
record_format Article
series BMC Biology
spelling doaj.art-402f5011b83a418da6494e4ad564b2082022-12-21T23:05:48ZengBMCBMC Biology1741-70072021-02-0119111710.1186/s12915-021-00949-xCharacterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq dataAmit Blumberg0Yixin Zhao1Yi-Fei Huang2Noah Dukler3Edward J. Rice4Alexandra G. Chivu5Katie Krumholz6Charles G. Danko7Adam Siepel8Simons Center for Quantitative Biology, Cold Spring Harbor LaboratorySimons Center for Quantitative Biology, Cold Spring Harbor LaboratorySimons Center for Quantitative Biology, Cold Spring Harbor LaboratorySimons Center for Quantitative Biology, Cold Spring Harbor LaboratoryBaker Institute for Animal Health, College of Veterinary Medicine, Cornell UniversityBaker Institute for Animal Health, College of Veterinary Medicine, Cornell UniversitySimons Center for Quantitative Biology, Cold Spring Harbor LaboratoryBaker Institute for Animal Health, College of Veterinary Medicine, Cornell UniversitySimons Center for Quantitative Biology, Cold Spring Harbor LaboratoryAbstract 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.https://doi.org/10.1186/s12915-021-00949-xRNA half-lifeRNA splicingEpigenomicsPRO-seqStructural equation modeling
spellingShingle Amit Blumberg
Yixin Zhao
Yi-Fei Huang
Noah Dukler
Edward J. Rice
Alexandra G. Chivu
Katie Krumholz
Charles G. Danko
Adam Siepel
Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data
BMC Biology
RNA half-life
RNA splicing
Epigenomics
PRO-seq
Structural equation modeling
title Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data
title_full Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data
title_fullStr Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data
title_full_unstemmed Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data
title_short Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data
title_sort characterizing rna stability genome wide through combined analysis of pro seq and rna seq data
topic RNA half-life
RNA splicing
Epigenomics
PRO-seq
Structural equation modeling
url https://doi.org/10.1186/s12915-021-00949-x
work_keys_str_mv AT amitblumberg characterizingrnastabilitygenomewidethroughcombinedanalysisofproseqandrnaseqdata
AT yixinzhao characterizingrnastabilitygenomewidethroughcombinedanalysisofproseqandrnaseqdata
AT yifeihuang characterizingrnastabilitygenomewidethroughcombinedanalysisofproseqandrnaseqdata
AT noahdukler characterizingrnastabilitygenomewidethroughcombinedanalysisofproseqandrnaseqdata
AT edwardjrice characterizingrnastabilitygenomewidethroughcombinedanalysisofproseqandrnaseqdata
AT alexandragchivu characterizingrnastabilitygenomewidethroughcombinedanalysisofproseqandrnaseqdata
AT katiekrumholz characterizingrnastabilitygenomewidethroughcombinedanalysisofproseqandrnaseqdata
AT charlesgdanko characterizingrnastabilitygenomewidethroughcombinedanalysisofproseqandrnaseqdata
AT adamsiepel characterizingrnastabilitygenomewidethroughcombinedanalysisofproseqandrnaseqdata