Mutational pathway determines whether drug gradients accelerate evolution of drug-resistant cells

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, λ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 λ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.

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 rescaling approach (Gjini and Wood, 2021). We show how the resistance to drugs in space, and the consequent adaptation of growth rate is governed by a Price equation with diffusion. The covariance terms in this equation integrate 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. Applying concepts from perturbation theory and reaction-diffusion equations, we propose an analytical characterization of average mutant fitness in the spatial system based on the principal eigenvalue of our linear problem. This enables an accurate translation from drug spatial gradients and mutant antibiotic susceptibility traits, to the relative advantage of each mutant across the environment. Such a mathematical understanding allows to predict the precise outcomes of selection over space, ultimately from the fundamental balance between growth and movement traits, and their diversity in a population.

Assisted percolation of slow-spreading mutants in heterogeneous environments

Environmental heterogeneity can drive genetic heterogeneity in expanding populations; mutant strains may emerge that trade overall growth rate for an improved ability to survive in patches that are hostile to the wild type. This evolutionary dynamic is of practical importance when seeking to prevent the emergence of damaging traits. We show that a subcritical slow-spreading mutant can attain dominance even when the density of patches is below their percolation threshold and predict this transition using geometrical arguments. This work demonstrates a phenomenon of "assisted percolation", where one subcritical process assists another to achieve supercriticality.

Bacterial motility can govern the dynamics of antibiotic resistance evolution

Spatial heterogeneity in antibiotic concentrations is thought to accelerate the evolution of antibiotic resistance, but current theory and experiments have overlooked the effect of cell motility on bacterial adaptation. Here, we study bacterial evolution in antibiotic landscapes with a quantitative model where bacteria evolve under the stochastic processes of proliferation, death, mutation and migration. Numerical and analytical results show that cell motility can both accelerate and decelerate bacterial adaptation by affecting the degree of genotypic mixing and ecological competition. Moreover, we find that for sufficiently high rates, cell motility can limit bacterial survival, and we derive conditions for all these regimes. Similar patterns are observed in more complex scenarios, namely where bacteria can bias their motion in chemical gradients (chemotaxis) or switch between motility phenotypes either stochastically or in a density-dependent manner. Overall, our work reveals limits to bacterial adaptation in antibiotic landscapes that are set by cell motility.

Model for Quorum-Sensing Mediated Stochastic Biofilm Nucleation

Surface-attached bacterial biofilms cause disease and industrial biofouling, as well as being widespread in the natural environment. Density-dependent quorum sensing is one of the mechanisms implicated in biofilm initiation. Here we present and analyze a model for quorum-sensing triggered biofilm initiation. In our model, individual, planktonic bacteria adhere to a surface, proliferate, and undergo a collective transition to a biofilm phenotype. This model predicts a stochastic transition between a loosely attached, finite layer of bacteria near the surface and a growing biofilm. The transition is governed by two key parameters: the collective transition density relative to the carrying capacity and the immigration rate relative to the detachment rate. Biofilm initiation is complex, but our model suggests that stochastic nucleation phenomena may be relevant.

A computational model for microbial colonization of an antifouling surface

Biofouling of marine surfaces such as ship hulls is a major industrial problem. Antifouling (AF) paints delay the onset of biofouling by releasing biocidal chemicals. We present a computational model for microbial colonization of a biocide-releasing AF surface. Our model accounts for random arrival from the ocean of microorganisms with different biocide resistance levels, biocide-dependent proliferation or killing, and a transition to a biofilm state. Our computer simulations support a picture in which biocide-resistant microorganisms initially form a loosely attached layer that eventually transitions to a growing biofilm. Once the growing biofilm is established, immigrating microorganisms are shielded from the biocide, allowing more biocide-susceptible strains to proliferate. In our model, colonization of the AF surface is highly stochastic. The waiting time before the biofilm establishes is exponentially distributed, suggesting a Poisson process. The waiting time depends exponentially on both the concentration of biocide at the surface and the rate of arrival of resistant microorganisms from the ocean. Taken together our results suggest that biofouling of AF surfaces may be intrinsically stochastic and hence unpredictable, but immigration of more biocide-resistant species, as well as the biological transition to biofilm physiology, may be important factors controlling the time to biofilm establishment.

Lineage dynamics in growing biofilms: Spatial patterns of standing vs. de novo diversity

Microbial biofilms show high phenotypic and genetic diversity, yet the mechanisms underlying diversity generation and maintenance remain unclear. Here, we investigate how spatial patterns of growth activity within a biofilm lead to spatial patterns of genetic diversity. Using individual-based computer simulations, we show that the active layer of growing cells at the biofilm interface controls the distribution of lineages within the biofilm, and therefore the patterns of standing and de novo diversity. Comparing biofilms of equal size, those with a thick active layer retain more standing diversity, while de novo diversity is more evenly distributed within the biofilm. In contrast, equal-sized biofilms with a thin active layer retain less standing diversity, and their de novo diversity is concentrated at the top of the biofilm, and in fewer lineages. In the context of antimicrobial resistance, biofilms with a thin active layer may be more prone to generate lineages with multiple resistance mutations, and to seed new resistant biofilms via sloughing of resistant cells from the upper layers. Our study reveals fundamental "baseline" mechanisms underlying the patterning of diversity within biofilms.

Emergence of Resistant Escherichia coli Mutants in Microfluidic On-Chip Antibiotic Gradients

Spatiotemporal structures and heterogeneities are common in natural habitats, yet their role in the evolution of antibiotic resistance is still to be uncovered. We applied a microfluidic gradient generator device to study the emergence of resistant bacteria in spatial ciprofloxacin gradients. We observed biofilm formation in regions with sub-inhibitory concentrations of antibiotics, which quickly expanded into the high antibiotic regions. In the absence of an explicit structure of the habitat, this multicellular formation led to a spatial structure of the population with local competition and limited migration. Therefore, such structures can function as amplifiers of selection and aid the spread of beneficial mutations. We found that the physical environment itself induces stress-related mutations that later prove beneficial when cells are exposed to antibiotics. This shift in function suggests that exaptation occurs in such experimental scenarios. The above two processes pave the way for the subsequent emergence of highly resistant specific mutations.

Hydrodynamic flow and concentration gradients in the gut enhance neutral bacterial diversity

