Inférence en génétique des populations IV.

Size: px
Start display at page:

Download "Inférence en génétique des populations IV."

Transcription

1 Inférence en génétique des populations IV. François Rousset & Raphaël Leblois M2 Biostatistiques FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

2 Modeling and analysis of covariances in allele frequencies Spatial Intra-family FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

3 Modeling and analysis of covariances in allele frequencies Modeling: important parameters in selection processes FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

4 Modeling and analysis of covariances in allele frequencies Modeling: important parameters in selection processes Artificial selection, or how to increase milk production by selecting bulls Wright s shifting balance theory General formal theory of selection with interactions between relatives FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

5 Modeling and analysis of covariances in allele frequencies Modeling: important parameters in selection processes Artificial selection, or how to increase milk production by selecting bulls Wright s shifting balance theory General formal theory of selection with interactions between relatives Statistical analysis: relatively simple (and routine) methods of inference FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

6 Modeling and analysis of covariances in allele frequencies Modeling: important parameters in selection processes Artificial selection, or how to increase milk production by selecting bulls Wright s shifting balance theory General formal theory of selection with interactions between relatives Statistical analysis: relatively simple (and routine) methods of inference Two main spatial models (Wright): island model, and isolation-by-distance Estimate Nm or the neighborhood size 4Dπσ 2 A few statistically desirable properties Analytical understanding FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

7 Modeling and analysis of covariances in allele frequencies Modeling: important parameters in selection processes Artificial selection, or how to increase milk production by selecting bulls Wright s shifting balance theory General formal theory of selection with interactions between relatives Statistical analysis: relatively simple (and routine) methods of inference Two main spatial models (Wright): island model, and isolation-by-distance Estimate Nm or the neighborhood size 4Dπσ 2 A few statistically desirable properties Analytical understanding Huge literature, often based on very naive ideas FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

8 Outline An example of the role of spatial structure in evolution Modeling of the role of spatial structure Formal modelling of covariances: diffusion vs. coalescence Genealogical interpretations of covariances and regressions An example of statistical inference based on covariances FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

9 Outline An example of the role of spatial structure in evolution Modeling of the role of spatial structure Formal modelling of covariances: diffusion vs. coalescence Genealogical interpretations of covariances and regressions An example of statistical inference based on covariances FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

10 Example: evolution of social behaviour Siderophore production as public good D Pseudomonas aeruginosa C D A Defector that does not produce the public good still benefits from those produced by the Cooperators. Public good production should be counter-selected. FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

11 Biologists test theories about the importance of dispersal Controlling bacterial dispersal Different agar concentrations: FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

12 Public goods production mean relative mutant fitness ± s.e. (v) % 0.25% 0.5% 1% agar concentration Viscosity Siderophore-defective strain wins Siderophore-producing strain wins Kümmerli et al, PRSB, 2009 FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

13 The effect of dispersal almost explained FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

14 The effect of dispersal almost explained FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

15 The effect of dispersal almost explained FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

16 The effect of dispersal almost explained FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

17 The effect of dispersal almost explained FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

18 The effect of dispersal almost explained FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

19 The effect of dispersal almost explained FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

20 The effect of dispersal almost explained FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

21 The effect of dispersal almost explained FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

22 The effect of dispersal almost explained FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

23 The effect of dispersal almost explained Same logic important for virtually any trait evolving in a population spatially subdivided in (small) subpopulations. FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

24 Outline An example of the role of spatial structure in evolution Modeling of the role of spatial structure Formal modelling of covariances: diffusion vs. coalescence Genealogical interpretations of covariances and regressions An example of statistical inference based on covariances FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

25 Formulation in terms of partial regression We have seen that (for deterministic models) the change in allele frequency can be written in terms of allelic fitnesses w and w, as p p =(w w )p (1 p ) =β w,1 Var(1 ) = Cov[w i, 1 (i))] where β w,1 is the regression coefficient of an individual s allelic fitness w i on her genotype 1 (i) FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

26 Formulation in terms of partial regression We have seen that (for deterministic models) the change in allele frequency can be written in terms of allelic fitnesses w and w, as p p =(w w )p (1 p ) =β w,1 Var(1 ) = Cov[w i, 1 (i))] where β w,1 is the regression coefficient of an individual s allelic fitness w i on her genotype 1 (i) We can apply the standard recursive relationship for regression coefficients between a response variable W and two predictors F and N: β W,F = β W,F.N + β N,F β W,N.F. here for F := 1, the indicator variable for the (focal) individual s genotype, and N := the mean value of the indicator variable for a group of social partners FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

27 Formulation in terms of partial regression We have seen that (for deterministic models) the change in allele frequency can be written in terms of allelic fitnesses w and w, as p p =(w w )p (1 p ) =β w,1 Var(1 ) = Cov[w i, 1 (i))] where β w,1 is the regression coefficient of an individual s allelic fitness w i on her genotype 1 (i) This entails the two-predictor regression expression for the change in allele frequency: p p = (β W,F.N + β N,F β W,N.F ) Var(1 ) where β N,F is the regression coefficient of the allele frequency among social partners to allele presence in the focal individual. β N,F is (one definition of) relatedness in the genetic literature. FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

28 Outline An example of the role of spatial structure in evolution Modeling of the role of spatial structure Formal modelling of covariances: diffusion vs. coalescence Genealogical interpretations of covariances and regressions An example of statistical inference based on covariances FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

29 Probability of identity Probability that two gene copies are of same allelic state = Identity in state FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

30 Probability of identity Probability that two gene copies are of same allelic state = Identity in state Cf Wright-Fisher model with two alleles a u A v FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

31 Probability of identity Probability that two gene copies are of same allelic state = Identity in state Cf Wright-Fisher model with two alleles a u A v pdf of allele frequency f (p) =Beta(α = 2Nu, β = 2Nv) FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

32 Probability of identity Probability that two gene copies are of same allelic state = Identity in state Cf Wright-Fisher model with two alleles a u A v pdf of allele frequency f (p) =Beta(α = 2Nu, β = 2Nv) 1 P(identity) = 0 f (x) dx FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

33 Probability of identity Probability that two gene copies are of same allelic state = Identity in state Cf Wright-Fisher model with two alleles a u A v pdf of allele frequency f (p) =Beta(α = 2Nu, β = 2Nv) 1 P(identity) = [x 2 + (1 x) 2 ]f (x) dx = u + 2Nu2 + v + 2Nv 2 0 (u + v)[1 + 2N(u + v)] Nu if u = v FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

