Email updates

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

Open Access Highly Accessed Research article

Selection for altruism through random drift in variable size populations

Bahram Houchmandzadeh* and Marcel Vallade

Author affiliations

Univ. Grenoble 1/CNRS, LIPhy UMR 5588, Grenoble, F-38401, France

For all author emails, please log on.

Citation and License

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

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


Received:9 November 2011
Accepted:14 March 2012
Published:10 May 2012

© 2012 Houchmandzadeh and Vallade; 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

Altruistic behavior is defined as helping others at a cost to oneself and a lowered fitness. The lower fitness implies that altruists should be selected against, which is in contradiction with their widespread presence is nature. Present models of selection for altruism (kin or multilevel) show that altruistic behaviors can have ‘hidden’ advantages if the ‘common good’ produced by altruists is restricted to some related or unrelated groups. These models are mostly deterministic, or assume a frequency dependent fitness.

Results

Evolutionary dynamics is a competition between deterministic selection pressure and stochastic events due to random sampling from one generation to the next. We show here that an altruistic allele extending the carrying capacity of the habitat can win by increasing the random drift of “selfish” alleles. In other terms, the fixation probability of altruistic genes can be higher than those of a selfish ones, even though altruists have a smaller fitness. Moreover when populations are geographically structured, the altruists advantage can be highly amplified and the fixation probability of selfish genes can tend toward zero. The above results are obtained both by numerical and analytical calculations. Analytical results are obtained in the limit of large populations.

Conclusions

The theory we present does not involve kin or multilevel selection, but is based on the existence of random drift in variable size populations. The model is a generalization of the original Fisher-Wright and Moran models where the carrying capacity depends on the number of altruists.

Keywords:
Frequency independent fitness; Genetic drif; Fixation probabilities; Non-structured populations

Background

Light production in Vibrio fischeri[1,2], siderophore production in Pseudomonas aeruginosa[3], invertase enzyme production in Saccharomyces cerevisiae[4], stalk formation by Dictyostelium discoideum, [2,5] are but a few examples of individuals in a community who help others at their own cost by devoting part of their resources to this task. This behavior has been termed “altruistic”. From the evolutionary point of view, altruists have a lower fitness than other individuals in the community who don’t help, but are recipient of the benefits produced by altruists. Through this paper, we will call these latter individuals ‘selfish’.

From the inception of evolution theory, the problem of the existence of altruists has been puzzling: how can a mutant with lower fitness prevail? And how does a community of altruists resist the spread of selfish allele (see [6] for a historical perspective)? In the last 40 years many models have emerged to explain the apparent contradiction between the smaller fitness of altruists and their widespread presence in various communities (for a review, see [7,8]). It is shown in these models that the actual fitness of an altruistic gene can be increased by other factors such as ‘common good’ restricted to kin (inclusive fitness [9,10]), or advantages conferred at another level of selection (group or multilevel selection [11,12]). These models which can be formulated through the Price equation have seen various generalizations and they are sometimes widely debated (see [13] and the numerous replies it has elicited).

The above models are either deterministic, i.e. populations change their size exactly according to their relative fitness, or involve frequency dependent fitness [14,15]. We show here that another possibility exists: an altruistic individual can produce a common good benefiting everybody in the community regardless of its nature (altruistic or selfish) and therefore increasing the carrying capacity of the habitat. Even though selfish individuals have always a higher fitness, genetic drift effects can favor the altruists.

It was established by the founding fathers of Population Genetics that a mutation that confers a relative fitness 1 + s does not automatically spread and take over the whole community, but has only a higher probability, called the fixation probability, to do so [16-18]. For a community of fixed size N of haploid individuals, the fixation probability π of a mutant appearing at one copy, for small selection pressure Ns < <1, is

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

(1)

The fixation probability is composed of two terms: even in the absence of selection, the population will become homogenic; in this neutral case, all individuals at generation zero have an equal probability 1/N of becoming fixed. When a beneficial mutation is present, the fixation probability of its carrier is increased by the relative excess fitness.

For populations of fixed size, as can be seen from expression (1) or the more precise expression (10) obtained by Kimura [19] and Moran [20], the fixation probability is a monotonically increasing function of the sole relative fitness. In the competition between alleles, arguments based on fitness parameter alone or the fixation probability lead to the same conclusions . However, if population size is not fixed, the fixation probability π, which takes into account both randomness due to finite size and selection, can lead to other conclusions than the fitness parameter alone.

