WhatsHap: Haplotype assembly for future-generation sequencing reads

التفاصيل البيبلوغرافية
العنوان: WhatsHap: Haplotype assembly for future-generation sequencing reads
المؤلفون: M. Patterson, T. Marschall, PISANTI, NADIA, L. Van Iersel, L. Stougie, G. W. Klau, A. Schönhuth
المساهمون: M., Patterson, T., Marschall, Pisanti, Nadia, L., Van Iersel, L., Stougie, G. W., Klau, A., Schönhuth
بيانات النشر: Springer
DEU
Dordrecht
سنة النشر: 2014
المجموعة: ARPI - Archivio della Ricerca dell'Università di Pisa
مصطلحات موضوعية: Haplotype Assembly, Fixed Parameter Tractable, Dynamic Programming
الوصف: The human genome is diploid, that is each of its chromosomes comes in two copies. This requires to phase the single nucleotide polymorphisms (SNPs), that is, to assign them to the two copies, beyond just detecting them. The resulting haplotypes, lists of SNPs belonging to each copy, are crucial for downstream analyses in population genetics. Currently, statistical approaches, which avoid making use of direct read information, constitute the state-of-the-art. Haplotype assembly, which addresses phasing directly from sequencing reads, suffers from the fact that sequencing reads of the current generation are too short to serve the purposes of genome-wide phasing. Future sequencing technologies, however, bear the promise to generate reads of lengths and error rates that allow to bridge all SNP positions in the genome at sufficient amounts of SNPs per read. Existing haplotype assembly approaches, however, profit precisely, in terms of computational complexity, from the limited length of current-generation reads, because their runtime is usually exponential in the number of SNPs per sequencing read. This implies that such approaches will not be able to exploit the benefits of long enough, future-generation reads. Here, we suggest WhatsHap, a novel dynamic programming approach to haplotype assembly. It is the first approach that yields provably optimal solutions to the weighted minimum error correction (wMEC) problem in runtime linear in the number of SNPs per sequencing read, making it suitable for future-generation reads. WhatsHap is a fixed parameter tractable (FPT) approach with coverage as the parameter. We demonstrate that WhatsHap can handle datasets of coverage up to 20x, processing chromosomes on standard workstations in only 1-2 hours. Our simulation study shows that the quality of haplotypes assembled by WhatsHap significantly improves with increasing read length, both in terms of genome coverage as well as in terms of switch errors. The switch error rates we achieve in our simulations are superior to those ...
نوع الوثيقة: conference object
وصف الملف: STAMPA
اللغة: English
العلاقة: info:eu-repo/semantics/altIdentifier/isbn/9783319052687; info:eu-repo/semantics/altIdentifier/wos/WOS:000345119100019; ispartofbook:18th ACM Annual International Conference on Research in Computational Molecular Biology (RECOMB); 18th ACM International Conference on Research in COMputational Biology (RECOMB); volume:8394; firstpage:237; lastpage:249; numberofpages:13; serie:LECTURE NOTES IN COMPUTER SCIENCE; http://hdl.handle.net/11568/460292Test; info:eu-repo/semantics/altIdentifier/scopus/2-s2.0-84898978697
DOI: 10.1007/978-3-319-05269-4-19
الإتاحة: https://doi.org/10.1007/978-3-319-05269-4-19Test
http://hdl.handle.net/11568/460292Test
حقوق: info:eu-repo/semantics/closedAccess
رقم الانضمام: edsbas.F9F4FA8
قاعدة البيانات: BASE