34 Probability of identity Probability that two gene copies are of same allelic state = Identity in state Cf Wright-Fisher model with two alleles a u A v pdf of allele frequency f (p) =Beta(α = 2Nu, β = 2Nv) 1 P(identity) = [x 2 + (1 x) 2 ]f (x) dx = u + 2Nu2 + v + 2Nv 2 0 (u + v)[1 + 2N(u + v)] Nu if u = v Mathematical overkill! FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

35 Probability of identity: coalescent argument Infinite allele model identity by descent FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

36 Probability of identity: coalescent argument Infinite allele model identity by descent For c t the probability that coalescence of a a given pair of genes occurred t generations ago, P(identity) = Q = c t...? t=0 FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

37 Probability of identity: coalescent argument Infinite allele model identity by descent For c t the probability that coalescence of a a given pair of genes occurred t generations ago, P(identity) = Q = c t (1 µ) 2t = t=0 c t γ t for γ (1 µ) 2 t=0 FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

38 Probability of identity: coalescent argument Infinite allele model identity by descent For c t the probability that coalescence of a a given pair of genes occurred t generations ago, P(identity) = Q = c t (1 µ) 2t = t=0 c t γ t for γ (1 µ) 2 t=0 Generating function FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

39 Probability of identity: coalescent argument Infinite allele model identity by descent For c t the probability that coalescence of a a given pair of genes occurred t generations ago, P(identity) = Q = Generating function c t (1 µ) 2t = t=0 For Wright-Fisher model: Q = 1 c t γ t for γ (1 µ) 2 t=0 ( ) t 1 γ t γ = N N γ + N(1 γ) Nµ when N and µ 0 FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

40 Infinite island model; Wright s F ST A very simple, easily solved model Immigration m Frequency of some focal allele locally, p = p i for given island i, and globally, p. FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

41 Infinite island model; Wright s F ST Diffusion approximation for frequency in each island i: P i Const p 2Nm p 1 i (1 p i ) 2Nm(1 p) 1 =Beta(2Nm p, 2Nm(1 p)) FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

42 Infinite island model; Wright s F ST Diffusion approximation for frequency in each island i: P i Const p 2Nm p 1 i (1 p i ) 2Nm(1 p) 1 =Beta(2Nm p, 2Nm(1 p)) We have seen that the stationarity distribution satisfies: 0 = f (p, t) t a(p)f (p, t) = + 2 b(p)f (p, t) p 2 p 2 f (p) exp ( 2 p a(x)/b(x) dx ) b(p) FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

43 Infinite island model; Wright s F ST Diffusion approximation for frequency in each island i: P i Const p 2Nm p 1 i (1 p i ) 2Nm(1 p) 1 =Beta(2Nm p, 2Nm(1 p)) Wright-Fisher model: b(p) = p(1 p) and a(p) = N( vp + u(1 p)) = N[(u + v)(p u/(u + v))] Frequency P Const p 2Nu 1 (1 p) 2Nv 1 =Beta(α = 2Nu, β = 2Nv) FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

44 Infinite island model; Wright s F ST Diffusion approximation for frequency in each island i: P i Const p 2Nm p 1 i (1 p i ) 2Nm(1 p) 1 =Beta(2Nm p, 2Nm(1 p)) Wright-Fisher model: b(p) = p(1 p) and a(p) = N( vp + u(1 p)) = N[(u + v)(p u/(u + v))] Frequency P Const p 2Nu 1 (1 p) 2Nv 1 =Beta(α = 2Nu, β = 2Nv) Island model b(p i ) = p i (1 p i ) and a(p i ) = N(m(p i p)) Frequency P Beta[α = 2Nm p, β = 2Nm(1 p)] FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

45 Infinite island model; Wright s F ST Diffusion approximation for frequency in each island i: P i Const p 2Nm p 1 i (1 p i ) 2Nm(1 p) 1 =Beta(2Nm p, 2Nm(1 p)) F ST Var p i p(1 p) = E(p2 i ) p2 p(1 p) = Nm More generally variances and covariances of allele frequencies can be expressed in terms of probabilities of identity of pair of genes. Here locally, P( p i ; p) = F ST p i + (1 F ST ) p 2 (with mutation at a low enough rate: proven in next slides); hence globally. P( ; p) = F ST p + (1 F ST ) p 2 FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

46 Outline An example of the role of spatial structure in evolution Modeling of the role of spatial structure Formal modelling of covariances: diffusion vs. coalescence Genealogical interpretations of covariances and regressions Infinite island model F ST, separation of time scales Finite island model coalescence times, separation of time scales Isolation by distance separation of time scales Robustness of regression coefficients An example of statistical inference based on covariances FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

47 Genealogical interpretation of F ST R R R R R R R R R R R R R R R R R R R R R R R R R FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

48 Genealogical interpretation of F ST past generations -2-1 R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

49 Genealogical interpretation of F ST past generations p (or p i ) -2-1 R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

50 Genealogical interpretation of F ST past generations p (or p i ) p R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

51 Genealogical interpretation of F ST past generations F ST p (or p i ) (1 F ST ) p R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

52 Wright s F ST redefined (For any genetic marker) One can describe the correlation of allele types within some given sub-population relative to the allele types of random genes in the total population as F ST = Q w Q b 1 Q b in terms of probabilities of identity Q w and Q b within and between subpopulations. This leads to relatively simple analysis. FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

53 General facts about Wright s F -statistics ˆF ST = ˆQ w ˆQ b 1 ˆQ b F IS = Q 0 Q w 1 Q w FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

54 General facts about Wright s F -statistics Naive but straightforward estimators ˆF ST = ˆQ w ˆQ b 1 ˆQ b F IS = Q 0 Q w 1 Q w FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

55 General facts about Wright s F -statistics ˆF ST = ˆQ w ˆQ b 1 ˆQ b F IS = Q 0 Q w 1 Q w Relationship with average coalescence times Q = c t γ t for γ (1 µ) 2 = 1 2µ E(T ) + O(µ 2 ) t=0 F ST E(T b) E(T w ) E(T b ) when µ 0. FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

