Comprehensive evaluation of RNA-seq quantification methods for linearity

BMC Bioinformatics

Research

Ying-Wooi Wan¹  Zhandong Liu²  Haijing Jin³ 

[1] Department of Molecular and Human Genetics, Baylor College of Medicine, One Baylor Plaza, 77030, Houston, TX, USA;Department of Pediatrics-Neurology, Jan and Dan Duncan Neurological Research Institute, Baylor College of Medicine, 1250 Moursund St., Suite 1325, 77030, Houston, TX, USA;Graduate Program in Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, One Baylor Plaza, 77030, Houston, TX, USA;

关键词: RNA-seq;    Deconvolution;    Linearity;   

DOI  :  10.1186/s12859-017-1526-y

来源: Springer

【 摘 要 】

BackgroundDeconvolution is a mathematical process of resolving an observed function into its constituent elements. In the field of biomedical research, deconvolution analysis is applied to obtain single cell-type or tissue specific signatures from a mixed signal and most of them follow the linearity assumption. Although recent development of next generation sequencing technology suggests RNA-seq as a fast and accurate method for obtaining transcriptomic profiles, few studies have been conducted to investigate best RNA-seq quantification methods that yield the optimum linear space for deconvolution analysis.ResultsUsing a benchmark RNA-seq dataset, we investigated the linearity of abundance estimated from seven most popular RNA-seq quantification methods both at the gene and isoform levels. Linearity is evaluated through parameter estimation, concordance analysis and residual analysis based on a multiple linear regression model. Results show that count data gives poor parameter estimations, large intercepts and high inter-sample variability; while TPM value from Kallisto and Salmon shows high linearity in all analyses.ConclusionsSalmon and Kallisto TPM data gives the best fit to the linear model studied. This suggests that TPM values estimated from Salmon and Kallisto are the ideal RNA-seq measurements for deconvolution studies.

【 授权许可】

CC BY   
© The Author(s) 2017

【 图 表 】

Fig. 4

Fig. 2

Scheme 1

12951_2016_223_Article_IEq1.gif

12951_2017_255_Article_IEq50.gif

12951_2017_255_Article_IEq49.gif

12936_2017_1904_Article_IEq1.gif

Fig. 1