Range expansions transition from pulled to pushed waves as growth becomes more cooperative in an experimental microbial population

Biotic resistance predictably shifts microbial invasion regimes

Invading new territory is a central aspect of the microbial lifestyle. However, invading microbes rarely find novel territories uninhabited; resident microbes can interact with the newcomers and, in many cases, impede their invasion - an effect known as 'biotic resistance'. Accordingly, invasions are shaped by the interplay between dispersal and resistance. However, these two factors are difficult to disentangle or manipulate in natural systems, making their interplay challenging to understand. To address this challenge, we track microbial invasions in the lab over space and time - first in a model system of two interacting microbes, then in a multi-strain system involving a pathogen invading resident communities. In the presence of biotic resistance, we observe three qualitatively different invasion regimes: 'consistent', 'pulsed', and 'pinned', where, in the third regime, strong biotic resistance stalls the invasion entirely despite ongoing invader dispersal. These rich invasion dynamics could be qualitatively predicted with a simple, parameter-free framework that ignores individual species interactions, even for rather complex communities. Moreover, we show that this simple framework could accurately predict simulated invasions from different mechanistic models, indicating its broad applicability. Our work offers an understanding of how biotic resistance impacts invasions and introduces a predictive tool to identify invasion-resistant communities.

Linking spatial drug heterogeneity to microbial growth dynamics in theory and experiment

Diffusion and migration play pivotal roles in microbial communities - shaping, for example, colonization in new environments and the maintenance of spatial structures of biodiversity. While previous research has extensively studied free diffusion, such as range expansion, there remains a gap in understanding the effects of biologically or physically deleterious confined environments. In this study, we examine the interplay between migration and spatial drug heterogeneity within an experimental meta-community of E. faecalis, a Gram-positive opportunistic pathogen. When the community is confined to spatially-extended habitats ('islands') bordered by deleterious conditions, we find that the population level response depends on the trade-off between the growth rate within the island and the rate of transfer into regions with harsher conditions, a phenomenon we explore by modulating antibiotic concentration within the island. In heterogeneous islands, composed of spatially patterned patches that support varying levels of growth, the population's fate depends critically on the specific spatial arrangement of these patches - the same spatially averaged growth rate leads to diverging responses. These results are qualitatively captured by simple simulations, and analytical expressions which we derive using first-order perturbation approximations to reaction-diffusion models with explicit spatial dependence. Among all possible spatial arrangements, our theoretical and experimental findings reveal that the arrangement with the highest growth rates at the center most effectively mitigates population decline, while the arrangement with the lowest growth rates at the center is the least effective. Extending this approach to more complex experimental communities with varied spatial structures, such as a ring-structured community, further validates the impact of spatial drug arrangement. Our findings suggest new approaches to interpreting diverging clinical outcomes when applying identical drug doses and inform the possible optimization of spatially-explicit dosing strategies.

Recolonization of secondary forests by a locally extinct Caribbean anole through the lens of range expansion theory

Disturbance and recovery dynamics are characteristic features of many ecosystems. Disturbance dynamics are widely studied in ecology and conservation biology. Still, we know less about the ecological processes that drive ecosystem recovery. The ecological processes that mediate ecosystem recovery stand at the intersection of many theoretical frameworks. Range expansion theory is one of these complementary frameworks that can provide unique insights into the population-level processes that mediate ecosystem recovery, particularly fauna recolonization. Although the biodiversity patterns that follow the fauna recolonization of recovering forests have been well described in the literature, the ecological processes at the population level that drive these patterns remain conspicuously unknown. In this study, we tested three fundamental predictions of range expansion theory during the recolonization of recovering forests in Puerto Rico by a shade specialist anole, Anolis gundlachi. Range expansion theory predicts that individuals at the early stages of recolonization (i.e., younger forests) would have a high prevalence of dispersive traits, experience less density dependence, and suffer less parasitism. To test these predictions, we conducted a chronosequence study applying space-for-time substitution where we compared phenotypic traits (i.e., body size, body condition, and relative limb size), population density, population growth rates, and Plasmodium parasitism rates among lizard populations living in young (<30 years), mid (~40-70 years), and old-growth forests (>75 years). Lizard populations in younger forests had lower densities, higher population growth rates, and lower rates of Plasmodium parasitism compared with old-growth forests. Still, while we found that individuals had larger body sizes, and longer forelimbs in young forests in one site, this result was not consistent among sites. This suggests a potential trade-off between the traits that provide a dispersal advantage during the initial stages of recolonization and those that are advantageous to establish in novel environmental conditions. Overall, our study emphasizes the suitability of range expansion theory to describe fauna recolonization but also highlights that the ecological processes that drive recolonization are time-dependent, complex, and nuanced.

Tipping points emerge from weak mutualism in metacommunities

The coexistence of obligate mutualists is often precariously close to tipping points where small environmental changes can drive catastrophic shifts in species composition. For example, microbial ecosystems can collapse by the decline of a strain that provides an essential resource on which other strains cross-feed. Here, we show that tipping points, ecosystem collapse, bistability and hysteresis arise even with very weak (non-obligate) mutualism provided the population is spatially structured. Based on numeric solutions of a metacommunity model and mean-field analyses, we demonstrate that weak mutualism lowers the minimal dispersal rate necessary to avoid stochastic extinction, while species need to overcome a mean threshold density to survive in this low dispersal rate regime. Our results allow us to make numerous predictions for mutualistic metacommunities regarding tipping points, hysteresis effects, and recovery from external perturbations, and let us draw general conclusions for ecosystems even with random, not necessarily mutualistic, interactions and systems with density-dependent dispersal rather than direct mutualistic interactions.

Mutualism at the leading edge: insights into the eco-evolutionary dynamics of host-symbiont communities during range expansion

