Revisiting Xia' s Index of Substitution Saturation from a Mathematical Perspective

技術的要約：数学的観点からの Xia の置換飽和指数（Index of Substitution Saturation）の再検討

問題提起
Xiaら [2003] によって導入された置換飽和指数 $I_{ss} = \bar{H}/H_{FSS}$ は、系統樹構築前に配列アライメントが系統学的シグナルを保持しているかどうかを判断するために広く用いられている。この指標は、サイト組成の平均エントロピー（$\bar{H}$）を完全飽和下で期待されるエントロピー（$H_{FSS}$）と比較するものである。この統計量は計算効率が高い（樹形やモデルの適合を必要としない）一方で、その期待値の理論的特性化や、臨界閾値（$I_{ss.c}$）の妥当性は、第一原理から厳密に導出されていない。本論文は、この統計量の期待値の数学的正当性、それがより単純な指標へと崩壊する条件、およびレート不均一性と配列長に対する臨界値の堅牢性について取り組むものである。

手法
本研究では、スペクトル分解理論、厳密な確率論的導出、および広範なシミュレーションを組み合わせて用いる。

1. スペクトル分解: 著者は、対称化された生成行列 $S = \Pi^{1/2}Q\Pi^{-1/2}$ の固有基底上で遷移行列 $P(t)$ を分解する。これにより、$E[I_{ss}]$ の期待値を、非ゼロのレート行列 $Q$ の固有値 $\mu$ に対応するスペクトル座標 $\xi_\mu = \langle e^{\mu d} \rangle$ の線形汎関数として表現することが可能となる。
2. 厳密な閉形式: 本論文は、スターツリーにおける $E[I_{ss}]$ の厳密な閉形式（定理7）、および一般の樹形における二次の漸近展開（定理5）を導出する。これらの導出において、2つの正規化慣例を区別している。すなわち、XiaらおよびDAMBEで使用されている、経験的な組成による正規化を行う「慣例(A)」と、プラグイン・エントロピー・アプローチである「慣例(B)」である。
3. シミュレーション設計: シミュレーションには IQ-TREE 内の AliSim を使用した。重要な設計上の特徴は、統計量が感度を失う飽和プラトーを回避するため、指数がその範囲内で一様に分布するようにする一様 $I_{ss}$ サンプリング戦略（式21）を用いていることである。また、本研究では「フェルゼンスタイン–ファリス転送」（式22）を利用して、指数を一定に保ちながら樹形トポロジーと枝長を変化させ、統計量の系統学的シグナルに対する感度を分離して検証している。
4. 経験的検証: 著者は、脊椎動物のミトコンドリア COI データセット（Xia and Lemey, 2009）およびゲノムスケールの鳥類アライメント（Jarvis et al., 2014）を再解析し、実データに対する指数の挙動をテストした。

主要な貢献と結果

1. 統計量の特性化と「崩壊」条件
本論文は、$E[I_{ss}]$ がスペクトル座標のベクトル $\xi$ を通じてのみ樹形に依存することを証明している。決定的なことに、代数的重複度（$\tau_\mu$）と組成ウェイト（$w_\mu$）の比が、すべての異なる非ゼロ固有値にわたって一定である場合に限り、指数は平均ペアワイズ $p$-距離（補正前の相違サイトの割合）へと崩壊する（定理13）。

• 一様組成: 一様な塩基組成を持つあらゆる可逆モデルにおいて、この条件は正確に成立する。
• 退化スペクトル: 単一の異なる非ゼロ固有値を持つモデル（例：JC69, F81）では、組成に関わらずこの条件は成立する。
• 失敗レジーム: 非一様な組成と非退化スペクトルを持つモデル（例：偏った頻度を持つ HKY や GTR）では、この条件は成立しない。これは、元の臨界値を較正するために使用されたまさにそのレジームである。このレジームにおいて、$I_{ss}$ は置換の実現数ではなく、モデルの組成に関する情報を追加してしまう。

2. 臨界値とレート不均一性
本研究は、表形式の臨界値（$I_{ss.c}$）がレート不均一性に対して誤って較正されていることを示している。

