COVID-Align: Accurate online alignment of hCoV-19 genomes using a profile HMM - Archive ouverte HAL Access content directly
Journal Articles Bioinformatics Year : 2021

COVID-Align: Accurate online alignment of hCoV-19 genomes using a profile HMM

(1, 2) , (1, 3) , (1, 4) , (1)
1
2
3
4

Abstract

Motivation The first cases of the COVID-19 pandemic emerged in December 2019. Until the end of February 2020, the number of available genomes was below 1,000, and their multiple alignment was easily achieved using standard approaches. Subsequently, the availability of genomes has grown dramatically. Moreover, some genomes are of low quality with sequencing/assembly errors, making accurate re-alignment of all genomes nearly impossible on a daily basis. A more efficient, yet accurate approach was clearly required to pursue all subsequent bioinformatics analyses of this crucial data. Results hCoV-19 genomes are highly conserved, with very few indels and no recombination. This makes the profile HMM approach particularly well suited to align new genomes, add them to an existing alignment and filter problematic ones. Using a core of ∼2,500 high quality genomes, we estimated a profile using HMMER, and implemented this profile in COVID-Align, a user-friendly interface to be used online or as standalone via Docker. The alignment of 1,000 genomes requires ∼50mn on our cluster. Moreover, COVID-Align provides summary statistics, which can be used to determine the sequencing quality and evolutionary novelty of input genomes (e.g. number of new mutations and indels).
Fichier principal
Vignette du fichier
btaa871.pdf (511.09 Ko) Télécharger le fichier
Origin : Publication funded by an institution
Loading...

Dates and versions

pasteur-02995703 , version 1 (09-11-2020)

Licence

Attribution - NonCommercial - CC BY 4.0

Identifiers

Cite

Frédéric Lemoine, Luc Blassel, Jakub Voznica, Olivier Gascuel. COVID-Align: Accurate online alignment of hCoV-19 genomes using a profile HMM. Bioinformatics, 2021, pp.btaa871. ⟨10.1093/bioinformatics/btaa871⟩. ⟨pasteur-02995703⟩
61 View
76 Download

Altmetric

Share

Gmail Facebook Twitter LinkedIn More