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

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. 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 SEquencing Quality Control (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.

IDENTIFICATION OF SHORTENED 3′ UNTRANSLATED REGIONS FROM EXPRESSION ARRAYS

Cancer cells have been recently shown to express high level of short 3′UTR isoforms that can escape miRNA-mediated regulation. We present here a computational procedure for systematically identifying shortened 3′UTRs by Affymetrix 3′ microarrays. The advantage of this technology compared to more recent and promising ones such as exon arrays and RNA-Seq is that, giving the relatively small cost, already existing datasets in public databases include a considerably higher number of experiments. Moreover, the design of Affymetrix Gene Chips is well-suited for 3′UTR analysis of a large number of genes. Initially, Affymetrix individual probes are regrouped into customized probesets mapping specifically the CDS or the 3′UTR of the transcript, according to RefSeq annotation. Then, candidate 3′UTR shortening events are identified by statistical differential expression analysis of customized probesets in different biological conditions. The procedure has been applied to expression data from two ovarian adenocarcinoma datasets. Selected gene sets are significantly enriched for annotated splice variant genes as well as genes involved in estrogen dependent cancer mechanisms, confirming the validity of the proposed procedure.