The evolution of mutualism between host and symbiont communities plays an essential role in maintaining ecosystem function and should therefore have a profound effect on their range expansion dynamics. In particular, the presence of mutualistic symbionts at the leading edge of a host-symbiont community should enhance its propagation in space. We develop a theoretical framework that captures the eco-evolutionary dynamics of host-symbiont communities, to investigate how the evolution of resource exchange may shape community structure during range expansion. We consider a community with symbionts that are mutualistic or parasitic to various degrees, where parasitic symbionts receive the same amount of resource from the host as mutualistic symbionts, but at a lower cost. The selective advantage of parasitic symbionts over mutualistic ones is increased with resource availability (i.e. with host density), promoting mutualism at the range edges, where host density is low, and parasitism at the population core, where host density is higher. This spatial selection also influences the speed of spread. We find that the host growth rate (which depends on the average benefit provided by the symbionts) is maximal at the range edges, where symbionts are more mutualistic, and that host-symbiont communities with high symbiont density at their core (e.g. resulting from more mutualistic hosts) spread faster into new territories. These results indicate that the expansion of host-symbiont communities is pulled by the hosts but pushed by the symbionts, in a unique push-pull dynamic where both the host and symbionts are active and tightly-linked players.

Imaging with spatio-temporal modelling to characterize the dynamics of plant-pathogen lesions

Within-host spread of pathogens is an important process for the study of plant-pathogen interactions. However, the development of plant-pathogen lesions remains practically difficult to characterize beyond the common traits such as lesion area. Here, we address this question by combining image-based phenotyping with mathematical modelling. We consider the spread of Peyronellaea pinodes on pea stipules that were monitored daily with visible imaging. We assume that pathogen propagation on host-tissues can be described by the Fisher-KPP model where lesion spread depends on both a logistic growth and an homogeneous diffusion. Model parameters are estimated using a variational data assimilation approach on sets of registered images. This modelling framework is used to compare the spread of an aggressive isolate on two pea cultivars with contrasted levels of partial resistance. We show that the expected slower spread on the most resistant cultivar is actually due to a significantly lower diffusion coefficient. This study shows that combining imaging with spatial mechanistic models can offer a mean to disentangle some processes involved in host-pathogen interactions and further development may allow a better identification of quantitative traits thereafter used in genetics and ecological studies.

Genetic diversity and structure of a recent fish invasion: Tench (Tinca tinca) in eastern North America

Introduced and geographically expanding populations experience similar eco-evolutionary challenges, including founder events, genetic bottlenecks, and novel environments. Theory predicts that reduced genetic diversity resulting from such phenomena limits the success of introduced populations. Using 1900 SNPs obtained from restriction-site-associated DNA sequencing, we evaluated hypotheses related to the invasion history and connectivity of an invasive population of Tench (Tinca tinca), a Eurasian freshwater fish that has been expanding geographically in eastern North America for three decades. Consistent with the reported history of a single introduction event, our findings suggest that multiple introductions from distinct genetic sources are unlikely as Tench had a small effective population size (~114 [95% CI = 106-123] individuals), no strong population subdivision across time and space, and evidence of a recent genetic bottleneck. The large genetic neighbourhood size (220 km) and weak within-population genetic substructure suggested high connectivity across the invaded range, despite the relatively large area occupied. There was some evidence for a small decay in genetic diversity as the species expanded northward, but not southward, into new habitats. As eradicating the species within a ~112 km radius would be necessary to prevent recolonization, eradicating Tench is likely not feasible at watershed-and possibly local-scales. Management should instead focus on reducing abundance in priority conservation areas to mitigate adverse impacts. Our study indicates that introduced populations can thrive and exhibit relatively high levels of genetic diversity despite severe bottlenecks (<1.5% of the ancestral effective population size) and suggests that landscape heterogeneity and population demographics can generate variability in spatial patterns of genetic diversity within a single range expansion.

Exploiting spatial dimensions to enable parallelized continuous directed evolution

Current strategies to improve the throughput of continuous directed evolution technologies often involve complex mechanical fluid-controlling system or robotic platforms, which limits their popularization and application in general laboratories. Inspired by our previous study on bacterial range expansion, in this study, we report a system termed SPACE for rapid and extensively parallelizable evolution of biomolecules by introducing spatial dimensions into the landmark phage-assisted continuous evolution system. Specifically, M13 phages and chemotactic Escherichia coli cells were closely inoculated onto a semisolid agar. The phages came into contact with the expanding front of the bacterial range, and then comigrated with the bacteria. This system leverages competition over space, wherein evolutionary progress is closely associated with the production of spatial patterns, allowing the emergence of improved or new protein functions. In a prototypical problem, SPACE remarkably simplified the process and evolved the promoter recognition of T7 RNA polymerase (RNAP) to a library of 96 random sequences in parallel. These results establish SPACE as a simple, easy to implement, and massively parallelizable platform for continuous directed evolution in general laboratories.

The role of pathogen-mediated insect superabundance in the East African emergence of a plant virus

One of the major crops for food security is cassava. Superabundant Bemisia tabaci whitefly, comprising unusually high landscape populations of the insect, have been implicated in cassava virus emergence. Studies have been unable to select from several hypotheses, however, as to the dynamic drivers of superabundant whitefly associated with the emergence in East Africa of severe cassava mosaic disease. One possibility is that pathogenic modification of infected plants can itself increase the growth of insect vector colonies on infected plants.Through the modelling of population processes at the landscape scale we introduce a framework for analysing patterns in the association of disease and insect waves.Our analyses demonstrate the role of pathogen-mediated insect superabundance in a plant disease invasion. Synthesis. An elevated abundance of insects at the landscape scale is frequently implicated in invasions of the plant pathogens that they carry. We advance ecological understanding of plant disease invasions by showing how landscape data can be used to investigate the causes of insect vector superabundance.

Predicting spatially heterogeneous invasive spread: Pyracantha angustifolia invading a dry Andean valley in northern Argentina

