Metapopulation-level adaptation of insect host plant preference and extinction-colonization dynamics in heterogeneous landscapes

Species living in highly fragmented landscapes typically occur as metapopulations with frequent turnover of local populations. The turnover rate depends on population sizes and connectivities, but it may also depend on the phenotypic and genotypic composition of populations. The Glanville fritillary butterfly (Melitaea cinxia) in Finland uses two host plant species, which show variation in their relative abundances at two spatial scales: locally among individual habitat patches and regionally among networks of patches. Female butterflies in turn exhibit spatial variation in genetically determined host plant preference within and among patch networks. Emigration, immigration and establishment of new populations have all been shown to be strongly influenced by the match between the host plant composition of otherwise suitable habitat patches and the host plant preference of migrating butterflies. The evolutionary consequences of such biased migration and colonization with respect to butterfly phenotypes might differ depending on spatial configuration and plant species composition of the patches in heterogeneous patch networks. Using a spatially realistic individual-based model we show that the model-predicted evolution of host plant preference due to biased migration explains a significant amount of the observed variation in host plant use among metapopulations living in dissimilar networks. This example illustrates how the ecological extinction-colonization dynamics may be linked with the evolutionary dynamics of life history traits in metapopulations.     

SPATIALLY CORRELATED EXTINCTIONS SELECT FOR LESS EMIGRATION BUT LARGER DISPERSAL DISTANCES IN THE SPIDER MITE TETRANYCHUS URTICAE

Dispersal is a central process to almost all species on earth, as it connects spatially structured populations and thereby increases population persistence. Dispersal is subject to (rapid) evolution and local patch extinctions are an important selective force in this context. In contrast to the randomly distributed local extinctions considered in most theoretical studies, habitat fragmentation or other anthropogenic interventions will lead to spatially correlated extinction patterns. Under such conditions natural selection is thought to lead to more long‐distance dispersal, but this theoretical prediction has not yet been verified empirically. We test this prediction in experimental spatially structured populations of the spider mite Tetranychus urticae and supplement these empirical results with insights from an individual‐based evolutionary model. We demonstrate that the spatial correlation of local extinctions changes the entire distribution of dispersal distances (dispersal kernel) and selects for overall less emigration but more long‐distance dispersal.     

Eco-evolutionary dynamics of dispersal in spatially heterogeneous environments

Ecology Letters (2011) 14: 1025-1034 ABSTRACT: Evolutionary changes in natural populations are often so fast that the evolutionary dynamics may influence ecological population dynamics and vice versa. Here we construct an eco-evolutionary model for dispersal by combining a stochastic patch occupancy metapopulation model with a model for changes in the frequency of fast-dispersing individuals in local populations. We test the model using data on allelic variation in the gene phosphoglucose isomerase (Pgi), which is strongly associated with dispersal rate in the Glanville fritillary butterfly. Population-specific measures of immigration and extinction rates and the frequency of fast-dispersing individuals among the immigrants explained 40% of spatial variation in Pgi allele frequency among 97 local populations. The model clarifies the roles of founder events and gene flow in dispersal evolution and resolves a controversy in the literature about the consequences of habitat loss and fragmentation on the evolution of dispersal.     

Finite metapopulation models with density-dependent migration and stochastic local dynamics

The effects of small density-dependent migration on the dynamics of a metapopulation are studied in a model with stochastic local dynamics. We use a diffusion approximation to study how changes in the migration rate and habitat occupancy affect the rates of local colonization and extinction. If the emigration rate increases or if the immigration rate decreases with local population size, a positive expected rate of change in habitat occupancy is found for a greater range of habitat occupancies than when the migration is density-independent. In contrast, the reverse patterns of density dependence in respective emigration and immigration reduce the range of habitat occupancies where the metapopulation will be viable. This occurs because density-dependent migration strongly influences both the establishment and rescue effects in the local dynamics of metapopulations.     

Matrix quality and disturbance frequency drive evolution of species behavior at habitat boundaries

