Bowtie (sequence analysis) explained

Bowtie
Author:Ben Langmead
Cole Trapnell
Mihai Pop
Steven Salzberg
Developer:Ben Langmead et al.,
Repo:Bowtie:
Bowtie 2:
Programming Language:C++
Operating System:Linux
macOS
Windows
Genre:Bioinformatics

Bowtie is a software package commonly used for sequence alignment and sequence analysis in bioinformatics.[1] The source code for the package is distributed freely and compiled binaries are available for Linux, macOS and Windows platforms. As of 2017, the Genome Biology paper describing the original Bowtie method has been cited more than 11,000 times. Bowtie is open-source software and is currently maintained by Johns Hopkins University.

History

The Bowtie sequence aligner was originally developed by Ben Langmead et al. at the University of Maryland in 2009. The aligner is typically used with short reads and a large reference genome, or for whole genome analysis. Bowtie is promoted as "an ultrafast, memory-efficient short aligner for short DNA sequences." The speed increase of Bowtie is partly due to implementing the Burrows–Wheeler transform for aligning,[2] which reduces the memory footprint (typically to around 2.2GB for the human genome);[3] a similar method is used by the BWA[4] and SOAP2[5] alignment methods.

Bowtie conducts a quality-aware, greedy, randomized, depth-first search through the space of possible alignments. Because the search is greedy, the first valid alignment encountered by Bowtie will not necessarily be the 'best' in terms of the number of mismatches or in terms of quality.

Bowtie is used as a sequence aligner by a number of other related bioinformatics algorithms, including TopHat,[6] Cufflinks[7] and the CummeRbund Bioconductor package.[8]

Bowtie 2

On 16 October 2011, the developers released a beta fork of the project called Bowtie 2.[9] In addition to the Burrows-Wheeler transform, Bowtie 2 also uses an FM-index (similar to a suffix array) to keep its memory footprint small. Due to its implementation, Bowtie 2 is more suited to finding longer, gapped alignments in comparison with the original Bowtie method. There is no upper limit on read length in Bowtie 2 and it allows alignments to overlap ambiguous characters in the reference.

External links

Notes and References

  1. Langmead. Ben. Cole Trapnell. Mihai Pop. Steven L Salzberg. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biology. 4 March 2009. 10. 3. 10:R25. 10.1186/gb-2009-10-3-r25. 19261174. 2690996. free .
  2. Ferragina. Paolo. Manzini. Giovanni. Indexing compressed text. Journal of the ACM. 2005 . 52. 4. 552–581. 10.1145/1082036.1082039. 6200428 .
  3. Web site: Bowtie: An ultrafast, memory-efficient short read aligner - SourceForge. 29 November 2013.
  4. Li. H.. Durbin. R.. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 18 May 2009. 25. 14. 1754–1760. 10.1093/bioinformatics/btp324. 19451168. 2705234.
  5. Li. R.. Yu. C.. Li. Y.. Lam. T.-W.. Yiu. S.-M.. Kristiansen. K.. Wang. J.. SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 3 June 2009. 25. 15. 1966–1967. 10.1093/bioinformatics/btp336. 19497933. free.
  6. Trapnell. C.. Pachter. L.. Lior Pachter . Salzberg. S. L.. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 16 March 2009. 25. 9. 1105–1111. 10.1093/bioinformatics/btp120. 19289445. 2672628.
  7. Trapnell. Cole. Roberts. Adam. Goff. Loyal. Pertea. Geo. Kim. Daehwan. Kelley. David R. Pimentel. Harold. Salzberg. Steven L. Rinn. John L. Pachter. Lior. Lior Pachter . Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nature Protocols. 1 March 2012. 7. 3. 562–578. 10.1038/nprot.2012.016. 22383036. 3334321.
  8. Web site: CummeRbund - An R package for persistent storage, analysis, and visualization of RNA-Seq from cufflinks output. 11 August 2015.
  9. Langmead. Ben. Salzberg. Steven L. Fast gapped-read alignment with Bowtie 2. Nature Methods. 4 March 2012. 9. 4. 357–359. 10.1038/nmeth.1923. 22388286. 3322381.