Abstract:BACKGROUND:The high-throughput sequencing technology, RNA-Seq, has been widely used to quantify gene and isoform expression in the study of transcriptome in recent years. Accurate expression measurement from the millions or billions of short generated reads is obstructed by difficulties. One is ambiguous mapping of reads to reference transcriptome caused by alternative splicing. This increases the uncertainty in estimating isoform expression. The other is non-uniformity of read distribution along the reference transcriptome due to positional, sequencing, mappability and other undiscovered sources of biases. This violates the uniform assumption of read distribution for many expression calculation approaches, such as the direct RPKM calculation and Poisson-based models. Many methods have been proposed to address these difficulties. Some approaches employ latent variable models to discover the underlying pattern of read sequencing. However, most of these methods make bias correction based on surrounding sequence contents and share the bias models by all genes. They therefore cannot estimate gene- and isoform-specific biases as revealed by recent studies.RESULTS:We propose a latent variable model, NLDMseq, to estimate gene and isoform expression. Our method adopts latent variables to model the unknown isoforms, from which reads originate, and the underlying percentage of multiple spliced variants. The isoform- and exon-specific read sequencing biases are modeled to account for the non-uniformity of read distribution, and are identified by utilizing the replicate information of multiple lanes of a single library run. We employ simulation and real data to verify the performance of our method in terms of accuracy in the calculation of gene and isoform expression. Results show that NLDMseq obtains competitive gene and isoform expression compared to popular alternatives. Finally, the proposed method is applied to the detection of differential expression (DE) to show its usefulness in the downstream analysis.CONCLUSIONS:The proposed NLDMseq method provides an approach to accurately estimate gene and isoform expression from RNA-Seq data by modeling the isoform- and exon-specific read sequencing biases. It makes use of a latent variable model to discover the hidden pattern of read sequencing. We have shown that it works well in both simulations and real datasets, and has competitive performance compared to popular methods. The method has been implemented as a freely available software which can be found at https://github.com/PUGEA/NLDMseq.

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

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

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

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

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

Estimation of Isoform Expression in Rna-Seq Data Using A Hierarchical Bayesian Model

BIISQ: Bayesian nonparametric discovery of Isoforms and Individual Specific Quantification

Statistical Inferences for Isoform Expression in RNA-Seq

PennSeq: Accurate Isoform-Specific Gene Expression Quantification in RNA-Seq by Modeling Non-Uniform Read Distribution

Statistical modeling of isoform splicing dynamics from RNA-seq time series data

DEIsoM: a Hierarchical Bayesian Model for Identifying Differentially Expressed Isoforms Using Biological Replicates

Statistical Modeling of RNA-Seq Data

Improving RNA-Seq Expression Estimation by Modeling Isoform- and Exon-Specific Read Sequencing Rate

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

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

NSMAP: A method for spliced isoforms identification and quantification from RNA-Seq

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

LIQA: long-read isoform quantification and analysis

Robust estimation of isoform expression with RNA-Seq data

A mixed model approach for joint genetic analysis of alternatively spliced transcript isoforms using RNA-Seq data

QuickMIRSeq: a Pipeline for Quick and Accurate Quantification of Both Known Mirnas and Isomirs by Jointly Processing Multiple Samples from Microrna Sequencing.