Genetic demixing and evolution in linear stepping stone models

Stability and gene strand bias of lambda prophages and chromosome organization in Escherichia coli

Temperate phage-mediated horizontal gene transfer is a potent driver of genetic diversity in the evolution of bacteria. Most lambdoid prophages in Escherichia coli are integrated into the chromosome with the same orientation with respect to the direction of chromosomal replication, and their location on the chromosome is far from homogeneous. To better understand these features, we studied the interplay between lysogenic and lytic states of phage lambda in both native and inverted integration orientations at the wild-type integration site as well as at other sites on the bacterial chromosome. Measurements of free phage released by spontaneous induction showed that the stability of lysogenic states is affected by location and orientation along the chromosome, with stronger effects near the origin of replication. Competition experiments and range expansions between lysogenic strains with opposite orientations and insertion loci indicated that there are no major differences in growth. Moreover, measurements of the level of transcriptional bursts of the cI gene coding for the lambda phage repressor using single-molecule fluorescence in situ hybridization resulted in similar levels of transcription for both orientations and prophage location. We postulate that the preference for a given orientation and location is a result of a balance between the maintenance of lysogeny and the ability to lyse.IMPORTANCEThe integration of genetic material of temperate bacterial viruses (phages) into the chromosomes of bacteria is a potent evolutionary force, allowing bacteria to acquire in one stroke new traits and restructure the information in their chromosomes. Puzzlingly, this genetic material is preferentially integrated in a particular orientation and at non-random sites on the bacterial chromosome. The work described here reveals that the interplay between the maintenance of the stability of the integrated phage, its ability to excise, and its localization along the chromosome plays a key role in setting chromosomal organization in Escherichia coli.

Mechanically-driven growth and competition in a Voronoi model of tissues

The mechanisms leading cells to acquire a fitness advantage and establish themselves in a population are paramount to understanding the development and growth of cancer. Although there are many works that study separately either the evolutionary dynamics or the mechanics of cancer, little has been done to couple evolutionary dynamics to mechanics. To address this question, we study a confluent model of tissue using a Self-Propelled Voronoi (SPV) model with stochastic growth rates that depend on the mechanical variables of the system. The SPV model is an out-of-equilibrium model of tissue derived from an energy functional that has a jamming/unjamming transition between solid-like and liquid-like states. By considering several scenarios of mutants invading a resident population in both phases, we determine the range of parameters that confer a fitness advantage and show that the preferred area and perimeter are the most relevant ones. We find that the liquid-like state is more resistant to invasion and show that the outcome of the competition can be determined from the simulation of a non-growing mixture. Moreover, a mean-field approximation can accurately predict the fate of a mutation affecting mechanical properties of a cell. Our results can be used to infer evolutionary dynamics from tissue images, understand cancer-suppressing effects of tissue mechanics, and even search for mechanics-based therapies.

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.

Fitness effects of a demography-dispersal trade-off in expandingSaccharomyces cerevisiaemats

Fungi expand in space and time to form complex multicellular communities. The mechanisms by which they do so can vary dramatically and determine the life-history and dispersal traits of expanding populations. These traits influence deterministic and stochastic components of evolution, resulting in complex eco-evolutionary dynamics during colony expansion. We perform experiments on budding yeast strains genetically engineered to display rough-surface and smooth-surface phenotypes in colony-like structures called 'mats'. Previously, it was shown that the rough-surface strain has a competitive advantage over the smooth-surface strain when grown on semi-solid media. We experimentally observe the emergence and expansion of segments with a distinct smooth-surface phenotype during rough-surface mat development. We propose a trade-off between dispersal and local carrying capacity to explain the relative fitness of these two phenotypes. Using a modified stepping-stone model, we demonstrate that this trade-off gives the high-dispersing, rough-surface phenotype a competitive advantage from standing variation, but that it inhibits this phenotype's ability to invade a resident smooth-surface population via mutation. However, the trade-off improves the ability of the smooth-surface phenotype to invade in rough-surface mats, replicating the frequent emergence of smooth-surface segments in experiments. Together, these computational and experimental findings advance our understanding of the complex eco-evolutionary dynamics of fungal mat expansion.

Genetic mixing and demixing on expanding spherical frontiers

