Abstract:Motivation: RNA-Seq technique has been demonstrated as a revolutionary means for exploring transcriptome because it provides deep coverage and base pair-level resolution. RNA-Seq quantification is proven to be an efficient alternative to Microarray technique in gene expression study, and it is a critical component in RNA-Seq differential expression analysis. Most existing RNA-Seq quantification tools require the alignments of fragments to either a genome or a transcriptome, entailing a time-consuming and intricate alignment step. To improve the performance of RNA-Seq quantification, an alignment-free method, Sailfish, has been recently proposed to quantify transcript abundances using all k-mers in the transcriptome, demonstrating the feasibility of designing an efficient alignment-free method for transcriptome quantification. Even though Sailfish is substantially faster than alternative alignment-dependent methods such as Cufflinks, using all k-mers in the transcriptome quantification impedes the scalability of the method. Results: We propose a novel RNA-Seq quantification method, RNA-Skim, which partitions the transcriptome into disjoint transcript clusters based on sequence similarity, and introduces the notion of sig-mers, which are a special type of k-mers uniquely associated with each cluster. We demonstrate that the sig-mer counts within a cluster are sufficient for estimating transcript abundances with accuracy comparable with any state-of-the-art method. This enables RNA-Skim to perform transcript quantification on each cluster independently, reducing a complex optimization problem into smaller optimization tasks that can be run in parallel. As a result, RNA-Skim uses <4% of the k-mers and <10% of the CPU time required by Sailfish. It is able to finish transcriptome quantification in <10 min per sample by using just a single thread on a commodity computer, which represents >100 speedup over the state-of-the-art alignment-based methods, while delivering comparable or higher accuracy. Availability and implementation: The software is available at http://www.csbio.unc.edu/rs. Contact: weiwang@cs.ucla.edu Supplementary information: Supplementary data are available at Bioinformatics online.

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

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

A Novel Analytical Strategy To Identify Fusion Transcripts Between Repetitive Elements And Protein Coding-Exons Using Rna-Seq

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

RNA-Skim: a Rapid Method for RNA-Seq Quantification at Transcript Level.

Accurate, Fast and Lightweight Clustering of de novo Transcriptomes using Fragment Equivalence Classes

Evaluating Optimal Safe Flows Decomposition for RNA Assembly

GeneScissors: a comprehensive approach to detecting and correcting spurious transcriptome inference owing to RNA-seq reads misalignment.

On de novo Bridging Paired-end RNA-seq Data

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

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

STAR: ultrafast universal RNA-seq aligner

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

Union Exon Based Approach for Rna-Seq Gene Quantification: to Be or Not to Be?

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

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

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

On the complexity of Minimum Path Cover with Subpath Constraints for multi-assembly

Forseti : a mechanistic and predictive model of the splicing status of scRNA-seq reads

Recursive Indirect-Paths Modularity (RIP-M) for Detecting Community Structure in RNA-Seq Co-expression Networks

Count-based differential expression analysis of RNA sequencing data using R and Bioconductor