Spatiotemporal microbial evolution on antibiotic landscapes

See allHide authors and affiliations

Science  09 Sep 2016:
Vol. 353, Issue 6304, pp. 1147-1151
DOI: 10.1126/science.aag0822


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.

The worldwide increase in antibiotic resistance has motivated numerous studies aimed at understanding the phenotypic and genotypic evolution of antibiotic resistance (17). These experiments have shed light on the trade-offs constraining adaptive evolution in single- and multidrug environments (5, 6, 8, 9). However, most of our current knowledge about the evolution of resistance is based on laboratory setups with well-mixed environments (17, 10, 11).

In natural and clinical settings, bacteria migrate between spatially distinct regions of selection (5, 6, 8, 9, 12, 13). Theoretical models show that spatially structured pressures change the nature of selection: Instead of competing with its neighbors for limited resources, an adapted individual needs only to be the first with the capability to venture and survive in a new region (14, 15). A pioneering study focusing on small population sizes showed that structured microenvironments increase the rate of adaptation to antibiotics through highly reproducible genetic changes (9). It is unknown how evolution is shaped by the diversification potential and differences in adaptive constraints of large populations in spatial environments.

Here, we present a device for the evolution of bacteria that allows migration and adaptation in a large, spatially structured environment. The microbial evolution and growth arena (MEGA)–plate consists of a rectangular acrylic dish, 120 × 60 cm, in which successive regions of black-colored agar containing different concentrations of antibiotics are overlaid by soft agar allowing bacterial motility (Fig. 1A). Motile bacteria inoculated at one location on the plate deplete nutrients locally and then spread by chemotaxis to other regions (16). Only increasingly resistant mutants can spread into sections containing higher levels of antibiotic. The large size of the plate serves two purposes: It provides for a large population and mutational supply, and it maintains the antibiotic gradient despite diffusion (drug diffusion time scales quadratically with distance while the bacterial front advances linearly; thus, the large plate size prevents the antibiotic gradient from equilibrating over the duration of the experiment). Once a mutant has exhausted the resources of a region of the plate, other mutants do not meaningfully migrate by chemotaxis to that region (because they move diffusively without a nutrient gradient). In this manner, mutational lineages can block each other physically—a phenomenon notably observed in biofilm formation (17). This partitioning of mutants into stable spatial domains also enables sampling of individual mutants for later analysis. Using periodic photography of the plate, we constructed time-lapse movies of evolution (movie S1). Combining these with analysis of isolates, this system allows reconstruction of the phenotypic and genotypic evolutionary histories of evolving bacteria.

Fig. 1 An experimental device for studying microbial evolution in a spatially structured environment.

(A) Setup of the four-step gradient of trimethoprim (TMP). Antibiotic is added in sections to make an exponential gradient rising inward. (B) The four-step TMP MEGA-plate after 12 days. E. coli appear as white on the black background. The 182 sampled points of clones are indicated by circles, colored by their measured MIC. Lines indicate video-imputed ancestry. (C) Time-lapse images of a section of the MEGA-plate. Repeated mutation and selection can be seen at each step. Images have been aligned and linearly contrast-enhanced but are otherwise unedited.