Previous theoretical studies suggest that a species' landscape should influence the evolution of its dispersal characteristics, because landscape structure affects the costs and benefits of dispersal. However, these studies have not considered the evolution of boundary crossing, that is, the tendency of animals to cross from habitat to nonhabitat (“matrix”). It is important to understand this dispersal behavior, because of its effects on the probability of population persistence. Boundary‐crossing behavior drives the rate of interaction with matrix, and thus, it influences the rate of movement among populations and the risk of dispersal mortality. We used an individual‐based, spatially explicit model to simulate the evolution of boundary crossing in response to landscape structure. Our simulations predict higher evolved probabilities of boundary crossing in landscapes with more habitat, less fragmented habitat, higher‐quality matrix, and more frequent disturbances (i.e., fewer generations between local population extinction events). Unexpectedly, our simulations also suggest that matrix quality and disturbance frequency have much stronger effects on the evolution of boundary crossing than either habitat amount or habitat fragmentation. Our results suggest that boundary‐crossing responses are most affected by the costs of dispersal through matrix and the benefits of escaping local extinction events. Evolution of optimal behavior at habitat boundaries in response to the landscape may have implications for species in human‐altered landscapes, because this behavior may become suboptimal if the landscape changes faster than the species' evolutionary response to that change. Understanding how matrix quality and habitat disturbance drive evolution of behavior at boundaries, and how this in turn influences the extinction risk of species in human‐altered landscapes should help us identify species of conservation concern and target them for management.     

Eco-evolutionary metapopulation dynamics and the spatial scale of adaptation

We construct a model that combines extinction-colonization dynamics with the dynamics of local adaptation in a network of habitat patches of dissimilar qualities. We derive a deterministic approximation for the stochastic model that allows the calculation of patch-specific incidences of occupancy and levels of adaptation at steady state. Depending on (i) the strength of local selection, (ii) the amount of genetic variance, (iii) the demographic cost of maladaptation, (iv) the spatial scale of gene flow, and (v) the amount of habitat heterogeneity, the model predicts adaptation at different spatial scales. Local adaptation is predicted when there is much genetic variance and strong selection, while network-level adaptation occurs when the demographic cost of maladaptation is low. For little genetic variance and high cost of maladaptation, the model predicts network-level habitat specialization in species with long-range migration but an intermediate scale of adaptation (mosaic specialization) in species with short-range migration. In fragmented landscapes, the evolutionary dynamics of adaptation may both decrease and enhance metapopulation viability in comparison with no evolution. The model can be applied to real patch networks with given sizes, qualities, and spatial positions of habitat patches.     

The influence of patch demographics on metapopulations, with particular reference to successional landscapes

The classical view of metapopulations relates the regional abundance of a species to the balance between the extinction and colonization dynamics of identical local populations. Species in successional landscapes may represent the most appropriate examples of classical metapopulations. However, Levins‐type metapopulation models do not explicitly separate population loss due to successional habitat change from other causes of extinction. A further complication is that the chance of population loss due to successional habitat change may be related to the age of a patch. I developed simple patch occupancy models to include succession and included consideration of patch age structure to address two related questions: what are the implications of changes in patch demographic rates and when is a move to a structured patch occupancy model justified? Age‐related variation in patch demography could increase or decrease the equilibrium fraction of the available habitat occupied by a species when compared to the predictions of an unstructured model. Metapopulation persistence was enhanced when the age class of patches with the highest species occupancy suffered relatively low losses to habitat succession. Conversely, when the age class of patches with the highest species occupancy also had relatively high successional loss rates, extinction thresholds were higher that would be predicted by a simple unstructured model. Hence age‐related variation in patch successional rate introduces biases into the predictions of simple unstructured models. Such biases can be detected from field surveys of the fraction of occupied and unoccupied patches in each age class. Where a bias is demonstrated, unstructured models will not be adequate for making predictions about the effects of changing parameters on metapopulation size. Thinking in successional terms emphasizes how landscapes might be managed to enhance or reduce the patch occupancy by any particular metapopulation     

Evolutionary suicide and evolution of dispersal in structured metapopulations

 We study the evolution of dispersal in a structured metapopulation model. The metapopulation consists of a large (infinite) number of local populations living in patches of habitable environment. Dispersal between patches is modelled by a disperser pool and individuals in transit between patches are exposed to a risk of mortality. Occasionally, local catastrophes eradicate a local population: all individuals in the affected patch die, yet the patch remains habitable. We prove that, in the absence of catastrophes, the strategy not to migrate is evolutionarily stable. Under a given set of environmental conditions, a metapopulation may be viable and yet selection may favor dispersal rates that drive the metapopulation to extinction. This phenomenon is known as evolutionary suicide. We show that in our model evolutionary suicide can occur for catastrophe rates that increase with decreasing local population size. Evolutionary suicide can also happen for constant catastrophe rates, if local growth within patches shows an Allee effect. We study the evolutionary bifurcation towards evolutionary suicide and show that a discontinuous transition to extinction is a necessary condition for evolutionary suicide to occur. In other words, if population size smoothly approaches zero at a boundary of viability in parameter space, this boundary is evolutionarily repelling and no suicide can occur. Received: 10 November 2000 / Revised version: 13 February 2002 / Published online: 17 July 2002     

