- 著者: Sonya Hessey, Abigail Bunkum, Ariana Huebner, Kerstin Haase, Kristiana Grigoriadis, Cristina Naceur-Lombardelli, Wing Kin Liu, Simone Zaccaria, Nicholas McGranahan, Charles Swanton, Mariam Jamal-Hanjani
- Corresponding author: Simone Zaccaria (UCL); Nicholas McGranahan (UCL); Charles Swanton (Francis Crick Institute); Mariam Jamal-Hanjani (UCL)
- 雑誌: Nature
- 発行年: 2026
- Epub日: 2026-03-17
- Article種別: Original Article
- PMID: 42056508
背景
転移は癌関連死亡の主要な原因であり、肺癌は転移性疾患の最大の割合を占める。癌細胞が原発部位から二次部位へ拡散する生物学的プロセスの理解は、その予防と治療戦略に不可欠である (Pich et al. 2022)。しかし、患者における転移播種のメカニズムに関する包括的かつ縦断的な解析は、臨床的制約によりこれまで困難であった。先行研究では、限られた転移サンプリングのため、転移の不均一性や播種パターンが過小評価されてきた可能性が指摘されている (McPherson et al. 2016; Al Bakir et al. 2023)。特に、非小細胞肺癌 (NSCLC) においては、その高い発生率、de novo転移性疾患としての頻繁な発現、および根治的意図の手術後の高い再発率により、転移性癌症例の大部分を占める (Sung et al. 2021; Ganti et al. 2021)。
癌細胞の移動は一過性の事象であり、患者でリアルタイムに観察することは困難である。しかし、DNAシーケンスデータを用いて原発腫瘍と転移における癌細胞のサブクローン集団の進化史を追跡することで、このプロセスの遡及的な理解が得られる (Reiter et al. 2017; Naxerova 2025)。卵巣癌 (McPherson et al. 2016)、乳癌 (Brown et al. 2017)、前立腺癌 (Gundem et al. 2015) の研究では、原発腫瘍と転移間、および解剖学的に異なる転移部位間で双方向的に移動する複数の腫瘍サブクローンを含む複雑な転移移動パターンが明らかになっている。この転移に関する見解の正確性と完全性は、サンプリングされた転移が患者の疾患負荷をどの程度代表しているか、および比較のための原発腫瘍の利用可能性に依存する (Hessey et al. 2023)。
このような包括的かつ縦断的なサンプリングは、生きた患者では臨床的に実現不可能であるため稀にしか行われないが、TRACERx (Tracking Non-small Cell Lung Cancer Evolution through Therapy (Rx); ClinicalTrials.gov identifier NCT01888601) のような前向き臨床研究とPEACE (Posthumous Evaluation of Advanced Cancer Environment; NCT03004755) のような研究剖検プログラムを統合することで達成可能となる。TRACERxは早期の切除可能な原発NSCLCの多領域プロファイリングを実施し、PEACEは診断から死亡までの完全な疾患経過をカバーする腫瘍組織リソースの構築を可能にする。これらの統合されたリソースは、転移の遺伝的不均一性、播種パターン、および転移能に関わるゲノム的・時間的・解剖学的特性を包括的に解明するための貴重な機会を提供する。しかし、転移の多様性と播種におけるゲノム的、時間的、解剖学的要因の包括的な相互作用については、依然として未解明な点が多く、特に広範な転移性疾患におけるサブクローン進化の全容を捉えるためのデータが不足している。
目的
本研究の目的は、TRACERxおよびPEACEコホートに登録されたNSCLC患者24例から得られた、診断から死亡までの縦断的多領域WESデータを用いて、転移の遺伝的不均一性、播種パターン、および転移能に関わるゲノム的、時間的、解剖学的特性を包括的に解明することである。具体的には、転移内および転移間のサブクローン多様性の程度、変異プロセスの時間的変化、転移固有のドライバー変異の発生、転移から転移への播種 (met-met seeding) の広範性、播種能力と転移のin situ期間との関連、および染色体不安定性 (CIN) が胸腔外播種に与える影響を詳細に解析することを目的とした。
結果
転移の遺伝的不均一性: 個々の転移内では中央値で7個 (範囲3–37) のサブクローンが検出され、79%の転移が他の転移では検出されないサブクローンを含んでいた。患者あたりの転移固有サブクローンの中央値は28個であり、これは先行TRACERx研究 (Al Bakir et al. 2023) の中央値2.5個と比較して11-foldの増加であった。これは、サンプリング数の増加が転移のサブクローン多様性の検出に大きく寄与することを示唆する。同一患者由来の転移同士は、対応する原発腫瘍よりも相互に類似していた (変異多様性 p = 0.004、SCNA多様性 p = 6.0 × 10⁻⁵、Wilcoxon signed-rank test)。しかし、原発腫瘍と転移間の遺伝的乖離は大きく、切除された原発NSCLC腫瘍は、臨床経過中に検出される転移を代表しない可能性が示唆された (Fig. 1c)。
変異プロセスの時間的変化: 喫煙およびAPOBECシグネチャは原発サブクローンで優位であり、転移固有サブクローンでは減少した。特に、喫煙シグネチャSBS4に起因する全変異の83.4%はtruncalで発生しており、NSCLC進化の初期プロセスであることが示された。転移固有サブクローンでは、clock-like (SBS1, SBS5)、APOBEC (SBS2, SBS13)、白金 (SBS31, SBS35)、喫煙 (SBS4) シグネチャの順に多く検出された。白金化学療法を受けた14例中9例 (64%) の転移固有サブクローンで白金関連シグネチャが検出され、治療関連変異誘発の可能性が示唆された。これらの変異プロセスは部位特異的ではなく、複数の転移部位で発生していた。
転移固有ドライバー変異の発生: 同定された196個のドライバー変異のうち、27% (53個) が転移固有であり、83% (20/24) の患者が少なくとも1つの転移固有ドライバー変異を有していた (中央値2個、範囲0–6個)。転移固有ドライバー変異の数は、転移固有変異の総数と相関し (Pearson’s R = 0.8, p = 3.2 × 10⁻⁶)、両者とも化学療法期間と関連していた。白金シグネチャ活性を持つ転移固有サブクローンは、そうでないサブクローンよりも転移固有ドライバー変異を含む可能性が高かった (p = 0.016, chi-squared test)。最も頻繁に変異していた腫瘍抑制遺伝子 (TSG) はTP53であり (75% (18/24) の患者)、89% (16/18) でbiallelic不活性化が認められた。全ゲノム重複 (WGD) は92% (22/24) の患者で発生し、24% (10/41) は転移内で発生した。臨床的にactionableな転移固有変異はSMARCA4 (p.E1083K) の1例のみであった (Fig. 2)。
転移から転移への播種 (met-met seeding) の広範性: 88% (21/24) の患者でmet-met seedingが検出された。サンプリングされた全転移の60% (156/258) が他の転移からの播種であり、38% (98/258) が原発腫瘍からの播種であった。62.5% (15/24) の患者では、複数の原発腫瘍サブクローンが播種し、それぞれが異なる転移を樹立した。しかし、原発から播種したサブクローンのうち、二次的なmet-met播種後代を生じたのは45% (30/67) に過ぎなかった。MACHINAアルゴリズムによる転移移動の確率的評価でも、met-met seedingが優勢であることが示された (Fig. 3c)。
播種能力とin situ期間: Met-met seedingが優勢な患者は、そうでない患者と比較して、より長い無病生存期間 (DFS中央値14ヶ月 vs 5ヶ月、p = 0.036) と全生存期間 (OS中央値32ヶ月 vs 15ヶ月、p = 0.019) を示した。初回再発スキャンで検出された転移の36.8% (14/38) が二次転移を播種したのに対し、最終スキャン後 (剖検時) に検出された転移では16.4% (10/61) であった (p = 0.03, Fisher’s exact test)。播種した転移は、画像上の最大体積が有意に大きく (p = 0.025)、サブクローン数が多く (p = 0.028)、分岐系統樹を持つ傾向が強かった (p = 9.8 × 10⁻⁸)。これらの結果は、転移がin situで存在する期間が長いほど、met-met seedingの可能性が高まることを示唆する (Fig. 4d)。
解剖学的制約と染色体不安定性 (CIN): Met-met seedingの71.3% (92/129) は同一体腔内に留まった (p = 2.40 × 10⁻⁵, Fisher’s exact test)。胸腔内転移は胸腔外転移よりも高頻度に二次播種した (26.4% (32/121) vs 13.9% (19/137), p = 0.013)。原発からの播種サブクローンは非播種サブクローンよりもSCNA負荷が有意に高く (p = 9.6 × 10⁻⁵, LME model)、met-met播種サブクローンでも同様であった (p = 3.0 × 10⁻⁷, LME model)。特に、胸腔外播種サブクローンのSCNA負荷は胸腔内播種サブクローンよりも有意に高く (p = 7.2 × 10⁻⁴, LME model)、CINが胸腔外播種に特異的に関与する可能性が示唆された (Fig. 5e,f)。
考察/結論
本研究は、TRACERxおよびPEACEコホートのNSCLC患者24例における転移進化の最も包括的な記述であり、従来の限定的サンプリングでは過小評価されていた転移の遺伝的多様性と播種パターンの複雑さを明らかにした。
先行研究との違い: これまでの研究では単一の原発播種サブクローンがほとんどの転移を形成するとされてきた (Al Bakir et al. 2023) のに対し、本研究では62.5%の患者で複数の原発サブクローンが播種し、それぞれが異なる転移を樹立したことを示した点でこれまでと異なり、転移播種の多様性を強調する。また、サンプリングされた転移の60%が他の転移から播種されるというmet-met seedingの広範性は、先行研究で報告されたよりも転移間播種の寄与が大きいことを示唆する。
新規性: 本研究で初めて、転移のin situ期間が二次転移播種能力に影響を与えることを明らかにした。早期に検出された転移ほど二次播種能力が高いという知見は、転移の進化における時間的要素の重要性を新規に提示する。さらに、染色体不安定性 (CIN) が胸腔外播種と特異的に関連するという発見は、循環播種や非胸腔微小環境への適応にCINが寄与する可能性を示唆し、この経路の転移プロセスを促進するゲノム的メカニズムを本研究で初めて特定した。
臨床応用: Met-met seedingの広範性は、既存転移の局所制御が転移進行予防に有効である可能性を示唆し、局所固形療法 (LCT) の生物学的根拠を提供する。Met-met播種にはin situ期間が必要であることから、治療介入の時間的窓が存在する。CINが高い転移が胸腔外播種能力を持つという知見は、LCTの対象患者を適切に選択する上で、特に高CINの転移を標的とすることの臨床的意義を示唆する。しかし、臨床的にactionableな転移固有ドライバー変異は稀であったことから、単一転移生検に基づく治療選択の限界が改めて示され、転移のゲノムプロファイリングの重要性が強調される。
残された課題: 本研究の限界として、88%の患者が全身療法を受けていたため、未治療患者やde novo転移性疾患患者への結果の一般化可能性は不明である。また、コホートサイズが小さいため、治療が転移進化やゲノムイベントに与える影響を詳細に評価する能力には限界があった。バルクWESは高深度で実施されたものの、全ゲノムシーケンスやシングルセルシーケンスと比較してサブクローン多様性を過小評価する可能性があり、構造変異、染色体外DNA、腫瘍微小環境、その他の非遺伝的プロセスの役割については調査できなかった。これらの残された課題は、今後のTRACERx EVO研究で高深度全ゲノムシーケンスを用いて対処される予定である。
方法
TRACERxおよびPEACEに共登録されたNSCLC患者24例から、原発腫瘍108領域、生前転移41領域、剖検転移352領域の計501検体に対し、高深度WES (中央値深度401.2×、四分位範囲 360.2–441.5) を実施した。このサンプリングは、死亡前に放射線学的に検出された転移の70% (112/160) をカバーした。
サブクローン構造と系統樹の再構築: 体細胞変異、体細胞コピー数異常 (SCNA)、および全ゲノム重複 (WGD) の詳細なゲノム解析を実施し、各患者の疾患のサブクローンアーキテクチャと系統発生史を再構築した。CONIPHER (Leshchiner et al. 2023) を用いて体細胞変異クラスターを同定し、腫瘍系統樹を構築した。変異クラスターは、転移播種に対する発生時期に基づいて、truncal、primary-unique、metastasis-unique、shared subclonalの4つのグループに分類された。ALPACA (Pawlik et al. 2025) を用いてサブクローン特異的なコピー数プロファイルを推測し、ParallelGDDetect (Frankell et al. 2023) を用いてWGDイベントを同定した。
転移播種経路の推定: MACHINAアルゴリズム (El-Kebir et al. 2018) を各患者の系統樹に適用し、転移を播種したサブクローンとその対応する移動経路を推定した。MACHINAは、系統樹のトポロジーと各腫瘍領域におけるサブクローンの存在に基づいて、最もパースィモニウスな転移経路を推測する。転移間播種の信頼性を評価するため、各転移移動に確率的アプローチを開発し、LOHイベントの保存性を定量化することで移動方向性を検証した。
変異シグネチャ解析: CONIPHERで定義された変異クラスターに対し、SigProfilerAssignment (Alexandrov et al. 2020) を用いてCOSMIC変異シグネチャ活性を推定した。シグネチャは、clock-like (SBS1, SBS5)、喫煙 (SBS4)、APOBEC (SBS2, SBS13)、白金 (SBS31, SBS35) など、NSCLCで既知の病因学に基づいて分類された。
ドライバー変異のタイミング解析: OncoKB (Chakravarty et al. 2017)、openCRAVAT、Ensembl Variant Effect Predictorを用いて体細胞変異をアノテーションし、既知の腫瘍形成ドライバーを同定した。TSGのbiallelic不活性化は、変異とLOHの組み合わせによって評価された。
SCNA負荷と播種能力の関連評価: サブクローンごとのSCNA負荷は、検出されたブレークポイントの総数として計算された。SCNA負荷と播種能力の関連を評価するため、線形混合効果モデル (LME model) を用いて、播種サブクローンと非播種サブクローンのSCNA負荷を比較した。
放射線学的データとの統合: 23例の患者から得られた匿名化された臨床画像スキャンとレポート (CT, PET-CT, MRI) を用いて、転移の初回検出時期と最大体積を評価した。放射線学的に可視で測定可能な転移は、ITK-SNAP v.4.2.0を用いて輪郭が描かれ、3次元腫瘍体積が算出された。統計解析には、Wilcoxon signed-rank test、Mann-Whitney U test、Fisher’s exact test、Pearson’s R、chi-squared testが用いられた。