The gut microbiota features important genetic diversity, and the specific spatial features of the gut may shape evolution within this environment. We investigate the fixation probability of neutral bacterial mutants within a minimal model of the gut that includes hydrodynamic flow and resulting gradients of food and bacterial concentrations. We find that this fixation probability is substantially increased, compared with an equivalent well-mixed system, in the regime where the profiles of food and bacterial concentration are strongly spatially dependent. Fixation probability then becomes independent of total population size. We show that our results can be rationalized by introducing an active population, which consists of those bacteria that are actively consuming food and dividing. The active population size yields an effective population size for neutral mutant fixation probability in the gut.

Lanza V, Rodríguez-Beltrán J, Galán JC, San Millán A, Cantón R, Coque TM. Evolutionary Pathways and Trajectories in Antibiotic Resistance

Evolution is the hallmark of life. Descriptions of the evolution of microorganisms have provided a wealth of information, but knowledge regarding "what happened" has precluded a deeper understanding of "how" evolution has proceeded, as in the case of antimicrobial resistance. The difficulty in answering the "how" question lies in the multihierarchical dimensions of evolutionary processes, nested in complex networks, encompassing all units of selection, from genes to communities and ecosystems. At the simplest ontological level (as resistance genes), evolution proceeds by random (mutation and drift) and directional (natural selection) processes; however, sequential pathways of adaptive variation can occasionally be observed, and under fixed circumstances (particular fitness landscapes), evolution is predictable. At the highest level (such as that of plasmids, clones, species, microbiotas), the systems' degrees of freedom increase dramatically, related to the variable dispersal, fragmentation, relatedness, or coalescence of bacterial populations, depending on heterogeneous and changing niches and selective gradients in complex environments. Evolutionary trajectories of antibiotic resistance find their way in these changing landscapes subjected to random variations, becoming highly entropic and therefore unpredictable. However, experimental, phylogenetic, and ecogenetic analyses reveal preferential frequented paths (highways) where antibiotic resistance flows and propagates, allowing some understanding of evolutionary dynamics, modeling and designing interventions. Studies on antibiotic resistance have an applied aspect in improving individual health, One Health, and Global Health, as well as an academic value for understanding evolution. Most importantly, they have a heuristic significance as a model to reduce the negative influence of anthropogenic effects on the environment.

Spatial segregation and cooperation in radially expanding microbial colonies under antibiotic stress

Antibiotic resistance in microbial communities reflects a combination of processes operating at different scales. In this work, we investigate the spatiotemporal dynamics of bacterial colonies comprised of drug-resistant and drug-sensitive cells undergoing range expansion under antibiotic stress. Using the opportunistic pathogen Enterococcus faecalis with plasmid-encoded β-lactamase, we track colony expansion dynamics and visualize spatial patterns in fluorescently labeled populations exposed to antibiotics. We find that the radial expansion rate of mixed communities is approximately constant over a wide range of drug concentrations and initial population compositions. Imaging of the final populations shows that resistance to ampicillin is cooperative, with sensitive cells surviving in the presence of resistant cells at otherwise lethal concentrations. The populations exhibit a diverse range of spatial segregation patterns that depend on drug concentration and initial conditions. Mathematical models indicate that the observed dynamics are consistent with global cooperation, despite the fact that β-lactamase remains cell-associated. Experiments confirm that resistant colonies provide a protective effect to sensitive cells on length scales multiple times the size of a single colony, and populations seeded with (on average) no more than a single resistant cell can produce mixed communities in the presence of the drug. While biophysical models of drug degradation suggest that individual resistant cells offer only short-range protection to neighboring cells, we show that long-range protection may arise from synergistic effects of multiple resistant cells, providing surprisingly large protection zones even at small population fractions.

Price equation captures the role of drug interactions and collateral effects in the evolution of multidrug resistance

Bacterial adaptation to antibiotic combinations depends on the joint inhibitory effects of the two drugs (drug interaction [DI]) and how resistance to one drug impacts resistance to the other (collateral effects [CE]). Here we model these evolutionary dynamics on two-dimensional phenotype spaces that leverage scaling relations between the drug-response surfaces of drug-sensitive (ancestral) and drug-resistant (mutant) populations. We show that evolved resistance to the component drugs - and in turn, the adaptation of growth rate - is governed by a Price equation whose covariance terms encode geometric features of both the two-drug-response surface (DI) in ancestral cells and the correlations between resistance levels to those drugs (CE). Within this framework, mean evolutionary trajectories reduce to a type of weighted gradient dynamics, with the drug interaction dictating the shape of the underlying landscape and the collateral effects constraining the motion on those landscapes. We also demonstrate how constraints on available mutational pathways can be incorporated into the framework, adding a third key driver of evolution. Our results clarify the complex relationship between drug interactions and collateral effects in multidrug environments and illustrate how specific dosage combinations can shift the weighting of these two effects, leading to different and temporally explicit selective outcomes.

Evolution in alternating environments with tunable interlandscape correlations

Natural populations are often exposed to temporally varying environments. Evolutionary dynamics in varying environments have been extensively studied, although understanding the effects of varying selection pressures remains challenging. Here, we investigate how cycling between a pair of statistically related fitness landscapes affects the evolved fitness of an asexually reproducing population. We construct pairs of fitness landscapes that share global fitness features but are correlated with one another in a tunable way, resulting in landscape pairs with specific correlations. We find that switching between these landscape pairs, depending on the ruggedness of the landscape and the interlandscape correlation, can either increase or decrease steady-state fitness relative to evolution in single environments. In addition, we show that switching between rugged landscapes often selects for increased fitness in both landscapes, even in situations where the landscapes themselves are anticorrelated. We demonstrate that positively correlated landscapes often possess a shared maximum in both landscapes that allows the population to step through sub-optimal local fitness maxima that often trap single landscape evolution trajectories. Finally, we demonstrate that switching between anticorrelated paired landscapes leads to ergodic-like dynamics where each genotype is populated with nonzero probability, dramatically lowering the steady-state fitness in comparison to single landscape evolution.

The Moran process on 2-chromatic graphs

Resources are rarely distributed uniformly within a population. Heterogeneity in the concentration of a drug, the quality of breeding sites, or wealth can all affect evolutionary dynamics. In this study, we represent a collection of properties affecting the fitness at a given location using a color. A green node is rich in resources while a red node is poorer. More colors can represent a broader spectrum of resource qualities. For a population evolving according to the birth-death Moran model, the first question we address is which structures, identified by graph connectivity and graph coloring, are evolutionarily equivalent. We prove that all properly two-colored, undirected, regular graphs are evolutionarily equivalent (where “properly colored” means that no two neighbors have the same color). We then compare the effects of background heterogeneity on properly two-colored graphs to those with alternative schemes in which the colors are permuted. Finally, we discuss dynamic coloring as a model for spatiotemporal resource fluctuations, and we illustrate that random dynamic colorings often diminish the effects of background heterogeneity relative to a proper two-coloring.

