1 / 12

Two-locus selection

Two-locus selection. p 1 ’ =p 1 2 +p 1 p 2 +p 1 p 3 +(1-r)p 1 p 4 +rp 2 p 3 p 2 ’ =p 2 2 +p 1 p 2 +p 2 p 4 +rp 1 p 4 +(1-r)p 2 p 3 p 3 ’ =p 3 2 +p 3 p 4 +p 1 p 3 +rp 1 p 4 +(1-r)p 2 p 3 p 4 ’ =p 4 2 +p 3 p 4 +p 2 p 4 +(1-r)p 1 p 4 +rp 2 p 3.

quang
Télécharger la présentation

Two-locus selection

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Two-locus selection

  2. p1’=p12+p1p2+p1p3+(1-r)p1p4+rp2p3 p2’=p22+p1p2+p2p4+rp1p4+(1-r)p2p3 p3’=p32+p3p4+p1p3+rp1p4+(1-r)p2p3 p4’=p42+p3p4+p2p4+(1-r)p1p4+rp2p3 r is probabilities of cross-over (coefficient of recombination). Usually 0 r  0.5. If r=0.5 then loci are called unlinked (or independent). If r=0 then population transform to one loci population with four alleles. AB Ab aB ab p1 p2 p3 p4

  3. AB Ab aB ab p1 p2 p3 p4 p1’=p12+p1p2+p1p3+p1p4 - rD p2’=p22+p1p2+p2p4 +p1p4+ rD p3’=p32+p3p4+p1p3 +p1p4+ rD p4’=p42+p3p4+p2p4 +p1p4- rD The miniumum and the maximum for D

  4. Selection coefficients Haploid genotypes AB, Ab, aB, ab Zygote genotypes ABAB, ABAb, ABaB ABab AbAB, AbAb, AbaB Abab aBAB, aBAb, aBaB aBab abAB, abAb, abaB abab Fitness coefficients BB Bb bb AA w11 w12 w22 Aa w13 w14 w24 aa w33 w34 w44

  5. p1’=p12+p1p2+p1p3+p1p4 - rD p2’=p22+p1p2+p2p4 +p1p4+ rD p3’=p32+p3p4+p1p3 +p1p4+ rD p4’=p42+p3p4+p2p4 +p1p4- rD p1’=( w11p12+w12p1p2+w13p1p3+w14 p1p4 –w14rD)/W p2’=( w22p22+w12p1p2+w24p2p4+w33p2p3 +w14rD)/W p3’=( w33p32+w34p3p4+w13p1p3+w33p2p3 +w14rD)/W p4’=( w44p42+w34p3p4+w24p2p4+w14p1p4 –w14rD)/W Measure of disequilibria D= p1p4-p2p3

  6. One Selected Locus, One Neutral Locus

  7. Haploid genotypes AB, Ab, aB, ab p1’=( w11p12+w12p1p2+w13p1p3+w14 p1p4 –w14rD)/W p2’=( w22p22+w12p1p2+w24p2p4+w33p2p3 +w14rD)/W p3’=( w33p32+w34p3p4+w13p1p3+w33p2p3 +w14rD)/W p4’=( w44p42+w34p3p4+w24p2p4+w14p1p4 –w14rD)/W Fitness coefficients BB Bb bb AA w11 w12 w22 Aa w13 w14 w24 aa w33 w34 w44 p1’=( p12+(1+hs)p1p2+p1p3+(1+hs)p1p4 –(1+hs)rD)/W p2’=((1+s)p22+p1p2+(1+s)p2p4+p2p3 +(1+hs)rD)/W p3’=( p32+(1+hs)p3p4+p1p3+p2p3 +(1+hs)rD)/W p4’=((1+s)p42+(1+hs)p3p4+(1+s)p2p4+(1+hs)p1p4 –(1+hs)rD)/W

  8. p1’=( p12+(1+hs)p1p2+p1p3+(1+hs)p1p4 –(1+hs)rD)/W p2’=((1+s)p22+p1p2+(1+s)p2p4+p2p3 +(1+hs)rD)/W p3’=( p32+(1+hs)p3p4+p1p3+p2p3 +(1+hs)rD)/W p4’=((1+s)p42+(1+hs)p3p4+(1+s)p2p4+(1+hs)p1p4 –(1+hs)rD)/W Evolutionary operator

  9. Only when D=0 will the frequency of a remain constant over time, as expected from the one-locus model. Otherwise, the a allele will change in frequency because of its association (linkage disequilibrium) with another selected locus: Genetic Hitchhiking. When s is positive (b favored) and D is positive (AB and ab are more common than expected), the a allele will increase in frequency over time.         When s is positive (b favored) and D is negative (Ab and aB are more common than expected), the a allele will decrease in frequency over time. If recombination is loose (r > s), disequilibrium decays rapidly and the a allele will change little in frequency. If recombination is tight (r < s), the extent of hitchhiking can b

  10. Result 2 • Linkage disequilibrium does not always decay to zero but can be maintained indefinitely within a population. • In fact, very special conditions are required for linkage equilibrium (D=0) to be maintained when both loci are • polymorphic and under selection. • Result 3 • Continuous cycling can occur (Hastings 1981): "apparent cases of directional selection may be due to stable cycling". • Result 4 • Even when the double heterozygote has the highest fitness, there may be no stable equilibrium maintaining • polymorphisms at both loci! • Claim 1 • If there is pure directional selection favoring a genotype, say A2B2 (such that wi,j+1 > wi,j and wi+1,j > wi,j), it has • been claimed but not proven • that selection will lead eventually to the fixation of A2B2. • Nordborg (pers. comm.) ran simulations with millions of fitness sets, which all confirmed this claim, but it has • not been shown analytically. • More results are known for special cases of the two-locus model: • Special Case 1: r=0 • When r=0, the system is dynamically the same as a one-locus four allele model. •        Mean fitness never decreases. • Special Case 2: Additive fitnesses • If the fitnesses of the A and B loci add together to give two-locus fitnesses, then: • There is only one equilibrium with both loci polymorphc for r>0, which has x1 = pA1 pB1, x2 = pA1 pB2, x3 = pA2 pB1, • and x4 = pA2 pB2, where pA1 and pB1 equal the allele frequencies found at equilibrium in the one-locus model with • heterozygote advantage. • This "internal" equilibrium is globally stable if there is heterozygote advantage at both loci. • D=0 at equilibrium. • Mean fitness never decreases (Ewens 1969). • Away from equilibrium, selection WILL generate linkage disequilibrium within a population that is initially in linkage equilibrium. • Special Case 3: Multiplicative fitnesses • If the fitnesses of the A and B loci multiply together to give two-locus fitnesses, then: • The same internal polymorphism exists with x1 = pA1 pB1, x2 = pA1 pB2, x3 = pA2 pB1, and x4 = pA2 pB2. • This "internal" equilibrium is only stable if there is heterozygote advantage at both loci AND if r is large enough. • For tight linkage, the population goes to one of two stable equilibria at which D ≠ 0

More Related