EMSAR: Estimation of transcript abundance from RNA-seq data by mappability-based segmentation and reclustering

Soohyun Lee, Chae Hwa Seo, Burak Han Alver, Sanghyuk Lee, Peter J. Park

Research output: Contribution to journalArticlepeer-review

11 Scopus citations

Abstract

Background: RNA-seq has been widely used for genome-wide expression profiling. RNA-seq data typically consists of tens of millions of short sequenced reads from different transcripts. However, due to sequence similarity among genes and among isoforms, the source of a given read is often ambiguous. Existing approaches for estimating expression levels from RNA-seq reads tend to compromise between accuracy and computational cost. Results: We introduce a new approach for quantifying transcript abundance from RNA-seq data. EMSAR (Estimation by Mappability-based Segmentation And Reclustering) groups reads according to the set of transcripts to which they are mapped and finds maximum likelihood estimates using a joint Poisson model for each optimal set of segments of transcripts. The method uses nearly all mapped reads, including those mapped to multiple genes. With an efficient transcriptome indexing based on modified suffix arrays, EMSAR minimizes the use of CPU time and memory while achieving accuracy comparable to the best existing methods. Conclusions: EMSAR is a method for quantifying transcripts from RNA-seq data with high accuracy and low computational cost. EMSAR is available at https://github.com/parklab/emsar

Original languageEnglish
Article number278
JournalBMC Bioinformatics
Volume16
Issue number1
DOIs
StatePublished - 3 Sep 2015

Bibliographical note

Publisher Copyright:
© 2015 Lee et al.

Keywords

  • Expression quantification
  • Isoforms
  • Multi-reads
  • Optimization
  • Suffix array

Fingerprint

Dive into the research topics of 'EMSAR: Estimation of transcript abundance from RNA-seq data by mappability-based segmentation and reclustering'. Together they form a unique fingerprint.

Cite this