Abstract:MOTIVATION: A critical task in high-throughput sequencing is aligning millions of short reads to a reference genome. Alignment is especially complicated for RNA sequencing (RNA-Seq) because of RNA splicing. A number of RNA-Seq algorithms are available, and claim to align reads with high accuracy and efficiency while detecting splice junctions. RNA-Seq data are discrete in nature; therefore, with reasonable gene models and comparative metrics RNA-Seq data can be simulated to sufficient accuracy to enable meaningful benchmarking of alignment algorithms. The exercise to rigorously compare all viable published RNA-Seq algorithms has not been performed previously.RESULTS: We developed an RNA-Seq simulator that models the main impediments to RNA alignment, including alternative splicing, insertions, deletions, substitutions, sequencing errors and intron signal. We used this simulator to measure the accuracy and robustness of available algorithms at the base and junction levels. Additionally, we used reverse transcription-polymerase chain reaction (RT-PCR) and Sanger sequencing to validate the ability of the algorithms to detect novel transcript features such as novel exons and alternative splicing in RNA-Seq data from mouse retina. A pipeline based on BLAT was developed to explore the performance of established tools for this problem, and to compare it to the recently developed methods. This pipeline, the RNA-Seq Unified Mapper (RUM), performs comparably to the best current aligners and provides an advantageous combination of accuracy, speed and usability.AVAILABILITY: The RUM pipeline is distributed via the Amazon Cloud and for computing clusters using the Sun Grid Engine (http://cbil.upenn.edu/RUM).CONTACT: ggrant@pcbi.upenn.edu; epierce@mail.med.upenn.eduSUPPLEMENTARY INFORMATION: The RNA-Seq sequence reads described in the article are deposited at GEO, accession GSE26248.

A polynomial delay algorithm for the enumeration of bubbles with length constraints in directed graphs and its application to the detection of alternative splicing in RNA-seq data

A Novel Combinatorial Method for Estimating Transcript Expression with RNA-Seq: Bounding the Number of Paths

A novel splicing graph allows a direct comparison between exon-based and splice junction-based approaches to alternative splicing detection

A comprehensive benchmarking of differential splicing tools for RNA-seq analysis at the event level

Linear-time superbubble identification algorithm for genome assembly

A survey of computational methods in transcriptome-wide alternative splicing analysis

Distribution of graph-distances in Boltzmann ensembles of RNA secondary structures

Observations on Novel Splice Junctions from RNA Sequencing Data

Improved methods for RNAseq-based alternative splicing analysis

Machine learning-optimized targeted detection of alternative splicing

Comparative analysis of RNA-Seq alignment algorithms and the RNA-Seq unified mapper (RUM)

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

: intuitive analysis and visualisation of differential alternative splicing using beta distributions

Safe Paths and Sequences for Scalable ILPs in RNA Transcript Assembly Problems

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

Assembly Arena: Benchmarking RNA isoform reconstruction algorithms for nanopore sequencing

An optimal method to segment piecewise poisson distributed signals with application to sequencing data.

Approximation Algorithms for Predicting RNA Secondary Structures with Arbitrary Pseudoknots

Annotation-free quantification of RNA splicing using LeafCutter

Using synthetic RNA to benchmark poly(A) length inference from direct RNA sequencing.

Nucleotide-level distance metrics to quantify alternative splicing implemented in TranD