[논문 리뷰] The continuous limit of the Moran process and the diffusion of mutant genes in infinite populations

Fabio A. C. C. Chalub, Max O. Souza|ArXiv.org|2006. 02. 23.

Stochastic processes and statistical mechanics참고 문헌 37인용 수 7

한 줄 요약

이 논문은 무한 모집단에서의 빈도 의존 Moran 과정의 연속 근사(limit)를 유도하며, 이산 확률적 역학이 이완-확산 스케일링 하에서 퇴화된 포물형 편미분방정식(PDE)로 수렴함을 보여준다. 주요 결과는 이 PDE 가 키무라의 고정 방정식을 일반화하며, 빈도 독립적 적합도의 경우 형식적으로 이를 복원함을 보이며, 큰 모집단에서 잘 정의되어 있고 수렴함을 증명한다.

ABSTRACT

We consider the so called Moran process with frequency dependent fitness given by a certain pay-off matrix. For finite populations, we show that the final state must be homogeneous, and show how to compute the fixation probabilities. Next, we consider the infinite population limit, and discuss the appropriate scalings for the drift-diffusion limit. In this case, a degenerated parabolic PDE is formally obtained that, in the special case of frequency independent fitness, recovers the celebrated Kimura equation in population genetics. We then show that the corresponding initial value problem is well posed and that the discrete model converges to the PDE model as the population size goes to infinity. We also study some game-theoretic aspects of the dynamics and characterize the best strategies, in an appropriate sense.

연구 동기 및 목표

• 빈도 의존적 적합도를 갖는 유한 모집단에서의 Moran 과정의 열역학적 극한을 수립하기 위해.
• 큰 모집단에 대해 연속 PDE 근사를 얻기 위한 적절한 스케일링(이완 및 이완-확산)을 규명하기 위해.
• 유도된 퇴화된 포물형 PDE의 잘 정의됨(well-posedness)을 증명하고, 이산 모델이 무한 모집단 극한에서 PDE로 수렴함을 보장하기 위해.
• 큰 유한 모집단에서 진화적으로 안정된 전략(ESS)과 고정 역학 사이의 관계를 명확히 하기 위해.
• 빈도 독립적 적합도 조건 하에서 유도된 PDE 가 키무라의 고전적 방정식으로 복원됨을 보장하기 위해.

제안 방법

• 크기 $N$ 의 유한 모집단에서 빈도 의존적 선택을 모델링하기 위해 지급 행렬을 사용한 Moran 과정을 사용한다.
• 두 가지 시간 스케일링 체제를 적용한다: $\Delta t = 1/N$ (이완 극한) 및 $\Delta t = 1/N^2$ (이완-확산 극한), 후자의 경우 약한 선택을 적용한다.
• 모집단이 무한해지는 극한에서 퇴화된 포물형 PDE 를 도출하며, 이 때 확산 계수는 경계 $x=0$ 과 $x=1$ 에서 0이 된다.
• 형식적 점근 전개와 약한 공식화를 사용하여 PDE 를 유도하고 경계 조건을 확인한다.
• 부분 적분과 시험 함수 분석을 통해 경계 항의 구조를 규명하고 PDE 의 일관성을 확인한다.
• 초기값 문제의 잘 정의됨을 확립하고, $N \to \infty$ 일 때 이산 과정이 PDE 로 수렴함을 증명한다.

실험 결과

연구 질문

• RQ1 빈도 의존적 적합도를 갖는 유한 모집단에서의 Moran 과정이 무한 모집단 극한에서 어떻게 연속 PDE 로 수렴하는가?
• RQ2 비어 있지 않은 확산 극한을 얻기 위해 필요한 시간 스케일링은 무엇이며, 이는 유도된 PDE 에 어떤 영향을 미치는가?
• RQ3 유도된 PDE 가 빈도 독립적 적합도 조건에서 어떻게 키무라의 고전적 방정식으로 축소되는가?
• RQ4 연속 복제자 방정식의 진화적으로 안정된 전략(ESS)은 큰 유한 모집단에서의 고정 확률과 어떻게 관련이 있는가?
• RQ5 경계 확산은 고정 또는 상실을 보장하는 데 어떤 역할을 하는가? 그리고 이는 어떻게 퇴화된 PDE 에서 포착되는가?

주요 결과

• 무한 모집단 극한에서의 Moran 과정은 시간에 따른 돌연변이 유전자 빈도의 확률 분포를 모델링하는 퇴화된 포물형 PDE 를 도출한다.
• 이완-확산 스케일링($\Delta t = 1/N^2$)과 약한 선택 하에서 PDE 는 이완과 확산 효과를 모두 포괄하며, 확산 계수는 $x=0$ 과 $x=1$ 에서 0이 된다.
• 빈도 독립적 적합도의 경우 PDE 는 형식적으로 키무라의 방정식으로 축소되며, 고정 확률에 대한 잘 알려진 전진-후진 쌍의 방정식을 복원한다.
• 유도된 PDE 의 초기값 문제는 잘 정의되어 있으며, 적절한 조건 하에서 해의 존재성과 유일성을 보장한다.
• 이산 Moran 과정은 모집단 크기 $N \to \infty$ 일 때 PDE 모델로 수렴함을 보여, 연속 근사의 타당성을 확립한다.
• 큰 유한 모집단에서 최적의 전략은 진화적으로 안정된 전략(ESS)에 해당하며, 이는 복제자 방정식의 역학적 행동과 일치한다.