Deleterious mutations can surf to high densities on the wave front of an expanding population

There is an increasing recognition that evolutionary processes play a key role in determining the dynamics of range expansion. Recent work demonstrates that neutral mutations arising near the edge of a range expansion sometimes surf on the expanding front leading them rather than that leads to reach much greater spatial distribution and frequency than expected in stationary populations. Here, we extend this work and examine the surfing behavior of nonneutral mutations. Using an individual-based coupled-map lattice model, we confirm that, regardless of its fitness effects, the probability of survival of a new mutation depends strongly upon where it arises in relation to the expanding wave front. We demonstrate that the surfing effect can lead to deleterious mutations reaching high densities at an expanding front, even when they have substantial negative effects on fitness. Additionally, we highlight that this surfing phenomenon can occur for mutations that impact reproductive rate (i.e., number of offspring produced) as well as mutations that modify juvenile competitive ability. We suggest that these effects are likely to have important consequences for rates of spread and the evolution of spatially expanding populations.     

The frequency of fitness peak shifts is increased at expanding range margins due to mutation surfing

Dynamic species' ranges, those that are either invasive or shifting in response to environmental change, are the focus of much recent interest in ecology, evolution, and genetics. Understanding how range expansions can shape evolutionary trajectories requires the consideration of nonneutral variability and genetic architecture, yet the majority of empirical and theoretical work to date has explored patterns of neutral variability. Here we use forward computer simulations of population growth, dispersal, and mutation to explore how range-shifting dynamics can influence evolution on rugged fitness landscapes. We employ a two-locus model, incorporating sign epistasis, and find that there is an increased likelihood of fitness peak shifts during a period of range expansion. Maladapted valley genotypes can accumulate at an expanding range front through a phenomenon called mutation surfing, which increases the likelihood that a mutation leading to a higher peak will occur. Our results indicate that most peak shifts occur close to the expanding front. We also demonstrate that periods of range shifting are especially important for peak shifting in species with narrow geographic distributions. Our results imply that trajectories on rugged fitness landscapes can be modified substantially when ranges are dynamic.     

Long distance dispersal and the fate of a gene from the colonization front

There is an increasing recognition that long distance dispersal (LDD) plays a key role in establishing spatial genetic structure during colonization. Recent works, focused on short distance dispersal, demonstrated that a neutral mutation arising at the colonization front can either ‘surf’ with the wave front and reach high frequencies or stay near its place of origin at low frequencies. Here, we examine how LDD, and more generally the shape of the dispersal kernel, modifies this phenomenon and how the width of the colonization corridor affects the fate of the mutation. We demonstrate that when LDD events are more frequent, the ‘surfing phenomenon’ is less frequent, probably because any alleles can get far ahead from the colonization front and preclude the invasion by others alleles, thus leading to an attenuation of the diversity loss. We also demonstrate that the width of the colonization corridor influences the fate of the mutation, wide spaces decreasing the probability of invasion. Overall, the genetic structure of diversity resulted not only from LDD but also particularly from the shape of the dispersal kernel.     

The rate of beneficial mutations surfing on the wave of a range expansion

Many theoretical and experimental studies suggest that range expansions can have severe consequences for the gene pool of the expanding population. Due to strongly enhanced genetic drift at the advancing frontier, neutral and weakly deleterious mutations can reach large frequencies in the newly colonized regions, as if they were surfing the front of the range expansion. These findings raise the question of how frequently beneficial mutations successfully surf at shifting range margins, thereby promoting adaptation towards a range-expansion phenotype. Here, we use individual-based simulations to study the surfing statistics of recurrent beneficial mutations on wave-like range expansions in linear habitats. We show that the rate of surfing depends on two strongly antagonistic factors, the probability of surfing given the spatial location of a novel mutation and the rate of occurrence of mutations at that location. The surfing probability strongly increases towards the tip of the wave. Novel mutations are unlikely to surf unless they enjoy a spatial head start compared to the bulk of the population. The needed head start is shown to be proportional to the inverse fitness of the mutant type, and only weakly dependent on the carrying capacity. The precise location dependence of surfing probabilities is derived from the non-extinction probability of a branching process within a moving field of growth rates. The second factor is the mutation occurrence which strongly decreases towards the tip of the wave. Thus, most successful mutations arise at an intermediate position in the front of the wave. We present an analytic theory for the tradeoff between these factors that allows to predict how frequently substitutions by beneficial mutations occur at invasion fronts. We find that small amounts of genetic drift increase the fixation rate of beneficial mutations at the advancing front, and thus could be important for adaptation during species invasions.     

