PAN-GWES: Pangenome-spanning epistasis and co-selection analysis via de Bruijn graphs

Abstract

Studies of bacterial adaptation and evolution are hampered by the difficulty of measuring traits such as virulence, drug resistance and transmissibility in large populations. In contrast, it is now feasible to obtain high-quality complete assemblies of many bacterial genomes thanks to scalable and affordable long-read sequencing technology. To exploit this opportunity we introduce a phenotype- and alignment-free method for discovering co-selected and epistatically interacting genomic variation from genome assemblies covering both core and accessory parts of genomes. Our approach uses a compact coloured de Bruijn graph to approximate the intra-genome distances between pairs of loci for a collection of bacterial genomes to account for the impacts of linkage disequilibrium (LD). We demonstrate the versatility of our approach to efficiently identify associations between loci linked with drug resistance and adaptation to the hospital niche in the major human bacterial pathogens Streptococcus pneumoniae and Enterococcus faecalis.