Heterogeneity in environmental conditions can have profound effects on long-term evolutionary outcomes in structured populations. We consider a population evolving on a colored graph, wherein the color of a node represents the resources at that location. Using a combination of analytical and numerical methods, we quantify the effects of background heterogeneity on a population’s dynamics. In addition to considering the notion of an “optimal” coloring with respect to mutant invasion, we also study the effects of dynamic spatial redistribution of resources as the population evolves. Although the effects of static background heterogeneity can be quite striking, these effects are often attenuated by the movement (or “flow”) of the underlying resources.

Adapt or Perish: Evolutionary Rescue in a Gradually Deteriorating Environment

We investigate the evolutionary rescue of a microbial population in a gradually deteriorating environment, through a combination of analytical calculations and stochastic simulations. We consider a population destined for extinction in the absence of mutants, which can survive only if mutants sufficiently adapted to the new environment arise and fix. We show that mutants that appear later during the environment deterioration have a higher probability to fix. The rescue probability of the population increases with a sigmoidal shape when the product of the carrying capacity and of the mutation probability increases. Furthermore, we find that rescue becomes more likely for smaller population sizes and/or mutation probabilities if the environment degradation is slower, which illustrates the key impact of the rapidity of environment degradation on the fate of a population. We also show that our main conclusions are robust across various types of adaptive mutants, including specialist and generalist ones, as well as mutants modeling antimicrobial resistance evolution. We further express the average time of appearance of the mutants that do rescue the population and the average extinction time of those that do not. Our methods can be applied to other situations with continuously variable fitnesses and population sizes, and our analytical predictions are valid in the weak-to-moderate mutation regime.

A Roadblock-and-Kill Mechanism of Action Model for the DNA-Targeting Antibiotic Ciprofloxacin

Fluoroquinolones, antibiotics that cause DNA damage by inhibiting DNA topoisomerases, are clinically important, but their mechanism of action is not yet fully understood. In particular, the dynamical response of bacterial cells to fluoroquinolone exposure has hardly been investigated, although the SOS response, triggered by DNA damage, is often thought to play a key role. Here, we investigated the growth inhibition of the bacterium Escherichia coli by the fluoroquinolone ciprofloxacin at low concentrations. We measured the long-term and short-term dynamical response of the growth rate and DNA production rate to ciprofloxacin at both the population and single-cell levels. We show that, despite the molecular complexity of DNA metabolism, a simple roadblock-and-kill model focusing on replication fork blockage and DNA damage by ciprofloxacin-poisoned DNA topoisomerase II (gyrase) quantitatively reproduces long-term growth rates in the presence of ciprofloxacin. The model also predicts dynamical changes in the DNA production rate in wild-type E. coli and in a recombination-deficient mutant following a step-up of ciprofloxacin. Our work highlights that bacterial cells show a delayed growth rate response following fluoroquinolone exposure. Most importantly, our model explains why the response is delayed: it takes many doubling times to fragment the DNA sufficiently to inhibit gene expression. We also show that the dynamical response is controlled by the timescale of DNA replication and gyrase binding/unbinding to the DNA rather than by the SOS response, challenging the accepted view. Our work highlights the importance of including detailed biophysical processes in biochemical-systems models to quantitatively predict the bacterial response to antibiotics.

Population structure across scales facilitates coexistence and spatial heterogeneity of antibiotic-resistant infections

Antibiotic-resistant infections are a growing threat to human health, but basic features of the eco-evolutionary dynamics remain unexplained. Most prominently, there is no clear mechanism for the long-term coexistence of both drug-sensitive and resistant strains at intermediate levels, a ubiquitous pattern seen in surveillance data. Here we show that accounting for structured or spatially-heterogeneous host populations and variability in antibiotic consumption can lead to persistent coexistence over a wide range of treatment coverages, drug efficacies, costs of resistance, and mixing patterns. Moreover, this mechanism can explain other puzzling spatiotemporal features of drug-resistance epidemiology that have received less attention, such as large differences in the prevalence of resistance between geographical regions with similar antibiotic consumption or that neighbor one another. We find that the same amount of antibiotic use can lead to very different levels of resistance depending on how treatment is distributed in a transmission network. We also identify parameter regimes in which population structure alone cannot support coexistence, suggesting the need for other mechanisms to explain the epidemiology of antibiotic resistance. Our analysis identifies key features of host population structure that can be used to assess resistance risk and highlights the need to include spatial or demographic heterogeneity in models to guide resistance management.

Using Selection by Nonantibiotic Stressors to Sensitize Bacteria to Antibiotics

Evolutionary adaptation of bacteria to nonantibiotic selective forces, such as osmotic stress, has been previously associated with increased antibiotic resistance, but much less is known about potentially sensitizing effects of nonantibiotic stressors. In this study, we use laboratory evolution to investigate adaptation of Enterococcus faecalis, an opportunistic bacterial pathogen, to a broad collection of environmental agents, ranging from antibiotics and biocides to extreme pH and osmotic stress. We find that nonantibiotic selection frequently leads to increased sensitivity to other conditions, including multiple antibiotics. Using population sequencing and whole-genome sequencing of single isolates from the evolved populations, we identify multiple mutations in genes previously linked with resistance to the selecting conditions, including genes corresponding to known drug targets or multidrug efflux systems previously tied to collateral sensitivity. Finally, we hypothesized based on the measured sensitivity profiles that sequential rounds of antibiotic and nonantibiotic selection may lead to hypersensitive populations by harnessing the orthogonal collateral effects of particular pairs of selective forces. To test this hypothesis, we show experimentally that populations evolved to a sequence of linezolid (an oxazolidinone antibiotic) and sodium benzoate (a common preservative) exhibit increased sensitivity to more stressors than adaptation to either condition alone. The results demonstrate how sequential adaptation to drug and nondrug environments can be used to sensitize bacteria to antibiotics and highlight new potential strategies for exploiting shared constraints governing adaptation to diverse environmental challenges.

Resist or perish: Fate of a microbial population subjected to a periodic presence of antimicrobial

