How do you simulate genetic drift?
Wright-Fisher population:
Genetic drift:
Decay in heterozygosity by drift:
H=1-G, where
H = frequency of the heterozygotes in the population
G = frequency of the homozygotes in a population
Ht = Ho(1-1/2N)to the power of t
Ht = Ho(1-1/2N)to the power of t
If N is really high (eg, 1,000, 000)
Mutation (u) and drift:
A problem in Neutral Theory:
Effective population size (Ne)
Harmonic mean:
Factors effecting Ne:
Different number of males than females: Ne = 4NmNf/Nm+Nf
eg) Elephant seals have harems, 5 males and 20 females (ie. Nc=25, Ne = 16)
Ne for Y: Ne=Nem/2
Ne for mitochondria: Ne(y)=Ne(mitochondria)
The X chromosome or haplodiploidy: Ne=9NmNf/4Nm+2nf
eg) a bee hive may have 20,000 bess, a singel queen and a dozen drones, Ne=9x12/4x12+2 = 2.16
Migration:
Gene flow:
Continent-Island model is one way flow from a very large stable population:
H=4Nu/(1+4Nu)