Understanding the drivers of invasive species spread is key to designing optimal management programmes for controlling them. Population models, parameterized from demographic and dispersal data, are useful for simulating invasion when long-term observations are lacking, which is particularly important when invasions are occurring rapidly. We aimed to understand how the spread of Pyracantha angustifolia (an invasive ornamental shrub) is influenced by heterogeneity in local conditions within a dry inter-montane valley in northwestern Argentina that is currently experiencing rapid urban expansion and concomitant reduction in cattle farming. Field surveys demonstrated that individuals growing in shrublands have a higher fecundity than those in grassland and rocky habitats. Furthermore, grazing pressure of cattle and horses substantially reduces fecundity and impacts the relationship between age and fecundity. We incorporated our field-estimated habitat and age-specific fecundities into both analytical integro-difference equation models and individual-based models to predict rates of spread across each habitat type in the presence or absence of grazing pressure. Results indicate that the rate of spread would be substantially higher in shrublands (by up to 33%) and highlight that grazing pressure can, through direct consumption of seeds by livestock, substantially depress the rate of spread (by up to 53%). These results suggest that shrubland areas in the valley are most vulnerable to invasion and that grazing may help to reduce the impact of Pyracantha. This suggests that a synergism between urban expansion and reduction in domestic livestock density contributes to the spread of an invasive species in this environment.

The ecological roles of bacterial chemotaxis

How bacterial chemotaxis is performed is much better understood than why. Traditionally, chemotaxis has been understood as a foraging strategy by which bacteria enhance their uptake of nutrients and energy, yet it has remained puzzling why certain less nutritious compounds are strong chemoattractants and vice versa. Recently, we have gained increased understanding of alternative ecological roles of chemotaxis, such as navigational guidance in colony expansion, localization of hosts or symbiotic partners and contribution to microbial diversity by the generation of spatial segregation in bacterial communities. Although bacterial chemotaxis has been observed in a wide range of environmental settings, insights into the phenomenon are mostly based on laboratory studies of model organisms. In this Review, we highlight how observing individual and collective migratory behaviour of bacteria in different settings informs the quantification of trade-offs, including between chemotaxis and growth. We argue that systematically mapping when and where bacteria are motile, in particular by transgenerational bacterial tracking in dynamic environments and in situ approaches from guts to oceans, will open the door to understanding the rich interplay between metabolism and growth and the contribution of chemotaxis to microbial life.

Methods to Evaluate Bacterial Motility and Its Role in Bacterial–Host Interactions

Bacterial motility is a widespread characteristic that can provide several advantages for the cell, allowing it to move towards more favorable conditions and enabling host-associated processes such as colonization. There are different bacterial motility types, and their expression is highly regulated by the environmental conditions. Because of this, methods for studying motility under realistic experimental conditions are required. A wide variety of approaches have been developed to study bacterial motility. Here, we present the most common techniques and recent advances and discuss their strengths as well as their limitations. We classify them as macroscopic or microscopic and highlight the advantages of three-dimensional imaging in microscopic approaches. Lastly, we discuss methods suited for studying motility in bacterial–host interactions, including the use of the zebrafish model.

A traveling-wave solution for bacterial chemotaxis with growth

Bacterial cells navigate their environment by directing their movement along chemical gradients. This process, known as chemotaxis, can promote the rapid expansion of bacterial populations into previously unoccupied territories. However, despite numerous experimental and theoretical studies on this classical topic, chemotaxis-driven population expansion is not understood in quantitative terms. Building on recent experimental progress, we here present a detailed analytical study that provides a quantitative understanding of how chemotaxis and cell growth lead to rapid and stable expansion of bacterial populations. We provide analytical relations that accurately describe the dependence of the expansion speed and density profile of the expanding population on important molecular, cellular, and environmental parameters. In particular, expansion speeds can be boosted by orders of magnitude when the environmental availability of chemicals relative to the cellular limits of chemical sensing is high. Analytical understanding of such complex spatiotemporal dynamic processes is rare. Our analytical results and the methods employed to attain them provide a mathematical framework for investigations of the roles of taxis in diverse ecological contexts across broad parameter regimes.

Spatial modulation of individual behaviors enables an ordered structure of diverse phenotypes during bacterial group migration

Coordination of diverse individuals often requires sophisticated communications and high-order computational abilities. Microbial populations can exhibit diverse individualistic behaviors, and yet can engage in collective migratory patterns with a spatially sorted arrangement of phenotypes. However, it is unclear how such spatially sorted patterns emerge from diverse individuals without complex computational abilities. Here, by investigating the single-cell trajectories during group migration, we discovered that, despite the constant migrating speed of a group, the drift velocities of individual bacteria decrease from the back to the front. With a Langevin-type modeling framework, we showed that this decreasing profile of drift velocities implies the spatial modulation of individual run-and-tumble random motions, and enables the bacterial population to migrate as a pushed wave front. Theoretical analysis and stochastic simulations further predicted that the pushed wave front can help a diverse population to stay in a tight group, while diverse individuals perform the same type of mean reverting processes around centers orderly aligned by their chemotactic abilities. This mechanism about the emergence of orderly collective migration from diverse individuals is experimentally demonstrated by titration of bacterial chemoreceptor abundance. These results reveal a simple computational principle for emergent ordered behaviors from heterogeneous individuals.

Organisms living in large groups often have to move together in order to navigate, forage for food, and increase their roaming range. Such groups are often made up of distinct individuals that must integrate their different behaviors in order to migrate in the same direction at a similar pace. For instance, for the bacteria Escherichia coli to travel as a condensed group, they must coordinate their response to a set of chemical signals called chemoattractants that tell them where to go.

The chemoattractants surrounding the bacteria are unequally distributed so that there is more of them at the front than the back of the group. During migration, each bacterium moves towards this concentration gradient in a distinct way, spontaneously rotating its direction in a ‘run-and-tumble’ motion that guides it towards areas where there are high levels of these chemical signals. In addition to this variability, how well individual bacteria are able to swim up the gradient also differs within the population. Bacteria that are better at sensing the chemoattractant gradient are placed at the front of the group, while those that are worst are shifted towards the back. This spatial arrangement is thought to help the bacteria migrate together. But how E. coli organize themselves in to this pattern is unclear, especially as they cannot communicate directly with one another and display such diverse, randomized behaviors.

