Abstract:Abstract Motivation High-throughput RNA sequencing has revolutionized the scope and depth of transcriptome analysis. Accurate reconstruction of a phenotype-specific transcriptome is challenging due to the noise and variability of RNA-seq data. This requires computational identification of transcripts from multiple samples of the same phenotype, given the underlying consensus transcript structure. Results We present a Bayesian method, integrated assembly of phenotype-specific transcripts (IntAPT), that identifies phenotype-specific isoforms from multiple RNA-seq profiles. IntAPT features a novel two-layer Bayesian model to capture the presence of isoforms at the group layer and to quantify the abundance of isoforms at the sample layer. A spike-and-slab prior is used to model the isoform expression and to enforce the sparsity of expressed isoforms. Dependencies between the existence of isoforms and their expression are modeled explicitly to facilitate parameter estimation. Model parameters are estimated iteratively using Gibbs sampling to infer the joint posterior distribution, from which the presence and abundance of isoforms can reliably be determined. Studies using both simulations and real datasets show that IntAPT consistently outperforms existing methods for the IntAPT. Experimental results demonstrate that, despite sequencing errors, IntAPT exhibits a robust performance among multiple samples, resulting in notably improved identification of expressed isoforms of low abundance. Availability and implementation The IntAPT package is available at http://github.com/henryxushi/IntAPT. Supplementary information Supplementary data are available at Bioinformatics online.

Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown

Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks

Analysis of RNA-Seq Data Using TopHat and Cufflinks

Transcriptome assembly from long-read RNA-seq alignments with StringTie2

StringTie enables improved reconstruction of a transcriptome from RNA-seq reads

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

Enhancing transcriptome expression quantification through accurate assignment of long RNA sequencing reads with TranSigner

Agtar: A Novel Approach For Transcriptome Assembly And Abundance Estimation Using An Adapted Genetic Algorithm From Rna-Seq Data

UTAP: User-friendly Transcriptome Analysis Pipeline

RNA-combine: a Toolkit for Comprehensive Analyses on Transcriptome Data from Different Sequencing Platforms.

RADPRE: A Computational Program for Identification of Differential Mrna Processing Including Alternative Polyadenylation

OneStopRNAseq: A Web Application for Comprehensive and Efficient Analyses of RNA-Seq Data

Comprehensive benchmark of differential transcript usage analysis for static and dynamic conditions

RNA-seq analysis is easy as 1-2-3 with limma, Glimma and edgeR

IntAPT: integrated assembly of phenotype-specific transcripts from multiple RNA-seq profiles

Transient transcriptome sequencing: computational pipeline to quantify genome-wide RNA kinetic parameters and transcriptional enhancer activity

Transcriptomics in the RNA-seq era

TransRef Enables Accurate Transcriptome Assembly by Redefining Accurate Neo-Splicing Graphs.

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

RNA-eXpress annotates novel transcript features in RNA-seq data

hppRNA-a Snakemake-based handy parameter-free pipeline for RNA-Seq analysis of numerous samples