Spatio‐temporal dynamics of density‐dependent dispersal during a population colonisation

Predicting population colonisations requires understanding how spatio‐temporal changes in density affect dispersal. Density can inform on fitness prospects, acting as a cue for either habitat quality, or competition over resources. However, when escaping competition, high local density should only increase emigration if lower‐density patches are available elsewhere. Few empirical studies on dispersal have considered the effects of density at the local and landscape scale simultaneously. To explore this, we analyze 5 years of individual‐based data from an experimental introduction of wild guppies Poecilia reticulata. Natal dispersal showed a decrease in local density dependence as density at the landscape level increased. Landscape density did not affect dispersal among adults, but local density‐dependent dispersal switched from negative (conspecific attraction) to positive (conspecific avoidance), as the colonisation progressed. This study demonstrates that densities at various scales interact to determine dispersal, and suggests that dispersal trade‐offs differ across life stages.     

LOCAL ADAPTATION WHEN COMPETITION DEPENDS ON PHENOTYPIC SIMILARITY

Recent work incorporating demographic–genetic interactions indicates the importance of population size, gene flow, and selection in influencing local adaptation. This work typically assumes that density‐dependent survival affects individuals equally, but individuals in natural population rarely compete equally. Among‐individual differences in resource use generate stronger competition between more similar phenotypes (frequency‐dependent competition) but it remains unclear how this additional form of selection changes the interactions between population size, gene flow, and local stabilizing selection. Here, we integrate migration–selection dynamics with frequency‐dependent competition. We developed a coupled demographic‐quantitative genetic model consisting of two patches connected by dispersal and subject to local stabilizing selection and competition. Our model shows that frequency‐dependent competition slightly increases local adaptation, greatly increases genetic variance within patches, and reduces the amount that migration depresses population size, despite the increased genetic variance load. The effects of frequency‐dependence depend on the strength of divergent selection, trait heritability, and when mortality occurs in the life cycle in relation to migration and reproduction. Essentially, frequency‐dependent competition reduces the density‐dependent interactions between migrants and residents, the extent to which depends on how different and common immigrants are compared to residents. Our results add new dynamics that illustrate how competition can alter the effects of gene flow and divergent selection on local adaptation and population carrying capacities.     

Density-dependent dispersal in birds and mammals

Density‐dependent dispersal can be caused by various mechanisms, from competition inducing individuals to emigrate (positive density‐dependence) to social crowding effects impeding free movement (negative density‐dependence). Various spatial population models have incorporated positively density‐dependent dispersal algorithms, and recent theoretical models have explored the conditions for density‐dependent dispersal (DD) to evolve. However, while the existence of DD is well documented in some taxa such as insects, there is no clear picture on its generality in vertebrates. Here I review the available empirical data on DD in birds and mammals, focusing mainly on variation in dispersal between years and on experimental density manipulations. Surprisingly few studies have explicitly focused on DD, and interpretation of the available data is often hampered by differences in approach, small sample sizes and/or statistical shortcomings. Positive DD was reported in 50 and 33% of the selected mammal and bird studies, respectively, while two studies on mammals (out of eight) reported negative DD. Among bird studies, DD was more often reported for emigration rates or long‐distance recoveries than for average distances within finite study areas. Experimental studies manipulating densities (mainly on mammals) have consistently generated positive DD, typically showing reduced emigration in response to partial population removal. Studies that examined dispersal in relation to seasonal changes in density (small mammals only) have more often reported negative DD. Studies that compared dispersal between sites differing in density, also show a mixture of positive and negative DD. This suggests that dispersal changes in a more complex way with seasonal and spatial density variation than with annual densities, and/or that these results are confounded by other factors differing between seasons and sites, such as habitat quality. I conclude that both correlational and experimental studies support the existence of positive, rather than negative, density‐dependent dispersal in birds and mammals.     

Selfing ability and dispersal are positively related,but not affected by range position: a multispecies study on southern African Asteraceae

Dispersal and breeding system traits are thought to affect colonization success. As species have attained their present distribution ranges through colonization, these traits may vary geographically. Although several theories predict associations between dispersal ability, selfing ability and the relative position of a population within its geographic range, there is little theoretical or empirical consensus on exactly how these three variables are related. We investigated relationships between dispersal ability, selfing ability and range position across 28 populations of 13 annual, wind‐dispersed Asteraceae species from the Namaqualand region of South Africa. Controlling for phylogeny, relative dispersal ability – assessed from vertical fall time of fruits – was positively related to an index of autofertility – determined from hand‐pollination experiments. These findings support the existence of two discrete syndromes: high selfing ability associated with good dispersal and obligate outcrossing associated with lower dispersal ability. This is consistent with the hypothesis that selection for colonization success drives the evolution of an association between these traits. However, no general effect of range position on dispersal or breeding system traits was evident. This suggests selection on both breeding system and dispersal traits acts consistently across distribution ranges.     