Landscape structure and boundary effects determine the fate of mutations occurring during range expansions

The interplay between the spatial dynamics of range expansion and evolutionary processes is receiving considerable attention. Recent theory has demonstrated that mutations occurring towards the front of a spatially expanding population can sometimes 'surf' to high frequency and spatial extent. Here, we extend this work to consider how the fate of a novel mutation is influenced by where and when it occurs. Specifically, we are interested in establishing how the origin of a mutation relative to a habitat edge influences its dynamics, and in understanding how this is mediated by the behaviour of individuals at those boundaries. Using a coupled-map lattice model, we demonstrate that the survival probability, abundance and spatial extent of surviving mutants can depend on their origin. An edge effect is often observed and can be quite different both qualitatively and quantitatively depending on the behavioural rules assumed. Mutations, especially those that are deleterious, that arise at a habitat edge with reflective boundary conditions can be many more times likely to survive for substantial periods of time than those that arise away from the edge. Conversely, with absorbing boundary conditions, their survival is greater when they arise well away from the edge. Our results clearly illustrate that landscape structure, habitat edges and boundary conditions have a considerable influence on the likely fate of mutations that occur during a period of range expansion.     

Accelerating invasion rates result from the evolution of density-dependent dispersal

Evolutionary processes play an important role in shaping the dynamics of range expansions, and selection on dispersal propensity has been demonstrated to accelerate rates of advance. Previous theory has considered only the evolution of unconditional dispersal rates, but dispersal is often more complex. For example, many species emigrate in response to crowding. Here, we use an individual-based model to investigate the evolution of density dependent dispersal into empty habitat, such as during an invasion. The landscape is represented as a lattice and dispersal between populations follows a stepping-stone pattern. Individuals carry three ‘genes’ that determine their dispersal strategy when experiencing different population densities. For a stationary range we obtain results consistent with previous theoretical studies: few individuals emigrate from patches that are below equilibrium density. However, during the range expansion of a previously stationary population, we observe evolution towards dispersal strategies where considerable emigration occurs well below equilibrium density. This is true even for moderate costs to dispersal, and always results in accelerating rates of range expansion. Importantly, the evolution we observe at an expanding front depends upon fitness integrated over several generations and cannot be predicted by a consideration of lifetime reproductive success alone. We argue that a better understanding of the role of density dependent dispersal, and its evolution, in driving population dynamics is required especially within the context of range expansions.     

Risky movement increases the rate of range expansion

The movement rules used by an individual determine both its survival and dispersal success. Here, we develop a simple model that links inter-patch movement behaviour with population dynamics in order to explore how individual dispersal behaviour influences not only its dispersal and survival, but also the population's rate of range expansion. Whereas dispersers are most likely to survive when they follow nearly straight lines and rapidly orient movement towards a non-natal patch, the most rapid rates of range expansion are obtained for trajectories in which individuals delay biasing their movement towards a non-natal patch. This result is robust to the spatial structure of the landscape. Importantly, in a set of evolutionary simulations, we also demonstrate that the movement strategy that evolves at an expanding front is much closer to that maximizing the rate of range expansion than that which maximizes the survival of dispersers. Our results suggest that if one of our conservation goals is the facilitation of range-shifting, then current indices of connectivity need to be complemented by the development and utilization of new indices providing a measure of the ease with which a species spreads across a landscape.     

The fate of mutations surfing on the wave of a range expansion