To help answer this question, Bai, He et al. discovered a general principle that describes how single bacterial cells move within a group. The results showed that E. coli alter their run-and-tumble motion depending on where they reside within the population: individuals at the rear drift faster so they can catch up with the group, while those leading the group drift slower to draw themselves back. This ‘reversion behavior’ allows the migrating bacteria to travel at a constant speed around a mean position relative to the group.

A cell’s drifting speed is determined by how well it moves towards the chemoattractant and its response to the concentration gradient. As a result, the mean position around which the bacterium accelerates or deaccelerates will vary depending on how sensitive it is to the chemoattractant gradient. The E. coli therefore spatially arrange themselves so that the more sensitive bacteria are located at the front of the group where the gradient is shallower; and cells that are less sensitive are located towards the back where the gradient is steeper.

These findings suggest a general principle for how bacteria form ordered patterns whilst migrating as a collective group. This behavior could also apply to other populations of distinct individuals, such as ants following a trail or flocks of birds migrating in between seasons.

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.

Diffusive wave dynamics beyond the continuum limit

Scientists have observed and studied diffusive waves in contexts as disparate as population genetics and cell signaling. Often, these waves are propagated by discrete entities or agents, such as individual cells in the case of cell signaling. For a broad class of diffusive waves, we characterize the transition between the collective propagation of diffusive waves, in which the wave speed is well described by continuum theory, and the propagation of diffusive waves by individual agents. We show that this transition depends heavily on the dimensionality of the system in which the wave propagates and that disordered systems yield dynamics largely consistent with lattice systems. In some system dimensionalities, the intuition that closely packed sources more accurately mimic a continuum can be grossly violated.

The spatial organization of microbial communities during range expansion

Microbes in nature often live in dense and diverse communities exhibiting a variety of spatial structures. Microbial range expansion is a universal ecological process that enables populations to form spatial patterns. It can be driven by both passive and active processes, for example, mechanical forces from cell growth and bacterial motility. In this review, we provide a taste of recent creative and sophisticated efforts being made to address basic questions in spatial ecology and pattern formation during range expansion. We especially highlight the role of motility to shape community structures, and discuss the research challenges and future directions.

Dynamics of diffusive cell signaling relays

In biological contexts as diverse as development, apoptosis, and synthetic microbial consortia, collections of cells or subcellular components have been shown to overcome the slow signaling speed of simple diffusion by utilizing diffusive relays, in which the presence of one type of diffusible signaling molecule triggers participation in the emission of the same type of molecule. This collective effect gives rise to fast-traveling diffusive waves. Here, in the context of cell signaling, we show that system dimensionality – the shape of the extracellular medium and the distribution of cells within it – can dramatically affect the wave dynamics, but that these dynamics are insensitive to details of cellular activation. As an example, we show that neutrophil swarming experiments exhibit dynamical signatures consistent with the proposed signaling motif. We further show that cell signaling relays generate much steeper concentration profiles than does simple diffusion, which may facilitate neutrophil chemotaxis.

Eco-evolutionary dynamics of range expansion

Understanding the movement of species' ranges is a classic ecological problem that takes on urgency in this era of global change. Historically treated as a purely ecological process, range expansion is now understood to involve eco-evolutionary feedbacks due to spatial genetic structure that emerges as populations spread. We synthesize empirical and theoretical work on the eco-evolutionary dynamics of range expansion, with emphasis on bridging directional, deterministic processes that favor evolved increases in dispersal and demographic traits with stochastic processes that lead to the random fixation of alleles and traits. We develop a framework for understanding the joint influence of these processes in changing the mean and variance of expansion speed and its underlying traits. Our synthesis of recent laboratory experiments supports the consistent role of evolution in accelerating expansion speed on average, and highlights unexpected diversity in how evolution can influence variability in speed: results not well predicted by current theory. We discuss and evaluate support for three classes of modifiers of eco-evolutionary range dynamics (landscape context, trait genetics, and biotic interactions), identify emerging themes, and suggest new directions for future work in a field that stands to increase in relevance as populations move in response to global change.

Dispersal evolution diminishes the negative density dependence in dispersal

In many organisms, dispersal varies with the local population density. Such patterns of density-dependent dispersal (DDD) are expected to shape the dynamics, spatial spread, and invasiveness of populations. Despite their ecological importance, empirical evidence for the evolution of DDD patterns remains extremely scarce. This is especially relevant because rapid evolution of dispersal traits has now been empirically confirmed in several taxa. Changes in DDD of dispersing populations could help clarify not only the role of DDD in dispersal evolution, but also the possible pattern of subsequent range expansion. Here, we investigate the relationship between dispersal evolution and DDD using a long-term experimental evolution study on Drosophila melanogaster. We compared the DDD patterns of four dispersal-selected populations and their non-selected controls. The control populations showed negative DDD, which was stronger in females than in males. In contrast, the dispersal-selected populations showed DDD, where neither males nor females exhibited DDD. We compare our results with previous evolutionary predictions that focused largely on positive DDD, and highlight how the direction of evolutionary change depends on the initial DDD pattern of a population. Finally, we discuss the implications of DDD evolution for spatial ecology and evolution.

Mutualistic cross-feeding in microbial systems generates bistability via an Allee effect

In microbial ecosystems, species not only compete for common resources but may also display mutualistic interactions as a result from metabolic cross-feeding. Such mutualism can lead to bistability. Depending on the initial population sizes, species will either survive or go extinct. Various phenomenological models have been suggested to describe bistability in mutualistic systems. However, these models do not account for interaction mediators such as nutrients. In contrast, nutrient-explicit models do not provide an intuitive understanding of what causes bistability. Here, we reduce a theoretical nutrient-explicit model of two mutualistic cross-feeders in a chemostat, uncovering an explicit relation to a growth model with an Allee effect. We show that the dilution rate in the chemostat leads to bistability by turning a weak Allee effect into a strong Allee effect. This happens as long as there is more production than consumption of cross-fed nutrients. Thanks to the explicit relationship of the reduced model with the underlying experimental parameters, these results allow to predict the biological conditions that sustain or prevent the survival of mutualistic species.

