Background: Haplotype assembly is the process of assigning the different alleles of the variants covered by mapped sequencing reads to the two haplotypes of the genome of a human individual. Long reads, which are nowadays cheaper to produce and more widely available than ever before, have been used to reduce the fragmentation of the assembled haplotypes since their ability to span several variants along the genome. These long reads are also characterized by a high error rate, an issue which may be mitigated, however, with larger sets of reads, when this error rate is uniform across genome positions. Unfortunately, current state-of-the-art dynamic programming approaches designed for long reads deal only with limited coverages. Results: Here, we propose a new method for assembling haplotypes which combines and extends the features of previous approaches to deal with long reads and higher coverages. In particular, our algorithm is able to dynamically adapt the estimated number of errors at each variant site, while minimizing the total number of error corrections necessary for finding a feasible solution. This allows our method to significantly reduce the required computational resources, allowing to consider datasets composed of higher coverages. The algorithm has been implemented in a freely available tool, HapCHAT: Haplotype Assembly Coverage Handling by Adapting Thresholds. An experimental analysis on sequencing reads with up to 60 × coverage reveals improvements in accuracy and recall achieved by considering a higher coverage with lower runtimes. Conclusions: Our method leverages the long-range information of sequencing reads that allows to obtain assembled haplotypes fragmented in a lower number of unphased haplotype blocks. At the same time, our method is also able to deal with higher coverages to better correct the errors in the original reads and to obtain more accurate haplotypes as a result. Availability: HapCHAT is available at http://hapchat.algolab.eu under the GNU Public License (GPL).

Beretta, S., Patterson, M., Zaccaria, S., Della Vedova, G., Bonizzoni, P. (2018). HapCHAT: Adaptive haplotype assembly for efficiently leveraging high coverage in long reads. BMC BIOINFORMATICS, 19(1) [10.1186/s12859-018-2253-8].

HapCHAT: Adaptive haplotype assembly for efficiently leveraging high coverage in long reads

Beretta, Stefano;Patterson, Murray D
;
Zaccaria, Simone;Della Vedova, Gianluca;Bonizzoni, Paola
2018

Abstract

Background: Haplotype assembly is the process of assigning the different alleles of the variants covered by mapped sequencing reads to the two haplotypes of the genome of a human individual. Long reads, which are nowadays cheaper to produce and more widely available than ever before, have been used to reduce the fragmentation of the assembled haplotypes since their ability to span several variants along the genome. These long reads are also characterized by a high error rate, an issue which may be mitigated, however, with larger sets of reads, when this error rate is uniform across genome positions. Unfortunately, current state-of-the-art dynamic programming approaches designed for long reads deal only with limited coverages. Results: Here, we propose a new method for assembling haplotypes which combines and extends the features of previous approaches to deal with long reads and higher coverages. In particular, our algorithm is able to dynamically adapt the estimated number of errors at each variant site, while minimizing the total number of error corrections necessary for finding a feasible solution. This allows our method to significantly reduce the required computational resources, allowing to consider datasets composed of higher coverages. The algorithm has been implemented in a freely available tool, HapCHAT: Haplotype Assembly Coverage Handling by Adapting Thresholds. An experimental analysis on sequencing reads with up to 60 × coverage reveals improvements in accuracy and recall achieved by considering a higher coverage with lower runtimes. Conclusions: Our method leverages the long-range information of sequencing reads that allows to obtain assembled haplotypes fragmented in a lower number of unphased haplotype blocks. At the same time, our method is also able to deal with higher coverages to better correct the errors in the original reads and to obtain more accurate haplotypes as a result. Availability: HapCHAT is available at http://hapchat.algolab.eu under the GNU Public License (GPL).
Articolo in rivista - Articolo scientifico
Haplotype assembly; High coverage; Long reads; Minimum error correction; Single individual haplotyping;
Haplotype assembly; High coverage; Long reads; Minimum error correction; Single individual haplotyping
English
2018
19
1
252
partially_open
Beretta, S., Patterson, M., Zaccaria, S., Della Vedova, G., Bonizzoni, P. (2018). HapCHAT: Adaptive haplotype assembly for efficiently leveraging high coverage in long reads. BMC BIOINFORMATICS, 19(1) [10.1186/s12859-018-2253-8].
File in questo prodotto:
File Dimensione Formato  
2018-BMCBioinfo-HapCHAT.pdf

accesso aperto

Tipologia di allegato: Publisher’s Version (Version of Record, VoR)
Dimensione 989.92 kB
Formato Adobe PDF
989.92 kB Adobe PDF Visualizza/Apri
hapchat.pdf

Solo gestori archivio

Tipologia di allegato: Publisher’s Version (Version of Record, VoR)
Dimensione 989.92 kB
Formato Adobe PDF
989.92 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/10281/204865
Citazioni
  • Scopus 7
  • ???jsp.display-item.citation.isi??? 6
Social impact