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.

RSQ: a Statistical Method for Quantification of Isoform-Specific Structurome Using Transcriptome-Wide Structural Profiling Data

Toward a Next-Generation Atlas of RNA Secondary Structure

Improving the Diversity of Captured Full-Length Isoforms Using a Normalized Single-Molecule RNA-sequencing Method

Statistical Inferences for Isoform Expression in RNA-Seq

Differential Analysis of RNA Structure Probing Experiments at Nucleotide Resolution: Uncovering Regulatory Functions of RNA Structure

RNA Regulations and Functions Decoded by Transcriptome-wide RNA Structure Probing

Accurate isoform quantification by joint short- and long-read RNA-sequencing

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

RNAsmc: A Integrated Tool for Comparing RNA Secondary Structure and Evaluating Allosteric Effects

Network-Based Isoform Quantification with RNA-Seq Data for Cancer Transcriptome Analysis

Joint Estimation of Isoform Expression and Isoform-Specific Read Distribution Using Multisample RNA-Seq Data.

Comprehensive assessment of mRNA isoform detection methods for long-read sequencing data

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

RIC-seq for global in situ profiling of RNA–RNA spatial interactions

Single-cell RNA Cap and Tail Sequencing (Scrcat-Seq) Reveals Subtype-Specific Isoforms Differing in Transcript Demarcation

IQSeq: Integrated Isoform Quantification Analysis Based on Next-Generation Sequencing.

Evaluation and comparison of computational tools for RNA-seq isoform quantification

AStruct: detection of allele-specific RNA secondary structure in structuromic probing data

Isoform Abundance Inference Provides a More Accurate Estimation of Gene Expression Levels in RNA-seq.

STRUCT: a statistical approach to identify RNA secondary structures from raw sequencing data, bypassing multiple sequence alignment

In vivo genome-wide profiling of RNA secondary structure reveals novel regulatory features