Patch quality and context, but not patch number, drive multi-scale colonization dynamics in experimental aquatic landscapes

Colonization and extinction are primary drivers of local population dynamics, community structure, and spatial patterns of biological diversity. Existing paradigms of island biogeography, metapopulation biology, and metacommunity ecology, as well as habitat management and conservation biology based on those paradigms, emphasize patch size, number, and isolation as primary characteristics influencing colonization and extinction. Habitat selection theory suggests that patch quality could rival size, number, and isolation in determining rates of colonization and resulting community structure. We used naturally colonized experimental landscapes to address four issues: (a) how do colonizing aquatic beetles respond to variation in patch number, (b) how do they respond to variation in patch quality, (c) does patch context affect colonization dynamics, and (d) at what spatial scales do beetles respond to habitat variation? Increasing patch number had no effect on per patch colonization rates, while patch quality and context were critical in determining colonization rates and resulting patterns of abundance and species richness at multiple spatial scales. We graphically illustrate how variation in immigration rates driven by perceived predation risk (habitat quality) can further modify dynamics of the equilibrium theory of island biogeography beyond predator-driven effects on extinction rates. Our data support the importance of patch quality and context as primary determinants of colonization rate, occupancy, abundance, and resulting patterns of species richness, and reinforce the idea that management of metapopulations for species preservation, and metacommunities for local and regional diversity, should incorporate habitat quality into the predictive equation.     

EVOLUTION OF DISPERSAL RATES IN METAPOPULATION MODELS: BRANCHING AND CYCLIC DYNAMICS IN PHENOTYPE SPACE

We study the evolution of dispersal rates in a two patch metapopulation model. The local dynamics in each patch are given by difference equations, which, together with the rate of dispersal between the patches, determine the ecological dynamics of the metapopulation. We assume that phenotypes are given by their dispersal rate. The evolutionary dynamics in phenotype space are determined by invasion exponents, which describe whether a mutant can invade a given resident population. If the resident metapopulation is at a stable equilibrium, then selection on dispersal rates is neutral if the population sizes in the two patches are the same, while selection drives dispersal rates to zero if the local abundances are different. With non-equilibrium metapopulation dynamics, non-zero dispersal rates can be maintained by selection. In this case, and if the patches are ecologically identical, dispersal rates always evolve to values which induce synchronized metapopulation dynamics. If the patches are ecologically different, evolutionary branching into two coexisting dispersal phenotypes can be observed. Such branching can happen repeatedly, leading to polymorphisms with more than two phenotypes. If there is a cost to dispersal, evolutionary cycling in phenotype space can occur due to the dependence of selection pressures on the ecological attractor of the resident population, or because phenotypic branching alternates with the extinction of one of the branches. Our results extend those of Holt and McPeek (1996), and suggest that phenotypic branching is an important evolutionary process. This process may be relevant for sympatric speciation.     

Dispersal, migration, and offspring retention in saturated habitats

We examine the evolutionary stability of year-round residency in territorial populations, where breeding sites are a limiting resource. The model links individual life histories to the population-wide competition for territories and includes spatial variation in habitat quality as well as a potential parent-offspring conflict over territory ownership. The general form of the model makes it applicable to the evolution of dispersal, migration, partial migration, and delayed dispersal (offspring retention). We show that migration can be evolutionarily stable only if year-round residency in a given area would produce a sink population, where mortality exceeds reproduction. If this applies to a fraction of the breeding habitat only, partial migration is expected to evolve. In the context of delayed dispersal, habitat saturation has been argued to form an ecological constraint on independent breeding, which favors offspring retention and cooperative breeding. We show that habitat saturation must be considered as a dynamic outcome of birth, death, and dispersal rates in the population, rather than an externally determined constraint. Although delayed dispersal often associates with intense competition for territories, life-history traits have direct effects on stable dispersal strategies, which can often override the effect of habitat saturation. As an example, high survival of floaters selects against delayed dispersal, even though it increases the number of competitors for each breeding vacancy (the "habitat saturation factor"). High survival of territory owners, by contrast, generally favors natal philopatry. We also conclude that spatial variation in habitat quality only rarely selects for delayed dispersal. Within a population, however, offspring retention is more likely in high-quality territories.     

The evolution of patch selection in stochastic environments

