Union Exon Based Approach for RNA-Seq Gene Quantification: To Be or Not to Be?

• Published in: PloS one 2015-11, Vol.10 (11), p.e0141910-e0141910
• Main Authors: Zhao, Shanrong, Xi, Li, Zhang, Baohong
• Format: Article
• Language: English
• Subjects: Accuracy; Algorithms; Annotations; Bioinformatics; Consortia; Data analysis; Exons; Exons (Molecular genetics); Gene expression; Gene loci; Genetic research; Genomes; Genomics; Isoforms; Methods; Ontology; Properties; R&D; Relative abundance; Research & development; Ribonucleic acid; RNA; RNA sequencing; Sequence Analysis, RNA - methods; Switches; Transcription; Transcription (Genetics)
• ISSN: 1932-6203; 1932-6203
• DOI: 10.1371/journal.pone.0141910

In recent years, RNA-seq is emerging as a powerful technology in estimation of gene and/or transcript expression, and RPKM (Reads Per Kilobase per Million reads) is widely used to represent the relative abundance of mRNAs for a gene. In general, the methods for gene quantification can be largely divided into two categories: transcript-based approach and 'union exon'-based approach. Transcript-based approach is intrinsically more difficult because different isoforms of the gene typically have a high proportion of genomic overlap. On the other hand, 'union exon'-based approach method is much simpler and thus widely used in RNA-seq gene quantification. Biologically, a gene is expressed in one or more transcript isoforms. Therefore, transcript-based approach is logistically more meaningful than 'union exon'-based approach. Despite the fact that gene quantification is a fundamental task in most RNA-seq studies, however, it remains unclear whether 'union exon'-based approach for RNA-seq gene quantification is a good practice or not. In this paper, we carried out a side-by-side comparison of 'union exon'-based approach and transcript-based method in RNA-seq gene quantification. It was found that the gene expression levels are significantly underestimated by 'union exon'-based approach, and the average of RPKM from 'union exons'-based method is less than 50% of the mean expression obtained from transcript-based approach. The difference between the two approaches is primarily affected by the number of transcripts in a gene. We performed differential analysis at both gene and transcript levels, respectively, and found more insights, such as isoform switches, are gained from isoform differential analysis. The accuracy of isoform quantification would improve if the read coverage pattern and exon-exon spanning reads are taken into account and incorporated into EM (Expectation Maximization) algorithm. Our investigation discourages the use of 'union exons'-based approach in gene quantification despite its simplicity.