[논문 리뷰] A stochastic model for phylogenetic trees

T. M. Liggett, R. B. Schinazi|arXiv (Cornell University)|2009. 05. 14.

Evolution and Genetic Dynamics인용 수 5

한 줄 요약

이 논문은 적합도 기반 선택을 갖는 출생-사망 과정을 기반으로 한 단순한 확률 모델을 제안한다: 새로운 유형(돌연변이)은 $\lambda n$의 비율로 추가되고, 가장 적합도가 낮은 유형은 $n$의 비율로 제거된다. $\lambda < 1$일 때 모델은 장기간에 걸쳐 단일 우세 유형을 가진 가늘고 길쭉한 나무를 생성한다—인플루엔자와 일치한다. $\lambda > 1$일 때는 많은 공존 유형을 가진 덤불 같은 나무를 생성한다—HIV와 일치한다. 핵심 결과는, $\lambda < 1$일 때 우세 유형이 시간 $\alpha t$에서 $t$까지 지속될 확률이 $\alpha$로 수렴하는 반면, $\lambda > 1$일 때는 0으로 수렴한다는 것이다.

ABSTRACT

We propose the following simple stochastic model for phylogenetic trees. New types are born and die according to a birth and death chain. At each birth we associate a fitness to the new type sampled from a fixed distribution. At each death the type with the smallest fitness is killed. We show that if the birth (i.e. mutation) rate is subcritical we get a phylogenetic tree consistent with an influenza tree (few types at any given time and one dominating type lasting a long time). When the birth rate is supercritical we get a phylogenetic tree consistent with an HIV tree (many types at any given time, none lasting very long).

연구 동기 및 목표

• 돌연변이, 사망, 적합도 선택 이외의 요소 없이도 인플루엔자와 HIV의 대조적인 계통수 나무 구조를 재현할 수 있는 최소한의 확률 모델을 개발하기 위해.
• 동일한 기반이 되는 메커니즘—돌연변이—가 어떻게 한 유형이 우세한 계통수(인플루엔자)와 여러 유형이 공존하는 계통수(HIV)를 만들어내는지 설명하기 위해.
• 다양한 돌연변이 비율에서 우세 유형의 장기적 지속 행동을 분석하기 위해.
• 돌연변이 비율 $\lambda$와 결과 계통수의 구조적 특성 사이의 수학적 연결 고리를 설정하기 위해.

제안 방법

• 유형 수에 대한 연속시간 출생-사망 과정으로 모델의 진화를 기술: 출생 비율 $n\lambda$, 사망 비율 $n$ (단, $n \geq 2$).
• 각 신규 유형에 대해 $[0,1]$에서 균일 분포를 따르는 독립적인 적합도 값을 할당한다.
• 각 사망 사건에서 가장 낮은 적합도를 가진 유형을 제거함으로써 적합도 순위 기반 선택을 보장한다.
• 임의의 시점에서 가장 높은 적합도를 가진 유형을 최대 유형으로 정의하고, 시간 간격 동안의 지속성을 추적한다.
• 재생 이론과 재생-보상 과정을 사용하여, 단일 유형으로 복귀하는 시간 간격($\lambda < 1$) 또는 새로운 최대치에 도달하는 첫 번째 통과 시간($\lambda > 1$)을 분석한다.
• 강한 대수법칙, 분포 수렴, 안정한 극한 정리(예: 인덱스 $1/2$의 안정적 하위조직)를 적용하여 渐近적 행동을 도출한다.

실험 결과

연구 질문

• RQ1 단순한 확률 모델(돌연변이, 사망, 적합도 기반 선택 포함)이 인플루엔자와 HIV의 상이한 계통수 나무 형태를 재현할 수 있는가?
• RQ2 돌연변이 비율 $\lambda$는 나무 내에서 최우세(가장 적합한) 유형의 지속성에 어떤 영향을 미치는가?
• RQ3 시간 $t \to \infty$일 때, 시간 $\alpha t$의 최우세 유형이 시간 $t$까지 지속될 한계 확률은 무엇인가?
• RQ4 모델의 행동은 $\lambda < 1$, $\lambda = 1$, $\lambda > 1$일 때 어떻게 정성적으로 변화하는가?
• RQ5 인플루엔자 나무는 가늘고, HIV 나무는 덤불 같은데, 둘 다 돌연변이에서 기인함에도 불구하고 이러한 차이를 설명할 수학적 메커니즘이 존재하는가?

주요 결과

• $\lambda < 1$일 때, 시간 $\alpha t$의 최우세 유형이 시간 $t$까지 지속될 확률은 $t \to \infty$일 때 $\alpha$로 수렴하며, 이는 가장 적합한 계통의 장기적 지속성을 나타낸다.
• $\lambda > 1$일 때, 시간 $\alpha t$에서 $t$까지 최우세 유형이 지속될 확률은 0으로 수렴하며, 이는 빠른 교체와 다수의 공존 유형을 나타낸다.
• $\lambda = 1$일 때, 지속 확률의 극한도 역시 $\alpha$이며, 이는 $S_{N(\alpha t)} / S_{N(t)}$의 수렴을 통해 유도된다.
• $\lambda = 1$일 때, 유형 수 $S_n$은 분포적으로 $n^2$만큼 증가하며, 이는 인덱스 $1/2$의 안정 분포와 일치한다.
• $\lambda < 1$일 때 지속 확률이 $\alpha$로 수렴하는 것은 재생 이론과 유계 수렴 정리를 통해 입증된다.
• $\lambda > 1$일 때 $S_n / n$이 유한한 극한으로 수렴하는 것은 랜덤 워크의 쌍용 및 지배 수렴 정리를 사용하여 증명되며, 이는 지속 확률이 극한에서 0이 되도록 보장한다.