Bacterial coexistence driven by motility and spatial competition

Elucidating elementary mechanisms that underlie bacterial diversity is central to ecology^1,2 and microbiome research³. Bacteria are known to coexist by metabolic specialization⁴, cooperation⁵ and cyclic warfare^6-8. Many species are also motile⁹, which is studied in terms of mechanism^10,11, benefit^12,13, strategy^14,15, evolution^16,17 and ecology^18,19. Indeed, bacteria often compete for nutrient patches that become available periodically or by random disturbances^2,20,21. However, the role of bacterial motility in coexistence remains unexplored experimentally. Here we show that-for mixed bacterial populations that colonize nutrient patches-either population outcompetes the other when low in relative abundance. This inversion of the competitive hierarchy is caused by active segregation and spatial exclusion within the patch: a small fast-moving population can outcompete a large fast-growing population by impeding its migration into the patch, while a small fast-growing population can outcompete a large fast-moving population by expelling it from the initial contact area. The resulting spatial segregation is lost for weak growth-migration trade-offs and a lack of virgin space, but is robust to population ratio, density and chemotactic ability, and is observed in both laboratory and wild strains. These findings show that motility differences and their trade-offs with growth are sufficient to promote diversity, and suggest previously undescribed roles for motility in niche formation and collective expulsion-containment strategies beyond individual search and survival.

Evolutionary Rescue from a Wave of Biological Invasion

Evolution can potentially rescue populations from being driven extinct by biological invasions, but predictions for this occurrence are generally lacking. Here I derive theoretical predictions for evolutionary rescue of a resident population experiencing invasion from an introduced competitor that spreads over its introduced range as a traveling spatial wave that displaces residents. I compare the likelihood of evolutionary rescue from invasion for two modes of competition: exploitation and interference competition. I find that, all else equal, evolutionary rescue is less effective at preventing extinction caused by interference-driven invasions than by exploitation-driven invasions. Rescue from interference-driven invasions is, surprisingly, independent of invader dispersal rate or the speed of invasion and is more weakly dependent on range size than in the exploitation-driven case. In contrast, rescue from exploitation-driven invasions strongly depends on range size and is less likely during fast invasions. The results presented here have potential applications for conserving endemic species from nonnative invaders and for ensuring extinction of pests using intentionally introduced biocontrol agents.

Cooperation mitigates diversity loss in a spatially expanding microbial population

Spatially expanding populations lose genetic diversity rapidly because of repeated bottlenecks formed at the colonization front. However, the rate of diversity loss depends on the specifics of the expanding population, such as its growth and dispersal dynamics. We have previously demonstrated that changing the amount of within-species cooperation leads to a qualitative transition in the nature of expansion from pulled (driven by growth at the low-density tip) to pushed (driven by migration from the high-density region at the front, but behind the tip). Here we demonstrate experimentally that pushed waves result in strongly reduced genetic drift during range expansions, thus preserving genetic diversity in the newly colonized region.

The evolution and potentially even the survival of a spatially expanding population depends on its genetic diversity, which can decrease rapidly due to a serial founder effect. The strength of the founder effect is predicted to depend strongly on the details of the growth dynamics. Here, we probe this dependence experimentally using a single microbial species, Saccharomyces cerevisiae, expanding in multiple environments that induce varying levels of cooperativity during growth. We observe a drastic reduction in diversity during expansions when yeast grows noncooperatively on simple sugars, but almost no loss of diversity when cooperation is required to digest complex metabolites. These results are consistent with theoretical expectations: When cells grow independently from each other, the expansion proceeds as a pulled wave driven by growth at the low-density tip of the expansion front. Such populations lose diversity rapidly because of the strong genetic drift at the expansion edge. In contrast, diversity loss is substantially reduced in pushed waves that arise due to cooperative growth. In such expansions, the low-density tip of the front grows much more slowly and is often reseeded from the genetically diverse population core. Additionally, in both pulled and pushed expansions, we observe a few instances of abrupt changes in allele fractions due to rare fluctuations of the expansion front and show how to distinguish such rapid genetic drift from selective sweeps.

Time scales and wave formation in non-linear spatial public goods games

The co-evolutionary dynamics of competing populations can be strongly affected by frequency-dependent selection and spatial population structure. As co-evolving populations grow into a spatial domain, their initial spatial arrangement and their growth rate differences are important factors that determine the long-term outcome. We here model producer and free-rider co-evolution in the context of a diffusive public good (PG) that is produced by the producers at a cost but evokes local concentration-dependent growth benefits to all. The benefit of the PG can be non-linearly dependent on public good concentration. We consider the spatial growth dynamics of producers and free-riders in one, two and three dimensions by modeling producer cell, free-rider cell and public good densities in space, driven by the processes of birth, death and diffusion (cell movement and public good distribution). Typically, one population goes extinct, but the time-scale of this process varies with initial conditions and the growth rate functions. We establish that spatial variation is transient regardless of dimensionality, and that structured initial conditions lead to increasing times to get close to an extinction state, called ε-extinction time. Further, we find that uncorrelated initial spatial structures do not influence this ε-extinction time in comparison to a corresponding well-mixed (non-spatial) system. In order to estimate the ε-extinction time of either free-riders or producers we derive a slow manifold solution. For invading populations, i.e. for populations that are initially highly segregated, we observe a traveling wave, whose speed can be calculated. Our results provide quantitative predictions for the transient spatial dynamics of cooperative traits under pressure of extinction.