Genetic fluctuation during range expansion is a key process driving evolution. When a bacterial population is expanding on a 2D surface, random fluctuations in the growth of the pioneers at the front line cause a strong demixing of genotypes. Even when there is no selective advantage, sectors of low genetic diversity are formed. Experimental studies of range expansions in surface-attached colonies of fluorescently labelled micro-organisms have contributed significantly to our understanding of fundamental evolutionary dynamics. However, experimental studies on genetic fluctuations in 3D range expansions have been sparse, despite their importance for tumour or biofilm development. We encapsulated populations of two fluorescent Escherichia coli strains in inoculation droplets (volumes [Formula: see text] nl). The confined ensemble of cells grew when embedded in a hydrogel-with nutrients-and developed 3D colonies with well-defined, sector-like regions. Using confocal laser scanning microscopy, we imaged the development of 3D colonies and the emergence of sectors. We characterized how cell concentration in the inoculation droplet controls sectors, growth rate, and the transition from branched colonies to quasi-spherical colonies. We further analysed how sectors on the surface change over time. We complement these experimental results with a modified 3D Eden growth model. The model in 3D spherical growth predicts a phase, where sectors are merging, followed by a steady increase (constant rate), and the experimentally analysed sectors were consistent with this prediction. Therefore, our results demonstrate qualitative differences between radial (2D) and spherical (3D) range expansions and their importance in gene fixation processes.

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.

Interfacial morphodynamics of proliferating microbial communities

In microbial communities, various cell types often coexist by occupying distinct spatial domains. What determines the shape of the interface between such domains-which in turn influences the interactions between cells and overall community function? Here, we address this question by developing a continuum model of a 2D spatially-structured microbial community with two distinct cell types. We find that, depending on the balance of the different cell proliferation rates and substrate friction coefficients, the interface between domains is either stable and smooth, or unstable and develops finger-like protrusions. We establish quantitative principles describing when these different interfacial behaviors arise, and find good agreement both with the results of previous experimental reports as well as new experiments performed here. Our work thus helps to provide a biophysical basis for understanding the interfacial morphodynamics of proliferating microbial communities, as well as a broader range of proliferating active systems.

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.

Specialization at an expanding front

As a population grows, spreading to new environments may favor specialization. In this paper, we introduce and explore a model for specialization at the front of a colony expanding synchronously into new territory. We show through numerical simulations that, by gaining fitness through accumulating mutations, progeny of the initial seed population can differentiate into distinct specialists. With competition and selection limited to the growth front, the emerging specialists first segregate into sectors, which then expand to dominate the entire population. We quantify the scaling of the fixation time with the size of the population and observe different behaviors corresponding to distinct universality classes: unbounded and bounded gains in fitness lead to superdiffusive (z=3/2) and diffusive (z=2) stochastic wanderings of the sector boundaries, respectively.

The Impact of Genetic Surfing on Neutral Genomic Diversity

Range expansions have been common in the history of most species. Serial founder effects and subsequent population growth at expansion fronts typically lead to a loss of genomic diversity along the expansion axis. A frequent consequence is the phenomenon of "gene surfing," where variants located near the expanding front can reach high frequencies or even fix in newly colonized territories. Although gene surfing events have been characterized thoroughly for a specific locus, their effects on linked genomic regions and the overall patterns of genomic diversity have been little investigated. In this study, we simulated the evolution of whole genomes during several types of 1D and 2D range expansions differing by the extent of migration, founder events, and recombination rates. We focused on the characterization of local dips of diversity, or "troughs," taken as a proxy for surfing events. We find that, for a given recombination rate, once we consider the amount of diversity lost since the beginning of the expansion, it is possible to predict the initial evolution of trough density and their average width irrespective of the expansion condition. Furthermore, when recombination rates vary across the genome, we find that troughs are over-represented in regions of low recombination. Therefore, range expansions can leave local and global genomic signatures often interpreted as evidence of past selective events. Given the generality of our results, they could be used as a null model for species having gone through recent expansions, and thus be helpful to correctly interpret many evolutionary biology studies.

Morphological instability and roughening of growing 3D bacterial colonies

The morphogenesis of two-dimensional bacterial colonies has been well studied. However, little is known about the colony morphologies of bacteria growing in three dimensions, despite the prevalence of three-dimensional environments (e.g., soil, inside hosts) as natural bacterial habitats. Using experiments on bacteria in granular hydrogel matrices, we find that dense multicellular colonies growing in three dimensions undergo a common morphological instability and roughen, adopting a characteristic broccoli-like morphology when they exceed a critical size. Analysis of a continuum “active fluid” model of the expanding colony reveals that this behavior originates from an interplay of competition for nutrients with growth-driven colony expansion, both of which vary spatially. These results shed light on the fundamental biophysical principles underlying growth in three dimensions.

