Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Research article

The impact of random frequency-dependent mutations on the average population fitness

Weini Huang, Benjamin Werner and Arne Traulsen*

Author affiliations

Evolutionary Theory Group, Max Planck Institute for Evolutionary Biology Plön, 24306 Plön, August-Thienemann-Straße 2, Germany

For all author emails, please log on.

Citation and License

BMC Evolutionary Biology 2012, 12:160  doi:10.1186/1471-2148-12-160


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/12/160


Received:30 March 2012
Accepted:22 August 2012
Published:30 August 2012

© 2012 Huang et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

In addition to selection, the process of evolution is accompanied by stochastic effects, such as changing environmental conditions, genetic drift and mutations. Commonly it is believed that without genetic drift, advantageous mutations quickly fixate in a halpoid population due to strong selection and lead to a continuous increase of the average fitness. This conclusion is based on the assumption of constant fitness. However, for frequency dependent fitness, where the fitness of an individual depends on the interactions with other individuals in the population, this does not hold.

Results

We propose a mathematical model that allows to understand the consequences of random frequency dependent mutations on the dynamics of an infinite large population. The frequencies of different types change according to the replicator equations and the fitness of a mutant is random and frequency dependent. To capture the interactions of different types, we employ a payoff matrix of variable size and thus are able to accommodate an arbitrary number of mutations. We assume that at most one mutant type arises at a time. The payoff entries to describe the mutant type are random variables obeying a probability distribution which is related to the fitness of the parent type.

Conclusions

We show that a random mutant can decrease the average fitness under frequency dependent selection, based on analytical results for two types and simulations for n types. Interestingly, in the case of at most two types the probabilities to increase or decrease the average fitness are independent of the concrete probability density function. Instead, they only depend on the probability that the payoff entries of the mutant are larger than the payoff entries of the parent type.

Background

Mutations provide a continuous source of variation in natural populations, on which natural selection can act. When fitness is assumed to be constant, only those mutations with higher fitness values will be fixed in a haploid population under strong selection and negligible random drift. Thus, the average fitness of the population would monotonically increase in evolutionary time. There have been numerous hypotheses why this is not what is observed in nature: for instance, environmental changes require new adaptions [1,2] or coevolution can imply continuous adaptation without increasing the average fitness [3-5]. However, these are not aspects that we intend to include here. Instead, we focus on a haploid population in a constant environment, and explore frequency dependent fitness, which can be described by evolutionary game theory [6-11]. In this framework, the fitness of a type depends on the frequencies of other types of individuals in the population. We address the very general question of how the average fitness changes when it is driven by random mutations under frequency dependent selection.

The fitness effects of new mutations have gained significant attention both in experimental research and theoretical work [12,13]. In experiments, the distribution of fitness effects depends on several aspect of the experimental setup, e.g. how well adapted the organism is to the environment and whether only single mutants or also double mutants (mutants differing from the wild type by two mutations) are considered. Different shaped distributions were proposed to capture the fitness distributions of random mutants under constant selection [14-17]. The concrete shape of fitness distributions of spontaneous mutations varies between species and even within the same species on different parts of DNA [18]. Although no common conclusion on this has been obtained yet – and a universal fitness distribution may as well not exist – it is often possible to estimate some general properties, such as the proportion of advantageous mutations and the mean value of the fitness of the mutations [19,20].

The concept of random distributed and frequency dependent fitness of mutations can be addressed by evolutionary game theory [21], which considers evolutionary processes under frequency dependent selection [22]. In this framework, a population of interacting individuals is considered. In the simplest case of linear frequency dependence, the interactions of different types of individuals are captured by a payoff matrix for a game. Those types which are more successful in the game will have a higher reproduction rate. We introduce a payoff matrix with variable size to capture mutations and extinctions. The new payoff entries introduced by mutations are independently drawn from a probability distribution, which corresponds to the concept of randomly distributed fitness. By tracking the dynamics of the payoff matrix and the compositions of the population, we are able to investigate several aspects of an evolving system, such as the average fitness changes of the population, the impact of the fitness distribution on these changes and the expected level of diversity.

Results

Dynamics for populations with two types

Let us start with a population of a resident wild type (R) and a mutant type (M). Suppose the fitness of a wild type in a homogenous population is d. For constant selection, the fitness distribution of a mutant is simply a one dimensional distribution around d. For frequency dependent selection, the fitness of a mutant must be defined based on more than a single number. We can write it as an evolutionary game based on a 2×2 payoff matrix with three new payoff entries, a, b and c

