Improved sensitivity and reliability of anchor based genome alignment

Uricaru , Raluca
        Université de Montpellier 2Centre National de la Recherche Scientifique, MontpellierMontpellierMontpellier.; Michotey , Célia
        INRA
        , Jouy-En-Josas .
         UR 1077 Mathématique, Informatique et Génome; Noé , Laurent
        .; Chiapello , Helene
        INRA
         .
         UR 1077 Mathématique, Informatique et Génome; Rivals , Eric
        INRA
         .

Improved sensitivity and reliability of anchor based genome alignment

2009

Uricaru , Raluca (Université Montpellier 2 (Sciences et Techniques)Université de Montpellier 2Centre National de la Recherche Scientifique, MontpellierMontpellierMontpellier(France).) | Michotey , Célia (INRA , Jouy-En-Josas (France). UR 1077 Mathématique, Informatique et Génome) | Noé , Laurent(collaborateur) (Université de Lille, Villeuneuve d'Ascq(France).) | Chiapello , Helene(auteur de correspondance) (INRA (France). UR 1077 Mathématique, Informatique et Génome) | Rivals , Eric(collaborateur) (INRA (France). )

Whole genome alignment is a challenging problem in computational comparative genomics. It is essential for the functional annotation of genomes, the understanding of their evolution, and for phylogenomics. Many global alignment programs are heuristic variations on the anchor based strategy, which relies on the initial detection of similarities and their selection in an ordered chain. Considering that alignment tools fail to align some pairs of bacterial strains, we investigate whether this is intrinsically due to the strategy or to a lack of sensitivity of the similarity detection method. For this, we implement and compare 6 programs based on three different detection methods (from exact matches to local alignments) on a large benchmark set. Our results suggest that the sensitivity of well known methods, like MGA or Mauve, can be greatly improved in the case of divergent genomes if one exploits spaced seeds at the detection phase. In other cases, such methods yield alignments that cover nearly the whole genome. Then, we focus on global reliability of alignments: should an aligned pair of segments be included in the global genome alignment? We investigate this reliability according to both the segment ”alignability” and to inclusion of orthologs. Again, we provide evidence that for both close and divergent genomes, one of our programs, YH, achieves alignments with sometimes a lower coverage, but a higher inclusion of orthologs. It opens the way to the first reliable alignments for some highly divergent species like Buchnera aphidicola or Prochlorococcus marinus.

Показать больше [+]