Consider an altruistic gene that by some means (production of a ‘common good’, limited grazing of natural resources, …) allows the carrying capacity to increase: if the community were composed only of altruists its population size would be Nf; if it were composed only of selfish individuals the population size would be Ni (Ni < Nf) (Figure 1a). The production of common good decreases the relative fitness of altruists by s.

thumbnailFigure 1. Variable carrying capacity.

(a): A community where the carrying capacity is an increasing function of altruist number, varying from Ni when the population is composed only of selfish individuals to Nf (Nf > Ni) when only altruists are present. (b) Two examples of random walks describing the stochastic behavior a such a system (transition probabilities 4–7), where m,n are the number of selfish and altruistic individuals. Red line: loss of S’s; Blue line: loss of A. A Moran process in this scheme corresponds to a random walk constrained to remain on an anti-diagonal line.

Consider now the fixation probability πA of one altruist mutant appearing in a community of Ni selfish individuals. A crude use of expression (1) shows that <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M2">View MathML</a>. On the other hand, the fixation probability πs of one selfish individual appearing in a community of Nf  altruists is <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M3">View MathML</a>. We see that if

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

(2)

i.e. the cost to the altruist is smaller than the benefits in term of relative population increase, then an altruist has a larger fixation probability than a selfish one, even though its relative fitness is smaller. The relative advantage of a selfish mutant is compensated by the increased ‘random noise’ to which it is exposed. Note that in a deterministic model of the above process, the A always lose, since S individuals always increase their proportion.

The above argument will be refined in the following. In the next section, we formulate precisely the stochastic process of altruism outlined above by generalizing the Moran model for non-structured, well mixed populations and we show that altruists can indeed be favored in their competition with selfish individuals. We outline the amplification of this advantage in geographically structured, viscous populations in the third section. The final section is dedicated to concluding remarks.

Results and discussion

Stochastic model for altruism

The fundamental aspects of population genetics were clarified in the framework of the classical Fisher-Wright (FW) stochastic model of non-overlapping generations or its continuous time alternative introduced by Moran [20]. Moran and FW are equivalent in the limit of large populations, where both are well approximated by the same diffusion equation [21]. These are the simplest models that capture the key elements of population genetics (genetic drift, fixation probability, fixation time,…) with the fewest possible ingredients.

In the Moran model, a population of size N is composed of two types of individual, say A and S. Empty spots are created randomly with fixed rate α, increasing the carrying capacity by unity. Once an empty spot has been created, it will be colonized by the progeny of either an A or an S individual according to their proportion in the community. In order to keep the population constant, Moran added the constraint that the colonization of a new spot be followed immediately by the death of an individual in the community, restoring the population size to N. Moran is therefore a simultaneous model of duplication and annihilation; the transition probability densities for the A to increase or decrease their number n by one individual are

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

(3)

where m is the number of S individuals and c is the ‘cost’: 1/ c is the relative fitness of the A and c > 1 indicates a selective disadvantage. W+ stands for the probability density that the new spot is colonized by an A and death occurs among the S. In principle, a similar set of equations must be written for the S individuals; however, as the population size is fixed, <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M6">View MathML</a>, the quantity m in eq.(3) can be replaced by N-n and the whole stochastic process treated as a one dimensional random walk for the A.

We generalize this model by including two ingredients. First, the fixed size constraint can be relaxed and we let N vary between two bounds Ni and Nf: empty spots are created-colonized and individuals die, without these two events necessarily succeeding each other. More importantly, in order to include the effects of altruists, we suppose that the rate of creation of empty spots is proportional to the number of altruists and is equal to αn; in contrast, the death rate is proportional to the number of S individuals and is equal to α m. This is the simplest hypothesis that implies that the increase in the carrying capacity of the habitat is proportional to the number of altruists (see also Methods, mean field approximation).

The stochastic model that captures all these features is a two dimensional random walk with the following transition probability densities (Figure 1b):

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

(4)

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

(5)

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

(6)

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

(7)

Consider for example the first two lines of the above equations, which are about birth events: the factor <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M11">View MathML</a> is the relaxation of Moran constraints and insures that population size remains below Nf; the factor αn accounts for the fact that empty spot creations are proportional to the number of A; finally, once a birth event has occurred, the probability for it to be an A or an S is proportional to the number of the corresponding sub-populations present at this time. The last two lines, which govern population decrease, are similar: the factor <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M12">View MathML</a> ensures that population size remains above Ni; the factor αm is the death rate (population decrease) for everybody due to the presence of selfish individuals. The cost of altruism is included in these equations: the proportion of A is <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M13">View MathML</a>, but once a death event has occurred, the probability for it to be an A is:

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