Challenging bacteria in spatial gradients of antibiotics leads to large increases in resistance through sequential adaptive steps across competing lineages (Fig. 1 and movie S1). We first set up the MEGA-plate with symmetric four-step gradients of trimethoprim (TMP) or ciprofloxacin (CPR) proceeding inward with order-of-magnitude increases in concentration per step [Fig. 1A; TMP: 0, 3, 30, 300, and 3000 × wild-type minimum inhibitory concentration (MIC); CPR: 0, 20, 200, 2000, and 20,000 × MIC] and inoculated the drug-free regions with Escherichia coli. Bacteria swim and spread until they reach a concentration in which they can no longer grow (TMP, Fig. 1C and movies S1 and S2; CPR, movie S3). As resistant mutants arise in the population, their descendants migrate into the next step of drug concentration and fan out (Fig. 1C, 88 hours). Adjacent mutant lineages exclude each other and compete for limited space, resulting in some lineages entirely blocking off growth of others (Fig. 1C). When the winning lineages reach a further increased level of drug concentration at which they too are unable to grow, secondary mutations arise and the process repeats. Ultimately, the bacteria reach and overspread the highest drug concentration, showing marked increases in drug resistance: Phenotyping of sampled mutants from the highest-concentration region showed a factor of 104 increase in MIC for TMP (Fig. 1B) and a factor of 105 increase in MIC for CPR (fig. S1). The adaption time (10 days in TMP, 12 days in CPR) is consistent with evolution in well-mixed environments (4), yet is slower than reported adaptation rates in microspatial environments, likely because of the additional time required to swim between concentration steps (9). It is possible that at different dimensions, the MEGA-plate will yield different evolutionary dynamics; a wider front would increase the effective population size and thus the mutational supply, whereas a longer run between steps would increase selection among adjacent lineages.

To test the importance of the size of intermediate steps in the evolution of high-level resistance, we set up a variant of the MEGA-plate in which bacteria go from no drug to a high level directly or through one middle region of variable magnitude (Fig. 2; TMP: high step 3000 × MIC, middle step 0, 3, 30, or 300 × MIC; CPR: high step 2000 × MIC, middle step 0, 2, 20, or 200 × MIC). Bacteria were unable to adapt directly from zero to the highest concentration of either drug. Diffusive smoothing of these large steps enabled the appearance of partially resistant mutants, but their lineages did not advance (Fig. 2A, left). The addition of an intermediate concentration step enabled adaptation, although this was impeded when this middle step was too high (Fig. 2B). Even across the permissive intermediate steps, evolution often proceeded through multiple mutations taking advantage of the local gradients formed by diffusion (TMP, movie S4; CPR, movie S5). Thus, by progressing through colonization of regions with moderately challenging selective pressures, intermediate-resistance mutants can expand to sufficient numbers to facilitate the rise of high-resistance mutants. Analogous to evolutionary rescue in temporal selective gradients (1820), a gradual spatial gradient allows adaptation to previously inhospitable environments. However, unlike in a temporal gradient, a spatial gradient does not impose a minimal time for a mutant’s appearance and spread; at any time, a mutant appearing on the stalled front can expand and evolve further, provided it is sufficiently resistant to colonize the next step. Thus, concordant with theoretical predictions (21, 22), access to intermediate regions of moderate selection is critical for enabling a range of evolutionary paths to high-level resistance.

Fig. 2 Initial adaptation to low drug concentrations facilitates later adaptation to high concentrations.

(A) Frames from a section of the TMP intermediate-step MEGA-plate over time (TMP, movie S4; CPR, movie S5). The first frame showing a mutant in the highest band is indicated by a blue box. (B) Rates of adaptation in the intermediate-step experiments across TMP and CPR, showing the necessity of intermediate adaptation for the evolution of high levels of resistance. Error bars show the appearance times of multiple lineages in the highest concentration. Because the intermediate step with no drug puts the highest and lowest concentrations adjacent, it serves as both the highest and lowest intermediate steps (dashed line).

