[論文レビュー] Incomplete Lineage Sorting: Consistent Phylogeny Estimation From Multiple Loci

Elchanan Mossel, Sébastien Roch|arXiv (Cornell University)|Oct 1, 2007

Genomics and Phylogenetic Studies参考文献 14被引用数 6

ひとこと要約

この論文は、不完全な系統選択下で複数の遺伝子系統木から系統樹を推定する統計的に整合的なアルゴリズムであるGLASS（Global LAteSt Split）を紹介する。GLASSは、loci間の最小共析時刻を活用して正しい系統樹トポロジーを推定し、標準的な共析モデル下で一貫性を示し、推定誤差が中程度の範囲にある場合にも耐性を示す。明示的な収束速度が保証されている。

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.

研究の動機と目的

• 不完全な系統選択により遺伝子系統木が真の系統樹と食い違うという、不正確な系統樹推定の課題に対処すること。
• 共析の確率的ばらつきによって遺伝子系統木が不一致になる状況でも、計算的に効率的でありながら統計的に一貫性を保つ手法を開発すること。
• 共析モデル下での系統樹推定に向けた明示的な収束速度を提示し、再構築品質の実用的評価を可能にすること。
• 距離ベースの再構築に一般化可能な、複数lociへのトポロジー整合法の拡張を図り、分子時計の仮定を必要としないこと。

提案手法

• GLASSは、各系統対から得られる個体群間の最小分岐時刻を用いて、個体群間の共析時刻を推定する。
• 分岐長のしきい値に基づく組み合わせ的条件を用いて、これらの最小共析時刻から正しい系統樹トポロジーを同定する。
• 共析モデル下では、2つの系統クラスタ間の最短共析時刻が、系統樹におけるそれらの最も最近い共通祖先に対応することを利用している。
• 共析時刻の推定誤差が制限された範囲内にある場合に一般化され、誤差が最短分岐長に関連するしきい値未満であれば正しさが保たれる。
• アルゴリズムは、分子時計を仮定しない任意の距離ベースの系統樹再構築手法と互換性がある。
• Rosenbergのトポロジー整合法を複数lociに拡張し、遺伝子間で共析時刻情報を集約することで、系統レベルの関係を推定する。

実験結果

リサーチクエスチョン

• RQ1 遺伝子系統木が不完全な系統選択によって不一致になる状況下でも、単純かつ計算的に効率的な手法が、系統樹推定において統計的整合性を達成できるか。
• RQ2 複数lociにわたる最小共析時刻が、真の系統樹トポロジーを正しく同定するための条件は何か。
• RQ3 共析時刻測定における有界な推定誤差が生じる場合、この手法はどのように動作するか。
• RQ4 分子時計の仮定を必要としない距離ベースの手法へ一般化できるか。
• RQ5 共析モデル下での系統樹再構築に向けた明示的な収束速度は、どのように導出できるか。

主な発見

• GLASSは標準的な共析モデル下で統計的に一貫しており、十分に多くの非連鎖lociが与えられれば、高確率で正しい系統樹を再構築する。
• 誤差が最短分岐長の半分未満であれば、共析時刻の推定誤差が中程度の範囲にある場合でも、この手法は一貫性を保つ。
• 明示的な収束速度が導出された：誤った再構築の確率は、loci数とともに指数関数的に減少し、$ \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) $ 個のlociを使用することで、高確率で正しい再構築が達成される。
• GLASSは任意の距離行列手法と組み合わせ可能であり、分子時計の仮定を必要としない。ただし、系統樹の各枝内で変異率が一貫していることが条件である。
• この手法はトポロジー整合法を複数lociに一般化し、遺伝子系統木の不一致からの系統樹推定に向けた堅牢なフレームワークを提供する。