Abstract:Motivation: Assigning RNA-seq reads to their transcript of origin is a fundamental task in transcript expression estimation. Where ambiguities in assignments exist due to transcripts sharing sequence, e.g. alternative isoforms or alleles, the problem can be solved through probabilistic inference. Bayesian methods have been shown to provide accurate transcript abundance estimates compared to competing methods. However, exact Bayesian inference is intractable and approximate methods such as Markov chain Monte Carlo (MCMC) and Variational Bayes (VB) are typically used. While providing a high degree of accuracy and modelling flexibility, standard implementations can be prohibitively slow for large datasets and complex transcriptome annotations. Results: We propose a novel approximate inference scheme based on VB and apply it to an existing model of transcript expression inference from RNA-seq data. Recent advances in VB algorithmics are used to improve the convergence of the algorithm beyond the standard Variational Bayes Expectation Maximisation (VBEM) algorithm. We apply our algorithm to simulated and biological datasets, demonstrating a significant increase in speed with only very small loss in accuracy of expression level estimation. We carry out a comparative study against seven popular alternative methods and demonstrate that our new algorithm provides excellent accuracy and inter-replicate consistency while remaining competitive in computation time. Availability: The methods were implemented in R and C++, and are available as part of the BitSeq project at \url{<a class="link-external link-https" href="https://github.com/BitSeq" rel="external noopener nofollow">this https URL</a>}. The method is also available through the BitSeq Bioconductor package. The source code to reproduce all simulation results can be accessed via \url{<a class="link-external link-https" href="https://github.com/BitSeq/BitSeqVB_benchmarking" rel="external noopener nofollow">this https URL</a>}.

Identifying differentially expressed transcripts from RNA-seq data with biological variation

A Bayesian model selection approach for identifying differentially expressed transcripts from RNA-Seq data

Fast and accurate approximate inference of transcript expression from RNA-seq data

Bayesian estimation of Differential Transcript Usage from RNA-seq data

Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation

Detecting Differential Expression from RNA-seq Data with Expression Measurement Uncertainty

EBSeq-HMM: a Bayesian Approach for Identifying Gene-Expression Changes in Ordered RNA-seq Experiments

DifferentialRegulation: a Bayesian hierarchical approach to identify differentially regulated genes

Differential expression analysis for paired RNA-seq data

Rseqdiff: Detecting Differential Isoform Expression from RNA-Seq Data Using Hierarchical Likelihood Ratio Test.

Degps is a Powerful Tool for Detecting Differential Expression in RNA-sequencing Studies

Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences

A comparison of methods for differential expression analysis of RNA-seq data

Non-parametric Bayesian modelling of digital gene expression data

Identifying Differentially Spliced Genes from Two Groups of Rna-Seq Samples

A deep generative model for single-cell RNA sequencing with application to detecting differentially expressed genes

Detecting differential usage of exons from RNA-seq data

BADER: Bayesian analysis of differential expression in RNA sequencing data

3t-seq: automatic gene expression analysis of single-copy genes, transposable elements, and tRNAs from RNA-seq data

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

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