Natal habitat effects drive density—dependent scaling of dispersal decisions

The dispersal behavior of a species is critical for the stability and persistence of its populations across a landscape. How population density affects dispersal decisions is important for predicting these dynamics, as the form of density‐dependent dispersal influences the stability and persistence of populations. Natal habitat experience often has strong impacts on individual dispersal behavior as well, but its influence on density‐dependent dispersal behaviors remains unexplored. Here we address this conceptual gap in two experiments separately examining habitat selection and emigration from recently colonized patches for two species of flour beetle Tribolium sp. We found that interactions between the quality of habitat experienced during natal development and current habitat for dispersal capable adults can strongly affect the form of density dependence, including reversing the direction of nonlinearities (accelerating to decelerating), or even negating the influence of population density for individual dispersal decisions. Across heterogeneous landscapes, where individuals from different populations may experience different natal habitats, this altering of density‐dependent relationships is predicted by theory to fundamentally influence regional population dynamics. Our results indicate that species which occur across heterogeneous environments, such as during conservation reintroductions, or as invasive species spread, have much potential for natal experience to interact with density dependence and influence local and regional population dynamics.     

Source–sink dynamics: how sinks affect demography of sources

Models of source–sink population dynamics have to make assumptions about whether, and eventually how, demographic parameters in source habitats are dependent on the demography in sink habitats. However, the empirical basis for making such assumptions has been weak. Here we report a study on experimental root vole populations, where estimates of demographic parameters were contrasted between source patches in source–sink (treatment) and source–source systems (control). In the presence of a sink patch (simulated by a pulsed removal of immigrants), source‐patch populations failed to increase over the breeding season, mainly due to a high spatially density‐dependent dispersal rate from source to sink patches. The per capita recruitment rate was almost two times higher in source–sink than in the source–source systems, but this did not compensate for the loss rate due to dispersal from source to sink patches. Sex ratio in the source–sink systems became less female biased, probably as a result of an enhanced frequency of dispersal movements in females. Good knowledge of the degree of density‐and habitat‐dependent dispersal is critical for predicting the dynamics of source–sink populations.     

Joint evolution of differential seed dispersal and self‐fertilization

Differential seed dispersal, in which selfed and outcrossed seeds possess different dispersal propensities, represents a potentially important individual‐level association. A variety of traits can mediate differential seed dispersal, including inflorescence and seed size variation. However, how natural selection shapes such associations is poorly known. Here, we developed theoretical models for the evolution of mating system and differential seed dispersal in metapopulations, incorporating heterogeneous pollination, dispersal cost, cost of outcrossing and environment‐dependent inbreeding depression. We considered three models. In the ‘fixed dispersal model’, only selfing rate is allowed to evolve. In the ‘fixed selfing model’, in which selfing is fixed but differential seed dispersal can evolve, we showed that natural selection favours a higher, equal or lower dispersal rate for selfed seeds to that for outcrossed seeds. However, in the ‘joint evolution model’, in which selfing and dispersal can evolve together, evolution necessarily leads to higher or equal dispersal rate for selfed seeds compared to that for outcrossed. Further comparison revealed that outcrossed seed dispersal is selected against by the evolution of mixed mating or selfing, whereas the evolution of selfed seed dispersal undergoes independent processes. We discuss the adaptive significance and constraints for mating system/dispersal association.     

The magnitude of local adaptation under genotype‐dependent dispersal

Dispersal moves individuals from patches where their immediate ancestors were successful to sites where their genotypes are untested. As a result, dispersal generally reduces fitness, a phenomenon known as “migration load.” The strength of migration load depends on the pattern of dispersal and can be dramatically lessened or reversed when individuals move preferentially toward patches conferring higher fitness. Evolutionary ecologists have long modeled nonrandom dispersal, focusing primarily on its effects on population density over space, the maintenance of genetic variation, and reproductive isolation. Here, we build upon previous work by calculating how the extent of local adaptation and the migration load are affected when individuals differ in their dispersal rate in a genotype‐dependent manner that alters their match to their environment. Examining a one‐locus, two‐patch model, we show that local adaptation occurs through a combination of natural selection and adaptive dispersal. For a substantial portion of parameter space, adaptive dispersal can be the predominant force generating local adaptation. Furthermore, genetic load may be largely averted with adaptive dispersal whenever individuals move before selective deaths occur. Thus, to understand the mechanisms driving local adaptation, biologists must account for the extent and nature of nonrandom, genotype‐dependent dispersal, and the potential for adaptation via spatial sorting of genotypes.     

