Dynamics of Fst for the island model

F(st) is a measure of genetic differentiation in a subdivided population. Sewall Wright observed that F(st)=1/1+2Nm in a haploid diallelic infinite island model, where N is the effective population size of each deme and m is the migration rate. In demonstrating this result, Wright relied on the infinite size of the population. Natural populations are not infinite and therefore they change over time due to genetic drift. In a finite population, F(st) becomes a random variable that evolves over time. In this work we ask, given an initial population state, what are the dynamics of the mean and variance of F(st) under the finite island model? In application both of these quantities are critical in the evaluation of F(st) data. We show that after a time of order N generations the mean of F(st) is slightly biased below 1/1+2Nm. Further we show that the variance of F(st) is of order 1/d where d is the number of demes in the population. We introduce several new mathematical techniques to analyze coalescent genealogies in a dynamic setting.