Evolutionary public good (PG) games capture the essence of production of growth-beneficial factors that are vulnerable to exploitation by free-riders who do not carry the cost of production. PGs emerge in cellular populations, for example in growing bacteria and cancer cells. We study the eco-evolutionary dynamics of a PG in populations that grow in space. In our model, PG-producer cells and free-rider cells can grow according to their own birth and death rates. Co-evolution occurs due to public good-driven surplus in the intrinsic growth rates at a cost to producers. A net growth rate-benefit to free-riders leads to the well-known tragedy of the commons in which producers go extinct. What is often omitted from discussions is the time scale on which this extinction can occur, especially in spatial populations. Here, we derive analytical estimates of the ε-extinction time in different spatial settings. As we do not consider a stochastic process, the ε-extinction time captures the time needed to approach an extinction state. We identify spatial scenarios in which extinction takes long enough such that the tragedy of the commons never occurs within a meaningful lifetime of the system. Using numerical simulations we analyze the deviations from our analytical predictions.

Evolution at the Edge of Expanding Populations

Predicting the evolution of expanding populations is critical to controlling biological threats such as invasive species and cancer metastasis. Expansion is primarily driven by reproduction and dispersal, but nature abounds with examples of evolution where organisms pay a reproductive cost to disperse faster. When does selection favor this "survival of the fastest"? We searched for a simple rule, motivated by evolution experiments where swarming bacteria evolved into a hyperswarmer mutant that disperses ∼100% faster but pays a growth cost of ∼10% to make many copies of its flagellum. We analyzed a two-species model based on the Fisher equation to explain this observation: the population expansion rate (v) results from an interplay of growth (r) and dispersal (D) and is independent of the carrying capacity: v = 2 ( rD ) 1 / 2 . A mutant can take over the edge only if its expansion rate (v2) exceeds the expansion rate of the established species (v1); this simple condition ( v 2 > v 1 ) determines the maximum cost in slower growth that a faster mutant can pay and still be able to take over. Numerical simulations and time-course experiments where we tracked evolution by imaging bacteria suggest that our findings are general: less favorable conditions delay but do not entirely prevent the success of the fastest. Thus, the expansion rate defines a traveling wave fitness, which could be combined with trade-offs to predict evolution of expanding populations.

Spatial coordination in a mutually beneficial bacterial community enhances its antibiotic resistance

Microbial communities can survive in complex and variable environments by using different cooperative strategies. However, the behaviors of these mutuality formed communities remain poorly understood, particularly with regard to the characteristics of spatial cooperation. Here, we selected two Escherichia coli strains, designated as the nutrition provider and the antibiotic protector, respectively, for construction of a mutually beneficial bacterial community that could be used to study these behaviors. We found that in addition to the functional mutualism, the two strains also cooperated through their spatial distribution. Under antibiotic pressure, the bacterial distribution changed to yield different spatial distributions, which resulted in community growth advantages beyond functional cooperation. The mutualistic behavior of these two strains suggested that similar communities could also use variations in spatial distribution to improve their survival rates in a natural environment or under the action of antibiotics.

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.

Bet-hedging strategies in expanding populations

In ecology, species can mitigate their extinction risks in uncertain environments by diversifying individual phenotypes. This observation is quantified by the theory of bet-hedging, which provides a reason for the degree of phenotypic diversity observed even in clonal populations. Bet-hedging in well-mixed populations is rather well understood. However, many species underwent range expansions during their evolutionary history, and the importance of phenotypic diversity in such scenarios still needs to be understood. In this paper, we develop a theory of bet-hedging for populations colonizing new, unknown environments that fluctuate either in space or time. In this case, we find that bet-hedging is a more favorable strategy than in well-mixed populations. For slow rates of variation, temporal and spatial fluctuations lead to different outcomes. In spatially fluctuating environments, bet-hedging is favored compared to temporally fluctuating environments. In the limit of frequent environmental variation, no opportunity for bet-hedging exists, regardless of the nature of the environmental fluctuations. For the same model, bet-hedging is never an advantageous strategy in the well-mixed case, supporting the view that range expansions strongly promote diversification. These conclusions are robust against stochasticity induced by finite population sizes. Our findings shed light on the importance of phenotypic heterogeneity in range expansions, paving the way to novel approaches to understand how biodiversity emerges and is maintained.

Ecological populations are often exposed to unpredictable and variable environmental conditions. A number of strategies have evolved to cope with such uncertainty. One of them is stochastic phenotypic switching, by which some individuals in the community are enabled to tackle adverse conditions, even at the price of reducing overall growth in the short term. In this paper, we study the effectiveness of these “bet-hedging” strategies for a population in the process of colonizing new territory. We show that bet-hedging is more advantageous when the environment varies spatially rather than temporally, and infrequently rather than frequently.

The three regimes of spatial recovery

An enduring challenge for ecology is identifying the drivers of ecosystem and population stability. In a spatially explicit context, key features to consider are landscape spatial structure, local interactions, and dispersal. Substantial work has been done on each of these features as a driver of stability, but little is known on the interplay between them. Missing has been a more integrative approach, able to map and identify different dynamical regimes, predicting a system's response to perturbations. Here we first consider a simple scenario, i.e., the recovery of a homogeneous metapopulation from a single localized pulse disturbance. The analysis of this scenario reveals three fundamental recovery regimes: Isolated Regime when dispersal is not significant, Rescue Regime when dispersal mediates recovery, and Mixing Regime when perturbations spread throughout the system. Despite its simplicity, our approach leads to remarkably general predictions. These include the qualitatively different outcomes of various scenarios of habitat fragmentation, the surprising benefits of local extinctions on population persistence at the transition between regimes, and the productivity shifts of metacommunities in a changing environment. This study thus provides context to known results and insight into future directions of research.

Delays in Fitness Adjustment Can Lead to Coexistence of Hierarchically Interacting Species

