I saw a concept on the Internet that says "the strength of genetic drift is inversely proportional to the population size". I don't know why they are inversely proportional? Can somebody explain? Thank you all!
Answer
Plane Crash Analogy
4 people in a plane crash
In a small aeroplane, there are 2 people that wear a blue shirt and 2 people that wear a green shirt. The plane crashes, half of the people died. The 2 survivors are those wearing the green shirt… well, nothing so surprising!
400 people in a plane crash
In a very big aeroplane, there are 200 people that wear a blue shirt and 200 people that wear a green shirt. The plane crashes, half of the people died. The 200 survivors are those wearing green shirt… This is quite surprising!
Genetic Drift
The same logic applies to genetic drift. Genetic drift is caused by events that modify the reproductive success of individuals in a random way (independently of their genotype). We usually referred to this as random sampling. At each generation, individuals are randomly chosen to reproduce and some genotypes might just happen to be chosen more often than others at a given draw (=at a given generation). Genetic drift pushes the frequency of allele slightly away from what would be predicted. According to the Wright-Fisher model, the frequency of alleles (of a bi-allelic gene) in a haploid population (to make it easier) in the next time step is given by:
$$p' = \frac{p \cdot WA}{p \cdot WA + (1-p) \cdot Wa}$$
where $p$ is there frequency of an allele at time = $t$ and $p'$ is the frequency at time = $t+1$. $WA$ is the fitness of the genotype which frequency is $p$ and $Wa$ is the fitness of the genotype which frequency is $1-p$. If the population is infinite, the predictions of this equation are exactly correct.
Now if we say that meteorites fall and half of the individual get killed. The probability of getting killed by a meteorite obviously does not depend on genetic predisposition, it is a question of chance! If you look at a population of 1 million individuals, half of them having the genotype $A$, the other half having the genotype $a$. It is very unlikely that more than 60% of all individuals that get killed are of the same genotype. Therefore, the meteorites won't change much the frequency of the genotypes. If you look at a population of 4 individuals 2 are $a$ and 2 are $A$, 2 of them get killed by a meteorite. Well, you have a probability of one half that the two survivors are of the same genotype and that the genotypes frequency would have changed drastically.
Genetic drift refers to these changes in allele frequency which are due to random events (such as meteorites) and the strength of genetic drift indeed depends on the population size for probabilistic reasons. The greatest the population size, the lowest is the strength (or the relative importance) of genetic drift.
How to model genetic drift
There are three famous models of genetic drift that all lead to very similar expectations. I shall just name them here but I will not develop the underlying mathematics.
Wright-Fisher model of genetic drift
- Not to confuse with the Wright-Fisher model of selection written above
- It models genetic drift as a random sample of the previous generation and hence uses a binomial distribution.
Moran model
- It is based on a birth-death model (a type of Markov model).
Kimura's diffusion equation model
- It is an extension of the above two models for a case of continuous time.
No comments:
Post a Comment