High-Performance Haplotype Assembly

Marco Aldinucci (1) , Andrea Bracciali (2) , Tobias Marschall (3, 4) , Murray Patterson (5) , Nadia Pisanti (6, 7) , Massimo Torquati (6)

1 University of Torino
2 University of Stirling
3 Max Planck Institute for Informatics [Saarbrücken]
4 Saarland University [Saarbrücken]
5 BPGE - Bioinformatique, phylogénie et génomique évolutive
6 DI - Dipartimento di Informatica [Pisa]
7 ERABLE - Equipe de recherche européenne en algorithmique et biologie formelle et expérimentale

The problem of Haplotype Assembly is an essential step in human genome analysis. It is typically formalised as the Minimum Error Correction (MEC) problem which is NP-hard. MEC has been approached using heuristics, integer linear programming, and fixed-parameter tractability (FPT), including approaches whose runtime is exponential in the length of the DNA fragments obtained by the sequencing process. Technological improvements are currently increasing fragment length, which drastically elevates computational costs for such methods. We present pWhatsHap, a multi-core parallelisation of WhatsHap, a recent FPT optimal approach to MEC. WhatsHap moves complexity from fragment length to fragment overlap and is hence of particular interest when considering sequencing technology's current trends. pWhat-sHap further improves the efficiency in solving the MEC problem, as shown by experiments performed on datasets with high coverage.

Domaines: Bio-informatique [q-bio.QM]

Marco Aldinucci, Andrea Bracciali, Tobias Marschall, Murray Patterson, Nadia Pisanti, et al.. High-Performance Haplotype Assembly. Computational Intelligence Methods for Bioinformatics and Biostatistics (CIBB'15), 2015, Cambridge, United Kingdom. pp.281 - 258, ⟨10.1016/j.compbiolchem.2005.05.001⟩. ⟨hal-01526652⟩