Deme Size Research Articles

Evolutionary games in deme structured, finite populations

We describe a fairly general model for the evolutionary dynamics in a sub-divided (or deme structured) population with migration and mutation. The number and size of demes are finite and fixed. The fitness of each individual is determined by pairwise interactions with other members of the same deme. The dynamics within demes can be modeled according to a broad range of evolutionary processes. With a probability proportional to fitness, individuals migrate to another deme. Mutations occur randomly. In the limit of few migrations and even rarer mutations we derive a simple analytic condition for selection to favor one strategic type over another. In particular, we show that the Pareto efficient type is favored when competition within demes is sufficiently weak. We then apply the general results to the prisoner's dilemma game and discuss selected dynamics and the conditions for cooperation to prevail.

Demographic and genetic factors shaping contemporary metapopulation effective size and its empirical estimation in salmonid fish.

The preservation of biodiversity requires an understanding of the maintenance of its components, including genetic diversity. Effective population size determines the amount of genetic variance maintained in populations, but its estimation can be complex, especially when populations are interconnected in a metapopulation. Theory predicts that the effective size of a metapopulation (meta-N(e)) can be decreased or increased by population subdivision, but little empirical work has evaluated these predictions. Here, we use neutral genetic markers and simulations to estimate the effective size of a putative metapopulation in Atlantic salmon (Salmo salar). For a weakly structured set of rivers, we find that meta-N(e) is similar to the sum of local deme sizes, whereas higher genetic differentiation among demes dramatically reduces meta-N(e) estimates. Interdemic demographic processes, such as asymmetrical gene flow, may explain this pattern. However, simulations also suggest that unrecognized population subdivision can also introduce downward bias into empirical estimation, emphasizing the importance of identifying the proper scale of distinct demographic and genetic processes. Under natural patterns of connectivity, evolutionary potential may generally be maintained at higher levels than the local population, with implications for conservation given ongoing species declines and habitat fragmentation.

Weighted local sharing and local clearing for multimodal optimisation

Local sharing is a method designed for efficient multimodal optimisation that combines fitness sharing with spatially structured populations and elitist replacement. In local sharing, the bias towards sharing and the influence of spatial structure is controlled by the deme (neighbourhood) size. This introduces an undesirable trade-off; to maximise the sharing effect large deme sizes must be used, but the opposite must be true if one wishes to maximise the influence of spatial population structure. This paper introduces two modifications to the local sharing method. The first alters local sharing so that parent selection and fitness sharing operate at two different spatial levels; parent selection is performed within small demes, while the effect of fitness sharing is weighted according to the distance between individuals in the entire population structure. The second method replaces fitness sharing within demes with clearing to produce a method that we call local clearing. The proposed methods, as tested on several benchmark problems, demonstrate a level of efficiency that surpasses that of traditional fitness sharing and standard local sharing. Additionally, they offer a level of parameter robustness that surpasses other elitist niching methods, such as clearing. Through analysis of the local clearing method, we show that this parameter robustness is a result of the isolated nature of the demes in a spatially structured population being able to independently concentrate on subsets of the desired optima in a fitness landscape.

Local Extinction and Recolonization, Species Effective Population Size, and Modern Human Origins

A primary objection from a population genetics perspective to a multiregional model of modern human origins is that the model posits a large census size, whereas genetic data suggest a small effective population size. The relationship between census size and effective size is complex, but arguments based on an island model of migration show that if the effective population size reflects the number of breeding individuals and the effects of population subdivision, then an effective population size of 10,000 is inconsistent with the census size of 500,000 to 1,000,000 that has been suggested by archeological evidence. However, these models have ignored the effects of population extinction and recolonization, which increase the expected variance among demes and reduce the inbreeding effective population size. Using models developed for population extinction and recolonization, we show that a large census size consistent with the multiregional model can be reconciled with an effective population size of 10,000, but genetic variation among demes must be high, reflecting low interdeme migration rates and a colonization process that involves a small number of colonists or kin-structured colonization. Ethnographic and archeological evidence is insufficient to determine whether such demographic conditions existed among Pleistocene human populations, and further work needs to be done. More realistic models that incorporate isolation by distance and heterogeneity in extinction rates and effective deme sizes also need to be developed. However, if true, a process of population extinction and recolonization has interesting implications for human demographic history.

Joint Effects of Inbreeding and Local Adaptation on the Evolution of Genetic Load after Fragmentation