56 Modelling in terms of probabilities of identity Q Example: finite island model Relationship of Q w with Q w, Q b one generation earlier? where and FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

57 Modelling in terms of probabilities of identity Q Example: finite island model Relationship of Q w with Q w, Q b one generation earlier? where Q(t + 1) = γa[q(t) + c(t)] and FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

58 Modelling in terms of probabilities of identity Q Example: finite island model Relationship of Q w with Q w, Q b one generation earlier? Q(t + 1) = γa[q(t) + c(t)] where and Q(t) = ( Qw (t) Q b (t) ) and c(t) = ( 1 Qw(t) N 0 ) FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

59 Modelling in terms of probabilities of identity Q Example: finite island model Relationship of Q w with Q w, Q b one generation earlier? where and Q(t) = A = ( Q(t + 1) = γa[q(t) + c(t)] ( Qw (t) Q b (t) ) and c(t) = ( 1 Qw(t) N 0 ) a 11 = (1 m) 2 + m2 n d 1 1 a 11 1 a 11 a 11 n d 1 n d 1 is the transition matrix for the random walk of two gene lineages between the state within same deme and the complementary state. ) FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

60 The stationary solution in an informative form Q = γ(i γa) 1 Ac Eigensystem of A λ 1 = 1, e 1 = (1, 1) ( ) 2, λ 2 = 1 m n d n d 1 e2 = (n d 1, 1) ( 1 Qw ) c = N = 1 Q w 1 (e e 2 ) N n d FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

61 The stationary solution in an informative form Q = γ(i γa) 1 Ac Eigensystem of A λ 1 = 1, e 1 = (1, 1) ( ) 2, λ 2 = 1 m n d n d 1 e2 = (n d 1, 1) ( 1 Qw ) c = N = 1 Q w 1 (e e 2 ) N n d Q = 1 Q w N 1 n d ( γ 1 γ e 1 + γλ ) 2 e 2 1 γλ 2 FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

62 The stationary solution in an informative form Q = γ(i γa) 1 Ac Eigensystem of A λ 1 = 1, e 1 = (1, 1) ( ) 2, λ 2 = 1 m n d n d 1 e2 = (n d 1, 1) ( 1 Qw ) c = N = 1 Q w 1 (e e 2 ) N n d Q = 1 Q w N 1 n d ( γ 1 γ e 1 + γλ ) 2 e 2 1 γλ 2 Q w Q b 1 Q w = 1 N γλ 2 1 γλ 2 FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

63 The stationary solution in an informative form Q = γ(i γa) 1 Ac Eigensystem of A λ 1 = 1, e 1 = (1, 1) ( ) 2, λ 2 = 1 m n d n d 1 e2 = (n d 1, 1) ( 1 Qw ) c = N = 1 Q w 1 (e e 2 ) N n d Q = 1 Q w N 1 n d ( γ 1 γ e 1 + γλ ) 2 e 2 1 γλ 2 F ST Q w Q b 1 Q b = γλ 2 γλ 2 + N(1 γλ 2 ) Nm FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

64 The stationary solution in another informative form in terms of eigenvalues of another matrix G: Q = γa(q + c) = γ[gq + (A G)1] = γ[gq + (I G)1] where G is the matrix whose ijth element is the coefficient of the jth element of Q in the ith element of A[Q + c]. Denote l j and u j the eigenvalues and eigenvectors of G; and write (I G)1 as j x ju j. Then FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

65 The stationary solution in another informative form in terms of eigenvalues of another matrix G: Q = γa(q + c) = γ[gq + (A G)1] = γ[gq + (I G)1] where G is the matrix whose ijth element is the coefficient of the jth element of Q in the ith element of A[Q + c]. Denote l j and u j the eigenvalues and eigenvectors of G; and write (I G)1 as j x ju j. Then Q = γ(i γg) 1 j x j u j = j γx j u j = 1 γl j t γ t j l t 1 j x j u j i.e. for ith element of Q, c i,t = j lt 1 j x j u ji. FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

66 Genealogical interpretation of F -statistics (N = 100, m = 10 3, n d = 10) Probability of coalescence FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

67 Genealogical interpretation of F -statistics Interpretation of F ST from comparison of two distributions of coalescence times FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

68 Genealogical interpretation of F -statistics robustness: ancient events have no effect on F ST FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

69 Genealogical interpretation of F -statistics (N = 100, m = 10 3, n d = 10) Probability of coalescence FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

70 Genealogical interpretation of F -statistics (N = 100, m = 10 3, n d = 1000) Probability of coalescence FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

71 Genealogical interpretation of F -statistics (N = 100, m = 10 3, n d = ) Probability of coalescence FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

72 Genealogical interpretation of F ST past generations F ST p (1 F ST ) p R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R R FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

73 Island model with self-fertilization ( selfing ) Three eigenvalues of G = different time scales dependent on magnitude of n d N vs. N vs selfing rate. FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

74 Isolation by distance Pattern of isolation by distance Genealogical interpretation FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

75 Isolation by distance: a quantitative approximation Generating function of one-lineage dispersal ψ(z) := k m kz k ; FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

76 Isolation by distance: a quantitative approximation Generating function of one-lineage dispersal ψ(z) := k m kz k ; Generating function of two-lineage increments in distance dispersal ψ 2 (z); FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

77 Isolation by distance: a quantitative approximation Generating function of one-lineage dispersal ψ(z) := k m kz k ; Generating function of two-lineage increments in distance dispersal ψ 2 (z); In terms of the eigensystem of A for any distribution of dispersal distance: Q = 1 1 Q w Nn d j γλ j 1 γλ j e j λ k = ψ 2 ( e ı2πk/nx ) and e k = (cos(2πki/n x )) for (i = 0,..., n x 1); FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

78 Isolation by distance: a quantitative approximation Generating function of one-lineage dispersal ψ(z) := k m kz k ; Generating function of two-lineage increments in distance dispersal ψ 2 (z); In terms of the eigensystem of A for any distribution of dispersal distance: Q = 1 1 Q w Nn d j γλ j 1 γλ j e j λ k = ψ 2 ( e ı2πk/nx ) and e k = (cos(2πki/n x )) for (i = 0,..., n x 1); Two dimensional λ kl = ψ 2 ( e ı2πk/nx ) ψ 2 ( e ı2πl/ny ) ; FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

