The different limits of weak selection and the evolutionary dynamics of finite populations

Evolutionary theory often resorts to weak selection, where different individuals have very similar fitness. Here, we relate two ways to introduce weak selection. The first considers evolutionary games described by payoff matrices with similar entries. This approach has recently attracted a lot of interest in the context of evolutionary game dynamics in finite populations. The second way to introduce weak selection is based on small distances in phenotype space and is a standard approach in kin-selection theory. Whereas both frameworks are interchangeable for constant fitness, frequency-dependent selection shows significant differences between them. We point out the difference between both limits of weak selection and discuss the condition under which the differences vanish. It turns out that this condition is fulfilled by the popular parametrization of the prisoner's dilemma in benefits and costs. However, for general payoff matrices differences between the two frameworks prevail.     

Evolutionary dynamics of continuous strategy games on graphs and social networks under weak selection

Understanding the emergence of cooperation among selfish individuals has been a long-standing puzzle, which has been studied by a variety of game models. Most previous studies presumed that interactions between individuals are discrete, but it seems unrealistic in real systems. Recently, there are increasing interests in studying game models with a continuous strategy space. Existing research work on continuous strategy games mainly focuses on well-mixed populations. Especially, little theoretical work has been conducted on their evolutionary dynamics in a structured population. In the previous work (Zhong et al., BioSystems, 2012), we showed that under strong selection, continuous and discrete strategies have significantly different equilibrium and game dynamics in spatially structured populations. In this paper, we further study evolutionary dynamics of continuous strategy games under weak selection in structured populations. By using the fixation probability based stochastic dynamics, we derive exact conditions of natural selection favoring cooperation for the death–birth updating scheme. We also present a network gain decomposition of the game equilibrium, which might provide a new view of the network reciprocity in a quantitative way. Finally, we make a detailed comparison between games using discrete and continuous strategies. As compared to the former, we find that for the latter (i) the same selection conditions are derived for the general 2 × 2 game; especially, the rule b/c > k in a simplified Prisoner's Dilemma is valid as well; however, (ii) for a coordination game, interestingly, the risk-dominant strategy is disfavored. Numerical simulations have also been conducted to validate our results.     

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 probability of fixation of a single mutant in an exchangeable selection model

The Cannings exchangeable model for a finite population in discrete time is extended to incorporate selection. The probability of fixation of a mutant type is studied under the assumption of weak selection. An exact formula for the derivative of this probability with respect to the intensity of selection is deduced, and developed in the case of a single mutant. This formula is expressed in terms of mean coalescence times under neutrality assuming that the coefficient of selection for the mutant type has a derivative with respect to the intensity of selection that takes a polynomial form with respect to the frequency of the mutant type. An approximation is obtained in the case where this derivative is a continuous function of the mutant frequency and the population size is large. This approximation is consistent with a diffusion approximation under moment conditions on the number of descendants of a single individual in one time step. Applications to evolutionary game theory in finite populations are presented.      

A model of weak selection in the infinite alleles framework

Ewens (1972) proposed a model in the infinite allele framework for populations with neutrality of all alleles at a particular locus. This paper proposes a generalisation of Ewens' result for situations where there is a form of weak selection. The models considered here are continuous time, discrete state space Markov processes.     

A branching process, its application in biology: influence of demographic parameters on the social structure in mammal groups

Branching processes are widely used in biology. This theoretical tool is used in cell dynamics, epidemics and population dynamics. In population dynamics, branching processes are mainly used to access extinction probabilities of populations, groups or families, with the Galton-Watson branching process. Many mammal species live in socially-structured groups, and the smallest units of these groups are lineages (or families) of kin-related individuals. In many primate species, these lineages are matrilines, as females remain in their natal groups most of the time, whereas males generally disperse. Lineage parameters, such as numbers of matrilines, size of each matriline and average degree of relatedness, could strongly influence the genetic composition of groups. Evidence indicates that division along matrilines could induce substantial differentiation among fission groups. Here, we develop a novel mathematical model based on the branching process theory describing demographic dynamics of groups. The main result of this model is an explicit analytical expression of the joint distribution of numbers of lineages and sizes of socially-structured groups. We investigated the influence of parameters such as natality and mortality on the outcome of the process, including extinction probability. Finally, we discuss this theoretical result with respect to biological significance.     

Diffusion approximations for one-locus multi-allele kin selection,mutation and random drift in group-structured populations: a unifying approach to selection models in population genetics

Diffusion approximations are ascertained from a two-time-scale argument in the case of a group-structured diploid population with scaled viability parameters depending on the individual genotype and the group type at a single multi-allelic locus under recurrent mutation, and applied to the case of random pairwise interactions within groups. The main step consists in proving global and uniform convergence of the distribution of the group types in an infinite population in the absence of selection and mutation, using a coalescent approach. An inclusive fitness formulation with coefficient of relatedness between a focal individual J affecting the reproductive success of an individual I, defined as the expected fraction of genes in I that are identical by descent to one or more genes in J in a neutral infinite population, given that J is allozygous or autozygous, yields the correct selection drift functions. These are analogous to the selection drift functions obtained with pure viability selection in a population with inbreeding. They give the changes of the allele frequencies in an infinite population without mutation that correspond to the replicator equation with fitness matrix expressed as a linear combination of a symmetric matrix for allozygous individuals and a rank-one matrix for autozygous individuals. In the case of no inbreeding, the mean inclusive fitness is a strict Lyapunov function with respect to this deterministic dynamics. Connections are made between dispersal with exact replacement (proportional dispersal), uniform dispersal, and local extinction and recolonization. The timing of dispersal (before or after selection, before or after mating) is shown to have an effect on group competition and the effective population size. In memory of Sam Karlin.     

On the survival probability of a slightly advantageous mutant gene in a multitype population: A multidimensional branching process model

The estimated survival probability of a slightly supercritical Galton-Watson process is generalized to a multitype branching process. The result is used to estimate the probability of initial success of a mutant gene whose effect on the individual carrier depends on the carrier's sex, class, etc. The probability of initial success is also estimated in a case where the effect of the mutation is manifested in terms of the distribution of types within one's progeny, e.g. in a case of a change in the sex ratio.