Disruption of gene flow among demes after landscape fragmentation can facilitate local adaptation but increase the effect of genetic drift and inbreeding. The joint effects of these conflicting forces on the mean fitness of individuals in a population are unknown. Through simulations, we explored the effect of increased isolation on the evolution of genetic load over the short and long term when fitness depends in part on local adaptation. We ignored genetic effects on demography. We modeled complex genomes, where a subset of the loci were under divergent selection in different localities. When a fraction of the loci were under heterogeneous selection, isolation increased mean fitness in larger demes made up of hundreds of individuals because of improved local adaptation. In smaller demes of tens of individuals, increased isolation improved local adaptation very little and reduced overall fitness. Short-term improvement of mean fitness after fragmentation may not be indicative of the long-term evolution of fitness. Whatever the deme size and potential for local adaptation, migration of one or two individuals per generation minimized the genetic load in general. The slow dynamics of mean fitness following fragmentation suggests that conservation measures should be implemented before the consequences of isolation on the genetic load become of concern.

A STRUCTURED COALESCENT PROCESS FOR SEASONALLY FLUCTUATING POPULATIONS

Many short-lived organisms pass through several generations during favorable growing seasons, separated by inhospitable periods during which only small hibernating or estivating refugia remain. This induces pronounced seasonal fluctuations in population size and metapopulation structure. The first generations in the growing season will be characterized by small, relatively isolated demes whereas the later generations will experience larger deme sizes with more extensive gene flow. Fluctuations of this sort can induce changes in the amount of genetic variation in early season samples compared to late season samples, a classical example being the observations of seasonal variation in allelism in New England Drosophila populations by P.T. Ives. In this article, we study the properties of a structured coalescent process under seasonal fluctuations using numerical analysis of exact state equations, analytical approximations that rely on a separation of timescales between intrademic versus interdemic processes, and individual-based simulations. We show that although an increase in genetic variation during each favorable growing season is observed, it is not as pronounced as in the empirical observations. This suggests that some of the temporal patterns of variation seen by Ives may be due to selection against deleterious lethals rather than neutral processes.

Local Extinction and Recolonization, Species Effective Population Size, and Modern Human Origins

A primary objection from a population genetics perspective to a multiregional model of modern human origins is that the model posits a large census size, whereas genetic data suggest a small effective population size. The relationship between census size and effective size is complex, but arguments based on an island model of migration show that if the effective population size reflects the number of breeding individuals and the effects of population subdivision, then an effective population size of 10,000 is inconsistent with the census size of 500,000 to 1,000,000 that has been suggested by archeological evidence. However, these models have ignored the effects of population extinction and recolonization, which increase the expected variance among demes and reduce the inbreeding effective population size. Using models developed for population extinction and recolonization, we show that a large census size consistent with the multiregional model can be reconciled with an effective population size of 10,000, but genetic variation among demes must be high, reflecting low interdeme migration rates and a colonization process that involves a small number of colonists or kin-structured colonization. Ethnographic and archeological evidence is insufficient to determine whether such demographic conditions existed among Pleistocene human populations, and further work needs to be done. More realistic models that incorporate isolation by distance and heterogeneity in extinction rates and effective deme sizes also need to be developed. However, if true, a process of population extinction and recolonization has interesting implications for human demographic history.

Genetic diversity in a crop metapopulation

The need to protect crop genetic resources has sparked a growing interest in the genetic diversity maintained in traditional farming systems worldwide. Although traditional seed management has been proposed as an important determinant of genetic diversity and structure in crops, no models exist that can adequately describe the genetic effects of seed management. We present a metapopulation model that accounts for several features unique to managed crop populations. Using traditional maize agriculture as an example, we develop a coalescence-based model of a crop metapopulation undergoing pollen and seed flow as well as seed replacement. In contrast to metapopulation work on natural systems, we model seed migration as episodic and originating from a single source per population rather than as a constant immigration from the entire metapopulation. We find that the correlated origin of migrants leads to surprising results, including a loss of invariance of within-deme diversity and a parabolic relationship between F(ST) and migration quantity. In contrast, the effects of migration frequency on diversity and structure are more similar to classical predictions, suggesting that seed migration in managed crop populations cannot be described by a single parameter. In addition to migration, we investigate the effects of deme size and extinction rates on genetic structure, and show that high levels of pollen migration may mask the effects of seed management on structure. Our results highlight the importance of analytically evaluating the effects of deviations from classical metapopulation models, especially in systems for which data are available to estimate specific model parameters.

Environmental variation, fluctuating selection and genetic drift in subdivided populations

Although there have many studies of the population genetical consequences of environmental variation, little is known about the combined effects of genetic drift and fluctuating selection in structured populations. Here we use diffusion theory to investigate the effects of temporally and spatially varying selection on a population of haploid individuals subdivided into a large number of demes. Using a perturbation method for processes with multiple time scales, we show that as the number of demes tends to infinity, the overall frequency converges to a diffusion process that is also the diffusion approximation for a finite, panmictic population subject to temporally fluctuating selection. We find that the coefficients of this process have a complicated dependence on deme size and migration rate, and that changes in these demographic parameters can determine both the balance between the dispersive and stabilizing effects of environmental variation and whether selection favors alleles with lower or higher fitness variance.

