Evolutionary rescue of resistant mutants is governed by a balance between radial expansion and selection in compact populations

Selective sweep probabilities in spatially expanding populations

Evolution during range expansions shapes biological systems from microbial communities and tumours up to invasive species. A fundamental question is whether, when a beneficial mutation arises during a range expansion, it will evade clonal interference and sweep through the population to fixation. However, most theoretical investigations of range expansions have been confined to regimes in which selective sweeps are effectively impossible, while studies of selective sweeps have either assumed constant population size or have ignored spatial structure. Here we use mathematical modelling and analysis to investigate selective sweep probabilities in the alternative yet biologically relevant scenario in which mutants can outcompete and displace a slowly spreading wildtype. Assuming constant radial expansion speed, we derive probability distributions for the arrival time and location of the first surviving mutant and hence find surprisingly simple approximate and exact expressions for selective sweep probabilities in one, two and three dimensions, which are independent of mutation rate. Namely, the selective sweep probability is approximately 1 - c w t / c m d , where c w t and c m are the wildtype and mutant radial expansion speeds, and d the spatial dimension. Using agent-based simulations, we show that our analytical results accurately predict selective sweep frequencies in the two-dimensional spatial Moran process. We further compare our results with those obtained for alternative growth laws. Parameterizing our model for human tumours, we find that selective sweeps are predicted to be rare except during very early solid tumour growth, thus providing a general, pan-cancer explanation for findings from recent sequencing studies.

A seven-step guide to spatial, agent-based modelling of tumour evolution

Spatial agent-based models are frequently used to investigate the evolution of solid tumours subject to localized cell-cell interactions and microenvironmental heterogeneity. As spatial genomic, transcriptomic and proteomic technologies gain traction, spatial computational models are predicted to become ever more necessary for making sense of complex clinical and experimental data sets, for predicting clinical outcomes, and for optimizing treatment strategies. Here we present a non-technical step by step guide to developing such a model from first principles. Stressing the importance of tailoring the model structure to that of the biological system, we describe methods of increasing complexity, from the basic Eden growth model up to off-lattice simulations with diffusible factors. We examine choices that unavoidably arise in model design, such as implementation, parameterization, visualization and reproducibility. Each topic is illustrated with examples drawn from recent research studies and state of the art modelling platforms. We emphasize the benefits of simpler models that aim to match the complexity of the phenomena of interest, rather than that of the entire biological system. Our guide is aimed at both aspiring modellers and other biologists and oncologists who wish to understand the assumptions and limitations of the models on which major cancer studies now so often depend.

Quantifying stochastic establishment of mutants in microbial adaptation

Studies of microbial evolution, especially in applied contexts, have focused on the role of selection in shaping predictable, adaptive responses to the environment. However, chance events - the appearance of novel genetic variants and their establishment, i.e. outgrowth from a single cell to a sizeable population - also play critical initiating roles in adaptation. Stochasticity in establishment has received little attention in microbiology, potentially due to lack of awareness as well as practical challenges in quantification. However, methods for high-replicate culturing, mutant labelling and detection, and statistical inference now make it feasible to experimentally quantify the establishment probability of specific adaptive genotypes. I review methods that have emerged over the past decade, including experimental design and mathematical formulas to estimate establishment probability from data. Quantifying establishment in further biological settings and comparing empirical estimates to theoretical predictions represent exciting future directions. More broadly, recognition that adaptive genotypes may be stochastically lost while rare is significant both for interpreting common lab assays and for designing interventions to promote or inhibit microbial evolution.