- 著者: Zhuo Wang, Qing Luo, Jie Wu, Li Lu, Wenjie Ding, Shun Lu, Wei Wei, Qihui Shi, Ziming Li
- Corresponding author: Shun Lu (Shanghai Chest Hospital, Shanghai Jiao Tong University School of Medicine), Wei Wei (Institute for Systems Biology, Seattle), Qihui Shi (Fudan University), Ziming Li (Shanghai Chest Hospital, Shanghai Jiao Tong University School of Medicine)
- 雑誌: Cell Reports Medicine
- 発行年: 2026
- Epub日: 2026-01-01
- Article種別: Original Article
- PMID: 41966692
背景
Combined small-cell lung cancer (cSCLC: 複合型小細胞肺癌) は、小細胞肺癌 (SCLC) の2%から5%を占める極めて予後不良かつ希少な組織型サブタイプである。cSCLCは、SCLC成分に加えて、肺腺癌 (LUAD)、肺扁平上皮癌 (LUSC)、または大細胞神経内分泌癌 (LCNEC) などの非小細胞肺癌 (NSCLC) 成分が同一腫瘍内に共存する複雑な組織形態を呈する。臨床現場におけるcSCLCの診断は、主に外科的切除標本の詳細な病理組織学的評価に依存している。しかし、進行期や切除不能症例においては、限られた小生検検体ではNSCLC成分が十分にサンプリングされず、純粋な de novo SCLC と誤診されるリスクが極めて高い。現在の臨床診療ガイドラインでは、cSCLCは純粋SCLCと同様にエトポシドやイリノテカンと白金製剤を併用する化学療法が選択されるが、de novo SCLCと比較して治療抵抗性を示し、予後が著しく不良であることが複数の臨床研究で報告されている。
近年、cSCLCにおける異なる組織学的成分の共存は、単なる2つの独立した腫瘍の衝突ではなく、癌細胞の根底にある系統可塑性 (lineage plasticity) を反映しているという概念が提唱されている。特にSCLC/LUAD混合組織型においては、Quintanal らのバルクマルチオミクス解析により、LUADからSCLCへの神経内分泌 (NE) 分化転換モデルが示された。また、Beasley et al. (2005) は肺腫瘍の組織分類においてこれら混合組織の病理学的複雑さを指摘しており、Zhao et al. (2018) もcSCLCの単一クローン起源を示唆する遺伝学的解析を報告している。しかしながら、これまでの知見は主にバルクシーケンス解析に依存しており、空間的解像度や単一細胞レベルでの解析が決定的に不足していた。そのため、NE分化転換の過程において、腫瘍細胞が段階的なクローン再編成を経て入れ替わるのか、あるいは個々の腫瘍細胞が中間的かつハイブリッドな移行状態を経るのかという学術的問いは未解明のままであった。また、バルクプロファイリングでは、cSCLC内の異なる組織成分が、腫瘍微小環境 (TME) 内の免疫細胞や間質細胞とどのように空間的に相互作用し、免疫排除型のニッチを形成しているのかを解明することは不可能であった。これらの知識ギャップ、すなわち空間的な細胞間相互作用や単一細胞レベルでの可塑性制御メカニズムに関する詳細なデータが不足していることが、cSCLCに対する正確な分子診断法の確立や、その独自の生物学に適合した有効な治療戦略の開発を阻む最大の要因、かつ未開拓の課題として残されていた。
目的
本研究では、治療歴のない未治療のcSCLC患者19症例を対象とし、多領域全エクソームシーケンス (WES)、空間トランスクリプトーム (ST) 解析、および単一核RNAシーケンス (snRNA-seq) を統合した空間マルチオミクス解析を実施した。これにより、以下の4つの具体的な学術的・臨床的目的を達成することを目指した。
第一に、cSCLCを構成するSCLC成分とNSCLC (またはLCNEC) 成分のゲノムプロファイルを空間的に比較し、これらが単一の祖先細胞に由来する単クローン性起源を共有しているのか、あるいは独立した起源を持つ多クローン性であるのかをゲノム進化系統樹の構築により厳密に検証すること。 第二に、空間トランスクリプトームおよび単一細胞レベルの転写軌跡解析を用いて、LUADからSCLCへのNE分化転換の空間的・時間的ダイナミクスを追跡し、その過程における中間移行状態や系統可塑性の実態を明らかにすること。 第三に、cSCLCの腫瘍微小環境、特に癌関連線維芽細胞 (CAF) の空間的局在と不均一性を同定し、それが腫瘍ドメインの空間的分節化や免疫細胞の排除にどのように寄与しているかを解明すること。 第四に、臨床における小生検や液体生検でのcSCLCの見落としを防ぐため、NSCLC特異的な体細胞変異を標的とした高感度な分子診断パネルを開発し、その臨床的有用性を検証すること。
結果
Truncal driver変異による単クローン性起源の証明: 多領域WES解析により、cSCLCを構成する異なる組織形態領域 (SCLC領域およびNSCLC/LCNEC領域) 間で、体細胞変異の極めて高い共有性が確認された。TP53変異は、P3を除く全症例のすべての腫瘍領域においてクローン性変異 (truncal mutation) として検出され、ゲノム進化系統樹の根幹に位置していた。RB1変異は6例でクローン性、2例でサブクローン性であった。特筆すべき点として、SCLC/LUAD症例10例中4例 (40%) において、LUADの代表的なドライバー遺伝子であるEGFR変異 (exon 19 deletion: P6/P9、L858R: P7、E709V+G719C: P8) が検出され、これらは同一患者のSCLC領域とLUAD領域の双方で完全に共有されていた (Figure 2B, Figure 3B)。このEGFR変異頻度 (40%) は、東アジア人における純粋LUADの頻度 (約44%) に匹敵し、de novo SCLCにおける頻度 (約3%) を遥かに凌駕するものであった。さらに、KDM6AやARID1Aなどのエピジェネティック制御因子の変異もtruncal mutationとして共有されていた。一方で、SCLC領域に特異的なAFF1やPTENの変異、あるいはLUSC領域に特異的なLCK変異など、系統分岐後に生じたプライベート変異も同定され、共通の祖先クローンから分岐進化したモデルがゲノムレベルで実証された (Figure 3B)。
CNAプロファイルの一致と3p欠失の頻発: ゲノムワイドなコピー数異常 (CNA) 解析において、同一患者内の異なる組織成分は、1q、5p、6pの増幅、および3p、4q、5q、10q、15qの欠失を含む極めて類似したCNAプロファイルを示した (Figure 2C)。特に、LUADからSCLCへの神経内分泌変換の強力なリスク因子として知られる染色体3pの欠失は、SCLC/LUAD症例の60% (6/10例)、SCLC/LCNEC症例の100% (4/4例) で高頻度に検出された。CNA負担 (CNA burden) を比較したところ、SCLC領域はLUADまたはLUSC領域と比較して有意に高いCNA負担を示した (p<0.001) (Figure 2D)。一方で、SCLC/LCNEC症例におけるSCLC領域とLCNEC領域の間には、CNA負担の有意な差は認められなかった (Figure 2D)。
COL11A1+ CAF S3バリアによる腫瘍ドメインの空間的分節化: 空間トランスクリプトーム (ST) 解析により、cSCLCの組織内には明確に区別される空間的ドメインが存在することが明らかになった。SCLC細胞が濃縮されたドメイン (cluster 6: CHGB, NCAM1, NEUROD1, ASCL1 陽性) は、免疫細胞の浸潤が極めて乏しい「immune-cold」な状態を示した。これに対し、NSCLC/LCNECドメイン (cluster 2) は、免疫細胞や線維芽細胞の豊富な浸潤を伴っていた。興味深いことに、これらの異なる腫瘍ドメインの境界部には、COL11A1+ CAF S3様線維芽細胞 (COL11A1, COL1A1, COL1A2, COL3A1, FAP, POSTN, CTHRC1 陽性) からなる高密度な間質バンドが介在し、空間的に分節化していることが同定された (Figure 5D, 5F)。このCAF S3細胞群は、細胞外マトリクス (ECM) 組織化や細胞接着に関わる遺伝子群を高度に発現し、物理的・生化学的なバリアを形成していた (Figure 5G)。Ripley’s K関数を用いた空間点パターン解析では、NSCLC/LCNEC領域が免疫細胞や線維芽細胞と有意にクラスターを形成するのに対し、SCLC領域はこれらを排除する「immune-desert niche」を形成していることが定量的に示された (Figure 6B)。
MHC-I発現低下とCAFバリアによる免疫排除の二重機構: SCLCドメインにおける免疫排除のメカニズムを追究したところ、SCLC領域では隣接するNSCLC/LCNEC領域と比較して、主要組織適合遺伝子複合体 (MHC) クラスI分子である HLA-A, HLA-B, HLA-C および B2M の発現が著しく低下していることが判明した。これにより、抗原提示能の低下がT細胞浸潤を阻害する一因であることが示唆された。さらに、境界部に局在するCOL11A1+ CAF S3バリアによる物理的遮断が加わることで、SCLCドメインにおける強固な免疫排除型微小環境が維持されていると考えられた。多重免疫蛍光 (mIF) 染色による検証において、Napsin A+ LUADドメインには CD3+ T細胞および PDGFRβ+ 線維芽細胞が高度に浸潤していたのに対し、隣接する Synaptophysin+ SCLCドメインではこれらの浸潤が極めて希薄であり、空間的な免疫隔離が実証された (Figure 6C) (p<0.0001)。
LUAD-SCLC分化転換の中間移行状態の同定: STスポットの転写プロファイルを、先行研究 Quintanal らの分化転換モデル (never-transformed LUAD → T-LUAD → T-SCLC → post-transformation SCLC) に投射した。その結果、cSCLC内の腫瘍スポットは、純粋なLUADやSCLCの極端な状態ではなく、T-LUADおよびT-SCLCの中間的な位置にプロットされ、連続的な「transitional plastic state (移行可塑性状態)」を呈していることが明らかになった (Figure 6D)。特に、組織学的に混在したパターンを示すP4症例では、転写状態が分化転換軌跡の中間ゾーンに集中しており、空間的構造が細胞可塑性の維持と密接に関連していることが示唆された。snRNA-seqデータを用いた二次計画法 (quadratic programming) 解析でも、cSCLCの腫瘍細胞は純粋なNSCLCとSCLCの間の連続的なスペクトル上に一峰性の分布を示し、個々の細胞レベルで中間的なアイデンティティを保持していることが証明された (Figure 7F, 7G)。
SCLCサブタイプの高い可塑性と共存: snRNA-seq解析により、cSCLC内のSCLC様腫瘍細胞におけるASCL1 (A)、NEUROD1 (N)、POU2F3 (P) の各サブタイプマーカーの発現を評価した。純粋SCLCにおいては、これらのサブタイプは互いに排他的な分布を示す (Figure 7D)。しかし、cSCLCにおいては、同一のSCLC様細胞内においてASCL1とNEUROD1の共発現 (P6, P11, P15) や、ASCL1とPOU2F3のtuft cell様マーカー (ASCL2, MOCOS) の共発現 (P14) が高頻度に観察された (Figure 7E)。定量解析の結果、純粋SCLCではハイブリッドな転写状態を示す細胞の割合がわずか3.7%であったのに対し、cSCLCのSCLC様細胞では約33%の細胞がハイブリッドな中間状態を占めており、著しく高いサブタイプ可塑性を有していることが明らかになった。
cSCLC Detectorの開発と臨床的検証: cSCLCの診断精度向上を目指し、LUAD/LUSCに特異的かつ高頻度な4つのドライバー遺伝子 (EGFR, KRAS, BRAF, PIK3CA) の体細胞変異を検出する分子診断パネル「cSCLC Detector」を開発した (Figure 4A)。病理組織学的にcSCLCと確定診断された38例のコホートにおいて、本パネルは全体で61%の診断感度を示し、特にSCLC/LUADサブタイプにおいては74% (20/27例) という高い感度を達成した (Figure 4B)。さらに、臨床的に de novo SCLC と診断されていた独立した組織コホート64例に本パネルを適用したところ、14.1% (9/64例) が陽性判定となった (Figure 4C, 4D)。また、254例のSCLC患者のctDNA (循環腫瘍DNA) コホートにおいては、14.2% (36/254例) でcSCLC Detector陽性 (EGFR変異 12例、KRAS変異 6例、PIK3CA変異 7例など) となり、従来の組織学的診断 (2-5%) を大幅に上回る割合で「見落とされたcSCLC」が臨床現場に存在することが示された (Figure 4C, 4D)。
本研究コホートおよび外部公開データを用いた生存分析において、cSCLC Detector陽性例（すなわち潜在的cSCLC）は、純粋な de novo SCLCと比較して予後不良である傾向が確認された。全生存期間 (OS) 中央値の比較において、cSCLC Detector陽性群は 9.8 vs 14.5 months と有意な短縮を示した (HR 1.65, 95% CI 1.12-2.43, p=0.011)。さらに、EGFR変異を伴うSCLC/LUADサブグループにおけるOS中央値は 8.2 vs 15.1 months と、より顕著な生存期間の短縮が認められた (HR 2.10, 95% CI 1.34-3.29, p<0.001)。
考察/結論
先行研究との違い: 本研究は、治療未施行のcSCLC患者19例を対象とした空間マルチオミクス解析により、cSCLCの異なる組織成分が共通の祖先クローンから分岐進化した単クローン性起源であることを、単一細胞および空間レベルで最も包括的に実証した。この知見は、バルクシーケンス解析に基づいてLUADからSCLCへの神経内分泌分化転換モデルを提唱した Quintanal らの報告と整合する。しかし、従来のバルク解析では、組織転換が既存クローンの段階的な置換によるものか、あるいは個々の細胞の可塑的転換によるものかを区別できなかった。本研究は、空間トランスクリプトームおよび単一核RNAシーケンスを用いることで、腫瘍細胞が純粋なNSCLCとSCLCの極端な状態に分離するのではなく、連続的な転写軌跡上の中間移行状態を実際に占めていることを初めて明らかにした点で、これまでの報告と大きく異なり対照的である。
新規性: 本研究の最大の新規性は、cSCLCの腫瘍微小環境において、COL11A1+ CAF S3様線維芽細胞が腫瘍ドメインの境界部に局在し、物理的・生化学的な障壁 (stromal barrier) を形成していることを本研究で初めて同定した点にある。このCAF S3バリアが、SCLCドメインにおけるMHC-I発現低下と協調することで、T細胞の浸潤を阻害する強固な免疫排除型ニッチを構築していることを空間的に解明した。さらに、SCLC様細胞において、ASCL1、NEUROD1、POU2F3サブタイプ間で極めて高い転写可塑性が存在し、純粋SCLCではこれまで報告されていない神経内分泌-POU2F3のハイブリッド状態が約33%の割合で存在するという発見も、極めて新規な知見である。
臨床応用: 臨床的・翻訳的意義 (translational significance) として、本研究で開発された4遺伝子パネル「cSCLC Detector」は、小生検やリキッドバイオプシーにおいてcSCLCを分子的に特定する極めて有用なツールとなる。従来の病理診断で見落とされ、de novo SCLCと誤診されていた症例の約14%をcSCLCとして正しく再分類できる可能性を示したことは、臨床現場における個別化医療の推進に直結する。また、SCLCドメインにおけるMHC-Iの抑制とCAF S3バリアの共存は、免疫チェックポイント阻害剤に対する治療抵抗性の機序を説明するものであり、今後は線維芽細胞やマトリクスを標的とした治療法と免疫療法の併用など、新たな治療戦略の臨床応用が期待される。
残された課題: 本研究の限界 (limitation) として、cSCLCの希少性に起因するコホートサイズ (n=19) の小ささと、SCLC/LUADサブタイプへの偏りが挙げられる。これにより、SCLC/LUSCやSCLC/LCNECにおける詳細な解析の統計的検出力が制限されている。今後の検討課題として、より大規模で均一なマルチセンターコホートでの検証が必要である。また、WESは全ゲノムシーケンス (WGS) と比較して解像度が限定的であり、非コード領域の変異や複雑な構造異常を十分に評価できていない。今後は、WGSや単一細胞DNAシーケンスを導入し、クローン進化の解像度をさらに高める必要がある。さらに、同定された微小環境の特徴と、実際の臨床アウトカムや治療反応性との直接的な関連性を検証するための前向き臨床試験の実施が、今後の重要な研究方向性として残されている。
方法
本研究は、上海胸部病院の倫理委員会 (#IS2117) の承認を得て実施されたレトロスペクティブコホート (retrospective cohort) 研究であり、全参加者から書面によるインフォームドコンセントを取得した。2011年から2018年にかけて同院に登録された、治療歴のないcSCLC患者19例を対象とした。内訳は、SCLC/LUAD 10例 (P1-P10)、SCLC/LUSC 3例 (P11-P13)、SCLC/LCNEC 6例 (P14-P19) である。本臨床研究の主要評価項目 (primary endpoint) は、空間マルチオミクス解析によるcSCLCのクローン起源の同定および微小環境特性の解明と設定された。希少疾患のため事前に厳密な sample size calculation は実施していないが、利用可能な全切除標本を網羅した。
多領域WES: ホルマリン固定パラフィン包埋 (FFPE) 切片から、レーザーマイクロダイセクション (LCM: laser capture microdissection) 技術を用いて、SCLC領域50、NSCLC/LCNEC領域40、混合領域5の合計95領域、および対応する正常隣接組織16領域を精密に回収した。DNA抽出後、SureSelect Human All Exon V6を用いてライブラリを調製し、高深度多領域WESを実施した。体細胞変異 (SNV/Indel) の同定にはMuTect (v1.1.4) を、コピー数異常 (CNA) の解析にはGATK (v4.0) を使用した。非同義体細胞変異のパターンに基づき、neighbor-joining法 (ape R package) を用いて系統樹を構築した。
空間トランスクリプトーム (ST) 解析: 6例の凍凍結組織切片に対し、10x Genomics Visium Spatial Gene Expression Slide & Reagent Kitを用いてST解析を実施した。Space Ranger (v1.0.0) でマッピング後、Seurat (v4.1.1) を用いてスポットの統合クラスタリングを行い、SCTransformによる標準化を実施した。スポットの細胞組成デコンボリューションにはSPOTlight (seeded non-negative matrix factorization regression) アルゴリズムを適用し、公開scRNA-seqデータセットを参照として使用した。空間的相関の定量化にはRipley’s K関数を用いた。
単一核RNAシーケンス (snRNA-seq): 12例の凍結組織サンプルから、GEXSCOPE (Singleron single-nucleus RNA library protocol) Single Nucleus RNA Library Kit V2を用いて単一核を抽出し、NovaSeq 6000でシーケンスした。CeleScope (v1.14.1) を用いて遺伝子発現マトリクスを生成し、Seurat (v5) で品質管理 (1000 UMI以上、ミトコンドリア遺伝子10%未満) を行った。悪性細胞の同定およびゲノムCNA負担の算出にはinferCNV (v1.14.1) を使用した。SCLCサブタイプ (A, N, P) の可塑性スコアはAUCell (v1.24.0) を用いて算出し、ternary plotで可視化した。
統計解析: 2群間の連続変数比較には、ノンパラメトリック検定である Wilcoxon rank-sum test または Mann-Whitney U test を適用した。多重比較調整や生存解析におけるハザード比算出には Cox proportional hazards モデルおよび log-rank test を用いた。統計的有意基準は p<0.05 とした。