Modeling spatial evolution of multi-drug resistance under drug environmental gradients

Multi-drug combinations to treat bacterial populations are at the forefront of approaches for infection control and prevention of antibiotic resistance. Although the evolution of antibiotic resistance has been theoretically studied with mathematical population dynamics models, extensions to spatial dynamics remain rare in the literature, including in particular spatial evolution of multi-drug resistance. In this study, we propose a reaction-diffusion system that describes the multi-drug evolution of bacteria based on a drug-concentration rescaling approach. We show how the resistance to drugs in space, and the consequent adaptation of growth rate, is governed by a Price equation with diffusion, integrating features of drug interactions and collateral resistances or sensitivities to the drugs. We study spatial versions of the model where the distribution of drugs is homogeneous across space, and where the drugs vary environmentally in a piecewise-constant, linear and nonlinear manner. Although in many evolution models, per capita growth rate is a natural surrogate for fitness, in spatially-extended, potentially heterogeneous habitats, fitness is an emergent property that potentially reflects additional complexities, from boundary conditions to the specific spatial variation of growth rates. Applying concepts from perturbation theory and reaction-diffusion equations, we propose an analytical metric for characterization of average mutant fitness in the spatial system based on the principal eigenvalue of our linear problem, lambda 1. This enables an accurate translation from drug spatial gradients and mutant antibiotic susceptibility traits to the relative advantage of each mutant across the environment. Our approach allows one to predict the precise outcomes of selection among mutants over space, ultimately from comparing their lambda 1 values, which encode a critical interplay between growth functions, movement traits, habitat size and boundary conditions. Such mathematical understanding opens new avenues for multi-drug therapeutic optimization. In this study we develop a framework to model multi-drug resistance evolution in space by combining drug-rescaling arguments with a reaction-diffusion type model. In response to multi-drug environmental gradients, each independent mutant can grow and diffuse following an individual spatially-varying growth function and diffusion rate. Applying concepts from perturbation theory and reaction-diffusion models, we propose an analytical metric to quantify average mutant fitness in the spatial system and to predict the outcome of selection. Our findings highlight that in spatially-extended habitats fitness is an emergent property that potentially integrates many complexities, from boundary conditions to environmental variation, as well as individual growth and diffusion traits.