Body condition dependent dispersal in a heterogeneous environment

We find the evolutionarily stable dispersal behaviour of a population that inhabits a heterogeneous environment where patches differ in safety (the probability that a juvenile individual survives until reproduction) and productivity (the total competitive weight of offspring produced by the local individual), assuming that these characteristics do not change over time. The body condition of clonally produced offspring varies within and between families. Offspring compete for patches in a weighted lottery, and dispersal is driven by kin competition. Survival during dispersal may depend on body condition, and competitive ability increases with increasing body condition. The evolutionarily stable strategy predicts that families abandon patches which are too unsafe or do not produce enough successful dispersers. From families that invest in retaining their natal patches, individuals stay in the patch that are less suitable for dispersal whereas the better dispersers disperse. However, this clear within-family pattern is often not reflected in the population-wide body condition distribution of dispersers or non-dispersers. This may be an explanation why empirical data do not show any general relationship between body condition and dispersal. When all individuals are equally good dispersers, then there exist equivalence classes defined by the competitive weight that remains in a patch. An equivalence class consists of infinitely many dispersal strategies that are selectively neutral. This provides an explanation why very diverse patterns found in body condition dependent dispersal data can all be equally evolutionarily stable.     

Metapopulation dynamics of the bog fritillary butterfly: experimental changes in habitat quality induced negative density‐dependent dispersal

Habitat quality and habitat geometry are two crucial factors driving metapopulation dynamics. However, their intricacy has prevented so far a reliable test of their relative impact on local population dynamics and persistence. Here we report on a long‐term study in which we manipulated habitat quality within a butterfly metapopulation, whereas habitat geometry was kept constant. The treatment consisted in lowering the quality of certain habitat patches while others were kept untreated, using the same spatial design over years. The effect of the treatment on metapopulation dynamics was assessed by comparing residence probability and dispersal rates within the same habitat network on 11 and 6 independent butterfly generations before and after treatment, respectively. Results showed that the experimental decrease in habitat quality generated significantly higher emigration rates from treated patches. This increase was associated with a significant decrease in dispersal rates out of untreated patches, and a significant higher residence probability in these patches. The direct relation between lower habitat quality and higher dispersal propensity in treated patches was expected. However, the lower dispersal from untreated patches after treatment was opposite to the expectation of positive density dependent dispersal generally observed in butterflies. Such negative density‐dependent dispersal would allow a rapid fine‐tuning of dispersal rates to changes in habitat quality, particularly when the spatial autocorrelation of the environmental is low. Accordingly, dispersal would promote an ideal free distribution of individuals in the landscape according to their fitness expectation.     

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.     

Temporal autocorrelation can enhance the persistence and abundance of metapopulations comprised of coupled sinks

In spatially heterogeneous landscapes, some habitats may be persistent sources, providing immigrants to sustain populations in unfavorable sink habitats (where extinction is inevitable without immigration). Recent theoretical and empirical studies of source-sink systems demonstrate that temporally variable local growth rates in sinks can substantially increase average abundance of a persisting population, provided that the variation is positively autocorrelated--in effect, temporal variation inflates average abundance. Here we extend these results to a metapopulation in which all habitat patches are sinks. Using numerical studies of a population with discrete generations (buttressed by analytic results), we show that temporal variation and moderate dispersal can jointly permit indefinite persistence of the metapopulation and that positive autocorrelation both lowers the magnitude of variation required for persistence and increases the average abundance of persisting metapopulations. These effects are weakened--but not destroyed--if variation in local growth rates is spatially synchronized and dispersal is localized. We show that the inflationary effect is robust to a number of extensions of the basic model, including demographic stochasticity and density dependence. Because ecological and environmental processes contributing to temporally variable growth rates in natural populations are typically autocorrelated, these observations may have important implications for species persistence.     

Allee effect and self-fertilization in hermaphrodites: reproductive assurance in a structured metapopulation

