Sex-specific spatio-temporal variability in reproductive success promotes the evolution of sex-biased dispersal

Inbreeding depression, asymmetries in costs or benefits of dispersal, and the mating system have been identified as potential factors underlying the evolution of sex-biased dispersal. We use individual-based simulations to explore how the mating system and demographic stochasticity influence the evolution of sex-specific dispersal in a metapopulation with females competing over breeding sites, and males over mating opportunities. Comparison of simulation results for random mating with those for a harem system (locally, a single male sires all offspring) reveal that even extreme variance in local male reproductive success (extreme male competition) does not induce male-biased dispersal. The latter evolves if the between-patch variance in reproductive success is larger for males than females. This can emerge due to demographic stochasticity if the habitat patches are small. More generally, members of a group of individuals experiencing higher spatio-temporal variance in fitness expectations may evolve to disperse with greater probability than others.     

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.     

Evolution of dispersal in metapopulations with local density dependence and demographic stochasticity

In this paper, we predict the outcome of dispersal evolution in metapopulations based on the following assumptions: (i) population dynamics within patches are density-regulated by realistic growth functions; (ii) demographic stochasticity resulting from finite population sizes within patches is accounted for; and (iii) the transition of individuals between patches is explicitly modelled by a disperser pool. We show, first, that evolutionarily stable dispersal rates do not necessarily increase with rates for the local extinction of populations due to external disturbances in habitable patches. Second, we describe how demographic stochasticity affects the evolution of dispersal rates: evolutionarily stable dispersal rates remain high even when disturbance-related rates of local extinction are low, and a variety of qualitatively different responses of adapted dispersal rates to varied levels of disturbance become possible. This paper shows, for the first time, that evolution of dispersal rates may give rise to monotonically increasing or decreasing responses, as well as to intermediate maxima or minima.     

Population demography and the evolution of helping behaviors

Limited dispersal may favor the evolution of helping behaviors between relatives as it increases their relatedness, and it may inhibit such evolution as it increases local competition between these relatives. Here, we explore one way out of this dilemma: if the helping behavior allows groups to expand in size, then the kin-competition pressure opposing its evolution can be greatly reduced. We explore the effects of two kinds of stochasticity allowing for such deme expansion. First, we study the evolution of helping under environmental stochasticity that may induce complete patch extinction. Helping evolves if it results in a decrease in the probability of extinction or if it enhances the rate of patch recolonization through propagules formed by fission of nonextinct groups. This mode of dispersal is indeed commonly found in social species. Second, we consider the evolution of helping in the presence of demographic stochasticity. When fecundity is below its value maximizing deme size (undersaturation), helping evolves, but under stringent conditions unless positive density dependence (Allee effect) interferes with demographic stochasticity. When fecundity is above its value maximizing deme size (oversaturation), helping may also evolve, but only if it reduces negative density-dependent competition.     

The evolution of dispersal in reserve networks

The fragmentation of an environment into developed and protected areas may influence selection pressure on dispersal by increasing the chance of moving from a favorable to an unfavorable habitat. We theoretically explore this possibility through two cases: (1) marine systems in which reduced predation and/or increased feeding drive the evolution of planktonic larval duration and (2) more generally, where stochasticity in reproductive yield drives the evolution of the proportion of offspring dispersing. Model results indicate that habitat fragmentation generally shifts selection pressure toward reduced dispersal, particularly when areas outside reserves are uninhabitable. However, shifts to increased dispersal may occur when temporal heterogeneity is the primary selective force and constant-quota harvest occurs outside reserves. In addition, model results suggest the potential for changes in the genetic variability in dispersal after habitat fragmentation. The predicted evolutionary changes in dispersal will depend on factors such as the relative genetic and environmental contributions to dispersal-related traits and the extent of anthropogenic impacts outside reserves. If the predicted evolutionary changes are biologically attainable, they may suggest altering current guidelines for the appropriate size and spacing of marine reserves necessary to achieve conservation and fisheries goals.     

Genetic structuring in fluctuating populations

