A Lauren-associated mitotic programme in gastric cancer: integrative bulk, single-cell and spatial analyses

技術的要約：胃癌におけるローレン分類に関連する有糸分裂プログラム

問題提起
胃癌（GC）は顕著な生物学的異質性を示し、歴史的にローレン分類による腸型および弥漫型癌へと分類されてきた。これらのサブタイプは組織学的構造や臨床転帰において異なるものの、再現可能なサブタイプ特異的な有糸分裂転写プログラムは完全には解明されていない。既存の研究は、統計的な関連性と因果メカニズムを混同しているか、あるいはサブタイプ特有のシグナルと一般的な増殖または細胞組成の変化を区別できていないことが多い。さらに、単一のデータセットから導出された候補遺伝子リストは、モデル選択バイアスやネットワークの不確実性に脆弱である場合が多い。本研究は、各データレイヤー（バルク、シングルセル、空間、および実験）の推論限界を厳格に維持し、因果関係の主張を避けつつ、ローレン分類に関連する堅牢で再現可能な有糸分裂プログラムを定義することを目的とする。

方法論
本研究では、バルク転写物、ネットワーク解析、シングルセル解像度、空間トランスクリプトミクス、および予備的なウェットラボ検証を組み合わせた、統合的かつ多層的なワークフローを採用した。

1. 発見および発現変動解析：

   • コホート： GSE62254（ローレン分類が付与された266検体：腸型137、弥漫型129）。
   • モデル： 年齢と性別で調整した一次線形モデル（~ Lauren + age + sex）およびフル臨床感度モデル（~ Lauren + age + sex + T + N + M）。
   • 基準： FDR < 0.05 かつ |標準化効果量| ≥ 0.30 を両方のモデルで満たし、方向性が一致した遺伝子を選定し、2,978個の堅牢な遺伝子を得た。

2. ネットワーク構築と優先順位付け：

   • WGCNA： 共変量（年齢、性別、T、N、Mの影響を除去し、ローレンの変動は保持した）を考慮した重み付き遺伝子共発現ネットワーク解析を実施した。4つのゲート（ローレンとの関連性、有糸分裂プログラムとの関連性、有糸分裂リファレンスの濃縮、およびスプリットハーフの安定性）を通過した唯一のモジュールとして、449遺伝子の「イエロー」モジュールを特定した。
   • 候補の交差： 2,978個の堅牢なDEGs、449個のイエローモジュール遺伝子、および153個のReactome有糸分裂遺伝子の交差から、39個の候補を導出した。
   • PPIおよび中心性： STRING物理相互作用ネットワーク（スコア ≥ 400）を構築した。次数（Degree）、最大クリーク中心性（MCC）、および媒介中心性（Betweenness）から得られた上位15遺伝子を交差させた。
   • 安定性テスト： 2,000回の固定シード・エッジ削除摂動（エッジの10%を削除）により、優先順位付けされた候補の安定性を確認した（非候補の安定性が0であるのに対し、候補のメディアン安定性は0.981であった）。

3. クロスコホート検証：

   • コホート： GSE15459 および GSE113255（合計 n=570）。
   • メタ解析： ランダム効果メタ解析（REML）を用いて、Hedges' g（弥漫型 － 腸型）を算出した。
   • ロック基準： 最終的な「ハブセット」は、方向の一致、BH補正 P < 0.05、I² ≤ 50%、|プールされたHedges' g| ≥ 0.30、およびleave-one-cohort-out解析における堅牢性に基づき、凍結された。

4. マルチオミクスおよび実験的裏付け：

   • シングルセル（GSE183904）： 患者ごとの擬似バルク集計を用い、特定の細胞系譜へのプログラムの局在化を行った（個々の細胞を独立した反復試行として扱うことを回避）。
   • 免疫コンテキスト： デコンボリューション（quanTIseq, MCP-counter, ESTIMATE）により、免疫共変動を評価した。
   • 空間（GSE270678）： ドナーレベルの解析（ドナー内のセクションをネスト）を行い、ANGPTL4関連のストローマシグナルとの空間的結合を評価した。
   • RT-qPCR： 10例の新鮮な胃癌検体（腸型5、弥漫型5）を用いた独立した検証を実施した。
   • 探索的解析： 仮説生成のために、Connectivity Map (CMap) スクリーニングおよび分子ドッキング（CDC20-Apcinポケット）を実施した。CellChat/NicheNet解析は、偽造テスト（falsification testing）に従事した。

