Spatial heterogeneity plays an important role in the evolution of drug resistance. While recent studies have indicated that spatial gradients of selection pressure can accelerate resistance evolution, much less is known about evolution in more complex spatial profiles. Here we use a stochastic toy model of drug resistance to investigate how different spatial profiles of selection pressure impact the time to fixation of a resistant allele. Using mean first passage time calculations, we show that spatial heterogeneity accelerates resistance evolution when the rate of spatial migration far exceeds that of mutation but slows fixation when mutation dominates. Interestingly, there exists an intermediate regime-characterized by comparable rates of migration and mutation-in which the rate of fixation can be either accelerated or decelerated depending on the spatial profile, even when spatially averaged selection pressure remains constant. Finally, we demonstrate that optimal tuning of the spatial profile can dramatically slow the spread and fixation of resistant subpopulations, which may lay the groundwork for optimized, spatially-resolved drug dosing strategies for mitigating the effects of drug resistance.Drug resistance is a rapidly growing public health threat and a central impediment to the treatment of cancer, viruses, and microbial infections [1][2][3][4]. The battle against resistance has been largely fought at the molecular level, leading to an increasingly mature understanding of its underlying biochemical and genetic roots. At the same time, evolutionary biologists have long recognized resistance as a fundamentally stochastic process governed by the complex interplay between microbial ecology and evolutionary selection. The last decade, in particular, has seen a significant surge in efforts to develop and understand evolution-based treatment strategies to forestall resistance [5][6][7][8][9][10][11][12][13][14][15][16]. While the vast majority of this work focuses on spatially homogeneous environments, a number of recent studies, both theoretical and experimental, have demonstrated that spatial heterogeneity in drug concentration can dramatically alter the evolutionary dynamics leading to resistance [16][17][18][19][20][21][22][23][24]. On a practical level, the picture that emerges is somewhat bleak, as resistance evolution is dramatically accelerated in the presence of spatial gradients in drug concentration [18][19][20][22][23][24] or heterogeneous drug penetration [17,21]. Interestingly, however, recent work shows that this acceleration can be tempered or even reversed when the mutational pathway (i.e. the genotypic fitness landscape) leading to resistance contains fitness valleys [18], which are known to inhibit evolution [25][26][27][28]. Unfortunately, because the fitness landscape is a genetic property of the cells themselves, the potential for accelerated evolution appears to be "built in", making it difficult to combat in a treatment setting. However, these results raise the question of whether non-...