Genomic constraints shape the evolution of alternative routes to drug resistance in prokaryotes.

BackgroundVariation within the prokaryotic pangenome is not random, and natural selection that favours particular combinations of genes appears to dominate over random drift. What is less clear is whether similar phenotypes tend to arise independently within pangenomes, or whether there is a single route through which they evolve. We have previously shown that the wider genomic context is vital to understanding specific phenotypes, such as antimicrobial resistance. Considering the amount of genomic variation we see in many pathogens, the question of routes to resistance and how they arise is particularly urgent. MethodsWe hypothesised that mutually exclusive routes to multidrug antimicrobial resistance evolve through distinct yet predictable pathways dependent upon the species-specific pangenomic background. To test this, we developed an integrated pangenomic and machine learning framework to analyse the emergence of multidrug resistance (MDR) in two clinically significant pathogens, Pseudomonas aeruginosa and Escherichia coli. Machine learning models uncovered genes statistically linked to particular MDR phenotypes. We then examined whether these genes were significantly associated or dissociated with one another in gene co-occurrence networks and evaluated the importance of genomic context on pairs of genes. ResultsWe demonstrate three key findings on the evolution of multidrug resistance in our dataset. First, antimicrobial resistance genes (ARGs) were distributed across both the core and accessory genomes, challenging the prevailing view that ARGs are typically confined to the accessory genome. Second, we demonstrate that E. coli possesses mutually exclusive pathways to MDR, with several resistance genes showing distinct dissociation patterns, suggesting alternative evolutionary routes to identical MDR phenotypes. Third, the network of ARG coassociations differed significantly between the P. aeruginosa and E. coli pangenomes with 33 gene pairs showing opposite association patterns between species, revealing species-specific genomic constraints on MDR emergence. These findings emphasise that bacterial evolution is more constrained than previously thought and that designing intervention therapies requires consideration of both historical evolutionary trajectories and the specific genomic context in which the resistance phenotypes arose.