[논문 리뷰] Stochastic Models for Speciation Events in Phylogenetic trees

Tanja Gernhard|arXiv (Cornell University)|2006. 10. 30.

Evolution and Paleontology Studies참고 문헌 11인용 수 4

한 줄 요약

이 학위논문은 뿌리 박힌 계통수에서의 종 분화 사건을 위한 확률 모델을 개발하며, 클러백-라이블러 발산과 마팅게일 이론을 사용하여 유일 모델과 균일 모델을 비교한다. 정점의 확률 분포를 위한 랭크 함수를 도입하고, 이러한 모델을 이용해 진화 가설을 검증하며, 특히 유일 나무에서 간선 길이를 추정하고 영장류 계통수의 진화가 유일 과정을 따르는지 평가한다.

ABSTRACT

In a phylogenetic tree, we often don't have information about the time a speciation event (inner node) occured. Under a neutral model for speciation, I develop fast algorithms for calculating the probability that an inner node i is the k-th speciation event. For the Yule and the coalescent model, I develop an edge length estimation as well. Various properties of the Yule model are discussed throughout the thesis.

연구 동기 및 목표

• 확률 과정, 특히 유일 모델과 균일 모델을 사용하여 계통수에서의 종 분화 사건을 모델링하는 것.
• 가설 검정과 정보 이론적 측정치를 사용하여 관측된 계통수의 유일 모델에 대한 적합도를 평가하는 것.
• 다른 모델 하에서 나무의 내부 정점의 확률을 정량화하기 위한 랭크 함수를 개발하는 것.
• 랭크 함수 기반의 분석적 및 알고리즘적 방법을 사용하여 유일 나무에서 간선 길이를 추정하는 것.
• 실제 생물학적 데이터, 예를 들어 영장류 슈퍼트리에 이러한 모델을 적용하여 진화 가설을 검증하는 것.

제안 방법

• 일정한 종 분화 비율을 가진 순수 출생 과정인 유일 모델과 모든 나무 형태가 동일하게 확률을 가지는 균일 모델을 사용한다.
• 유일 모델과 균일 모델 간의 비교를 위해 클러백-라이블러 발산을 적용하여 모델 식별 능력을 정량화한다.
• 균일 모델 하에서 나무 형태 통계량의 편차를 유한하게 제한하기 위해 마팅게일 이론과 아즈마 부등식을 활용한다.
• 정점이 나무에서 주어진 랭크에 있을 확률을 계산하기 위한 랭크 함수를 도입하여 나무 구조에 대한 가설 검정을 가능하게 한다.
• 이진 및 비이진 나무에서 정점 랭크 분포를 계산하기 위한 다항 시간 알고리즘(예: RankProb, RankCount)을 개발한다.
• 유사도 비율 검정과 네이만-피어슨 보조정리를 사용하여 주어진 나무가 유일 모델이 아닌지 균일 모델에서 유래했는지 검증한다.

실험 결과

연구 질문

• RQ1 클러백-라이블러 발산이 계통수에서 유일 모델과 균일 모델을 효과적으로 구분할 수 있는가?
• RQ2 유일 또는 균일 나무에서 내부 정점의 랭크 확률 분포는 어떻게 되는가?
• RQ3 마팅게일 과정을 사용하여 균일 모델 하에서 나무 형태 통계량의 편차를 어떻게 제한할 수 있는가?
• RQ4 영장류 슈퍼트리는 어느 정도까지 유일 모델의 종 분화 메커니즘을 지지하는가?
• RQ5 랭크 함수와 관련 확률을 사용하여 유일 나무에서 간선 길이를 분석적으로 추정할 수 있는가?

주요 결과

• 유일 모델과 균일 모델 간의 클러백-라이블러 발산은 유계이면서 계산 가능하며, 유일 모델은 균일 모델보다 더 높은 엔트로피를 보인다.
• 유일 모델 하에서 나무의 확률은 각 내부 정점의 서브트리 크기의 역수의 곱에 비례하며, 이는 정리 2.2.1에 의해 정밀하게 기술된다.
• 랭크 함수를 통해 정점의 확률 분포를 효율적으로 계산할 수 있으며, RankProb 및 RankProbGen과 같은 알고리즘이 다항 시간 내에 실행된다.
• 랭크 함수 기반의 유사도 비율 검정은 주어진 나무가 유일 모델에 부합하는지 여부를 통계적으로 타당한 방법으로 검증할 수 있다.
• 영장류 슈퍼트리는 균일 모델과는 상당히 다른 나무 형태 분포를 보이며, 이는 유일 모델이 종 분화의 타당한 메커니즘임을 지지한다.
• 내부 정점의 기대 랭크를 사용하여 유일 나무에서 간선 길이를 분석적으로 추정할 수 있으며, 시뮬레이션 및 알고리즘 구현을 통한 검증 결과가 확인되었다.