技術要約：指数関数的に増殖する集団における選択を伴うサイト頻度スペクトル

問題提起

本研究は、**ドライバー変異（driver mutations）**によって駆動される指数関数的増殖下にある集団における、サイト頻度スペクトル（SFS）、特に変異アレル頻度の分布について調査するものである。中立進化（変異が適応度上の利点を持たない場合）におけるSFSは、期待値として典型的には $1/j^2$ の法則に従うことがよく知られているが、**選択（selection）**下でのSFSの挙動については、まだ十分に理解されていない。

著者らは、以下の条件下での**二類型連続時間線形出生死滅過程（two-type continuous-time linear birth-death process）**に焦点を当てている：

• **タイプ0（野生型）**細胞は、正の純マルサス成長率 $\lambda_0 > 0$ で増殖する。
• **タイプ1（変異型／ドライバー）**細胞は、タイプ0細胞から変異によって発生し、より高い純成長率 $\lambda_1 = \lambda_0 + b$ を持つ。ここで $b > 0$ は選択的利点を表す。
• 本モデルは、複製依存的な変異（細胞分裂時に発生）と、複製に依存しない変異（一定の割合で発生）の両方を許容する。

主な目的は、ドライバーSFS $S_j(t)$ を特徴付けることである。これは、時刻 $t$ においてサイズ $j$ を持つ個別のドライバークローンの数をカウントするものである。本研究は、中立進化や半決定論的近似に関する先行研究を拡張し、正確なモーメント、漸近挙動（時間および頻度における）、およびほとんど確実な収束（almost sure convergence）の結果を導出することを目指している。

手法

著者らは、**分岐過程理論（branching process theory）**に基づいた厳密な確率論的枠組みを用いている。

1. モデルの定式化:

   • システムは、タイプ0細胞がタイプ1クローンを播種する、可約分岐過程（reducible branching process）としてモデル化される。
   • 各タイプ1クローンは、出生率 $a+b$ および死亡率 $d$ を持つ独立した線形出生死滅過程として進化する。
   • 解析では、決定論的適応度（固定された $b$）と、ランダム適応度（各クローンの適応度利点 $B$ が分布 $Q$ から抽出される）の2つのシナリオを検討する。

2. 解析手法:

   • 正確なモーメント: 著者らは、クローンの開始時刻とタイプ0集団の軌跡に条件付けることで、 $S_j(t)$ の一次および二次モーメントの積分表示を導出する。
   • 漸近解析: 大時間極限（$t \to \infty$）および大頻度極限（$j \to \infty$）は、分岐過程の極限となる確率変数の性質（例：$Y = \lim e^{-\lambda_0 t} Z_0(t)$ および $Y^{(i)}_1 = \lim e^{-\lambda_1 t} Z^{(i)}_1(t)$）を用いて解析される。
   • ほとんど確実な収束: 強法則の数（strong laws of large numbers）を証明するために、著者らは $L^2$ 近似を構築する。彼らはSFSを単調プロセス（サイズ $j$ への進入と退出をカウントするもの）へと分解し、マルチンゲール引数とボレル・カンテリの補題を利用して収束を確立する。
   • 固定検出サイズ: 解析は、ランダムな検出時刻 $\tau_n$（総集団サイズが $n$ に達した時）へと拡張される。これは、超臨界増殖における時間とサイズの間の関係（$\tau_n \sim \frac{1}{\lambda_1} \ln n$）を利用している。

主要な貢献と結果

1. 正確なモーメントと平均挙動

本論文は、ドライバーSFSの平均および二次モーメントの正確な積分公式を提供している。

• 選択の効果: 著者らは、選択強度 $b$ を大きくすると、平均SFSにおいて小さなクローン（$j$ が1に近いもの）の頻度が減少することを証明している。これは、強い選択がクローンを小さなサイズからより急速に通過させるため、小さな頻度のビンに滞在する時間が短くなるからである。
• 漸近的冪乗則: 大きな $j$ および大きな $t$ に対して、平均SFSは以下の冪乗則に従う：
  $$\lim_{t \to \infty} e^{-\lambda_0 t} E[S_j(t)] \sim C \cdot j^{-\lambda_0/\lambda_1 - 1}$$
  ランダム適応度モデルでは、テール挙動は最大適応度利点 $b_{max}$ を持つクローンによって支配され、対数補正を伴う冪乗則を示す。