The evolution of antimicrobial resistance can be strongly affected by variations of antimicrobial concentration. Here, we study the impact of periodic alternations of absence and presence of antimicrobial on resistance evolution in a microbial population, using a stochastic model that includes variations of both population composition and size, and fully incorporates stochastic population extinctions. We show that fast alternations of presence and absence of antimicrobial are inefficient to eradicate the microbial population and strongly favor the establishment of resistance, unless the antimicrobial increases enough the death rate. We further demonstrate that if the period of alternations is longer than a threshold value, the microbial population goes extinct upon the first addition of antimicrobial, if it is not rescued by resistance. We express the probability that the population is eradicated upon the first addition of antimicrobial, assuming rare mutations. Rescue by resistance can happen either if resistant mutants preexist, or if they appear after antimicrobial is added to the environment. Importantly, the latter case is fully prevented by perfect biostatic antimicrobials that completely stop division of sensitive microorganisms. By contrast, we show that the parameter regime where treatment is efficient is larger for biocidal drugs than for biostatic drugs. This sheds light on the respective merits of different antimicrobial modes of action.

Delayed antibiotic exposure induces population collapse in enterococcal communities with drug-resistant subpopulations

The molecular underpinnings of antibiotic resistance are increasingly understood, but less is known about how these molecular events influence microbial dynamics on the population scale. Here, we show that the dynamics of E. faecalis communities exposed to antibiotics can be surprisingly rich, revealing scenarios where increasing population size or delaying drug exposure can promote population collapse. Specifically, we demonstrate how density-dependent feedback loops couple population growth and antibiotic efficacy when communities include drug-resistant subpopulations, leading to a wide range of behavior, including population survival, collapse, or one of two qualitatively distinct bistable behaviors where survival is favored in either small or large populations. These dynamics reflect competing density-dependent effects of different subpopulations, with growth of drug-sensitive cells increasing but growth of drug-resistant cells decreasing effective drug inhibition. Finally, we demonstrate how populations receiving immediate drug influx may sometimes thrive, while identical populations exposed to delayed drug influx collapse.

Antibiotic interactions shape short-term evolution of resistance in E. faecalis

Antibiotic combinations are increasingly used to combat bacterial infections. Multidrug therapies are a particularly important treatment option for E. faecalis, an opportunistic pathogen that contributes to high-inoculum infections such as infective endocarditis. While numerous synergistic drug combinations for E. faecalis have been identified, much less is known about how different combinations impact the rate of resistance evolution. In this work, we use high-throughput laboratory evolution experiments to quantify adaptation in growth rate and drug resistance of E. faecalis exposed to drug combinations exhibiting different classes of interactions, ranging from synergistic to suppressive. We identify a wide range of evolutionary behavior, including both increased and decreased rates of growth adaptation, depending on the specific interplay between drug interaction and drug resistance profiles. For example, selection in a dual β-lactam combination leads to accelerated growth adaptation compared to selection with the individual drugs, even though the resulting resistance profiles are nearly identical. On the other hand, populations evolved in an aminoglycoside and β-lactam combination exhibit decreased growth adaptation and resistant profiles that depend on the specific drug concentrations. We show that the main qualitative features of these evolutionary trajectories can be explained by simple rescaling arguments that correspond to geometric transformations of the two-drug growth response surfaces measured in ancestral cells. The analysis also reveals multiple examples where resistance profiles selected by drug combinations are nearly growth-optimized along a contour connecting profiles selected by the component drugs. Our results highlight trade-offs between drug interactions and resistance profiles during the evolution of multi-drug resistance and emphasize evolutionary benefits and disadvantages of particular drug pairs targeting enterococci.

The effect of spatiotemporal antibiotic inhomogeneities on the evolution of resistance

Combating the evolution of widespread antibiotic resistance is one of the most pressing challenges facing modern medicine. Recent research has demonstrated that the evolution of pathogens with high levels of resistance can be accelerated by spatial and temporal inhomogeneities in antibiotic concentration, which frequently arise in patients and the environment. Strategies to predict and counteract the effects of such inhomogeneities will be critical in the fight against resistance. In this paper we develop a mechanistic framework for modelling the adaptive evolution of resistance in the presence of spatiotemporal antibiotic concentrations, which treats the adaptive process as an interaction between two mutually orthogonal forces; the first returns cells to their wild-type state in the absence of antibiotic selection, and the second selects for increased coping ability in the presence of an antibiotic. We apply our model to investigate laboratory adaptation experiments, and then extend it to consider the case in which multiple strategies for resistance undergo competitive evolution.

Diffusion-driven enhancement of the antibiotic resistance selection window

The current crisis of antimicrobial resistance in clinically relevant pathogens has highlighted our limited understanding of the ecological and evolutionary forces that drive drug resistance adaptation. For instance, although human tissues are highly heterogeneous, most of our mechanistic understanding about antibiotic resistance evolution is based on constant and well-mixed environmental conditions. A consequence of considering spatial heterogeneity is that, even if antibiotics are prescribed at high dosages, the penetration of drug molecules through tissues inevitably produces antibiotic gradients, exposing bacterial populations to a range of selective pressures and generating a dynamic fitness landscape that changes in space and time. In this paper, we will use a combination of mathematical modelling and computer simulations to study the population dynamics of susceptible and resistant strains competing for resources in a network of micro-environments with varying degrees of connectivity. Our main result is that highly connected environments increase diffusion of drug molecules, enabling resistant phenotypes to colonize a larger number of spatial locations. We validated this theoretical result by culturing fluorescently labelled Escherichia coli in 3D-printed devices that allow us to control the rate of diffusion of antibiotics between neighbouring compartments and quantify the spatio-temporal distribution of resistant and susceptible bacterial cells.

Growth-dependent drug susceptibility can prevent or enhance spatial expansion of a bacterial population

As a population wave expands, organisms at the tip typically experience plentiful nutrients while those behind the front become nutrient-depleted. If the environment also contains a gradient of some inhibitor (e.g. a toxic drug), a tradeoff exists: the nutrient-rich tip is more exposed to the inhibitor, while the nutrient-starved region behind the front is less exposed. Here we show that this can lead to complex dynamics when the organism's response to the inhibitory substance is coupled to nutrient availability. We model a bacterial population which expands in a spatial gradient of antibiotic, under conditions where either fast-growing bacteria at the wave's tip, or slow-growing, resource-limited bacteria behind the front are more susceptible to the antibiotic. We find that growth-rate dependent susceptibility can have strong effects on the dynamics of the expanding population. If slow-growing bacteria are more susceptible, the population wave advances far into the inhibitory zone, leaving a trail of dead bacteria in its wake. In contrast, if fast-growing bacteria are more susceptible, the wave is blocked at a much lower concentration of antibiotic, but a large population of live bacteria remains behind the front. Our results may contribute to understanding the efficacy of different antimicrobials for spatially structured microbial populations such as biofilms, as well as the dynamics of ecological population expansions more generally.