How do growing bacterial colonies get their shapes? While colony morphogenesis is well studied in two dimensions, many bacteria grow as large colonies in three-dimensional (3D) environments, such as gels and tissues in the body or subsurface soils and sediments. Here, we describe the morphodynamics of large colonies of bacteria growing in three dimensions. Using experiments in transparent 3D granular hydrogel matrices, we show that dense colonies of four different species of bacteria generically become morphologically unstable and roughen as they consume nutrients and grow beyond a critical size—eventually adopting a characteristic branched, broccoli-like morphology independent of variations in the cell type and environmental conditions. This behavior reflects a key difference between two-dimensional (2D) and 3D colonies; while a 2D colony may access the nutrients needed for growth from the third dimension, a 3D colony inevitably becomes nutrient limited in its interior, driving a transition to unstable growth at its surface. We elucidate the onset of the instability using linear stability analysis and numerical simulations of a continuum model that treats the colony as an “active fluid” whose dynamics are driven by nutrient-dependent cellular growth. We find that when all dimensions of the colony substantially exceed the nutrient penetration length, nutrient-limited growth drives a 3D morphological instability that recapitulates essential features of the experimental observations. Our work thus provides a framework to predict and control the organization of growing colonies—as well as other forms of growing active matter, such as tumors and engineered living materials—in 3D environments.

Growth instabilities shape morphology and genetic diversity of microbial colonies

Cellular populations assume an incredible variety of shapes ranging from circular molds to irregular tumors. While we understand many of the mechanisms responsible for these spatial patterns, little is known about how the shape of a population influences its ecology and evolution. Here, we investigate this relationship in the context of microbial colonies grown on hard agar plates. This a well-studied system that exhibits a transition from smooth circular disks to more irregular and rugged shapes as either the nutrient concentration or cellular motility is decreased. Starting from a mechanistic model of colony growth, we identify two dimensionless quantities that determine how morphology and genetic diversity of the population depend on the model parameters. Our simulations further reveal that population dynamics cannot be accurately described by the commonly-used surface growth models. Instead, one has to explicitly account for the emergent growth instabilities and demographic fluctuations. Overall, our work links together environmental conditions, colony morphology, and evolution. This link is essential for a rational design of concrete, biophysical perturbations to steer evolution in the desired direction.

Spatial population genetics with fluid flow

The growth and evolution of microbial populations is often subjected to advection by fluid flows in spatially extended environments, with immediate consequences for questions of spatial population genetics in marine ecology, planktonic diversity and origin of life scenarios. Here, we review recent progress made in understanding this rich problem in the simplified setting of two competing genetic microbial strains subjected to fluid flows. As a pedagogical example we focus on antagonsim, i.e., two killer microorganism strains, each secreting toxins that impede the growth of their competitors (competitive exclusion), in the presence of stationary fluid flows. By solving two coupled reaction-diffusion equations that include advection by simple steady cellular flows composed of characteristic flow motifs in two dimensions (2D), we show how local flow shear and compressibility effects can interact with selective advantage to have a dramatic influence on genetic competition and fixation in spatially distributed populations. We analyze several 1D and 2D flow geometries including sources, sinks, vortices and saddles, and show how simple analytical models of the dynamics of the genetic interface can be used to shed light on the nucleation, coexistence and flow-driven instabilities of genetic drops. By exploiting an analogy with phase separation with nonconserved order parameters, we uncover how thesegeneticdrops harness fluid flows for novel evolutionary strategies, even in the presence of number fluctuations, as confirmed by agent-based simulations as well.

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.

Branching structure of genealogies in spatially growing populations and its implications for population genetics inference

Spatial models where growth is limited to the population edge have been instrumental to understanding the population dynamics and the clone size distribution in growing cellular populations, such as microbial colonies and avascular tumours. A complete characterization of the coalescence process generated by spatial growth is still lacking, limiting our ability to apply classic population genetics inference to spatially growing populations. Here, we start filling this gap by investigating the statistical properties of the cell lineages generated by the two dimensional Eden model, leveraging their physical analogy with directed polymers. Our analysis provides quantitative estimates for population measurements that can easily be assessed via sequencing, such as the average number of segregating sites and the clone size distribution of a subsample of the population. Our results not only reveal remarkable features of the genealogies generated during growth, but also highlight new properties that can be misinterpreted as signs of selection if non-spatial models are inappropriately applied.

Population genetics in microchannels

Many microbial populations proliferate in small channels. In such environments, reproducing cells organize in parallel lanes. Reproducing cells shift these lanes, potentially expelling other cells from the channel. In this paper, we combine theory and experiments to understand how these dynamics affects the diversity of a microbial population. We theoretically predict that genetic diversity is quickly lost along lanes of cells. Our experiments confirm that a population of proliferating Escherichia coli in a microchannel organizes into lanes of genetically identical cells within a few generations. Our findings elucidate the effect of lane formation on populations evolution, with potential applications ranging from microbial ecology in soil to dynamics of epithelial tissues in higher organisms.