Many species, including humans, have dramatically expanded their range in the past, and such range expansions had certainly an impact on their genetic diversity. For example, mutations arising in populations at the edge of a range expansion can sometimes surf on the wave of advance and thus reach a larger spatial distribution and a much higher frequency than would be expected in stationary populations. We study here this surfing phenomenon in more detail, by performing extensive computer simulations under a two-dimensional stepping-stone model. We find that the probability of survival of a new mutation depends to a large degree on its proximity to the edge of the wave. Demographic factors such as deme size, migration rate, and local growth rate also influence the fate of these new mutations. We also find that the final spatial and frequency distributions depend on the local deme size of a subdivided population. This latter result is discussed in the light of human expansions in Europe as it should allow one to distinguish between mutations having spread with Paleolithic or Neolithic expansions. By favoring the spread of new mutations, a consequence of the surfing phenomenon is to increase the rate of evolution of spatially expanding populations.     

On the evolution of dispersal via heterogeneity in spatial connectivity

Dispersal has long been recognized as a mechanism that shapes many observed ecological and evolutionary processes. Thus, understanding the factors that promote its evolution remains a major goal in evolutionary ecology. Landscape connectivity may mediate the trade-off between the forces in favour of dispersal propensity (e.g. kin-competition, local extinction probability) and those against it (e.g. energetic or survival costs of dispersal). It remains, however, an open question how differing degrees of landscape connectivity may select for different dispersal strategies. We implemented an individual-based model to study the evolution of dispersal on landscapes that differed in the variance of connectivity across patches ranging from networks with all patches equally connected to highly heterogeneous networks. The parthenogenetic individuals dispersed based on a flexible logistic function of local abundance. Our results suggest, all else being equal, that landscapes differing in their connectivity patterns will select for different dispersal strategies and that these strategies confer a long-term fitness advantage to individuals at the regional scale. The strength of the selection will, however, vary across network types, being stronger on heterogeneous landscapes compared with the ones where all patches have equal connectivity. Our findings highlight how landscape connectivity can determine the evolution of dispersal strategies, which in turn affects how we think about important ecological dynamics such as metapopulation persistence and range expansion.     

Survival of GacS/GacA mutants of the biological control bacterium Pseudomonas aureofaciens 30-84 in the wheat rhizosphere

GacS/GacA comprises a two-component regulatory system that controls the expression of secondary metabolites required for the control of plant diseases in many pseudomonads. High mutation frequencies of gacS and gacA have been observed in liquid culture. We examined whether gacS/gacA mutants could competitively displace the wild-type populations on roots and thus pose a threat to the efficacy of biological control. The survival of a gac mutant alone and in competition with the wild type on roots was examined in the biological control strain Pseudomonas aureofaciens 30-84. In this bacterium, GacS/GacA controls the expression of phenazine antibiotics that are inhibitory to plant pathogenic fungi and enhance the competitive survival of the bacterium. Wheat seedlings were inoculated with strain 30-84, and bacteria were recovered from roots after 21 days in sterile or nonsterile soil to check for the presence of gacS or gacA mutants. Although no mutants were detected in the inoculum, gacS/gacA mutants were recovered from 29 out of 31 roots and comprised up to 36% of the total bacterial populations. Southern hybridization analysis of the recovered gacA mutants did not indicate a conserved mutational mechanism. Replacement series analysis on roots utilizing strain 30-84 and a gacA mutant (30-84.gacA) or a gacS mutant (30-84.A2) demonstrated that although the mutant population partially displaced the wild type in sterile soil, it did not do so in natural soil. In fact, in natural soil final rhizosphere populations of wild-type strain 30-84 starting from mixtures were at least 1.5 times larger than would be predicted from their inoculation ratio and generally were greater than or equal to the population of wild type alone despite lower inoculation rates. These results indicate that although gacS/gacA mutants survive in natural rhizosphere populations, they do not displace wild-type populations. Better survival of wild-type populations in mixtures with mutants suggests that mutants arising de novo or introduced within the inoculum may be beneficial for the survival of wild-type populations in the rhizosphere.     

Reid's paradox revisited: the evolution of dispersal kernels during range expansion

