April 2009 Archives
When a novel, adaptive mutation arises in a population it is more likely to go extinct than go to fixation. This is because rare alleles can by chance be lost from the population through accidental deaths, chromosomal segregation, or other forces behind genetic drift.
In 1955 Motoo Kimura used what’s known as diffusion theory to find formulas for calculating the probability that a novel allele with selection coefficient goes to fixation. For example, take a haploid, Wright-Fisher population of size that consists of only A individuals with fitness . Assume that one of those individuals mutates to B, such that the . From this setup, Kimura found that the the probability that B will eventually become fixed in the population and A goes extinct is approximately
under the right assumptions. For diploids with fitnesses , , and . This is
Sella and Hirch (2005) use intuition to modify Kimura’s results and derive some better and more useful alternatives for the above equations, which we use in our work for their nice properties. Sella and Hirch give the following equation in their paper:
where if the population is haploid and if the population is diploid “with multiplicative fitness within loci. …” Unfortunately, they did not specify the diploid model that they were using because there are two different ones used in the literature.
The first one is used above, , , and .
The second one is , , and .
As you can see, is not measuring the same thing in both models. In the first approach, ; in the second . This is important because accidentally mixing up the models will lead to erroneous results.
The difference between these two diploid models came up last week when we were applying Sella and Hirch (2005) to one of our projects. For diploids, we just couldn’t get their approximation to work (), and we suspected that they were using the second model, while we had read their paper to imply the first model. We looked back at their paper and realized that they neglect to specify exactly how they calculate and for diploids; we weren’t sure what they used.
I ended up working through their unspecified math to verify that in fact Sella and Hirch (2005) used the second model of diploid fitnesses, while we wanted to work with the first model. So here is a clarification of their equation, based on our way of thinking:
where if the population is haploid with fitnesses and , and if the population is diploid with fitnesses , , and .
- Sella and Hirch (2005) The application of statistical physics to evolutionary biology. PNAS 102:27 9541–9546.
The Panda’s Thumb is down right now, and as far as I can tell other servers at our Texas location are out as well. If Dallas had any storms over the night, then PT could have been kicked off the internet. Hopefully, they’ll be back up this evening.
Of course, as I post this it comes back up.