Spatial constraints, such as rigid barriers, affect the dynamics of cell populations, potentially altering the course of natural evolution. In this paper, we investigate the population genetics of Escherichia coli proliferating in microchannels with open ends. Our analysis is based on a population model, in which reproducing cells shift entire lanes of cells toward the open ends of the channel. The model predicts that diversity is lost very rapidly within lanes but at a much slower pace among lanes. As a consequence, two mixed, neutral E. coli strains competing in a microchannel must organize into an ordered regular stripe pattern in the course of a few generations. These predictions are in quantitative agreement with our experiments. We also demonstrate that random mutations appearing in the middle of the channel are much more likely to reach fixation than those occurring elsewhere. Our results illustrate fundamental mechanisms of microbial evolution in spatially confined space.

Survival in branching cellular populations

We analyze evolutionary dynamics in a confluent, branching cellular population, such as in a growing duct, vasculature, or in a branching microbial colony. We focus on the coarse-grained features of the evolution and build a statistical model that captures the essential features of the dynamics. Using simulations and analytic approaches, we show that the survival probability of strains within the growing population is sensitive to the branching geometry: Branch bifurcations enhance survival probability due to an overall population growth (i.e., "inflation"), while branch termination and the small effective population size at the growing branch tips increase the probability of strain extinction. We show that the evolutionary dynamics may be captured on a wide range of branch geometries parameterized just by the branch diameter N₀ and branching rate b. We find that the survival probability of neutral cell strains is largest at an "optimal" branching rate, which balances the effects of inflation and branch termination. We find that increasing the selective advantage s of the cell strain mitigates the inflationary effect by decreasing the average time at which the mutant cell fate is determined. For sufficiently large selective advantages, the survival probability of the advantageous mutant decreases monotonically with the branching rate.

Slow expanders invade by forming dented fronts in microbial colonies

Living organisms never cease to evolve, so there is a significant interest in predicting and controlling evolution in all branches of life sciences. The most basic question is whether a trait should increase or decrease in a given environment. The answer seems to be trivial for traits such as the growth rate in a bioreactor or the expansion rate of a tumor. Yet, it has been suggested that such traits can decrease, rather than increase, during evolution. Here, we report a mutant that outcompeted the ancestor despite having a slower expansion velocity when in isolation. To explain this observation, we developed and validated a theory that describes spatial competition between organisms with different expansion rates and arbitrary competitive interactions.

Most organisms grow in space, whether they are viruses spreading within a host tissue or invasive species colonizing a new continent. Evolution typically selects for higher expansion rates during spatial growth, but it has been suggested that slower expanders can take over under certain conditions. Here, we report an experimental observation of such population dynamics. We demonstrate that mutants that grow slower in isolation nevertheless win in competition, not only when the two types are intermixed, but also when they are spatially segregated into sectors. The latter was thought to be impossible because previous studies focused exclusively on the global competitions mediated by expansion velocities, but overlooked the local competitions at sector boundaries. Local competition, however, can enhance the velocity of either type at the sector boundary and thus alter expansion dynamics. We developed a theory that accounts for both local and global competitions and describes all possible sector shapes. In particular, the theory predicted that a slower on its own, but more competitive, mutant forms a dented V-shaped sector as it takes over the expansion front. Such sectors were indeed observed experimentally, and their shapes matched quantitatively with the theory. In simulations, we further explored several mechanisms that could provide slow expanders with a local competitive advantage and showed that they are all well-described by our theory. Taken together, our results shed light on previously unexplored outcomes of spatial competition and establish a universal framework to understand evolutionary and ecological dynamics in expanding populations.

Seascape origin of Richards growth

First proposed as an empirical rule over half a century ago, the Richards growth equation has been frequently invoked in population modeling and pandemic forecasting. Central to this model is the advent of a fractional exponent γ, typically fitted to the data. While various motivations for this nonanalytical form have been proposed, it is still considered foremost an empirical fitting procedure. Here, we find that Richards-like growth laws emerge naturally from generic analytical growth rules in a distributed population, upon inclusion of (i) migration (spatial diffusion) among different locales, and (ii) stochasticity in the growth rate, also known as "seascape noise." The latter leads to a wide (power law) distribution in local population number that, while smoothened through the former, can still result in a fractional growth law for the overall population. This justification of the Richards growth law thus provides a testable connection to the distribution of constituents of the population.

