Selection for mutational robustness in finite populations

We investigate the evolutionary dynamics of a finite population of RNA sequences replicating on a neutral network. Despite the lack of differential fitness between viable sequences, we observe typical properties of adaptive evolution, such as increase of mean fitness over time and punctuated-equilibrium transitions, after initial mutation-selection balance has been reached. We find that a product of population size and mutation rate of approximately 30 or larger is sufficient to generate selection pressure for mutational robustness, even if the population size is orders of magnitude smaller than the neutral network on which the population resides. Our results show that quasispecies effects and neutral drift can occur concurrently, and that the relative importance of each is determined by the product of population size and mutation rate.     

Language dynamics in finite populations

Any mechanism of language acquisition can only learn a restricted set of grammars. The human brain contains a mechanism for language acquisition which can learn a restricted set of grammars. The theory of this restricted set is universal grammar (UG). UG has to be sufficiently specific to induce linguistic coherence in a population. This phenomenon is known as "coherence threshold". Previously, we have calculated the coherence threshold for deterministic dynamics and infinitely large populations. Here, we extend the framework to stochastic processes and finite populations. If there is selection for communicative function (selective language dynamics), then the analytic results for infinite populations are excellent approximations for finite populations; as expected, finite populations need a slightly higher accuracy of language acquisition to maintain coherence. If there is no selection for communicative function (neutral language dynamics), then linguistic coherence is only possible for finite populations.     

On selection in finite populations

Journal of Mathematical Biology - Two major forces shaping evolution are drift and selection. The standard models of neutral drift—the Wright–Fisher (WF) and Moran processes—can...     

Fitness-dependent mutation rates in finite populations

Mutation rate may be condition dependent, whereby individuals in poor condition, perhaps from high mutation load, have higher mutation rates than individuals in good condition. Agrawal (J. Evol. Biol.15, 2002, 1004) explored the basic properties of fitness-dependent mutation rate (FDMR) in infinite populations and reported some heuristic results for finite populations. The key parameter governing how infinite populations evolve under FDMR is the curvature (k) of the relationship between fitness and mutation rate. We extend Agrawal's analysis to finite populations and consider dominance and epistasis. In finite populations, the probability of long-term existence depends on k. In sexual populations, positive curvature leads to low equilibrium mutation rate, whereas negative curvature results in high mutation rate. In asexual populations, negative curvature results in rapid extinction via 'mutational meltdown', whereas positive curvature sometimes allows persistence. We speculate that fitness-dependent mutation rate may provide the conditions for genetic architecture to diverge between sexual and asexual taxa.     

Evolutionary game dynamics in finite populations

We introduce a model of stochastic evolutionary game dynamics in finite populations which is similar to the familiar replicator dynamics for infinite populations. Our focus is on the conditions for selection favoring the invasion and/or fixation of new phenotypes. For infinite populations, there are three generic selection scenarios describing evolutionary game dynamics among two strategies. For finite populations, there are eight selection scenarios. For a fixed payoff matrix a number of these scenarios can occur for different population sizes. We discuss several examples with unexpected behavior.     

Multistage index selection in finite populations

Multistage selection with fixed proportions and selection indices based on covariates of the target variable is studied. Assuming a multivariate normal distribution before the selection, expressions are presented for the expectation and the variance of the target variable in the retained subpopulation. As the numerical evaluation for finite populations requires lengthy computations, some approximations using methods for infinite populations are proposed. Numerical illustrations are given for selections in up to three stages.     

Concepts of substitutional load in finite populations

In defining the substitutional load in a finite population, one can use a definition relating to the total load suffered until substitution of a favoured allele occurs (this will normally involve a contribution to the load from a number of cases where the favoured allele is lost), or one can use a definition using the load only for a single substitution process in which successful substitution of the favoured allele does occur. It is not evident at the moment which definition will be the more useful and meaningful in discussing how load arguments might limit the rate of selectively controlled gene substitutions. In this paper relatively simple methods are used to find formulae for both loads (the formula for the first confirming a previous result of Kimura and Maruyama, 1969). Numerical comparisons of the two loads are given, together with the results of a Monte Carlo experiment confirming the theoretical value for the second load. The latter theoretical value is found by considering a conditional diffusion process, the condition being that the favoured allele is eventually fixed. This diffusion process can be used, amongst other things, to arrive at the formulae for the mean and variance of the conditional fixation time found by other methods by  and .     

A theoretical basis for measures of kin selection in subdivided populations: finite populations and localized dispersal

