[论文解读] Complete intersection for equivariant models

Marta Casanellas, Jesús Fernández-Sánchez+1 more|arXiv (Cornell University)|Dec 22, 2015

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-等变模型类，包括链对称模型和一般马尔可夫模型。

提出的方法

• 基于将较小的树连接起来的归纳方法，从更小的树结构中构建系统发育不变量。
• 通过引入‘细扁平化’——一种专门的不变量类——来补充诱导不变量，以在内部节点处细化描述。
• 利用群 G⊂𝔖κ 的表示理论来计算簇的余维数和局部结构。
• 识别出一组构成正则序列的边不变量，形成完整交集，从而保证簇在局部由恰好余维数个方程定义。
• 应用傅里叶基底显式表示不变量，从而实现定义方程的算法计算。
• 证明在通用无进化点处簇是光滑的，这保证了完整交集足以实现局部描述。

实验结果

研究问题

• RQ1 对于任意 G-等变模型和树拓扑，是否能用完整交集在通用无进化点附近局部描述系统发育簇？
• RQ2 给定树和等变模型，系统发育簇的余维数是多少？它与群作用和状态空间有何关系？
• RQ3 无进化点在簇上是否为非奇异点，从而允许使用极小的定义方程组？
• RQ4 能否以低次方程算法计算出定义方程，使其适用于生物学中的实际推断？
• RQ5 所构造的方程是否会引入不对应于真实系统发育簇的虚假不可约分支？

主要发现

• 对于任意 G-等变模型和任意树，系统发育簇在通用无进化点处是光滑的，从而保证了局部完整交集描述的可行性。
• 簇的余维数可显式表示为群 G、叶数以及 G 在状态空间上的置换表示的重数的函数。
• 通过一组已知为系统发育不变量的正则序列边不变量构造了完整交集，从而确保未引入虚假分支。
• 对于生物学上相关的模型（如链对称模型和一般马尔可夫模型），定义方程的次数最高为13，相比指数次构造方法有巨大改进。
• 该算法是有效且显式的，可从树结构和群作用计算出局部定义方程。
• 该构造覆盖了除可能一个低维子集外的所有生物学上有意义的点，因此适用于实际推断。