Competing evolutionary paths in growing populations with applications to multidrug resistance

Investigating the emergence of a particular cell type is a recurring theme in models of growing cellular populations. The evolution of resistance to therapy is a classic example. Common questions are: when does the cell type first occur, and via which sequence of steps is it most likely to emerge? For growing populations, these questions can be formulated in a general framework of branching processes spreading through a graph from a root to a target vertex. Cells have a particular fitness value on each vertex and can transition along edges at specific rates. Vertices represent cell states, say genotypes or physical locations, while possible transitions are acquiring a mutation or cell migration. We focus on the setting where cells at the root vertex have the highest fitness and transition rates are small. Simple formulas are derived for the time to reach the target vertex and for the probability that it is reached along a given path in the graph. We demonstrate our results on several scenarios relevant to the emergence of drug resistance, including: the orderings of resistance-conferring mutations in bacteria and the impact of imperfect drug penetration in cancer.

How long does it take for a treatment naive, growing bacterial colony to be able to survive exposure to a cocktail of antibiotics? En route to multidrug resistance, what order did the drugs become impotent in? Questions such as these that pertain to the emergence of a significant cell type in a growing population arise frequently. They are often investigated via mathematical modelling but biologically insightful results are challenging to obtain. Here we outline a general framework of a stochastically growing population spreading through a graph to study such questions and provide simple formulas as answers. The significant cell type appears upon the population reaching a target vertex. Due to their simplicity, the derived formulas are widely accessible and can be used to guide and develop intuition on a range of biological scenarios. We demonstrate this on several settings including: how a region where drugs cannot penetrate affects the emergence of resistance, and, the ordering of mutations that leads to drugs being ineffective.

Environmental heterogeneity can tip the population genetics of range expansions

The population genetics of most range expansions is thought to be shaped by the competition between Darwinian selection and random genetic drift at the range margins. Here, we show that the evolutionary dynamics during range expansions is highly sensitive to additional fluctuations induced by environmental heterogeneities. Tracking mutant clones with a tunable fitness effect in bacterial colonies grown on randomly patterned surfaces we found that environmental heterogeneity can dramatically reduce the efficacy of selection. Time-lapse microscopy and computer simulations suggest that this effect arises generically from a local 'pinning' of the expansion front, whereby stretches of the front are slowed down on a length scale that depends on the structure of the environmental heterogeneity. This pinning focuses the range expansion into a small number of 'lucky' individuals with access to expansion paths, altering the neutral evolutionary dynamics and increasing the importance of chance relative to selection.

Bacterial growth: a statistical physicist's guide

Bacterial growth presents many beautiful phenomena that pose new theoretical challenges to statistical physicists, and are also amenable to laboratory experimentation. This review provides some of the essential biological background, discusses recent applications of statistical physics in this field, and highlights the potential for future research.

Access to high-impact mutations constrains the evolution of antibiotic resistance in soft agar

Despite widespread resistance to many important antibiotics, the factors that govern the emergence and prevalence of antibiotic-resistant bacteria are still unclear. When exposed to antibiotic gradients in soft agar plates measuring as little as 1.25 × 11 cm we found that Escherichia coli rapidly became resistant to representatives from every class of antibiotics active against Gram-negative bacteria. Evolution kinetics were independent of the frequency of spontaneous mutations that confer antibiotic resistance or antibiotic dose-response curves, and were only loosely correlated to maximal antibiotic concentrations. Instead, rapid evolution required unrealized mutations that could markedly decrease antibiotic susceptibility. When bacteria could not evolve through these “high-impact” mutations, populations frequently bottlenecked, reducing the number of cells from which mutants could arise and prolonging evolution times. This effect was independent of the antibiotic’s mechanism of action, and may affect the evolution of antibiotic resistance in clinical settings.

Quantifying the impact of a periodic presence of antimicrobial on resistance evolution in a homogeneous microbial population of fixed size

The evolution of antimicrobial resistance generally occurs in an environment where antimicrobial concentration is variable, which has dramatic consequences on the microorganisms' fitness landscape, and thus on the evolution of resistance. We investigate the effect of these time-varying patterns of selection within a stochastic model. We consider a homogeneous microbial population of fixed size subjected to periodic alternations of phases of absence and presence of an antimicrobial that stops growth. Combining analytical approaches and stochastic simulations, we quantify how the time necessary for fit resistant bacteria to take over the microbial population depends on the alternation period. We demonstrate that fast alternations strongly accelerate the evolution of resistance, reaching a plateau for sufficiently small periods. Furthermore, this acceleration is stronger in larger populations. For asymmetric alternations, featuring a different duration of the phases with and without antimicrobial, we shed light on the existence of a minimum for the time taken by the population to fully evolve resistance. The corresponding dramatic acceleration of the evolution of antimicrobial resistance likely occurs in realistic situations, and may have an important impact both in clinical and experimental situations.

Protein evolution speed depends on its stability and abundance and on chaperone concentrations

Some biological evolution is slow (millions of years), and some is fast (months to years). The speed at which a protein evolves depends on how stable a protein’s folded structure is, how well it avoids aggregation, and how well-chaperoned it is. What are the mechanisms? We compute fitness landscapes by combining a model of protein-folding equilibria with sequence-change dynamics. We find that adapting to a new environment is fastest for proteins that are least stably folded, because those sit on steep downhill parts of fitness potentials. The modeling shows that cells should adapt to warmer environments faster than to colder ones, explains why increasing a protein’s abundance slows cell evolution, and explains how chaperones accelerate evolution by mitigating this effect.

Proteins evolve at different rates. What drives the speed of protein sequence changes? Two main factors are a protein’s folding stability and aggregation propensity. By combining the hydrophobic–polar (HP) model with the Zwanzig–Szabo–Bagchi rate theory, we find that: (i) Adaptation is strongly accelerated by selection pressure, explaining the broad variation from days to thousands of years over which organisms adapt to new environments. (ii) The proteins that adapt fastest are those that are not very stably folded, because their fitness landscapes are steepest. And because heating destabilizes folded proteins, we predict that cells should adapt faster when put into warmer rather than cooler environments. (iii) Increasing protein abundance slows down evolution (the substitution rate of the sequence) because a typical protein is not perfectly fit, so increasing its number of copies reduces the cell’s fitness. (iv) However, chaperones can mitigate this abundance effect and accelerate evolution (also called evolutionary capacitance) by effectively enhancing protein stability. This model explains key observations about protein evolution rates.