(8)

if c > 1. The results below don’t change significantly if the cost of altruism is included in other rates. For example, a higher probability for an S to reproduce, or any combination that favors S over A. Note that if the increase/decrease rates were independent of m and n, we recover the Moran model by setting <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M15">View MathML</a>, in which case each birth/death is succeeded by a death/birth event (see Methods, relation to Moran model).

The above rates ensure that if A are lost ( n = 0), the population size tends toward Ni and if S are lost ( m = 0), it tends toward Nf. Note that in the mean field approximation of the above process where fluctuations are neglected and the deterministic limit is taken, the A are always eliminated if c > 1 (see Methods, mean field approximation).

In finite size populations however, fluctuations play an important role. The focus of this paper is the computation of the fixation probability of the above process and the probability that altruists or selfish mutants take over the community. The fixation probability π( k) of a general stochastic process beginning with the initial state k and fixing either to k0 or <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M16">View MathML</a> is the solution of Kolmogorov backward equation which is a linear set of equations [22]

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

(9)

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

(10)

where the sum is over all the states q attainable from the state k with transition probabilities <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M19">View MathML</a>. For one dimensional, one step processes such as Moran, k =n and the solution of the linear system is easily obtained [22]:

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

(11)

where <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M21">View MathML</a> is the proportion of the A. The approximation corresponds to the Kimura solution obtained through a backward diffusion equation [19] and <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M22">View MathML</a>. Expression (1) is the first order expansion of the above expression in s.

For the two dimensional process (4–7) where <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M23">View MathML</a> is the initial number of the S and A, no closed form solution can be obtained. We can however solve equation (8) numerically by standard linear solvers or else resort to a Gillespie algorithm [23] to solve the stochastic equations (4–7) directly. Both these methods are used in this paper and the analytical approximations obtained below are compared to them.

For large populations, we use the usual diffusion equation approximation of eq.(8) [19,22]. For weak selection pressure, the diffusion approximation error for the simple Moran process is <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M24">View MathML</a>[24]; for more general cases, the validity of the approximation has been discussed by Zhou and Qian [25]. Setting <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M25">View MathML</a> and denoting π( x, y) the fixation probability for the initial composition ( x, y), the diffusion equation reads:

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

(12)

where

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

(13)

and <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M28">View MathML</a><a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M29">View MathML</a>. This is a complicated elliptic partial differential equation. In the absence of selection (c = 1) however, the trivial neutral solution is <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M30">View MathML</a> which as expected, is just the proportion of altruists. Building upon this solution, and denoting <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M31">View MathML</a> for the proportion of altruists and <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M32">View MathML</a>, we can check that to the first order of perturbation <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M33">View MathML</a>, the solution reads

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

(14)

where

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

(15)

and γ is a numerical coefficient: <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M36">View MathML</a>. The first order perturbation solution (12), which was derived for small selection pressures <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M37">View MathML</a>, proves in fact to be an excellent approximation for selection pressure as high as <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M38">View MathML</a>, (Figure 2).

thumbnailFigure 2. Fixation probabilities.

Comparison of analytical solution (12) (solid lines) to numerical solution of eq.(8) for increasing selection pressure indicated by the arrows: <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M39">View MathML</a>. Nf = 100, Ni = 90.

The general solution (12) allows for the computation of the fixation probability of one individual introduced into a community of the other type. To the first order of perturbation in <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M40">View MathML</a>, the fixation probabilities πA of one A introduced in a community of S reads:

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

(16)

and the fixation probabilities πs of one S introduced in a community of A is

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

(17)

Figure (3a) shows the evolution of these probabilities as a function of selection pressure for various Ni and Nf. Equations (1314) show that the condition for the altruist to be favored, <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M43">View MathML</a>, is simply

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

(18)