A null model for habitat patch selection in spatially heterogeneous environments is the ideal free distribution (IFD), which assumes individuals have complete knowledge about the environment and can freely disperse. Under equilibrium conditions, the IFD predicts that local population growth rates are zero in all occupied patches, sink patches are unoccupied, and the fraction of the population selecting a patch is proportional to the patch's carrying capacity. Individuals, however, often experience stochastic fluctuations in environmental conditions and cannot respond to these fluctuations instantaneously. An evolutionary stability analysis for fixed patch-selection strategies reveals that environmental uncertainty disrupts the classical IFD predictions: individuals playing the evolutionarily stable strategy may occupy sink patches, local growth rates are negative and typically unequal in all patches, and individuals prefer higher-quality patches less than predicted by their carrying capacities. Spatial correlations in environmental fluctuations can enhance or marginalize these trends. The analysis predicts that continually increasing environmental variation first selects for range expansion, then selects for persisting coupled sink populations, and ultimately leads to regional extinction. In contrast, continually increasing habitat degradation first selects for range contraction and may select for persisting coupled sink populations before regional extinction. These results highlight the combined roles of spatial and temporal heterogeneity on the evolution of habitat selection.     

The evolution of dispersal – the importance of information about population density and habitat characteristics

The evolution of mobility patterns and dispersal strategies depend on different population, habitat and life history characteristics. The ability to perceive and make use of information about the surrounding environment for dispersal decisions will also differ between organisms. To investigate the evolutionary consequences of such differences, we have used a simulation model with nearest-neighbour dispersal in a metapopulation to study how variation in the ability to obtain and make use of information about habitat quality and conspecific density affects the evolution of dispersal strategies. We found a rather strong influence of variation in information on the overall rate of dispersal in a metapopulation. The highest emigration rate evolved in organisms with no information about either density or habitat quality and the lowest rate was found in organisms with information about both the natal and the neighbouring patches. For organisms that can make use of information about conspecific density, positively density-dependent dispersal evolved in the majority of cases, with the strongest density dependence occurring when an individual only has information about density in the natal patch. However, we also identified situations, involving strong local population fluctuations and frequent local extinctions, where negatively density-dependent dispersal evolved.     

Predicting which species will benefit from corridors in fragmented landscapes from population growth models

Connecting isolated patches of habitat in fragmented landscapes with corridors is a popular conservation strategy. This strategy is also controversial in large part because of uncertainty about what characteristics of a species and its environment promote corridor use. In this article we address the question, For what types of species will populations benefit from corridors? We asked this question using a model of two logistically growing populations connected by migration in which both emigration and migration success were determined by the presence or absence of a corridor. We found that in the short run (e.g., during recovery from disaster), corridors are most effective for species with fast-growing populations that have low survivorship when dispersing through unsuitable (matrix) habitat. We also found that emigration rates and habitat-specific mortality rates are key determinants of the effects of corridors on population size. In the long term, corridors are most likely to benefit species with slow-growing populations that have low survivorship when dispersing through matrix habitat. Our results confirm the major conclusions from previous empirical studies of corridor benefits. However, most studies fail to consider the most appropriate questions to determine the potential benefits of habitat corridors. First, what is the time scale of the conservation goal? Corridors have positive effects on different suites of species in the short and long term. Second, is the major threat of local extinction due to sustained population decline or boom-bust cycles? Third, what is the migration rate through the matrix? Fourth, what fraction of migrants dispersing through the matrix successfully immigrate to another patch?     

Kin selection and natal dispersal in an age-structured population

We examine the effect of iteroparity on the evolution of dispersal for a species living in a stable but fragmented habitat. We use a kin selection model that incorporates the effects of demographic stochasticity on the local age structure and age-specific genetic identities. We consider two cases: when the juvenile dispersal rate is allowed to change with maternal age and when it is not. In the latter case, we find that the unconditional evolutionarily stable dispersal rate increases when the adult survival rate increases. Two antagonistic forces act upon the evolution of age-specific dispersal rates. First, when the local age structure varies between patches of habitat, the intensity of competition between adults and juveniles in the natal patch is, on average, lower for offspring born to older senescent mothers. This selects for decreasing dispersal with maternal age. Second, offspring born to older parents are on average more related to other juveniles in the same patch and they experience a higher intensity of kin competition, which selects for increasing dispersal with maternal age. We show that the evolutionary outcome results from a balance between these two opposing forces, which depends on the amount of variance in age structure among sub-populations.     

The downward spiral: eco‐evolutionary feedback loops lead to the emergence of ‘elastic’ ranges