Current approaches to modeling range advance assume that the distribution describing dispersal distances in the population (the "dispersal kernel") is a static entity. We argue here that dispersal kernels are in fact highly dynamic during periods of range advance because density effects and spatial assortment by dispersal ability ("spatial selection") drive the evolution of increased dispersal on the expanding front. Using a spatially explicit individual-based model, we demonstrate this effect under a wide variety of population growth rates and dispersal costs. We then test the possibility of an evolved shift in dispersal kernels by measuring dispersal rates in individual cane toads (Bufo marinus) from invasive populations in Australia (historically, toads advanced their range at 10 km/year, but now they achieve >55 km/year in the northern part of their range). Under a common-garden design, we found a steady increase in dispersal tendency with distance from the invasion origin. Dispersal kernels on the invading front were less kurtotic and less skewed than those from origin populations. Thus, toads have increased their rate of range expansion partly through increased dispersal on the expanding front. For accurate long-range forecasts of range advance, we need to take into account the potential for dispersal kernels to be evolutionarily dynamic.     

Survival, Growth, and Localization of Epiphytic Fitness Mutants of Pseudomonas syringae on Leaves

Among 82 epiphytic fitness mutants of a Pseudomonas syringae pv. syringae strain that were characterized in a previous study, 4 mutants were particularly intolerant of the stresses associated with dry leaf surfaces. These four mutants each exhibited distinctive behaviors when inoculated onto and into plant leaves. For example, while none showed measurable growth on dry potato leaf surfaces, they grew to different population sizes in the intercellular spaces of bean leaves and on dry bean leaf surfaces, and one mutant appeared incapable of growth in both environments although it grew well on moist bean leaves. The presence of the parental strain did not influence the survival of the mutants immediately following exposure of leaves to dry, high-light incubation conditions, suggesting that the reduced survival of the mutants did not result from an inability to produce extracellular factors in planta. On moist bean leaves that were colonized by either a mutant or the wild type, the proportion of the total epiphytic population that was located in sites protected from a surface sterilant was smaller for the mutants than for the wild type, indicating that the mutants were reduced in their ability to locate, multiply in, and/or survive in such protected sites. This reduced ability was only one of possibly several factors contributing to the reduced epiphytic fitness of each mutant. Their reduced fitness was not specific to the host plant bean, since they also exhibited reduced fitness on the nonhost plant potato; the functions altered in these strains are thus of interest for their contribution to the general fitness of bacterial epiphytes.     

Neutral processes forming large clones during colonization of new areas

In species reproducing both sexually and asexually clones are often more common in recently established populations. Earlier studies have suggested that this pattern arises due to natural selection favouring generally or locally successful genotypes in new environments. Alternatively, as we show here, this pattern may result from neutral processes during species’ range expansions. We model a dioecious species expanding into a new area in which all individuals are capable of both sexual and asexual reproduction, and all individuals have equal survival rates and dispersal distances. Even under conditions that favour sexual recruitment in the long run, colonization starts with an asexual wave. After colonization is completed, a sexual wave erodes clonal dominance. If individuals reproduce more than one season, and with only local dispersal, a few large clones typically dominate for thousands of reproductive seasons. Adding occasional long‐distance dispersal, more dominant clones emerge, but they persist for a shorter period of time. The general mechanism involved is simple: edge effects at the expansion front favour asexual (uniparental) recruitment where potential mates are rare. Specifically, our model shows that neutral processes (with respect to genotype fitness) during the population expansion, such as random dispersal and demographic stochasticity, produce genotype patterns that differ from the patterns arising in a selection model. The comparison with empirical data from a post‐glacially established seaweed species (Fucus radicans) shows that in this case, a neutral mechanism is strongly supported.     

Reduced dispersal propensity in the wingless waterstrider Aquarius najas in a highly fragmented landscape

Dispersal behaviour of animals in fragmented habitats has generated intensive theoretical attention but empirical data on the evolution of dispersal are still relatively scarce. Theory predicts reduced dispersal propensity in small and isolated habitat patches. We tested these predictions in the waterstrider Aquarius najas, a wingless species with special habitat demands. Flightlessness constrains insect dispersal and as a stream specialist A. najas cannot survive on still water. Lakes therefore represent a dispersal barrier for this species. We measured dispersal propensity of the waterstrider A. najas which originated from ten fragmented populations. In the experiment, we transplanted laboratory-grown individuals to the field. We did not find differences between sex in dispersal propensity. However, we found that waterstriders that originated from small and isolated patches moved less than individuals from large and more continuous habitats. This suggests that the cost of dispersal over hostile surrounding habitats may be high. We conclude that a low likelihood of dispersal in A. najas is an adaptation to small and isolated stream habitats.     