主な結果

• 5遺伝子ハブセット： 本研究は、再現可能な5遺伝子プログラムを特定し、検証した：BIRC5, CDC20, CENPA, CKS1B, NUF2。
• 発現パターン： これら5つの遺伝子はすべて、弥漫型腫瘍と比較して腸型腫瘍において有意に高い発現を示した。プログラムのプールされたHedges' g は −0.827（95% CI −1.007 to −0.647; I² = 0）であり、強力かつ一貫した関連を示した。
• 堅牢性： この関連性は、一般的な増殖マーカーであるMKI67および臨床変数による調整後も方向性が一貫しており、このプログラムが一般的な細胞周期活性とは異なる、特定の有糸分裂組織化を捉えていることを示唆している。
• 細胞局在： シングルセル擬似バルク解析により、このプログラムは主に上皮およびT細胞コンパートメントに局在していることが判明した。ただし、当該データセットにはローレンのラベルが欠如していたため、シングルセルレベルでのサブタイプ特異性を確認することはできなかった。
• 免疫コンテキスト： 弥漫型腫瘍は、推定されるM2マクロファージおよびB細胞の豊富さが高い。この有糸分裂プログラムは、これら弥漫型に富む免疫特徴と負の共変動を示した。
• 空間的結合： 空間トランスクリプトミクスにおいて、このプログラムは、バルクの結果とは方向が逆転したストローマシグナルとのドナーレベルの空間的結合を示した。具体的には、平均ドナーレベルの空間的結合は、腸型ドナー（−0.0216）よりも弥漫型ドナー（0.0349）の方が高かった（差は0.0565）。これは、バルクの有糸分裂プログラムは腸型腫瘍で濃縮されている一方で、このプログラムの局所的なストローマシグナルとの空間的結合は弥漫型においてより強いことを示しており、単純なバルクシグナルの複製ではなく、文脈依存的な組織構成を浮き彫りにしている。
• 実験的検証： 10例の独立した検体を用いたRT-qPCRにより、方向性が確認された（すべての5遺伝子が腸型腫瘍において高い発現を示し、群間の完全な分離が見られた、P = 0.00794）。
• 陰性/無効の結果：
  • CellChat/NicheNet通信解析は、多重検定補正後、安定したサブタイプ特異的な通信軸を特定できなかった。
  • 分子ドッキングは、Tyrphostin 9、palbociclib、およびethacrynic acidのCDC20ポケットへの適合性を示唆したが、これらはあくまで計算上の仮説として提示されており、生化学的な検証は行われていない。

意義および主張
本論文は、胃癌のローレン分類に関連する、再現可能なクロスコホート有糸分裂転写プログラムを定義することを主張している。その意義は、以下の厳格な分析フレームワークにある：

1. 関連性と因果関係の分離： 著者らは、自らの知見が因果メカニズムや治療効果ではなく、あくまで関連性の枠組みを示すものであることを明示している。
2. 堅牢性の確保： 一次モデルおよびフル臨床モデル間での一致、WGCNAの安定性、PPIトポロジー、および3つの独立した検証コホートを要求することで、単一データセットのバイオインフォマティクスに共通する偽陽性のリスクを最小限に抑えている。
3. 推論境界の尊重： バルク組織シグナル、シングルセル局在（擬似バルクを使用）、およびドナーレベルの空間的結合を注意深く区別し、これらのレイヤーが相補的ではあるが等価ではない証拠を提供することを認めている。
4. 検証可能な枠組みの提供： 5遺伝子セット（BIRC5, CDC20, CENPA, CKS1B, NUF2）は、ローレン分類の生物学における機能的分離のための、優先順位付けされた、実験的にテスト可能な標的を提供する。

著者らは、これらの遺伝子が腸型胃癌において高い活性を持つ堅牢な有糸分裂プログラムを定義しているものの、本研究の知見は、現時点では検証された臨床ツール、因果メカニズム、または治療戦略を確立するものではないと結論付けている。本研究は、限定された、仮説生成のための基礎となるものである。