Aligning the unalignable: bacteriophage whole genome alignments

BMC Bioinformatics

Methodology Article

Sèverine Bérard¹  Nicolas Pompidor²  Annie Chateau³  Krister M. Swenson³  Anne Bergeron⁴  Paul Guertin⁵ 

[1] ISEM, CNRS - Univ. Montpellier, Montpellier, France;LIRMM, CNRS - Univ. Montpellier, 161 rue Ada, 34392, Montpellier, France;LIRMM, CNRS - Univ. Montpellier, 161 rue Ada, 34392, Montpellier, France;LIRMM, CNRS - Univ. Montpellier, 161 rue Ada, 34392, Montpellier, France;IBC Institut de Biologie Computationnelle, Montpellier, France;LaCIM, Université du Québec à Montréal, Montréal, Canada;LaCIM, Université du Québec à Montréal, Montréal, Canada;Département de mathématiques, Collège André-Grasset, Montréal, Canada;

关键词: Bacteriophages;    Whole genome alignments;    Comparative genomics;    Partial orders;   

DOI  :  10.1186/s12859-015-0869-5

 received in 2015-10-06, accepted in 2015-12-22,  发布年份 2016

来源: Springer

【 摘 要 】

BackgroundIn recent years, many studies focused on the description and comparison of large sets of related bacteriophage genomes. Due to the peculiar mosaic structure of these genomes, few informative approaches for comparing whole genomes exist: dot plots diagrams give a mostly qualitative assessment of the similarity/dissimilarity between two or more genomes, and clustering techniques are used to classify genomes. Multiple alignments are conspicuously absent from this scene. Indeed, whole genome aligners interpret lack of similarity between sequences as an indication of rearrangements, insertions, or losses. This behavior makes them ill-prepared to align bacteriophage genomes, where even closely related strains can accomplish the same biological function with highly dissimilar sequences.ResultsIn this paper, we propose a multiple alignment strategy that exploits functional collinearity shared by related strains of bacteriophages, and uses partial orders to capture mosaicism of sets of genomes. As classical alignments do, the computed alignments can be used to predict that genes have the same biological function, even in the absence of detectable similarity. The Alpha aligner implements these ideas in visual interactive displays, and is used to compute several examples of alignments of Staphylococcus aureus and Mycobacterium bacteriophages, involving up to 29 genomes. Using these datasets, we prove that Alpha alignments are at least as good as those computed by standard aligners. Comparison with the progressiveMauve aligner – which implements a partial order strategy, but whose alignments are linearized – shows a greatly improved interactive graphic display, while avoiding misalignments.ConclusionsMultiple alignments of whole bacteriophage genomes work, and will become an important conceptual and visual tool in comparative genomics of sets of related strains.A python implementation of Alpha, along with installation instructions for Ubuntu and OSX, is available on bitbucket (https://bitbucket.org/thekswenson/alpha).

【 授权许可】

CC BY   
© Bérard et al. 2015

【 图 表 】

Fig. 4

Fig. 5

Fig. 6

Fig. 2

Fig. 5

Fig. 1

12951_2017_255_Article_IEq39.gif

Fig. 1

Fig. 2

12951_2015_155_Article_IEq62.gif

12951_2017_292_Article_IEq1.gif

Fig. 3

Fig. 5