The effect of genetic drift in spatially distributed dispersal-linked and density-regulated populations is studied in a classical one-locus two-allele system. We analyse emergence of genetic differentiation assuming random drift only, where the noise-like variability is due to demographic stochasticity. We find emergence of clusters of sub-units with local allele fixation and persistence of both alleles in lengthy simulations. We demonstrate that local allele fixation (extending over a number of adjoining spatial sub-units) – without global loss of alleles – may occur when the carrying capacities of local patches are small, under a full range population dynamic regimes, when dispersal rate is small, and when redistribution (through dispersal) does not act as global mixer. These results are novel. The key to the observations is that drift is simultaneously influenced by distance-dependent dispersal, demographic stochasticity and autocorrelated population fluctuations due to delayed-density dependence. These are standard elements of contemporary population models in spatially structured context. With stable large populations, no stochasticity and dispersal limited to neighbours only, our model collapses to the stepping-stone model, while with dispersal being random and global, the model collapses to Wright's island model.     

Evolution of positive and negative density-dependent dispersal

Understanding the evolution of density-dependent dispersal strategies has been a major challenge for evolutionary ecologists. Some existing models suggest that selection should favour positive and others negative density-dependence in dispersal. Here, we develop a general model that shows how and why selection may shift from positive to negative density-dependence in response to key ecological factors, in particular the temporal stability of the environment. We find that in temporally stable environments, particularly with low dispersal costs and large group sizes, habitat heterogeneity selects for negative density-dependent dispersal, whereas in temporally variable environments, particularly with high dispersal costs and small group sizes, habitat heterogeneity selects for positive density-dependent dispersal. This shift reflects the changing balance between the greater competition for breeding opportunities in more productive patches, versus the greater long-term value of offspring that establish themselves there, the latter being very sensitive to the temporal stability of the environment. In general, dispersal of individuals out of low-density patches is much more sensitive to habitat heterogeneity than is dispersal out of high-density patches.     

Evolution of Complex Density-Dependent Dispersal Strategies

The question of how dispersal behavior is adaptive and how it responds to changes in selection pressure is more relevant than ever, as anthropogenic habitat alteration and climate change accelerate around the world. In metapopulation models where local populations are large, and thus local population size is measured in densities, density-dependent dispersal is expected to evolve to a single-threshold strategy, in which individuals stay in patches with local population density smaller than a threshold value and move immediately away from patches with local population density larger than the threshold. Fragmentation tends to convert continuous populations into metapopulations and also to decrease local population sizes. Therefore we analyze a metapopulation model, where each patch can support only a relatively small local population and thus experience demographic stochasticity. We investigated the evolution of density-dependent dispersal, emigration and immigration, in two scenarios: adult and natal dispersal. We show that density-dependent emigration can also evolve to a nonmonotone, “triple-threshold” strategy. This interesting phenomenon results from an interplay between the direct and indirect benefits of dispersal and the costs of dispersal. We also found that, compared to juveniles, dispersing adults may benefit more from density-dependent vs. density-independent dispersal strategies.     

Sex ratio evolution when fitness and dispersal vary

In a heterogeneous environment, when the fitness of males and females are differently influenced by habitat quality, habitat-dependent sex ratios may evolve to favor the production of the sex that benefits more (or loses less) from the local habitat. Similarly, sex-biased dispersal favors the evolution of habitat-dependent sex ratios. The present study documents the convergence stable sex ratios expected in the presence of sex-specific fitness gains when dispersal is partial, sex-biased or costly, using a simple model with patches of two qualities. Results show that partial dispersal reduces the sex ratio bias expected with sex-specific fitness gains. The direction of the sex ratio bias can be reversed by sex-biased dispersal or the existence of sex-specific dispersal costs, provided that fitness gains for the two sexes are not too different. The reversal of the sex ratio bias is more readily observed when sex-specific dispersal rates are opposite and extreme. Both dispersal and fitness gains, especially when they are sex-specific, should thus be considered when making predictions about sex ratio evolution in a heterogeneous environment.     

Metapopulation extinction risk: Dispersal’s duplicity

Metapopulation extinction risk is the probability that all local populations are simultaneously extinct during a fixed time frame. Dispersal may reduce a metapopulation’s extinction risk by raising its average per-capita growth rate. By contrast, dispersal may raise a metapopulation’s extinction risk by reducing its average population density. Which effect prevails is controlled by habitat fragmentation. Dispersal in mildly fragmented habitat reduces a metapopulation’s extinction risk by raising its average per-capita growth rate without causing any appreciable drop in its average population density. By contrast, dispersal in severely fragmented habitat raises a metapopulation’s extinction risk because the rise in its average per-capita growth rate is more than offset by the decline in its average population density. The metapopulation model used here shows several other interesting phenomena. Dispersal in sufficiently fragmented habitat reduces a metapopulation’s extinction risk to that of a constant environment. Dispersal between habitat fragments reduces a metapopulation’s extinction risk insofar as local environments are asynchronous. Grouped dispersal raises the effective habitat fragmentation level. Dispersal search barriers raise metapopulation extinction risk. Nonuniform dispersal may reduce the effective fraction of suitable habitat fragments below the extinction threshold. Nonuniform dispersal may make demographic stochasticity a more potent metapopulation extinction force than environmental stochasticity.     