LIFE AT THE FRONT OF AN EXPANDING POPULATION

Environmental changes have caused episodes of habitat expansions in the evolutionary history of many species. These range changes affect the dynamics of biological evolution in multiple ways. Recent microbial experiments as well as simulations suggest that enhanced genetic drift at the frontier of a two-dimensional range expansion can cause genetic sectoring patterns with fractal domain boundaries. Here, we propose and analyze a simple model of asexual biological evolution at expanding frontiers that explains these neutral patterns and predicts the effect of natural selection. We find that beneficial mutations give rise to sectors with an opening angle that depends sensitively on the selective advantage of the mutants. Deleterious mutations, on the other hand, are not able to establish a sector permanently. They can, however, temporarily "surf" on the population front, and thereby reach unusually high frequencies. As a consequence, expanding frontiers are loaded with a high fraction of mutants at mutation–selection balance. Numerically, we also determine the condition at which the wild type is lost in favor of deleterious mutants (genetic meltdown) at a growing front. Our prediction for this error threshold differs qualitatively from existing well-mixed theories, and sets tight constraints on sustainable mutation rates for populations that undergo frequent range expansions.     

Evolutionary cheating in Escherichia coli stationary phase cultures

Starved cultures of Escherichia coli are highly dynamic, undergoing frequent population shifts. The shifts result from the spread of mutants able to grow under conditions that impose growth arrest on the ancestral population. To analyze competitive interactions underlying this dynamic we measured the survival of a typical mutant and the wild type during such population shifts. Here we show that the survival advantage of the mutant at any given time during a takeover is inversely dependent on its frequency in the population, its growth adversely affects the survival of the wild type, and its ability to survive in stationary phase at fixation is lower than that of its ancestor. These mutants do not enter, or exit early, the nondividing stationary-phase state, cooperatively maintained by the wild type. Thus they end up overrepresented as compared to their initial frequency at the onset of the stationary phase, and subsequently they increase disproportionately their contribution in terms of progeny to the succeeding generation in the next growth cycle, which is a case of evolutionary cheating. If analyzed through the game theory framework, these results might be explained by the prisoner's dilemma type of conflict, which predicts that selfish defection is favored over cooperation.     

Ascorbate-deficient mutants of Arabidopsis grow in high light despite chronic photooxidative stress

Acclimation to changing environments, such as increases in light intensity, is necessary, especially for the survival of sedentary organisms like plants. To learn more about the importance of ascorbate in the acclimation of plants to high light (HL), vtc2, an ascorbate-deficient mutant of Arabidopsis, and the double mutants vtc2npq4 and vtc2npq1 were tested for growth in low light and HL and compared with the wild type. The vtc2 mutant has only 10% to 30% of wild-type levels of ascorbate, vtc2npq4 has lower ascorbate levels and lacks non-photochemical quenching of chlorophyll fluorescence (NPQ) because of the absence of the photosystem II protein PsbS, and vtc2npq1 is NPQ deficient and also lacks zeaxanthin in HL but has PsbS. All three genotypes were able to grow in HL and had wild-type levels of Lhcb1, cytochrome f, PsaF, and 2-cysteine peroxiredoxin. However, the mutants had lower electron transport and oxygen evolution rates and lower quantum efficiency of PSII compared with the wild type, implying that they experienced chronic photooxidative stress. The mutants lacking NPQ in addition to ascorbate were only slightly more affected than vtc2. All three mutants had higher glutathione levels than the wild type in HL, suggesting a possible compensation for the lower ascorbate content. These results demonstrate the importance of ascorbate for the long-term acclimation of plants to HL.     

Cell lineage dynamics in stratified shoot apical meristems