When a mutant and a wild type interact, the mutant obtains fitness a, and the wild type obtains c. When a mutant meets another mutant, it obtains b. Following the concept of randomly distributed fitness of mutations, the entries a, b and c are defined as random variables. We assume that a, b and c independently follow the same probability distribution given by a probability density function f(x). While this is the simplest possibility, it may be more realistic to assume correlations between the payoff entries characterizing each type, i.e. between a and b as well as between c and d (see below, section Games with equal gains from switching). However, in the extreme case of a = b and c = d, this would recover the case of constant selection, so we expect that such correlations would lead to results intermediate between constant and frequency dependent selection. We discuss how this distribution affects the changes in the average fitness during the evolutionary process. It turns out, the probability <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M2">View MathML</a> that a payoff entry is larger than the fitness of the wild type (the parent type in the case of n types) d, is of particular interest and determines the change in the average fitness. Remarkably, all other aspects of the fitness distribution turn out to be irrelevant for this observable.

The dynamics of evolving populations shows stochastic fluctuations when selection is weak and when populations are small. In addition, stochasticity can arise based on environmental changes or stochastic effects due to mutations. As we are interested in the effects of frequency dependent selection, we only consider stochasticity arising from random frequency dependent mutations and use the replicator equations to model evolutionary dynamics. The frequency of a certain type changes deterministically according to the difference of its own fitness to the average fitness in the population.

Suppose x is the frequency of the mutant type and 1 − x the frequency of the wild type, respectively. We can define the fitness of the mutant type, W1, and the fitness of the wild type, W2, as

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M3">View MathML</a>

(1)

where a, b, c, and d are the entries in the payoff matrix. The average fitness of the population <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M4">View MathML</a> is given by

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M5">View MathML</a>

(2)

If the fitness of the mutant type is larger than the average fitness, its frequency will increase. If the fitness of the mutant type is below the average fitness, its frequency will decrease. We follow the usual assumption that the change of the frequency of the mutant type is given by the replicator equation [23-25]

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M6">View MathML</a>

(3)

The change of the wild type frequency follows immediately as <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M7">View MathML</a>. This dynamics is fully determined by the entries of the payoff matrix. Different constellations of the payoff entries cause different dynamical patterns. In the following, we discuss all generic cases of two-type interactions and how the average fitness of the population changes under the different situations.

First, we analyze the case where the mutant has higher fitness than the wild type for all frequencies x. This is the case for a > c and b > d. The wild type goes extinct and the mutant type will be fixed in the population. Thus, the average fitness <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M8">View MathML</a> in the new equilibrium x = 1 is given by the payoff entry of the mutant type interacting with itself, a. We are interested in the probability, that the fitness of the population is increased after the fixation of the mutant. This becomes a conditional probability of a > d given that a > c and b > d. Applying Bayes Rule, this can be expressed as

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M9">View MathML</a>

(4)

We assume that the random variables a, b and c are independently derived from the same probability distribution. Hence, b does not depend on a or on c. Thus, the probability of b > d is independent from the probability that a > d, which is used in Eq. (4). Since a and c are sampled from the same distribution, we have p(a > c) = 1/2 in the denominator. For the numerator, we have

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M10">View MathML</a>

(5)

where F(x) is the cumulative distribution function of a random variable with probability density function f(x). The probability that one of the new payoff entries a,b,c is greater than the wild type fitness d is <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M11">View MathML</a>. Using this expression in Eq. (5), we arrive at

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M12">View MathML</a>

(6)

Strikingly, this only depends on θ, and is independent of the concrete choice of the probability density function f(x). In population genetics, beneficial mutation rates are measured based on the concept of constant fitness, where the fitness of the mutant and the fitness of the wild type are both constant numbers. However, if we consider frequency dependent fitness, a new parameter is needed to represent the proportion of beneficial mutations. One option arising from our approach is to compare the payoff value of the mutant with the payoff value of the wild type when they are confronted by the same opponent. Since θ is the probability that the new payoff value of the mutant is larger than the wild type’s payoff d, it corresponds to the probability that a mutation is beneficial under the constant selection scenario. If θ can be measured, the probability that the average fitness is increased by a random mutant is independent of the payoff distribution according to Eq. (5). But different choices of probability density functions f(x) will result in different values of θ, thus leading to different probabilities to increase the average fitness.

