- 著者: Rong Tang, Liquan Jin, Chen Chen, Yingna Liao, Dingru Li, Yangyi Li, Yueyue Chen, et al.
- Corresponding author: Si Shi, Xianjun Yu (Fudan University Shanghai Cancer Center / Shanghai Pancreatic Cancer Institute)
- 雑誌: Cancer Cell
- 発行年: 2026
- Epub日: 2026-05-19
- Article種別: Original Article
- PMID: 42276047
背景
膵管腺癌(PDAC: pancreatic ductal adenocarcinoma)は依然として最も予後不良の固形癌の一つであり、根治切除後でも5年生存率は約13%、術後中央生存期間は10~30か月にとどまる。これまでのバルクRNA-seq(RNA sequencing)を用いた大規模コホート解析(TCGA-PAAD等)により多数の予後関連遺伝子が同定されてきたが(Liu et al. 2018)、バルクレベルの測定は腫瘍微小環境(TME: tumor microenvironment)内の各細胞型の遺伝子発現を混在させるため、「どの細胞型でどの遺伝子の発現が予後に関連するか」というcell type解像度の情報を提供できないという本質的な限界があった。近年、scRNA-seqをはじめとする単細胞解像度技術がTME内の細胞型別トランスクリプトームプロファイリングを可能にしてきた一方(Roerden et al. 2025)、長期生存データとの統合はほとんど行われておらず、未解明のままであった。snRNA-seqは凍結アーカイブ組織からの解析が可能であり、長期追跡コホートの組み立てを大幅に簡略化できるという特長から、本研究ではPDAC術後長期生存データとsnRNA-seqの統合による予後マップ構築が試みられた。PDACのTMEはがん関連線維芽細胞(CAF: cancer-associated fibroblast)・内皮細胞・免疫細胞・膵外分泌細胞など多様な細胞型から構成される特異な構造を持つが(Bailey et al. 2016)、こうした細胞型特異的な予後関連シグナルはこれまで体系的に解析されておらず、知識の空白が残されていた。
目的
手術切除PDAC患者152例を対象にsnRNA-seqと長期全生存(OS)データを統合し、cell type解像度の予後関連遺伝子(crPRG)を網羅的に同定する。さらに、空間トランスクリプトミクスを用いてcrPRGの空間分布と治療反応性の関係を明らかにし、研究コミュニティ向けのインタラクティブwebプラットフォーム(ctPANDA)を構築する。また、pan-cell-typeのSS-PRGであるPLOD2を標的としたPROTAC(タンパク質分解誘導化合物)の開発とその抗腫瘍効果を実証する。
結果
所見1 - 大規模PDACコホートの構築と7,248個のcell type解像度予後遺伝子(crPRG)同定:
Fudan University Shanghai Cancer Center(FUSCC)手術PDAC患者3,103例のうちn=152例に対してsnRNA-seqを実施し、1,203,189個の単核を非教師なしクラスタリングにより上皮・間質・内皮・免疫・内分泌細胞の5大細胞群、さらに23種の主要細胞型に分類した(Fig. 1B)。全追跡患者の80.3%が死亡エンドポイントに到達しており、TCGA-PAAD(The Cancer Genome Atlas pancreatic adenocarcinoma、48.6%)を大幅に上回る予後解析の精度を確保した(Fig. S1D)。Cell type別平均発現量からCox比例ハザードモデルを用いて7,248個のcell type解像度の予後関連遺伝子(crPRG)を同定し(Table S2A・S2B)、そのうち67.6%は臨床病理学的変数(TNM病期・神経浸潤・リンパ管浸潤)から独立した予後因子であった(Figs. S2B・S2C)。代表的SS-PRGであるPLOD2はPDAC ductal細胞でHR 2.34(高発現群 vs. 低発現群、n=152例)と不良予後と顕著に関連し、短期生存(SS)関連PRG(SS-PRG)は全crPRGの41.9%、長期生存(LS)関連PRG(LS-PRG)は49.5%を占めた。
所見2 - 細胞型依存性の双方向性予後関連遺伝子(bi-PRG)の同定と機能検証:
619遺伝子が細胞型によって相反する予後関連性を示すbi-PRGとして同定された(Fig. 2B・S2D-E)。代表例としてCD44はPDAC細胞・CAFでは不良予後に関連する一方、T細胞では良好予後と関連し、ITGB1もPDAC細胞では不良予後、T細胞では良好予後と関連した。ITGB1 KOによりCAR-T細胞の腫瘍浸潤が有意に減少し、T細胞の抗腫瘍機能維持に必須であることが機能的に実証された(Figs. S2I-S2K)。ADAM17はPDAC細胞では不良予後関連であるがDCでは良好予後と関連し、ADAM17阻害剤処理によりOVA搭載DCのOT-1細胞プライミング能がELISpotおよびCTV追跡アッセイで有意に低下することが示された(Figs. S2M-S2P)。
所見3 - 低酸素-間質増殖軸がpan-cellular短期生存(SS)関連経路として浮上:
50種Hallmark遺伝子セットのうち90%(45/50)が少なくとも1種の細胞型で予後関連性を示した(Fig. 2C)。最も頻度高くSS-PRGに濃縮された経路は低酸素・デスモプラジア軸であり、mTORC1シグナリング・解糖系・HIF-1シグナリング・ECM受容体相互作用・フォーカルアドヒージョンが複数の細胞型で悪性予後と関連していた(Figs. 2D・S3D・S4A)。解糖系フラックス(グルコース→ピルビン酸→乳酸)は腫瘍細胞のみならず間質・免疫細胞でも不良予後と関連し(Fig. S3E)、TCAサイクル代謝(アセチルCoA・コハク酸生成)は良好予後と相関した(Fig. S3F)。LTS(OS>5年)vs. STS(OS<1年)比較では830個の差次発現crPRGが同定され、LTS上昇特徴としてCXCL17(PDAC ductal cellで良好予後、CD8+ T細胞との近接性と関連)、STS上昇特徴としてHLA-G(非古典的免疫抑制MHC分子)・SERPINA1・ANGPT2(毛細管状内皮細胞)が示された(Fig. 2G)。
所見4 - Xenium in situ 5K空間解析による治療反応性とSS-type細胞ニッチ同定:
29例のPDAC組織(未治療・術前化学療法/化学免疫療法)に対してXenium in situ 5K空間トランスクリプトミクスを適用し、3,110,831細胞の5,095遺伝子パネルデータを取得した(Fig. 4A)。19種の細胞ニッチ(CN)を同定し、CN0・CN4(免疫濃縮)はPR症例で増加、CN6(SS-PDAC ductal・SS-CAF・マクロファージ濃縮)はSD比較でPRで有意に減少した(Fig. 4D-F; ANOVA with Tukey’s test)。CN6は術前療法不応例と関連する三重ニッチ(SS-PDAC ductal / SS-Mono/Macro / SS-CAF)として同定され、その中のSS-CAF→SS-PDAC間のコラーゲン・CD99・CXCL・THY1・JAG1-NOTCH2を介したLR相互作用がPR症例で減少していた(Fig. 4G-H)。SS型腫瘍内神経はPR症例においてSS-PDAC ductal細胞により近接しており(paired t検定)、神経微小浸潤が治療抵抗性に寄与する可能性が示された(Figs. 4J-K)。
所見5 - PLOD2のpan-cell-type SS-PRG同定とPROTAC-G3による治療標的化:
全SS-PRGの中でPLOD2とP4HA1が最多の細胞型で不良予後と関連し(Fig. 6A)、PLOD2はPDAC ductal細胞でHR 2.34・CAFでHR 1.87と特に高い値を示した。PLOD2はHIF1Aと正の相関を示し(Spearman r=0.62、p<0.001)、低酸素条件下で腫瘍細胞・非腫瘍細胞を含む複数の細胞型で発現上昇し、免疫組織化学(IHC)コホート(FFPE、n=166例)で高PLOD2発現は腫瘍濃縮領域・間質濃縮領域の両方で不良予後と関連した(Figs. S9G-S9H)。基質スティフネスがPLOD2の発現と安定性を増加させ、これはF-アクチン重合阻害剤LatA・プロテアソーム阻害剤MG132・TRIM21 KO実験によりユビキチン-プロテアソーム経路依存的であることが示された(Figs. S10C-S10H)。1,635,125化合物の仮想スクリーニングと表面プラズモン共鳴(SPR)検証を経て、TRIM21を利用したPROTAC-G3が設計・合成され、Panc-1細胞において濃度依存的なPLOD2分解能が確認された(0.01-10 μM、Fig. 6C)。PDXマウスモデル(n=5)ではPROTAC-G3処理で腫瘍体積・重量が有意に減少し(two-way ANOVA + unpaired t検定)、ストロマ含量も低下した(Figs. 6E-F)。CD34+ヒト化マウスモデル(n=5)においてもPROTAC-G3はCFPAC-1異種移植腫瘍の増殖を有意に抑制し(two-way ANOVA)、体重変化は対照群と同等であった(Fig. S12I)。scRNA-seq解析ではPROTAC-G3処理後にT細胞浸潤増加・CD8+ T細胞の腫瘍反応性スコア上昇・腫瘍細胞の増殖能低下が確認された(Figs. S13C-S13G)。
考察/結論
本研究は、これまでバルクRNA-seqに依存していたPDAC予後研究を、本研究で初めてcell type解像度にまで拡張した点において画期的な意義を持つ。従来のTCGA等バルクデータは各細胞型の遺伝子発現を混在して測定するため、腫瘍細胞での不良予後関連遺伝子がT細胞では良好予後関連である(bi-PRGの概念)という重要な情報が隠蔽されていた。これと異なり、本研究ではcell type解像度の情報によりADAM17やITGB1がbi-PRGであることが示され、従来の阻害剤開発戦略に対して重要な注意点が提示された。先行研究では腫瘍微小環境の複雑性と免疫回避の機序が示されてきたが(Li et al. Nature 2025)、本研究は空間トランスクリプトミクスを用いてSS-type細胞の三重ニッチが治療抵抗性と関連することを示した新規な知見を提供する。
PLOD2はコラーゲン線維架橋を制御するリジルヒドロキシラーゼファミリーとして知られていたが、これまでPDACにおける8細胞型横断的な予後関連性とPROTAC標的としての可能性は報告されていなかった。PROTAC分野では近年様々な標的の分解誘導が試みられてきており(Pike et al. NatRevChem 2026)、本研究ではPLOD2のE3リガーゼパートナーとしてTRIM21を同定し、PROTAC-G3の設計・合成に成功した。ヒト化マウスモデルでの有効性と良好な忍容性は、PROTAC-G3のtranslational potentialを示す proof-of-conceptとして評価できる。
臨床応用の観点では、ctPANDA(https://pancreaticcancer.cn)は公開されており、研究者がcell type別に予後関連遺伝子・タンパク質活性・代謝フラックスを探索できる初のPDAC特化ツールとして機能する。STS関連特徴は術前化学療法の候補患者同定に、LTS関連特徴は根治切除適応患者の選別に活用できる可能性がある。また、LS-PRGをmRNA-LNP技術で直接腫瘍に送達したり、T細胞特異的LS-PRGをCAR-T細胞製造中に過剰発現させる戦略も提唱された。空間トランスクリプトミクス技術は本研究のように治療抵抗性ニッチの解析に有望であり、単細胞空間解析のがん生物学への応用が加速している(Frangieh et al. NatCancer 2026)。
残された課題として、crPRGの因果関係の確立はなされておらず(相関研究に留まる)、単細胞Perturb-seqによる機能解析との統合が今後の検討として必要である。また、本研究は切除可能PDACSに限定されており、転移性PDACSや他の癌腫へのcrPRGの適用可能性はlimitationとして挙げられている。snRNA-seqは免疫サブセットの捕捉効率に偏りがあるため、大規模scRNA-seqコホートとの相補的な活用が今後の検討課題である。
方法
手術切除PDACの凍結腫瘍組織からDounce均質化+遠心分離で核を単離しsnRNA-seqを実施した(FUSCC由来切除可能PDAC患者n=152例)。152例の高品質データを用いてCell Ranger・Seurat・Harmony等の標準パイプラインで処理し、23種の主要細胞型に分類した。Cox比例ハザードモデルで細胞型別平均発現量(正規化データおよびTPMの2種)からcrPRGを同定した。空間解析にはXenium in situ 5K(5,095遺伝子パネル)を29例に適用し、3,110,831細胞を解析した。空間距離定量化はselmaパッケージ、細胞ニッチ解析はCellChat、ligand-receptor解析はNicheNetを使用した。in vitro実験にはPanc-1・CFPAC-1・BxPC-3の各PDAC細胞株を用い、PLOD2候補化合物の同定には仮想スクリーニング(1,635,125化合物)→表面プラズモン共鳴(SPR)結合解析を実施した。PDXモデル(n=5、患者由来腫瘍のNOD-SCIDマウスへの移植)ではPROTAC-G3を腫瘍サイズが一定サイズに達した時点から投与し、2方向ANOVAと対応のないt検定で統計解析した。ヒト化マウスモデル(CD34+ human cord blood stem cell再構成BALB/c-Rag2-/-Il2rγ-/-マウス)ではhuCD45+含有率をベースラインに取り、腫瘍増殖を2方向ANOVAで評価した。統計検定にはlog rank検定(KM曲線)、ANOVAとTukey検定(細胞ニッチ解析)、Pearson相関係数(空間解析)、対応のあるt検定(放射距離解析)を適用した。