[논문 리뷰] Incomplete Lineage Sorting: Consistent Phylogeny Estimation From Multiple Loci

Elchanan Mossel, Sébastien Roch|arXiv (Cornell University)|2007. 10. 01.

Genomics and Phylogenetic Studies참고 문헌 14인용 수 6

한 줄 요약

이 논문은 불완전한 선형 분할에 기인한 유전자 트리의 충돌을 고려하여 종 트리 추정을 위한 통계적으로 일관된 알고리즘인 GLASS(Global LAteSt Split)를 소개한다. 이 알고리즘은 다양한 유전자 위치에서의 최소 공통 조상 형성 시간을 활용하여 올바른 종 트리 구조를 추론하며, 표준 공통 조상 모델 하에서 일관성을 입증하고, 중간 정도의 추정 오차에 대해서도 내성적 성능을 유지하며 명시적인 수렴 속도를 제공한다.

ABSTRACT

We introduce a simple algorithm for reconstructing phylogenies from multiple gene trees in the presence of incomplete lineage sorting, that is, when the topology of the gene trees may differ from that of the species tree. We show that our technique is statistically consistent under standard stochastic assumptions, that is, it returns the correct tree given sufficiently many unlinked loci. We also show that it can tolerate moderate estimation errors.

연구 동기 및 목표

• 불완전한 선형 분할로 인해 유전자 트리가 진정한 종 트리와 다를 수 있는 문제를 해결하기 위해.
• 유전자 트리가 공통 조상의 확률적 변동으로 인해 상충될 경우에도 통계적으로 일관되면서 계산 효율성이 높은 방법을 개발하기 위해.
• 공통 조상 모델 하에서 종 트리 추정의 명시적 수렴 속도를 제공하여 재구성 품질의 실용적 평가를 가능하게 하기 위해.
• 분자 시계를 가정하지 않고도 거리 기반 재구성에 일반화할 수 있도록, 다수의 유전자 위치에 대한 구조적 일치 방법을 확장하기 위해.

제안 방법

• GLASS는 각 종 집단 쌍에서의 모든 유전자 트리와 개체군의 최소 분리 시점을 취하여 집단 간 공통 조상 형성 시간을 추정한다.
• 이 방법은 분기 길이 기준 임계값을 활용한 조합 조건을 사용하여 이러한 최소 공통 조상 형성 시간에서 정확한 종 트리 구조를 식별한다.
• 이 방법은 공통 조상 모델 하에서, 두 종 집단 간의 가장 짧은 공통 조상 형성 시간이 그들의 종 트리에서 가장 최근의 공통 조상에 해당한다는 성질에 기반한다.
• 공통 조상 형성 시간의 추정 오차가 제한된 범위 내에 있을 경우, 이 방법은 노이즈가 있는 데이터에 대해서도 일반화되어 정확성을 유지한다. 오차가 가장 짧은 분기 길이와 관련된 임계값 이하일 경우에만 정확성이 유지된다.
• 이 알고리즘은 어떤 거리 기반 계통수 재구성 방법과도 호환되며, 분자 시계를 가정하지 않아도 된다.
• 이 방법은 다수의 유전자 위치로 확장된 로즈버그의 구조적 일치 방법을 활용하여, 여러 유전자를 통해 공통 조상 형성 시간 정보를 통합함으로써 종 수준의 관계를 추론한다.

실험 결과

연구 질문

• RQ1 불완전한 선형 분할로 인해 유전자 트리가 상충될 경우, 단순하고 계산 효율성이 높은 방법이 통계적으로 일관된 종 트리 추정을 달성할 수 있는가?
• RQ2 다수의 유전자 위치에서 최소 공통 조상 형성 시간이 진정한 종 트리의 구조를 올바르게 식별할 수 있는 조건은 무엇인가?
• RQ3 공통 조상 형성 시간 측정의 유한한 추정 오차 하에서 이 방법의 성능은 어떠한가?
• RQ4 이 방법은 분자 시계 가정 없이도 거리 기반 방법으로 일반화될 수 있는가?
• RQ5 공통 조상 모델 하에서 종 트리 재구성에 대해 유도할 수 있는 명시적 수렴 속도는 무엇인가?

주요 결과

• GLASS는 표준 공통 조상 모델 하에서 통계적으로 일관되며, 서로 연결되지 않은 유전자 위치가 충분히 많을 경우 높은 확률로 정확한 종 트리를 재구성한다.
• 오차가 종 트리에서 가장 짧은 분기 길이의 반 이하일 경우, 공통 조상 형성 시간의 중간 정도의 추정 오차가 있더라도 방법은 여전히 일관성을 유지한다.
• 명시적인 수렴 속도가 유도되었다: 잘못된 재구성 확률은 유전자 위치 수에 따라 지수적으로 감소하며, $ \mathbb{P}[\text{Incorrect Reconstruction}] \leq (n-1)[e^{-\mu\Lambda}]^k $로 유계된다. 여기서 $ \Lambda = \frac{m - 2\delta}{m} $이다.
• 임의의 원하는 오차 허용 범위 $ \varepsilon $에 대해, $ k = \frac{1}{\mu\Lambda} \ln\left(\frac{n-1}{\varepsilon}\right) $개의 유전자 위치를 사용함으로써 높은 확률로 정확한 재구성을 달성할 수 있다.
• GLASS는 어떤 거리 행렬 방법과도 조합 가능하며, 종 트리의 각 분기 내에서 돌연변이율이 일관될 경우 분자 시계 가정 없이도 사용 가능하다.
• 이 방법은 다수의 유전자 위치로 확장된 구조적 일치를 일반화하여, 유전자 트리의 상충성에서 기인한 종 트리 추정을 위한 강력한 프레임워크를 제공한다.