• 非識別性: 系統学的回復が失敗する臨界樹長（$T_{Lc}$）は、レート不均一性パラメータ $\alpha$ の厳密な単調関数である（命題11）。しかし、指数 $I_{ss}$ 自体は、近飽和領域において $\alpha$ に対してほぼ不変である。
• 解像度の崩壊: 指数の天井付近において、$I_{ss}$ の樹長変化に対する感度は消失する（式17）。樹長が15倍変化しても（回復が失敗する範囲を跨いでいるにもかかわらず）、指数の値の変化は約1.2%にとどまる。したがって、臨界値は、テストが最も必要とされる領域において、中程度の飽和と深刻な飽和を区別することができない。

3. 配列長への依存性
本論文は、指数の標準誤差が $L^{-1/2}$ に比例してスケールすることを確立している（定理19）。臨界値はシミュレーションから導出された固定の閾値であるため、テストの統計的検出力は配列長とともに増大する。短いアライメントにとって有意なマージンであっても、長いアライメントにとっては無視できる可能性がある。つまり、閾値は配列長と共同で解釈される場合にのみ意味を持つ。

4. 予測能と情報量
条件付き独立性テスト（セクション 2.13）を通じて、著者は $I_{ss}$ が $p$-距離と高い相関関係にある一方で、「遷移帯」（$p$-距離が 0.63 から 0.70 の間）においては $p$-距離単独よりも系統学的回復をより良く予測できることを示している。しかし、この予測能は「置換のカウント」を測定することから生じるのではない。むしろ、$I_{ss}$ は樹形の形状（特に内部枝の割合）と高次のサイトパターンの組成を追跡している。この指数は「情報フリー」ではないが、それが提供する情報は、その名称の由来である置換数とは異なるものである。

5. 経験的知見

• 脊椎動物 COI: 第3コドンポジションの分割（$I_{ss} \approx 0.715$）は、報告されている2つの臨界値（対称的 vs 非対称的トポロジー）の間に位置している。指数は天井に達しており、解像度の崩壊によりテストは決定不能となっている。
• 鳥類ゲノム: ゲノムスケールのアライメント（$I_{ss} \approx 0.119$）は、意図通り閾値を十分に下回っている。
• 他の指標との比較: 本論文は $I_{ss}$ を $\lambda$ 指数（Jermiin & Misof, 2020）および $C$ 値（Struck et al., 2008）と比較している。その結果、平均 $\lambda$ は組成によって再スケールされた平均 $p$-距離と数学的に等価であり、同様の限界を継承していることが判明した。一方、遷移/トランスバージョン比に基づく $C$ 値は、$I_{ss}$ を追跡せず、飽和した分割を正しく特定できない。

意義と主張
本論文は、置換飽和指数は特定のレジーム（一様組成または退化スペクトル）においては飽和を検出するための有効な統計量であるが、一般的なケースにおいては根本的な限界がある、と結論付けている。

• テストに合格すること（低い $I_{ss}$）は、飽和に対する証拠（反証）となる。
• テストに失敗すること（高い $I_{ss}$）は、飽和の弱い証拠であり、必ずしもデータが解決不能であることを意味しない。それは単に、指数が分岐に対する感度を失う天井に達していることを示している可能性がある。
• 臨界値は普遍的な定数ではなく、配列長、レート不均一性、および較正に使用される特定の推論エンジンに依存する。
• 著者は、$I_{ss}$ は置換数の直接的な尺度ではなく、むしろレート行列のスペクトル特性とアライメントの組成の関数であると断言している。したがって、これらへの依存性を理解せずに「系統学的シグナル」の代用として用いることは、特に組成が非一様であったりレート不均一性があったりする場合、誤解を招く可能性がある。

本研究は新しい指数を提案するものではなく、既存の指数の解釈のための数学的枠組みを提供し、ユーザーに対して、閾値は自身の配列長を伴う「パワー閾値」であること、そして飽和付近での統計量の挙動は「解像度の崩壊」現象によって支配されていることを警告するものである。