Snowdrift game induces pattern formation in systems of self-propelled particles

Evolutionary games between species are known to lead to intriguing spatiotemporal patterns in systems of diffusing agents. However, the role of interspecies interactions is hardly studied when agents are (self-)propelled, as is the case in many biological systems. Here, we combine aspects from active matter and evolutionary game theory and study a system of two species whose individuals are (self-)propelled and interact through a snowdrift game. We derive hydrodynamic equations for the density and velocity fields of both species from which we identify parameter regimes in which one or both species form macroscopic orientational order as well as regimes of propagating wave patterns. Interestingly, we find simultaneous wave patterns in both species that result from the interplay between alignment and snowdrift interactions-a feedback mechanism that we call game-induced pattern formation. We test these results in agent-based simulations and confirm the different regimes of order and spatiotemporal patterns as well as game-induced pattern formation.

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. This protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

Strong noise limit for population dynamics in incompressible advection

Genetic diversity is at the basis of the evolution process of populations and it is responsible for the populations' degree of fitness to a particular ecosystem. In marine environments many factors play a role in determining the dynamics of a population, including the amount of nutrients, the temperature, and many other stressing factors. An important and yet rather unexplored challenge is to figure out the role of individuals' dispersion, due to flow advection, on population genetics. In this paper we focus on two populations, one of which has a slight selective advantage, advanced by an incompressible two-dimensional flow. In particular, we want to understand how this advective flow can modify the dynamics of the advantageous allele. We generalize, through a theoretical analysis, previous evidence according to which the fixation probability is independent of diffusivity, showing that this is also independent of fluid advection. These findings may have important implications in the understanding of the dynamics of a population of microorganism, such as plankton or bacteria, in marine environments under the influence of (turbulent) currents.

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.

The large N c limit of QCD on the lattice

We review recent progress in the study of the large N c limit of gauge theories from lattice simulations. The focus is not only the planar limit but also the size of O ( N c - 1 ) corrections for values of N c ≳ 3 . Some concrete examples of the topics we include are tests of large- N c factorization, the topological susceptibility, the glueball, meson and baryon spectra, the chiral dependence of masses and decay constants, and weak matrix elements related to the Δ I = 1 / 2 rule in kaon decays.

Theoretical modelling of competitive microbial range expansion with heterogeneous mechanical interactions

Microbial range expansion experiments provide insight into the complex link between dynamic structure, pattern formation and evolutionary dynamics of growing populations. In this work, we develop a theoretical model in order to investigate the interplay of growth statistics and mechanical interactions which are implemented as division driven pushing and swapping of cells. For the case of the competitive growth of a strongly and a weakly interacting strain we investigate the influence of different mean division times, as well as different mechanical interactions on the development of the colony. Our results show that the susceptibility to cell division induced pushing has a much stronger influence on the structure of the colony than cell sorting towards the colony's perimeter. Motivated by microbial range expansion experiments of Neisseria gonorrhoeae bacteria, we also consider the influence of mutating cells on the structure of the colony. We show that the outgrowth of the three different strains is strongly influenced by the relative strengths of their mechanical interaction. The experimentally observed patterns are reproduced for mechanical interactions of the mutants, which range between those of the strongly and weakly interacting strain.

Geospatial distributions reflect temperatures of linguistic features

Quantifying the speed of linguistic change is challenging because the historical evolution of languages is sparsely documented. Consequently, traditional methods rely on phylogenetic reconstruction. Here, we propose a model-based approach to the problem through the analysis of language change as a stochastic process combining vertical descent, spatial interactions, and mutations in both dimensions. A notion of linguistic temperature emerges naturally from this analysis as a dimensionless measure of the propensity of a linguistic feature to undergo change. We demonstrate how temperatures of linguistic features can be inferred from their present-day geospatial distributions, without recourse to information about their phylogenies. Thus, the evolutionary dynamics of language, operating across thousands of years, leave a measurable geospatial signature. This signature licenses inferences about the historical evolution of languages even in the absence of longitudinal data.

Population extinction on a random fitness seascape

We explore the role of stochasticity and noise in the statistical outcomes of commonly studied population dynamics models within a space-independent (mean-field) perspective. Specifically, we consider a distributed population with logistic growth at each location, subject to "seascape" noise, wherein the population's fitness randomly varies with location and time. Despite its simplicity, the model actually incorporates variants of directed percolation, and directed polymers in random media, within a mean-field perspective. Probability distributions of the population can be computed self-consistently, and the extinction transition is shown to exhibit novel critical behavior with exponents dependent on the ratio of the strengths of migration and noise amplitudes. The results are compared and contrasted with the more conventional choice of demographic noise due to stochastic changes in reproduction.

