Sailfish: Alignment-free Isoform Quantification from RNA-seq Reads using Lightweight Algorithms

Sailfish: Alignment-free Isoform Quantification from RNA-seq Reads using Lightweight Algorithms
Rob Patro (1), Stephen M. Mount (2), Carl Kingsford (1) ((1) Lane Center for Computational Biology, School of Computer Science, Carnegie Mellon University, (2) Department of Cell Biology and Molecular Genetics and Center for Bioinformatics and Computational Biology, University of Maryland)
(Submitted on 16 Aug 2013)

RNA-seq has rapidly become the de facto technique to measure gene expression. However, the time required for analysis has not kept up with the pace of data generation. Here we introduce Sailfish, a novel computational method for quantifying the abundance of previously annotated RNA isoforms from RNA-seq data. Sailfish entirely avoids mapping reads, which is a time-consuming step in all current methods. Sailfish provides quantification estimates much faster than existing approaches (typically 20-times faster) without loss of accuracy.

8 thoughts on “Sailfish: Alignment-free Isoform Quantification from RNA-seq Reads using Lightweight Algorithms

  1. Naively, given that Sailfish does not explicitly use the added specificity from long reads or mates in paired-end data, it would be interesting to see its performance on the synthetic data as a function of read length +/- pairing.

  2. We suspect that the accuracy of Sailfish is due to a tradeoff between not using information about the co-occurrence of k-mers within a single read (or paired end information) and bypassing problems with read mapping (e.g. what to do about reads that map multiple places. and how much and what sort of errors to tolerate). As people start to use the tool, we are interested in hearing about performance on different data types. I’m thinking that important variables will be the quality of the transcriptome annotation and the extent (and range) of sequence divergence between the reference transcriptome and the sample, but there may be others. Also, we don’t see much sensitivity to the value of k, but that could also be different with real data. We have set up a user’s group – sailfish-users on googlegroups.com (short URL http://ongen.us/SForum) – where people can share results.

  3. Pingback: Most viewed on Haldane’s Sieve: August 2013 | Haldane's Sieve

  4. Pingback: Some preprint comment streams at Haldane’s sieve and related sites | Haldane's Sieve

  5. Pingback: Sifting through 2013 with Haldane’s Sieve | Haldane's Sieve

  6. Pingback: Sailfish: Alignment-free isoform quantification from RNA-seq reads

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google photo

You are commenting using your Google account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s