Organisms that exploit different environments may experience a stochastic delay in adjusting their fitness when they switch habitats. We study two such organisms whose fitness is determined by the species composition of the local environment, as they interact through a public good. We show that a delay in the fitness adjustment can lead to the coexistence of the two species in a metapopulation, although the faster-growing species always wins in well-mixed competition experiments. Coexistence is favored over wide parameter ranges and is independent of spatial clustering. It arises when species are heterogeneous in their fitness and can keep each other balanced.

Chirality provides a direct fitness advantage and facilitates intermixing in cellular aggregates

Chirality in shape and motility can evolve rapidly in microbes and cancer cells. To determine how chirality affects cell fitness, we developed a model of chiral growth in compact aggregates such as microbial colonies and solid tumors. Our model recapitulates previous experimental findings and shows that mutant cells can invade by increasing their chirality or switching their handedness. The invasion results either in a takeover or stable coexistence between the mutant and the ancestor depending on their relative chirality. For large chiralities, the coexistence is accompanied by strong intermixing between the cells, while spatial segregation occurs otherwise. We show that the competition within the aggregate is mediated by bulges in regions where the cells with different chiralities meet. The two-way coupling between aggregate shape and natural selection is described by the chiral Kardar-Parisi-Zhang equation coupled to the Burgers’ equation with multiplicative noise. We solve for the key features of this theory to explain the origin of selection on chirality. Overall, our work suggests that chirality could be an important ecological trait that mediates competition, invasion, and spatial structure in cellular populations.

Is it better to be left- or right-handed? The answer depends on whether the goal is making a handshake or winning a boxing match. The need for coordination favors the handedness of the majority, but being different could also provide an advantage. The same rules could apply to microbial colonies and cancer tumors. Like humans, cells often have handedness (chirality) that reflects the lack of mirror symmetry in their shapes or movement patterns. We find that cells gain a substantial fitness advantage by either increasing the magnitude of their chirality or switching to the opposite handedness. Selection for specific chirality can overcome differences in growth rate and is mediated by the formation of bulges along the colony edge in regions where cells with different chiralities meet.

Role of carbon allocation efficiency in the temperature dependence of autotroph growth rates

Relating the temperature dependence of photosynthetic biomass production to underlying metabolic rates in autotrophs is crucial for predicting the effects of climatic temperature fluctuations on the carbon balance of ecosystems. We present a mathematical model that links thermal performance curves (TPCs) of photosynthesis, respiration, and carbon allocation efficiency to the exponential growth rate of a population of photosynthetic autotroph cells. Using experiments with the green alga, Chlorella vulgaris, we apply the model to show that the temperature dependence of carbon allocation efficiency is key to understanding responses of growth rates to warming at both ecological and longer-term evolutionary timescales. Finally, we assemble a dataset of multiple terrestrial and aquatic autotroph species to show that the effects of temperature-dependent carbon allocation efficiency on potential growth rate TPCs are expected to be consistent across taxa. In particular, both the thermal sensitivity and the optimal temperature of growth rates are expected to change significantly due to temperature dependence of carbon allocation efficiency alone. Our study provides a foundation for understanding how the temperature dependence of carbon allocation determines how population growth rates respond to temperature.

Fluctuations uncover a distinct class of traveling waves

Epidemics, flame propagation, and cardiac rhythms are classic examples of reaction-diffusion waves that describe a switch from one alternative state to another. Only two types of waves are known: pulled, driven by the leading edge, and pushed, driven by the bulk of the wave. Here, we report a distinct class of semipushed waves for which both the bulk and the leading edge contribute to the dynamics. These hybrid waves have the kinetics of pushed waves, but exhibit giant fluctuations similar to pulled waves. The transitions between pulled, semipushed, and fully pushed waves occur at universal ratios of the wave velocity to the Fisher velocity. We derive these results in the context of a species invading a new habitat by examining front diffusion, rate of diversity loss, and fluctuation-induced corrections to the expansion velocity. All three quantities decrease as a power law of the population density with the same exponent. We analytically calculate this exponent, taking into account the fluctuations in the shape of the wave front. For fully pushed waves, the exponent is -1, consistent with the central limit theorem. In semipushed waves, however, the fluctuations average out much more slowly, and the exponent approaches 0 toward the transition to pulled waves. As a result, a rapid loss of genetic diversity and large fluctuations in the position of the front occur, even for populations with cooperative growth and other forms of an Allee effect. The evolutionary outcome of spatial spreading in such populations could therefore be less predictable than previously thought.

Multiple scales in metapopulations of public goods producers

Multiple scales in metapopulations can give rise to paradoxical behavior: in a conceptual model for a public goods game, the species associated with a fitness cost due to the public good production can be stabilized in the well-mixed limit due to the mere existence of these scales. The scales in this model involve a length scale corresponding to separate patches, coupled by mobility, and separate time scales for reproduction and interaction with a local environment. Contrary to the well-mixed high mobility limit, we find that for low mobilities, the interaction rate progressively stabilizes this species due to stochastic effects, and that the formation of spatial patterns is not crucial for this stabilization.

Experimental evidence that density dependence strongly influences plant invasions through fragmented landscapes

Populations of range expanding species encounter patches of both favorable and unfavorable habitat as they spread across landscapes. Theory shows that increasing patchiness slows the spread of populations modeled with continuously varying population density when dispersal is not influence by the environment or individual behavior. However, as is found in uniformly favorable landscapes, spread remains driven by fecundity and dispersal from low density individuals at the invasion front. In contrast, when modeled populations are composed of discrete individuals, patchiness causes populations to build up to high density before dispersing past unsuitable habitat, introducing an important influence of density dependence on spread velocity. To test the hypothesized interaction between habitat patchiness and density dependence, we simultaneously manipulated these factors in a greenhouse system of annual plants spreading through replicated experimental landscapes. We found that increasing the size of gaps and amplifying the strength of density dependence both slowed spread velocity, but contrary to predictions, the effect of amplified density dependence was similar across all landscape types. Our results demonstrate that the discrete nature of individuals in spreading populations has a strong influence on how both landscape patchiness and density dependence influence spread through demographic and dispersal stochasticity. Both finiteness and landscape structure should be critical components to theoretical predictions of future spread for range expanding native species or invasive species colonizing new habitat.