Invasion front dynamics in disordered environments

Invasion occurs in environments that are normally spatially disordered, however, the effect of such a randomness on the dynamics of the invasion front has remained less understood. Here, we study Fisher's equation in disordered environments both analytically and numerically. Using the Effective Medium Approximation, we show that disorder slows down invasion velocity and for ensemble average of invasion velocity in disordered environment we have [Formula: see text] where [Formula: see text] is the amplitude of disorder and [Formula: see text] is the invasion velocity in the corresponding homogeneous environment given by [Formula: see text]. Additionally, disorder imposes fluctuations on the invasion front. Using a perturbative approach, we show that these fluctuations are Brownian with a diffusion constant of: [Formula: see text]. These findings were approved by numerical analysis. Alongside this continuum model, we use the Stepping Stone Model to check how our findings change when we move from the continuum approach to a discrete approach. Our analysis suggests that individual-based models exhibit inherent fluctuations and the effect of environmental disorder becomes apparent for large disorder intensity and/or high carrying capacities.

The impact of long-range dispersal on gene surfing

Range expansions lead to distinctive patterns of genetic variation in populations, even in the absence of selection. These patterns and their genetic consequences have been well studied for populations advancing through successive short-ranged migration events. However, most populations harbor some degree of long-range dispersal, experiencing rare yet consequential migration events over arbitrarily long distances. Although dispersal is known to strongly affect spatial genetic structure during range expansions, the resulting patterns and their impact on neutral diversity remain poorly understood. Here, we systematically study the consequences of long-range dispersal on patterns of neutral variation during range expansion in a class of dispersal models which spans the extremes of local (effectively short-ranged) and global (effectively well-mixed) migration. We find that sufficiently long-ranged dispersal leaves behind a mosaic of monoallelic patches, whose number and size are highly sensitive to the distribution of dispersal distances. We develop a coarse-grained model which connects statistical features of these spatial patterns to the evolution of neutral diversity during the range expansion. We show that growth mechanisms that appear qualitatively similar can engender vastly different outcomes for diversity: Depending on the tail of the dispersal distance distribution, diversity can be either preserved (i.e., many variants survive) or lost (i.e., one variant dominates) at long times. Our results highlight the impact of spatial and migratory structure on genetic variation during processes as varied as range expansions, species invasions, epidemics, and the spread of beneficial mutations in established populations.

Discrete Eulerian model for population genetics and dynamics under flow

Marine species reproduce and compete while being advected by turbulent flows. It is largely unknown, both theoretically and experimentally, how population dynamics and genetics are changed by the presence of fluid flows. Discrete agent-based simulations in continuous space allow for accurate treatment of advection and number fluctuations, but can be computationally expensive for even modest organism densities. In this report, we propose an algorithm to overcome some of these challenges. We first provide a thorough validation of the algorithm in one and two dimensions without flow. Next, we focus on the case of weakly compressible flows in two dimensions. This models organisms such as phytoplankton living at a specific depth in the three-dimensional, incompressible ocean experiencing upwelling and/or downwelling events. We show that organisms born at sources in a two-dimensional time-independent flow experience an increase in fixation probability.

Cooperation in Microbial Populations: Theory and Experimental Model Systems

Cooperative behavior, the costly provision of benefits to others, is common across all domains of life. This review article discusses cooperative behavior in the microbial world, mediated by the exchange of extracellular products called public goods. We focus on model species for which the production of a public good and the related growth disadvantage for the producing cells are well described. To unveil the biological and ecological factors promoting the emergence and stability of cooperative traits we take an interdisciplinary perspective and review insights gained from both mathematical models and well-controlled experimental model systems. Ecologically, we include crucial aspects of the microbial life cycle into our analysis and particularly consider population structures where ensembles of local communities (subpopulations) continuously emerge, grow, and disappear again. Biologically, we explicitly consider the synthesis and regulation of public good production. The discussion of the theoretical approaches includes general evolutionary concepts, population dynamics, and evolutionary game theory. As a specific but generic biological example, we consider populations of Pseudomonas putida and its regulation and use of pyoverdines, iron scavenging molecules, as public goods. The review closes with an overview on cooperation in spatially extended systems and also provides a critical assessment of the insights gained from the experimental and theoretical studies discussed. Current challenges and important new research opportunities are discussed, including the biochemical regulation of public goods, more realistic ecological scenarios resembling native environments, cell-to-cell signaling, and multispecies communities.

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.