where <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M45">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M46">View MathML</a> and we have kept only the leading terms. <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M47">View MathML</a> is the equilibrium relative excess cost of altruism at which A and S individuals become equivalent. Figure 3b shows the excellent agreement between the above results and exact numerical results. Altruists have a selective advantage if the selection pressure against them, i.e. the combined effect of fitness and population size, is smaller than the relative increase in population size. Unlike a Hamilton rule, criterion (15) is a finite size effect and is of purely stochastic nature : because of the demographic effect, selfish mutants are submitted to a higher stochastic noise than altruist; this can be sufficient to prevent them from prevailing. Note that the above computations were performed for the limiting case of weak selection (Nfs < < 1), which is considered by most, but not all, scientists, to be the relevant limit of evolutionary dynamics [26,27]. Direct numerical resolution of eq. (8) shows however that an equilibrium excess fitness exists even at high selection pressure, given a high enough relative increase in population size.

thumbnailFigure 3. Criterion for Altruists selection. ( a) fixation probabilities πs (red squares) and πA (blue circles) as a function of selection pressure <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M48">View MathML</a>, for Nf = 100 and <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M49">View MathML</a>. Values are obtained by numerically solving eq.(8). Solid lines are theoretical values (eqs. 13,14). Increasing Ni are indicated by the arrow. (b) equilibrium selection pressure <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M50">View MathML</a> for which πA = πs for multiple combinations of <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M51">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M52">View MathML</a>, as a function of relative population increase, obtained numerically. The solid line is the theoretical value (eq.15).

Geographically structured populations

The altruists’ advantage can be enhanced for large structured populations [28-31]. Geographically structured populations can be modeled as divided into colonies that exchange migrants [32]. The Moran model on graph is a non trivial problem [33]; we restrict our treatment here to the simplest case where the migration time scale is small compared to fixation time of one mutant (viscous populations): a migrant is either lost or fixed before a new migration event happens. The argument we develop below is similar to the two level model of Traulsen and Nowak [34]. Consider a one dimensional community subdivided into M colonies (Figure 4), exchanging migrants with neighboring patches at rate m. As the migration event is rare, these colonies are fixed either into an A or S state. The probability density per unit time pSA for an S colony on the border to become an A colony is to receive one migrant from the neighboring A colony multiplied by the probability that this mutant gets fixed:

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

(19)

thumbnailFigure 4. Geographically structured populations.

Geographically structured population where patches can exchange migrants. For low migration rates, the border between A and S domains can be modeled as a biased random walk.

Similarly the probability density for an A colony on the border to become S is

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

(20)

Therefore, the movement of the border itself can be considered a biased random walk. The probability ΠA for an altruist mutant to take over the whole community is thus the probability for a mutant to take over one colony and then for this colony to take over the whole community:

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

(21)

where <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M56">View MathML</a>. If the criterion (15) is satisfied, then obviously r < 1 and for large number of communities M > >  1,

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

(22)

On the other hand, the probability Πs for a selfish mutant to be fixed is

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

(23)

and <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M59">View MathML</a> for M > >  1: once altruists dominate, the chances for a selfish mutant to invade the community is close to zero! Increased random noise due to production of common good and a small migration rate are an efficient way of keeping selfishness in check.

The above computation concerns the low migration limit. In the high migration limit, the community is non-structured and its effective size is <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M60">View MathML</a>. Criterion (15) shows that in this regime, altruists cannot emerge; this is indeed equivalent to the deterministic case where emergence of altruists calls for other mechanisms. Between these two regimes of high and low migration rate, there is a rich interval where migration rate is a key ingredient in the competition between altruists and cheaters.

Conclusions

The main concepts of Population Genetics were clarified in the framework of the original model of Fisher-Wright and Moran (FWM). These models introduced the key ingredient of population size and its role in the randomness of selection. It became clear in the 1920-30’s that a beneficial mutation does not spread automatically to the whole population, but has to overcome the “random noise” of population sampling over generations. The idea that random noise plays also a role for the selection of altruism has been introduced in two kind of models, which have a marked difference with the model we present here. The first class of models, formulated mostly through evolutionary game theory formalism, concerns fixed size populations, where the transition rates are frequency dependent [14]: the fitness of an A individual can be superior to the fitness of an S individual if the number of A individuals already present is high enough. It can then be shown, upon very general conditions, that the fixation probability of altruists can become superior to that of selfish ones. These models can be seen as the generalization of Hamilton’s original idea, where “altruistic” help is restricted to genetically related individuals, even though Traulsen [35] has argued that the underlying mathematics is fundamentally different. The second class of models concerns group (or multilevel) selection. It has been shown [34] that the fixation probability of altruists can be higher than those of selfish ones, if the population is structured into groups and the splitting of one group leads to the elimination of another. It has also recently been noticed that random noise in a growing population can favor altruists during a transient period [36].

