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.

RabbitTrim: Highly Optimized Trimming of Illumina Sequencing Data on Multi-core Platforms

RabbitFX: Efficient Framework for FASTA/Q File Parsing on Modern Multi-Core Platforms.

RabbitQCPlus 2.0: More efficient and versatile quality control for sequencing data

RabbitQCPlus: More Efficient Quality Control for Sequencing Data.

Trimmomatic: a flexible trimmer for Illumina sequence data

RabbitQC: high-speed scalable quality control for sequencing data

RabbitV: Fast Detection of Viruses and Microorganisms in Sequencing Data on Multi-Core Architectures

RabbitSAlign: Accelerating Short-Read Alignment for CPU-GPU Heterogeneous Platforms

StreamingTrim 1.0: a Java software for dynamic trimming of 16S rRNA sequence data from metagenetic studies

FIGARO: An efficient and objective tool for optimizing microbiome rRNA gene trimming parameters

RabbitKSSD: Accelerating Genome Distance Estimation on Modern Multi-Core Architectures.

RabbitMash: Accelerating Hash-Based Genome Analysis on Modern Multi-Core Architectures.

Minirmd: Accurate and Fast Duplicate Removal Tool for Short Reads Via Multiple Minimizers.

RabbitVar: Ultra-Fast and Accurate Somatic Small-Variant Calling on Multi-Core Architectures

AdapterRemoval v2: rapid adapter trimming, identification, and read merging

An efficient trimming algorithm based on multi-feature fusion scoring model for NGS data.

Sequence-matching adapter trimmers generate consistent quality and assembly metrics for Illumina sequencing of RNA viruses

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

InfoTrim: A DNA Read Quality Trimmer Using Entropy.

Sam2bam: High-Performance Framework for NGS Data Preprocessing Tools

Accelerating spliced alignment of long RNA sequencing reads using parallel maximal exact match retrieval