Next, we assume that a mutant type occurs with lower fitness than the wild type. With frequency dependence, there are two situations for such a mutant type. The mutant type can either have lower fitness than the wild type for all frequencies, or it can have a lower fitness only for small frequencies. In both cases, the mutant will go extinct and the average fitness will remain unchanged, since a mutant type is supposed to arise with a small amount.

Finally, a mutant type could be initially advantageous compared to wild types, but turn to be disadvantageous when it has reached a certain frequency. This occurs for a < c and b > d. In this case neither the wild type nor the mutant type can take over the population, but there exists a mixed equilibrium consisting of the mutant type at a frequency <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M13">View MathML</a> and the wild type at a frequency 1 − x. In this coexistence equilibrium, the fitness of the wild type subpopulation is equal to the fitness of the mutant type subpopulation. The average fitness of the system in the equilibrium is given by

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M14">View MathML</a>

(7)

Again, we ask for the probability of having a coexistence game that increases the average fitness. This is the conditional probability that <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M15">View MathML</a> given that a < c and b > d, which can be written as

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M16">View MathML</a>

(8)

This is identical to Eq. (4) if one exchanges a c. Since a and c have the same distribution, we recover the result from Eq. (6),

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M17">View MathML</a>

(9)

In other words, the probability to increase fitness is the same in a coexistence game as in a game where the mutant dominates the wild type.

Let us now combine the results and consider the changes of the average fitness over all types of interactions. The probability to increase the fitness due to a new mutation is given by

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M18">View MathML</a>

(10)

In a similar manner, we can calculate the probability to decrease the average fitness due to a new mutation. When the mutant dominates the wild type, the average fitness may still decrease. This is exactly what happens in the Prisoner’s Dilemma [26,27]. Equivalently to the calculation above, we have

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M19">View MathML</a>

(11)

For the probability to decrease the average fitness in a coexistence game, we find

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M20">View MathML</a>

(12)

Thus, using a calculation similar to Eq. (10), the overall probability to decrease the average fitness is given by

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M21">View MathML</a>

(13)

Also the probability to maintain a constant average fitness can be calculated in this way. For continuous fitness distributions, there are no strictly neutral mutations. As the fitness of the wild type is a specific value of the continuous random variable, the probability of having a strict neutral mutation, the fitness of which is equal to the fitness of the wild type, is 0. Thus, the average fitness is only maintained when the mutant goes extinct, which occurs with probability

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M22">View MathML</a>

(14)

We discussed the changes of the average fitness in a two-type population under frequency dependent selection above. Under constant selection, the average fitness will increase with probability θ and decrease with probability 0. In the same way as for frequency dependent selection, it will remain constant with probability 1 − θ. Figure 1 illustrates these results and compares frequency dependent selection to constant selection for all values of θ. For frequency dependent selection, there is an intersection point θ, where the probability to increase the average fitness and to decrease the average fitness are equal. Using Eq. (10) and Eq. (13), this becomes <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M23">View MathML</a>, and we have <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M24">View MathML</a>. Small values of θ are typically considered to be of biological relevance. In this case, frequency dependent selection tends to decrease the average fitness: for <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M25">View MathML</a>, it is more likely that the average fitness of the population is decreased by a single random frequency dependent mutation; for <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M26">View MathML</a>, it is more likely that it is increased.

thumbnailFigure 1. Probability of increasing or decreasing the average fitness in the new equilibrium after one mutation event in an initial homogenous population. θ is the probability that a random payoff entry of the mutant, a, b or c is larger than wild type initial fitness d. Blue symbols and lines are simulation and analytical results under frequency dependent selection (average over 106 runs). Red lines are analytical results under constant selection. For constant selection, the average fitness either increases or is unchanged by a new mutation, where the fraction of mutants that increases fitness is determined by θ. However, under frequency dependent selection, the average fitness of the population in the new equilibrium after a mutation can also decrease. The probability to increase, decrease the average fitness or maintain the same average fitness, depends on θ, for <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M27">View MathML</a> the probability to increase the average fitness in the new equlibrium is larger than the probability to decrease it.