Evolution in range expansions with competition at rough boundaries

When a biological population expands into new territory, genetic drift develops an enormous influence on evolution at the propagating front. In such range expansion processes, fluctuations in allele frequencies occur through stochastic spatial wandering of both genetic lineages and the boundaries between genetically segregated sectors. Laboratory experiments on microbial range expansions have shown that this stochastic wandering, transverse to the front, is superdiffusive due to the front's growing roughness, implying much faster loss of genetic diversity than predicted by simple flat front diffusive models. We study the evolutionary consequences of this superdiffusive wandering using two complementary numerical models of range expansions: the stepping stone model, and a new interpretation of the model of directed paths in random media, in the context of a roughening population front. Through these approaches we compute statistics for the times since common ancestry for pairs of individuals with a given spatial separation at the front, and we explore how environmental heterogeneities can locally suppress these superdiffusive fluctuations.

Nucleation of antagonistic organisms and cellular competitions on curved, inflating substrates

We consider the dynamics of spatially distributed, diffusing populations of organisms with antagonistic interactions. These interactions are found on many length scales, ranging from kilometer-scale animal range dynamics with selection against hybrids to micron-scale interactions between poison-secreting microbial populations. We find that the dynamical line tension at the interface between antagonistic organisms suppresses survival probabilities of small clonal clusters: the line tension introduces a critical cluster size that an organism with a selective advantage must achieve before deterministically spreading through the population. We calculate the survival probability as a function of selective advantage δ and antagonistic interaction strength σ. Unlike a simple Darwinian selective advantage, the survival probability depends strongly on the spatial diffusion constant D_{s} of the strains when σ>0, with suppressed survival when both species are more motile. Finally, we study the survival probability of a single mutant cell at the frontier of a growing spherical cluster of cells, such as the surface of an avascular spherical tumor. Both the inflation and curvature of the frontier significantly enhance the survival probability by changing the critical size of the nucleating cell cluster.

Lane formation and critical coarsening in a model of bacterial competition

We study competition of two nonmotile bacterial strains in a three-dimensional channel numerically and analyze how their configuration evolves in space and time. We construct a lattice model that takes into account self-replication, mutation, and killing of bacteria. When mutation is not significant, the two strains segregate and form stripe patterns along the channel. The formed lanes are gradually rearranged, with increasing length scales in the two-dimensional cross-sectional plane. We characterize it in terms of coarsening and phase ordering in statistical physics. In particular, for the simple model without mutation and killing, we find logarithmically slow coarsening, which is characteristic of the two-dimensional voter model. With mutation and killing, we find a phase transition from a monopolistic phase, in which lanes are formed and coarsened until the system is eventually dominated by one of the two strains, to an equally mixed and disordered phase without lane structure. Critical behavior at the transition point is also studied and compared with the generalized voter class and the Ising class. These results are accounted for by continuum equations, obtained by applying a mean-field approximation along the channel axis. Our findings indicate relevance of critical coarsening of two-dimensional systems in the problem of bacterial competition within anisotropic three-dimensional geometry.

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.

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 range expansions on a growing front: Roughness, fixation, and directed percolation

Directed percolation (DP) is a classic model for nonequilibrium phase transitions into a single absorbing state (fixation). It has been extensively studied by analytical and numerical techniques in diverse contexts. Recently, DP has appeared as a generic model for the evolutionary and ecological dynamics of competing bacterial populations. Range expansion-the stochastic reproduction of bacteria competing for space to be occupied by their progeny-leads to a fluctuating and rough growth front, which is known from experiment and simulation to affect the underlying critical behavior of the DP transition. In this work, we employ symmetry arguments to construct a pair of nonlinear stochastic partial differential equations describing the coevolution of surface roughness with the composition field of DP. Macroscopic manifestations (phenomenology) of these equations on growth patterns and genealogical tracks of range expansion are discussed; followed by a renormalization group analysis of possible scaling behaviors at the DP transition.

Emergence of evolutionary driving forces in pattern-forming microbial populations

Evolutionary dynamics are controlled by a number of driving forces, such as natural selection, random genetic drift and dispersal. In this perspective article, we aim to emphasize that these forces act at the population level, and that it is a challenge to understand how they emerge from the stochastic and deterministic behaviour of individual cells. Even the most basic steric interactions between neighbouring cells can couple evolutionary outcomes of otherwise unrelated individuals, thereby weakening natural selection and enhancing random genetic drift. Using microbial examples of varying degrees of complexity, we demonstrate how strongly cell-cell interactions influence evolutionary dynamics, especially in pattern-forming systems. As pattern formation itself is subject to evolution, we propose to study the feedback between pattern formation and evolutionary dynamics, which could be key to predicting and potentially steering evolutionary processes. Such an effort requires extending the systems biology approach from the cellular to the population scale.This article is part of the theme issue 'Self-organization in cell biology'.

