The impact of random frequency-dependent mutations on the average population fitness
© Huang et al.; licensee BioMed Central Ltd. 2012
Received: 30 March 2012
Accepted: 22 August 2012
Published: 30 August 2012
Skip to main content
© Huang et al.; licensee BioMed Central Ltd. 2012
Received: 30 March 2012
Accepted: 22 August 2012
Published: 30 August 2012
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.
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.
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.
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 . 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 , which considers evolutionary processes under frequency dependent selection . 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.
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 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.
The change of the wild type frequency follows immediately as . 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.
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.
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.
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.
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, A ij 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 , where j = 1,2,3…,n, and x j is the frequency of type j, such that .
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 a mj is a random variable which is drawn from a probability density function f(x) and is larger than a pj 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 , 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, .
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.
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.
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.
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).
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 . 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 . 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 . 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 . 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.
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.
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.
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.).
This article is published under license to 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.