Page 27 - Applied Probability
P. 27

1. Basic Principles of Population Genetics
                              10
                                At a fixed point p ∞ ∈ [0, 1], we have ∆p ∞ = 0. In view of equation
                                                                                   s
                                                                                     . The third
                              (1.4), this can occur only when p ∞ equals 0, 1, or possibly
                                                                                  r+s
                              point is a legitimate fixed point if and only if r and s have the same sign.
                              In the case r> 0 and s ≤ 0, the linear function g(p)= s − (r + s)p satisfies
                              g(0) ≤ 0 and g(1) < 0. It is therefore negative throughout the open interval
                              (0, 1), and equation (1.4) implies that ∆p n < 0 for all p n ∈ (0, 1). It follows
                              that the decreasing sequence p n has a limit p ∞ < 1 when p 0 < 1. Equation
                              (1.4) shows that p ∞ > 0 is inconsistent with ∆p ∞ = lim n→∞ ∆p n =0.
                              Hence, we arrive at the intuitively obvious conclusion that the A 1 allele is
                              driven to extinction. In the opposite case r ≤ 0 and s> 0, the A 2 allele is
                              driven to extinction.
                                When r and s have the same sign, it is helpful to consider the difference
                                             s                    s
                                    p n+1 −      =∆p n + p n −
                                           r + s                r + s
                                                                       s

                                                       (r + s)p n q n p n −        s
                                                                      r+s
                                                 = −            2    2     + p n −
                                                           1 − rp − sq            r + s
                                                                n    n
                                                           2    2
                                                     1 − rp − sq − (r + s)p n q n    s
                                                           n
                                                                n
                                                 =                             p n −
                                                                 2
                                                           1 − rp − sq 2            r + s
                                                                 n    n
                                                                         s
                                                     1 − rp n − sq n
                                                 =                p n −      .
                                                           2
                                                     1 − rp − sq 2 n   r + s
                                                           n
                              If both r and s are negative, then the factor
                                                              1 − rp n − sq n
                                                   λ(p n )=
                                                                   2
                                                              1 − rp − sq 2 n
                                                                   n
                                                          > 1,
                                         s
                              and p n −    has constant sign and grows in magnitude. Therefore, ar-
                                        r+s
                              guments similar to those given in the r> 0 and s ≤ 0 case imply that
                                                     s
                                                                                   s
                              lim n→∞ p n = 0 for p 0 <  r+s  and lim n→∞ p n = 1 for p 0 >  r+s . The point
                                s  is an unstable equilibrium.
                               r+s
                                                                                         s
                                If both r and s are positive, then 0 ≤ λ(p n ) < 1, and p n −  has
                                                                                        r+s
                              constant sign and declines in magnitude. In this case, lim n→∞ p n =  s  ,
                                                                                           r+s
                              and the point  s  is a stable equilibrium.For p 0 ≈  s  ,
                                           r+s                                r+s
                                                                   s

                                                    λ(p n )  ≈ λ
                                                                 r + s
                                                               r + s − 2rs
                                                           =             ,
                                                               r + s − rs
                                         s        s  n      s
                              and p n −     ≈ λ(    ) (p 0 −  ). In other words, p n approaches its
                                        r+s      r+s       r+s
                              equilibrium value locally at the geometric rate λ(  s  ).
                                                                          r+s
                                The rate of convergence of p n to 0 or 1 depends on whether there
                              is selection against the heterozygous genotype A 1 /A 2 . Consider the case
   22   23   24   25   26   27   28   29   30   31   32