[논문 리뷰] Complete intersection for equivariant models

Marta Casanellas, Jesús Fernández-Sánchez외 1명|arXiv (Cornell University)|2015. 12. 22.

Geometric and Algebraic Topology참고 문헌 33인용 수 4

한 줄 요약

이 논문은 생물학적 추론에 핵심적인 무진화 점 주변에서 G-등변 모델의 계통발생 다양체를 국소적으로 기술하는 완전 교차를 명시적으로 구성하는 알고리즘을 제시한다. 이 다양체가 이러한 점에서 매끄럽다는 것을 증명하여, 최대 차수 13 이하의 저차수 다항방정식으로 구성된 최소 시스템이 제로지 열린 근방에서 다양체를 정확히 기술할 수 있음을 보장하며, 이는 이전의 지수적 차수의 구성 방식보다 크게 향상된 결과이다.

ABSTRACT

Phylogenetic varieties related to equivariant substitution models have been studied largely in the last years. One of the main objectives has been finding a set of generators of the ideal of these varieties, but this has not yet been achieved in some cases (for example, for the general Markov model this involves the open "salmon conjecture") and it is not clear how to use all generators in practice. Motivated by applications in biology, we tackle the problem from another point of view. The elements of the ideal that could be useful for applications in phylogenetics only need to describe the variety around certain "points of no evolution". We produce a collection of explicit equations that describe the variety on a Zariski open neighborhood of these points. Namely, for any tree on any number of leaves (and any degrees at the interior nodes) and for any equivariant model on any set of states, we compute the codimension of the corresponding phylogenetic variety. We prove that this variety is smooth at general points of no evolution, and provide an algorithm to produce a complete intersection that describes the variety around these points.

연구 동기 및 목표

• 생물학적 응용을 위해 실용적이고 작은 다항식 생성자 집합으로 계통발생 다양체를 기술하는 데 도전하는 것.
• 유전자가 변화하지 않는 점 주변의 국소적 구조에 집중하는 것—즉, 조상에서 잎으로 이르기까지 서열이 변하지 않는 생물학적으로 의미 있는 구성이다.
• 이러한 점에서 다양체의 완전 교차 기술을 제공하여 필요한 방정식 수를 최소화하는 것.
• 구성된 방정식이 저차수이면서도 불필요한 성분을 포함하지 않고 진정한 계통발생 다양체만을 기술하도록 보장하는 것.
• 아벨 군 기반 모델에 대한 이전 결과를 더 넓은 G-등변 모델 클래스로 일반화하는 것—이에 포함되는 것은 스트랜드 대칭 모델 및 일반 마르코프 모델 등이다.

제안 방법

• 더 작은 나무들을 연결하는 방식으로 기반한 귀납적 접근을 사용하여, 더 작은 나무의 구조에서 계통발생 불변량을 구성한다.
• 유도된 불변량을 보완하기 위해 '박막 평탄화'(thin flattenings)를 도입한다—이것은 내부 노드에서의 기술을 정밀화하는 특수한 종류의 불변량이다.
• G ⊂ 𝔖κ의 표현 이론을 활용하여 다양체의 여피호모지와 국소적 구조를 계산한다.
• 완전 교차를 이루는 엣지 불변량의 정규열을 식별하여, 다양체가 정확히 여피호모지 수만큼의 방정식으로 국소적으로 정의됨을 보장한다.
• 푸리에 기저를 적용하여 불변량을 명시적으로 표현함으로써, 정의 방정식의 알고리즘적 계산을 가능하게 한다.
• 일반적인 무진화 점에서의 매끄러움을 증명하여, 국소적 기술을 위해 완전 교차가 충분함을 보장한다.

실험 결과

연구 질문

• RQ1 임의의 G-등변 모델과 나무 구조에 대해, 일반적인 무진화 점 주변에서 계통발생 다양체를 완전 교차로 국소적으로 기술할 수 있는가?
• RQ2 주어진 나무와 등변 모델에 대해 계통발생 다양체의 여피호모지는 무엇이며, 이는 군 작용과 상태 공간과 어떻게 관련되는가?
• RQ3 무진화 점들이 다양체에서 비특이점인가? 이는 최소 정의 방정식 시스템이 가능함을 의미한다.
• RQ4 정의 방정식은 저차수로 알고리즘적으로 계산 가능하여 생물학적 추론에 실용적으로 적용 가능한가?
• RQ5 구성된 방정식은 진정한 계통발생 다양체에 해당하지 않는 불필요한 기약성 성분을 도입하지 않는가?

주요 결과

• 임의의 G-등변 모델과 임의의 나무에 대해, 일반적인 무진화 점에서 계통발생 다양체는 매끄럽다. 이는 국소적 완전 교차 기술이 가능함을 보장한다.
• 다양체의 여피호모지는 군 G, 잎의 수, 그리고 상태 공간 위에서 G의 순열 표현의 중복도를 이용해 명시적으로 계산된다.
• 완전 교차는 엣지 불변량의 정규열을 사용하여 구성되며, 이 불변량들이 계통발생 불변량임이 알려져 있어 불필요한 성분이 도입되지 않는다.
• 생물학적으로 의미 있는 모델들—예를 들어 스트랜드 대칭 모델 및 일반 마르코프 모델—에 대해 정의 방정식의 차수는 최대 13 이하이며, 이는 이전의 지수적 차수의 구성 방식에 비해 극적으로 향상된 결과이다.
• 알고리즘이 효과적이고 명시적이며, 나무 구조와 군 작용으로부터 국소 정의 방정식을 계산할 수 있다.
• 구성은 가능한 한 낮은 차원의 부분집합을 제외하고 모든 생물학적으로 의미 있는 점을 커버하므로 실용적 추론에 적합하다.