Tuning Spatial Profiles of Selection Pressure to Modulate the Evolution of Drug Resistance

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 is sufficiently large relative to mutation but slows fixation for small migration rates. 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, even in the absence of a fitness cost for resistance. Our results may lay the groundwork for optimized, spatially resolved drug dosing strategies for mitigating the effects of drug resistance.

Predicting the dynamics of bacterial growth inhibition by ribosome-targeting antibiotics

Understanding how antibiotics inhibit bacteria can help to reduce antibiotic use and hence avoid antimicrobial resistance—yet few theoretical models exist for bacterial growth inhibition by a clinically relevant antibiotic treatment regimen. In particular, in the clinic, antibiotic treatment is time-dependent. Here, we use a theoretical model, previously applied to steady-state bacterial growth, to predict the dynamical response of a bacterial cell to a time-dependent dose of ribosome-targeting antibiotic. Our results depend strongly on whether the antibiotic shows reversible transport and/or low-affinity ribosome binding (‘low-affinity antibiotic’) or, in contrast, irreversible transport and/or high affinity ribosome binding (‘high-affinity antibiotic’). For low-affinity antibiotics, our model predicts that growth inhibition depends on the duration of the antibiotic pulse, and can show a transient period of very fast growth following removal of the antibiotic. For high-affinity antibiotics, growth inhibition depends on peak dosage rather than dose duration, and the model predicts a pronounced post-antibiotic effect, due to hysteresis, in which growth can be suppressed for long times after the antibiotic dose has ended. These predictions are experimentally testable and may be of clinical significance.

Exploiting ecology in drug pulse sequences in favour of population reduction

A deterministic population dynamics model involving birth and death for a two-species system, comprising a wild-type and more resistant species competing via logistic growth, is subjected to two distinct stress environments designed to mimic those that would typically be induced by temporal variation in the concentration of a drug (antibiotic or chemotherapeutic) as it permeates through the population and is progressively degraded. Different treatment regimes, involving single or periodical doses, are evaluated in terms of the minimal population size (a measure of the extinction probability), and the population composition (a measure of the selection pressure for resistance or tolerance during the treatment). We show that there exist timescales over which the low-stress regime is as effective as the high-stress regime, due to the competition between the two species. For multiple periodic treatments, competition can ensure that the minimal population size is attained during the first pulse when the high-stress regime is short, which implies that a single short pulse can be more effective than a more protracted regime. Our results suggest that when the duration of the high-stress environment is restricted, a treatment with one or multiple shorter pulses can produce better outcomes than a single long treatment. If ecological competition is to be exploited for treatments, it is crucial to determine these timescales, and estimate for the minimal population threshold that suffices for extinction. These parameters can be quantified by experiment.

Convection shapes the trade-off between antibiotic efficacy and the selection for resistance in spatial gradients

Since penicillin was discovered about 90 years ago, we have become used to using drugs to eradicate unwanted pathogenic cells. However, using drugs to kill bacteria, viruses or cancer cells has the serious side effect of selecting for mutant types that survive the drug attack. A crucial question therefore is how one could eradicate as many cells as possible for a given acceptable risk of drug resistance evolution. We address this general question in a model of drug resistance evolution in spatial drug gradients, which recent experiments and theories have suggested as key drivers of drug resistance. Importantly, our model takes into account the influence of convection, resulting for instance from blood flow. Using stochastic simulations, we study the fates of individual resistance mutations and quantify the trade-off between the killing of wild-type cells and the rise of resistance mutations: shallow gradients and convection into the antibiotic region promote wild-type death, at the cost of increasing the establishment probability of resistance mutations. We can explain these observed trends by modeling the adaptation process as a branching random walk. Our analysis reveals that the trade-off between death and adaptation depends on the relative length scales of the spatial drug gradient and random dispersal, and the strength of convection. Our results show that convection can have a momentous effect on the rate of establishment of new mutations, and may heavily impact the efficiency of antibiotic treatment.

Assessing the Risk of Probiotic Dietary Supplements in the Context of Antibiotic Resistance

Probiotic bacteria are known to harbor intrinsic and mobile genetic elements that confer resistance to a wide variety of antibiotics. Their high amounts in dietary supplements can establish a reservoir of antibiotic resistant genes in the human gut. These resistant genes can be transferred to pathogens that share the same intestinal habitat thus resulting in serious clinical ramifications. While antibiotic resistance of probiotic bacteria from food, human and animal sources have been well-documented, the resistant profiles of probiotics from dietary supplements have only been recently studied. These products are consumed with increasing regularity due to their health claims that include the improvement of intestinal health and immune response as well as prevention of acute and antibiotic-associated diarrhea and cancer; but, a comprehensive risk assessment on the spread of resistant genes to human health is lacking. Here, we highlight recent reports of antibiotic resistance of probiotic bacteria isolated from dietary supplements, and propose complementary strategies that can shed light on the risks of consuming such products in the context of a global widespread of antibiotic resistance. In concomitant with a broader screening of antibiotic resistance in probiotic supplements is the use of computational simulations, live imaging and functional genomics to harvest knowledge on the evolutionary behavior, adaptations and dynamics of probiotics studied in conditions that best represent the human gut including in the presence of antibiotics. The underlying goal is to enable the health benefits of probiotics to be exploited in a responsible manner and with minimal risk to human health.

Drug delivery in a tumour cord model: a computational simulation

The tumour vasculature and microenvironment is complex and heterogeneous, contributing to reduced delivery of cancer drugs to the tumour. We have developed an in silico model of drug transport in a tumour cord to explore the effect of different drug regimes over a 72 h period and how changes in pharmacokinetic parameters affect tumour exposure to the cytotoxic drug doxorubicin. We used the model to describe the radial and axial distribution of drug in the tumour cord as a function of changes in the transport rate across the cell membrane, blood vessel and intercellular permeability, flow rate, and the binding and unbinding ratio of drug within the cancer cells. We explored how changes in these parameters may affect cellular exposure to drug. The model demonstrates the extent to which distance from the supplying vessel influences drug levels and the effect of dosing schedule in relation to saturation of drug-binding sites. It also shows the likely impact on drug distribution of the aberrant vasculature seen within tumours. The model can be adapted for other drugs and extended to include other parameters. The analysis confirms that computational models can play a role in understanding novel cancer therapies to optimize drug administration and delivery.