Reproductive assurance through selfing during colonization events or when population densities are low has often been put forward as a mechanism selecting for the evolution of self-fertilization. Such arguments emphasize on the role of both local demography and metapopulation processes. We developed a model for the evolution of self-fertilization in a structured metapopulation in which local densities are not steady because of population growth. Reproduction by selfing is density-independent (reproductive assurance) but selfed seeds endure inbreeding depression, whereas reproduction by outcrossing is density-dependent (Allee effect). First, we derived an analytical criterion for metapopulation viability as a function of the selfing rate and metapopulation parameters. We show that outcrossers can develop a viable metapopulation when they produce a high amount of dispersal seeds that counterbalances their incapacity to found new populations from low densities. Second, the model shows there is a positive feedback between demography and outcrossing rates, leading to either complete outcrossing or selfing. Specifically, we illustrate that inbreeding depression can paradoxically favor the evolution of selfing because of its negative effect on density. Also, complete outcrossing can be selected despite pollen limitation, although it does not provide a full seed set. This model underlines the influence of the mating system both on demography and gene dynamics in a metapopulation context.     

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.     

Negative density-dependent dispersal emerges from the joint evolution of density- and body condition-dependent dispersal strategies

Empirical studies have documented both positive and negative density-dependent dispersal, yet most theoretical models predict positive density dependence as a mechanism to avoid competition. Several hypotheses have been proposed to explain the occurrence of negative density-dependent dispersal, but few of these have been formally modeled. Here, we developed an individual-based model of the evolution of density-dependent dispersal. This model is novel in that it considers the effects of density on dispersal directly, and indirectly through effects on individual condition. Body condition is determined mechanistically, by having juveniles compete for resources in their natal patch. We found that the evolved dispersal strategy was a steep, increasing function of both density and condition. Interestingly, although populations evolved a positive density-dependent dispersal strategy, the simulated metapopulations exhibited negative density-dependent dispersal. This occurred because of the negative relationship between density and body condition: high density sites produced low-condition individuals that lacked the resources required for dispersal. Our model, therefore, generates the novel hypothesis that observed negative density-dependent dispersal can occur when high density limits the ability of organisms to disperse. We suggest that future studies consider how phenotype is linked to the environment when investigating the evolution of dispersal.     

Evidence for positive density-dependent emigration in butterfly metapopulations

A positive effect of (meta)population density on emigration has been predicted by many theoretical models and confirmed empirically in various organisms. However, in butterflies, the most popular species for dispersal studies, the evidence for its existence has so far been equivocal, with negative relationships between density and emigration being reported more frequently. We analysed dispersal in sympatric metapopulations of two Maculinea butterflies, intensively surveyed with mark–release–recapture methods for 7 years. Dispersal parameters, derived using the virtual migration model, were assessed against butterfly densities, which fluctuated strongly over the study period. Emigration was positively correlated with density, and this effect was particularly strong at densities above carrying capacity, when emigration increased up to threefold in females and twofold in males compared with the normal levels. In turn, density had little impact on other dispersal parameters analysed. Our findings provide good evidence for positive density-dependence of emigration in butterflies. Emigrating at high densities is particularly beneficial for females, because it gives them a chance to lay part of their egg-load in less crowded patches, where offspring survival is higher due to lower intraspecific competition. Even though the rise in emigration becomes considerable at densities exceeding carrying capacity, i.e. relatively infrequently, it still has serious implications for many ecological phenomena, such as species range expansions, gene flow, and metapopulation persistence. Consequently, instead of treating emigration as a fixed trait, it is worth allowing for its density-dependence in applications such as population viability analyses, genetic models or metapopulation models.     

Regional and local scale metapopulation dynamics in the interaction between Callosobruchus maculatus and Anisopteromalus calandrae

Habitat structure increases the persistence of many extinction‐prone resource–consumer interactions. Metapopulation theory is one of the leading approaches currently used to explain why local, ephemeral populations persist at a regional scale. Central to the metapopulation concept is the amount of dispersal occurring between patches, too much or too little can result in regional extinction. In this study, the role of dispersal on the metapopulation dynamics of an over‐exploitative host–parasitoid interaction is assessed. In the absence of the parasitoid the highly vagile bruchid, Callosobruchus maculatus, can maintain a similar population size regardless of the permeability of the inter‐patch matrix and exhibits strong negative density‐dependence. After the introduction of the parasitoid the size of the bruchid population decreases with a corresponding increase in the occurrence of empty patches. In this case, limiting the dispersal of both species decouples the interaction to a greater extent and results in larger regional bruchid populations. Given the disparity between the dispersal rates of the two species, it is proposed that the more dispersive host benefits from the reduction in landscape permeability by increasing the opportunity to colonise empty patches and rescue extinction prone populations. Associated with the introduction of the parasitoid is a shift in the strength of density‐dependence as the population moves from bottom–up towards top–down regulation. The importance of local and regional scale measurements is apparent when the role of individual patches on regional dynamics is considered. By only taking regional dynamics into account the importance of dispersal regime on local dynamics is overlooked. Similarly, when local dynamics were examined, patches were found to have different influences on regional dynamics depending on dispersal regime and patch location.     

