RNA sequencing, or RNA-seq for short, is widely applied for extracting gene and transcript expression from biological samples. Given numerous quantification pipelines for RNA-seq data, one fundamental challenge is to identify a pipeline that can produce the most accurate gene and/or transcript expression. Exploring all pipelines requires extensive computational resources, so we propose to use a subsampling approach that can speed up the pipeline evaluation and selection process. We applied our approach to one simulated and two real RNA-seq datasets and found that expression estimates derived from subsampled data are close surrogates for those derived from original data. In addition, the ranking of quantification pipelines based on the subsampled data was highly concordant with that based on the original data. Therefore, we conclude that subsampling is a valid approach to facilitating efficient quantification pipeline selection using RNA-seq data.
['Po-Yen Wu', 'May D. Wang']
The selection of quantification pipelines for Illumina RNA-seq data using a subsampling approach