- 著者: Somnath Tagore, Lindsay Caprio, Amit Dipak Amin, Kresimir Bestak, Karan Luthria, Edridge D’Souza, Irving Barrera, Johannes C. Melms, Sharon Wu, Sinan Abuzaid, Yiping Wang, Viktoria Jakubikova, Peter Koch, D. Zack Brodtman, Banpreet Bawa, Sachin K. Deshmukh, Leon Ebel, Miguel A. Ibarra-Arellano, Abhinav Jaiswal, Carino Gurjao, Jana Biermann, Neha Shaikh, Priyanka Ramaradj, Yohanna Georgis, Galina G. Lagos, Matthew I. Ehrlich, Patricia Ho, Zachary H. Walsh, Meri Rogava, Michelle Garlin Politis, Devanik Biswas, Azzurra Cottarelli, Nikhil Rizvi, Catherine A. Shu, Benjamin Herzberg, Niroshana Anandasabapathy, George Sledge, Emmanuel Zorn, Peter Canoll, Jeffrey N. Bruce, Naiyer A. Rizvi, Alison M. Taylor, Anjali Saqi, Hanina Hibshoosh, Gary K. Schwartz, Brian S. Henick, Fei Chen, Denis Schapiro, Parin Shah, Benjamin Izar
- Corresponding author: Benjamin Izar (Columbia University Irving Medical Center, New York, NY, USA)
- 雑誌: Nature Medicine
- 発行年: 2025
- Epub日: 2025-02-27
- Article種別: Original Article
- PMID: 40016452
背景
非小細胞肺がん (NSCLC) は脳転移 (BM) の主要な原因であり、米国では年間約23万人のNSCLC新規診断患者のうち、診断時に約10%、経過中に最大40%がBMを発症すると報告されている。免疫チェックポイント阻害薬 (ICB) やEGFR/ALK標的薬の導入により一部の患者の予後は改善したものの、多くのBM患者の予後は依然として不良である。BMの分子基盤に関する理解は未だ不十分であり、その複雑な生態系の解明が喫緊の課題である。
これまでの遺伝子解析研究は、主に治療後の検体を用いたバルクプロファイリングに依存していた。このアプローチでは、治療による耐性クローンの選択という交絡因子が存在し、腫瘍内の不均一性や非悪性細胞との相互作用を十分に捉えることができなかった。例えば、Tirosh et al. Science 2016 や vanDerLeun et al. NatRevCancer 2020 などの先行研究は、単一細胞解析の重要性を示唆しているが、NSCLC BMにおける未治療検体の包括的な解析は不足していた。
メラノーマBMに関する最近の単一細胞解析では、未治療検体において頭蓋外転移とは異なる特徴が示されたが、肺腺癌 (LUAD) BMの生態系に関する高解像度の理解は不十分なままであった。凍結組織からアーティファクトの少ない単一核マルチオミクス解析を可能にする技術進歩は、この知識ギャップを埋める上で重要な役割を果たす。本研究は、これらの技術的進歩を活用し、未治療NSCLC BMの包括的な分子・空間ランドスケープを初めて構築することを目的とした。これにより、BMの発生と進展に関わる生物学的メカニズムの解明に貢献し、新たな治療標的の同定に繋がる知見を提供することが期待される。
目的
本研究の目的は、未治療NSCLC (LUAD) 患者の原発腫瘍 (PT) と脳転移 (BM) 組織に対し、単一核RNAシーケンス (snRNA-seq)、T細胞受容体 (TCR) シーケンス、空間トランスクリプトミクス、および低カバレッジ全ゲノムシーケンス (lp-WGS) を統合的に解析することである。これにより、BMに特異的なゲノム不安定性、腫瘍細胞状態、および腫瘍微小環境 (TME) の再構成を同定する。さらに、これらの発見を独立した大規模コホート (合計約17,000例) と多重免疫蛍光ペアコホート (18組) で検証し、BMの治療標的となりうる生物学的特徴を抽出することを目指す。最終的には、NSCLC BMの生物学に関する基礎的な洞察を提供し、今後の研究のための重要なリソースを構築する。
結果
BMにおける染色体不安定性 (CIN) の顕著な濃縮: lp-WGSとsnRNA-seqによるCNA検出の一致度はPearson R=0.73 (p=7.8×10⁻⁸) であった。FGAはBMでPTより有意に高く (Wilcoxon p=0.023)、CIN70スコアもBMで有意に上昇した (Wilcoxon p=2.1×10⁻⁶)。4,710例のNSCLCコホート、38例および58例のマッチWES 2コホート、Caris Life Sciencesの12,575例 (LUAD/LUSC) のLOHベース解析においても、BMが最もCINが高いことを示した (Fig 2a-j)。整合して、cGAS (cyclic GMP-AMP synthase)+小核の免疫蛍光定量では、BMはPTより有意に増加しており (p=3.04×10⁻²)、CIN高BMではSTING (Stimulator of interferon genes) 経路が抑制されていることが示された。CINはNSCLC患者の全生存期間の不良と関連していた (Extended Data Fig. 2c)。
神経様転写プログラムの出現: 非負行列因子分解 (NMF) により8つの再現性あるメタプログラム (MP) を同定した。このうち4つ (MP4, MP5, MP7, MP8) はBMに排他的に存在し、CIN70シグネチャが濃縮されていた (Fig 3c)。MP7/8は神経前駆細胞やオリゴデンドロサイト前駆細胞の転写ハルマークと強く相関し、原発性脳腫瘍と共通の特徴を示した。MP4は既存の脳転移を含まない参照における「unassigned」プログラムに対応するBM特異的プログラムであった (Fig 3e)。外部scRNA-seqデータセット Kim et al. NatCommun 2020 を用いた検証では、BM特異的MPはBMでのみ有意に濃縮されており、一般的な転移プログラムではないことが示された (Fig 3f)。
Cluster 21 — 患者横断的に共有されるBM濃縮集団: 悪性細胞154,771個のunsupervised clusteringにより、全患者由来の希少な悪性細胞集団「cluster 21」を同定した (全悪性細胞の1.81%、n=2,812 cells) (Fig 4a-c)。このcluster 21はBMで有意に濃縮されており (Fisher p<0.001)、CIN70スコアが高値であった (Extended Data Fig. 4j)。また、EMT (上皮間葉転換) 制御因子 (ZEB1, ZEB2, VIM, SPARCL1) および神経様分化マーカー (PLXDC2, PCDH9, GFAP, MBP, TCF4) を共発現する一方、肺上皮lineageマーカー (NKX2-1, SFTPB, EPCAM) は抑制されていた (Fig 4d)。18組のペアPT/BM検体を用いた多重免疫蛍光染色では、EpCAM+GFAP+細胞は全癌細胞の約0.12%と希少であったが、その大部分がBMに存在した (p=1×10⁻⁴) (Fig 4e,f)。cluster 21のシグネチャスコアが高い患者は全生存予後が不良であり (Extended Data Fig. 4o)、この集団がBM特異的な高度な可塑性を持つ脱分化集団として機能することが示唆された。Markov吸収確率を用いた解析では、悪性細胞の82.95%がmixed-lineageプログラムを示し、cluster 21は特に正常参照細胞から離れた状態にあることが明らかになった (Fig 4j)。
TME — T細胞、骨髄系細胞、CNS細胞の再構成: T/NK細胞 (n=28,646 cells) の解析では、CD8+T細胞の活性化・分化・記憶サブセットがBMで有意に減少し、Treg細胞、ナイーブT細胞、NK細胞が増加していた (Fig 5c)。多重免疫蛍光染色でも、BMにおけるCD8+TIL浸潤はPTより低いことが確認された (p=4×10⁻³) (Fig 5e)。TCR解析ではクローン拡大は限定的であった (GINI指数=0.013-0.08)。骨髄系細胞 (n=25,741 cells) の解析では、PTは主に肺組織常在性の肺胞マクロファージ (alveolar macrophage) と肥満細胞で構成されるのに対し、BMは単球/単球由来マクロファージ (MDM)、IFN応答性マクロファージ (cluster 1)、プロ腫瘍性マーカー (CD163, CD163L, MERTK, SIGLEC1/10) を発現するマクロファージ、好中球、ミクログリアを含む多様な骨髄コンパートメントを呈した (Fig 5h)。BMではCD163タンパク質の発現強度がPTより高かった (Fig 5j)。
CNS細胞・神経ニッチ相互作用: CNS細胞 (n=17,941 cells) は6つの主要な亜集団に分類された。ContactTracingを用いた細胞間相互作用解析では、骨髄系細胞由来のPECAM1とCNS細胞由来のCD38の相互作用がBMニッチで示唆され、これは血管透過性や非カノニカルプリン作動性シグナル伝達を介した免疫抑制的なTMEを駆動する可能性がある (Fig 5l)。CNS細胞と悪性細胞間の相互作用では、EMT、血管新生、低酸素経路の濃縮が認められた (Extended Data Fig. 5g)。また、ニッチ細胞とT細胞間のIL-6/JAK/STAT3シグナル軸の濃縮も確認され、これはプロ腫瘍性TMEに寄与する可能性がある。骨髄系細胞由来のCCL3およびCCL5とBM癌細胞との相互作用も同定され、これらは血管新生、EMT、増殖を促進し、ミクログリアを免疫抑制表現型に極性化させることが示唆された。
空間的変動性: Slide-seqV2プラットフォームを用いた空間トランスクリプトミクス解析では、cluster 21シグネチャを発現する細胞が血管周囲空間に濃縮されていることが示された (Fig 6d-g)。これらの細胞はCIN70シグネチャの発現と強く相関しており (Pearson R=0.68, p=7.2×10⁻³)、高度に転移性の細胞集団であるという考えをさらに裏付けた (Fig 6p)。骨髄系細胞密度の高いBM領域では、MYC標的遺伝子、TGF-βシグナル伝達、EMTの濃縮が認められた一方、PTではインターフェロンγ応答、インターフェロンα応答、抗原提示、テルペノイド代謝の濃縮が認められた (Fig 6j,k,n,o)。
考察/結論
本研究は、NSCLC脳転移 (BM) の包括的なマルチモーダル単一細胞・空間・ゲノムリソースを提供し、BMが単なる原発巣の遠隔コピーではなく、独自の生物学的特徴を持つ生態系であることを明らかにした。
先行研究との違い: これまでの研究は治療済み検体のバルクプロファイルに依存し、腫瘍内不均一性や非悪性細胞との相互作用を捉えきれていなかった。本研究は、未治療検体を用いた単一細胞・空間解析により、Achrol et al. NatRevDisPrimers 2019 が指摘するBMの複雑性を高解像度で解明した点で、これまでの知見と対照的である。特に、Bakhoum et al. Nature 2018 が報告したCINと転移の関連性を、NSCLC BMにおいて大規模コホートで検証し、cGAS+小核の免疫蛍光染色により機能的にも裏付けた点が強みである。
新規性: 本研究で初めて、BMにおいて染色体不安定性 (CIN) が顕著に高く、神経様メタプログラム (MP7/8) とEMT・神経分化マーカーを共発現する希少な悪性細胞集団「cluster 21」が脳転移で強く濃縮されていることを新規に同定した。このcluster 21は、全悪性細胞のわずか1.81%を占めるに過ぎないが、患者横断的に共有され、BMの発生と進展に重要な役割を果たす可能性が示唆された。また、BMの腫瘍微小環境 (TME) がCD8+T細胞の減少、CD163高発現M2マクロファージの存在、および神経系細胞と骨髄系細胞の相互作用を伴う免疫抑制的な特徴を持つことを明らかにした点も新規性である。
臨床応用: 本研究の知見は、NSCLC BMの診断、予後予測、および治療戦略の開発に臨床応用される可能性がある。(1) CINはBMのリスクおよび予後バイオマーカーとして患者層別化に活用できる。(2) CIN高BMで抑制されているSTING経路を再活性化する治療法 (STINGアゴニストとICBの併用など) が有効な可能性がある。(3) cluster 21の神経様集団を標的とするEMT/神経分化阻害剤 (TCF4, ZEB1/2, GFAP+集団を標的) の開発が期待される。(4) PT主体の肺胞マクロファージとBMのプロ腫瘍性マクロファージ・好中球の違いを標的とした骨髄系細胞の再プログラミングが治療選択肢となりうる。(5) IL-6/JAK/STAT3軸やCCL3/CCL5–骨髄系細胞–癌細胞相互作用の遮断がBMの免疫抑制的なTMEを改善する可能性がある。
残された課題: 今後の検討課題として、cluster 21の原因となる遺伝子ドライバーの特定、治療後のBMにおける動態変化の解析、LUSCやSCLCなどの他の肺がんタイプへの一般化、および前向きコホートでの治療応答予測性能の検証が挙げられる。また、本研究は時間的なスナップショットであり、BMの全進展過程における癌細胞内在性およびTMEの特徴をより深く理解するためには、ヒトBMを忠実に再現する堅牢な免疫能を有する前臨床モデルの開発が必要である。
方法
本研究では、未治療NSCLC患者43検体 (PT 12例、BM 31例、うち2例はペアPT/BM) を対象に、snRNA-seq、matched TCR-seq、single-cell spatial transcriptomics、およびpool-matched lp-WGSを適用した。厳格な品質管理 (QC) 後、274,817個の単一細胞トランスクリプトームを解析に含めた。悪性細胞の特定には、lp-WGS、inferCNV、および遺伝子発現を統合した分類器を使用し、全細胞の55.83%が悪性細胞として同定された。
主要な8つの細胞型を同定した。これには悪性細胞 (n=154,771 cells)、T/NK細胞 (n=28,646 cells)、骨髄系細胞 (n=25,741 cells)、B/形質細胞 (n=12,110 cells)、中枢神経系 (CNS) 細胞 (n=17,941 cells)、内皮細胞 (n=5,164 cells)、間質細胞 (n=15,412 cells)、およびその他の細胞 (n=15,032 cells) が含まれる。染色体不安定性 (CIN) の評価は、FGA (fraction of genome altered) とCIN70 (chromosomal instability 70-gene signature) mRNAシグネチャを用いて行った。これらの評価結果は、Caris Life Sciencesの12,575例 (PT 8,508例、BM 708例、ECM 3,359例)、TCGAの4,710例、およびWESの86例からなる独立コホートで検証された。
転写プログラムの同定には、非負行列因子分解 (NMF) を用いてメタプログラム (MP) を抽出し、Markov吸収確率を用いて健常肺細胞参照との近接性を算出した。細胞間相互作用の推定にはContactTracingツールを使用し、TCR-seqデータを用いてT細胞のクローン多様性を解析した。さらに、18組のペアPT/BM検体に対し、EpCAM、cGAS、GFAP、CD8、CD163の多重免疫蛍光染色を行い、これらのマーカーの発現を定量的に評価した。統計解析には、Wilcoxon検定、Kruskal-Wallis検定、Fisherの正確確率検定、Pearson相関係数、log-rank検定などが用いられた。