Spatiotemporal microbial evolution on antibiotic landscapes

A key aspect of bacterial survival is the ability to evolve while migrating across spatially varying environmental challenges. Laboratory experiments, however, often study evolution in well-mixed systems. Here, we introduce an experimental device, the microbial evolution and growth arena (MEGA)-plate, in which bacteria spread and evolved on a large antibiotic landscape (120 × 60 centimeters) that allowed visual observation of mutation and selection in a migrating bacterial front. While resistance increased consistently, multiple coexisting lineages diversified both phenotypically and genotypically. Analyzing mutants at and behind the propagating front, we found that evolution is not always led by the most resistant mutants; highly resistant mutants may be trapped behind more sensitive lineages. The MEGA-plate provides a versatile platform for studying microbial adaption and directly visualizing evolutionary dynamics.

Toward a quantitative understanding of antibiotic resistance evolution

The rising prevalence of antibiotic resistant bacteria is an increasingly serious public health challenge. To address this problem, recent work ranging from clinical studies to theoretical modeling has provided valuable insights into the mechanisms of resistance, its emergence and spread, and ways to counteract it. A deeper understanding of the underlying dynamics of resistance evolution will require a combination of experimental and theoretical expertise from different disciplines and new technology for studying evolution in the laboratory. Here, we review recent advances in the quantitative understanding of the mechanisms and evolution of antibiotic resistance. We focus on key theoretical concepts and new technology that enables well-controlled experiments. We further highlight key challenges that can be met in the near future to ultimately develop effective strategies for combating resistance.

Beyond the Hypercube: Evolutionary Accessibility of Fitness Landscapes with Realistic Mutational Networks

Evolutionary pathways describe trajectories of biological evolution in the space of different variants of organisms (genotypes). The probability of existence and the number of evolutionary pathways that lead from a given genotype to a better-adapted genotype are important measures of accessibility of local fitness optima and the reproducibility of evolution. Both quantities have been studied in simple mathematical models where genotypes are represented as binary sequences of two types of basic units, and the network of permitted mutations between the genotypes is a hypercube graph. However, it is unclear how these results translate to the biologically relevant case in which genotypes are represented by sequences of more than two units, for example four nucleotides (DNA) or 20 amino acids (proteins), and the mutational graph is not the hypercube. Here we investigate accessibility of the best-adapted genotype in the general case of K > 2 units. Using computer generated and experimental fitness landscapes we show that accessibility of the global fitness maximum increases with K and can be much higher than for binary sequences. The increase in accessibility comes from the increase in the number of indirect trajectories exploited by evolution for higher K. As one of the consequences, the fraction of genotypes that are accessible increases by three orders of magnitude when the number of units K increases from 2 to 16 for landscapes of size N ∼ 10⁶ genotypes. This suggests that evolution can follow many different trajectories on such landscapes and the reconstruction of evolutionary pathways from experimental data might be an extremely difficult task.

Biological evolution is driven by heritable, genetic alterations that affect the fitness of organisms. However, the pool of “fitter” variants (genotypes) is often restricted and it is not at all obvious how evolution finds its way from low-fitness to high-fitness genotypes in a complex, multidimensional “fitness landscapes” with many peaks (fit organisms) and valleys (unfit ones). To address this question we investigate how likely it is for biological evolution to find a way “uphill” from a lower-fitness organism to the best adapted organism. We discover that the accessibility of the fittest organism depends on the number of types of basic “units” used to encode genotypes. These units can be, for example, the four DNA nucleotides A,T,C,G, or the ∼20 amino acids used for synthesizing proteins, and the choice of the most appropriate unit is dictated by how the genotypes and the fitnesses are related—a relationship that researchers have begun to unveil only recently. We find that increasing the number of units strongly increases the probability that there will be at least one uphill path to the best-adapted genotype, and the number of evolutionary pathways leading to it. Our findings suggest that biological evolution can follow many more pathways than previously thought.

Evolutionary scalpels for dissecting tumor ecosystems

Amidst the growing literature on cancer genomics and intratumor heterogeneity, essential principles in evolutionary biology recur time and time again. Here we use these principles to guide the reader through major advances in cancer research, highlighting issues of "hit hard, hit early" treatment strategies, drug resistance, and metastasis. We distinguish between two frameworks for understanding heterogeneous tumors, both of which can inform treatment strategies: (1) The tumor as diverse ecosystem, a Darwinian population of sometimes-competing, sometimes-cooperating cells; (2) The tumor as tightly integrated, self-regulating organ, which may hijack developmental signals to restore functional heterogeneity after treatment. While the first framework dominates literature on cancer evolution, the second framework enjoys support as well. Throughout this review, we illustrate how mathematical models inform understanding of tumor progression and treatment outcomes. Connecting models to genomic data faces computational and technical hurdles, but high-throughput single-cell technologies show promise to clear these hurdles. This article is part of a Special Issue entitled: Evolutionary principles - heterogeneity in cancer?, edited by Dr. Robert A. Gatenby.

The adaptation rate of a quantitative trait in an environmental gradient

The spatial range of a species habitat is generally determined by the ability of the species to cope with biotic and abiotic variables that vary in space. Therefore, the species range is itself an evolvable property. Indeed, environmental gradients permit a mode of evolution in which range expansion and adaptation go hand in hand. This process can contribute to rapid evolution of drug resistant bacteria and viruses, because drug concentrations in humans and livestock treated with antibiotics are far from uniform. Here, we use a minimal stochastic model of discrete, interacting organisms evolving in continuous space to study how the rate of adaptation of a quantitative trait depends on the steepness of the gradient and various population parameters. We discuss analytical results for the mean-field limit as well as extensive stochastic simulations. These simulations were performed using an exact, event-driven simulation scheme that can deal with continuous time-, density- and coordinate-dependent reaction rates and could be used for a wide variety of stochastic systems. The results reveal two qualitative regimes. If the gradient is shallow, the rate of adaptation is limited by dispersion and increases linearly with the gradient slope. If the gradient is steep, the adaptation rate is limited by mutation. In this regime, the mean-field result is highly misleading: it predicts that the adaptation rate continues to increase with the gradient slope, whereas stochastic simulations show that it in fact decreases with the square root of the slope. This discrepancy underscores the importance of discreteness and stochasticity even at high population densities; mean-field results, including those routinely used in quantitative genetics, should be interpreted with care.

Evolution of Drug Resistance in Bacteria