79 Isolation by distance: a quantitative approximation Q w Q ij 1 Q w = 1 Nn d 1 Nπ 2 k l π π Q w Q d 1 Q w log(d) 2Nπσ 2 + constant 0 0 γλ kl 1 γλ kl (e kl,0 e kl,ij ) γψ 2 (e ıx ) ψ 2 (e ıy ) 1 γψ 2 (e ıx ) ψ 2 (e ıy [1 cos(kx) cos(ly)] dx dy ) FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

80 Isolation by distance: a quantitative approximation Q w Q ij 1 Q w = 1 Nn d 1 Nπ 2 k l π π Q w Q d 1 Q w log(d) 2Nπσ 2 + constant 0 0 γλ kl 1 γλ kl (e kl,0 e kl,ij ) γψ 2 (e ıx ) ψ 2 (e ıy ) 1 γψ 2 (e ıx ) ψ 2 (e ıy [1 cos(kx) cos(ly)] dx dy ) 2Nπσ 2 Q w Q d 1 Q w FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

81 Outline An example of the role of spatial structure in evolution Modeling of the role of spatial structure Formal modelling of covariances: diffusion vs. coalescence Genealogical interpretations of covariances and regressions An example of statistical inference based on covariances FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

82 A simple inference under isolation by distance 2Nπσ 2 Q w Q d 1 Q w regress ˆQ w ˆQ d 1 ˆQ w against log(distance) to obtain an estimate of 1/2Nπσ 2. FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

83 A simple inference under isolation by distance 2Nπσ 2 Q w Q d 1 Q w regress ˆQ w ˆQ d 1 ˆQ w against log(distance) to obtain an estimate of 1/2Nπσ 2. FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

84 Statistical inference: an illustration FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

85 Statistical inference: an illustration FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

86 Statistical inference: an illustration FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

87 Statistical inference: an illustration 1/regression slope = estimate of neighborhood size Nb 4Dπσ 2 Here one finds Nb ˆ=753 with CI obtained by bootstrapping over loci. (mark-recapture experiments yielded Nb ˆ=555) FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

88 Moment vs. likelihood method Nb Regression PAC likelihood FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

89 Summary Various evolutionary processes can be understood in terms of covariances in allele frequencies; Formal analysis in terms of probabilities of identity and distributions of coalescence times; Wright s F ST and its variants can characterize recent events; They can provide some reasonable data analyses, although not the most precise ones. FR & RL Inférence en génétique des populations IV. M2 Biostatistiques / 33

Lecture 18 : Ewens sampling formula

Lecture 18 : Ewens sampling formula Lecture 8 : Ewens sampling formula MATH85K - Spring 00 Lecturer: Sebastien Roch References: [Dur08, Chapter.3]. Previous class In the previous lecture, we introduced Kingman s coalescent as a limit of

More information

Solutions to Even-Numbered Exercises to accompany An Introduction to Population Genetics: Theory and Applications Rasmus Nielsen Montgomery Slatkin

Solutions to Even-Numbered Exercises to accompany An Introduction to Population Genetics: Theory and Applications Rasmus Nielsen Montgomery Slatkin Solutions to Even-Numbered Exercises to accompany An Introduction to Population Genetics: Theory and Applications Rasmus Nielsen Montgomery Slatkin CHAPTER 1 1.2 The expected homozygosity, given allele

More information

How robust are the predictions of the W-F Model?

How 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 information

Population Genetics I. Bio

Population Genetics I. Bio Population Genetics I. Bio5488-2018 Don Conrad dconrad@genetics.wustl.edu Why study population genetics? Functional Inference Demographic inference: History of mankind is written in our DNA. We can learn

More information

Mathematical Population Genetics II

Mathematical 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 information

Population Genetics: a tutorial

Population 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 information

STAT 536: Genetic Statistics

STAT 536: Genetic Statistics STAT 536: Genetic Statistics Frequency Estimation Karin S. Dorman Department of Statistics Iowa State University August 28, 2006 Fundamental rules of genetics Law of Segregation a diploid parent is equally

More information

The Λ-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 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 information

Evolutionary Theory. Sinauer Associates, Inc. Publishers Sunderland, Massachusetts U.S.A.

Evolutionary Theory. Sinauer Associates, Inc. Publishers Sunderland, Massachusetts U.S.A. Evolutionary Theory Mathematical and Conceptual Foundations Sean H. Rice Sinauer Associates, Inc. Publishers Sunderland, Massachusetts U.S.A. Contents Preface ix Introduction 1 CHAPTER 1 Selection on One

More information

Computational Systems Biology: Biology X

Computational 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 information

Frequency Spectra and Inference in Population Genetics

Frequency 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 information

Mathematical Population Genetics II

Mathematical Population Genetics II Mathematical Population Genetics II Lecture Notes Joachim Hermisson June 9, 2018 University of Vienna Mathematics Department Oskar-Morgenstern-Platz 1 1090 Vienna, Austria Copyright (c) 2013/14/15/18 Joachim

More information

Population Structure

Population Structure Ch 4: Population Subdivision Population Structure v most natural populations exist across a landscape (or seascape) that is more or less divided into areas of suitable habitat v to the extent that populations

More information

Gene Genealogies Coalescence Theory. Annabelle Haudry Glasgow, July 2009

Gene Genealogies Coalescence Theory. Annabelle Haudry Glasgow, July 2009 Gene Genealogies Coalescence Theory Annabelle Haudry Glasgow, July 2009 What could tell a gene genealogy? How much diversity in the population? Has the demographic size of the population changed? How?

More information

Genetic Variation in Finite Populations

Genetic 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 information

Methods for Cryptic Structure. Methods for Cryptic Structure

Methods for Cryptic Structure. Methods for Cryptic Structure Case-Control Association Testing Review Consider testing for association between a disease and a genetic marker Idea is to look for an association by comparing allele/genotype frequencies between the cases

More information

Selection and Population Genetics

Selection 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 information

THE CENTRAL CONCEPTS OF INCLUSIVE FITNESS 3000 word article in the Oxford University Press Encyclopedia of Evolution, January 2002.

THE CENTRAL CONCEPTS OF INCLUSIVE FITNESS 3000 word article in the Oxford University Press Encyclopedia of Evolution, January 2002. 1 TH CNTRAL CONCPTS OF INCLUSIV FITNSS 3 word article in the Oxford University Press ncyclopedia of volution, January. Peter D. Taylor Dept. of Mathematics and Statistics Queen's University Kingston ON

