<p>Population genetic processes, such as the adaptation of a quantitative trait to directional selection, may occur on longer time scales than the sweep of a single advantageous mutation. To study such processes in finite populations, approximations for the time course of the distribution of a beneficial mutation were derived previously by branching process methods. The application to the evolution of a quantitative trait requires bounds for the probability of survival <InlineEquation ID="IEq1"> <EquationSource Format="TEX">\(S^{(n)}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>S</mi> <mrow> <mo stretchy="false">(</mo> <mi>n</mi> <mo stretchy="false">)</mo> </mrow> </msup> </math></EquationSource> </InlineEquation> up to generation <i>n</i> of a single beneficial mutation. Here, we present a method to obtain a simple, analytically explicit, either upper or lower, bound for <InlineEquation ID="IEq2"> <EquationSource Format="TEX">\(S^{(n)}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>S</mi> <mrow> <mo stretchy="false">(</mo> <mi>n</mi> <mo stretchy="false">)</mo> </mrow> </msup> </math></EquationSource> </InlineEquation> in a supercritical Galton-Watson process. We prove the existence of an upper bound for offspring distributions including Poisson, binomial, and negative binomial. They are constructed by bounding the given generating function, <InlineEquation ID="IEq3"> <EquationSource Format="TEX">\(\varphi \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>φ</mi> </math></EquationSource> </InlineEquation>, by a fractional linear one that has the same survival probability <InlineEquation ID="IEq4"> <EquationSource Format="TEX">\(S^\infty \)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>S</mi> <mi>∞</mi> </msup> </math></EquationSource> </InlineEquation> and yields the same rate of convergence of <InlineEquation ID="IEq5"> <EquationSource Format="TEX">\(S^{(n)}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>S</mi> <mrow> <mo stretchy="false">(</mo> <mi>n</mi> <mo stretchy="false">)</mo> </mrow> </msup> </math></EquationSource> </InlineEquation> to <InlineEquation ID="IEq6"> <EquationSource Format="TEX">\(S^\infty \)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>S</mi> <mi>∞</mi> </msup> </math></EquationSource> </InlineEquation> as <InlineEquation ID="IEq7"> <EquationSource Format="TEX">\(\varphi \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>φ</mi> </math></EquationSource> </InlineEquation>. For distributions with at most three offspring, we characterize when this method yields an upper bound, a lower bound, or only an approximation. Because for many distributions it is difficult to get a handle on <InlineEquation ID="IEq8"> <EquationSource Format="TEX">\(S^\infty \)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>S</mi> <mi>∞</mi> </msup> </math></EquationSource> </InlineEquation>, we derive an approximation by series expansion in <i>s</i>, where <i>s</i> is the selective advantage of the mutant. We briefly review well-known asymptotic results that generalize Haldane’s approximation 2<i>s</i> for <InlineEquation ID="IEq9"> <EquationSource Format="TEX">\(S^\infty \)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>S</mi> <mi>∞</mi> </msup> </math></EquationSource> </InlineEquation>, as well as less well-known results on sharp bounds for <InlineEquation ID="IEq10"> <EquationSource Format="TEX">\(S^\infty \)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>S</mi> <mi>∞</mi> </msup> </math></EquationSource> </InlineEquation>. We apply them to explore when bounds for <InlineEquation ID="IEq11"> <EquationSource Format="TEX">\(S^{(n)}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>S</mi> <mrow> <mo stretchy="false">(</mo> <mi>n</mi> <mo stretchy="false">)</mo> </mrow> </msup> </math></EquationSource> </InlineEquation> exist for a family of generalized Poisson distributions. Numerical results demonstrate the accuracy of our and of previously derived bounds for <InlineEquation ID="IEq12"> <EquationSource Format="TEX">\(S^\infty \)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>S</mi> <mi>∞</mi> </msup> </math></EquationSource> </InlineEquation> and <InlineEquation ID="IEq13"> <EquationSource Format="TEX">\(S^{(n)}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>S</mi> <mrow> <mo stretchy="false">(</mo> <mi>n</mi> <mo stretchy="false">)</mo> </mrow> </msup> </math></EquationSource> </InlineEquation>. Finally, we treat an application of these results to determine the response of a quantitative trait to prolonged directional selection.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Bounds for survival probabilities in supercritical Galton-Watson processes and applications to population genetics

  • Reinhard Bürger

摘要

Population genetic processes, such as the adaptation of a quantitative trait to directional selection, may occur on longer time scales than the sweep of a single advantageous mutation. To study such processes in finite populations, approximations for the time course of the distribution of a beneficial mutation were derived previously by branching process methods. The application to the evolution of a quantitative trait requires bounds for the probability of survival \(S^{(n)}\) S ( n ) up to generation n of a single beneficial mutation. Here, we present a method to obtain a simple, analytically explicit, either upper or lower, bound for \(S^{(n)}\) S ( n ) in a supercritical Galton-Watson process. We prove the existence of an upper bound for offspring distributions including Poisson, binomial, and negative binomial. They are constructed by bounding the given generating function, \(\varphi \) φ , by a fractional linear one that has the same survival probability \(S^\infty \) S and yields the same rate of convergence of \(S^{(n)}\) S ( n ) to \(S^\infty \) S as \(\varphi \) φ . For distributions with at most three offspring, we characterize when this method yields an upper bound, a lower bound, or only an approximation. Because for many distributions it is difficult to get a handle on \(S^\infty \) S , we derive an approximation by series expansion in s, where s is the selective advantage of the mutant. We briefly review well-known asymptotic results that generalize Haldane’s approximation 2s for \(S^\infty \) S , as well as less well-known results on sharp bounds for \(S^\infty \) S . We apply them to explore when bounds for \(S^{(n)}\) S ( n ) exist for a family of generalized Poisson distributions. Numerical results demonstrate the accuracy of our and of previously derived bounds for \(S^\infty \) S and \(S^{(n)}\) S ( n ) . Finally, we treat an application of these results to determine the response of a quantitative trait to prolonged directional selection.