Quantifying alternative splicing from paired-end RNA-sequencing data. RNA-sequencing has revolutionized biomedical research and, in particular, our ability to study gene alternative splicing. The problem has important implications for human health, as alternative splicing may be involved in malfunctions at the cellular level and multiple diseases. However, the high-dimensional nature of the data and the existence of experimental biases pose serious data analysis challenges. We find that the standard data summaries used to study alternative splicing are severely limited, as they ignore a substantial amount of valuable information. Current data analysis methods are based on such summaries and are hence suboptimal. Further, they have limited flexibility in accounting for technical biases. We propose novel data summaries and a Bayesian modeling framework that overcome these limitations and determine biases in a nonparametric, highly flexible manner. These summaries adapt naturally to the rapid improvements in sequencing technology. We provide efficient point estimates and uncertainty assessments. The approach allows to study alternative splicing patterns for individual samples and can also be the basis for downstream analyses. We found a severalfold improvement in estimation mean square error compared popular approaches in simulations, and substantially higher consistency between replicates in experimental data. Our findings indicate the need for adjusting the routine summarization and analysis of alternative splicing RNA-seq studies. We provide a software implementation in the R package casper.
Keywords for this software
References in zbMATH (referenced in 5 articles )
Showing results 1 to 5 of 5.
- Li, Wei Vivian; Zhao, Anqi; Zhang, Shihua; Li, Jingyi Jessica: MSIQ: joint modeling of multiple RNA-seq samples for accurate isoform quantification (2018)
- Papastamoulis, Panagiotis; Rattray, Magnus: Bayesian estimation of differential transcript usage from RNA-seq data (2017)
- Lin, Zhixiang; Li, Mingfeng; Sestan, Nenad; Zhao, Hongyu: A Markov random field-based approach for joint estimation of differentially expressed genes in mouse transcriptome data (2016)
- Rossell, David; Attolini, Camille Stephan-Otto; Kroiss, Manuel; Stöcker, Almond: Corrigendum: quantifying alternative splicing from paired-end RNA-seq data (2015)
- Rossell, David; Attolini, Camille Stephan-Otto; Kroiss, Manuel; Stöcker, Almond: Quantifying alternative splicing from paired-end RNA-sequencing data (2014)