Frequency-dependent selection can arise from different mechanisms. In a haploid population, frequency-dependent selection is caused by the interactions of different types. In this case, the fitness of a particular type depends on the frequency of its own and other types in the population. However, in a diploid population, frequency dependent selection on alleles can arise also from the interactions of two alleles at one locus [8,28,29]. Our model can be easily extended to a diploid population in such a case, which leads to different results for the average change in fitness, see Appendix.

Games with n types

So far, we have discussed the change of the average fitness of a population consisting of at most two types. However, when two types coexist in a stable polymorphism, an additional type can enter the population and persist. To describe the interaction of individuals in a population with more than two types, we extend the 2 × 2 payoff matrix to a n × n payoff matrix A, where n is the number of types in the population. The entry in the i-th row and the j-th column, Aij represents the fitness of an i-type individual interacting with a j-type individual. The fitness of type i on average can be written as <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M28">View MathML</a>, where j = 1,2,3…,n, and xj is the frequency of type j, such that <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M29">View MathML</a>.

In our model, n is not a fixed number. When a type goes extinct, the corresponding row and column are deleted in the payoff matrix. Thus, the value of n decreases by one. When a mutation occurs, one row and one column are added to describe the interactions of the mutant type and resident types, which increases the size of the payoff matrix by one. The new entries introduced by a mutation are generated based on the assumption that the interactions between the mutant type m and any resident type i are similar to those between the parent type p and the resident type i. In our case, we assume amj is a random variable which is drawn from a probability density function f(x) and is larger than apj with probability θ.

Since the complexity of the population dynamics increases considerably with the number of types, it would be difficult to obtain the changes of the average fitness in a polymorphic population of n > 2 types analytically. Therefore, we use the replicator equations to simulate the dynamics of the system with several types. We start the simulation from a homogenous population. However, since we are interested in the average fitness changes and other stationary quantities averaged over a long time period, the initial number of types has no effects on the results. The time intervals are sufficiently small that at most one mutant type can appear during one time interval. The probability that a resident type i produces a mutant type is <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M30">View MathML</a>, where i = 1,2,3,…,n. Thus the probability that a mutant arises from a resident type i increases with the fitness of this type. However, for the whole population, the probability that a mutant type appears is just the mutation rate, <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M31">View MathML</a>.

We can chose arbitrary mutation rates in our simulations. However, when the mutation rate is very high, a population might experience a new mutation when it is still in a non-equilibrium state triggered by the previous mutation. In this case, the fate of a mutant is not only driven by selection, but also by the interplay of mutations. Since we are interested in the fitness consequences of frequency dependent selection, we choose the mutation rate small enough such that a population disturbed by a mutation reaches the new equilibrium before the next mutation arises.

We first look at the transition probability between different levels of diversity under mutation and selection. Once a mutation occurs it can coexist with all resident types, replace one resident type, outcompete some resident types, or go extinct. The transition matrix T describes this dynamics. Suppose the number of types in the current population is n. The element Tni denotes the transition probability from n to i coexisting types, where i = 1,2,3,…,n + 1, see Figure 2. We obtain the values in the transition matrix from numerical simulations. Every transition event triggered by a mutation is recorded and the probability to go from a certain number of types to another number of types is averaged over many realizations. These transition probabilities show some interesting properties. The probability to keep the current diversity (the element in the main diagonal in a row) is always higher than the probabilities to decrease or increase the diversity (all the other elements in the same row), see Figure 2 and Ref. [30]. Interestingly, for a population consisting of less than 4 types, the probability to increase the diversity Tii + 1 is higher than the probability to decrease the diversity <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M32">View MathML</a> in the parameter regime of Figure 2. Once the population reaches the threshold of 4 types, this pattern reverses. Thus in the long run the population tends towards an intermediate level of diversity. Furthermore, we observe the ranking, T12 > T23 > T34 > T45. This suggests that the probability to reach higher levels of diversity decreases with increasing diversity even for larger number of initial types. The transition probability from one type to a two-type coexistence can be calculated analytically based on the comparison of payoff entries, see above. Thus, T12 = p(a < c)p(b > d) = θ/2, which is confirmed by our simulation results of T12 under different θ for the n-type model.

thumbnailFigure 2. Transition probabilities between different levels of diversity. The entry in row i and column j is the transition probability from a stable coexistence of i types to a stable coexistence of j types, numbers are also color coded. The mutation rate is so low that the transitions between different states are caused by the appearance of a single mutation. The higher the number of coexisting types is, the more difficult the state is to be reached. Here we show the transition for up to six co-existing types (θ = 0.5, averages obtained over 500 independent realizations and 20000 mutations per realization).