The model we present here is not frequency dependent: an A individual has always a lesser chance of reproducing than an S individual; the mean field description of this model has only one stable fixed point which corresponds to the disappearance of altruists. Moreover, The mechanism we propose is for non-structured populations, even though the altruist effect can be amplified when the population is structured into groups with small migration rate between groups. Imagine a group of M islands composed only of altruists and another group of M islands composed only of selfish individuals. Introduce one S mutant in each island of the first group and one A mutant in each islands of the second group. After some time, the number of islands in the first group is increased if the criterion (15) is satisfied.

In summary, we have shown, by a slight generalization of the Moran model, that in finite size populations, the fixation probability of altruists can be higher than that of selfish ones, even though their fitness is lower and their emergence is ‘forbidden’ by a Hamilton rule. We have also shown that in large, structured populations, and in the limit of small migration rate, the same arguments hold. Production of the ‘common good’ and increase in the carrying capacity of the habitat increase the random noise for selfish individuals and can therefore favor altruists.

The aim of the present article is not to contest the merits of kin/group selection models which have been investigated during the last forty years with a large number of case studies. We believe we are providing an alternative way of thinking about altruism which is complementary to the above models and which restores the key ingredients of population genetic to this topic.

Methods

Diffusion equation derivation

In the discrete backward Kolmogorov eq. (8) set <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M61">View MathML</a> and q all the states reachable from k, i.e. all states of the form ( m ± 1, n) and m, n ±  1. The equations read

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

(24)

For large populations <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M63">View MathML</a>, we set <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M64">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M65">View MathML</a> and develop the above expression to the second order in <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M66">View MathML</a> (Kramers-Moyal expansion). Combining all the resulting terms leads to the partial differential equation (11). It is fruitful to express this equation in terms of total relative population <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M67">View MathML</a> and proportion of altruists <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M68">View MathML</a>; the inside domain shown in Figure 1 then maps into the <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M69">View MathML</a> rectangle, where <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M70">View MathML</a>. In these coordinates, the diffusion equation reads:

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

(25)

where

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

(26)

Mean field approximation

In the deterministic approximation, fluctuations are neglected. Denoting by m and n the ensemble average of the number of S and A individuals, their deterministic evolution equation reads:

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

(27)

It is more fruitful to write directly the evolution of the proportion of A-individuals <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M74">View MathML</a>. Using the expression for transition probabilities (4–7), we have

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

(28)

where <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M76">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M77">View MathML</a>. It is then obvious that for <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M78','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M78">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M79','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M79">View MathML</a>. In the deterministic model, A-individuals always disappear.

The equation for total population reads

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

(29)

for <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M81">View MathML</a>, the stationary solution of this equation, assuming that μ is held constant is

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

(30)

which shows that the increase in carrying capacity of the habitat η - k, at small selection pressure, is mostly proportional to the number of A-individuals. A closer look at the above equations (16,18) shows that <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M83">View MathML</a> is the only stable fixed point when c >  1.

Relation to Moran model

In a simple model where population size is variable, but birth and death rates are independent of the number of altruists and selfish individuals, a constant α will replace (αn) and (α m) in equations (4–7). In the case where <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M84">View MathML</a>, the stochastic movement pictured in Figure 1b reduces to a movement on the anti-diagonal staircase: births and deaths occur only when the total population <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M85">View MathML</a> is respectively equal to Ni and Nf. The analog of the Moran process is obtained by computing the two steps transition probabilities <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M86">View MathML</a>. If <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M87">View MathML</a>, this implies first the birth of one individual of one type and then the death of an individual of the other type. Combining the rates given by eqs.(4–7) where birth and death rates are constant, we obtain

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

(31)

The same expression is obtained if <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M89','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M89">View MathML</a>.

Numerical resolution of fixation probabilities

Two different kinds of numerical resolution were used to check the validity of our analytical results on the fixation probabilities: A Gillespie stochastic algorithm and direct resolution of eq. (8).

Gillespie algorithm

The stochastic equations given by the rates (4–7) can be seen as 2 chemical reactions for the species <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M90','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M90">View MathML</a>:

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

(32)

which we solve by the classical Gillespie algorithm [23] written in C++. We are interested here only in the fixation probability and not in the fixation time; the program can therefore be accelerated by computing only the nature of the event that occurs at each turn (and not its time of occurrence). In general, to solve for the fixation probability, R = 106 stochastic trajectories are generated.

