Alexander Kanitz1, Foivos Gypas2, Andreas J Gruber3, Andreas R Gruber4, Georges Martin5, Mihaela Zavolan6. 1. Biozentrum, University of Basel and Swiss Institute of Bioinformatics, Basel, Switzerland. alexander.kanitz@unibas.ch. 2. Biozentrum, University of Basel and Swiss Institute of Bioinformatics, Basel, Switzerland. foivos.gypas@unibas.ch. 3. Biozentrum, University of Basel and Swiss Institute of Bioinformatics, Basel, Switzerland. aj.gruber@unibas.ch. 4. Biozentrum, University of Basel and Swiss Institute of Bioinformatics, Basel, Switzerland. andreas.gruber@unibas.ch. 5. Biozentrum, University of Basel and Swiss Institute of Bioinformatics, Basel, Switzerland. georges.martin@unibas.ch. 6. Biozentrum, University of Basel and Swiss Institute of Bioinformatics, Basel, Switzerland. mihaela.zavolan@unibas.ch.
Abstract
BACKGROUND: Understanding the regulation of gene expression, including transcription start site usage, alternative splicing, and polyadenylation, requires accurate quantification of expression levels down to the level of individual transcript isoforms. To comparatively evaluate the accuracy of the many methods that have been proposed for estimating transcript isoform abundance from RNA sequencing data, we have used both synthetic data as well as an independent experimental method for quantifying the abundance of transcript ends at the genome-wide level. RESULTS: We found that many tools have good accuracy and yield better estimates of gene-level expression compared to commonly used count-based approaches, but they vary widely in memory and runtime requirements. Nucleotide composition and intron/exon structure have comparatively little influence on the accuracy of expression estimates, which correlates most strongly with transcript/gene expression levels. To facilitate the reproduction and further extension of our study, we provide datasets, source code, and an online analysis tool on a companion website, where developers can upload expression estimates obtained with their own tool to compare them to those inferred by the methods assessed here. CONCLUSIONS: As many methods for quantifying isoform abundance with comparable accuracy are available, a user's choice will likely be determined by factors such as the memory and runtime requirements, as well as the availability of methods for downstream analyses. Sequencing-based methods to quantify the abundance of specific transcript regions could complement validation schemes based on synthetic data and quantitative PCR in future or ongoing assessments of RNA-seq analysis methods.
BACKGROUND: Understanding the regulation of gene expression, including transcription start site usage, alternative splicing, and polyadenylation, requires accurate quantification of expression levels down to the level of individual transcript isoforms. To comparatively evaluate the accuracy of the many methods that have been proposed for estimating transcript isoform abundance from RNA sequencing data, we have used both synthetic data as well as an independent experimental method for quantifying the abundance of transcript ends at the genome-wide level. RESULTS: We found that many tools have good accuracy and yield better estimates of gene-level expression compared to commonly used count-based approaches, but they vary widely in memory and runtime requirements. Nucleotide composition and intron/exon structure have comparatively little influence on the accuracy of expression estimates, which correlates most strongly with transcript/gene expression levels. To facilitate the reproduction and further extension of our study, we provide datasets, source code, and an online analysis tool on a companion website, where developers can upload expression estimates obtained with their own tool to compare them to those inferred by the methods assessed here. CONCLUSIONS: As many methods for quantifying isoform abundance with comparable accuracy are available, a user's choice will likely be determined by factors such as the memory and runtime requirements, as well as the availability of methods for downstream analyses. Sequencing-based methods to quantify the abundance of specific transcript regions could complement validation schemes based on synthetic data and quantitative PCR in future or ongoing assessments of RNA-seq analysis methods.
Authors: Mihaela Zavolan; Shinji Kondo; Christian Schonbach; Jun Adachi; David A Hume; Yoshihide Hayashizaki; Terry Gaasterland Journal: Genome Res Date: 2003-06 Impact factor: 9.043
Authors: Ning Leng; John A Dawson; James A Thomson; Victor Ruotti; Anna I Rissman; Bart M G Smits; Jill D Haag; Michael N Gould; Ron M Stewart; Christina Kendziorski Journal: Bioinformatics Date: 2013-02-21 Impact factor: 6.937
Authors: Josep F Abril; Pär G Engström; Felix Kokocinski; Tamara Steijger; Tim J Hubbard; Roderic Guigó; Jennifer Harrow; Paul Bertone Journal: Nat Methods Date: 2013-11-03 Impact factor: 28.547
Authors: Narayanan Raghupathy; Kwangbom Choi; Matthew J Vincent; Glen L Beane; Keith S Sheppard; Steven C Munger; Ron Korstanje; Fernando Pardo-Manual de Villena; Gary A Churchill Journal: Bioinformatics Date: 2018-07-01 Impact factor: 6.937
Authors: Arfa Mehmood; Asta Laiho; Mikko S Venäläinen; Aidan J McGlinchey; Ning Wang; Laura L Elo Journal: Brief Bioinform Date: 2020-12-01 Impact factor: 11.622
Authors: Becky C Carlyle; Robert R Kitchen; Jing Zhang; Rashaun S Wilson; Tukiet T Lam; Joel S Rozowsky; Kenneth R Williams; Nenad Sestan; Mark B Gerstein; Angus C Nairn Journal: J Proteome Res Date: 2018-09-06 Impact factor: 4.466
Authors: Sheryl E Arambula; Scott M Belcher; Antonio Planchart; Stephen D Turner; Heather B Patisaul Journal: Endocrinology Date: 2016-08-29 Impact factor: 4.736