Diverse mutant selection windows shape spatial heterogeneity in evolving populations

Mutant selection windows (MSWs), the range of drug concentrations that select for drug-resistant mutants, have long been used as a model for predicting drug resistance and designing optimal dosing strategies in infectious disease. The canonical MSW model offers comparisons between two subtypes at a time: drug-sensitive and drug-resistant. In contrast, the fitness landscape model with N alleles, which maps genotype to fitness, allows comparisons between N genotypes simultaneously, but does not encode continuous drug response data. In clinical settings, there may be a wide range of drug concentrations selecting for a variety of genotypes. Therefore, there is a need for a more robust model of the pathogen response to therapy to predict resistance and design new therapeutic approaches. Fitness seascapes, which model genotype-by-environment interactions, permit multiple MSW comparisons simultaneously by encoding genotype-specific dose-response data. By comparing dose-response curves, one can visualize the range of drug concentrations where one genotype is selected over another. In this work, we show how N-allele fitness seascapes allow for N ∗2N−1 unique MSW comparisons. In spatial drug diffusion models, we demonstrate how fitness seascapes reveal spatially heterogeneous MSWs, extending the MSW model to more accurately reflect the selection fo drug resistant genotypes. Furthermore, we find that the spatial structure of MSWs shapes the evolution of drug resistance in an agent-based model. Our work highlights the importance and utility of considering dose-dependent fitness seascapes in evolutionary medicine. Author Summary Drug resistance in infectious disease and cancer is a major driver of mortality. While undergoing treatment, the population of cells in a tumor or infection may evolve the ability to grow despite the use of previously effective drugs. Researchers hypothesize that the spatial organization of these disease populations may contribute to drug resistance. In this work, we analyze how spatial gradients of drug concentration impact the evolution of drug resistance. We consider a decades-old model called the mutant selection window (MSW), which describes the drug concentration range that selects for drug-resistant cells. We show how extending this model with continuous dose-response data, which describes how different types of cells respond to drug, improves the ability of MSWs to predict evolution. This work helps us understand how the spatial organization of cells, such as the organization of blood vessels within a tumor, may promote drug resistance. In the future, we may use these methods to optimize drug dosing to prevent resistance or leverage known vulnerabilities of drug-resistant cells.