The coalescent, Fisher-Wright-Moran model, and Lyapunov Equation

Marek Kimmel

Dept. of Statistics, Rice University, Houston, USA

with Adam Bobrowski, Ranajit Chakraborty and Ovide Arino


Abstract

We provide new results for the time-continuous Moran model with mutations of the general Markov-chain form. The matrix R(t) (possibly infinite) of the joint distributions of the types of a pair of alleles sampled from the population at time t, satisfies a matrix differential equation of the form dR(t)/dt=[Q^T R(t)+R(t)Q]-(1/(2N))R(t)+(1/(2N))P(t), where Q is the intensity matrix of the Markov chain, P(t) is its diagonalized probability distribution, and N is the effective population size. This is the Lyapunov differential equation, known in the control theory. Investigation of behavior of its solutions leads to consideration of tensor products of transition (Markov Semigroups). Semigroup theory methods allow to prove asymptotic results for the model, also in the case when the population size did not stay constant. Special cases of the model described above include stepwise mutation models with and without allele size constraints, and directional bias of mutations. Allele state changes caused by recombinatorial misalignment and more complex sequence conversion patterns also can be incorporated in this model. As an application, we investigate the influence of size constraints on evolutionary dynamics of microsatellite loci. We conclude that observed differences of variability between neutral and selective microsatellite loci are less likely due to allele-size constraints than to different mutation rates. Consequences for population dynamics of certain genetic diseases (such as fragile X or myotonic dystrophy) are discussed.