- 著者: Ahmed M. Elhossiny, Padma Kadiyala, Jude Ogechukwu Okoye, Harrison L. Hiraki, Megan C. Procario, Thejaswini Giridharan, Hannah R. Watkoske, Mariana Tannus Ruckert, Jiayue Wang, Brian D. Griffith, Alexander W. Bray, Jamie N. Mills, Carlos E. Espinoza, Jörg Zeller, Nicole Peterson, Filip Bednar, Yaqing Zhang, Arvind Rao, Costas A. Lyssiotis, Julianne M. Szczepanski, Jiaqi Shi, Atul Deshpande, Anirban Maitra, Elana J. Fertig, Eileen S. Carpenter, Timothy L. Frankel, Marina Pasca di Magliano
- Corresponding author: Marina Pasca di Magliano (University of Michigan), Timothy L. Frankel (University of Michigan), Eileen S. Carpenter (University of Michigan), Elana J. Fertig (University of Maryland)
- 雑誌: Cancer Discovery
- 発行年: 2026
- Epub日: N/A
- Article種別: Original Article
- PMID: 42165710
背景
膵管腺癌 (PDAC: pancreatic ductal adenocarcinoma) は、極めて予後不良な悪性腫瘍であり、腫瘍微小環境 (TME: tumor microenvironment) が腫瘍組織全体の最大90%を占めるという特徴的な間質豊富な構造を有している。膵癌の主要な前駆病変である膵管内上皮内腫瘍性病変 (PanIN: pancreatic intraepithelial neoplasia) は、顕微鏡レベルの微小な非侵襲性病変であり、臨床的に非侵襲的な手法で検出することが困難であるため、その生物学的特性や微小環境の進化プロセスは長らく未解明のままであった。
先行研究において、Hruban et al. (2019) はPanINが膵癌の主要な前駆病変であることを病理組織学的な総説で解説し、さらに Kanda et al. (2012) は早期段階のPanINの大部分にすでに oncogenic KRAS 遺伝子変異が存在することを示した。しかし、Carpenter et al. (2023) が健常ドナー膵臓30例の解析において60%以上にPanIN病変が存在することを見出した事実は、大多数のPanINが実際には浸潤癌へと進行しないというパラドックスを提示した。
従来のシングルセル解析や空間トランスクリプトミクス研究 (Bell et al., 2024; Fiorini et al., 2025; Shiau et al., 2024) は、主に切除された腫瘍組織を対象としており、癌に隣接する「正常対照」組織を用いていた。しかし、癌隣接組織は腫瘍による炎症やフィールドエフェクト (field defect) の影響を強く受けており、真の健常状態を反映していないという課題があった。すなわち、健常なドナー膵臓における孤発性 (sporadic) PanINから浸潤癌への進行過程で、上皮と周囲の間質成分がどのように協調、あるいは乖離して変化するのかを直接比較したヒト組織研究は決定的に不足していた。特に、マウスモデルでは胚性期からのKRAS活性化などヒトの自然発症プロセスとは異なる発生様式をたどるため、ヒト組織を用いた直接的なエビデンスの欠如が、膵癌の早期介入や予防戦略における最大の知識ギャップ (knowledge gap) となっていた。
目的
健常ドナー膵臓由来の孤発性PanINと、膵管腺癌 (PDAC) 患者由来の浸潤癌組織における上皮成分および間質成分の転写プロファイルを、空間トランスクリプトミクス (spatial transcriptomics) と単一細胞RNAシーケンシング (scRNA-seq: single-cell RNA sequencing) を統合した高解像度解析によって直接比較し、前駆病変から浸潤癌への進行に伴う微小環境の再プログラミング機構と、上皮・間質の非同期的進化プロセスを解明することを目的とする。
結果
上皮ドメインにおけるPanINから浸潤癌への転写的連続体の形成: 空間ドメイン解析により、健常ドナー由来のPanIN上皮は、正常な腺房や膵管ドメインとは明確に区別され、分化型腫瘍である腺管状腫瘍 (GLTumor: glandular tumor) ドメインと転写的に極めて類似した連続的スペクトラムを形成していることが示された (Fig. 2D, G)。Palantir擬時間解析では、PanIN → GLTumor → PDTumor (poorly differentiated tumor: 低分化型腫瘍) へと至る一方向性の明確な進化軌跡が描出され、進行に伴ってKRAS/MAPKシグナルが連続的に上昇した。DEG解析において、PanIN上皮でGLTumorと比較して有意に高発現する遺伝子として、細胞外マトリックス (ECM) リモデリングに関与するカリクレイン関連ペプチダーゼ群である KLK6 (kallikrein-related peptidase 6)、KLK7 (kallikrein-related peptidase 7)、KLK8 (kallikrein-related peptidase 8) や、補体活性化抑制因子である CD55、およびグルコーストランスポーターである SLC2A1 (GLUT1) が同定された (Fig. 3E)。これらの遺伝子は、GLTumorにおいて 2.5-fold 以上の有意な発現上昇を示した (p<0.001)。一方、GLTumorおよびPDTumorでは、GSEAにおいてインターフェロンαおよびγ応答経路が有意に濃縮されていた (p=0.002)。免疫染色による検証では、MAPK活性化指標であるp-ERKおよびJAK/STAT経路の指標であるp-STAT3陽性細胞の比率が、PanINと比較して浸潤癌ドメインで有意に増加していることが確認された (Fig. 3I, K)。
PanIN周囲微小環境の独自性と間質成分の非同期的進化: 各上皮構造から150μm以内のスポットを対象とした空間的近傍解析 (neighborhood analysis) を実施したところ、n=11のドナーサンプルにおけるPanIN周囲の細胞ドメイン組成は、正常な膵管周囲と極めて類似しており、浸潤癌のTMEとは劇的に異なることが明らかになった (Fig. 4A)。PanIN周囲には JCHAIN 陽性および CD27 陽性の形質細胞 (plasma cells) が高頻度に集積している一方、マクロファージはPanIN上皮から離れた位置に排除されていた (Fig. 4B, E)。しかし、浸潤癌組織においてはこれが完全に逆転し、形質細胞は散在化し、免疫抑制性のマクロファージドメインが腫瘍上皮の直近に集積する空間的再配置が生じていた。また、三次リンパ組織 (TLS: tertiary lymphoid structure) ドメインは、n=7のPDACサンプルにおいて腫瘍細胞から地理的に離れた領域に特異的に形成されており、CCL19、CCL21、CCL5 の高発現およびそれらの受容体である CCR7、CCR4 との強い相互作用が検出された (Fig. 4H, I)。この形質細胞のPanIN近傍への選択的集積は、n=6の独立したドナーバリデーションコホートでも再現された。
腫瘍特異的線維芽細胞サブタイプ Fibro2 の同定と臨床的予後相関: scRNA-seqアトラスの階層的デコンボリューションにより、7つの線維芽細胞サブタイプが同定された。そのうち、LRRC15 (leucine-rich repeat containing 15) [High] かつ ACTA2 (α-SMA) [High]、さらにWNTシグナルの転写因子である LEF1 (lymphoid enhancer binding factor 1) を高発現するサブタイプ (Fibro2) は、PDAC症例の腫瘍組織にのみ特異的に出現し、正常膵管やPanINの周囲には一切存在しなかった (Fig. 5B, G)。TCGA-PAADデータセットを用いた生存解析において、この LRRC15 [High] (Fibro2) シグニチャーの高値群は、他の6つの線維芽細胞シグニチャーと異なり、唯一有意に不良な生存予後と相関していた (p=0.003)。さらに、Fibro2シグニチャーは代替活性化 (M2型) マクロファージシグニチャーと極めて高い空間的相関を示した (Fig. 5N)。XeniumおよびCosMxの独立した単一細胞空間データセット (n=9 PDAC) を用いた検証でも、Fibro2スコアが癌細胞の直近で特異的に濃縮され、距離の増加に伴って減衰する現象が確認された (Fig. 6B, D)。
癌細胞によるCAFのWNTシグナル活性化と免疫抑制極性化の誘導: 空間転写解析において、Fibro2ドメインでは LEF1、WNT2、WNT5A が高発現しており、腫瘍ドメインでは WNT7A、WNT7B、WNT10A が高発現していた (Fig. 7A, B)。PDAC患者由来オルガノイド (n=3 lines) とCAF (n=3 lines) の3D共培養実験において、CAFにおける LEF1 の発現が共培養によって 3.2-fold の著しい誘導を示す一方 (p<0.001)、正常線維芽細胞ではこの誘導が最小限にとどまることが示された (Fig. 7F)。また、CAFとの共培養により、腫瘍オルガノイドにおいて基底型 (basal-like) マーカー (KRT17, S100A2, KRT6A) の発現が上昇し、古典型 (classical) マーカー (MUC5AC, TFF1) が減少した。さらに、腫瘍オルガノイドの条件培地をヒトPBMC由来の M0 マクロファージ (n=3 replicates) に曝露した実験では、代替活性化マーカーである SPP1 (secreted phosphoprotein 1) の発現が複数のオルガノイド株において一貫して 4.5-fold 以上の有意な上昇を示すことが確認された (Fig. 7I, p<0.001)。
考察/結論
先行研究との違い: 従来の膵癌空間転写解析研究は、主に癌組織に隣接する「正常対照」組織を用いてPanINを解析していたが、これらは腫瘍による広範な炎症性フィールドエフェクトの影響を受けていた。これに対し、本研究はGift of Life Michiganプログラムを介して得られた「真の健常ドナー」由来の孤発性PanIN組織を直接比較対象として用いた点で、これまでの研究と大きく異なっている。これにより、癌の存在によるバイアスを排除し、純粋な前駆病変段階における微小環境の特性を正確に描き出すことに成功した。
新規性: 本研究は、PanINの上皮成分が転写プロファイルにおいてすでに浸潤癌と連続するスペクトラムを形成している一方で、その周囲の微小環境は正常な膵管周囲と類似した状態を維持しているという「上皮と間質の非同期的進化 (asynchronous evolution)」を本研究で初めて明らかにした。特に、腫瘍特異的な LRRC15 [High]/LEF1 [High] 活性化線維芽細胞 (Fibro2) や、免疫抑制性の SPP1 陽性マクロファージの空間的再配置が、前駆病変段階では存在せず、浸潤癌への進行に伴って初めて出現することを新規に同定した点は極めて新規性が高い。
臨床応用: 本研究の知見は、膵癌の早期介入 (cancer interception) 戦略に新たなパラダイムをもたらす。大多数のPanINが癌に進行しない理由が「間質リプログラミングの欠如」にある可能性が示されたため、上皮の遺伝子変異そのものを標的にするのではなく、Fibro2線維芽細胞の活性化やマクロファージの空間的再配置といった「間質許容性環境」の形成を阻止することが、膵癌発症を予防する画期的な治療標的となる。また、同定された LRRC15 [High] シグニチャーは、PDAC患者の予後不良と有意に相関しており、臨床現場における新たな予後予測バイオマーカーとしての臨床的有用性が期待される。
残された課題: 本研究の限界 (limitation) および今後の検討課題として、解析が横断的なヒト組織サンプルの観察に基づいているため、PanINから癌への進行を同一被験者内で縦断的に追跡できていない点が挙げられる。また、微小なPanIN病変のサンプルサイズの制約から、制御性T細胞や骨髄由来抑制細胞 (MDSC) などの希少な免疫サブセットの空間的動態を完全には解明できていない。今後は、WNTシグナル阻害や線維芽細胞リプログラミング阻止が、実際に前駆病変から浸潤癌への進行を抑制するかを検証する機能的な介入試験や、より大規模な前向きコホートでの検証が必要である。
本研究の成果は、腫瘍のクローン進化 (ITH et al) や、上皮間葉転換 (EMT)、および癌の生物学的特徴である免疫微小環境の再プログラミング (cancer et al. Hallmarks of) の理解を深める上で、極めて重要な学術的価値を有している。
方法
本研究は、Gift of Life Michiganドナープログラムを通じて取得した健常ドナー膵臓 (空間トランスクリプトミクス解析用 n=11 donors、バリデーションコホート用 n=6 donors) および治療前PDAC患者由来の膵癌組織 (n=7 patients)、腫瘍隣接正常膵組織 (n=2 patients) を対象とした。
10x Genomics CytAssist Visiumプラットフォームを用いて空間トランスクリプトミクス解析を実施した。また、健常ドナー (n=24)、隣接正常 (n=3)、PDAC (n=18) から得られた計20万細胞超のscRNA-seqアトラスを統合参照データとして構築した。組織学的アノテーションは、消化器専門病理医3名が独立して実施した。スポット内の細胞組成を推定するため、RCTD (Robust Cell Type Deconvolution) アルゴリズムによるデコンボリューションを行い、BayesSpaceアルゴリズムを用いて17の空間ドメインを定義した。さらに、Palantirアルゴリズムを用いた擬時間 (pseudotime) 解析により上皮の進化軌跡を推定し、Numbatアルゴリズムを用いてscRNA-seqデータからコピー数多型 (CNV: copy number variation) を予測した。
空間的なリガンド-受容体相互作用の解析には LIANA+ (Ligand-Receptor Interaction Analysis+) および CellChat を使用した。独立したバリデーションとして、公開されているXenium単一細胞空間データセット2件 (Bell et al., Fiorini et al.) およびCosMxデータセット (Shiau et al.、164視野、9 PDAC患者) を用い、Proseg unsupervised segmentationアルゴリズムにより細胞境界を再定義して解析した。
機能的検証として、患者由来PDACオルガノイド (n=3 lines) と癌関連線維芽細胞 (CAF: cancer-associated fibroblast) (n=3 lines) またはドナー由来正常線維芽細胞 (n=3 lines) を用いた3D Cultrex共培養実験、およびヒト末梢血単核球 (PBMC) 由来M0マクロファージに対するオルガノイド条件培地曝露実験を実施した。統計解析には、DESeq2を用いた差次的発現遺伝子 (DEG: differential expression gene) 解析、およびMSigDB Hallmark遺伝子セットを用いた遺伝子セット濃縮解析 (GSEA: gene set enrichment analysis) を適用した。多群比較には one-way ANOVA を、2群間比較には Student t-test または Mann-Whitney U test を用いた。