For a population with n types, the changes of the average fitness are more complicated, as the interactions between different types are much more diverse than in a two-type population. Even a classification of different types of interactions in such a population is difficult and of limited value to understand the change in average fitness. Instead, we evaluate the changes of the average fitness between these states numerically.

A mutation can increase, maintain, or decrease the diversity level of the population. We present the changes of the average fitness in these three scenarios, see Figure 3, for those transitions which happen most frequently (see Figure 2). For small θ, mutants are more likely to obtain lower fitness than their parents type does, in the interactions with the same resident type. This can cause the decrease of the average fitness in all three situations. If θ is sufficiently small, the average fitness will decrease all the time. When θ becomes larger, the average fitness can increase. The larger θ is, the larger the increase is. Thus, our results under the replicator dynamics provide not only the change of the average fitness under a constant θ, but also the direction and magnitude of the average fitness changes. In real systems, one may expect that θdecreases during the adaption of the population. However, e.g. environmental changes could also increase it.

thumbnailFigure 3. Changes in the average fitness when a population evolves between different levels of diversity under various probabilities that a mutant payoff values is larger than the parent’s θ. The symbols are simulation results based on replicator dynamics. The number of different types can either stay the same, increase by one or decrease by any number, because at most a single mutation enters the population. Note that the average fitness of the population in the new equilibrium decreases for small θ in all three cases after a transition. Thus even if a mutant takes over a population, the average fitness can decrease. With increasing θ, the average fitness will increase over time, but the fitness gain reduces with increasing diversity. The difference among results under Gaussian distribution and uniform distribution with the same variance, shows that the absolute changes of the average fitness also depends the concrete shapes of the probability distribution (every symbol is averaged over 500 independent realizations and 20000 mutations per realization. The probability distribution f(x) is Gaussian (left) or uniform (right) with variance 1).

Games with equal gains from switching

So far, we have assumed that the payoff of the mutant interacting with another resident type is derived from the payoff of its parent interacting with the same resident type. In a population with only two types, this leads to the case where the three random payoff entries, a, b and c, are all related to d. As a null model, we have assumed that a, b and c are uncorrelated. While this is the simplest possibility, it may not be the case for concrete biological systems. Therefore, we analyze a different case here which focuses on particular cases of frequency dependence, but includes such correlations.

We focus on an evolutionary game with the payoff matrix

where ε and δ are independent random variables with probability distributions fε(x) and fδ(x) respectively. ε represents the effect of a mutation on the mutant type, and δ represents the effect of a mutation on those who interact with the mutant type. This game has the property of “equal gains from switching”, where the sum of the payoff values in the main diagonal is equal to the sum of the payoff values in the other diagonal [31]. It can arise from the assumption that the two types are close to each other in a continuous phenotype space [32]. The case of δ = 0 corresponds to constant selection. Note that there are no coexistence games when we assume such payoff matrices. If ε > 0, the mutant will take over the population (d + ε + δ > d + δ and d + ε > d), and the new average fitness becomes <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M34">View MathML</a>. Compared with the former average fitness d, the average fitness increases if ε + δ > 0, and decreases if ε + δ < 0. If ε < 0, the mutant will be outcompeted by the wild type (d + ε + δ < d + δ and d + ε < d), and the average fitness of the population remains the same. The probability to increase the average fitness becomes <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M35">View MathML</a>, where θε is the probability that ε is larger than 0, and p(ε + δ > 0 ∣ε > 0) is the conditional probability that the sum of ε and δ is larger than 0 given ε is larger than 0. This conditional probability can be written as

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M36">View MathML</a>

(15)

The values of θε and p(ε + δ > 0 ∣ ε > 0), which determine the probability that the average fitness increases, depend on the concrete choice of fε(x) and fδ(x). The integrals can only be carried out in special cases.

It is worth to mention there is a difference between games with equal gains from switching and games with independent random payoff entires on the population dynamics. In an infinite population, where genetic drift has no effect on the population dynamics, the resulting dynamics under positive frequency dependent selection and under constant selection are similar, as there are no stable coexistences. Successful mutants will invade and take over the population sequentially. The diversity will only increase if the mutation rate is high enough. On the contrary, when different kinds of interactions, especially negative frequency dependent selection, are allowed (for example, the case with independent random payoff entires), diversity can increase even for lower mutation rates (see above).