We next focused on the genotypic and phenotypic paths leading to high levels of resistance. We sequenced 21 isolates from the four-step TMP gradient experiment and 230 isolates from the multiple intermediate-step TMP experiment above. The samples separated into minimally and highly mutated (i.e., mutator phenotype) groups [>60 single-nucleotide polymorphisms (SNPs) and indels for high, <12 for low; Fig. 3A]. Similar separation was also seen when restricting the analysis to synonymous mutations to minimize differential effects of selection (fig. S2). All highly mutated isolates, but none of the others, had mutations in dnaQ (also called mutD), which was the only gene for which the presence of a mutation correlated perfectly with the mutator phenotype. This gene encodes DNA polymerase III, which is critical to proofreading (23, 24). Isolates carrying mutated dnaQ alleles showed increased rates of mutations on rifampin disk diffusion assays (fig. S2) (25). These mutators appeared repeatedly in distinct locations on the plate and across experiments. On the basis of lineage reconstructions from the time-lapse video as well as genotypic relationships, the mutator phenotype emerged at least six times independently between the four-step and intermediate-step TMP experiments above [fig. S2; four different alleles of dnaQ were observed: Val96 → Glu (V96E), Ile97 → Asn (I97N), Ile97 → Ser (I97S), and Ile97 → Thr (I97T), where I97T appeared three independent times]. Although these mutator lineages accumulated mutations more rapidly, their rate of phenotypic adaption was similar to that of the less mutated isolates, reaching the highest level of resistance at roughly the same time (fig. S3 and movie S4). Indeed, the highly mutated lineages had a close to neutral ratio of nonsynonymous to synonymous substitutions (Fig. 3B). In contrast, the less mutated isolates showed a high bias toward coding mutations, indicating that most of these mutations were likely adaptive (Fig. 3B).

Fig. 3 Diverse genotypic strategies for adaptation to trimethoprim.

(A) Numbers of observed mutations across individual isolates. Samples with a dnaQ mutation (solid symbols) consistently carried more mutations than those sampled with the wild-type dnaQ allele (crosses). Data points are horizontally jittered for clarity. (B) The normalized ratio of nonsynonymous to synonymous substitutions of isolates compared with the ancestor for samples with normal and highly mutated phenotypes. Error bars are the standard deviation of the Bayesian posterior estimate for the binomial parameter. (C) Numbers of distinct mutational events in genes that were mutated at least twice independently. Genes are colored by pathway per EcoCyc (37). Nonsynonymous, synonymous, and loss-of-function mutations (including indel and nonsense) are indicated. Genes that only had one mutation across all samples were combined into the “unique” column. Individual mutation events were inferred through ancestry (movies S1 and S4). Inset: The multistep MEGA-plate with samples containing the mutation soxR R20C (yellow) tracing mutational events from multiple samples and video.

Focusing on the nonmutator isolates, we identified a wide spectrum of putatively adaptive mutations for TMP resistance. The most frequently mutated gene was the primary target of TMP, folA (26), which encodes dihydrofolate reductase (DHFR) (Fig. 3C), with more mutations appearing as resistance increased. We also observed several genes that were repeatedly mutated yet are not involved in the folate biosynthesis pathway, and thus are not primarily associated with TMP resistance. These included stress response genes, such as those of the mar and sox operons, known to be important in general antibiotic and toxin resistance (27), as well as genes involved in transcription and translation, which have been shown to affect TMP resistance (28). We also found that three genes not classically associated with TMP resistance were repeatedly mutated, often with a probable loss of function (frameshift or nonsense): a phosphate transporter (pitA), shikimate kinase I (aroK), and a negative regulator of the PhoQP system (mgrB). Knockout of these genes in the ancestral strain confirmed their resistance phenotypes (fig. S4).

Mutations that increased resistance often came with a cost of reduced growth, which was subsequently restored by additional compensatory mutations (2931). Although some resistance-conferring mutations allowed colonization of regions of high drug concentration without affecting growth, many lineages capable of growing in these regions were deficient in yield, particularly during CPR resistance evolution (as measured by optical density: Fig. 4, A and B, for CPR; fig. S3 for TMP). These yield-deficient mutations were followed by compensatory mutations allowing growth to full density (Fig. 4, A and B, and movie S3; number of compensatory mutants observed in a single run: >50 for TMP, >500 for CPR). In the absence of a chemotaxis-inducing nutrient gradient, the compensatory mutants stayed localized behind the front, appearing in a characteristic pattern of localized spots spreading from single points (Fig. 4A and movie S3).

Fig. 4 Compensatory mutations can be spatially trapped.

