[논문 리뷰] A Hamilton-Jacobi approach to characterize the evolutionary equilibria in heterogeneous environments

Sepideh Mirrahimi|arXiv (Cornell University)|2016. 12. 19.

Evolution and Genetic Dynamics참고 문헌 32인용 수 26

한 줄 요약

이 논문은 변이, 선택 및 이동이 존재하는 두 개의 이질적 서식지에서 생체형 구조를 가진 집단의 진화적 평형을 기술하기 위해 보정항을 포함한 해밀턴-자코비 접근법을 개발한다. 단형성 또는 이형성 조건을 명시적으로 제시하고, 가우시안 가정을 초월한 고차항 근사치를 계산함으로써, 고전적 양적 유전학 모델에 비해 더 엄밀한 수학적 기반을 제공한다.

ABSTRACT

In this work, we characterize the solution of a system of elliptic integro-differential equations describing a phenotypically structured population subject to mutation, selection and migration between two habitats. Assuming that the effects of the mutations are small but nonzero, we show that the population's distribution has at most two peaks and we give explicit conditions under which the population will be monomorphic (unimodal distribution) or dimorphic (bimodal distribution). More importantly, we provide a general method to determine the dominant terms of the population's distribution in each case. Our work, which is based on Hamilton-Jacobi equations with constraint, goes further than previous works where such tools were used, for different problems from evolutionary biology, to identify the asymptotic solutions, while the mutations vanish, as a sum of Dirac masses. In order to extend such results to the case with non-vanishing effects of mutations, the main elements are a uniqueness property and the computation of the correctors. This method allows indeed to go further than the Gaussian approximation commonly used by biologists and makes a connection between the theories of adaptive dynamics and quantitative genetics. Our work being motivated by biological questions, the objective of this article is to provide the mathematical details which are necessary for our biological results [16].

연구 동기 및 목표

• 변이, 선택 및 이동 조건 하에서 이질적 환경에서의 생체형 분포를 특성화하는 수학적 프레임워크를 제공하는 것.
• 소멸하는 돌연변이의 극한에 국한되지 않고, 고유성과 다음 계수 보정항을 계산함으로써 해밀턴-자코비 접근법을 확장하는 것.
• 양적 유전학에서 널리 사용되는 가우시안 근사치에 대한 엄밀한 대체 모델을 제공하여, 진화적 결과의 더 정확한 예측을 가능하게 하는 것.
• 두 서식지 모델에서 집단 분포가 단형성(단봉) 또는 이형성(이봉)이 되는 조건을 명시적으로 규명하는 것.
• 비점근적 해를 위한 비점근적 해법을 도출함으로써 적응역학과 양적 유전학 사이의 격차를 메우는 것.

제안 방법

• 두 서식지 간의 인구 밀도, 성장률 및 이동을 모델링하는 타원형 적분미분방정식계를 수립한다.
• 인구 밀도 방정식을 해밀턴-자코비 방정식으로 변환하기 위해 로그(홉프-콜) 변환을 적용한다.
• 점성해의 프레임워크를 수립하고, 소멸하지 않는 돌연변이 효과 하에서 점성해의 고유성을 증명한다.
• 로그 인구 밀도 전개에서 1차 및 2차 보정항(v_i, w_i)을 계산하여 고차항 보정을 포착한다.
• 평형점 주변의 테일러 전개를 사용하여 인구 크기 보정(K_i) 및 보정 함수의 명시적 표현을 유도한다.
• 한계 분포의 지지집합 분석 및 해밀턴-자코비 해의 임계점에서의 행동 분석을 통해 단형성 및 이형성 조건을 유도한다.

실험 결과

연구 질문

• RQ1 두 개의 이질적 서식지에서 생체형 구조를 가진 집단이 어떤 조건에서 단형성 또는 이형성 분포로 진화하는가?
• RQ2 소멸하는 돌연변이 극한을 초월하여 비점근적 진화적 평형을 기술할 수 있도록 해밀턴-자코비 접근법을 어떻게 확장할 수 있는가?
• RQ3 돌연변이 효과가 비영일 경우, 인구 분포의 다음 계수 보정항(보정항)의 수학적 구조는 어떠한가?
• RQ4 이동률과 서식지별 선택 강도는 이형성의 발생에 어떻게 영향을 미치는가?
• RQ5 이 해밀턴-자코비 프레임워크를 통해 양적 유전학에서 사용되는 가우시안 근사치를 체계적으로 보정할 수 있는가?

주요 결과

• 인구 분포는 최대 두 개의 피크를 가지며, 단형성 또는 이형성 조건은 이동률, 선택 강도 및 경쟁 강도를 포함한 부등식을 통해 명시적으로 특성화된다.
• 제약 조건이 있는 해밀턴-자코비 방정식의 점성해에 대한 고유성 결과를 제공하여 엄밀한 고차항 근사치를 가능하게 한다.
• 보정항 v_i 및 w_i가 명시적으로 계산되어, O(ε²)까지의 인구 크기 보정 K_i를 도출할 수 있다.
• 단형성의 경우, 인구 크기 보정 K_2는 K_2 = N_2^{M*} ( (E_2 + 0.5 D_2²)/√g_1 + F_2 )로 유도되며, F_2는 방정식계를 통해 결정된다.
• 이 방법은 이론적 가우시안 가정 없이도 인구 분포의 주요 항을 성공적으로 계산하여, 이전에 진화 생물학에서 사용된 근사치를 보완한다.
• 분석 결과, ESS가 두 개의 서로 다른 특성을 지지할 경우 이형성이 발생하며, 임계점(z = -θ)에서의 일致성 검증을 통해 해 구조가 검증된다.