Computational method dramatically speeds up estimates of gene expression

Apr 20, 2014

With gene expression analysis growing in importance for both basic researchers and medical practitioners, researchers at Carnegie Mellon University and the University of Maryland have developed a new computational method that dramatically speeds up estimates of gene activity from RNA sequencing (RNA-seq) data.

With the new method, dubbed Sailfish after the famously speedy fish, estimates of that previously took many hours can be completed in a few minutes, with accuracy that equals or exceeds previous methods. The researchers' report on their new method is being published online April 20 by the journal Nature Biotechnology.

Gigantic repositories of RNA-seq data now exist, making it possible to re-analyze experiments in light of new discoveries. "But 15 hours a pop really starts to add up, particularly if you want to look at 100 experiments," said Carl Kingsford, an associate professor in CMU's Lane Center for Computational Biology. "With Sailfish, we can give researchers everything they got from previous methods, but faster."

Though an organism's genetic makeup is static, the activity of individual genes varies greatly over time, making gene expression an important factor in understanding how organisms work and what occurs during disease processes. Gene activity can't be measured directly, but can be inferred by monitoring RNA, the molecules that carry information from the genes for producing proteins and other cellular activities. RNA-seq is a leading method for producing these snapshots of gene expression; in genomic medicine, it has proven particularly useful in analyzing certain cancers.

The RNA-seq process results in short sequences of RNA, called "reads." In previous methods, the RNA molecules from which they originated could be identified and measured only by painstakingly mapping these reads to their original positions in the larger molecules.

But Kingsford, working with Rob Patro, a post-doctoral researcher in the Lane Center, and Stephen M. Mount, an associate professor in Maryland's Department of Cell Biology and Molecular Genetics and its Center for Bioinformatics and Computational Biology, found that the time-consuming mapping step could be eliminated. Instead, they found they could allocate parts of the reads to different types of RNA molecules, much as if each read acted as several votes for one molecule or another.

Without the mapping step, Sailfish can complete its RNA analysis 20-30 times faster than previous methods.

This numerical approach might not be as intuitive as a map to a biologist, but it makes perfect sense to a computer scientist, Kingsford said. Moreover, the Sailfish method is more robust—better able to tolerate errors in the reads or differences between individuals' genomes. These errors can prevent some reads from being mapped, he explained, but the Sailfish method can make use of all the RNA read "votes," which improves the method's accuracy.

Explore further: Rapid and accurate mRNA detection in plant tissues

More information: The Sailfish code has been released and is available for download at

Paper: Sailfish enables alignment-free isoform quantification from RNA-seq reads using lightweight algorithms, Nature Biotechnology, DOI: 10.1038/nbt.2862

add to favorites email to friend print save as pdf

Related Stories

Rapid and accurate mRNA detection in plant tissues

Apr 17, 2014

Gene expression is the process whereby the genetic information of DNA is used to manufacture functional products, such as proteins, which have numerous different functions in living organisms. Messenger RNA (mRNA) serves ...

New bioinformatics tool to visualize transcriptomes

Mar 09, 2014

ZENBU, a new, freely available bioinformatics tool developed at the RIKEN Center for Life Science Technology in Japan, enables researchers to quickly and easily integrate, visualize and compare large amounts of genomic information ...

Recommended for you

For legume plants, a new route from shoot to root

19 hours ago

A new study shows that legume plants regulate their symbiotic relationship with soil bacteria by using cytokinins—signaling molecules— that are transmitted through the plant structure from leaves into ...

Controlling the transition between generations

Sep 18, 2014

Rafal Ciosk and his group at the FMI have identified an important regulator of the transition from germ cell to embryonic cell. LIN-41 prevents the premature onset of embryonic transcription in oocytes poised ...

User comments : 0