In times of severe environmental changes and resulting shifts in the geographical distribution of animal and plant species it is crucial to unravel the mechanisms responsible for the dynamics of species’ ranges. Without such a mechanistic understanding, reliable projections of future species distributions are difficult to derive. Species’ ranges may be highly dynamic. One particularly interesting phenomenon is range contraction following a period of expansion, referred to as ‘elastic’ behaviour. It has been proposed that this phenomenon occurs in habitat gradients, which are characterized by a negative cline in selection for dispersal from the range core towards the margin, as one may find, for example, with increasing patch isolation. Using individual‐based simulations and numerical analyses we show that Allee effects are an important determinant of range border elasticity. If only intra‐specific processes are considered, Allee effects are even a necessary condition for ranges to exhibit elastic behavior. The eco‐evolutionary interplay between dispersal evolution, Allee effects and habitat isolation leads to lower colonization probability and higher local extinction risk after range expansions, which result in an increasing amount of marginal sink patches and consequently, range contraction. We also demonstrate that the nature of the gradient is crucial for range elasticity. Gradients which do not select for lower dispersal at the margin than in the core (especially gradients in patch size, demographic stochasticity and extinction rate) do not lead to elastic range behavior. Thus, we predict that range contractions are likely to occur after periods of expansion for species living in gradients of increasing patch isolation, which suffer from Allee effects.     

Climate change in metacommunities: dispersal gives double-sided effects on persistence

Climate change is increasingly affecting the structure and dynamics of ecological communities both at local and at regional scales, and this can be expected to have important consequences for their robustness and long-term persistence. The aim of the present work is to analyse how the spatial structure of the landscape and dispersal patterns of species (dispersal rate and average dispersal distance) affects metacommunity response to two disturbances: (i) increased mortality during dispersal and (ii) local species extinction. We analyse the disturbances both in isolation and in combination. Using a spatially and dynamically explicit metacommunity model, we find that the effect of dispersal on metacommunity persistence is two-sided: on the one hand, high dispersal significantly reduces the risk of bottom-up extinction cascades following the local removal of a species; on the other hand, when dispersal imposes a risk to the dispersing individuals, high dispersal increases extinction risks, especially when dispersal is global. Large-bodied species with long generation times at the highest trophic level are particularly vulnerable to extinction when dispersal involves a risk. This suggests that decreasing the mortality risk of dispersing individuals by improving the quality of the habitat matrix may greatly increase the robustness of metacommunities.     

Non-equilibria in small metapopulations: comparing the deterministic Levins model with its stochastic counterpart

In this paper, we examine, for small metapopulations, the stochastic analog of the classical Levins metapopulation model. We study its basic model output, the expected time to metapopulation extinction, for systems which are brought out of equilibrium by imposing sudden changes in patch number and the colonization and extinction parameters. We find that the expected metapopulation extinction time shows different behavior from the relaxation time of the original, deterministic, Levins model. This relaxation time is therefore limited in value for predicting the behavior of the stochastic model. However, predictions about the extinction time for deterministically unviable cases remain qualitatively the same. Our results further suggest that, if we want to counteract the effects of habitat loss or increased dispersal resistance, the optimal conservation strategy is not to restore the original situation, that is, to create habitat or decrease resistance against dispersal. As long as the costs for different management options are not too dissimilar, it is better to improve the quality of the remaining habitat in order to decrease the local extinction rate.     

Small-world networks decrease the speed of Muller's ratchet

Muller's ratchet is an evolutionary process that has been implicated in the extinction of asexual species, the evolution of non-recombining genomes, such as the mitochondria, the degeneration of the Y chromosome, and the evolution of sex and recombination. Here we study the speed of Muller's ratchet in a spatially structured population which is subdivided into many small populations (demes) connected by migration, and distributed on a graph. We studied different types of networks: regular networks (similar to the stepping-stone model), small-world networks and completely random graphs. We show that at the onset of the small-world network - which is characterized by high local connectivity among the demes but low average path length - the speed of the ratchet starts to decrease dramatically. This result is independent of the number of demes considered, but is more pronounced the larger the network and the stronger the deleterious effect of mutations. Furthermore, although the ratchet slows down with increasing migration between demes, the observed decrease in speed is smaller in the stepping-stone model than in small-world networks. As migration rate increases, the structured populations approach, but never reach, the result in the corresponding panmictic population with the same number of individuals. Since small-world networks have been shown to describe well the real contact networks among people, we discuss our results in the light of the evolution of microbes and disease epidemics.