DTI-based mesoscale connectomics detects survival-relevant white matter substructures in glioblastoma patients

技術要約：生存解析におけるDTIベースのメソスケール・コネクトミクス

問題提起
膠芽腫（GB）は、腫瘍細胞が白質（WM）路に沿って移動するという極めて浸潤性の高い性質を特徴としており、これが完全な外科的切除を困難にし、予後不良（全生存期間［OS］の中央値は12〜18ヶ月）を招いている。術後切除範囲（EOR）、年齢、分子プロファイルなどは既知の予後指標であるが、個別の線維サブコンポーネント内における微細構造の完全性が生存に与える具体的な影響については、依然として不明である。従来の画像アプローチは、白質路全体の平均値や手動による関心領域（ROI）解析に依存することが多いが、これらは腫瘍の浸潤や患者の予後に関連する、微細構造の完全性における微妙かつ空間的に限定された変化を検出できない可能性がある。

方法論
本レトロスペクティブ研究では、UPenn-GBMコホート（The Cancer Imaging Archive）から得られた498名のIDH野生型膠芽腫患者の拡散テンソル画像（DTI）および臨床データを分析した。本研究では、OSと相関する白質路セグメントを特定するために、**コネクトメトリー（connectometry）と相関トラクトグラフィー（correlational tractography）**を組み合わせたメソスケール・コネクトミクス・フレームワークを採用した。

1. データ前処理： 術前のDTIデータ（30拡散方向、b=1000 s/mm²）をDSI Studioを用いて前処理した。データは、スピン分布関数（SDF）を生成するために、q空間微分同型再構成（QSDR）を用いたMNI空間への空間正規化を行った。
2. コネクトメトリー： HCP-842アトラスに沿った線維方向にSDFをサンプリングすることで、局所的なコネクトームベクトルを推定した。これにより、隣接するボクセルにおける拡散スピンの密度を表す局所コネクトーム行列が作成され、腫瘍組織を明示的に除外することなく（腫瘍コアの固有の低異方性に依存）、メソスケールな結合性を捉えた。
3. 相関トラクトグラフィー： 共変量（年齢、性別、EOR）を制御した上で、局所コネクトーム指標（特に異方性比［FA］）とOSとの間の偏相関を計算した。統計的有意性は置換テストによって評価された。
4. 線維追跡： 統計マップ（t値閾値 4.2–5.2）から決定論的線維追跡を行い、OSと有意に関連する連続的な線維セグメントを生成した。線維は、トポロジーに基づく枝刈り、20ボクセルの長さ閾値、および偽発見率（FDR）閾値0.01を用いてフィルタリングされた。
5. 統計解析： 再構成された線維はサブコンポーネントに分割された。平均FA値は、健康な対照群および最大選択ランク統計から導出された閾値に基づき、二値化（正常 vs 低下）された。生存差は、Kaplan-Meier曲線（ログランク検定）およびCox比例ハザードモデルを用いて分析された。線形および指数回帰モデルを用いて、FAとOSの相関を評価し、モデルの適合度は赤池情報量基準（AIC）を通じて評価された。

主要な貢献

• メソスケール解像度： 本研究は、メソスケール・コネクトミクスが、生存に関連する特定の白質サブコンポーネント（白質路全体の平均ではなく）を特定できることを示している。
• 生存に関連する線維路の特定： 本研究は、OSと相関する主要な線維路として、主に**視床放射（投影線維）および脳梁（交連線維）**内の特定の線維クラスターをマッピングした。
• 葉特異的なパターン： 本研究は、相関する主要な線維タイプが腫瘍の位置によって異なることを明らかにした。これは規範的な解剖学的構造を反映しており、中心部腫瘍では投影線維、前頭部腫瘍では交連線維、側頭部腫瘍では連合線維がそれぞれ関連している。
• 定量的予後指標： 本研究は、これらの特定のサブコンポーネントにおける正常なFA値が、標準的な臨床因子を制御した後でも、独立して有意に長い生存期間と関連していることを確立した。

結果

• 生存優位性： 特定された線維サブコンポーネントにおいて正常なFA値を持つ患者は、低下したFA値を持つ患者と比較して、OSの中央値が有意に長かった。
  • 全脳： 487日 vs 207日 (p < 0.001)。
  • 左半球： 488日 vs 241日。
  • 右半球： 493日 vs 216日。
  • 1年生存率は、正常FA群の方が2倍以上高かった（全脳で68.7% vs 32.8%）。
• ハザード比： Coxモデルは、FAが独立した生存予測因子であることを確認した。全脳において、FAの増加（ΔFA = 0.1）は、ハザード比（HR）0.555（95% CI: 0.467–0.658）に関連しており、死亡リスクの減少を示した。左半球（HR = 0.589）および右半球（HR = 0.589）の腫瘍についても同様の独立した関連が見られた。
• モデル適合度： 指数回帰モデルは、線形モデルよりもデータへの適合性が高かった（ほとんどのグループでAICの差 > 10）。これは、FAの悪化がある臨界レベルを超えると急速な臨床的低下を招くという、非線形的で閾値に基づいた関係を示唆している。
• 葉特異的な有意性： FAに基づく統計的に有意な生存差は、左前頭葉、右中央領域、および右後頭葉・頭頂葉の腫瘍において観察された。側頭葉では、Coxモデルにおいて有意な関連は見られなかった（p = 0.151）。

意義と主張
本論文は、コネクトメトリーと相関トラクトグラフィーが、膠芽腫患者の生存と密接に関連する微細構造の完全性を持つ白質サブコンポーネントを特定することに成功したと主張している。著者らは、これらの知見が、将来の予後モデリングのためのバイオマーカーとしての価値を浮き彫りにしていると考えている。

本研究は、膠芽腫をネットワークレベルの疾患として理解すべきというパラダイムを支持している。すなわち、アウトカムは局所的な腫瘍負荷だけでなく、重要な経路（特に投影線維および交連線維）における構造的結合性の保持にも影響される。著者らは、術前の評価に線維特異的なFA指標を組み込むことで、患者のリスク層別化や治療の優先順位付けを容易にできる可能性を示唆している。

しかし、著者らは臨床実装に関して慎重なトーンを維持している。彼らは、これらのカットオフ値は探索的であり、仮説生成的なものであり、独立したコホートでの検証が必要であると明示している。また、レトロスペクティブな設計、分子データ（MGMTステータスなど）の欠如、治療レジメンの詳細の欠如、および交差線維の解像度を制限するシングルシェルDTIの使用といった限界も認めている。したがって、本論文は、メソスケール・コネクトミクスが膠芽腫の生物学を理解するための強力なツールである一方で、臨床実装の前に、調和されたプロトコルと縦断的デザインを用いた前向きな多施設共同検証が不可欠であると結論付けている。