(A) Ciprofloxacin experiment still frame with locations of 24 isolates showing a full-fitness mutation (cyan) or yield-deficient mutation followed by a compensatory mutation (purple). (B) Optical density at the marked points in (A) over the course of the experiment. The two example traces (indicated by a square and triangle) correspond to the points marked by the same glyph in (A). (C) Mutants isolated behind the front can have markedly higher resistance than the front at the time it passed the same location. Resistance of the front was measured by the concentration at which front progression stopped; isolate MICs were measured in vitro. (D) The front (marked F) and three compensatory mutants (marked 1 to 3) were sampled at 162 hours, and immediately inoculated ahead of the front as indicated by the arrows. Growth of the moved mutants is evident for the three compensatory mutants, despite being inoculated at a CPR concentration much higher than where they emerged, but not for the front. (E) Measured CPR MICs of the mutants from (D).

Focusing on evolution in CPR, we sampled and phenotyped compensatory mutants. We found that many of them had not only compensated for growth but had also increased in resistance, often beyond the resistance levels of the propagating front (Fig. 4C). Yet, as these mutants were engulfed by their parental lineage, they stayed constrained to the immediate vicinity in which they appeared and were unable to overtake the moving front. To test whether these compensatory mutants were capable of outcompeting the propagating front, we conducted an additional evolution experiment in which we sampled the trapped compensatory mutants and moved them forward, reinoculating them ahead of the still-moving front. These compensatory mutants were able to grow in a region where the front could not (Fig. 4D). Similarly, some trapped compensatory mutants were able to outcompete their parent when placed side-by-side on a fresh gradient plate (fig. S5). Hence, as compensatory mutations often occur behind the front, they are spatially restricted from contributing to the ultimate evolutionary course of the population. Indeed, in the rare cases when these compensatory mutations appeared at the front and were not physically blocked, they accelerated the adaptive process (fig. S3 and movie S3, 00:53). Thus, the fitness of the population is not driven by the fittest mutants (3234), but rather by those that are both sufficiently fit and arise sufficiently close to the advancing front.

The MEGA-plate is not intended to directly simulate natural or clinical settings, but it does capture unique aspects of evolution during range expansion. Evolution of high levels of resistance is enabled by intermediate regions of moderate selective pressure. Furthermore, as multiple lineages evolve in parallel, the propagating front can be led by lineages less fit than those trapped behind it. It will be interesting to explore how adaptation rates and mutational diversity depend on other spatiotemporal parameters, including population density, mutation rate, and the relative expansion speed and spatial dimensions. Owing to the relaxed evolutionary constraints in range expansion dynamics, the MEGA-plate is likely to reveal novel mutational pathways to high-level multiantibiotic resistance. Further, the MEGA-plate can be adapted to a range of organisms and challenges beyond antibiotics. Differences in evolutionary dynamics between evolution under different selection pressures appear visually, simplifying both hypothesis generation and testing. Owing to this flexibility, the MEGA-plate is a platform for exploring the interplay of spatial constraints and evolutionary pressures. The MEGA-plate provides a physical analog of the otherwise abstract Muller plots of population genetics (35, 36) and of other elusive aspects of evolution, including diversification, compensatory mutations, and clonal interference. Its relative simplicity and ability to visually demonstrate evolution makes the MEGA-plate a useful tool for science education and outreach.

Supplementary Materials

Materials and Methods

Figs. S1 to S6

Tables S1 and S2

Movies S1 to S5

References (3843)

References and Notes

Acknowledgments: Sequence data are available on NCBI SRA under accession number SRP077287. We thank X. R. Bao and A. C. Palmer for helpful discussions and the National BioResource Project (NIG, Japan) for providing the Keio collection. Supported by National Defense Science and Engineering Graduate fellowship 32 CFR 168a (E.D.K.), NIH grant R01-GM081617 (R.K.), and European Research Council FP7 ERC grant 281891 (R.K.).
View Abstract

Stay Connected to Science

Navigate This Article