RecGraph: recombination-aware alignment of sequences to variation graphs
Jorge Avila Cartes,Paola Bonizzoni,Simone Ciccolella,Gianluca Della Vedova,Luca Denti,Xavier Didelot,Davide Monti,Yuri Pirola
DOI: https://doi.org/10.1093/bioinformatics/btae292
IF: 5.8
2024-04-27
Bioinformatics
Abstract:Abstract Motivation Bacterial genomes present more variability than human genomes, which requires important adjustments in computational tools that are developed for human data. In particular, bacteria exhibit a mosaic structure due to homologous recombinations, but this fact is not sufficiently captured by standard read mappers that align against linear reference genomes. The recent introduction of pangenomics provides some insights in that context, as a pangenome graph can represent the variability within a species. However, the concept of sequence-to-graph alignment that captures the presence of recombinations has not been previously investigated. Results In this paper, we present the extension of the notion of sequence-to-graph alignment to a variation graph that incorporates a recombination, so that the latter are explicitly represented and evaluated in an alignment. Moreover, we present a dynamic programming approach for the special case where there is at most a recombinationâwe implement this case as RecGraph. From a modeling point of view, a recombination corresponds to identifying a new path of the variation graph, where the new arc is composed of two halves, each extracted from an original path, possibly joined by a new arc. Our experiments show that RecGraph accurately aligns simulated recombinant bacterial sequences that have at most a recombination, providing evidence for the presence of recombination events. Availability Our implementation is open source and available at https://github.com/AlgoLab/RecGraph Supplementary information Supplementary data are available at Bioinformatics online.
biochemical research methods,biotechnology & applied microbiology,mathematical & computational biology