More information

Evolutionary dynamics of populations with genotype-phenotype map

Evolutionary dynamics of populations with genotype-phenotype map Evolutionary dynamics of populations with genotype-phenotype map Esther Ibáñez Marcelo, Tomas Alarcon Cor Biomat 2013: Mathematics of Planet Earth CENTRE DE RECERCA MATEMÀTICA 17-21 June 2013 Esther Ibáñez

More information

Quantitative trait evolution with mutations of large effect

Quantitative 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 information

6 Introduction to Population Genetics

6 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 information

Mathematical models in population genetics II

Mathematical 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 information

Supplemental Information Likelihood-based inference in isolation-by-distance models using the spatial distribution of low-frequency alleles

Supplemental Information Likelihood-based inference in isolation-by-distance models using the spatial distribution of low-frequency alleles Supplemental Information Likelihood-based inference in isolation-by-distance models using the spatial distribution of low-frequency alleles John Novembre and Montgomery Slatkin Supplementary Methods To

More information

The Wright-Fisher Model and Genetic Drift

The 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 information

Demography April 10, 2015

Demography April 10, 2015 Demography April 0, 205 Effective Population Size The Wright-Fisher model makes a number of strong assumptions which are clearly violated in many populations. For example, it is unlikely that any population

More information

Stochastic Demography, Coalescents, and Effective Population Size

Stochastic 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 information

Statistics & Data Sciences: First Year Prelim Exam May 2018

Statistics & 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 information

I of a gene sampled from a randomly mating popdation,

I of a gene sampled from a randomly mating popdation, Copyright 0 1987 by the Genetics Society of America Average Number of Nucleotide Differences in a From a Single Subpopulation: A Test for Population Subdivision Curtis Strobeck Department of Zoology, University

More information

Diffusion Models in Population Genetics

Diffusion 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 information

Lecture 9. QTL Mapping 2: Outbred Populations

Lecture 9. QTL Mapping 2: Outbred Populations Lecture 9 QTL Mapping 2: Outbred Populations Bruce Walsh. Aug 2004. Royal Veterinary and Agricultural University, Denmark The major difference between QTL analysis using inbred-line crosses vs. outbred

More information

6 Introduction to Population Genetics

6 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 information

Cooperation Achieved by Migration and Evolution in a Multilevel Selection Context

Cooperation Achieved by Migration and Evolution in a Multilevel Selection Context Proceedings of the 27 IEEE Symposium on Artificial Life (CI-ALife 27) Cooperation Achieved by Migration and Evolution in a Multilevel Selection Context Genki Ichinose Graduate School of Information Science

More information

Stationary Distribution of the Linkage Disequilibrium Coefficient r 2

Stationary Distribution of the Linkage Disequilibrium Coefficient r 2 Stationary Distribution of the Linkage Disequilibrium Coefficient r 2 Wei Zhang, Jing Liu, Rachel Fewster and Jesse Goodman Department of Statistics, The University of Auckland December 1, 2015 Overview

More information

Evolution of cooperation in finite populations. Sabin Lessard Université de Montréal

Evolution of cooperation in finite populations. Sabin Lessard Université de Montréal Evolution of cooperation in finite populations Sabin Lessard Université de Montréal Title, contents and acknowledgements 2011 AMS Short Course on Evolutionary Game Dynamics Contents 1. Examples of cooperation

More information

Evolution of quantitative traits

Evolution of quantitative traits Evolution of quantitative traits Introduction Let s stop and review quickly where we ve come and where we re going We started our survey of quantitative genetics by pointing out that our objective was

More information

Mathematical Population Genetics. Introduction to the Stochastic Theory. Lecture Notes. Guanajuato, March Warren J Ewens

Mathematical Population Genetics. Introduction to the Stochastic Theory. Lecture Notes. Guanajuato, March Warren J Ewens 1 Mathematical Population Genetics Introduction to the Stochastic Theory Lecture Notes Guanajuato, March 2009 Warren J Ewens 2 Preface These notes should, ideally, be read before the Guanajuato meeting

More information

Lecture 32: Infinite-dimensional/Functionvalued. Functions and Random Regressions. Bruce Walsh lecture notes Synbreed course version 11 July 2013

Lecture 32: Infinite-dimensional/Functionvalued. Functions and Random Regressions. Bruce Walsh lecture notes Synbreed course version 11 July 2013 Lecture 32: Infinite-dimensional/Functionvalued Traits: Covariance Functions and Random Regressions Bruce Walsh lecture notes Synbreed course version 11 July 2013 1 Longitudinal traits Many classic quantitative

More information

BIOL Evolution. Lecture 9

BIOL Evolution. Lecture 9 BIOL 432 - Evolution Lecture 9 J Krause et al. Nature 000, 1-4 (2010) doi:10.1038/nature08976 Selection http://www.youtube.com/watch?v=a38k mj0amhc&feature=playlist&p=61e033 F110013706&index=0&playnext=1

More information

Contrasts for a within-species comparative method

Contrasts for a within-species comparative method Contrasts for a within-species comparative method Joseph Felsenstein, Department of Genetics, University of Washington, Box 357360, Seattle, Washington 98195-7360, USA email address: joe@genetics.washington.edu

More information

Breeding Values and Inbreeding. Breeding Values and Inbreeding

Breeding Values and Inbreeding. Breeding Values and Inbreeding Breeding Values and Inbreeding Genotypic Values For the bi-allelic single locus case, we previously defined the mean genotypic (or equivalently the mean phenotypic values) to be a if genotype is A 2 A

More information

19. Genetic Drift. The biological context. There are four basic consequences of genetic drift:

19. Genetic Drift. The biological context. There are four basic consequences of genetic drift: 9. Genetic Drift Genetic drift is the alteration of gene frequencies due to sampling variation from one generation to the next. It operates to some degree in all finite populations, but can be significant

More information

Association studies and regression

Association studies and regression Association studies and regression CM226: Machine Learning for Bioinformatics. Fall 2016 Sriram Sankararaman Acknowledgments: Fei Sha, Ameet Talwalkar Association studies and regression 1 / 104 Administration

More information

Problems on Evolutionary dynamics

Problems 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 information

Overview. Background

