- 著者: Oriol Pich, Sophia Ward, Andrew Rowan, Cristina Naceur-Lombardelli, Oliver Shutkever, Carlos Martinez-Ruiz, Siân Harries, Sonya Hessey, Babu Naidu, James D. Brenton, et al.
- Corresponding author: Oriol Pich、Charles Swanton、Mariam Jamal-Hanjani、Nicholas McGranahan (The Francis Crick Institute / UCL Cancer Institute)
- 雑誌: Nature
- 発行年: 2026
- Epub日: 2025-12-10
- Article種別: Original Article
- PMID: 41372419
背景
がん患者の治療成績の向上に伴い、長期生存者が増加するなか、がん治療が将来の二次がんリスクや正常臓器に与える影響への関心が高まっている。これまで高深度変異解析により、組織学的に正常な組織にも多数のドライバー変異を持つクローンが存在することが明らかにされてきた (Alexandrov et al. Nature 2013)。加齢と正常組織のクローン拡大の関係は皮膚・膀胱・結腸・肝臓・食道・脳など多くの組織で示されており、SBS5 (substitution-base signature 5; age-related) シグネチャーが主要な加齢変異プロセスとして同定されている (Alexandrov et al. Nature 2020)。
しかしながら、これらの先行研究の大部分は健常加齢個体を対象としており、がん治療という人工的な変異誘発源が正常組織の体細胞進化に与える影響については研究が不足しており、系統的な解析は手薄であった。腫瘍においては白金製剤・アルキル化剤・放射線療法が特徴的な変異シグネチャーを残すことは知られているが、これらの治療が血液・脳・肝臓などの非腫瘍性正常組織に与える変異量の組織差、さらには免疫チェックポイント阻害薬のような非変異誘発性治療が正常組織の選択圧として機能するか否かは未解明であった (McGranahan et al. Cell 2017)。正常組織での治療誘発変異の定量化と組織間差の理解は、将来の二次がんリスクの予測や治療選択の最適化に不可欠であり、この領域の知見の不足が課題であった。
目的
転移がん患者の剖検時に採取された正常組織を対象に、(1) がん治療を含む外因性要因が正常細胞の変異獲得に与える影響とその組織特異性、(2) 治療誘発変異シグネチャーの特定と定量、(3) 正常組織における組織特異的正の選択とドライバー変異の帰属、および (4) 免疫療法の非変異誘発性選択圧としての役割を検証する。
結果
所見1：正常組織の体細胞変異全体像と組織間差：
コホート全体で166,732個のSBS変異、16,108個の挿入欠失変異 (ID) および2,399個のDBS変異を同定した。変異の大部分は低VAFで検出され (中央値VAF 0.0000323、IQR 0.000027-0.000042)、少数のクローン集団に限局していることが示された (Fig. 1b)。サンプルあたりの変異数は305-2,854と個体差・組織差が大きく、変異数とシーケンシング深度には有意な正の相関があった (r=0.26、P=1.16×10⁻³)。変異率の分散は組織種別が50% (P=4.7×10⁻²⁴)、患者が20% (P=3.26×10⁻¹⁰)、年齢が6%を説明した (Fig. 1、Extended Data Fig. 2e)。肝臓サンプルは最高の変異量 (1.83、IQR 1.71-2.14 変異/Mb/細胞) を示し、未治療肝細胞がんの変異量と同等であった。一方、脳 (0.58、IQR 0.47-0.64) および下垂体組織 (0.50、IQR 0.47-0.54) は最低値を示した (Fig. 1b)。de novo変異シグネチャー解析により16種類のSBSシグネチャーを同定し、うち6種はCOSMIC (Catalogue of Somatic Mutations in Cancer) データベースの既知シグネチャーと一致した。一致したシグネチャーにはSBS4 (substitution-base signature; 喫煙・加齢関連) 等のSBS4/5/40/31/35/25が含まれ、10種のSBS-A〜SBS-J (novel substitution-base signatures A-to-J) が新規として同定された (Fig. 2a)。
所見2：加齢・生活習慣関連変異プロセスの組織特異性：
コホートを通じて最も広く認められた変異シグネチャーはSBS5 (加齢時計様プロセス) であり、サンプルあたり平均40% (IQR 26-52%) の変異に寄与していた。組織ごとの年間SBS5変異獲得率は肺で約43変異/細胞/年と最も高く、血液 (約24)、脳 (約21) と組織間で有意差が認められた (P=1.252×10⁻¹³) (Fig. 2d,e)。喫煙関連シグネチャーSBS4は肺・肝臓・心臓組織に同定され、全3組織でpack years (1日の喫煙箱数×喫煙年数) と有意な相関を示した。肺では1 pack yearあたり約20変異/細胞が蓄積され、心臓・肝臓では約5変異/細胞と推定された (Fig. 2f,g)。禁煙者のみに限定した解析でも有意な相関が維持され (R²=0.55、P=0.0005)、喫煙関連変異が禁煙後も組織内に蓄積したまま残存することが示された。アルコール関連シグネチャーSBS-B (novel substitution-base signature B; alcohol-associated) は16/16肝臓サンプルで同定され (1サンプル最大38%)、アルコール摂取量と有意に相関した (P=0.006、Jonckheere-Terpstra検定; Fig. 2h)。
所見3：がん治療による変異フットプリント：
白金製剤化学療法関連のSBS31/SBS35 (substitution-base; platinum-induced) は投与コース数と有意な相関を示し、腎臓サンプルで特に顕著であった (r²=0.94、P=0.00026) (Fig. 3a,b)。組織別の白金変異誘発率は血液が最高 (89変異/細胞/白金コース) であり、脳・副腎では最低値 (<27変異/細胞/コース) であった。血液では6コースの白金化学療法が約27加齢年相当の変異蓄積を引き起こすと推定され、同じ6コースが肺では4.47加齢年相当にとどまるという顕著な組織差が示された (Fig. 4a)。テモゾロミド関連シグネチャーSBS-Fは治療を受けた4患者の39サンプル中18サンプルで同定され、変異量は組織種別 (P=0.03) およびテモゾロミド投与日数 (P=2.64×10⁻⁵) と有意に相関した (Fig. 3d)。脾臓・下垂体・心臓でテモゾロミド関連変異の負荷が最大 (中央値20%、範囲5-43%) であった。放射線関連挿入欠失シグネチャーID-E (insertion-deletion signature E; radiation-induced) は放射線治療サンプルで有意に富化され (平均0.04挿入欠失/細胞/Mb、Mann-Whitney U検定)、主に脳・下垂体・肺組織で検出された (Fig. 3e)。免疫療法 (抗CTLA4・抗PD1・併用) については、いずれも変異率の有意な増加と関連しなかった (抗CTLA4: P=0.14; 抗PD1: P=0.049で負の相関; 併用: P=0.49) (Fig. 5e)。
所見4：組織特異的正の選択とドライバー変異：
dNdSc解析によりコホート全体で最も広く正の選択シグナルを示したのはTP53であり、156サンプル中34サンプル (21%) で有意なミスセンス変異過剰が認められた。注目すべきことに、コホート内で同定されたTP53ミスセンス変異の総数は2,439個であり、TCGA (The Cancer Genome Atlas) の全がん種10,295腫瘍サンプルにおけるTP53ミスセンス変異数2,314個に匹敵した (Fig. 5a,b、Extended Data Fig. 7a)。肺組織では肺がんドライバーに対応するPTENおよびPIK3CAの正の選択が認められ、一部サンプルではB2MおよびEGFR変異も観察された。肝臓ではNFE2L2 (nuclear factor erythroid 2-like 2; 肝細胞がんドライバー) とCTNNB1 (cadherin transducer; β-カテニン遺伝子) の優先的拡大が示され、45番コドンのホットスポット変異も確認された。脾臓ではNOTCH2ナンセンス変異とBRAFミスセンス変異 (特にD594、n=17; BRAF V600Eより高頻度) が同定された。血液ではCHEK2 (checkpoint-homolog effector kinase-2)・DNMT3A・PPM1D (protein phosphatase Mg-dependent 1D) 等のクローン性造血関連遺伝子における選択が認められたが、その変異VAFは標準シーケンシングの検出限界以下であった (Extended Data Fig. 9)。
所見5：免疫療法による非変異誘発性選択圧とドライバー変異帰属：
免疫療法は変異率を増加させないにもかかわらず、PPM1Dドライバー変異との有意な正の関連を示した (抗CTLA4: P=0.0022; 抗PD1: P=0.020) (Fig. 5e,f)。TP53ドライバー変異も抗CTLA4治療と有意に関連した (P=0.00957)。免疫療法治療群でTP53・CHEK2・PPM1Dに対するdN/dS値の上昇傾向が認められ、免疫療法が変異誘発を伴わずに特定クローンの選択圧として機能することが示された。治療誘発ドライバー変異の帰属解析では、白金製剤曝露サンプルの正常肺組織TP53ドライバー変異の26% (50件中13件)、血液のTP53ドライバーの20% (234件中46件) がSBS31 (白金シグネチャー) に帰属された。テモゾロミド関連シグネチャーSBS-Fはテモゾロミド曝露血液サンプルのTP53ドライバーの45% (51件中23件)、脾臓の35% (37件中13件) に帰属された。システム全体として、正常組織のドライバー変異の25%超が治療に起因することが示唆された (Fig. 5g)。
考察/結論
本研究はがん患者22名から採取した168正常組織の高深度duplex sequencingという新規な大規模リソースを活用し、がん治療が正常組織の体細胞進化に与える影響を初めて体系的に解析した。得られた知見は、正常組織の変異蓄積が単純な加齢プロセスの産物ではなく、複数の外因性要因 (治療・喫煙・アルコール) と内因性プロセスの複合的相互作用によって規定されることを示す。
① 先行研究との違い: これまでの変異蓄積研究は健常個体の加齢変異解析に集中しており、がん治療を受けた患者の正常組織における治療誘発変異の系統的定量は行われていなかった (Alexandrov et al. Nature 2013)。本研究は治療後患者の正常組織を対象とし、白金製剤・テモゾロミドなど治療特異的変異シグネチャーを正常組織で初めて定量化した点がこれまでと異なる。特に血液での6コース白金化学療法が27加齢年相当の変異を誘発するという定量的推定は、がん治療の長期毒性評価に新たな枠組みを提供する。
② 新規性: 免疫チェックポイント阻害薬が変異誘発を伴わずに正常組織のTP53・PPM1Dドライバー変異の選択圧として機能するという発見は新規な概念である (Lopez-Bigas et al. Nature 2026)。変異誘発なき選択圧 (non-mutagenic selection pressure) として免疫療法が機能することは、本研究で初めて正常組織レベルで直接示された。PEACE研究というユニークな剖検コホートによって、通常では解析不可能な複数臓器の同時解析が実現したことも新規な研究基盤である。
③ 臨床応用: 治療誘発変異の定量化は臨床的意義が大きい。若年男性への白金製剤投与は将来の固形がんリスクを高める (ハザード比2.4) ことが知られているが、本研究は血液細胞での高変異誘発を分子レベルで裏付けた。また、免疫療法によるPPM1Dクローン選択は将来の二次がん・血液悪性腫瘍リスクとの関連が懸念されることから、臨床応用として治療後の液体生検を用いたクローン性造血モニタリングが提唱される (McGranahan et al. Cell 2017)。肺での白金6サイクル=4.47加齢年相当は喫煙40 pack years=20加齢年相当を大幅に下回り、喫煙の影響が変異誘発に留まらない非変異誘発性発がん促進機序を含むことを示唆する。
④ 残された課題: 本研究最大の技術的限界はduplex sequencingが細胞種分解能を持たない点である。組織間の変異率差には異なる細胞種組成 (cell-type composition) が寄与し得るため、今後の研究では単一細胞レベルの変異解析による細胞種別評価が必要である。また、コホートサイズ22名という規模は統計的検出力に限界があり、特にまれな変異シグネチャーの組織分布については更なる検討を要する。免疫療法によるPPM1D/TP53クローン選択の機序 (免疫活性化が特定クローンに選択優位をもたらす機構) についても今後の検討が必要である。さらに、正常組織での治療誘発変異蓄積と将来の二次がん発症リスクの定量的関連を確立するための前向き longitudinal研究の構築が今後の課題である。
方法
PEACE (Posthumous Evaluation of Advanced Cancer Environment) 研究 (ClinicalTrials.gov NCT03004755; 13/Lo/0972) の剖検サンプルを用いた横断的基礎研究。転移がん患者22名 (中央値年齢65歳、範囲29-78歳) から168の正常組織サンプル (16臓器種別) を採取。2名については治療前の正常組織サンプルが TRACERx 研究 (NCT01888601) から補完的に入手可能であった。
変異解析法: 82.5-kb カスタム duplex sequencing (二重鎖シーケンシング) パネルを用い、30種の既知がん遺伝子 (TP53、EGFR、KRAS、BRAFを含む) と変異頻度の高いゲノム領域を標的とした。中央値シーケンシングカバレッジ 30,504× (IQR 27,684-36,450)。Duplex sequencingはユニーク分子識別子 (UMI) と二本鎖情報を用いてアーティファクト変異を除去し、低変異アレル頻度 (variant allele frequency, VAF) での高感度変異検出を可能とする。171サンプルに標準パネルを使用し、12サンプルには変異シグネチャー解析専用の小型パネルを使用 (後者は変異シグネチャー解析のみに供した)。品質管理により10サンプルを除外。
統計解析: de novo変異シグネチャー解析にはベイズ階層的ディリクレ過程 (Bayesian hierarchical-Dirichlet process) を適用し、一塩基置換 (SBS)・挿入欠失 (ID)・二塩基置換 (DBS) を独立に解析。正の選択解析には dNdScv (カバレッジ調整参照パネル) を用い、非同義/同義変異率 (dN/dS) を算出。治療効果の定量化には線形混合効果モデルを用い、患者・組織・腫瘍種を混合効果として補正した。Mann-Whitney U検定、Jonckheere-Terpstra検定、Pearson/Spearman相関を適宜使用。