2. ほとんど確実な収束（強法則）

ドライバーSFSに関する強法則の数の証明は、主要な貢献である。

• 固定時間: 非絶滅を条件として、$e^{-\lambda_0 t} S_j(t)$ は、極限となるタイプ0集団のサイズ $Y$ に比例する確率変数にほとんど確実に収束する。
• 固定サイズ: 同様に、$n^{-\lambda_0/\lambda_1} S_j(\tau_n)$ は、$Y$ と極限となるタイプ1集団のサイズ $W$ を含む確率変数にほとんど確実に収束する。
• 正規化されたSFS: 全変異数に対する頻度 $j$ の変異の割合は、成長率の比 $\lambda_0/\lambda_1$ と絶滅確率 $p_1$ のみに依存する決定論的な極限に収束する。これは、SFSの形状を用いて進化パラメータを推定できることを示唆している。

3. 希薄性とカットオフ周波数

著者らは、頻度 $j$ が時間とともに変化する場合、特に $j(t)$ の場合の $S_j(t)$ の挙動を調査している。彼らは、以下のスケールを持つ決定的なカットオフ周波数を特定している：
$$j(t) \sim \exp\left( \frac{\lambda_0 \lambda_1}{\lambda_0 + \lambda_1} t \right)$$

• 三分則（Trichotomy）:
  • $j(t)$ がこの速度よりも速く成長する場合、そのサイズのクローンの数は0に近づく。
  • $j(t)$ がこの速度で成長する場合、クローンの数は（非絶滅を条件として）ポアソン確率変数に分布収束する。
  • $j(t)$ がより遅く成長する場合、クローンの数は無限大に発散する。
    この閾値は、新しいクローンの生成と、それらがサイズ $j(t)$ まで成長することのバランスを表している。

4. 相対的SFSと大きなクローン

本論文は、全変異集団の一定の割合 $> f$ を占めるクローンをカウントする相対的SFS $R_f(t)$ を分析している。

• 小頻度極限: $f \to 0$ のとき、相対的SFSは $f^{-\lambda_0/\lambda_1}$ という冪乗則に従う。
• 半決定論的モデルとの接続: 微小変異極限において、著者らは先行文献（DurrettおよびTungなど）で使用されている半決定論的近似を厳密に正当化している。彼らは、確率的モデルが、野生型集団が決定論的に成長する半決定論的モデルに収束することを示している。
• 最大のクローン: 総集団に対する最大変異クローンのサイズ（正規化されたもの）は、小変異極限において、ポアソン・ディリクレ分布（具体的には $\alpha = \lambda_0/\lambda_1$ である $PD(\alpha, 0)$）の上限に収束する。

5. 固定検出サイズにおけるテール挙動

固定検出サイズ $n$ において、SFSは $j=n$ における点質量（サイズ $n$ に達した単一の最大のクローンを表す）を示す。この事象の確率は、最初の変異が発生した時点における野生型集団と変異型集団の間の競争に依存する。

重要性と主張

本論文は、選択を伴う確率的かつ指数関数的に増殖する集団における、ドライバーSFSの包括的な理論的特性付けを提供すると主張している。

• パラメータ推定: 著者らは、導出された極限定理により、進化パラメータ、特に変異細胞の適応度増加（野生型に対する $b$）の推定が可能であると断言している。正規化されたSFSは $\lambda_0/\lambda_1$ の決定論的関数に収束するため、観測されたSFSを反転させることで、理論的に選択強度を推定することができる。
• 厳密な正当化: 本研究は、以前に使用されていた半決定論的モデルに対して厳密な確率論的基礎を提供し、小変異極限におけるそれらの妥当性を確認している。
• 中立進化との区別: 結果は、選択がいかに中立モデルと比較してSFSの形状を根本的に変えるか（特に、小さな頻度のクローンを抑制し、成長率の比によって決定される特定の冪乗則のテールを作り出すこと）を強調している。

著者らは、現在の研究が単一のドライバー変異とドライバーSFSに焦ellしていることに注意している。彼らは、今後の研究において、多類型集団における「ヒッチハイカー（随伴）」変異や連続的なドライバー獲得を含む解析へと拡張することを認めている。