Discussion and conclusion

Mutants with high individual fitness do not necessarily increase the average fitness of the population under frequency dependent selection. Similarly, the mutants which maximize the average fitness of a population are not necessarily those leading to a stable equilibrium in this scenario. An example for a two-type population is that a mutant interacts with the wild type in a game like the Prisoners’ Dilemma [7,26]. This is a special case of a dominance game, where the defector (the mutant) outcompetes the cooperator (the wild type) and causes a reduction in the average fitness. For example, in the RNA phage ϕ6, the competitive interactions among the high multiplicities-of-infection phage (the defector) and the low multiplicities-of-infection phage (the cooperator) in the same host cell are studied, which conforms to the Prisoners’ Dilemma [33]. In this experiment, when the defector invades with a low frequency, it has higher fitness than the residents (c > a), but the average fitness decreases when the defector becomes fixed (d > a).

Since natural selection works on an individual level rather than a population level, it does not always lead to an increase of the average fitness. Our random mutant games model accommodates mutations under frequency-dependent selection, which can result in an increase or decrease in the average fitness, not only for the simplest case of two types but also for an arbitrary number of mutant types. An interesting aspect of our model is that even though it allows for an infinite number of mutant types, it does not result in a continuous growth of diversity in a population, but leads to an intermediate level of diversity [30]. We assume that the payoffs are constant in time and identical for individuals of the same type. If individuals vary in their payoffs despite being of the same type, the results are altered by this additional source of randomness [34,35]. In a population with two types, we calculate a particular value θ, where the probability that the average fitness increases is equal to the probability it decreases. The exact value of θ depends on the concrete implementation of the payoff matrix. An interesting result of our model is that the probability to decrease or increase fitness depends only on a particularly simple property of the fitness distribution. While this may not be of direct relevance to a concrete biological system, it illustrates conceptually that a decreasing fitness may not be counterintuitive even under the simplest possible assumptions of frequency dependence.

We have discussed the changes in the average fitness for an infinite asexual population under mutation and selection. Additional effects occur when the population size becomes finite and genetic drift is not negligible [30]. However, our main observation is that the average fitness at equilibrium can only increase or remain constant by random mutations under constant selection, but also decrease under frequency-dependent selection. This can shed new light on problems in evolutionary biology and leads to the exciting question on the dynamics of the average population fitness in real biological populations. In an asexual finite population, random genetic drift leads to the accumulation of deleterious mutations and an continuous decrease in the average fitness, which is well known as Muller’s ratchet [36]. Without any forms of recombination and epistasis, beneficial mutations are the only source to compensate the average fitness decline. Since the probability of increasing the average fitness by random mutations is lower under frequency-dependent selection (see Figure 1), we must conclude that asexual populations face an even bigger challenge to maintain their average fitnesses under frequency dependent selection than under constant selection in a finite population. This is particularly striking when θ is small, a case that is typically thought of as the biologically most relevant case.

In population genetics, the change of the average fitness has also been studied in diploid systems [37,38]. However, our approach starting from a different point of view, not only allows the interplay of mutation and selection, but also a wider interpretation of the fitness of heterzygotes. Suppose A and B are two alleles at the same locus. In population genetics, the fitness of genotype AB and BA is usually considered to be identical, which is a special case in our model called symmetric diploids. However, this does not hold in asymmetric diploids where the maternal allele and paternal allele are not equally expressed. Our model and our analysis allow both cases. In the framework of a well-mixed symmetric diploid population (corresponding to random mating), our result that the average fitness never decreases is consistent with the former statement in population genetics (see Appendix).

Frequency dependent interactions can lead to a decrease of the average fitness of a population during the process of evolution despite natural selection. This is because natural selection works on individual fitness instead of the average fitness of a population.

Methods

We explore the population dynamics driven by random mutations under frequency dependent selection based on the replicator dynamics. For populations with only two types, we obtain analytical results by analyzing the change of average fitness between two equilibria of the replicator equations. For a population with more than two types, we simulate the evolutionary dynamics numerically. The current group of replicator equations and the current population composition determine the equilibrium that the population moves to. Once the population reaches this equilibrium, a new mutation occurs. In our model, every mutation brings a new game and consequently an additional replicator equation. Our approach corresponds to low mutation rates, where the time a population needs to reach an equilibrium is shorter than the waiting time for the next mutation.