Overview. Background Overview Implementation of robust methods for locating quantitative trait loci in R Introduction to QTL mapping Andreas Baierl and Andreas Futschik Institute of Statistics and Decision Support Systems

More information

Population genetics snippets for genepop

Population genetics snippets for genepop Population genetics snippets for genepop Peter Beerli August 0, 205 Contents 0.Basics 0.2Exact test 2 0.Fixation indices 4 0.4Isolation by Distance 5 0.5Further Reading 8 0.6References 8 0.7Disclaimer

More information

This 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 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 information

Challenges when applying stochastic models to reconstruct the demographic history of populations.

Challenges when applying stochastic models to reconstruct the demographic history of populations. Challenges when applying stochastic models to reconstruct the demographic history of populations. Willy Rodríguez Institut de Mathématiques de Toulouse October 11, 2017 Outline 1 Introduction 2 Inverse

More information

arxiv: v1 [q-bio.pe] 14 Jan 2016

arxiv: v1 [q-bio.pe] 14 Jan 2016 Efficient Maximum-Likelihood Inference For The Isolation-With-Initial-Migration Model With Potentially Asymmetric Gene Flow Rui J. Costa and Hilde Wilkinson-Herbots arxiv:1601.03684v1 [q-bio.pe] 14 Jan

More information

Supplementary Figures.

Supplementary Figures. Supplementary Figures. Supplementary Figure 1 The extended compartment model. Sub-compartment C (blue) and 1-C (yellow) represent the fractions of allele carriers and non-carriers in the focal patch, respectively,

More information

Microevolution 2 mutation & migration

Microevolution 2 mutation & migration Microevolution 2 mutation & migration Assumptions of Hardy-Weinberg equilibrium 1. Mating is random 2. Population size is infinite (i.e., no genetic drift) 3. No migration 4. No mutation 5. No selection

More information

Introductory seminar on mathematical population genetics