Habitat heterogeneities versus spatial type frequency variances as driving forces of dispersal evolution

Understanding the evolution of dispersal is essential for understanding and predicting the dynamics of natural populations. Two main factors are known to influence dispersal evolution: spatio‐temporal variation in the environment and relatedness between individuals. However, the relation between these factors is still poorly understood, and they are usually treated separately. In this article, I present a theoretical framework that contains and connects effects of both environmental variation and relatedness, and reproduces and extends their known features. Spatial habitat variation selects for balanced dispersal strategies, whereby the population is kept at an ideal free distribution. Within this class of dispersal strategies, I explain how increased dispersal is promoted by perturbations to the dispersal type frequencies. An explicit formula shows the magnitude of the selective advantage of increased dispersal in terms of the spatial variability in the frequencies of the different dispersal strategies present. These variances are capable of capturing various sources of stochasticity and hence establish a common scale for their effects on the evolution of dispersal. The results furthermore indicate an alternative approach to identifying effects of relatedness on dispersal evolution.     

Local and Regional Determinants of Colonisation-Extinction Dynamics of a Riparian Mainland-Island Root Vole Metapopulation

The role of local habitat geometry (habitat area and isolation) in predicting species distribution has become an increasingly more important issue, because habitat loss and fragmentation cause species range contraction and extinction. However, it has also become clear that other factors, in particular regional factors (environmental stochasticity and regional population dynamics), should be taken into account when predicting colonisation and extinction. In a live trapping study of a mainland-island metapopulation of the root vole (Microtus oeconomus) we found extensive occupancy dynamics across 15 riparian islands, but yet an overall balance between colonisation and extinction over 4 years. The 54 live trapping surveys conducted over 13 seasons revealed imperfect detection and proxies of population density had to be included in robust design, multi-season occupancy models to achieve unbiased rate estimates. Island colonisation probability was parsimoniously predicted by the multi-annual density fluctuations of the regional mainland population and local island habitat quality, while extinction probability was predicted by island population density and the level of the recent flooding events (the latter being the main regionalized disturbance regime in the study system). Island size and isolation had no additional predictive power and thus such local geometric habitat characteristics may be overrated as predictors of vole habitat occupancy relative to measures of local habitat quality. Our results suggest also that dynamic features of the larger region and/or the metapopulation as a whole, owing to spatially correlated environmental stochasticity and/or biotic interactions, may rule the colonisation – extinction dynamics of boreal vole metapopulations. Due to high capacities for dispersal and habitat tracking voles originating from large source populations can rapidly colonise remote and small high quality habitat patches and re-establish populations that have gone extinct due to demographic (small population size) and environmental stochasticity (e.g. extreme climate events).     

The role of colonization in the dynamics of patchy populations of a cyclic vole species

The crash phase of vole populations with cyclic dynamics regularly leads to vast areas of uninhabited habitats. Yet although the capacity for cyclic voles to re-colonize such empty space is likely to be large and predicted to have become evolved as a distinct life history trait, the processes of colonization and its effect on the spatio-temporal dynamics have been little studied. Here we report from an experiment with root voles (Microtus oeconomus) specifically targeted at quantifying the process of colonization of empty patches from distant source patches and its resultant effect on local vole deme size variation in a patchy landscape. Three experimental factors: habitat quality, predation risk and inter-patch distance were employed among 24 habitat patches in a 100 × 300-m experimental area. The first-born cohort in the spring efficiently colonized almost all empty patches irrespective of the degree of patch isolation and predation risk, but this was dependent on habitat quality. Just after the initial colonization wave the deme sizes in patches of the same quality were underdispersed relative to Poisson variance, indicating regulated (density-dependent) settlement. Towards the end of the breeding season local demographic processes acted to smooth out the initial post-colonization differences among source and colonization patches, and among patches of initially different quality. However, at this time demographic stochasticity had also given rise to a large (overdispersed) variation in deme sizes that may have contributed to an overshadowing of the effect of other factors. The results of this experiment confirmed our expectation that the space-filling capacity of voles is large. The costs associated with transience appeared to be so low, at least at the spatial scale considered in this experiment, that such costs are not likely to substantially constrain habitat selection and colonization in the increase phase of cyclic patchy populations.     