Evolutionary game dynamics in a finite asymmetric two-deme population and emergence of cooperation

We consider evolutionary game dynamics in a finite population subdivided into two demes with both unequal deme sizes and different migration rates. Assuming viability differences in the population according to a linear game within each deme as a result of pairwise interactions, we specify conditions for weak selection favoring a mutant strategy to go to fixation, under the structured-coalescent assumptions, and their connections with evolutionary stability concepts. In the framework of the Iterated Prisoner's Dilemma with strategy ‘tit-for-tat’ as mutant strategy and ‘always defect’ as resident strategy, we deduce a condition under which the emergence of cooperation is favored by selection, when the game matrix is the same in both demes. We show how this condition extends the one-third law for a panmictic population and when an asymmetry in the spatial structure of a two-deme population facilitates the emergence of the cooperative tit-for-tat strategy in comparison with both its symmetric and panmictic population structure counterparts. We find that the condition is less stringent in the asymmetric scenario versus the symmetric scenario if both the fraction of the population in the deme where the mutant was initially introduced, and the expected proportion of migrant offspring in this deme among all migrant offspring after population regulation, are smaller than, or equal to, 1 2 , provided they are not too small. On the other hand, the condition is less stringent than the one-third law, which holds in the panmictic case, if the latter proportion remains not too close to 1.

Multilocus models in the infinite island model of population structure

Different methods have been developed to consider the effects of statistical associations among genes that arise in population genetics models: kin selection models deal with associations among genes present in different interacting individuals, while multilocus models deal with associations among genes at different loci. It was pointed out recently that these two types of models are very similar in essence. In this paper, we present a method to construct multilocus models in the infinite island model of population structure (where deme size may be arbitrarily small). This method allows one to compute recursions on allele frequencies, and different types of genetic associations (including associations between different individuals from the same deme), and incorporates selection. Recursions can be simplified using quasi-equilibrium approximations; however, we show that quasi-equilibrium calculations for associations that are different from zero under neutrality must include a term that has not been previously considered. The method is illustrated using simple examples.

The spread of incompatibility-inducing parasites in sub-divided host populations.

BackgroundMaternally transmitted symbionts have evolved a variety of ways to promote their spread through host populations. One strategy is to hamper the reproduction of uninfected females by a mechanism called cytoplasmic incompatibility (CI). CI occurs in crosses between infected males and uninfected females and leads to partial to near-complete infertility. CI-infections are under positive frequency-dependent selection and require genetic drift to overcome the range of low frequencies where they are counter-selected. Given the importance of drift, population sub-division would be expected to facilitate the spread of CI. Nevertheless, a previous model concluded that variance in infection between competing groups of breeding individuals impedes the spread of CI.ResultsIn this paper we derive a model on the spread of CI-infections in populations composed of demes linked by restricted migration. Our model shows that population sub-division facilitates the invasion of CI. While host philopatry (low migration) favours the spread of infection, deme size has a non-monotonous effect, with CI-invasion being most likely at intermediate deme size. Individual-based simulations confirm these predictions and show that high levels of local drift speed up invasion but prevent high levels of prevalence across the entire population. Additional simulations with sex-specific migration rates further show that low migration rates of both sexes are required to facilitate the spread of CI.ConclusionOur analyses show that population structure facilitates the invasion of CI-infections. Since some level of sub-division is likely to occur in most natural populations, our results help to explain the high incidence of CI-infections across species of arthropods. Furthermore, our work has important implications for the use of CI-systems in order to genetically modify natural populations of disease vectors.

Gene surfing in expanding populations

Large scale genomic surveys are partly motivated by the idea that the neutral genetic variation of a population may be used to reconstruct its migration history. However, our ability to trace back the colonization pathways of a species from their genetic footprints is limited by our understanding of the genetic consequences of a range expansion. Here, we study, by means of simulations and analytical methods, the neutral dynamics of gene frequencies in an asexual population undergoing a continual range expansion in one dimension. During such a colonization period, lineages can fix at the wave front by means of a “surfing” mechanism [Edmonds, C.A., Lillie, A.S., Cavalli-Sforza, L.L., 2004. Mutations arising in the wave front of an expanding population. Proc. Natl. Acad. Sci. 101, 975–979]. We quantify this phenomenon in terms of (i) the spatial distribution of lineages that reach fixation and, closely related, (ii) the continual loss of genetic diversity (heterozygosity) at the wave front, characterizing the approach to fixation. Our stochastic simulations show that an effective population size can be assigned to the wave that controls the (observable) gradient in heterozygosity left behind the colonization process. This effective population size is markedly higher in the presence of cooperation between individuals (“ pushed waves”) than when individuals proliferate independently (“ pulled waves”), and increases only sub-linearly with deme size. To explain these and other findings, we develop a versatile analytical approach, based on the physics of reaction–diffusion systems, that yields simple predictions for any deterministic population dynamics. Our analytical theory compares well with the simulation results for pushed waves, but is less accurate in the case of pulled waves when stochastic fluctuations in the tip of the wave are important.