The spatial and metabolic basis of colony size variation

Spatial structure impacts microbial growth and interactions, with ecological and evolutionary consequences. It is therefore important to quantitatively understand how spatial proximity affects interactions in different environments. We tested how proximity influences colony size when either Escherichia coli or Salmonella enterica are grown on various carbon sources. The importance of colony location changed with species and carbon source. Spatially explicit, genome-scale metabolic modeling recapitulated observed colony size variation. Competitors that determine territory size, according to Voronoi diagrams, were the most important drivers of variation in colony size. However, the relative importance of different competitors changed through time. Further, the effect of location increased when colonies took up resources quickly relative to the diffusion of limiting resources. These analyses made it apparent that the importance of location was smaller than expected for experiments with S. enterica growing on glucose. The accumulation of toxic byproducts appeared to limit the growth of large colonies and reduced variation in colony size. Our work provides an experimentally and theoretically grounded understanding of how location interacts with metabolism and diffusion to influence microbial interactions.

Microbial expansion-collision dynamics promote cooperation and coexistence on surfaces

Microbes colonizing a surface often experience colony growth dynamics characterized by an initial phase of spatial clonal expansion followed by collision between neighboring colonies to form potentially genetically heterogeneous boundaries. For species with life cycles consisting of repeated surface colonization and dispersal, these spatially explicit "expansion-collision dynamics" generate periodic transitions between two distinct selective regimes, "expansion competition" and "boundary competition," each one favoring a different growth strategy. We hypothesized that this dynamic could promote stable coexistence of expansion- and boundary-competition specialists by generating time-varying, negative frequency-dependent selection that insulates both types from extinction. We tested this experimentally in budding yeast by competing an exoenzyme secreting "cooperator" strain (expansion-competition specialists) against nonsecreting "defectors" (boundary-competition specialists). As predicted, we observed cooperator-defector coexistence or cooperator dominance with expansion-collision dynamics, but only defector dominance otherwise. Also as predicted, the steady-state frequency of cooperators was determined by colonization density (the average initial cell-cell distance) and cost of cooperation. Lattice-based spatial simulations give good qualitative agreement with experiments, supporting our hypothesis that expansion-collision dynamics with costly public goods production is sufficient to generate stable cooperator-defector coexistence. This mechanism may be important for maintaining public-goods cooperation and conflict in microbial pioneer species living on surfaces.

Spatial dynamics of synthetic microbial mutualists and their parasites

A major force contributing to the emergence of novelty in nature is the presence of cooperative interactions, where two or more components of a system act in synergy, sometimes leading to higher-order, emergent phenomena. Within molecular evolution, the so called hypercycle defines the simplest model of an autocatalytic cycle, providing major theoretical insights on the evolution of cooperation in the early biosphere. These closed cooperative loops have also inspired our understanding of how catalytic loops appear in ecological systems. In both cases, hypercycle and ecological cooperative loops, the role played by space seems to be crucial for their stability and resilience against parasites. However, it is difficult to test these ideas in natural ecosystems, where time and spatial scales introduce considerable limitations. Here, we use engineered bacteria as a model system to a variety of environmental scenarios identifying trends that transcend the specific model system, such an enhanced genetic diversity in environments requiring mutualistic interactions. Interestingly, we show that improved environments can slow down mutualistic range expansions as a result of genetic drift effects preceding local resource depletion. Moreover, we show that a parasitic strain is excluded from the population during range expansions (which acknowledges a classical prediction). Nevertheless, environmental deterioration can reshape population interactions, this same strain becoming part of a three-species mutualistic web in scenarios in which the two-strain mutualism becomes non functional. The evolutionary and ecological implications for the design of synthetic ecosystems are outlined.

Spatial gene drives and pushed genetic waves

Gene drives have the potential to rapidly replace a harmful wild-type allele with a gene drive allele engineered to have desired functionalities. However, an accidental or premature release of a gene drive construct to the natural environment could damage an ecosystem irreversibly. Thus, it is important to understand the spatiotemporal consequences of the super-Mendelian population genetics before potential applications. Here, we use a reaction-diffusion model for sexually reproducing diploid organisms to study how a locally introduced gene drive allele spreads to replace the wild-type allele, although it possesses a selective disadvantage s > 0. Using methods developed by Barton and collaborators, we show that socially responsible gene drives require 0.5 < s < 0.697, a rather narrow range. In this "pushed wave" regime, the spatial spreading of gene drives will be initiated only when the initial frequency distribution is above a threshold profile called "critical propagule," which acts as a safeguard against accidental release. We also study how the spatial spread of the pushed wave can be stopped by making gene drives uniquely vulnerable ("sensitizing drive") in a way that is harmless for a wild-type allele. Finally, we show that appropriately sensitized drives in two dimensions can be stopped, even by imperfect barriers perforated by a series of gaps.

Density dependence in demography and dispersal generates fluctuating invasion speeds

Density dependence plays an important role in population regulation and is known to generate temporal fluctuations in population density. However, the ways in which density dependence affects spatial population processes, such as species invasions, are less understood. Although classical ecological theory suggests that invasions should advance at a constant speed, empirical work is illuminating the highly variable nature of biological invasions, which often exhibit nonconstant spreading speeds, even in simple, controlled settings. Here, we explore endogenous density dependence as a mechanism for inducing variability in biological invasions with a set of population models that incorporate density dependence in demographic and dispersal parameters. We show that density dependence in demography at low population densities-i.e., an Allee effect-combined with spatiotemporal variability in population density behind the invasion front can produce fluctuations in spreading speed. The density fluctuations behind the front can arise from either overcompensatory population growth or density-dependent dispersal, both of which are common in nature. Our results show that simple rules can generate complex spread dynamics and highlight a source of variability in biological invasions that may aid in ecological forecasting.