We devise a stochastic and spatially explicit model for the dynamics of the initials cells in a stratified shoot apical meristem (SAM). The meristem is composed of three layers with seven initials per layer. We investigate the probability and number of divisions for a mutant lineage to either reach fixation or becoming purged through selection or drift. In contrast to previous studies our results show that the functional organization of the initials in stratified SAMs acts as an efficient purging mechanism particularly of deleterious mutations. All mutants are rapidly purged when deleterious. The probability of fixation for mutants with a higher fitness than the wild type increases linearly up to 70%. The median number of divisions to fixation of both genotypes is insensitive to the mutant's fitness. The median number of divisions to wild-type fixation is less than 100, with the upper quartile below 200. The largest number of divisions to wild-type fixation are in the order of 100 000 divisions. Our results indicate that the spatial organization of SAM enables the efficient purging of mutant lineages, particularly if they are deleterious. On the other hand, long-lived chimeric stages are common when mutant lineages succeed to overcome the initial numerical disadvantage.     

Survival of GacS/GacA Mutants of the Biological Control Bacterium Pseudomonas aureofaciens 30-84 in the Wheat Rhizosphere

GacS/GacA comprises a two-component regulatory system that controls the expression of secondary metabolites required for the control of plant diseases in many pseudomonads. High mutation frequencies of gacS and gacA have been observed in liquid culture. We examined whether gacS/gacA mutants could competitively displace the wild-type populations on roots and thus pose a threat to the efficacy of biological control. The survival of a gac mutant alone and in competition with the wild type on roots was examined in the biological control strain Pseudomonas aureofaciens 30-84. In this bacterium, GacS/GacA controls the expression of phenazine antibiotics that are inhibitory to plant pathogenic fungi and enhance the competitive survival of the bacterium. Wheat seedlings were inoculated with strain 30-84, and bacteria were recovered from roots after 21 days in sterile or nonsterile soil to check for the presence of gacS or gacA mutants. Although no mutants were detected in the inoculum, gacS/gacA mutants were recovered from 29 out of 31 roots and comprised up to 36% of the total bacterial populations. Southern hybridization analysis of the recovered gacA mutants did not indicate a conserved mutational mechanism. Replacement series analysis on roots utilizing strain 30-84 and a gacA mutant (30-84.gacA) or a gacS mutant (30-84.A2) demonstrated that although the mutant population partially displaced the wild type in sterile soil, it did not do so in natural soil. In fact, in natural soil final rhizosphere populations of wild-type strain 30-84 starting from mixtures were at least 1.5 times larger than would be predicted from their inoculation ratio and generally were greater than or equal to the population of wild type alone despite lower inoculation rates. These results indicate that although gacS/gacA mutants survive in natural rhizosphere populations, they do not displace wild-type populations. Better survival of wild-type populations in mixtures with mutants suggests that mutants arising de novo or introduced within the inoculum may be beneficial for the survival of wild-type populations in the rhizosphere.     

Genetic admixture between captive-bred and wild individuals affects patterns of dispersal in a brown trout (<Emphasis Type="Italic">Salmo trutta</Emphasis>) population

Genetic admixture between captive-bred and wild individuals has been demonstrated to affect many individual traits, although little is known about its potential influence on dispersal, an important trait governing the eco-evolutionary dynamics of populations. Here, we quantified and described the spatial distribution of genetic admixture in a brown trout (Salmo trutta) population from a small watershed that was stocked until 1999, and then tested whether or not individual dispersal parameters were related to admixture between wild and captive-bred fish. We genotyped 715 fish at 17 microsatellite loci sampled from both the mainstream and all populated tributaries, as well as 48 fish from the hatchery used to stock the study area. First, we used Bayesian clustering to infer local genetic structure and to quantify genetic admixture. We inferred first generation migrants to identify dispersal events and test which features (genetic admixture, sex and body length) affected dispersal parameters (i.e. probability to disperse, distance of dispersal and direction of the dispersal event). We identified two genetic clusters in the river basin, corresponding to wild fish on the one hand and to fish derived from the captive strain on the other hand, allowing us to define an individual gradient of admixture. Individuals with a strong assignment to the captive strain occurred almost exclusively in some tributaries, and were more likely to disperse towards a tributary than towards a site of the mainstream. Furthermore, dispersal probability increased as the probability of assignment to the captive strain increased, and individuals with an intermediate level of admixture exhibited the lowest dispersal distances. These findings show that various dispersal parameters may be biased by admixture with captive-bred genotypes, and that management policies should take into account the differential spread of captive-bred individuals in wild populations.