Towards reliable isoform quantification using RNA-SEQ data

BMC Bioinformatics

Proceedings

Steffen Heber¹  Brian E Howard¹ 

[1] Bioinformatics Research Center, North Carolina State University, 27606, Raleigh, NC, USA;

关键词: Design Matrix;    Read Length;    Differential Splice;    Good Linear Unbiased Estimator;    Good Linear Unbiased Estimator;   

DOI  :  10.1186/1471-2105-11-S3-S6

来源: Springer

【 摘 要 】

BackgroundIn eukaryotes, alternative splicing often generates multiple splice variants from a single gene. Here weexplore the use of RNA sequencing (RNA-Seq) datasets to address the isoform quantification problem. Given a set of known splice variants, the goal is to estimate the relative abundance of the individual variants.MethodsOur method employs a linear models framework to estimate the ratios of known isoforms in a sample. A key feature of our method is that it takes into account the non-uniformity of RNA-Seq read positions along the targeted transcripts.ResultsPreliminary tests indicate that the model performs well on both simulated and real data. In two publicly available RNA-Seq datasets, we identified several alternatively-spliced genes with switch-like, on/off expression properties, as well as a number of other genes that varied more subtly in isoform expression. In many cases, genes exhibiting differential expression of alternatively spliced transcripts were not differentially expressed at the gene level.ConclusionsGiven that changes in isoform expression level frequently involve a continuum of isoform ratios, rather than all-or-nothing expression, and that they are often independent of general gene expression changes, we anticipate that our research will contribute to revealing a so far uninvestigated layer of the transcriptome. We believe that, in the future, researchers will prioritize genes for functional analysis based not only on observed changes in gene expression levels, but also on changes in alternative splicing.

【 授权许可】

CC BY   
© Howard and Heber; licensee BioMed Central Ltd. 2010