Available energy fluxes drive a transition in the diversity, stability, and functional structure of microbial communities

A fundamental goal of microbial ecology is to understand what determines the diversity, stability, and structure of microbial ecosystems. The microbial context poses special conceptual challenges because of the strong mutual influences between the microbes and their chemical environment through the consumption and production of metabolites. By analyzing a generalized consumer resource model that explicitly includes cross-feeding, stochastic colonization, and thermodynamics, we show that complex microbial communities generically exhibit a transition as a function of available energy fluxes from a “resource-limited” regime where community structure and stability is shaped by energetic and metabolic considerations to a diverse regime where the dominant force shaping microbial communities is the overlap between species’ consumption preferences. These two regimes have distinct species abundance patterns, different functional profiles, and respond differently to environmental perturbations. Our model reproduces large-scale ecological patterns observed across multiple experimental settings such as nestedness and differential beta diversity patterns along energy gradients. We discuss the experimental implications of our results and possible connections with disorder-induced phase transitions in statistical physics.

The diversity, stability and functional structure of microbial communities have dramatic effects on the health of humans and of ecosystems. The complexity of these communities has so far precluded the development of a general predictive model that would capture the dependence of these features on environmental conditions. We confronted this challenge by constructing a flexible simulation framework, and randomly sampling parameters under a variety of modeling assumptions to identify generic patterns. We found two qualitatively distinct regimes of community structure, which reproduce observed patterns of biodiversity, and make new predictions about stability and function.

Reaction fronts in persistent random walks with demographic stochasticity

Standard reaction-diffusion systems are characterized by infinite velocities and no persistence in the movement of individuals, two conditions that are violated when considering living organisms. Here we consider a discrete particle model in which individuals move following a persistent random walk with finite speed and grow with logistic dynamics. We show that, when the number of individuals is very large, the individual-based model is well described by the continuous reactive Cattaneo equation (RCE), but for smaller values of the carrying capacity important finite-population effects arise. The effects of fluctuations on the propagation speed are investigated both considering the RCE with a cutoff in the reaction term and by means of numerical simulations of the individual-based model. Finally, a more general Lévy walk process for the transport of individuals is examined and an expression for the front speed of the resulting traveling wave is proposed.

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.

Collective motion conceals fitness differences in crowded cellular populations

Many cellular populations are tightly-packed, such as microbial colonies and biofilms, or tissues and tumors in multicellular organisms. Movement of one cell in those crowded assemblages requires motion of others, so that cell displacements are correlated over many cell diameters. Whenever movement is important for survival or growth, these correlated rearrangements could couple the evolutionary fate of different lineages. Yet, little is known about the interplay between mechanical forces and evolution in dense cellular populations. Here, by tracking slower-growing clones at the expanding edge of yeast colonies, we show that the collective motion of cells prevents costly mutations from being weeded out rapidly. Joint pushing by neighboring cells generates correlated movements that suppress the differential displacements required for selection to act. This mechanical screening of fitness differences allows slower-growing mutants to leave more descendants than expected under non-mechanical models, thereby increasing their chance for evolutionary rescue. Our work suggests that, in crowded populations, cells cooperate with surrounding neighbors through inevitable mechanical interactions. This effect has to be considered when predicting evolutionary outcomes, such as the emergence of drug resistance or cancer evolution.

Physical interactions reduce the power of natural selection in growing yeast colonies

Microbial populations often assemble in dense populations in which proliferating individuals exert mechanical forces on the nearby cells. Here, we use yeast strains whose doubling times depend differently on temperature to show that physical interactions among cells affect the competition between different genotypes in growing yeast colonies. Our experiments demonstrate that these physical interactions have two related effects: they cause the prolonged survival of slower-growing strains at the actively-growing frontier of the colony and cause faster-growing strains to increase their frequency more slowly than expected in the absence of physical interactions. These effects also promote the survival of slower-growing strains and the maintenance of genetic diversity in colonies grown in time-varying environments. A continuum model inspired by overdamped hydrodynamics reproduces the experiments and predicts that the strength of natural selection depends on the width of the actively growing layer at the colony frontier. We verify these predictions experimentally. The reduced power of natural selection observed here may favor the maintenance of drug-resistant cells in microbial populations and could explain the apparent neutrality of interclone competition within tumors.

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.

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.