- 著者: Jose Reyes, Isabella Del Priore, Andrea C. Chaikovsky, Nikhita Pasnuri, Ahmed M. Elhossiny, Jin Park, Philipp Weiler, Tobias Krause, Andrew Moorman, Catherine Snopkowski, Meril Takizawa, Cassandra Burdziak, Nalin Ratnayeke, Ignas Masilionis, Yu-Jui Ho, Ronan Chaligne, Paul B. Romesser, Aveline Filliol, Tal Nawy, John P. Morris 4th, Zhen Zhao, Marina Pasca Di Magliano, Direna Alonso-Curbelo, Dana Pe’er, Scott W. Lowe
- Corresponding author: Dana Pe’er (Memorial Sloan Kettering Cancer Center / Howard Hughes Medical Institute); Scott W. Lowe (Memorial Sloan Kettering Cancer Center / Howard Hughes Medical Institute)
- 雑誌: Cell
- 発行年: 2026
- Epub日: 2026-05-14
- Article種別: Original Article
- PMID: 41990751
背景
膵管腺癌(pancreatic ductal adenocarcinoma; PDAC)では Kras 変異が膵管内上皮腫瘍(pancreatic intraepithelial neoplasia; PanIN)などの前悪性病変を誘発し、TP53(約70%の症例で変異/消失)・CDKN2A・SMAD4の不活化を経て悪性化する。しかし、組織学的に正常な組織にもKras変異を持つクローン拡大が高頻度に認められるにもかかわらず、実際にPDACへ進展するのはごく一部であり、組織レベルでの制約がどのように破綻するかは長らく不明であった。
著者らの先行研究では、炎症がKRAS変異と相乗的に作用し、高い可塑性とがん様転写プログラムを持つNes(Nestin)陽性の前駆細胞様(progenitor-like)集団が生じることが示されており、同集団は上昇したクロマチンアクセシビリティと細胞間コミュニケーション能を持つことが報告されていた(Burdziak et al. Science 2023)。また、炎症と遺伝子-環境相互作用が先行してエピゲノムリプログラミングを誘発し、前悪性状態の発生に関わることも示されていた(Alonso et al. Nature 2021)。さらに、自然発生的なp53 LOH(loss of heterozygosity)後に特定のコピー数変化が決定論的順序で蓄積し、悪性化の「タイムスタンプ」となることが判明していた(Baslan et al. Nature 2022)。
しかし、p53・CDKN2A・SMAD4という主要腫瘍抑制プログラムが premalignant 段階でどの細胞集団に・どのように選択的に作用するか、また細胞状態の変化と微小環境リモデリングがどのように連動して良性-悪性移行スイッチを引き起こすかについてはgap in knowledgeが残っていた。特に、前駆細胞様状態の出現と腫瘍形成の間に存在する数ヶ月の時間差がなぜ生じるか、その間に「どのようなバリアが機能しているか」の分子・組織レベルの理解が不足していた。
目的
KP^LOH マウスモデルと単細胞・空間トランスクリプトーム解析を組み合わせ、p53 自然喪失時に生じる premalignant サブポピュレーションの転写状態を捕捉し、KRAS と p53 が前駆細胞様ニッチの形成・維持・崩壊をどのように制御するかを明らかにすること。さらにヒト膵炎組織での相同状態の存在を検証し、がん介入(cancer interception)の治療標的としての生物学的根拠を提供する。
結果
p53 LOH 前駆細胞の捕捉と悪性転換中間体の同定:
pre-tumor KP^LOH マウス(3〜4 ヶ月齢)では p53 LOH 細胞(mKate2+/GFP-)が全 mKate2+ 細胞の 1%〜3% と希少に存在した(Fig 1B-1C)。FACS (fluorescence-activated cell sorting) ソーティング後の scRNA-seq 解析では、p53 proficient premalignant 細胞が腺房-管様化生(acinar-to-ductal metaplasia; ADM; Cpa1/Krt19 (keratin 19) 陽性)・神経内分泌系(Scg5 (secretogranin V)/Chga (chromogranin A) 陽性)・tuft 系(Pou2f3 陽性)・増殖性(Mki67 陽性)・胃腺様 PanIN(Dmbt1/Muc6/Tff1/Tff2/Anxa10 陽性)など多様な転写状態を示した。これらのうち Nes/Msn/Hmga2/Hmga1/Vim 陽性の前駆細胞様(progenitor-like)集団が、拡散距離解析でPDACに最も近い転写距離にあることが同定され(Fig 1F)、良性-悪性移行の中間体として位置づけられた。smFISH では個々の Msn+ 前駆細胞様細胞が腺構造内に散在して出現することが確認され、これらは独立して再発的に生じることが示唆された(Fig 1H-1J)。CNA (copy number alteration) の推定解析では、p53 LOH 細胞の多くが chr11(p53)欠失のみの静寂なゲノムと premalignant 転写状態を示す一方、一部は PDAC と共通する chr4(Cdkn2a)欠失・chr2 増幅などの部分的核型変化を伴う「移行型」であり、悪性化トラジェクトリ上の中間体と考えられた。
発がんと腫瘍抑制プログラムの希少前駆細胞様集団への同時収束:
p53 proficient 前駆細胞様集団は、p53 ターゲット遺伝子(細胞周期停止: Cdkn1a, Ccng1; DNA 修復: Mgmt; アポトーシス: Bbc3, Bax)を他の premalignant 状態と比較して最高レベルで発現しており、p53 LOH 時にこれらが特異的に消失した(Fig 1G)。スプライスジャンクション解析で CDKN2A の p19^ARF と p16^INK4a 両 isoform が同集団で共活性化されることが確認され、さらに TGFβ-SMAD4 依存性遺伝子も選択的に高発現していた(Fig 1L)。すなわち PDAC で高頻度に消失する p53/CDKN2A/SMAD4 の 3 経路すべてが前駆細胞様状態に収束していた。同時に、同集団は RAS シグナル・解糖系(Warburg 代謝)・上皮間葉転換(EMT; epithelial-mesenchymal transition)という悪性化プログラムも選択的に高発現し、senescence 関連転写シグネチャも著明に濃縮されており、oncogene-induced senescence(OIS)に類似した発がん・腫瘍抑制力の同一細胞内拮抗構造が明らかとなった。この選択性は前駆細胞様状態以外の premalignant 上皮集団では観察されなかった。
Xenium 空間トランスクリプトームによる前駆細胞ニッチ段階的構築の解明:
KC マウス(n=10: 傷害後 1〜2 日、n=5: 3 週後)の合計 >3.5 百万細胞を 480 遺伝子パネルで Xenium 空間解析した。DC 解析で gastric-like 状態から progenitor-like 状態への連続軸(G-P DC)が転写変動の主軸として抽出され、Msn → Hmga2 → Vim の段階的誘導が確認された(Fig 2F)。G-P niche pseudotime 枠組み(60 µm 半径)で上皮転写状態を主軸に微小環境の動的変化を定量化したところ、G-P 軸の進行とともに管腔スコアの低下・上皮密度の減少・lesion サイズの拡大・間質浸潤が段階的かつ定量的に一致した(Fig 2H-2K)。ニッチ組成は連続的かつ秩序立った変化を示した: gastric-like 主体のニッチでは Gli1+ 筋線維芽細胞と Maf+ 骨髄系細胞が優勢であったが、前駆細胞様上皮への移行とともに Itgax+ 単球/マクロファージと POSTN+/TNC+ 筋線維芽細胞が段階的に集積した(Fig 3E)。最終的に Ccn1+/Ccn2+(YAP (yes-associated protein) シグナル関連)筋線維芽細胞と免疫抑制性 SPP1+/Arg1+/Il1b+ マクロファージが優勢な浸潤性 PDAC に類似した微小環境が、巨視的腫瘍形成前に形成された(Fig 3H)。創傷治癒応答遺伝子シグネチャも多 compartment で有意に誘導された(Wilcoxon rank-sum 検定; Fig 3G)。このニッチリモデリングはニッチ半径 30〜200 µm で頑健であり、単細胞解析のデータとも一致した。intercellular 通信解析では、前駆細胞ニッチへの移行とともに Il18-Il18rap(上皮-骨髄)・Csf2-Csf2rb(上皮-骨髄)・Tgfb1 シグナル(3 compartment 同時)・Jag1-Notch3(上皮-線維芽)・Pdgfb-Pdgfrb(上皮-線維芽)など多軸の双方向シグナルが段階的に活性化し、自己強化的組織回路が形成されることが示された。
ヒト膵炎組織での前駆細胞様ニッチの保存:
がん非罹患ヒト剖検膵臓 scRNA-seq の再解析で、MSN/HMGA1 陽性かつ KRAS シグナル・腫瘍抑制・senescence 関連プログラムを共活性化する希少 ADM/duct 様細胞クラスターが複数ドナーで同定された(Fig 4A-4B)。KRT17 をプロキシマーカーとし、膵炎患者の組織アレイ(n=67 コア; 非腫瘍性 27 例、PDAC 腫瘍隣接 40 例)を多重免疫蛍光染色で解析した(Fig 4C-4G)。KRT17+ 細胞の出現頻度は間質・免疫細胞分率と有意に正相関した(Spearman’s rank correlation; Fig 4D)。KRT17+ 細胞の周囲では、マウス前駆細胞様ニッチマーカー(HMGA1・TNC・ARG1)・senescence 関連マーカー(p53・p16^INK4a・CRYAB)・有糸分裂シグナル(p-ERK・KI67)が濃縮され、CD8a/GZMB(細胞傷害性 T 細胞マーカー)は顕著に欠乏していた(Fig 4G)。この空間的アーキテクチャはマウス前駆細胞様ニッチと構造的に一致した免疫特権的微小環境を示すものであった。
一過性 KRAS^G12D 阻害による前駆細胞ニッチ崩壊と悪性化遅延:
膵炎誘発後 48 時間の MRTX1133 パルス投与により、HMGA2+ 前駆細胞様細胞は 24-fold 枯渇したが(cycling 細胞は 10-fold、gastric pit 様細胞は 8-fold の枯渇にとどまった)、上皮全体は保たれた(Fig 6B)。MRTX1133 開始 4 時間後には、アポトーシスマーカー cleaved caspase-3(CC3)が前駆細胞様細胞の 44% に誘導されたのに対し、gastric 細胞等では <6% にとどまった。Xenium 解析(n=4 vs n=2、合計 2,686,667 細胞)では、前駆細胞様細胞の枯渇とともに Tnc+ 筋線維芽細胞と Itgax+ マクロファージ/単球も広範に枯渇しニッチが崩壊した(Fig 6C-6F)。6 週齢 KP^LOH マウスへの 72 時間 MRTX1133 パルス後の長期追跡(Fig 6G)では、vehicle 群との比較で有意な PDAC 発症遅延が確認された(log-rank 検定; Fig 6H)。最終的に発症した腫瘍は LGALS4(gastric 様)高発現の古典型 identity を持つ傾向を示した。また、抗 CSF1R 抗体によるマクロファージ枯渇実験(n=3 vs IgG n=4)でも、前駆細胞様集団と大型 progenitor lesion が有意に減少した(Kolmogorov-Smirnov 検定; Fig 5G-5H)。
p53 欠損による Progenitor 2 状態と免疫特権ニッチ形成:
KC^shp53 マウスでは傷害後 3 週時点で HMGA2+ 前駆細胞様集団が自然消退せず持続していた(KC^shCtrl では消退; Fig 7B)。scRNA-seq DC 解析では p53 欠損が G-P DC 軸上での上皮細胞進行を促し、Vim/VIM 高発現・インターフェロンシグナル・YAP シグナル・解糖系の亢進を特徴とする「Progenitor 2」状態を生じた(Fig 7C-7D)。Xenium 解析(KC^shp53 n=9、KC^shCtrl n=5、合計 4,611,972 細胞)では、p53 欠損組織で progenitor niche domains が焦点的小領域から連続的大領域へ拡大し、Ccn1+/Ccn2+ CAF(cancer-associated fibroblast)および Itgax+/Cd274^high(PD-L1^high)マクロファージが大量蓄積した(Fig 7H)。PD-L1^high マクロファージは progenitor-like 細胞が一定閾値を超えると急激に蓄積し、MHC (major histocompatibility complex) class II 構成要素を下方調節した免疫抑制状態を示した。sc multiome 解析(n=9,081 細胞)の統合により、前駆細胞様状態で p53 モチーフアクセシビリティと p53 依存的発現を両方示す 138 遺伝子が同定され、うち 84 遺伝子は in vitro p53 再活性化でも上方調節、54 遺伝子は前駆細胞様状態に特異的であった。これにより p53 が細胞周期停止・アポトーシス・senescence・DNA 修復・上皮接着の多軸プログラムを通じて前駆細胞様状態の解消を統括することが明らかとなった。
考察/結論
本研究は、PDAC の良性-悪性移行において、Kras 変異・p53/CDKN2A/SMAD4 腫瘍抑制・細胞状態(progenitor-like plasticity)・微小環境再構築(self-reinforcing niche)が単一の希少前駆細胞様サブポピュレーションに収束するという統合モデルを空間的に実証したものである。Progenitor-like 細胞は「腫瘍駆動と腫瘍抑制が拮抗する場」として機能し、p53 が機能している限りはニッチ形成と悪性化が制限される。
既存研究との相違: これまでの研究では、p53・CDKN2A・SMAD4 は膵がん進行を抑制するバリアとして前悪性上皮全体にわたって働くと考えられてきた。対照的に本研究は、これら 3 経路が KRAS 活性の高い前駆細胞様サブポピュレーションにのみ選択的に集中し、他の premalignant 状態(gastric・ADM・neuroendocrine など)はほとんど影響を受けないことを示した。この既報のパラダイムとは異なる「選択的収束」モデルは、p53/p16^INK4a が人工多能性幹細胞(induced pluripotent stem cell; iPSC)誘導を抑制するという長年の知見と整合し、p53 が「可塑性の守護者(guardian of plasticity)」として機能するという新概念を支持する。
本研究で新規に明らかにした点: 本研究で初めて、単一前駆細胞様集団への腫瘍抑制と発がん力の同時収束を空間トランスクリプトームで可視化し、巨視的 PDAC 形成前に self-reinforcing cancer-like niche が段階的に自律形成されることを新規に実証した。G-P niche pseudotime 法という新規の計算フレームワークにより、静的組織スナップショットから動的ニッチ形成の軌跡を再構築することに成功した点も本研究の独自の貢献である。さらに、がん非罹患ヒト膵炎組織にも KRT17+ バサル様細胞を中心とした類似ニッチが保存されていることを新規に示し、マウスモデルで得られた知見のヒトへの外挿可能性を実証した。
臨床的意義: 一過性の KRAS^G12D 選択的阻害(MRTX1133)が premalignant 段階でニッチを崩壊させ悪性化を有意に遅延できるという知見は、premalignant interception(前悪性介入)の生物学的根拠として重要な臨床的意義を持つ。bench-to-bedside 観点から、高リスク集団(家族性 PDAC・慢性膵炎患者)での KRAS 阻害剤や CSF1R 阻害剤を用いたニッチ指向的早期介入戦略が展望される。また、TP53 変異を持つ難治性がんでは、p53 機能の生化学的回復より、p53 が自然に排除する前駆細胞様状態とニッチ生態系を直接標的とする方が実現可能な治療アプローチとなる可能性がある。実際に KRAS 阻害剤は進行 PDAC でも basal/mesenchymal 様集団を優先的に排除し免疫療法への感受性を改善する報告があり、本研究の premalignant 結果はこれと一貫した臨床含意を持つ。
残された課題: 今後の検討として、(1) ヒト PanIN/IPMN (intraductal papillary mucinous neoplasm) 組織での前駆細胞様ニッチの大規模検証、(2) TP53 ミスセンス変異(特に gain-of-function 変異)を有する場合の前駆細胞様状態への修飾効果の解明、(3) 前駆細胞様状態の解消機序(アポトーシス・senescence 後免疫排除・分化)の単細胞レベルでの fate-mapping 追跡、(4) 肺・大腸など他臓器でのprogenitor-niche 共通性の検証が残された課題として挙げられる。limitation として、使用した p53 欠損モデルはノックアウト/shRNA ノックダウンであり実際の TP53 ミスセンス変異とは性質が異なること、MRTX1133 や抗 CSF1R 介入が前駆細胞ニッチ以外の上皮集団にも影響を及ぼすため、front-cell-circuit の厳密な分離には次世代摂動戦略が必要なことが挙げられる。
方法
動物モデル: KP^LOH モデル(Kras^G12D, Trp53^flox/+, Ptf1a-Cre (pancreatic transcription factor 1a-Cre), mKate2/GFP 蛍光レポーター)で p53 LOH 細胞(mKate2+/GFP-)を生体追跡。多アレル mESC(mouse embryonic stem cell)由来の同期 KC マウス(LSL-Kras^G12D; Ptf1a-Cre)に caerulein 腹腔内投与で膵炎を誘発し、2日後(n=10 マウス)・3週後(n=5 マウス)に組織採取。KC^shp53 / KC^shCtrl モデルはドキシサイクリン誘導型 p53 shRNA を用い1週間誘導後に膵炎処置。KRAS^G12D 選択的阻害には MRTX1133 を 48〜72 時間パルス投与、マクロファージ枯渇には抗 CSF1R (colony-stimulating factor 1 receptor) 抗体(n=3)または IgG コントロール(n=4)を 4 週間投与した。
単細胞解析: pre-tumor KP^LOH マウスからmKate2+/GFP- および mKate2+/GFP+ 細胞を FACS ソーティングし scRNA-seq を実施。傷害後 1 日の KC マウス上皮細胞に対し scRNA+ATAC マルチオーム解析(n=9,081 細胞)を行い、p53 応答エレメントのクロマチンアクセシビリティと転写依存性を統合評価した。拡散成分(diffusion component; DC)解析で gastric-like–progenitor-like 軸(G-P DC)を抽出し、inferCNV (copy number variation inference algorithm) でコピー数変化を推定した。
Xenium 空間トランスクリプトーム: 事前の scRNA-seq 知見に基づき設計した 480 遺伝子パネルを使用し、KC マウス組織の合計 >3.5 百万細胞を Xenium In Situ プラットフォームで解析。60 µm 半径のニッチ定義によるG-P niche pseudotime 枠組みを構築し、上皮転写状態を軸に微小環境変化を定量化した。MRTX1133 投与(n=4)対 vehicle(n=2)の合計 2,686,667 細胞、KC^shp53(n=9)対 KC^shCtrl(n=5)の合計 4,611,972 細胞を解析した。smFISH (single-molecule fluorescence in situ hybridization) で Msn・p53 ターゲット遺伝子を空間的に可視化し、共局在を確認した。
ヒト組織解析: がん非罹患剖検膵臓の公開 scRNA-seq データを再解析。膵炎患者の組織アレイ(67 コア; 非腫瘍性 27 例、PDAC 腫瘍隣接 40 例)を多重免疫蛍光染色で解析し、KRT17 をプロキシマーカーとして前駆細胞様状態を同定、KRT17+ 細胞中心の 200 µm 同心円放射状ビン解析を実施。
統計: 群間比較は two-tailed Wilcoxon rank-sum 検定、ペアデータは paired t 検定、KRT17+ 細胞と間質・免疫細胞の相関は Spearman’s rank correlation、生存曲線は log-rank 検定、progenitor 細胞クラスターサイズ比較は Kolmogorov-Smirnov 検定を使用した。