The analysis of kin selection in subdivided populations has been hampered by the lack of well‐defined measures of genealogical relatedness in the presence of localized dispersal. Furthermore, the usual arguments underlying the definition of game‐theoretical measures of inclusive fitness are not exact under localized dispersal. We define such measures to give the first‐order effects of selection on the probability of fixation of an allele. The derived measures of kin selection and relatedness are valid in finite populations and under localized dispersal. For the infinite island model, the resulting measure of kin selection is equivalent to a previously used measure. In other cases its definition is based on definitions of relatedness which are different from the usual ones. To illustrate the approach, we reanalyse a model with localized dispersal. We consider sex ratio evolution under sex‐specific dispersal behaviour, and the results confirm the earlier conclusion that the sex ratio is biased towards the sex with the dispersal rate closer to the optimal dispersal rate in the absence of sex‐specific dispersal behaviour.     

Stability properties of underdominance in finite subdivided populations

IN ISOLATED populations underdominance leads to bistable evolutionary dynamics: below a certain mutant allele frequency the wildtype succeeds. Above this point, the potentially underdominant mutant allele fixes. In subdivided populations with gene flow there can be stable states with coexistence of wildtypes and mutants: polymorphism can be maintained because of a migration-selection equilibrium, i.e., selection against rare recent immigrant alleles that tend to be heterozygous. We focus on the stochastic evolutionary dynamics of systems where demographic fluctuations in the coupled populations are the main source of internal noise. We discuss the influence of fitness, migration rate, and the relative sizes of two interacting populations on the mean extinction times of a group of potentially underdominant mutant alleles. We classify realistic initial conditions according to their impact on the stochastic extinction process. Even in small populations, where demographic fluctuations are large, stability properties predicted from deterministic dynamics show remarkable robustness. Fixation of the mutant allele becomes unlikely but the time to its extinction can be long.     

Linkage disequilibrium under genetic hitchhiking in finite populations

The model of genetic hitchhiking predicts a reduction in sequence diversity at a neutral locus closely linked to a beneficial allele. In addition, it has been shown that the same process results in a specific pattern of correlations (linkage disequilibrium) between neutral polymorphisms along the chromosome at the time of fixation of the beneficial allele. During the hitchhiking event, linkage disequilibrium on either side of the beneficial allele is built up whereas it is destroyed across the selected site. We derive explicit formulas for the expectation of the covariance measure D and standardized linkage disequilibrium sigma 2D between a pair of polymorphic sites. For our analysis we use the approximation of a star-like genealogy at the selected site. The resulting expressions are approximately correct in the limit of large selection coefficients. Using simulations we show that the resulting pattern of linkage disequilibrium is quickly-i.e., in <0.1N generations-destroyed after the fixation of the beneficial allele for moderately distant neutral loci, where N is the diploid population size.     

The distribution of lethal allelism in finite populations

It is shown that, under certain conditions, the distribution of lethal allelism in equilibrium populations is independent of the fitness of heterozygotes. For these conditions, an expression for the over-generation variance of allelism is given.     

Random mating and random union of gametes models for finite dioecious populations

This paper compares a random mating model for independent trials with two random union of gametes models in the case when the population is finite and dioecious. The first random union of gametes model is examined as a two-locus dioecious model with unequal recombination values in each sex, which is a generalization of previous work. The second model can be of some biological use under certain circumstances, and it is easier to analyse (when appropriate) than the first. It is examined with and without mutation to obtain both old and new results about the fixation probabilities (in the multilocus case) and the related rates in the autosomal locus case.     

Public goods games with reward in finite populations

Public goods games paraphrase the problem of cooperation in game theoretical terms. Cooperators contribute to a public good and thereby increase the welfare of others at a cost to themselves. Defectors consume the public good but do not pay its cost and therefore outperform cooperators. Hence, according to genetic or cultural evolution, defectors should be favored and the public good disappear – despite the fact that groups of cooperators are better off than groups of defectors. The maximization of short term individual profits causes the demise of the common resource to the detriment of all. This outcome can be averted by introducing incentives to cooperate. Negative incentives based on the punishment of defectors efficiently stabilize cooperation once established but cannot initiate cooperation. Here we consider the complementary case of positive incentives created by allowing individuals to reward those that contribute to the public good. The finite-population stochastic dynamics of the public goods game with reward demonstrate that reward initiates cooperation by providing an escape hatch out of states of mutual defection. However, in contrast to punishment, reward is unable to stabilize cooperation but, instead, gives rise to a persistent minority of cooperators.