Abstract:Next-generation RNA sequencing (RNA-seq) technology has been widely used to assess full-length RNA isoform abundance in a high-throughput manner. RNA-seq data offer insight into gene expression levels and transcriptome structures, enabling us to better understand the regulation of gene expression and fundamental biological processes. Accurate isoform quantification from RNA-seq data is challenging due to the information loss in sequencing experiments. A recent accumulation of multiple RNA-seq data sets from the same tissue or cell type provides new opportunities to improve the accuracy of isoform quantification. However, existing statistical or computational methods for multiple RNA-seq samples either pool the samples into one sample or assign equal weights to the samples when estimating isoform abundance. These methods ignore the possible heterogeneity in the quality of different samples and could result in biased and unrobust estimates. In this article, we develop a method, which we call "joint modeling of multiple RNA-seq samples for accurate isoform quantification" (MSIQ), for more accurate and robust isoform quantification by integrating multiple RNA-seq samples under a Bayesian framework. Our method aims to (1) identify a consistent group of samples with homogeneous quality and (2) improve isoform quantification accuracy by jointly modeling multiple RNA-seq samples by allowing for higher weights on the consistent group. We show that MSIQ provides a consistent estimator of isoform abundance, and we demonstrate the accuracy and effectiveness of MSIQ compared with alternative methods through simulation studies on <i>D. melanogaster</i> genes. We justify MSIQ's advantages over existing approaches via application studies on real RNA-seq data from human embryonic stem cells, brain tissues, and the HepG2 immortalized cell line. We also perform a comprehensive analysis of how the isoform quantification accuracy would be affected by RNA-seq sample heterogeneity and different experimental protocols.

MAJIQ V3 offers improvements in accuracy, performance, and usability for splicing analysis from RNA sequencing

Contrasting and combining transcriptome complexity captured by short and long RNA sequencing reads

MntJULiP and Jutils: Differential splicing analysis of RNA-seq data with covariates

A Deep Dive into Statistical Modeling of RNA Splicing QTLs Reveals New Variants that Explain Neurodegenerative Disease

Comprehensive Assessment of Isoform Detection Methods for Third-Generation Sequencing Data

ASJA: A Program for Assembling Splice Junctions Analysis

RNA-Seq Analysis of Differential Splice Junction Usage and Intron Retentions by DEXSeq

recount3: summaries and queries for large-scale RNA-seq expression and splicing

Mapping RNA splicing variations in clinically accessible and nonaccessible tissues to facilitate Mendelian disease diagnosis using RNA-seq

Observations on Novel Splice Junctions from RNA Sequencing Data

LIQA: long-read isoform quantification and analysis

Machine learning-optimized targeted detection of alternative splicing

Annotation-free quantification of RNA splicing using LeafCutter

Detection and Visualization of Differential Splicing in RNA-Seq Data with JunctionSeq

Selecting Differential Splicing Methods: Practical Considerations

SpliceWiz: interactive analysis and visualization of alternative splicing in R

Performance evaluation of computational methods for splice-disrupting variants and improving the performance using the machine learning-based framework

junctionCounts: comprehensive alternative splicing analysis and prediction of isoform-level impacts to the coding sequence

MSIQ: Joint modeling of multiple RNA-seq samples for accurate isoform quantification

In Situ Imaging of mRNA Splicing Variants by SpliceRCA.

An Empirical Bayes Change-Point Model for Identifying 3′ and 5′ Alternative Splicing by Next-Generation RNA Sequencing