Spatial and Temporal Genetic Structure of Anopheles arabiensis in Southern Zambia over Consecutive Wet and Drought Years

No studies have addressed the spatial complexity of Anopheles arabiensis populations in Zambia or the effects of drought on the genetic structure of this species. We genotyped approximately 420 An. arabiensis at 12 microsatellite loci representing 18 collections from the Southern Province of Zambia. Collections spanned three transmission seasons and covered a wet year-drought year-wet year cycle. Anopheles arabiensis within the 2,000 km(2) of the Macha study region were panmictic, with high gene flow between Macha and Namwala, Zambia, which are 80 km apart. There was little evidence for genetic structuring among years, with no significant shifts in allele frequency distributions or observed heterozygosity, and no evidence for a genetic bottleneck despite a drastic reduction in mosquito numbers during the drought year. Anopheles arabiensis in southern Zambia has a large deme size, and the regional genetic structure of this species was little affected by an extended drought period.

Fixation probability for a beneficial allele and a mutant strategy in a linear game under weak selection in a finite island model

The effect of population structure on the probability of fixation of a newly introduced mutant under weak selection is studied using a coalescent approach. Wright's island model in a framework of a finite number of demes is assumed and two selection regimes are considered: a beneficial allele model and a linear game among offspring. A first-order approximation of the fixation probability for a single mutant with respect to the intensity of selection is deduced. The approximation requires the calculation of expected coalescence times, under neutrality, for lineages starting from two or three sampled individuals. The results are obtained in a general setting without assumptions on the number of demes, the deme size or the migration rate, which allows for simultaneous coalescence or migration events in the genealogy of the sampled individuals. Comparisons are made with limit cases as the deme size or the number of demes goes to infinity or the migration rate goes to zero for which a diffusion approximation approach is possible. Conditions for selection to favor a mutant strategy replacing a resident strategy in the context of a linear game in a finite island population are addressed.

The Role of Life Cycle and Migration in Selection for Variance in Offspring Number

For two genotypes that have the same mean number of offspring but differ in the variance in offspring number, naturalselection will favor the genotype with lower variance. In such cases, the average growth rate is not sufficient as a measure of fitness or as a predictor of fixation probability. However, the effect of variance in offspring number on the fixationprobability of mutant strategies has been calculated under several scenarios with the general conclusion that variance in offspring number reduces fitness in proportion to the inverse of the population size [Gillespie, J., Genetics 76:601-606, 1974; Proulx, S.R., Theor. Popul. Biol. 58:33-47, 2000]. This relationship becomes more complicated under a metapopulation scenario where the "effective" population size depends on migration rate, population structure, and lifecycle. It is shown that in a life cycle where reproduction and migration (the birth-migration-regulation life cycle, or BMR)occur prior to density regulation within every deme, the fitness of a strategy depends on migration rate. When migration rates are near zero, the fitness of the strategy is determined by the size of individual demes, so that the strategy favoredin small populations tends to be fixed. As migration rate increases and approaches panmixis between demes, the fitness ofa reproductive strategy approaches what its value would be in a single, panmictic deme with a population size correspondingtothe census size of the metapopulation. Interestingly, when the life cycle is characterized by having density regulation in each deme prior to migration (the BRM life cycle) the fixation probability of a strategy is independent of migration rate. These results are found to be qualitatively consistent with the individual-based simulation results in Shpak [Theor. Biosci.124:65-85, 2005].

An evolutionary analysis of the relationship between spite and altruism

We investigate the selective pressures on a social trait when evolution occurs in a population of constant size. We show that any social trait that is spiteful simultaneously qualifies as altruistic. In other words, any trait that reduces the fitness of less related individuals necessarily increases that of related ones. Our analysis demonstrates that the distinction between "Hamiltonian spite" and "Wilsonian spite" is not justified on the basis of fitness effects. We illustrate this general result with an explicit model for the evolution of a social act that reduces the recipient's survival ("harming trait"). This model shows that the evolution of harming is favoured if local demes are of small size and migration is low (philopatry). Further, deme size and migration rate determine whether harming evolves as a selfish strategy by increasing the fitness of the actor, or as a spiteful/altruistic strategy through its positive effect on the fitness of close kin.

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.

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 Fate of Mutations Surfing on the Wave of a Range Expansion

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