Some mathematical models from population genetics
|
|
- Brian Moore
- 6 years ago
- Views:
Transcription
1 Some mathematical models from population genetics 5: Muller s ratchet and the rate of adaptation Alison Etheridge University of Oxford joint work with Peter Pfaffelhuber (Vienna), Anton Wakolbinger (Frankfurt) Charles Cuthbertson (Oxford), Feng Yu (Bristol) New York, Sept. 07 p.1
2 Muller s ratchet Asexually reproducing population. All mutations deleterious. Chromosomes passed down as indivisible blocks. New York, Sept. 07 p.2
3 Muller s ratchet Asexually reproducing population. All mutations deleterious. Chromosomes passed down as indivisible blocks. Number of deleterious mutations accumulated along any ancestral line can only increase. New York, Sept. 07 p.2
4 Muller s ratchet Asexually reproducing population. All mutations deleterious. Chromosomes passed down as indivisible blocks. Number of deleterious mutations accumulated along any ancestral line can only increase. When everyone in the current best class has accumulated at least one additional bad mutation the ratchet clicks. New York, Sept. 07 p.2
5 Muller s ratchet Asexually reproducing population. All mutations deleterious. Chromosomes passed down as indivisible blocks. Number of deleterious mutations accumulated along any ancestral line can only increase. When everyone in the current best class has accumulated at least one additional bad mutation the ratchet clicks. How many generations will it take for an asexually reproducing population to lose its best class? New York, Sept. 07 p.2
6 Haigh s model Wright-Fisher model: Individuals in (t + 1)st generation select parent at random from generation t. Probability individual which has accumulated k mutations is selected as parent proportional to relative fitness (1 s) k. Number of mutations carried by offspring then k + J, where J Poiss(λ) (independent). New York, Sept. 07 p.3
7 Type frequencies: x(t) = (x k (t)) k=0,1,... Nx k (t) = No. of individuals carrying exactly k mutations at time t. New York, Sept. 07 p.4
8 Type frequencies: x(t) = (x k (t)) k=0,1,... Nx k (t) = No. of individuals carrying exactly k mutations at time t. H an N 0 -valued random variable. P[H = k] (1 s) k x k (t), New York, Sept. 07 p.4
9 Type frequencies: x(t) = (x k (t)) k=0,1,... Nx k (t) = No. of individuals carrying exactly k mutations at time t. H an N 0 -valued random variable. J Poiss(λ) independent of H. P[H = k] (1 s) k x k (t), New York, Sept. 07 p.4
10 Type frequencies: x(t) = (x k (t)) k=0,1,... Nx k (t) = No. of individuals carrying exactly k mutations at time t. H an N 0 -valued random variable. P[H = k] (1 s) k x k (t), J Poiss(λ) independent of H. K 1, K 2,..., K N independent copies of H + J. Random type frequencies in next generation are X k (t + 1) = 1 N #{i : K i = k}. New York, Sept. 07 p.4
11 Infinite populations As N, LLN x(t + 1) = E x(t) [X(t + 1)]. New York, Sept. 07 p.5
12 Infinite populations As N, LLN x(t + 1) = E x(t) [X(t + 1)]. Suppose x(t) Poiss(α). P[H = k] (1 s) k x k = (1 s) k αk e α. k! Then H Poiss(α(1 s)), J Poiss(λ), so H + J Poiss(α(1 s) + λ). New York, Sept. 07 p.5
13 Infinite populations As N, LLN x(t + 1) = E x(t) [X(t + 1)]. Suppose x(t) Poiss(α). P[H = k] (1 s) k x k = (1 s) k αk e α. k! Then H Poiss(α(1 s)), J Poiss(λ), so H + J Poiss(α(1 s) + λ). Poisson weights Poisson weights. New York, Sept. 07 p.5
14 For every initial condition with x 0 > 0, the solution to the deterministic dynamics converges as t to the stationary point π := Poiss(λ/s). New York, Sept. 07 p.6
15 Back to finite populations Write k for the number of mutations carried by individuals in fittest class. New York, Sept. 07 p.7
16 Back to finite populations Write k for the number of mutations carried by individuals in fittest class. forms a recurrent Markov chain. Y := (Y k ) k=0,1,... := (X k +k ) k=0,1,2... New York, Sept. 07 p.7
17 Back to finite populations Write k for the number of mutations carried by individuals in fittest class. Y := (Y k ) k=0,1,... := (X k +k ) k=0,1,2... forms a recurrent Markov chain. Condition on Y(t) = y(t). Size of new best class, y 0 (t + 1) Binom(N, p 0 (t)), with p 0 (t) probability of sampling parent from best class and not acquiring any additional mutations: p 0 (t) = y 0(t) W (t) e λ, W (t) = y i (t)(1 s) i. i=0 New York, Sept. 07 p.7
18 Back to finite populations Write k for the number of mutations carried by individuals in fittest class. Y := (Y k ) k=0,1,... := (X k +k ) k=0,1,2... forms a recurrent Markov chain. Condition on Y(t) = y(t). Size of new best class, y 0 (t + 1) Binom(N, p 0 (t)), with p 0 (t) probability of sampling parent from best class and not acquiring any additional mutations: p 0 (t) = y 0(t) W (t) e λ, W (t) = y i (t)(1 s) i. i=0 Evolution of best class determined by W (t), the mean fitness in the population. New York, Sept. 07 p.7
19 Elements of Haigh s analysis Immediately after a click, the type frequencies are π := 1 1 π 0 (π 1, π 2,...). New York, Sept. 07 p.8
20 Elements of Haigh s analysis Immediately after a click, the type frequencies are π := 1 1 π 0 (π 1, π 2,...). Phase one: deterministic dynamical system dominates, decaying exponentially fast towards its equilibrium New York, Sept. 07 p.8
21 Elements of Haigh s analysis Immediately after a click, the type frequencies are π := 1 1 π 0 (π 1, π 2,...). Phase one: deterministic dynamical system dominates, decaying exponentially fast towards its equilibrium Phase two: the bulk of the population changes only slowly. Mean fitness assumed constant and then No. of individuals in best class approximated by Galton-Watson branching process with Poisson offspring distribution. New York, Sept. 07 p.8
22 The Fleming Viot diffusion Relevant for applications are: large N, small s, small λ. New York, Sept. 07 p.9
23 The Fleming Viot diffusion Relevant for applications are: large N, small s, small λ. Dynamics then captured by dx k = j s(j k)x j X k + λ(x k 1 X k ) dt + j k 1 N X jx k dw jk, k = 0, 1, 2,... where X 1 = 0 and (W jk ) j>k array of independent Brownian motions, W kj := W jk. New York, Sept. 07 p.9
24 The Fleming Viot diffusion Relevant for applications are: large N, small s, small λ. Dynamics then captured by dx k = j s(j k)x j X k + λ(x k 1 X k ) dt + j k 1 N X jx k dw jk, k = 0, 1, 2,... where X 1 = 0 and (W jk ) j>k array of independent Brownian motions, W kj := W jk. As before Y k = X k +k, dy 0 = s(m 1 (Y) λ)y 0 (t)dt + 1 N Y 0(1 Y 0 )dw 0, M 1 (Y) = j jy j. New York, Sept. 07 p.9
25 Infinite population limit with x 1 = 0. dx k = (s(m 1 (x) k) λ)x k + λx k 1 ) dt, k = 0, 1, 2,... New York, Sept. 07 p.10
26 Infinite population limit dx k = (s(m 1 (x) k) λ)x k + λx k 1 ) dt, k = 0, 1, 2,... with x 1 = 0. Transform into system of equations for cumulants: log x k e ξk = k=0 k=1 κ k ( ξ) k k!. Assume x 0 > 0 and set κ 0 = log x 0. Then κ k = sκ k+1 + λ, k = 0, 1, 2,... New York, Sept. 07 p.10
27 System can be solved. In particular, κ 1 (t) = kx k (t) = ξ log k=0 k=0 x k (0)e ξk ξ=st + λ s (1 e st ). New York, Sept. 07 p.11
28 System can be solved. In particular, κ 1 (t) = kx k (t) = ξ log k=0 k=0 x k (0)e ξk ξ=st + λ s (1 e st ). Notice the exponential decay towards the equilibrium vaue of λ/s. The time τ = log(λ/s) s corresponds exactly to the end of Haigh s phase one. New York, Sept. 07 p.11
29 Approximations dy 0 = s(m 1 (Y) λ)y 0 (t)dt + 1 N Y 0(1 Y 0 )dw 0. Cannot solve for M 1 (Y). Instead seek a good approximation of M 1 given Y 0. New York, Sept. 07 p.12
30 Approximations dy 0 = s(m 1 (Y) λ)y 0 (t)dt + 1 N Y 0(1 Y 0 )dw 0. Cannot solve for M 1 (Y). Instead seek a good approximation of M 1 given Y 0. Simulations suggest a good fit to a linear relationship between Y 0 and M 1. New York, Sept. 07 p.12
31 Extending Haigh s approach Haigh assumes at click times π 0 distributed evenly over other classes. Suppose now that this holds in between click times too: given Y 0 approximate state of system by the PPA (Poisson Profile Approximation) Π(Y 0 ) = ( Y 0, 1 Y ) 0 (π 1, π 2,...). 1 π 0 New York, Sept. 07 p.13
32 Extending Haigh s approach Haigh assumes at click times π 0 distributed evenly over other classes. Suppose now that this holds in between click times too: given Y 0 approximate state of system by the PPA (Poisson Profile Approximation) Π(Y 0 ) = ( Y 0, 1 Y ) 0 (π 1, π 2,...). 1 π 0 Estimate M 1 not from PPA but from relaxed PPA obtained by evolving PPA according to the deterministic dynamical system for time Aτ := A log(λ/s)/s. This gives M 1 = θ + η e η 1 ( 1 Y ) 0, η = (λ/s) 1 A. π 0 New York, Sept. 07 p.13
33 Three one dimensional diffusions Substituting in the one-dimensional diffusion approximation for Y 0 gives: A small, dy 0 = λ(π 0 Y 0 )Y 0 dt + ( A = 1, dy 0 = 0.58s 1 Y ) 0 π 0 A large, ( dy 0 = s 1 Y ) 0 Y 0 dt + π 0 1 N Y 0dW, Y 0 dt + 1 N Y 0 dw, 1 N Y 0 dw, New York, Sept. 07 p.14
34 A rescaling Z(t) = 1 π 0 Y 0 (Nπ 0 t). Set γ = Nλ Ns log(nλ). New York, Sept. 07 p.15
35 A rescaling Z(t) = 1 π 0 Y 0 (Nπ 0 t). Set γ = Nλ Ns log(nλ). A small, dz = (Nλ) 1 2γ (1 Z)Zdt + ZdW, A = 1, 1 dz = 0.58 γ log(nλ) (Nλ)1 γ (1 Z)Zdt + Z dw, A large, dz = 1 γ log(nλ) (Nλ)1 γ (1 Z)Zdt + Z dw. New York, Sept. 07 p.15
36 Implications A small fast clicking: dz = (Nλ) 1 2γ (1 Z)Zdt + ZdW. The ratchet never clicks for γ < 1/2. New York, Sept. 07 p.16
37 Implications A small fast clicking: dz = (Nλ) 1 2γ (1 Z)Zdt + ZdW. The ratchet never clicks for γ < 1/2. A = 1 (A large) moderate clicking (slow clicking): dz = 0.58(1) 1 γ log(nλ) (Nλ)1 γ (1 Z)Zdt + Z dw. New York, Sept. 07 p.16
38 Implications A small fast clicking: dz = (Nλ) 1 2γ (1 Z)Zdt + ZdW. The ratchet never clicks for γ < 1/2. A = 1 (A large) moderate clicking (slow clicking): dz = 0.58(1) 1 γ log(nλ) (Nλ)1 γ (1 Z)Zdt + Z dw. In order for γ log(nλ) (Nλ)1 γ to be > 5, γ Nλ New York, Sept. 07 p.16
39 Rule of thumb For biologically realistic parameters, transition from no clicks to moderate clicks (on evolutionary timescale) around γ = 0.5. New York, Sept. 07 p.17
40 Rule of thumb For biologically realistic parameters, transition from no clicks to moderate clicks (on evolutionary timescale) around γ = 0.5. The rate of the ratchet is of the order N γ 1 λ γ for γ (1/2, 1), whereas it is exponentially slow in (Nλ) 1 γ for γ < 1/2. New York, Sept. 07 p.17
41 Simulations New York, Sept. 07 p.18
42 New York, Sept. 07 p.19
43 Purely beneficial mutations We saw before that the probability of an isolated selected sweep 2σ. New York, Sept. 07 p.20
44 Purely beneficial mutations We saw before that the probability of an isolated selected sweep 2σ. What about overlapping sweeps? New York, Sept. 07 p.20
45 Purely beneficial mutations We saw before that the probability of an isolated selected sweep 2σ. What about overlapping sweeps? Interference reduces the chance of fixation as beneficial mutations compete with one another. New York, Sept. 07 p.20
46 Purely beneficial mutations We saw before that the probability of an isolated selected sweep 2σ. What about overlapping sweeps? Interference reduces the chance of fixation as beneficial mutations compete with one another. Is there a limit to the rate of adaptation? New York, Sept. 07 p.20
47 A mixture of mutations All mutations equal strength so ith individual s fitness characterized by net number, X i, of beneficial mutations. New York, Sept. 07 p.21
48 A mixture of mutations All mutations equal strength so ith individual s fitness characterized by net number, X i, of beneficial mutations. Mutation: For each individual i a mutation event occurs at rate µ. With probability 1 q, X i changes to X i 1 and with probability q, X i changes to X i + 1. Selection: For each pair of individuals (i, j), at rate σ N (X i X j ) +, individual i replaces individual j. Resampling: For each pair of individuals (i, j), at rate 1 N, individual i replaces individual j. New York, Sept. 07 p.21
49 A mixture of mutations All mutations equal strength so ith individual s fitness characterized by net number, X i, of beneficial mutations. Mutation: For each individual i a mutation event occurs at rate µ. With probability 1 q, X i changes to X i 1 and with probability q, X i changes to X i + 1. Selection: For each pair of individuals (i, j), at rate σ N (X i X j ) +, individual i replaces individual j. Resampling: For each pair of individuals (i, j), at rate 1 N, individual i replaces individual j. Technical modification: suppress mutations which would make the width of the population > L N 1/4. New York, Sept. 07 p.21
50 ..., P 0 (t),..., P k (t),... - Proportion of individuals with k mutations at time t k max (k min ) - type of the fittest (least fit) individual New York, Sept. 07 p.22
51 ..., P 0 (t),..., P k (t),... - Proportion of individuals with k mutations at time t k max (k min ) - type of the fittest (least fit) individual dp k = µ k (P ) dt + σ l Z(k l)p k P l dt + dm P k = [ µ k (P ) + σ(k m(p ))P k ] dt + dm P k µ k (P ) µ (qp k 1 P k + (1 q)p k+1 ) New York, Sept. 07 p.22
52 ..., P 0 (t),..., P k (t),... - Proportion of individuals with k mutations at time t k max (k min ) - type of the fittest (least fit) individual dp k = µ k (P ) dt + σ l Z(k l)p k P l dt + dm P k = [ µ k (P ) + σ(k m(p ))P k ] dt + dm P k µ k (P ) µ (qp k 1 P k + (1 q)p k+1 ) M P is a martingale with [ ] M P 2µ k (t) N t + 1 N t 0 (2 + σ(k l) + + σ(l k) + )P k (s)p l (s) ds l Z New York, Sept. 07 p.22
53 Moments Mean fitness m(p ) = k kp k satisfies dm(p ) = ( µ(p ) + σc 2 (P )) dt + dm P,m µ(p ) µ(2q 1) New York, Sept. 07 p.23
54 Moments Mean fitness m(p ) = k kp k satisfies dm(p ) = ( µ(p ) + σc 2 (P )) dt + dm P,m µ(p ) µ(2q 1) Ignoring mutation terms, centred moments c k = k (k m(p))n P k satisfy dc 2 σc 3 dt + small noise terms dc 3 σ(c 4 3c 2 c 2 ) dt + small noise terms dc 4 σ(c 5 4c 3 c 2 ) dt + small noise terms dc 5 σ(c 6 5c 4 c 2 ) dt + small noise terms New York, Sept. 07 p.23
55 Heuristics Centred process has a stationary distribution. New York, Sept. 07 p.24
56 Heuristics Centred process has a stationary distribution. Approximate by a deterministic distribution. New York, Sept. 07 p.24
57 Heuristics Centred process has a stationary distribution. Approximate by a deterministic distribution. Equations for centred moments give: 0 σc 3 dt 0 σ(c 4 3c 2 c 2 ) 0 σ(c 5 4c 3 c 2 ) 0 σ(c 6 5c 4 c 2 )... New York, Sept. 07 p.24
58 Heuristics Centred process has a stationary distribution. Approximate by a deterministic distribution. Equations for centred moments give: 0 σc 3 dt 0 σ(c 4 3c 2 c 2 ) 0 σ(c 5 4c 3 c 2 ) 0 σ(c 6 5c 4 c 2 )... Stationary distribution approximately Gaussian. New York, Sept. 07 p.24
59 Suppose stationary distribution Gaussian with mean m(p ), variance b 2. New York, Sept. 07 p.25
60 Suppose stationary distribution Gaussian with mean m(p ), variance b 2. Front is at K where 1 2πb e K2 /2b 2 = 1 N K b 2 log N New York, Sept. 07 p.25
61 Suppose stationary distribution Gaussian with mean m(p ), variance b 2. Front is at K where 1 2πb e K2 /2b 2 = 1 N K b 2 log N If there is a single individual at m(p ) + K at time t = 0, how long until there is an individual at m(p ) + K + 1? New York, Sept. 07 p.25
62 Suppose stationary distribution Gaussian with mean m(p ), variance b 2. Front is at K where 1 2πb e K2 /2b 2 = 1 N K b 2 log N If there is a single individual at m(p ) + K at time t = 0, how long until there is an individual at m(p ) + K + 1? Ignoring beneficial mutations occurring to individuals at m(p ) + K 1, until the front advances Z(t) e (σk (1 q)µ)t. New York, Sept. 07 p.25
63 More heuristics Front advances as soon as individual of Z(t) accumulates a beneficial mutation. New York, Sept. 07 p.26
64 More heuristics Front advances as soon as individual of Z(t) accumulates a beneficial mutation. The probability that the front doesn t advance by time t is exp { t qµ 0 } Z(s) ds exp { } qµ σk (1 q)µ (e(σk (1 q)µ)t 1) New York, Sept. 07 p.26
65 More heuristics Front advances as soon as individual of Z(t) accumulates a beneficial mutation. The probability that the front doesn t advance by time t is exp { t qµ 0 } Z(s) ds exp Average time for front to advance by one is { } qµ σk (1 q)µ (e(σk (1 q)µ)t 1) 1 σk (1 q)µ log(σk (1 q)µ). New York, Sept. 07 p.26
66 More heuristics Front advances as soon as individual of Z(t) accumulates a beneficial mutation. The probability that the front doesn t advance by time t is exp { t qµ 0 } Z(s) ds exp Average time for front to advance by one is { } qµ σk (1 q)µ (e(σk (1 q)µ)t 1) 1 σk (1 q)µ log(σk (1 q)µ). But wave speed is µ(2q 1) + σc 2 (P ) µ(2q 1) + σb 2 µ(2q 1) + σk2 2 log N. New York, Sept. 07 p.26
67 Conclusion Consistency condition: σk (1 q)µ log(σk (1 q)µ) = µ(2q 1) + σk2 2 log N. New York, Sept. 07 p.27
68 Conclusion Consistency condition: σk (1 q)µ log(σk (1 q)µ) = µ(2q 1) + σk2 2 log N. For large K, this approximately reduces to K log(σk) = 2 log N. New York, Sept. 07 p.27
69 Conclusion Consistency condition: σk (1 q)µ log(σk (1 q)µ) = µ(2q 1) + σk2 2 log N. For large K, this approximately reduces to K log(σk) = 2 log N. If K = log N, then LHS > RHS; K = log 1 δ N, δ > 0, then LHS < RHS. New York, Sept. 07 p.27
70 Conclusion Consistency condition: σk (1 q)µ log(σk (1 q)µ) = µ(2q 1) + σk2 2 log N. For large K, this approximately reduces to K log(σk) = 2 log N. If K = log N, then LHS > RHS; K = log 1 δ N, δ > 0, then LHS < RHS. So K between log N and any fractional power of log N rate of adaptation, of order between log N and any fractional power of log N. New York, Sept. 07 p.27
71 Rigorous result Theorem. If q > 0, then for any β > 0, there exists a positive constant c µ,σ,q such that if N is sufficiently large. E π [m(1)] E π [c 2 ] c µ,σ,q log 1 β N New York, Sept. 07 p.28
72 Simulations With µ = 0.01, q = 0.01, σ = 0.01, N = 1000, 2000, 5000, 10000, First row: mean; second row: variance mean 2 mean mean 0.5 mean mean time time time time time var 1 var var var var time time time time time New York, Sept. 07 p.29
73 6 x Adaptation Rate N Adaptation rate against population size. From top to bottom, q = 4%, 2%, 1%, 0.2%, µ = 0.01, σ = New York, Sept. 07 p.30
arxiv: v1 [math.pr] 18 Sep 2007
arxiv:0709.2775v [math.pr] 8 Sep 2007 How often does the ratchet click? Facts, heuristics, asymptotics. by A. Etheridge, P. Pfaffelhuber and A. Wakolbinger University of Oxford, Ludwig-Maximilians University
More informationProperties of an infinite dimensional EDS system : the Muller s ratchet
Properties of an infinite dimensional EDS system : the Muller s ratchet LATP June 5, 2011 A ratchet source : wikipedia Plan 1 Introduction : The model of Haigh 2 3 Hypothesis (Biological) : The population
More informationEvolution in a spatial continuum
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
More informationThe Λ-Fleming-Viot process and a connection with Wright-Fisher diffusion. Bob Griffiths University of Oxford
The Λ-Fleming-Viot process and a connection with Wright-Fisher diffusion Bob Griffiths University of Oxford A d-dimensional Λ-Fleming-Viot process {X(t)} t 0 representing frequencies of d types of individuals
More informationModelling populations under fluctuating selection
Modelling populations under fluctuating selection Alison Etheridge With Aleksander Klimek (Oxford) and Niloy Biswas (Harvard) The simplest imaginable model of inheritance A population of fixed size, N,
More informationThe mathematical challenge. Evolution in a spatial continuum. The mathematical challenge. Other recruits... The mathematical challenge
The mathematical challenge What is the relative importance of mutation, selection, random drift and population subdivision for standing genetic variation? Evolution in a spatial continuum Al lison Etheridge
More informationThe tree-valued Fleming-Viot process with mutation and selection
The tree-valued Fleming-Viot process with mutation and selection Peter Pfaffelhuber University of Freiburg Joint work with Andrej Depperschmidt and Andreas Greven Population genetic models Populations
More informationThe effects of a weak selection pressure in a spatially structured population
The effects of a weak selection pressure in a spatially structured population A.M. Etheridge, A. Véber and F. Yu CNRS - École Polytechnique Motivations Aim : Model the evolution of the genetic composition
More informationThe Wright-Fisher Model and Genetic Drift
The Wright-Fisher Model and Genetic Drift January 22, 2015 1 1 Hardy-Weinberg Equilibrium Our goal is to understand the dynamics of allele and genotype frequencies in an infinite, randomlymating population
More informationDynamics of the evolving Bolthausen-Sznitman coalescent. by Jason Schweinsberg University of California at San Diego.
Dynamics of the evolving Bolthausen-Sznitman coalescent by Jason Schweinsberg University of California at San Diego Outline of Talk 1. The Moran model and Kingman s coalescent 2. The evolving Kingman s
More informationHow robust are the predictions of the W-F Model?
How robust are the predictions of the W-F Model? As simplistic as the Wright-Fisher model may be, it accurately describes the behavior of many other models incorporating additional complexity. Many population
More informationHomework # , Spring Due 14 May Convergence of the empirical CDF, uniform samples
Homework #3 36-754, Spring 27 Due 14 May 27 1 Convergence of the empirical CDF, uniform samples In this problem and the next, X i are IID samples on the real line, with cumulative distribution function
More informationProblems on Evolutionary dynamics
Problems on Evolutionary dynamics Doctoral Programme in Physics José A. Cuesta Lausanne, June 10 13, 2014 Replication 1. Consider the Galton-Watson process defined by the offspring distribution p 0 =
More informationarxiv: v4 [math.pr] 15 Oct 2010
The Annals of Applied Probability 21, Vol. 2, No. 3, 978 14 DOI: 1.1214/9-AAP645 c Institute of Mathematical Statistics, 21 ASYMPTOTIC BEHAVIOR OF THE RATE OF ADAPTATION arxiv:78.3453v4 [math.pr] 15 Oct
More informationPopulation Genetics: a tutorial
: a tutorial Institute for Science and Technology Austria ThRaSh 2014 provides the basic mathematical foundation of evolutionary theory allows a better understanding of experiments allows the development
More informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationA representation for the semigroup of a two-level Fleming Viot process in terms of the Kingman nested coalescent
A representation for the semigroup of a two-level Fleming Viot process in terms of the Kingman nested coalescent Airam Blancas June 11, 017 Abstract Simple nested coalescent were introduced in [1] to model
More informationConvergence to equilibrium for rough differential equations
Convergence to equilibrium for rough differential equations Samy Tindel Purdue University Barcelona GSE Summer Forum 2017 Joint work with Aurélien Deya (Nancy) and Fabien Panloup (Angers) Samy T. (Purdue)
More informationSDE Coefficients. March 4, 2008
SDE Coefficients March 4, 2008 The following is a summary of GARD sections 3.3 and 6., mainly as an overview of the two main approaches to creating a SDE model. Stochastic Differential Equations (SDE)
More informationFrequency Spectra and Inference in Population Genetics
Frequency Spectra and Inference in Population Genetics Although coalescent models have come to play a central role in population genetics, there are some situations where genealogies may not lead to efficient
More informationDiffusion Models in Population Genetics
Diffusion Models in Population Genetics Laura Kubatko kubatko.2@osu.edu MBI Workshop on Spatially-varying stochastic differential equations, with application to the biological sciences July 10, 2015 Laura
More informationTime Series 3. Robert Almgren. Sept. 28, 2009
Time Series 3 Robert Almgren Sept. 28, 2009 Last time we discussed two main categories of linear models, and their combination. Here w t denotes a white noise: a stationary process with E w t ) = 0, E
More informationMathematical models in population genetics II
Mathematical models in population genetics II Anand Bhaskar Evolutionary Biology and Theory of Computing Bootcamp January 1, 014 Quick recap Large discrete-time randomly mating Wright-Fisher population
More informationUNIVERSITY OF LONDON IMPERIAL COLLEGE LONDON
UNIVERSITY OF LONDON IMPERIAL COLLEGE LONDON BSc and MSci EXAMINATIONS (MATHEMATICS) MAY JUNE 23 This paper is also taken for the relevant examination for the Associateship. M3S4/M4S4 (SOLUTIONS) APPLIED
More informationON COMPOUND POISSON POPULATION MODELS
ON COMPOUND POISSON POPULATION MODELS Martin Möhle, University of Tübingen (joint work with Thierry Huillet, Université de Cergy-Pontoise) Workshop on Probability, Population Genetics and Evolution Centre
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationURN MODELS: the Ewens Sampling Lemma
Department of Computer Science Brown University, Providence sorin@cs.brown.edu October 3, 2014 1 2 3 4 Mutation Mutation: typical values for parameters Equilibrium Probability of fixation 5 6 Ewens Sampling
More informationUniversal examples. Chapter The Bernoulli process
Chapter 1 Universal examples 1.1 The Bernoulli process First description: Bernoulli random variables Y i for i = 1, 2, 3,... independent with P [Y i = 1] = p and P [Y i = ] = 1 p. Second description: Binomial
More informationLAN property for ergodic jump-diffusion processes with discrete observations
LAN property for ergodic jump-diffusion processes with discrete observations Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Arturo Kohatsu-Higa (Ritsumeikan University, Japan) &
More informationIntensive Course on Population Genetics. Ville Mustonen
Intensive Course on Population Genetics Ville Mustonen Population genetics: The study of the distribution of inherited variation among a group of organisms of the same species [Oxford Dictionary of Biology]
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationExample: physical systems. If the state space. Example: speech recognition. Context can be. Example: epidemics. Suppose each infected
4. Markov Chains A discrete time process {X n,n = 0,1,2,...} with discrete state space X n {0,1,2,...} is a Markov chain if it has the Markov property: P[X n+1 =j X n =i,x n 1 =i n 1,...,X 0 =i 0 ] = P[X
More informationBranching Processes II: Convergence of critical branching to Feller s CSB
Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied
More information6 Introduction to Population Genetics
Grundlagen der Bioinformatik, SoSe 14, D. Huson, May 18, 2014 67 6 Introduction to Population Genetics This chapter is based on: J. Hein, M.H. Schierup and C. Wuif, Gene genealogies, variation and evolution,
More informationComputational Systems Biology: Biology X
Bud Mishra Room 1002, 715 Broadway, Courant Institute, NYU, New York, USA Human Population Genomics Outline 1 2 Damn the Human Genomes. Small initial populations; genes too distant; pestered with transposons;
More informationMULLER S RATCHET WITH COMPENSATORY MUTATIONS
The Annals of Applied Probability 212, Vol. 22, No. 5, 218 2132 DOI: 1.1214/11-AAP836 Institute of Mathematical Statistics, 212 MULLER S RATCHET WITH COMPENSATORY MUTATIONS BY P. PFAFFELHUBER 1,P.R.STAAB
More informationStatistics & Data Sciences: First Year Prelim Exam May 2018
Statistics & Data Sciences: First Year Prelim Exam May 2018 Instructions: 1. Do not turn this page until instructed to do so. 2. Start each new question on a new sheet of paper. 3. This is a closed book
More informationThe Combinatorial Interpretation of Formulas in Coalescent Theory
The Combinatorial Interpretation of Formulas in Coalescent Theory John L. Spouge National Center for Biotechnology Information NLM, NIH, DHHS spouge@ncbi.nlm.nih.gov Bldg. A, Rm. N 0 NCBI, NLM, NIH Bethesda
More informationLogFeller et Ray Knight
LogFeller et Ray Knight Etienne Pardoux joint work with V. Le and A. Wakolbinger Etienne Pardoux (Marseille) MANEGE, 18/1/1 1 / 16 Feller s branching diffusion with logistic growth We consider the diffusion
More informationBrownian motion and the Central Limit Theorem
Brownian motion and the Central Limit Theorem Amir Bar January 4, 3 Based on Shang-Keng Ma, Statistical Mechanics, sections.,.7 and the course s notes section 6. Introduction In this tutorial we shall
More informationExcited random walks in cookie environments with Markovian cookie stacks
Excited random walks in cookie environments with Markovian cookie stacks Jonathon Peterson Department of Mathematics Purdue University Joint work with Elena Kosygina December 5, 2015 Jonathon Peterson
More informationThis is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or
Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects
More informationECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering
ECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering Lecturer: Nikolay Atanasov: natanasov@ucsd.edu Teaching Assistants: Siwei Guo: s9guo@eng.ucsd.edu Anwesan Pal:
More informationTyler Hofmeister. University of Calgary Mathematical and Computational Finance Laboratory
JUMP PROCESSES GENERALIZING STOCHASTIC INTEGRALS WITH JUMPS Tyler Hofmeister University of Calgary Mathematical and Computational Finance Laboratory Overview 1. General Method 2. Poisson Processes 3. Diffusion
More informationWright-Fisher Models, Approximations, and Minimum Increments of Evolution
Wright-Fisher Models, Approximations, and Minimum Increments of Evolution William H. Press The University of Texas at Austin January 10, 2011 1 Introduction Wright-Fisher models [1] are idealized models
More informationEndowed with an Extra Sense : Mathematics and Evolution
Endowed with an Extra Sense : Mathematics and Evolution Todd Parsons Laboratoire de Probabilités et Modèles Aléatoires - Université Pierre et Marie Curie Center for Interdisciplinary Research in Biology
More informationGenetic Variation in Finite Populations
Genetic Variation in Finite Populations The amount of genetic variation found in a population is influenced by two opposing forces: mutation and genetic drift. 1 Mutation tends to increase variation. 2
More informationSimulation methods for stochastic models in chemistry
Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation
More informationDynamics of net-baryon density correlations near the QCD critical point
Dynamics of net-baryon density correlations near the QCD critical point Marcus Bluhm and Marlene Nahrgang The work of M.B. is funded by the European Union s Horizon 22 research and innovation programme
More information6 Introduction to Population Genetics
70 Grundlagen der Bioinformatik, SoSe 11, D. Huson, May 19, 2011 6 Introduction to Population Genetics This chapter is based on: J. Hein, M.H. Schierup and C. Wuif, Gene genealogies, variation and evolution,
More informationOn critical branching processes with immigration in varying environment
On critical branching processes with immigration in varying environment Márton Ispány Faculty of Informatics, University of Debrecen Alfréd Rényi Institute of Mathematics, Hungarian Academy of Sciences
More informationINFORMATION-THEORETIC BOUNDS OF EVOLUTIONARY PROCESSES MODELED AS A PROTEIN COMMUNICATION SYSTEM. Liuling Gong, Nidhal Bouaynaya and Dan Schonfeld
INFORMATION-THEORETIC BOUNDS OF EVOLUTIONARY PROCESSES MODELED AS A PROTEIN COMMUNICATION SYSTEM Liuling Gong, Nidhal Bouaynaya and Dan Schonfeld University of Illinois at Chicago, Dept. of Electrical
More informationP 1.5 X 4.5 / X 2 and (iii) The smallest value of n for
DHANALAKSHMI COLLEGE OF ENEINEERING, CHENNAI DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING MA645 PROBABILITY AND RANDOM PROCESS UNIT I : RANDOM VARIABLES PART B (6 MARKS). A random variable X
More informationElementary Applications of Probability Theory
Elementary Applications of Probability Theory With an introduction to stochastic differential equations Second edition Henry C. Tuckwell Senior Research Fellow Stochastic Analysis Group of the Centre for
More informationStochastic Demography, Coalescents, and Effective Population Size
Demography Stochastic Demography, Coalescents, and Effective Population Size Steve Krone University of Idaho Department of Mathematics & IBEST Demographic effects bottlenecks, expansion, fluctuating population
More informationYaglom-type limit theorems for branching Brownian motion with absorption. by Jason Schweinsberg University of California San Diego
Yaglom-type limit theorems for branching Brownian motion with absorption by Jason Schweinsberg University of California San Diego (with Julien Berestycki, Nathanaël Berestycki, Pascal Maillard) Outline
More informationFinite Markov Information-Exchange processes
Finite Markov Information-Exchange processes David Aldous March 30, 2011 Note that previous FMIE models were non-equilibrium. We digress to a quite different FMIE model designed to have an equilibrium.
More informationWald Lecture 2 My Work in Genetics with Jason Schweinsbreg
Wald Lecture 2 My Work in Genetics with Jason Schweinsbreg Rick Durrett Rick Durrett (Cornell) Genetics with Jason 1 / 42 The Problem Given a population of size N, how long does it take until τ k the first
More informationQuantitative trait evolution with mutations of large effect
Quantitative trait evolution with mutations of large effect May 1, 2014 Quantitative traits Traits that vary continuously in populations - Mass - Height - Bristle number (approx) Adaption - Low oxygen
More informationThe genealogy of branching Brownian motion with absorption. by Jason Schweinsberg University of California at San Diego
The genealogy of branching Brownian motion with absorption by Jason Schweinsberg University of California at San Diego (joint work with Julien Berestycki and Nathanaël Berestycki) Outline 1. A population
More informationTHE vast majority of mutations are neutral or deleterious.
Copyright Ó 2007 by the Genetics Society of America DOI: 10.1534/genetics.106.067678 Beneficial Mutation Selection Balance and the Effect of Linkage on Positive Selection Michael M. Desai*,,1 and Daniel
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationProcesses of Evolution
15 Processes of Evolution Forces of Evolution Concept 15.4 Selection Can Be Stabilizing, Directional, or Disruptive Natural selection can act on quantitative traits in three ways: Stabilizing selection
More informationBayesian Inference and the Symbolic Dynamics of Deterministic Chaos. Christopher C. Strelioff 1,2 Dr. James P. Crutchfield 2
How Random Bayesian Inference and the Symbolic Dynamics of Deterministic Chaos Christopher C. Strelioff 1,2 Dr. James P. Crutchfield 2 1 Center for Complex Systems Research and Department of Physics University
More informationMean field simulation for Monte Carlo integration. Part II : Feynman-Kac models. P. Del Moral
Mean field simulation for Monte Carlo integration Part II : Feynman-Kac models P. Del Moral INRIA Bordeaux & Inst. Maths. Bordeaux & CMAP Polytechnique Lectures, INLN CNRS & Nice Sophia Antipolis Univ.
More informationMutation-Selection Systems
Chapter 12 Mutation-Selection Systems The basic mechanisms of population biology are mutation, selection, recombination and genetic drift. In the previous chapter we concentrated on mutation and genetic
More informationExample 4.1 Let X be a random variable and f(t) a given function of time. Then. Y (t) = f(t)x. Y (t) = X sin(ωt + δ)
Chapter 4 Stochastic Processes 4. Definition In the previous chapter we studied random variables as functions on a sample space X(ω), ω Ω, without regard to how these might depend on parameters. We now
More informationSelection and Population Genetics
Selection and Population Genetics Evolution by natural selection can occur when three conditions are satisfied: Variation within populations - individuals have different traits (phenotypes). height and
More informationUncertainty quantification and systemic risk
Uncertainty quantification and systemic risk Josselin Garnier (Université Paris Diderot) with George Papanicolaou and Tzu-Wei Yang (Stanford University) February 3, 2016 Modeling systemic risk We consider
More informationGillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation
More informationThe fate of beneficial mutations in a range expansion
The fate of beneficial mutations in a range expansion R. Lehe (Paris) O. Hallatschek (Göttingen) L. Peliti (Naples) February 19, 2015 / Rutgers Outline Introduction The model Results Analysis Range expansion
More informationDynamics of Coexistence of Asexual and Sexual Reproduction in Adaptive Landscape
Dynamics of Coexistence of Asexual and Sexual Reproduction in Adaptive Landscape Shuyun Jiao,Yanbo Wang, Ping Ao Shanghai Center for Systems Biomedicine, Key Laboratory of Systems Biomedicine of Ministry
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationMacroscopic quasi-stationary distribution and microscopic particle systems
Macroscopic quasi-stationary distribution and microscopic particle systems Matthieu Jonckheere, UBA-CONICET, BCAM visiting fellow Coauthors: A. Asselah, P. Ferrari, P. Groisman, J. Martinez, S. Saglietti.
More informationGENETICS. On the Fixation Process of a Beneficial Mutation in a Variable Environment
GENETICS Supporting Information http://www.genetics.org/content/suppl/2/6/6/genetics..24297.dc On the Fixation Process of a Beneficial Mutation in a Variable Environment Hildegard Uecker and Joachim Hermisson
More informationA. Bovier () Branching Brownian motion: extremal process and ergodic theorems
Branching Brownian motion: extremal process and ergodic theorems Anton Bovier with Louis-Pierre Arguin and Nicola Kistler RCS&SM, Venezia, 06.05.2013 Plan 1 BBM 2 Maximum of BBM 3 The Lalley-Sellke conjecture
More informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 12, 215 Averaging and homogenization workshop, Luminy. Fast-slow systems
More informationThe Cameron-Martin-Girsanov (CMG) Theorem
The Cameron-Martin-Girsanov (CMG) Theorem There are many versions of the CMG Theorem. In some sense, there are many CMG Theorems. The first version appeared in ] in 944. Here we present a standard version,
More informationOptimization. The value x is called a maximizer of f and is written argmax X f. g(λx + (1 λ)y) < λg(x) + (1 λ)g(y) 0 < λ < 1; x, y X.
Optimization Background: Problem: given a function f(x) defined on X, find x such that f(x ) f(x) for all x X. The value x is called a maximizer of f and is written argmax X f. In general, argmax X f may
More informationState-Space Methods for Inferring Spike Trains from Calcium Imaging
State-Space Methods for Inferring Spike Trains from Calcium Imaging Joshua Vogelstein Johns Hopkins April 23, 2009 Joshua Vogelstein (Johns Hopkins) State-Space Calcium Imaging April 23, 2009 1 / 78 Outline
More informationAPM 541: Stochastic Modelling in Biology Diffusion Processes. Jay Taylor Fall Jay Taylor (ASU) APM 541 Fall / 61
APM 541: Stochastic Modelling in Biology Diffusion Processes Jay Taylor Fall 2013 Jay Taylor (ASU) APM 541 Fall 2013 1 / 61 Brownian Motion Brownian Motion and Random Walks Let (ξ n : n 0) be a collection
More informationThe implications of neutral evolution for neutral ecology. Daniel Lawson Bioinformatics and Statistics Scotland Macaulay Institute, Aberdeen
The implications of neutral evolution for neutral ecology Daniel Lawson Bioinformatics and Statistics Scotland Macaulay Institute, Aberdeen How is How is diversity Diversity maintained? maintained? Talk
More informationTwo viewpoints on measure valued processes
Two viewpoints on measure valued processes Olivier Hénard Université Paris-Est, Cermics Contents 1 The classical framework : from no particle to one particle 2 The lookdown framework : many particles.
More informationOn prediction and density estimation Peter McCullagh University of Chicago December 2004
On prediction and density estimation Peter McCullagh University of Chicago December 2004 Summary Having observed the initial segment of a random sequence, subsequent values may be predicted by calculating
More informationComputational statistics
Computational statistics Combinatorial optimization Thierry Denœux February 2017 Thierry Denœux Computational statistics February 2017 1 / 37 Combinatorial optimization Assume we seek the maximum of f
More informationLecture Notes 7 Random Processes. Markov Processes Markov Chains. Random Processes
Lecture Notes 7 Random Processes Definition IID Processes Bernoulli Process Binomial Counting Process Interarrival Time Process Markov Processes Markov Chains Classification of States Steady State Probabilities
More informationLecture 1. References
There follows an annotated bibliography. It is extensive, but still far from complete. I have separated the primary sources which form a less intimidating list. If I have the energy I ll add a section
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research education use, including for instruction at the authors institution
More informationSympatric Speciation
Sympatric Speciation Summary Speciation mechanisms: allopatry x sympatry The model of Dieckmann & Doebeli Sex Charles Darwin 1809-1882 Alfred R. Wallace 1823-1913 Seleção Natural Evolution by natural selection
More informationStochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
More informationPathwise construction of tree-valued Fleming-Viot processes
Pathwise construction of tree-valued Fleming-Viot processes Stephan Gufler November 9, 2018 arxiv:1404.3682v4 [math.pr] 27 Dec 2017 Abstract In a random complete and separable metric space that we call
More informationA Dynamic Contagion Process with Applications to Finance & Insurance
A Dynamic Contagion Process with Applications to Finance & Insurance Angelos Dassios Department of Statistics London School of Economics Angelos Dassios, Hongbiao Zhao (LSE) A Dynamic Contagion Process
More informationNonparametric Drift Estimation for Stochastic Differential Equations
Nonparametric Drift Estimation for Stochastic Differential Equations Gareth Roberts 1 Department of Statistics University of Warwick Brazilian Bayesian meeting, March 2010 Joint work with O. Papaspiliopoulos,
More informationThe number of accessible paths in the hypercube
Bernoulli 222), 2016, 653 680 DOI: 10.3150/14-BEJ641 The number of accessible paths in the hypercube JUIEN BERESTYCKI 1,*, ÉRIC BRUNET 2 and ZHAN SHI 1,** 1 aboratoire de Probabilités et Modèles Aléatoires,
More informationReading for Lecture 13 Release v10
Reading for Lecture 13 Release v10 Christopher Lee November 15, 2011 Contents 1 Evolutionary Trees i 1.1 Evolution as a Markov Process...................................... ii 1.2 Rooted vs. Unrooted Trees........................................
More informationA Bayesian method for the analysis of deterministic and stochastic time series
A Bayesian method for the analysis of deterministic and stochastic time series Coryn Bailer-Jones Max Planck Institute for Astronomy, Heidelberg DPG, Berlin, March 2015 Time series modelling heteroscedastic,
More informationMarkov Chain Monte Carlo (MCMC)
Markov Chain Monte Carlo (MCMC Dependent Sampling Suppose we wish to sample from a density π, and we can evaluate π as a function but have no means to directly generate a sample. Rejection sampling can
More informationAffine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009
Affine Processes Econometric specifications Eduardo Rossi University of Pavia March 17, 2009 Eduardo Rossi (University of Pavia) Affine Processes March 17, 2009 1 / 40 Outline 1 Affine Processes 2 Affine
More informationMathematical Population Genetics II
Mathematical Population Genetics II Lecture Notes Joachim Hermisson March 20, 2015 University of Vienna Mathematics Department Oskar-Morgenstern-Platz 1 1090 Vienna, Austria Copyright (c) 2013/14/15 Joachim
More informationTo appear in the American Naturalist
Unifying within- and between-generation bet-hedging theories: An ode to J.H. Gillespie Sebastian J. Schreiber Department of Evolution and Ecology, One Shields Avenue, University of California, Davis, California
More information