Cooperation‐mediated plasticity in dispersal and colonization

Kin selection theory predicts that costly cooperative behaviors evolve most readily when directed toward kin. Dispersal plays a controversial role in the evolution of cooperation: dispersal decreases local population relatedness and thus opposes the evolution of cooperation, but limited dispersal increases kin competition and can negate the benefits of cooperation. Theoretical work has suggested that plasticity of dispersal, where individuals can adjust their dispersal decisions according to the social context, might help resolve this paradox and promote the evolution of cooperation. Here, we experimentally tested the hypothesis that conditional dispersal decisions are mediated by a cooperative strategy: we quantified the density‐dependent dispersal decisions and subsequent colonization efficiency from single cells or groups of cells among six genetic strains of the unicellular Tetrahymena thermophila that differ in their aggregation level (high, medium, and low), a behavior associated with cooperation strategy. We found that the plastic reaction norms of dispersal rate relative to density differed according to aggregation level: highly aggregative genotypes showed negative density‐dependent dispersal, whereas low‐aggregation genotypes showed maximum dispersal rates at intermediate density, and medium‐aggregation genotypes showed density‐independent dispersal with intermediate dispersal rate. Dispersers from highly aggregative genotypes had specialized long‐distance dispersal phenotypes, contrary to low‐aggregation genotypes; medium‐aggregation genotypes showing intermediate dispersal phenotype. Moreover, highly aggregation genotypes showed evidence for beneficial kin‐cooperation during dispersal. Our experimental results should help to resolve the evolutionary conflict between cooperation and dispersal: cooperative individuals are expected to avoid kin‐competition by dispersing long distances, but maintain the benefits of cooperation by dispersing in small groups.     

LIMITED DISPERSAL, BUDDING DISPERSAL, AND COOPERATION: AN EXPERIMENTAL STUDY

Numerous theoretical studies have investigated how limited dispersal may provide an explanation for the evolution of cooperation, by leading to interactions between relatives. However, despite considerable theoretical attention, there has been a lack of empirical tests. In this article, we test how patterns of dispersal influence the evolution of cooperation, using iron-scavenging in the bacterium Pseudomonas aeruginosa as our cooperative trait. We found that relatively limited dispersal does not favor cooperation. The reason for this is that although limited dispersal increases the relatedness between interacting individuals, it also leads to increased local competition for resources between relatives. This result supports Taylor's prediction that in the simplest possible scenario, the effects of increased relatedness and local competition exactly cancel out. In contrast, we show that one way for cooperation to be favored is if individuals disperse in groups (budding dispersal), because this maintains high relatedness while reducing local competition between relatives (relatively global competition).     

Top-down and bottom-up effects on the spatiotemporal dynamics of cereal aphids: testing scaling theory for local density

The relationship between density and area depends on local growth rates and the area-dependence of migration rates. These rates vary among taxa due to dispersal behaviour, plot productivity and natural enemy impact. Previous studies in aphids suggest that aphid densities are highest in patches of intermediate sizes, and lower in small and large patches. The suggested mechanism causing these patterns is that the dispersal behaviour in aphids creates a mixture of area- and perimeter-dependent migration rates. In this paper, we used these predictions to examine the additional consequences of nutrient availability and natural enemies on the density-area relationship. The derived predictions were compared to data from a system with three aphid species, a set of aphid parasitoids and generalist natural enemies, and at two levels of plant nutrient availability. We find that predictions from the model based only on dispersal and local growth agree with the temporal dynamics of density-area relationships for aphids in high nutrient patches. In patches with low nutrients, high parasitism rates appeared to cause a negative density-area relationship for aphids, thereby deviating from predictions driven by the aphids' dispersal behavior. Hence, the dispersal model with scale-dependent migration rates can provide a useful tool for understanding insect distribution in patch size gradients, but the relative importance of top-down effects can completely change with plot productivity.