Abstract:Abstract Background RNA sequencing is currently the method of choice for genome-wide profiling of gene expression. A popular approach to quantify expression levels of genes from RNA-seq data is to map reads to a reference genome and then count mapped reads to each gene. Gene annotation data, which include chromosomal coordinates of exons for tens of thousands of genes, are required for this quantification process. There are several major sources of gene annotations that can be used for quantification, such as Ensembl and RefSeq databases. However, there is very little understanding of the effect that the choice of annotation has on the accuracy of gene expression quantification in an RNA-seq analysis. Results In this paper, we present results from our comparison of Ensembl and RefSeq human annotations on their impact on gene expression quantification using a benchmark RNA-seq dataset generated by the SEQC consortium. We show that the use of RefSeq gene annotation models led to better quantification accuracy, based on the correlation with ground truths including expression data from >800 real-time PCR validated genes, known titration ratios of gene expression and microarray expression data. We also found that the recent expansion of the RefSeq annotation has led to a decrease in its annotation accuracy. Finally, we demonstrated that the RNA-seq quantification differences observed between different annotations were not affected by the use of different normalization methods. Conclusion In conclusion, our study found that the use of the conservative RefSeq gene annotation yields better RNA-seq quantification results than the more comprehensive Ensembl annotation. We also found that, surprisingly, the recent expansion of the RefSeq database, which was primarily driven by the incorporation of sequencing data into the gene annotation process, resulted in a reduction in the accuracy of RNA-seq quantification.

Information transduction capacity reduces the uncertainties in annotation-free isoform discovery and quantification.

Improving the Diversity of Captured Full-Length Isoforms Using a Normalized Single-Molecule RNA-sequencing Method

Comprehensive Assessment of Isoform Detection Methods for Third-Generation Sequencing Data

Evaluation and comparison of computational tools for RNA-seq isoform quantification

Isoform Abundance Inference Provides a More Accurate Estimation of Gene Expression Levels in RNA-seq.

Enhancing novel isoform discovery: leveraging nanopore long-read sequencing and machine learning approaches

Accurate isoform quantification by joint short- and long-read RNA-sequencing

Statistical Inferences for Isoform Expression in RNA-Seq

Joint Estimation of Isoform Expression and Isoform-Specific Read Distribution Using Multisample RNA-Seq Data.

Network-Based Isoform Quantification with RNA-Seq Data for Cancer Transcriptome Analysis

Comprehensive assessment of mRNA isoform detection methods for long-read sequencing data

Alternative splicing: Human disease and quantitative analysis from high-throughput sequencing

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

PennSeq: Accurate Isoform-Specific Gene Expression Quantification in RNA-Seq by Modeling Non-Uniform Read Distribution

BIISQ: Bayesian nonparametric discovery of Isoforms and Individual Specific Quantification

RSQ: a Statistical Method for Quantification of Isoform-Specific Structurome Using Transcriptome-Wide Structural Profiling Data

TIF-Seq2 disentangles overlapping isoforms in complex human transcriptomes

Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences

Impact of gene annotation choice on the quantification of RNA-seq data

Statistical Modeling of RNA-Seq Data

MetaTX: deciphering the distribution of mRNA-related features in the presence of isoform ambiguity, with applications in epitranscriptome analysis