Evolution in a spatial continuum Drift, draft and structure Alison Etheridge University of Oxford Joint work with Nick Barton (Edinburgh) and Tom Kurtz (Wisconsin) New York, Sept. 2007 p.1
Kingman s Coalescent Kingman (1982) Neutral (haploid) population of constant size N Wright-Fisher model: new generation determined by multinomial sampling with equal weights New York, Sept. 2007 p.2
Kingman s Coalescent Kingman (1982) Neutral (haploid) population of constant size N Wright-Fisher model: new generation determined by multinomial sampling with equal weights Equivalently offspring choose their parent at random New York, Sept. 2007 p.2
Kingman s Coalescent Kingman (1982) Neutral (haploid) population of constant size N Wright-Fisher model: new generation determined by multinomial sampling with equal weights Equivalently offspring choose their parent at random Time in units of population size and let N New York, Sept. 2007 p.2
Kingman s Coalescent Kingman (1982) Neutral (haploid) population of constant size N Wright-Fisher model: new generation determined by multinomial sampling with equal weights MRCA Equivalently offspring choose their parent at random Time in units of population size and let N Coalescence rate ( ) k 2 Coalesent time Wright Fisher time New York, Sept. 2007 p.2
Forwards in time Two types a and A. p(t) = proportion of type a. New York, Sept. 2007 p.3
Forwards in time Two types a and A. p(t) = proportion of type a. Forwards in time, E[ p] = 0 (neutrality) E[( p) 2 ] = δtp(1 p) E[( p) 3 ] = O(δt) 2 dp t = p t (1 p t )dw t New York, Sept. 2007 p.3
Forwards in time Two types a and A. p(t) = proportion of type a. Forwards in time, E[ p] = 0 (neutrality) E[( p) 2 ] = δtp(1 p) E[( p) 3 ] = O(δt) 2 Backwards in time MRCA Coalesent time Wright Fisher time dp t = p t (1 p t )dw t New York, Sept. 2007 p.3
Forwards in time Two types a and A. p(t) = proportion of type a. Forwards in time, E[ p] = 0 (neutrality) E[( p) 2 ] = δtp(1 p) E[( p) 3 ] = O(δt) 2 Backwards in time MRCA Coalesent time Wright Fisher time dp t = p t (1 p t )dw t 1 dp τ = N p τ(1 p τ )dw τ, Coalescence rate 1 N ( ) k 2 New York, Sept. 2007 p.3
Basic observation Genetic diversity is orders of magnitude lower than expected from census numbers and genetic drift. Something else is going on... New York, Sept. 2007 p.4
Genetic hitchhiking When a selectively advantageous allele arises, it is either lost or it sweeps to fixation What about a neutral allele on the same chromosome? New York, Sept. 2007 p.5
Genetic hitchhiking When a selectively advantageous allele arises, it is either lost or it sweeps to fixation What about a neutral allele on the same chromosome? A D A C B D MEIOSIS B C Selected locus Neutral Locus Germ cell Gametes New York, Sept. 2007 p.5
... genetic draft Gillespie (2000) The frequency of the neutral allele linked to the new mutation will be boosted 1 selected frequency neutral u 1/N time New York, Sept. 2007 p.6
... genetic draft Gillespie (2000) The frequency of the neutral allele linked to the new mutation will be boosted Probability p, p u + p(1 u) Probability 1 p, p p(1 u) 1/N 1 frequency selected neutral u time New York, Sept. 2007 p.6
... genetic draft Gillespie (2000) The frequency of the neutral allele linked to the new mutation will be boosted Probability p, p u + p(1 u) Probability 1 p, p p(1 u) 1/N 1 frequency selected neutral u time E[ p] = 0, E[( p) 2 ] = p(1 p)e[u 2 ] New York, Sept. 2007 p.6
... genetic draft Gillespie (2000) The frequency of the neutral allele linked to the new mutation will be boosted Probability p, p u + p(1 u) Probability 1 p, p p(1 u) 1/N 1 frequency selected neutral u time E[ p] = 0, E[( p) 2 ] = p(1 p)e[u 2 ] E[( p) 3 ] = O(1) = multiple coalescences New York, Sept. 2007 p.6
Λ-coalescents Pitman (1999), Sagitov (1999) If there are currently p ancestral lineages, each transition involving j of them merging happens at rate β p,j = 1 Λ a finite measure on [0, 1] Kingman s coalescent, Λ = δ 0 0 u j 2 (1 u) p j Λ(du) New York, Sept. 2007 p.7
Forwards in time Bertoin & Le Gall (2003) The Λ-coalescent describes the genealogy of a sample from a population evolving according to a Λ-Fleming-Viot process. Poisson point process intensity dt u 2 Λ(du) individual sampled at random from population proportion u of population replaced by offspring of chosen individual New York, Sept. 2007 p.8
Spatial structure Kimura s stepping stone model dp i = j 1 m ij (p j p i )dt+ p i (1 p i )dw i N e System of interacting Wright-Fisher diffusions New York, Sept. 2007 p.9
Spatial structure Kimura s stepping stone model dp i = j 1 m ij (p j p i )dt+ p i (1 p i )dw i N e System of interacting Wright-Fisher diffusions Genealogy described by system of coalescing random walks New York, Sept. 2007 p.9
Evolution in a spatial continuum? For many biological populations it is more natural to consider a spatial continuum Can we replace the stepping stone model by a stochastic pde? New York, Sept. 2007 p.10
Evolution in a spatial continuum? For many biological populations it is more natural to consider a spatial continuum Can we replace the stepping stone model by a stochastic pde? dp = 1 2 pdt + p(1 p)dw W a space-time white noise New York, Sept. 2007 p.10
Evolution in a spatial continuum? For many biological populations it is more natural to consider a spatial continuum Can we replace the stepping stone model by a stochastic pde? dp = 1 2 pdt + p(1 p)dw W a space-time white noise In two dimensions the equation has no solution Diffusive rescaling leads to the heat equation But anyway local populations are finite New York, Sept. 2007 p.10
Another basic observation Real populations experience large scale fluctuations in which the movement and reproductive success of many individuals are correlated New York, Sept. 2007 p.11
An individual based model Start with Poisson point process intensity λdx New York, Sept. 2007 p.12
An individual based model Start with Poisson point process intensity λdx At rate µ(dr) dx dt throw down ball centre x, radius r. New York, Sept. 2007 p.12
An individual based model Start with Poisson point process intensity λdx At rate µ(dr) dx dt throw down ball centre x, radius r. Each individual in region dies with probability u New York, Sept. 2007 p.12
An individual based model Start with Poisson point process intensity λdx At rate µ(dr) dx dt throw down ball centre x, radius r. Each individual in region dies with probability u New individuals born according to a Poisson λudx New York, Sept. 2007 p.12
An individual based model Start with Poisson point process intensity λdx At rate µ(dr) dx dt throw down ball centre x, radius r. Each individual in region dies with probability u New individuals born according to a Poisson λudx New York, Sept. 2007 p.12
An individual based model Start with Poisson point process intensity λdx At rate µ(dr) dx dt throw down ball centre x, radius r. Each individual in region dies with probability u x r New individuals born according to a Poisson λudx New York, Sept. 2007 p.12
An individual based model Start with Poisson point process intensity λdx At rate µ(dr) dx dt throw down ball centre x, radius r. Each individual in region dies with probability u x r New individuals born according to a Poisson λudx New York, Sept. 2007 p.12
An individual based model Start with Poisson point process intensity λdx At rate µ(dr) dx dt throw down ball centre x, radius r. Each individual in region dies with probability u x r New individuals born according to a Poisson λudx New York, Sept. 2007 p.12
Two limits Let the local population density λ New York, Sept. 2007 p.13
Two limits Let the local population density λ Rescale space and time to investigate large scale effects New York, Sept. 2007 p.13
The spatial Λ-Fleming-Viot process State {ρ(t, x, ) M 1 (K), x R 2, t 0}. New York, Sept. 2007 p.14
The spatial Λ-Fleming-Viot process State {ρ(t, x, ) M 1 (K), x R 2, t 0}. π Poisson point process rate µ(dr) dx dt. New York, Sept. 2007 p.14
The spatial Λ-Fleming-Viot process State {ρ(t, x, ) M 1 (K), x R 2, t 0}. π Poisson point process rate µ(dr) dx dt. For each r > 0, ν r (du) M 1 ([0, 1]). New York, Sept. 2007 p.14
The spatial Λ-Fleming-Viot process State {ρ(t, x, ) M 1 (K), x R 2, t 0}. π Poisson point process rate µ(dr) dx dt. For each r > 0, ν r (du) M 1 ([0, 1]). Dynamics: for each (t, x, r) π, New York, Sept. 2007 p.14
The spatial Λ-Fleming-Viot process State {ρ(t, x, ) M 1 (K), x R 2, t 0}. π Poisson point process rate µ(dr) dx dt. For each r > 0, ν r (du) M 1 ([0, 1]). Dynamics: for each (t, x, r) π, u ν r (du) New York, Sept. 2007 p.14
The spatial Λ-Fleming-Viot process State {ρ(t, x, ) M 1 (K), x R 2, t 0}. π Poisson point process rate µ(dr) dx dt. For each r > 0, ν r (du) M 1 ([0, 1]). Dynamics: for each (t, x, r) π, u ν r (du) z U(B r (x)) New York, Sept. 2007 p.14
The spatial Λ-Fleming-Viot process State {ρ(t, x, ) M 1 (K), x R 2, t 0}. π Poisson point process rate µ(dr) dx dt. For each r > 0, ν r (du) M 1 ([0, 1]). Dynamics: for each (t, x, r) π, u ν r (du) z U(B r (x)) k ρ(t, z, ). New York, Sept. 2007 p.14
The spatial Λ-Fleming-Viot process State {ρ(t, x, ) M 1 (K), x R 2, t 0}. π Poisson point process rate µ(dr) dx dt. For each r > 0, ν r (du) M 1 ([0, 1]). Dynamics: for each (t, x, r) π, u ν r (du) z U(B r (x)) k ρ(t, z, ). For all y B r (x), ρ(t, y, ) = (1 u)ρ(t, y, ) + uδ k. New York, Sept. 2007 p.14
Conditions (1) ρ(t, x, ) experiences jump of size u A [0, 1] at rate πr 2 ν r (u)µ(dr). (0, ] A New York, Sept. 2007 p.15
Conditions (1) ρ(t, x, ) experiences jump of size u A [0, 1] at rate πr 2 ν r (u)µ(dr). (0, ] A Λ(du) = (0, ) u 2 r 2 ν r (du)µ(dr) M F ([0, 1]) New York, Sept. 2007 p.15
Conditions (1) ρ(t, x, ) experiences jump of size u A [0, 1] at rate πr 2 ν r (u)µ(dr). (0, ] A Λ(du) = (0, ) u 2 r 2 ν r (du)µ(dr) M F ([0, 1]) A single ancestral lineage evolves in series of jumps with intensity dt ( x /2, ) [0,1] L r (x) πr 2 on R + R 2 where L r (x) = B r (0) B r (x). u ν r (du)µ(dr)dx New York, Sept. 2007 p.15
Conditions (2) ( (1 x 2 ) R 2 ( x /2, ) [0,1] L r (x) πr 2 u ν r (du)µ(dr) ) dx <. New York, Sept. 2007 p.16
Conditions (2) ( (1 x 2 ) R 2 ( x /2, ) [0,1] L r (x) πr 2 u ν r (du)µ(dr) ) dx <. Two lineages currently at separation y R 2 coalesce at instantaneous rate ( ) L r (y) u 2 ν r (du) µ(dr). ( y /2, ) [0,1] New York, Sept. 2007 p.16
Generalisations Replace R 2 by an arbitrary Polish space New York, Sept. 2007 p.17
Generalisations Replace R 2 by an arbitrary Polish space Choose a Poisson number of parents at each reproduction event New York, Sept. 2007 p.17
Generalisations Replace R 2 by an arbitrary Polish space Choose a Poisson number of parents at each reproduction event Choose spatial position of parents non-uniformly. New York, Sept. 2007 p.17
Generalisations Replace R 2 by an arbitrary Polish space Choose a Poisson number of parents at each reproduction event Choose spatial position of parents non-uniformly. Impose spatial motion of individuals not linked directly to the reproduction events. New York, Sept. 2007 p.17
Generalisations Replace R 2 by an arbitrary Polish space Choose a Poisson number of parents at each reproduction event Choose spatial position of parents non-uniformly. Impose spatial motion of individuals not linked directly to the reproduction events. Instead of replacing a portion u of individuals from a ball centred on x, replace individuals sampled according to some distribution (e.g. Gaussian) centred on x. New York, Sept. 2007 p.17
Generalisations Replace R 2 by an arbitrary Polish space Choose a Poisson number of parents at each reproduction event Choose spatial position of parents non-uniformly. Impose spatial motion of individuals not linked directly to the reproduction events. Instead of replacing a portion u of individuals from a ball centred on x, replace individuals sampled according to some distribution (e.g. Gaussian) centred on x.... and many more. New York, Sept. 2007 p.17