[论文解读] Perfect Phylogeny Haplotyping is Complete for Logspace

Michael Elberfeld|ArXiv.org|May 5, 2009

Genetic Associations and Epidemiology参考文献 8被引用 3

一句话总结

本文证明了完美系谱单倍型分析（PPH）属于确定性对数空间（L）完全问题，解决了长期悬而未决的开放问题。论文提出一种基于图的新表征方法，利用从基因型数据导出的解析图，将PPH归约为二分图识别问题，而该问题已知属于L，从而通过一阶归约和空间高效计算证明了PPH的L完全性。

ABSTRACT

Haplotyping is the bioinformatics problem of predicting likely haplotypes based on given genotypes. It can be approached using Gusfield's perfect phylogeny haplotyping (PPH) method for which polynomial and linear time algorithms exist. These algorithm use sophisticated data structures or do a stepwise transformation of the genotype data into haplotype data and, therefore, need a linear amount of space. We are interested in the exact computational complexity of PPH and show that it can be solved space-efficiently by an algorithm that needs only a logarithmic amount of space. Together with the recently proved L-hardness of PPH, we establish L-completeness. Our algorithm relies on a new characterization for PPH in terms of bipartite graphs, which can be used both to decide and construct perfect phylogenies for genotypes efficiently.

研究动机与目标

• 解决关于完美系谱单倍型分析（PPH）是否属于确定性对数空间（L）的开放问题。
• 开发一种在O(log n)空间内运行的空间高效算法，优于以往的多项式空间方法。
• 通过编码基因型中2-位点间配对解析模式的解析图，提出PPH的新表征方式。
• 建立从PPH到二分图识别问题的一阶归约，利用解析图的结构特性。
• 证明PPH为L完全问题，结合新提出的L成员性结果与先前的L难结果。

提出的方法

• 为每个基因型构建解析图，其中顶点表示SNP位点，边表示2-位点间解析关系（相等或不相等）。
• 根据基因型是否以相等（0）或不相等（1）方式解析2-位点，为边分配权重（0或1）。
• 通过将权重为0的边替换为长度为2的路径，将加权图转换为无权图，同时保持路径长度的奇偶性。
• 将PPH归约为二分图问题：当且仅当变换后的图是二分图时，存在完美系谱。
• 使用一阶逻辑表达图的构建与二分性检查，确保归约属于L。
• 利用已知事实：二分性可在L中判定，从而将PPH置于L中。

实验结果

研究问题

• RQ1 完美系谱单倍型分析（PPH）是否可在确定性对数空间（L）中求解？
• RQ2 能否构建一种图论表征，避免猜测解析关系，而改用预计算的关系？
• RQ3 完美系谱的存在是否精确对应于变换后解析图的二分性？
• RQ4 能否通过一阶归约将PPH问题归约为二分图识别问题？
• RQ5 在已有L难性结果与新L成员性结果的基础上，PPH是否为L完全问题？

主要发现

• 证明了PPH属于确定性对数空间（L），解决了Elberfeld（2009）提出的主要开放问题。
• 提出了一种基于解析图的新PPH表征方式，用于编码基因型中2-位点的配对解析模式。
• 该问题可归约为检查变换后的无向图是否为二分图，而该问题已知属于L。
• 归约过程为一阶归约，保持对数空间可计算性，支持高效且空间最优的计算。
• 结合先前的L难性结果，本研究确立了PPH为L完全问题，意味着其属于对数空间中最具挑战性的问题之列。
• 该方法避免使用复杂数据结构与分步变换，转而依赖全局图性质与一阶逻辑。