Direct resolution

Equation (8) constitutes a linear system and can therefore be solved by standard numerical packages. For the present case however, the unknowns, i.e. the fixation probabilities π( m,n) don’t constitute a vector but a second rank tensor; the tensor formed from the rates W is of rank 4. To adapt our linear system to standard linear solvers, we have to re-index the unknowns and decrease their rank by one: <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M92','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M92">View MathML</a>. We have chosen the following scheme, which corresponds to a sequential scanning of the anti-diagonal lines (Figure 5a):

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

(33)

where <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M94','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M94">View MathML</a>. The <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M95','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M95">View MathML</a> points belong to the interior of the trapezoid <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M96">View MathML</a>, n ≥  1, m ≥  1.

thumbnailFigure 5. Tensor reindexation.

(a) To each 2 d index ( m,n), a new 1 d index k is associated by scanning sequentially the anti-diagonal lines. ( b) The re-indexation transforms the tensorial equation (8) into a normal linear system <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M97','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M97">View MathML</a>, where πk are the unknowns.

The re-indexation transforms the equation (8) into a normal linear system

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

(34)

where I(k) designates the 1 d indexes of the four nearest neighbors of the point ( m,n), where <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M99','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M99">View MathML</a>. The above equations can be written in standard matrix notation

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

(35)