Appendix

Diploid populations with two alleles

The impact of Mendelian inheritance on the population dynamics has been discussed in the framework of evolutionary game theory before [25,39-41]. In a diploid population, the combinations of two alleles at a given locus on a pair of homologous chromosomes, can be interpreted by a special two player game. Suppose there are allele A and allele B. The fitness of different genotypes, WAA,WAB and WBB can be described by a 2 × 2 matrix

This is mathematically identical to the game with two types discussed above. Here, WAA corresponds to a, WAB to c = b, and WBB to d. For a population initially only with homozygotes BB, the probability of increasing the average fitness <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M38">View MathML</a> caused by a random new allele A, can be calculated by setting c = b in Eq. (10). This becomes

<a onClick="popup('http://www.biomedcentral.com/1471-2148/12/160/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/160/mathml/M39">View MathML</a>

(16)

The probability that the average fitness decreases in such a population is 0, because the diploid AB and the diploid BA is indistinguishable, c = b. In asymmetric diploids, where the maternal alleles and paternal alleles are not equally expressed, the average fitness changes are exactly the same as shown in a general case of haploid populations.

Competing interests

The authors declare no competing interests.

Author’s contributions

W.H. and A.T. designed the model. W.H., B.W. and A.T. evaluated the model and wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank Bernhard Haubold, Christoph Hauert, Christian Hilbe, and Manfred Milinski for helpful discussions and suggestions. We also thank the reviewers for valuable and constructive comments. We gratefully acknowledge funding by the Chinese Scholarship Council (W.H.), the American Association of University Women (W.H.), and the Max Planck Society (W.H., B.W. and A.T.).