Resistance to antibiotics is an important and timely problem of contemporary medicine. Rapid evolution of resistant bacteria calls for new preventive measures to slow down this process, and a longer-term progress cannot be achieved without a good understanding of the mechanisms through which drug resistance is acquired and spreads in microbial populations. Here, we discuss recent experimental and theoretical advances in our knowledge how the dynamics of microbial populations affects the evolution of antibiotic resistance . We focus on the role of spatial and temporal drug gradients and show that in certain situations bacteria can evolve de novo resistance within hours. We identify factors that lead to such rapid onset of resistance and discuss their relevance for bacterial infections.

Antibiotic resistance: a physicist's view

The problem of antibiotic resistance poses challenges across many disciplines. One such challenge is to understand the fundamental science of how antibiotics work, and how resistance to them can emerge. This is an area where physicists can make important contributions. Here, we highlight cases where this is already happening, and suggest directions for further physics involvement in antimicrobial research.

Allele surfing promotes microbial adaptation from standing variation

The coupling of ecology and evolution during range expansions enables mutations to establish at expanding range margins and reach high frequencies. This phenomenon, called allele surfing, is thought to have caused revolutions in the gene pool of many species, most evidently in microbial communities. It has remained unclear, however, under which conditions allele surfing promotes or hinders adaptation. Here, using microbial experiments and simulations, we show that, starting with standing adaptive variation, range expansions generate a larger increase in mean fitness than spatially uniform population expansions. The adaptation gain results from 'soft' selective sweeps emerging from surfing beneficial mutations. The rate of these surfing events is shown to sensitively depend on the strength of genetic drift, which varies among strains and environmental conditions. More generally, allele surfing promotes the rate of adaptation per biomass produced, which could help developing biofilms and other resource-limited populations to cope with environmental challenges.

Competition for space during bacterial colonization of a surface

Competition for space is ubiquitous in the ecology of both microorganisms and macro-organisms. We introduce a bacterial model system in which the factors influencing competition for space during colonization of an initially empty habitat can be tracked directly. Using fluorescence microscopy, we follow the fate of individual Escherichia coli bacterial cell lineages as they undergo expansion competition (the race to be the first to colonize a previously empty territory), and as they later compete at boundaries between clonal territories. Our experiments are complemented by computer simulations of a lattice-based model. We find that both expansion competition, manifested as differences in individual cell lag times, and boundary competition, manifested as effects of neighbour cell geometry, can play a role in colonization success, particularly when lineages expand exponentially. This work provides a baseline for investigating how ecological interactions affect colonization of space by bacterial populations, and highlights the potential of bacterial model systems for the testing and development of ecological theory.

Growth-dependent bacterial susceptibility to ribosome-targeting antibiotics

Bacterial growth environment strongly influences the efficacy of antibiotic treatment, with slow growth often being associated with decreased susceptibility. Yet in many cases, the connection between antibiotic susceptibility and pathogen physiology remains unclear. We show that for ribosome-targeting antibiotics acting on Escherichia coli, a complex interplay exists between physiology and antibiotic action; for some antibiotics within this class, faster growth indeed increases susceptibility, but for other antibiotics, the opposite is true. Remarkably, these observations can be explained by a simple mathematical model that combines drug transport and binding with physiological constraints. Our model reveals that growth-dependent susceptibility is controlled by a single parameter characterizing the ‘reversibility’ of ribosome-targeting antibiotic transport and binding. This parameter provides a spectrum classification of antibiotic growth-dependent efficacy that appears to correspond at its extremes to existing binary classification schemes. In these limits, the model predicts universal, parameter-free limiting forms for growth inhibition curves. The model also leads to non-trivial predictions for the drug susceptibility of a translation mutant strain of E. coli, which we verify experimentally. Drug action and bacterial metabolism are mechanistically complex; nevertheless, this study illustrates how coarse-grained models can be used to integrate pathogen physiology into drug design and treatment strategies.

The spatial profiles and metabolic capabilities of microbial populations impact the growth of antibiotic-resistant mutants

Antibiotic resistance adversely affects clinical and public health on a global scale. Using the opportunistic human pathogen Pseudomonas aeruginosa, we show that increasing the number density of bacteria, on agar containing aminoglycoside antibiotics, can non-monotonically impact the survival of antibiotic-resistant mutants. Notably, at high cell densities, mutant survival is inhibited. A wide range of bacterial species can inhibit antibiotic-resistant mutants. Inhibition results from the metabolic breakdown of amino acids, which results in alkaline by-products. The consequent increase in pH acts in conjunction with aminoglycosides to mediate inhibition. Our work raises the possibility that the manipulation of microbial population structure and nutrient environment in conjunction with existing antibiotics could provide therapeutic approaches to combat antibiotic resistance.

Quantifying the Determinants of Evolutionary Dynamics Leading to Drug Resistance

The emergence of drug resistant pathogens is a serious public health problem. It is a long-standing goal to predict rates of resistance evolution and design optimal treatment strategies accordingly. To this end, it is crucial to reveal the underlying causes of drug-specific differences in the evolutionary dynamics leading to resistance. However, it remains largely unknown why the rates of resistance evolution via spontaneous mutations and the diversity of mutational paths vary substantially between drugs. Here we comprehensively quantify the distribution of fitness effects (DFE) of mutations, a key determinant of evolutionary dynamics, in the presence of eight antibiotics representing the main modes of action. Using precise high-throughput fitness measurements for genome-wide Escherichia coli gene deletion strains, we find that the width of the DFE varies dramatically between antibiotics and, contrary to conventional wisdom, for some drugs the DFE width is lower than in the absence of stress. We show that this previously underappreciated divergence in DFE width among antibiotics is largely caused by their distinct drug-specific dose-response characteristics. Unlike the DFE, the magnitude of the changes in tolerated drug concentration resulting from genome-wide mutations is similar for most drugs but exceptionally small for the antibiotic nitrofurantoin, i.e., mutations generally have considerably smaller resistance effects for nitrofurantoin than for other drugs. A population genetics model predicts that resistance evolution for drugs with this property is severely limited and confined to reproducible mutational paths. We tested this prediction in laboratory evolution experiments using the "morbidostat", a device for evolving bacteria in well-controlled drug environments. Nitrofurantoin resistance indeed evolved extremely slowly via reproducible mutations-an almost paradoxical behavior since this drug causes DNA damage and increases the mutation rate. Overall, we identified novel quantitative characteristics of the evolutionary landscape that provide the conceptual foundation for predicting the dynamics of drug resistance evolution.