where <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M101','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M101">View MathML</a> are the unknowns. <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M102','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M102">View MathML</a> is a sparse matrix, which apart from the diagonal elements, has at most four non-zero elements per line: if k is the image of element ( m,n), then <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M103','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M103">View MathML</a> only if k′ is the image of one of the four nearest neighbors of ( m,n), in which case its value is given according to rates (4–7). The right hand side vector Bk is a sparse vector provided by the limit conditions <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M104','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M104">View MathML</a>: if k′, one of the 4 nearest neighbors of the element k belongs to the border m = 0, then <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M105','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M105">View MathML</a> and the corresponding <a onClick="popup('http://www.biomedcentral.com/1471-2148/12/61/mathml/M106','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/12/61/mathml/M106">View MathML</a> is transferred to the right hand side to constitute the vector Bk. Note that because we index the interior of the trapezoid, the index k itself can never belong to the border.

Once the linear system (20) has been constituted, it can be solved by any linear solver. We have used the commercial package matlab for these manipulations.

Authors’ contributions

BH designed the reasearch, performed the numerical work and wrote the article. BH and MV performed the analytical computations. Both authors read and approved the final manuscript.

Acknowledgements

We are grateful to O. Rivoire and E. Geissler for fruitful discussions.

References

  1. Visick KL, Ruby EG: Vibrio fischeri and its host: it takes two to tango.

    Curr Opin Microbiol 2006, 9(6):632-638. PubMed Abstract | Publisher Full Text OpenURL

  2. Foster KR, Shaulsky G, Strassmann JE, Queller DC, Thompson CRL: Pleiotropy as a mechanism to stabilize cooperation.

    Nature 2004, 431(7009):693-696. PubMed Abstract | Publisher Full Text OpenURL

  3. West SA, Buckling A: Cooperation, virulence and siderophore production in bacterial parasites.

    Proc Biol Sci 2003, 270(1510):37-44. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Gore J, Youk H, van Oudenaarden A: Snowdrift game dynamics and facultative cheating in yeast.

    Nature 2009, 459(7244):253-256. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Kessin RH: Dictyostelium: Evolution, Cell Biology and the Development of Multicellularity. Press Cambridge: Cambridge Univ; 2001. OpenURL

  6. Dugatkin L: The Altruism Equation. Princeton University Press; 2006. OpenURL

  7. Nowak MA: Five rules for the evolution of cooperation.

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

  8. West SA, Griffin AS, Gardner A: Evolutionary explanations for cooperation.

    Curr Biol 2007, 17(16):R661-R672. PubMed Abstract | Publisher Full Text OpenURL

  9. Hamilton WD: The genetical evolution of social behaviour.

    I. J Theor Biol 1964, 7:1-16. Publisher Full Text OpenURL

  10. Maynard Smith J: Group selection and kin selection.

    Nature 1964, 201(4924):1145-1147. Publisher Full Text OpenURL

  11. Lewontin RC: The units of selection.

    Annu Rev Ecol Syst 1970, 1:1-18. Publisher Full Text OpenURL

  12. Wilson DS: The group selection controversy: history and current status.

    Annu Rev Ecol Syst 1983, 14:159-187. Publisher Full Text OpenURL

  13. Nowak MA, Tarnita CE, Wilson EO: The evolution of eusociality.

    Nature 2010, 466(7310):1057-1062. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Nowak MA, Sasaki A, Taylor C, Fudenberg D: Emergence of cooperation and evolutionary stability in finite populations.

    Nature 2004, 428(6983):646-650. PubMed Abstract | Publisher Full Text OpenURL

  15. Traulsen A, Claussen JC, Hauert C: Coevolutionary dynamics: from finite to infinite populations.

    Phys Rev Lett 2005, 95(23):238701. PubMed Abstract | Publisher Full Text OpenURL

  16. Fisher R: The genetical theory of natural selection, a complete variorum edition. Oxford University Press; 1999. OpenURL

  17. Wright S: Evolution in Mendelian populations.

    Genetics 1931, 16(2):97-159. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Patwa Z, Wahl LM: The fixation probability of beneficial mutations.

    J R Soc Interface 2008, 5(28):1279-1289. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Kimura M: On the probability of fixation of mutant genes in a population.

    Genetics 1962, 47:713-719. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Moran P: The Statistical processes of evolutionary theory. Oxford University Press; 1962. OpenURL

  21. Houchmandzadeh B, Vallade M: Alternative to the diffusion equation in population genetics.

    Phys Rev E Stat Nonlin Soft Matter Phys 2010, 82(5 Pt 1):051913. PubMed Abstract | Publisher Full Text OpenURL

  22. Ewens WJ: Mathematical Population Genetics. Springer-Verlag; 2004. OpenURL

  23. Gillespie DT: Exact stochastic simulation of coupled chemical reactions.

    J Phys Chem 1977, 81(25):2340-2361. Publisher Full Text OpenURL

  24. Ethier SN, Norman MF: Error estimate for the diffusion approximation of the Wright–Fisher model.

    Proc Natl Acad Sci U S A 1977, 74(11):5096-5098. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Zhou D, Qian H: Fixation, transient landscape, and diffusion dilemma in stochastic evolutionary game dynamics.

    Phys Rev E Stat Nonlin Soft Matter Phys 2011, 84(3–1):031907. PubMed Abstract | Publisher Full Text OpenURL

  26. Bell G: Fluctuating selection: the perpetual renewal of adaptation in variable environments.

    Philos Trans R Soc Lond B Biol Sci 2010, 365(1537):87-97. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Wu B, Altrock PM, Wang L, Traulsen A: Universality of weak selection.

    Phys Rev E Stat Nonlin Soft Matter Phys 2010, 82:046106. PubMed Abstract | Publisher Full Text OpenURL

  28. van Baalen M, Rand DA: The unit of selection in viscous populations and the evolution of altruism.

    J Theor Biol 1998, 193(4):631-648. PubMed Abstract | Publisher Full Text OpenURL

  29. Lion S, Gandon S: Habitat saturation and the spatial evolutionary ecology of altruism.

    J Evol Biol 2009, 22(7):1487-1502. PubMed Abstract | Publisher Full Text OpenURL

  30. Wakano JY, Nowak MA, Hauert C: Spatial dynamics of ecological public goods.

    Proc Natl Acad Sci U S A 2009, 106(19):7910-7914. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Hallatschek O: Noise driven evolutionary waves.

    PLoS Comput Biol 2011, 7(3):e1002005. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Maruyama T: A Markov process of gene frequency change in a geographically structured population.

    Genetics 1974, 76(2):367-377. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Houchmandzadeh B, Vallade M: The fixation probability of a beneficial mutation in a geographically structured population.

    New J Phys 2011, 13:073020. Publisher Full Text OpenURL

  34. Traulsen A, Nowak MA: Evolution of cooperation by multilevel selection.

    Proc Natl Acad Sci U S A 2006, 103(29):10952-10955. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Traulsen A: Mathematics of kin- and group-selection: formally equivalent?

    Evolution 2010, 64(2):316-323. PubMed Abstract | Publisher Full Text OpenURL

  36. Melbinger A, Cremer J, Frey E: Evolutionary game theory in growing populations.

    Phys Rev Lett 2010, 105(17):178101. PubMed Abstract | Publisher Full Text OpenURL