Introductory seminar on mathematical population genetics Exercises Sheets Introductory seminar on mathematical population genetics WS 20/202 Kristan Schneider, Ada Akerman Ex Assume a single locus with alleles A and A 2 Denote the frequencies of the three (unordered

More information

Unit roots in vector time series. Scalar autoregression True model: y t 1 y t1 2 y t2 p y tp t Estimated model: y t c y t1 1 y t1 2 y t2

Unit roots in vector time series. Scalar autoregression True model: y t 1 y t1 2 y t2 p y tp t Estimated model: y t c y t1 1 y t1 2 y t2 Unit roots in vector time series A. Vector autoregressions with unit roots Scalar autoregression True model: y t y t y t p y tp t Estimated model: y t c y t y t y t p y tp t Results: T j j is asymptotically

More information

Lecture 2: Genetic Association Testing with Quantitative Traits. Summer Institute in Statistical Genetics 2017

Lecture 2: Genetic Association Testing with Quantitative Traits. Summer Institute in Statistical Genetics 2017 Lecture 2: Genetic Association Testing with Quantitative Traits Instructors: Timothy Thornton and Michael Wu Summer Institute in Statistical Genetics 2017 1 / 29 Introduction to Quantitative Trait Mapping

More information

Index. Causality concept of, 128 selection and, 139, 298. Adaptation, 6, 7. See also Biotic adaptation. defining, 55, 133, 301

Index. Causality concept of, 128 selection and, 139, 298. Adaptation, 6, 7. See also Biotic adaptation. defining, 55, 133, 301 Index Adaptation, 6, 7. See also Biotic adaptation altruistic, 60-61 defining, 55, 133, 301 group, 136-137 group -related, 52, 53 at level of multicellular vehicle, 172 meaning of, 238 at organismallevel,

More information

COPYRIGHTED MATERIAL CONTENTS. Preface Preface to the First Edition

COPYRIGHTED MATERIAL CONTENTS. Preface Preface to the First Edition Preface Preface to the First Edition xi xiii 1 Basic Probability Theory 1 1.1 Introduction 1 1.2 Sample Spaces and Events 3 1.3 The Axioms of Probability 7 1.4 Finite Sample Spaces and Combinatorics 15

More information

Chapter 8: Introduction to Evolutionary Computation

Chapter 8: Introduction to Evolutionary Computation Computational Intelligence: Second Edition Contents Some Theories about Evolution Evolution is an optimization process: the aim is to improve the ability of an organism to survive in dynamically changing

More information

Computational statistics

Computational 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 information

Part IB Statistics. Theorems with proof. Based on lectures by D. Spiegelhalter Notes taken by Dexter Chua. Lent 2015

Part IB Statistics. Theorems with proof. Based on lectures by D. Spiegelhalter Notes taken by Dexter Chua. Lent 2015 Part IB Statistics Theorems with proof Based on lectures by D. Spiegelhalter Notes taken by Dexter Chua Lent 2015 These notes are not endorsed by the lecturers, and I have modified them (often significantly)

More information

DNA polymorphisms such as SNP and familial effects (additive genetic, common environment) to

DNA polymorphisms such as SNP and familial effects (additive genetic, common environment) to 1 1 1 1 1 1 1 1 0 SUPPLEMENTARY MATERIALS, B. BIVARIATE PEDIGREE-BASED ASSOCIATION ANALYSIS Introduction We propose here a statistical method of bivariate genetic analysis, designed to evaluate contribution

More information

Time Series: Theory and Methods

Time Series: Theory and Methods Peter J. Brockwell Richard A. Davis Time Series: Theory and Methods Second Edition With 124 Illustrations Springer Contents Preface to the Second Edition Preface to the First Edition vn ix CHAPTER 1 Stationary

More information

Endowed with an Extra Sense : Mathematics and Evolution

Endowed 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 information

URN MODELS: the Ewens Sampling Lemma

URN 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 information

Statistical population genetics

Statistical population genetics Statistical population genetics Lecture 2: Wright-Fisher model Xavier Didelot Dept of Statistics, Univ of Oxford didelot@stats.ox.ac.uk Slide 21 of 161 Heterozygosity One measure of the diversity of a

More information

Mathematical Biology. Sabin Lessard

Mathematical Biology. Sabin Lessard J. Math. Biol. DOI 10.1007/s00285-008-0248-1 Mathematical Biology Diffusion approximations for one-locus multi-allele kin selection, mutation and random drift in group-structured populations: a unifying

More information

Darwinian Selection. Chapter 7 Selection I 12/5/14. v evolution vs. natural selection? v evolution. v natural selection

Darwinian Selection. Chapter 7 Selection I 12/5/14. v evolution vs. natural selection? v evolution. v natural selection Chapter 7 Selection I Selection in Haploids Selection in Diploids Mutation-Selection Balance Darwinian Selection v evolution vs. natural selection? v evolution ² descent with modification ² change in allele

More information

Quantitative Genetics I: Traits controlled my many loci. Quantitative Genetics: Traits controlled my many loci

Quantitative Genetics I: Traits controlled my many loci. Quantitative Genetics: Traits controlled my many loci Quantitative Genetics: Traits controlled my many loci So far in our discussions, we have focused on understanding how selection works on a small number of loci (1 or 2). However in many cases, evolutionary

More information

Inbreeding depression due to stabilizing selection on a quantitative character. Emmanuelle Porcher & Russell Lande

Inbreeding depression due to stabilizing selection on a quantitative character. Emmanuelle Porcher & Russell Lande Inbreeding depression due to stabilizing selection on a quantitative character Emmanuelle Porcher & Russell Lande Inbreeding depression Reduction in fitness of inbred vs. outbred individuals Outcrossed

More information

Surfing genes. On the fate of neutral mutations in a spreading population

Surfing genes. On the fate of neutral mutations in a spreading population Surfing genes On the fate of neutral mutations in a spreading population Oskar Hallatschek David Nelson Harvard University ohallats@physics.harvard.edu Genetic impact of range expansions Population expansions

More information

Case-Control Association Testing. Case-Control Association Testing

Case-Control Association Testing. Case-Control Association Testing Introduction Association mapping is now routinely being used to identify loci that are involved with complex traits. Technological advances have made it feasible to perform case-control association studies

More information

Hierarchical Linear Models. Hierarchical Linear Models. Much of this material already seen in Chapters 5 and 14. Hyperprior on K parameters α:

Hierarchical Linear Models. Hierarchical Linear Models. Much of this material already seen in Chapters 5 and 14. Hyperprior on K parameters α: Hierarchical Linear Models Hierarchical Linear Models Much of this material already seen in Chapters 5 and 14 Hierarchical linear models combine regression framework with hierarchical framework Unified

More information

Lecture 4: Allelic Effects and Genetic Variances. Bruce Walsh lecture notes Tucson Winter Institute 7-9 Jan 2013

Lecture 4: Allelic Effects and Genetic Variances. Bruce Walsh lecture notes Tucson Winter Institute 7-9 Jan 2013 Lecture 4: Allelic Effects and Genetic Variances Bruce Walsh lecture notes Tucson Winter Institute 7-9 Jan 2013 1 Basic model of Quantitative Genetics Phenotypic value -- we will occasionally also use

More information

Proportional Variance Explained by QLT and Statistical Power. Proportional Variance Explained by QTL and Statistical Power

Proportional Variance Explained by QLT and Statistical Power. Proportional Variance Explained by QTL and Statistical Power Proportional Variance Explained by QTL and Statistical Power Partitioning the Genetic Variance We previously focused on obtaining variance components of a quantitative trait to determine the proportion

More information

A 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 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 information

(Genome-wide) association analysis

(Genome-wide) association analysis (Genome-wide) association analysis 1 Key concepts Mapping QTL by association relies on linkage disequilibrium in the population; LD can be caused by close linkage between a QTL and marker (= good) or by

More information

The co-evolution of culturally inherited altruistic helping and cultural transmission under random group formation

The co-evolution of culturally inherited altruistic helping and cultural transmission under random group formation Theoretical Population Biology 73 (2008) 506 516 Contents lists available at ScienceDirect Theoretical Population Biology journal homepage: www.elsevier.com/locate/tpb The co-evolution of culturally inherited

More information

Covariance function estimation in Gaussian process regression

Covariance function estimation in Gaussian process regression Covariance function estimation in Gaussian process regression François Bachoc Department of Statistics and Operations Research, University of Vienna WU Research Seminar - May 2015 François Bachoc Gaussian

More information

Closed-form sampling formulas for the coalescent with recombination

Closed-form sampling formulas for the coalescent with recombination 0 / 21 Closed-form sampling formulas for the coalescent with recombination Yun S. Song CS Division and Department of Statistics University of California, Berkeley September 7, 2009 Joint work with Paul

More information

Fitness landscapes and seascapes

Fitness landscapes and seascapes Fitness landscapes and seascapes Michael Lässig Institute for Theoretical Physics University of Cologne Thanks Ville Mustonen: Cross-species analysis of bacterial promoters, Nonequilibrium evolution of

More information

MIXED MODELS THE GENERAL MIXED MODEL

MIXED MODELS THE GENERAL MIXED MODEL MIXED MODELS This chapter introduces best linear unbiased prediction (BLUP), a general method for predicting random effects, while Chapter 27 is concerned with the estimation of variances by restricted

More information

... x. Variance NORMAL DISTRIBUTIONS OF PHENOTYPES. Mice. Fruit Flies CHARACTERIZING A NORMAL DISTRIBUTION MEAN VARIANCE

... x. Variance NORMAL DISTRIBUTIONS OF PHENOTYPES. Mice. Fruit Flies CHARACTERIZING A NORMAL DISTRIBUTION MEAN VARIANCE NORMAL DISTRIBUTIONS OF PHENOTYPES Mice Fruit Flies In:Introduction to Quantitative Genetics Falconer & Mackay 1996 CHARACTERIZING A NORMAL DISTRIBUTION MEAN VARIANCE Mean and variance are two quantities

More information

The Admixture Model in Linkage Analysis

The Admixture Model in Linkage Analysis The Admixture Model in Linkage Analysis Jie Peng D. Siegmund Department of Statistics, Stanford University, Stanford, CA 94305 SUMMARY We study an appropriate version of the score statistic to test the

More information

Recovery of a recessive allele in a Mendelian diploid model

Recovery of a recessive allele in a Mendelian diploid model Recovery of a recessive allele in a Mendelian diploid model Loren Coquille joint work with A. Bovier and R. Neukirch (University of Bonn) Young Women in Probability and Analysis 2016, Bonn Outline 1 Introduction

More information

STAT 730 Chapter 4: Estimation

STAT 730 Chapter 4: Estimation STAT 730 Chapter 4: Estimation Timothy Hanson Department of Statistics, University of South Carolina Stat 730: Multivariate Analysis 1 / 23 The likelihood We have iid data, at least initially. Each datum

More information

Inference For High Dimensional M-estimates. Fixed Design Results

Inference For High Dimensional M-estimates. Fixed Design Results : Fixed Design Results Lihua Lei Advisors: Peter J. Bickel, Michael I. Jordan joint work with Peter J. Bickel and Noureddine El Karoui Dec. 8, 2016 1/57 Table of Contents 1 Background 2 Main Results and

More information

EVOLUTIONARY DYNAMICS AND THE EVOLUTION OF MULTIPLAYER COOPERATION IN A SUBDIVIDED POPULATION

EVOLUTIONARY DYNAMICS AND THE EVOLUTION OF MULTIPLAYER COOPERATION IN A SUBDIVIDED POPULATION Friday, July 27th, 11:00 EVOLUTIONARY DYNAMICS AND THE EVOLUTION OF MULTIPLAYER COOPERATION IN A SUBDIVIDED POPULATION Karan Pattni karanp@liverpool.ac.uk University of Liverpool Joint work with Prof.

More information

Measurement Error and Linear Regression of Astronomical Data. Brandon Kelly Penn State Summer School in Astrostatistics, June 2007

Measurement Error and Linear Regression of Astronomical Data. Brandon Kelly Penn State Summer School in Astrostatistics, June 2007 Measurement Error and Linear Regression of Astronomical Data Brandon Kelly Penn State Summer School in Astrostatistics, June 2007 Classical Regression Model Collect n data points, denote i th pair as (η

More information

1.3 Convergence of Regular Markov Chains

1.3 Convergence of Regular Markov Chains Markov Chains and Random Walks on Graphs 3 Applying the same argument to A T, which has the same λ 0 as A, yields the row sum bounds Corollary 0 Let P 0 be the transition matrix of a regular Markov chain

More information

Brief history of The Prisoner s Dilemma (From Harman s The Price of Altruism)

Brief history of The Prisoner s Dilemma (From Harman s The Price of Altruism) Brief history of The Prisoner s Dilemma (From Harman s The Price of Altruism) 1948 The RAND Corporation. a civilian nonprofit think tank of an elite cadre of physicists, mathematicians, economists, and

More information

Coalescent based demographic inference. Daniel Wegmann University of Fribourg

Coalescent based demographic inference. Daniel Wegmann University of Fribourg Coalescent based demographic inference Daniel Wegmann University of Fribourg Introduction The current genetic diversity is the outcome of past evolutionary processes. Hence, we can use genetic diversity

More information

Selection on multiple characters

Selection on multiple characters Selection on multiple characters Introduction So far we ve studied only the evolution of a single trait, e.g., height or weight. But organisms have many traits, and they evolve at the same time. How can

More information

Lecture 2: Introduction to Quantitative Genetics

Lecture 2: Introduction to Quantitative Genetics Lecture 2: Introduction to Quantitative Genetics Bruce Walsh lecture notes Introduction to Quantitative Genetics SISG, Seattle 16 18 July 2018 1 Basic model of Quantitative Genetics Phenotypic value --

More information

DNA-based species delimitation

DNA-based species delimitation DNA-based species delimitation Phylogenetic species concept based on tree topologies Ø How to set species boundaries? Ø Automatic species delimitation? druhů? DNA barcoding Species boundaries recognized

More information

From Individual-based Population Models to Lineage-based Models of Phylogenies

From Individual-based Population Models to Lineage-based Models of Phylogenies From Individual-based Population Models to Lineage-based Models of Phylogenies Amaury Lambert (joint works with G. Achaz, H.K. Alexander, R.S. Etienne, N. Lartillot, H. Morlon, T.L. Parsons, T. Stadler)

More information

Chapter 3 - Temporal processes

Chapter 3 - Temporal processes STK4150 - Intro 1 Chapter 3 - Temporal processes Odd Kolbjørnsen and Geir Storvik January 23 2017 STK4150 - Intro 2 Temporal processes Data collected over time Past, present, future, change Temporal aspect

More information

Effective population size and patterns of molecular evolution and variation

Effective population size and patterns of molecular evolution and variation FunDamental concepts in genetics Effective population size and patterns of molecular evolution and variation Brian Charlesworth Abstract The effective size of a population,, determines the rate of change

More information

OVERVIEW. L5. Quantitative population genetics

OVERVIEW. L5. Quantitative population genetics L5. Quantitative population genetics OVERVIEW. L1. Approaches to ecological modelling. L2. Model parameterization and validation. L3. Stochastic models of population dynamics (math). L4. Animal movement

More information

Partitioning the Genetic Variance

Partitioning the Genetic Variance Partitioning the Genetic Variance 1 / 18 Partitioning the Genetic Variance In lecture 2, we showed how to partition genotypic values G into their expected values based on additivity (G A ) and deviations

More information

The Generalized Higher Criticism for Testing SNP-sets in Genetic Association Studies

The Generalized Higher Criticism for Testing SNP-sets in Genetic Association Studies The Generalized Higher Criticism for Testing SNP-sets in Genetic Association Studies Ian Barnett, Rajarshi Mukherjee & Xihong Lin Harvard University ibarnett@hsph.harvard.edu June 24, 2014 Ian Barnett

More information

REGRESSION WITH SPATIALLY MISALIGNED DATA. Lisa Madsen Oregon State University David Ruppert Cornell University

REGRESSION WITH SPATIALLY MISALIGNED DATA. Lisa Madsen Oregon State University David Ruppert Cornell University REGRESSION ITH SPATIALL MISALIGNED DATA Lisa Madsen Oregon State University David Ruppert Cornell University SPATIALL MISALIGNED DATA 10 X X X X X X X X 5 X X X X X 0 X 0 5 10 OUTLINE 1. Introduction 2.

More information

STAT 536: Migration. Karin S. Dorman. October 3, Department of Statistics Iowa State University

STAT 536: Migration. Karin S. Dorman. October 3, Department of Statistics Iowa State University STAT 536: Migration Karin S. Dorman Department of Statistics Iowa State University October 3, 2006 Migration Introduction Migration is the movement of individuals between populations. Until now we have

More information