Changes

Newer edit →

Sequence assembly

34,164 bytes added, 23:58, 18 December 2010

Created page with "In bioinformatics, sequence assembly refers to aligning and merging fragments of a much longer DNA sequence in order to reconstruct the ori..."

In bioinformatics, sequence assembly refers to aligning and merging fragments of a much longer DNA sequence in order to reconstruct the original sequence. This is needed as DNA sequencing technology cannot read whole genomes in one go, but rather small pieces between 20 and 1000 bases, depending on the technology used. Typically the short fragments, called reads, result from shotgun sequencing genomic DNA, or gene transcript (ESTs).
The problem of sequence assembly can be compared to taking many copies of a book, passing them all through a shredder, and piecing a copy of the book back together from only shredded pieces. The book may have many repeated paragraphs, and some shreds may be modified to have typos. Excerpts from another book may be added in, and some shreds may be completely unrecognizable.

<h2>Genome assemblers</h2>
The first sequence assemblers began to appear in the late 1980s and early 1990s as variants of simpler sequence alignment programs to piece together vast quantities of fragments generated by automated sequencing instruments called DNA sequencers. As the sequenced organisms grew in size and complexity (from small viruses over plasmids to bacteria and finally eukaryotes), the assembly programs needed to increasingly employ more and more sophisticated strategies to handle:
<ul>
<li>terabytes of sequencing data which need processing on computing clusters;</li>
<li>identical and nearly identical sequences (known as repeats) which can, in the worst case, increase the time and space complexity of algorithms exponentially;</li>
<li>and errors in the fragments from the sequencing instruments, which can confound assembly.</li>
</ul>
Faced with the challenge of assembling the first larger eukaryotic genomes, the fruit fly Drosophila melanogaster, in 2000 and the human genome just a year later, scientists developed assemblers like Celera Assembler[1] and Arachne[2] able to handle genomes of 100-300 million base pairs. Subsequent to these efforts, several other groups, mostly at the major genome sequencing centers, built large-scale assemblers, and an open source effort known as AMOS[3] was launched to bring together all the innovations in genome assembly technology under the open source framework.
<h2>EST assemblers</h2>
EST assembly differs from genome assembly in several ways. The sequences for EST assembly are the transcribed mRNA of a cell and represent only a subset of the whole genome. At a first glance, underlying algorithmical problems differ between genome and EST assembly. For instance, genomes often have large amounts of repetitive sequences, mainly in the inter-genic parts. Since ESTs represent gene transcripts, they will not contain these repeats. On the other hand, cells tend to have a certain number of genes that are constantly expressed in very high amounts (housekeeping genes), which again leads to the problem of similar sequences present in high amounts in the data set to be assembled.
Furthermore, genes sometimes overlap in the genome (sense-antisense transcription), and should ideally still be assembled separately. EST assembly is also complicated by features like (cis-) alternative splicing, trans-splicing, single-nucleotide polymorphism, recoding, and post-transcriptional modification.
<h2>De-novo vs. mapping assembly</h2>
In sequence assembly, two different types can be distinguished:
<ol>
<li>de-novo: assembling reads together so that they form a new, previously unknown sequence</li>
<li>mapping: assembling reads against an existing backbone sequence, building a sequence that is similar but not necessarily identical to the backbone sequence</li>
</ol>
In terms of complexity and time requirements, de-novo assemblies are orders of magnitude slower and more memory intensive than mapping assemblies. This is mostly due to the fact that the assembly algorithm need to compare every read with every other read (an operation that is has a complexity of O(<var>n</var>2) but can be reduced to O(<var>n</var> log(<var>n</var>)). Referring to the comparison drawn to shredded books in the introduction: while for mapping assemblies one would have a very similar book as template (perhaps with the names of the main characters and a few locations changed), the de-novo assemblies are more hardcore in a sense as one would not know beforehand whether this would become a science book, or a novel, or a catalogue etc.
<h2>Influence of technological changes</h2>
The complexity of sequence assembly is driven by two major factors: the number of fragments and their lengths. While more and longer fragments allow better identification of sequence overlaps, they also pose problems as the underlying algorithms show quadratic or even exponential complexity behaviour to both number of fragments and their length. And while shorter sequences are faster to align, they also complicate the layout phase of an assembly as shorter reads are more difficult to use with repeats or near identical repeats.
In the earliest days of DNA sequencing, scientists could only gain a few sequences of short length (some dozen bases) after weeks of work in laboratories. Hence, these sequences could be aligned in a few minutes by hand.
In 1975, the Dideoxy termination method (also known as Sanger sequencing) was invented and until shortly after 2000, the technology was improved up to a point were fully automated machines could churn out sequences in a highly parallelised mode 24 hours a day. Large genome centers around the world housed complete farms of these sequencing machines, which in turn led to the necessity of assemblers to be optimised for sequences from whole-genome shotgun sequencing projects where the reads
<ul>
<li>are about 800–900 bases long</li>
<li>contain sequencing artifacts like sequencing and cloning vectors</li>
<li>have error rates between 0.5 and 10%</li>
</ul>
With the Sanger technology, bacterial projects with 20,000 to 200,000 reads could easily be assembled on one computer. Larger ones like the human genome with approximately 35 million reads needed already large computing farms and distributed computing.
By 2004 / 2005, pyrosequencing had been brought to commercial viability by 454 Life Sciences. This new sequencing methods generated reads much shorter than from Sanger sequencing: initially about 100 bases, now 400 bases and expected to grow to 1000 bases by the end of 2010. However, due to the much higher throughput and lower cost than Sanger sequencing, the adoption of this technology by genome centers pushed development of sequence assemblers to deal with this new type of sequences. The sheer amount of data coupled with technology specific error patterns in the reads delayed development of assemblers, at the beginning in 2004 only the Newbler assembler from 454 was available. Presented in mid 2007[4], the hybrid version of the MIRA assembler by Chevreux et al. was the first freely available assembler who could assemble 454 reads and mixtures of 454 reads and Sanger reads; using sequences from different sequencing technologies was subsequently coined hybrid assembly.
Ironically, technological development of sequencing continued to improve in the wrong way (from a sequence assembly point of view). Since 2006, the Solexa technology is available and heavily used to generate roundabout 100 million reads per day on a single sequencing machine. Compare this to the 35 million reads of the human genome project which needed several years to be produced on hundreds of sequencing machines. The downside is that these reads have a length of only 36 bases (expected to grow to 50 bases by the end of 2008). This makes sequence alignment an even more daunting task. Presented by the end of 2007, the SHARCGS assembler[5] by Dohm et al. was the first published assembler that was used for an assembly with Solexa reads, quickly followed by a number of others.
<h2>Greedy algorithm</h2>
Given a set of sequence fragments the object is to find the Shortest common supersequence.
<ol>
<li>calculate pairwise alignments of all fragments</li>
<li>choose two fragments with the largest overlap</li>
<li>merge chosen fragments</li>
<li>repeat step 2. and 3. until only one fragment is left</li>
</ol>
The result is a suboptimal solution to the problem.
<h2>Available assemblers</h2>
The following table lists assemblers that have a de-novo assembly capability on at least one of the supported technologies.[6]

<table class="wikitable" border="1">
<tbody>
<tr>
<th>Name</th>
<th>Type</th>
<th>Technologies</th>
<th>Author</th>
<th>Presented / 
Last updated
</th>
<th>Licence*</th>
<th>Homepage</th>
</tr>
<tr>
<td>ABySS</td>
<td>genomes</td>
<td>Solexa, SOLiD</td>
<td>Simpson, J. et al.</td>
<td>2008 / 2010</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>AMOS</td>
<td>genomes</td>
<td>Sanger, 454</td>
<td>Salzberg, S. et al.</td>
<td>2002? / 2008?</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>Celera WGA Assembler / CABOG</td>
<td>(large) genomes</td>
<td>Sanger, 454, Solexa</td>
<td>Myers, G. et al.; Miller G. et al.</td>
<td>2004 / 2010</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>CLC Genomics Workbench</td>
<td>genomes</td>
<td>Sanger, 454, Solexa, SOLiD</td>
<td>CLC bio</td>
<td>2008 / 2010</td>
<td>C</td>
<td>link</td>
</tr>
<tr>
<td>Edena</td>
<td>genomes</td>
<td>Solexa</td>
<td>D. Hernandez, P. François, L. Farinelli, M. Osteras, and J. Schrenzel.</td>
<td>2008</td>
<td>C</td>
<td>link</td>
</tr>
<tr>
<td>Euler</td>
<td>genomes</td>
<td>Sanger, 454 (,Solexa ?)</td>
<td>Pevzner, P. et al.</td>
<td>2001 / 2006?</td>
<td>(C / NC-A?)</td>
<td>link</td>
</tr>
<tr>
<td>Euler-sr</td>
<td>genomes</td>
<td>454, Solexa</td>
<td>Chaisson, MJ. et al.</td>
<td>2008</td>
<td>NC-A</td>
<td>link</td>
</tr>
<tr>
<td>Forge</td>
<td>(large) genomes, EST, metagenomes</td>
<td>454, Solexa , SOLID, Sanger</td>
<td>Platt, DM, Evers, D.</td>
<td>2010</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>Geneious</td>
<td>genomes</td>
<td>Sanger, 454, Solexa</td>
<td>Biomatters Ltd</td>
<td>2009 / 2010</td>
<td>C</td>
<td>link</td>
</tr>
<tr>
<td>IDBA (Iterative De Bruijn graph short read Assembler)</td>
<td>(large) genomes</td>
<td>Sanger</td>
<td>Yu Peng, Henry C. M. Leung, Siu-Ming Yiu, Francis Y. L. Chin</td>
<td>2010</td>
<td>(C / NC-A?)</td>
<td>link</td>
</tr>
<tr>
<td>MIRA (Mimicking Intelligent Read Assembly)</td>
<td>genomes, ESTs</td>
<td>Sanger, 454, Solexa</td>
<td>Chevreux, B.</td>
<td>1998 / 2010</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>NextGENe</td>
<td>(small genomes?)</td>
<td>454, Solexa, SOLiD</td>
<td>Softgenetics</td>
<td>2008</td>
<td>C</td>
<td>link</td>
</tr>
<tr>
<td>Newbler</td>
<td>genomes, ESTs</td>
<td>454, Sanger</td>
<td>454/Roche</td>
<td>2009</td>
<td>C</td>
<td>link</td>
</tr>
<tr>
<td>Phrap</td>
<td>genomes</td>
<td>Sanger, 454</td>
<td>Green, P.</td>
<td>2002 / 2003 / 2008</td>
<td>C / NC-A</td>
<td>link</td>
</tr>
<tr>
<td>TIGR Assembler</td>
<td>genomic</td>
<td>Sanger</td>
<td>-</td>
<td>1995 / 2003</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>Ray[7]</td>
<td>genomes</td>
<td>Illumina, mix of Illumina and 454, paired or not</td>
<td>Sébastien Boisvert, François Laviolette & Jacques Corbeil.</td>
<td>2010</td>
<td>OS [GNU General Public License]</td>
<td>link</td>
</tr>
<tr>
<td>Sequencher</td>
<td>(small) genomes</td>
<td>Sanger</td>
<td>Gene Codes Corporation</td>
<td>1991 / 2009</td>
<td>C</td>
<td>link</td>
</tr>
<tr>
<td>SeqMan NGen</td>
<td>(small) genomes, ESTs</td>
<td>Sanger, 454, Solexa</td>
<td>DNASTAR</td>
<td> ? / 2008</td>
<td>C</td>
<td>link</td>
</tr>
<tr>
<td>SHARCGS</td>
<td>(small) genomes</td>
<td>Solexa</td>
<td>Dohm et al.</td>
<td>2007 / 2007</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>SOPRA</td>
<td>genomes</td>
<td>Solexa, SOLiD, Sanger, 454</td>
<td>Dayarian, A. et al.</td>
<td>2010 / 2010</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>SSAKE</td>
<td>(small) genomes</td>
<td>Solexa (SOLiD? Helicos?)</td>
<td>Warren, R. et al.</td>
<td>2007 / 2007</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>SOAPdenovo</td>
<td>genomes</td>
<td>Solexa</td>
<td>Li, R. et al.</td>
<td>2009 / 2009</td>
<td>Closed</td>
<td>link</td>
</tr>
<tr>
<td>Staden gap4 package</td>
<td>BACs (, small genomes?)</td>
<td>Sanger</td>
<td>Staden et al.</td>
<td>1991 / 2008</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>VCAKE</td>
<td>(small) genomes</td>
<td>Solexa (SOLiD?, Helicos?)</td>
<td>Jeck, W. et al.</td>
<td>2007 / 2007</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>Phusion assembler</td>
<td>(large) genomes</td>
<td>Sanger</td>
<td>Mullikin JC, et.al.</td>
<td>2003</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>Quality Value Guided SRA (QSRA)</td>
<td>genomes</td>
<td>Sanger, Solexa</td>
<td>Bryant DW, et.al.</td>
<td>2009</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td>Velvet (algorithm)</td>
<td>(small) genomes</td>
<td>Sanger, 454, Solexa, SOLiD</td>
<td>Zerbino, D. et al.</td>
<td>2007 / 2009</td>
<td>OS</td>
<td>link</td>
</tr>
<tr>
<td style="border-top: #333 1px solid" colspan="7">*Licences: OS = Open Source; C = Commercial; C / NC-A = Commercial but free for non-commercial and academics; Brackets = unclear, but most likely C / NC-A</td>
</tr>
</tbody>
</table>

<h2>See also</h2>
<ul>
<li>Sequence alignment</li>
<li>Genome assembly</li>
</ul>
<h2>References</h2>
<ol class="references">
<li id="cite_note-0"><a href="#cite_ref-0">^</a> Myers EW, Sutton GG, Delcher AL, et al. (March 2000). <a class="external text" href="http://www.sciencemag.org/cgi/pmidlookup?view=long&pmid=10731133" rel="nofollow">"A whole-genome assembly of Drosophila"</a>. Science 287 (5461): 2196–204. <a class="mw-redirect" title="PubMed Identifier" href="/wiki/PubMed_Identifier">PMID</a> <a class="external text" href="http://www.ncbi.nlm.nih.gov/pubmed/10731133" rel="nofollow">10731133</a>. <a class="external free" href="http://www.sciencemag.org/cgi/pmidlookup?view=long&pmid=10731133" rel="nofollow">http://www.sciencemag.org/cgi/pmidlookup?view=long&pmid=10731133</a>.<span class="Z3988" title="ctx_ver=Z39.88-2004&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&rft.genre=article&rft.atitle=A+whole-genome+assembly+of+Drosophila&rft.jtitle=Science&rft.aulast=Myers+EW%2C+Sutton+GG%2C+Delcher+AL%2C+%27%27et+al.%27%27&rft.au=Myers+EW%2C+Sutton+GG%2C+Delcher+AL%2C+%27%27et+al.%27%27&rft.date=March+2000&rft.volume=287&rft.issue=5461&rft.pages=2196%E2%80%93204&rft_id=info:pmid/10731133&rft_id=http%3A%2F%2Fwww.sciencemag.org%2Fcgi%2Fpmidlookup%3Fview%3Dlong%26pmid%3D10731133&rfr_id=info:sid/en.wikipedia.org:Sequence_assembly"> </li>
<li id="cite_note-1"><a href="#cite_ref-1">^</a> Batzoglou S, Jaffe DB, Stanley K, et al. (January 2002). <a class="external text" href="http://www.genome.org/cgi/pmidlookup?view=long&pmid=11779843" rel="nofollow">"ARACHNE: a whole-genome shotgun assembler"</a>. Genome Res. 12 (1): 177–89. <a title="Digital object identifier" href="/wiki/Digital_object_identifier">doi</a>:<a class="external text" href="http://dx.doi.org/10.1101%2Fgr.208902" rel="nofollow">10.1101/gr.208902</a>. <a class="mw-redirect" title="PubMed Identifier" href="/wiki/PubMed_Identifier">PMID</a> <a class="external text" href="http://www.ncbi.nlm.nih.gov/pubmed/11779843" rel="nofollow">11779843</a>. <a title="PubMed Central" href="/wiki/PubMed_Central">PMC</a> <a class="external text" href="http://www.pubmedcentral.gov/articlerender.fcgi?tool=pmcentrez&artid=155255" rel="nofollow">155255</a>. <a class="external free" href="http://www.genome.org/cgi/pmidlookup?view=long&pmid=11779843" rel="nofollow">http://www.genome.org/cgi/pmidlookup?view=long&pmid=11779843</a>.<span class="Z3988" title="ctx_ver=Z39.88-2004&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&rft.genre=article&rft.atitle=ARACHNE%3A+a+whole-genome+shotgun+assembler&rft.jtitle=Genome+Res.&rft.aulast=Batzoglou+S%2C+Jaffe+DB%2C+Stanley+K%2C+%27%27et+al.%27%27&rft.au=Batzoglou+S%2C+Jaffe+DB%2C+Stanley+K%2C+%27%27et+al.%27%27&rft.date=January+2002&rft.volume=12&rft.issue=1&rft.pages=177%E2%80%9389&rft_id=info:doi/10.1101%2Fgr.208902&rft_id=info:pmid/11779843&rft_id=http%3A%2F%2Fwww.genome.org%2Fcgi%2Fpmidlookup%3Fview%3Dlong%26pmid%3D11779843&rfr_id=info:sid/en.wikipedia.org:Sequence_assembly"> </li>
<li id="cite_note-2"><a href="#cite_ref-2">^</a> <a class="external text" href="http://amos.sourceforge.net/" rel="nofollow">AMOS page</a> with links to various papers</li>
<li id="cite_note-3"><a href="#cite_ref-3">^</a> Copy in Google groups of the <a class="external text" href="http://groups.google.com/group/bionet.software/browse_thread/thread/b34b348011d04f0e?fwc=1" rel="nofollow">post announcing MIRA 2.9.8 hybrid version</a> in the bionet.software Usenet group</li>
<li id="cite_note-4"><a href="#cite_ref-4">^</a> Dohm JC, Lottaz C, Borodina T, Himmelbauer H (November 2007). <a class="external text" href="http://www.genome.org/cgi/pmidlookup?view=long&pmid=17908823" rel="nofollow">"SHARCGS, a fast and highly accurate short-read assembly algorithm for de novo genomic sequencing"</a>. Genome Res. 17 (11): 1697–706. <a title="Digital object identifier" href="/wiki/Digital_object_identifier">doi</a>:<a class="external text" href="http://dx.doi.org/10.1101%2Fgr.6435207" rel="nofollow">10.1101/gr.6435207</a>. <a class="mw-redirect" title="PubMed Identifier" href="/wiki/PubMed_Identifier">PMID</a> <a class="external text" href="http://www.ncbi.nlm.nih.gov/pubmed/17908823" rel="nofollow">17908823</a>. <a title="PubMed Central" href="/wiki/PubMed_Central">PMC</a> <a class="external text" href="http://www.pubmedcentral.gov/articlerender.fcgi?tool=pmcentrez&artid=2045152" rel="nofollow">2045152</a>. <a class="external free" href="http://www.genome.org/cgi/pmidlookup?view=long&pmid=17908823" rel="nofollow">http://www.genome.org/cgi/pmidlookup?view=long&pmid=17908823</a>.<span class="Z3988" title="ctx_ver=Z39.88-2004&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&rft.genre=article&rft.atitle=SHARCGS%2C+a+fast+and+highly+accurate+short-read+assembly+algorithm+for+de+novo+genomic+sequencing&rft.jtitle=Genome+Res.&rft.aulast=Dohm+JC%2C+Lottaz+C%2C+Borodina+T%2C+Himmelbauer+H&rft.au=Dohm+JC%2C+Lottaz+C%2C+Borodina+T%2C+Himmelbauer+H&rft.date=November+2007&rft.volume=17&rft.issue=11&rft.pages=1697%E2%80%93706&rft_id=info:doi/10.1101%2Fgr.6435207&rft_id=info:pmid/17908823&rft_id=http%3A%2F%2Fwww.genome.org%2Fcgi%2Fpmidlookup%3Fview%3Dlong%26pmid%3D17908823&rfr_id=info:sid/en.wikipedia.org:Sequence_assembly"> </li>
<li id="cite_note-5"><a href="#cite_ref-5">^</a> <a class="external text" href="http://seqanswers.com/forums/showthread.php?t=43" rel="nofollow">list of software including mapping assemblers in the SeqAnswers discussion forum.</a></li>
<li id="cite_note-6"><a href="#cite_ref-6">^</a> Boisvert S, Laviolette F, Corbeil J. (October 2010). <a class="external text" href="http://www.liebertonline.com/doi/abs/10.1089/cmb.2009.0238" rel="nofollow">"Ray: simultaneous assembly of reads from a mix of high-throughput sequencing technologies."</a>. J Comput Biol. 17 (11): 1519-33. <a title="Digital object identifier" href="/wiki/Digital_object_identifier">doi</a>:<a class="external text" href="http://dx.doi.org/10.1089%2Fcmb.2009.0238" rel="nofollow">10.1089/cmb.2009.0238</a>. <a class="mw-redirect" title="PubMed Identifier" href="/wiki/PubMed_Identifier">PMID</a> <a class="external text" href="http://www.ncbi.nlm.nih.gov/pubmed/20958248" rel="nofollow">20958248</a>. <a class="external free" href="http://www.liebertonline.com/doi/abs/10.1089/cmb.2009.0238" rel="nofollow">http://www.liebertonline.com/doi/abs/10.1089/cmb.2009.0238</a>.</li>
</ol>

WikiSysop

Bureaucrats, Administrators

1,081

edits

kogic.kr β

Changes

Sequence assembly

kogic.kr ^β