Abstract: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.

Relative transcript quantification by Quantitative PCR: Roughly right or precisely wrong?

A statistical method for the conservative adjustment of false discovery rate (q-value)

Statistical significance of quantitative PCR

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

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

Quantification of mRNA using real-time RT-PCR

Models for transcript quantification from RNA-Seq

Benchmarking of RNA-sequencing analysis workflows using whole-transcriptome RT-qPCR expression data

Removal of between-run variation in a multi-plate qPCR experiment

Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays.

A PCR Primer Bank for Quantitative Gene Expression Analysis.

[Reverse transcriptase PCR (RT-PCR) and quantitative-competitive PCR (QC-PCR)]

Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples

Comprehensive algorithm for quantitative real-time polymerase chain reaction

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

Analyzing real-time PCR data by the comparative CT method

Identification and optimization of parameters for accurate quantification of RNA by RT-dPCR

Misuse of RPKM or TPM normalization when comparing across samples and sequencing protocols

Absolute quantification of real-time PCR data with stage signal difference analysis

Real‐time quantitative PCR: A tool for absolute and relative quantification

mRNA enrichment protocols determine the quantification characteristics of external RNA spike-in controls in RNA-Seq studies