Habitat destruction, environmental catastrophes, and metapopulation extinction

The extinction process of fragmented populations, characterized by a small number of conspecifics inhabiting each patch, is heavily affected by natural and human disturbance. To evaluate the risk of extinction we consider a network of identical patches connected by passive or active dispersal and hosting a finite, discrete number of individuals. We discuss three types of disturbance affecting the metapopulation: permanent loss of habitat patches, erosion of existing patches, and random catastrophes that wipe out the entire population of a patch. Starting from an infinite-dimensional Markov model that fully accounts for demographic stochasticity, we reduce it to finite dimension via moment closure with negative-binomial approximation. The compact models obtained in this way account for the dynamics of the fraction of empty patches, the average number of individuals in occupied patches, and the variance of their distribution. After comparing the performance of these compact models with that of the infinite-dimensional model in the case of no disturbances, we then proceed to computing persistence-extinction boundaries as bifurcation lines of the compact models in the space of demographic and disturbance parameters. We consider bifurcations with respect to demographic and environmental parameters and contrast our results with those of previous theories. We find out that environmental catastrophes increase the risk of extinction for both frequent and infrequent dispersers, while the random loss of patches has a much larger influence on frequent dispersers. This influence can be counterbalanced by active dispersal. Local erosion of habitat fragments has a larger influence on infrequent than on frequent dispersers. We finally discuss the important synergistic effects of disturbances acting simultaneously.     

Evolution of body condition‐dependent dispersal in metapopulations

Body condition‐dependent dispersal strategies are common in nature. Although it is obvious that environmental constraints may induce a positive relationship between body condition and dispersal, it is not clear whether positive body conditional dispersal strategies may evolve as a strategy in metapopulations. We have developed an individual‐based simulation model to investigate how body condition–dispersal reaction norms evolve in metapopulations that are characterized by different levels of environmental stochasticity and dispersal mortality. In the model, body condition is related to fecundity and determined either by environmental conditions during juvenile development (adult dispersal) or by those experienced by the mother (natal dispersal). Evolutionarily stable reaction norms strongly depend on metapopulation conditions: positive body condition dependency of dispersal evolved in metapopulation conditions with low levels of dispersal mortality and high levels of environmental stochasticity. Negative body condition‐dependent dispersal evolved in metapopulations with high dispersal mortality and low environmental stochasticity. The latter strategy is responsible for higher dispersal rates under kin competition when dispersal decisions are based on body condition reached at the adult life stage. The evolution of both positive and negative body condition‐dependent dispersal strategies is consequently likely in metapopulations and depends on the prevalent environmental conditions.     

The social evolution of dispersal with public goods cooperation

Selection can favour the evolution of individually costly dispersal if this alleviates competition between relatives. However, conditions that favour altruistic dispersal also mediate selection for other social behaviours, such as public goods cooperation, which in turn is likely to mediate dispersal evolution. Here, we investigate – both experimentally (using bacteria) and theoretically – how social habitat heterogeneity (i.e. the distribution of public goods cooperators and cheats) affects the evolution of dispersal. In addition to recovering the well‐known theoretical result that the optimal level of dispersal increases with genetic relatedness of patch mates, we find both mathematically and experimentally that dispersal is always favoured when average patch occupancy is low, but when average patch occupancy is high, the presence of public goods cheats greatly alters selection for dispersal. Specifically, when public goods cheats are localized to the home patch, higher dispersal rates are favoured, but when cheats are present throughout available patches, lower dispersal rates are favoured. These results highlight the importance of other social traits in driving dispersal evolution.     

Conditional dispersal, clines, and the evolution of dispersiveness

Conditional dispersal, in which an individual’s decision over whether to disperse is a response to environmental conditions, features prominently in studies of dispersal evolution. Using models of clines, I examine how one widely discussed cost of dispersal, namely, that dispersal impedes local adaptation, changes with conditional dispersal and what this implies for dispersal evolution. I examine the consequences for dispersal evolution of the responsiveness of dispersal to the environment, the accuracy of any proximal cues that individuals rely upon to assess habitat quality, and whether dispersal responds to fitness itself or only to some fitness components (juvenile survivorship). All of the conditional dispersal behaviors that I consider weaken the indirect cost of dispersal inhibiting local adaptation. However, if individuals rely on imprecise cues to assess habitat quality and base dispersal decisions on juvenile survivorship, then conditional dispersal can incur additional costs by exacerbating overcrowding. Conditional dispersal initially leads to steeper clines in traits under direct selection, but when dispersiveness can itself evolve, conditional dispersal allows sigmoidal clines to persist long after those obtained with unconditional movement would become stepped. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Single large AND several small habitat patches: A community perspective on their importance for biodiversity

