[논문 리뷰] The stationary distribution of a Wright-Fisher diffusion model with general small mutation rates

Conrad J. Burden, Robert Griffiths|arXiv (Cornell University)|2018. 02. 22.

Mathematical and Theoretical Epidemiology and Ecology Models인용 수 1

한 줄 요약

이 논문은 일반적인 작은 돌연변이 비율을 가진 Wright-Fisher 확산 모델의 정적 분포를 공통 조상 접근법을 사용하여 유도하며, 돌연변이 비율의 일차 항에서 최대 한 번의 돌연변이가 발생할 경우 근사가 성립함을 보여준다. 이 방법은 재조합이 가능한 이진 공통 조상 트리에서의 앨레르기 구성에 대해 명시적인 해를 제공하며, 복제율 λ를 가진 순수한 출생 과정도 포함된다.

ABSTRACT

The stationary distribution of a sample taken from a Wright-Fisher diffusion with general small mutation rates is found using a coalescent approach. The approximation is equivalent to having at most one mutation in the coalescent tree to the first order in the rates. The approach is different from Burden and Tang [1, 2] who use a probability flux argument to obtain the same results from a forward diffusion generator equation. The solution has interest because the solution is not known when rates are not small. An analogous solution is found for the configuration of alleles in a general exchangeable binary coalescent tree. In particular an explicit solution is found for a pure birth process tree when individuals reproduce at rate {\lambda}.

연구 동기 및 목표

• 일반적인 작은 돌연변이 비율 하에서 Wright-Fisher 확산 모델의 정적 분포를 결정하는 것.
• 확산 생성자 기반 기존 방법 외의 공통 조상 프레임워크를 사용한 새로운 유도 방법 제공.
• 일반적인 교환 가능한 이진 공통 조상 트리로의 해를 확장하고 순수한 출생 과정의 구성에 대해 명시적인 해를 구하는 것.
• 돌연변이 비율이 작지 않을 경우 알려진 해가 부족한 점을 메우는 것.

제안 방법

• 공통 조상 접근법을 사용하여 조상 과정에서의 선형 융합 및 돌연변이 사건을 모델링한다.
• 분석은 돌연변이 비율이 작다고 가정하므로, 공통 조상 트리에서 최대 한 번의 돌연변이가 발생하며 일차 항에서의 근사를 가정한다.
• 정적 분포는 공통 조상 트리에서 선형의 앨레르기 구성에 따라 유도된다.
• 이 방법은 재조합이 가능한 이진 공통 조상 트리에 적용되며, 복제율 λ를 가진 순수한 출생 과정의 경우도 포함된다.
• 정방향 확산 생성자로부터 알려진 결과와의 일관성을 통해 해를 검증하지만, 다른 메커니즘을 통해 유도된다.
• 이전 연구에서 사용된 확률 유량 논의를 피함으로써 정적 분포에 대한 새로운 시각을 제공한다.

실험 결과

연구 질문

• RQ1 돌연변이 비율이 작지만 일반적인 경우, Wright-Fisher 확산 모델의 정적 분포는 무엇인가?
• RQ2 공통 조상 접근법은 정적 분포 유도에서 확률 유량 방법과 어떻게 비교되는가?
• RQ3 일반적인 교환 가능한 이진 공통 조상 트리에서 앨레르기 구성에 대해 명시적인 해를 구할 수 있는가?
• RQ4 이 프레임워크 하에서 복제율 λ를 가진 순수한 출생 과정의 정적 분포는 무엇인가?
• RQ5 왜 돌연변이 비율이 작지 않을 경우 해가 알려져 있지 않은데, 이 작업은 어떻게 이 격차를 메우는가?

주요 결과

• 공통 조상 접근법을 사용하여 작은 돌연변이 비율의 일차 항에서 유효한 정적 분포가 유도되었다.
• 해는 공통 조상 트리에서 최대 한 번의 돌연변이가 발생한다는 가정과 동치이며, 분석을 단순화한다.
• 일반적인 교환 가능한 이진 공통 조상 트리에서의 앨레르기 구성에 대해 명시적인 해를 확보하였다.
• 복제율 λ를 가진 순수한 출생 과정의 경우, 정적 분포에 대해 닫힌 형태의 해를 도출하였다.
• 정방향 확산 생성자 기반 방법과 독립적인 새로운 유도를 제공하여 보완적인 통찰을 제공한다.
• 이 방법은 돌연변이 비율이 작지 않을 경우 기존 해의 한계를 부각하며, 현재 결과의 독창성을 강조한다.