[논문 리뷰] Rigorous results for a population model with selection II: genealogy of the population

Jason Schweinsberg|arXiv (Cornell University)|2015. 07. 01.

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$일 때 사라지는 항들로 유한하게 제한되며, 이는 약한 수렴을 보장한다.
• 증명은 다중 융합 사건—두 개 이상의 계통이 동시에 융합하는 사건—이 주어진 선택 조건 하에서 계통수의 주요 특징임을 확인한다.
• 고정된 시간 영역 $T$에 대해 수렴이 성립하며, 제시된 매개변수 스케일링 하에서 오차 한계는 $N \to \infty$일 때 사라진다.