The debate whether single large or several small (SLOSS) patches benefit biodiversity has existed for decades, but recent literature provides increasing evidence for the importance of small habitats. Possible beneficial mechanisms include reduced presence of predators and competitors in small habitat areas or specific functions such as stepping stones for dispersal. Given the increasing amount of studies highlighting individual behavioral differences that may influence these functions, we hypothesize that the advantage of small versus large habitat patches not only depends on patch functionality but also on the presence of animal personalities (i.e., risk-tolerant vs. risk-averse). Using an individual-based, spatially-explicit community model, we analyzed the diversity of mammal communities in landscapes consisting of a few large habitat islands interspersed with different amounts and sizes of small habitat patches. Within these heterogeneous environments, individuals compete for resources and form home-ranges, with only risk-tolerant individuals using habitat edges. Results show that when risk-tolerant individuals exist, small patches increase species diversity. A strong peak occurs at approximately 20% habitat cover in small patches when those small habitats are only used for foraging but not for breeding and home-range core position. Additional usage as stepping stones for juvenile dispersal further increases species persistence. Overall, our results reveal that a combination of a few large and several small habitat patches promotes biodiversity by enhancing landscape heterogeneity. Here, heterogeneity is created by pronounced differences in habitat functionality, increasing edge density, and variability in habitat use by different behavioral types. The finding that a combination of single large AND several small (SLASS) patches is needed for effective biodiversity preservation has implications for advancing landscape conservation. Particularly in structurally poor agricultural areas, modern technology enables precise management with the opportunity to create small foraging habitats by excluding less profitable agricultural land from cultivation.     

Will sympatric speciation fail due to stochastic competitive exclusion?

Sympatric speciation requires coexistence of the newly formed species. If divergence proceeds by small mutational steps, the new species utilize almost the same resources initially, and full speciation may be impeded by competitive exclusion in stochastic environments. We investigate this primarily ecological problem of sympatric speciation by studying the population dynamics of a diverging asexual population in a fluctuating environment. Correlation between species responses to environmental fluctuation is assumed to decrease with distance in trait space. Rapidly declining correlation in combination with high environmental variability may delay full speciation or even render it impossible. Stochastic extinctions impeding speciation are most likely when correlation decays faster than competition, for example, when demographic stochasticity is strong or when divergence is not accompanied by niche separation, such as in speciation driven entirely by sexual selection. Our general theoretical results show an interesting connection between short-term ecological dynamics and long-term, large-scale evolution.     

Multi-scale methods predict invasion speeds in variable landscapes

Spread rates of invasive plant species depend heavily on variable seed/seedling survivorships over various habitat types as well as on variability in seed dispersal induced by rapid transport of propagules in open areas and slow transport in vegetated areas. The ability to capture spatial variability in seed survivorship and dispersal is crucial to accurately predict the rate of spread of plants in real world landscapes. However, current analytic methods for predicting spread rates are not suited for arbitrary, spatially heterogeneous systems. Here, we analyze invasion rates of the invasive plant Phragmites australis (common reed) over variable wetland landscapes. Phragmites is one of the most pervasive perennial grasses, outcompeting native vegetation, providing poor wildlife habitat, and proving difficult to eradicate across its invasive range in North America. Phragmites spreads sexually via seeds and asexually via underground (rhizomes) and aboveground (stolons) stems. We construct a structured integrodifference equation model of the Phragmites life cycle capturing variable seed survivorship in a seed bank, sexual and asexual recruitment into a juvenile age class, and differential competition among all classes with adults. The demographic model is coupled with a homogenized ecological diffusion/settling seed dispersal model that allows for seed deposition that varies with habitat type. The dispersal kernel we develop does not require local normalization and can be implemented efficiently using standard computational techniques. The model generates a traveling wave of isolated patches, establishing only in suitable habitats. We use the method of multiple scales to predict invasion speed as a solvability condition at large scales and test the predictions numerically. Accurate predictions are generated for a wide range of landscape parameters, indicating that invasion speeds can be understood in landscapes of arbitrary structure using this approach.