References

  1. Remold SK, Lenski RE: Contribution of indiviual random mutations to genotype-by-environment interactions in Eschrichia coli.

    Proc Natl Acad Sci U S A 2001, 98:11388-11393. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Lalić J, Cuevas J, Elena SF: Effect of host species on the distribution of mutational fitness effects for an RNA virus.

    PLoS Genet 2011, 11:e1002378. OpenURL

  3. van Valen L: A new evolutionary law.

    Evolutionary Theory 1973, 1:1-30. OpenURL

  4. Brockhurst MA, Morgan AD, Rainey PB, Buckling A: Populating mixing accelerates coevolution.

    Ecol Lett 2003, 6:975-979. Publisher Full Text OpenURL

  5. Thompson J, Cunningham B: Geographic structure and dynamics of coevolutionary selection.

    Nature 2002, 417:735-738. PubMed Abstract | Publisher Full Text OpenURL

  6. Maynard SmithJ, Price GR: The logic of animal conflict.

    Nature 1973, 246:15-18. Publisher Full Text OpenURL

  7. Maynard SmithJ: Evolution and the Theory of Games. Cambridge: Cambridge University Press; 1982. OpenURL

  8. Hofbauer J, Sigmund K: Evolutionary Games and Population Dynamics. Cambridge: Cambridge University Press; 1998. OpenURL

  9. Cressman R: Evolutionary Dynamics and Extensive Form Games. Cambridge: MIT Press; 2003. OpenURL

  10. Nowak MA, Sigmund K: Evolutionary Dynamics of Biological Games.

    Science 2004, 303:793-799. PubMed Abstract | Publisher Full Text OpenURL

  11. Nowak MA: Evolutionary Dynamics. Cambridge: Harvard University Press; 2006. OpenURL

  12. Sanjuán R, Moya A, Elena SF: The distribution of fitness effects caused by single nucleotide substitutions in an rna virus.

    Proc Natl Acad Sci U S A 2004, 101:8396-8401. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Orr HA: The distribution of fitness effects among beneficial mutations.

    Genetics 2003, 163:1519-1526. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Gillespie JH: Some properties of finite populations experiencing strong selection and weak mutation.

    Am Naturalist 1983, 121:691-708. Publisher Full Text OpenURL

  15. Zeyl C, DeVisser JA: Estimates of the rate and distribution of fitness effects of spontaneous mutation in Saccharomyces cerevisiae.

    Genetics 2001, 157:53-61. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Cowperthwaite MC, Bull JJ, Meyers LA: Distributions of beneficial fitness effects in RNA.

    Genetics 2005, 170:1449-1457. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Orr HA: The distribution of fitness effects among beneficial mutations in fisher’s geometric model of adaptation.

    J Theor Biol 2006, 238:279-285. PubMed Abstract | Publisher Full Text OpenURL

  18. Eyre-Walker A, Keightley PD: The distribution of fitness effects of new mutations.

    Nature 2007, 8:610-618. OpenURL

  19. Lynch M, Blanchard J, Houle D, Kibota T, Schultz S, Vassilieva L, Willis J: Spontaneous deleterious mutation.

    Evolution 1999, 53:645-663. Publisher Full Text OpenURL

  20. Eyre-Walker A, Woolfit M, Phelps T: The distribution of fitness effects of new deleterious amino acid mutations in humans.

    Genetics 2006, 173:891-900. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Huang W, Traulsen A: Fixation probabilities of random mutants under frequency dependent selection.

    J Theor Biol 2010, 263:262-268. PubMed Abstract | Publisher Full Text OpenURL

  22. Levin SA, Grenfell B, Hastings A, Perelson AS: Mathematical and computational challenges in population biology and ecosystems science.

    Science 1997, 275:334-343. PubMed Abstract | Publisher Full Text OpenURL

  23. Taylor PD, Jonker B: Evolutionary stable strategies and game dynamics.

    Math Biosci 1978, 40:145-156. Publisher Full Text OpenURL

  24. Zeeman EC: Population dynamics from game theory.

    Lecture Notes in Mathematics 1980, 819:471-497. Publisher Full Text OpenURL

  25. Hofbauer J, Schuster P, Sigmund K: Game Dynamics in Mendelian Populations.

    Biol Cybernetics 1982, 43:51-57. Publisher Full Text OpenURL

  26. Rapoport A, Chammah AM: Prisoner’s Dilemma. Ann Arbor: University of Michigan Press; 1965. PubMed Abstract OpenURL

  27. Nowak MA: Five rules for the Evolution of Cooperation.

    Science 2006, 314:1560-1563. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Cressman R: The stability concept of evolutionary game theory: a dynamics approach.

    Lecture Notes in Biomathematics 1992, 94:1-128. Publisher Full Text OpenURL

  29. Traulsen A, Reed FA: From genes to games: Cooperation and cyclic dominance in meiotic drive.

    J Theor Biol 2012, 299:120-125. PubMed Abstract | Publisher Full Text OpenURL

  30. Huang W, Haubold B, Hauert C, Traulsen A: Emergence of stable polymorphism driven by evolutionary games between mutants.

    Nat Comm 2012, 3:919. OpenURL

  31. Nowak MA, Sigmund K: The evolution of stochastic strategies in the prisoner’s dilemma.

    Acta Applicandae Mathematicae 1990, 20:247-265. Publisher Full Text OpenURL

  32. Wild G, Traulsen A: The different limits of weak selection and the evolutionary dynamics of finite populations.

    J Theor Biol 2007, 247:382-390. PubMed Abstract | Publisher Full Text OpenURL

  33. Turner PE, Chao L: Prisoner’s Dilemma in an RNA virus.

    Nature 1999, 398:441-443. PubMed Abstract | Publisher Full Text OpenURL

  34. Fudenberg D, Harris C: Evolutionary Dynamics with aggregate shocks.

    J Econ Theory 1992, 57:420-441. Publisher Full Text OpenURL

  35. Johnson DDP, Stopka P, Bell J: Individual variation evades the Prisoner’s Dilemma.

    BMC Evolutionary Biol 2002, 2:15. BioMed Central Full Text OpenURL

  36. Muller HJ: The relation of recombination to mutational advance.

    Mutat Res 1964, 1:2-9. Publisher Full Text OpenURL

  37. Fisher RA: The Genetical Theory of Natural Selection. Oxford: Clarendon Press; 1930. OpenURL

  38. Fisher RA: Average excess and average effect of a gene substitution.

    Ann Hum Genet 1941, 11:53-63. Publisher Full Text OpenURL

  39. Hines WGS: An evolutionarily stable strategy model for randomly mating diploid populations.

    J Theor Biol 1980, 87:379-384. PubMed Abstract | Publisher Full Text OpenURL

  40. Eshel I: Evolutionarily stable strategies and viability selection in Mendelian populations.

    Theor Population Biol 1982, 22:204-217. Publisher Full Text OpenURL

  41. Brown RLW: Evolutionary game dynamics in diploid populations.

    Theor Population Biol 1983, 24:313-322. Publisher Full Text OpenURL