[論文レビュー] Rigorous results for a population model with selection II: genealogy of the population

Jason Schweinsberg|arXiv (Cornell University)|Jul 1, 2015

Evolution and Genetic Dynamics被引用数 7

ひとこと要約

この論文は、突然変異率と選択係数の特定のスケーリング極限のもとで、弱い正の選択が働く大規模で固定サイズの集団の系統が、ボルタウーゼン＝シュニットマン共析に収束することを厳密に確立している。拡散近似とカップリング技法を用いて、有益な突然変異の急速な固定によって複数の系統が同時に合流する現象が、強い選択下の集団の特徴であるという非厳密な予測を裏付けた。

ABSTRACT

We consider a model of a population of fixed size $N$ undergoing selection. Each individual acquires beneficial mutations at rate $μ_N$, and each beneficial mutation increases the individual's fitness by $s_N$. Each individual dies at rate one, and when a death occurs, an individual is chosen with probability proportional to the individual's fitness to give birth. Under certain conditions on the parameters $μ_N$ and $s_N$, we show that the genealogy of the population can be described by the Bolthausen-Sznitman coalescent. This result confirms predictions of Desai, Walczak, and Fisher (2013), and Neher and Hallatschek (2013).

研究の動機と目的

• 選択下の集団の系統がボルタウーゼン＝シュニットマン共析に収束するという非厳密な予測に対して、数学的に厳密な根拠を提供すること。
• 適応度分布と適応の速度に関する先行研究を拡張し、集団の祖先過程を分析すること。
• ボルタウーゼン＝シュニットマン過程が共析極限として得られるように、突然変異率 $\mu_N$ と選択強度 $s_N$ に課すべき条件を確立すること。
• 大規模集団において有益な突然変異の急速な固定が原因で生じる、同時に複数の系統が合流する現象が系統に与える役割を確認すること。

提案手法

• 時間の再スケーリングに $a_N = \frac{\log(s_N/\mu_N)}{s_N}$ を用いた拡散スケーリングを採用し、系統の合流に伴う自然な時間スケールを捉える。
• 系統的過程は、合流イベントをモデル化するポアソン点過程 $\Phi$ にカップリングされた共析過程 $\Pi_N$ の列によって近似される。
• 主な技術的道具は、有限集団の過程 $\Pi_N$ と極限のボルタウーゼン＝シュニットマン共析 $\Pi$ の間のカップリングであり、誤差項を小さく制御することで乖離を制御する。
• とくに、まれな同時合流のための確率的乖離を、モーメントの上限とポアソン尾部推定を用いて境界化することで、証明が成立する。
• 結果の構築には、先行研究 [27] における適応度のダイナミクスと突然変異蓄積の結果を活用し、特に $k_N = \frac{\log N}{\log(s_N/\mu_N)}$ として与えられる典型的な突然変異数に着目する。
• 時間ずらしの議論により、仮定 A1–A3 の下で、有限時間における系統的過程 $\Pi_N(t)$ が固定時間点で弱収束して $\Pi(t)$ に近づくことが保証される。

実験結果

リサーチクエスチョン

• RQ1 突然変異率 $\mu_N$ と選択係数 $s_N$ にどのような条件下で、選択下の集団の系統がボルタウーゼン＝シュニットマン共析に収束するか。
• RQ2 突然変異率と選択係数のスケーリング極限が、祖先過程における複数合流イベントの頻度と構造にどのように影響を与えるか。
• RQ3 デサイ、ウォルツァック、ファイザー（2013年）およびネーラーとハラーツィャック（2013年）の非厳密な予測を、確率過程とカップリング技法を用いて形式的に正当化できるか。
• RQ4 典型的な突然変異数 $k_N$ が、共析イベントの時間スケールと構造を決定づける役割を果たすか。

主な発見

• 仮定 A1–A3 の下で、$N \to \infty$ のとき、集団の系統は弱収束してボルタウーゼン＝シュニットマン共析に近づく。
• 時間スケール $a_N = \frac{\log(s_N/\mu_N)}{s_N}$ が、極限共析挙動を観察するのに適切なスケーリングであることが示された。
• 突然変異数 $k_N = \frac{\log N}{\log(s_N/\mu_N)}$ は、最も適応度の高い個体と平均個体との間の典型的な適応度差を特徴づけ、共析過程の分解能を制御する。
• 有限過程 $\Pi_N$ と極限過程 $\Pi$ 間の不一致合流イベントの確率は、$N \to \infty$ のとき消える項によって境界化されており、弱収束が保証される。
• 証明により、2つ以上の系統が同時に合流する複数合流イベントが、与えられた選択状態下での系統の支配的特徴であることが確認された。
• 固定時間区間 $T$ に対して収束が成り立ち、パrameterスケーリングが所定の条件下にある限り、誤差項は $N \to \infty$ のとき消える。