- 著者: Valentin Marteau, Niloofar Nemati, Kristina Handler, Deeksha Raju, Alexander Kirchmair, et al., Dominik Wolf, Isabelle C. Arnold, Stefan Salcher, Zlatko Trajanoski
- Corresponding author: Dominik Wolf, Isabelle C. Arnold, Stefan Salcher, Zlatko Trajanoski (Medical University of Innsbruck, Austria / University of Zurich, Switzerland)
- 雑誌: Cancer Cell
- 発行年: 2026
- Epub日: 2025-12-01
- Article種別: Original Article
- PMID: 41529665
背景
大腸癌 (colorectal cancer, CRC) は世界的に発生率が高く、GLOBOCANで癌関連死因の第2位であり、転移性CRC (metastatic CRC, mCRC) の全生存期間中央値は約30ヶ月に留まる。免疫チェックポイント阻害剤 (immune checkpoint inhibitor, ICI) はミスマッチ修復異常 (deficient mismatch repair, dMMR) / マイクロサテライト不安定性高 (microsatellite instability-high, MSI-H) の約4%のmCRC症例にのみ有効で、microsatellite stable (MSS) 大多数では奏効しない。先行研究では、Galon et al. Science 2006がCRCにおける免疫細胞の種類・密度・局在が予後を予測すること (Immunoscore) を示し、Khaliq et al. Cell 2024はCRCのscRNA-seqアトラスを構築したが10x Genomics platform限定で好中球は捕捉が困難であった。Salcher et al. Nature 2022はNSCLCで腫瘍関連好中球 (tumor-associated neutrophils, TANs) のヘテロジェナイティーと抗原提示能を持つ「hybrid neutrophil」を同定し、Gungabeesoon et al. Nature 2022はinterferon-stimulated neutrophilがICI応答に必須であることを示した。好中球はCRC腫瘍内・間質に豊富に存在しCD8+ T細胞と共局在し抗腫瘍/腫瘍促進両機能を持つが、CRCにおけるサブポピュレーション分類・KRAS変異依存的極性化・空間的組織化は未解明であった。何が足りなかったか：CRC全体規模で好中球を網羅し、抗原提示性TANの存在・予後意義・空間ニッチ・bone marrow-tumor axisを統合した大規模解析は存在せず、KRAS変異がneutrophilを腫瘍促進性へ極性化する分子機序も特定されていなかった。
目的
CRCの包括的シングルセルアトラス (4.27 million cells / 1,670 samples / 76 datasets) を構築して全主要細胞型を高解像度で同定・分類し、特に好中球サブポピュレーションの詳細な解析・KRAS変異と好中球極性化の関係・抗原提示能を持つTANの同定・空間的ニッチの解析・bone marrow axisを統合的に検証し、臨床的意義 (生存予後との関連) を確認すること。
結果
所見1：427万細胞CRCアトラスと4免疫表現型の同定 (best track B: cell counts + cohort scale)：48 studies / 76 datasetsから4,264,946 single cellsを統合 (1,146,196 epithelial / 1,552,108 T/NK / 804,042 B/Plasma / 385,589 myeloid / 339,767 stromal/endothelial、Figure 1D-E、n=487 CRC patients)。Cancer cellは600,740 cellsで5主要サブクラスタ (colonocyte-like FABP1/CEACAM7、BEST4+、Goblet-like MUC2/FCGBP、Crypt-like LGR5/SMOC2、transit-amplifying-like TOP2A/UBE2C) を同定。Leiden clusteringで4 immune phenotype (immune desert [ID]、T cell enriched [T]、B cell enriched [B]、myeloid enriched [M]) に層別化、MSI-H腫瘍はT phenotypeと有意相関 (adjusted OR=2.75、95% CI 1.03-7.33、p=0.043、Figure 2)。
所見2：骨髄球コンセンサスNMF プログラム (best track B: NMF programs + n datasets)：385,589 myeloid cellsに non-negative matrix factorization (NMF) を適用し、14 discovery datasetsで再現される consensus myeloid gene expression programsを同定。細胞活性 (cell cycle / epithelial interaction / metabolic / heat shock / unfolded protein response [UPR] / IFN response / hypoxia) に加え、4 immunomodulatory programs (immunosuppressive response / inflammatory response / B-plasma interaction / T cell interaction) を抽出。Neutrophil consensus programにはIL1A・IL1B (炎症性) と OLR1・VEGFA・PLAU (TAN特性) が含まれた (fold change > 2、FDR < 0.05、Figure 3)。
所見3：KRAS変異依存的TAN polarization (best track B: PDO conditioning + T cell suppression)：Bulk RNA-seq解析でKRAS変異腫瘍はWT腫瘍より好中球浸潤が約1.5-fold増加し、BRAF変異腫瘍では減少 (p<0.01、Figure 4)。KRAS-mutant PDO conditioned mediaで末梢血好中球 (n=6 healthy donors) を6時間曝露するとLOX-1 (OLR1) MFIが約2-fold上昇 (p<0.05)、CD181発現は低下、CD16+ mature neutrophil比率は保持、OLR1 mRNA・IL1A mRNAが各々2.5-fold・3-fold上昇 (qPCR、3 independent experiments)。このTAN-like polarized neutrophilをCD8+ T cellと共培養 (1:4比) するとCD69発現が約40%減少 (p<0.01、n=4 donors)、CD4+ T cellでは差なし (Figure 4)。
所見4：抗原提示HLA-DR+ TAN (TAN2/TAN3) の発見と予後意義 (best track A: cohort II survival)：好中球サブクラスタリングで blood neutrophil (BN1-BN3) と4 TAN clusters (TAN1-TAN4) を同定。TAN2/TAN3はHLA-DRA・HLA-DPA1・CD74高発現の抗原提示性 hybrid neutrophilであり、TAN4はIL1A/IL1B/CXCL1/CXCL2高発現の炎症-腫瘍促進型 (Figure 5)。Flow cytometryでTAN表面HLA-DR蛋白が血液好中球・正常腸粘膜好中球より約3-fold高発現 (p<0.001、n=12 cohort I)。Tissue microarray cohort II (n=210) で HLA-DR+ TAN高密度群 (median splitで層別、上位50% vs 下位50%) は腫瘍コアで有意に改善した5-yr OS (HR 0.62、95% CI 0.41-0.94、p=0.024、Cox回帰、Figure 5)、log-rank p=0.018。
所見5：空間的IL-1 signaling neutrophil-fibroblast ニッチ：Spatial transcriptomics 3.7 million cellsで好中球が間質線維芽細胞 fibroblast S1 / S3 と近接する特異的ニッチを同定。LIANA+細胞間コミュニケーション解析でこのニッチはIL-1 signaling (neutrophil IL1A→fibroblast IL1R1) が優勢、線維芽細胞のCXCL1/CXCL2/CXCL5産生と好中球生存・保持フィードバックループを形成 (Figure 6)。CXCR2 inhibitor SX-682処理 in orthotopic mouse modelでこのニッチが破綻し、tumor growth が約40%減少 (p<0.05、n=8/group、Figure 7)。Orthotopic MC38 mouse modelでbone marrow / blood / tumorのscRNA-seqを時系列解析、tumor-derived G-CSFが骨髄でgranulopoiesisを加速しimmature LY6G+ neutrophilを増加させた (約2.5-fold、Figure 7)。
考察/結論
本研究はCRCにおける最大規模 (427万細胞・650患者) のscRNA-seqアトラスを構築し、好中球サブポピュレーションの体系的分類と機能的多様性を解明した。先行研究との違い：Salcher et al. CancerCell 2022はNSCLCで hybrid neutrophil概念を確立したが、CRCでの存在・分布・予後意義は未確認であった点と異なり、本研究はHLA-DR+ TAN2/TAN3の存在・cohort II tissue microarrayでの良好生存との相関を実証した。またKhaliq et al. Cell 2024の先行CRC scRNA-seqアトラスとは対照的に、本研究はBD Rhapsody platform + 最適化protocolにより好中球を効率的に捕捉でき、myeloid compartmentの解像度が大幅に向上した。新規性：(1) CRCで HLA-DR+抗原提示性TAN (TAN2/TAN3) を本研究で初めて同定し予後意義を確認、(2) KRAS-mutant PDOが好中球をLOX-1+ TAN-likeに極性化しCD8+ T cell抑制を駆動するこれまで報告されていない機序を実証、(3) IL-1 signaling neutrophil-fibroblast ニッチというnovelな空間構造を空間トランスクリプトミクスで定義、(4) orthotopic mouse modelで tumor-bone marrow axisが granulopoiesisを再プログラムする新規な知見を統合した。臨床応用：(1) HLA-DR+ TAN密度を新規 prognostic biomarkerとしてImmunoscoreに追加可能で、MSS CRCの患者層別化にbench-to-bedsideで寄与、(2) KRAS変異腫瘍へのKRAS G12C / G12D 阻害薬 (sotorasib、adagrasib、MRTX1133) と免疫療法の組み合わせ戦略への臨床応用根拠を提供、(3) CXCR2 inhibitor (SX-682、AZD5069等) とICIの併用は本研究で同定されたIL-1 niche破綻という臨床的有用性を持つtranslational戦略となる。残された課題：(1) HLA-DR+ TAN誘導薬の開発と臨床試験 (CRC MSS患者対象、future Phase I)、(2) KRAS-IL-1 axis阻害剤のCRC適応開発、(3) bone marrow tumor-induced granulopoiesisを抑制する分子標的の同定、(4) 公開scRNA-seqデータの統合における batch effect / technical variability のlimitation評価、(5) prospective cohortでのHLA-DR+ TAN予後値の validation、が残された課題として残されている。本研究はCRC TMEを再定義する基盤データセットを提供し、好中球サブタイプベースの精密医療の科学的根拠を確立した。
方法
統合scRNA-seqアトラス構築：公開scRNA-seqデータセット48 studies / 76 datasetsを統合し、4,270,000 cells (1,670 samples、650 patients) を含む大規模CRCアトラスを構築した。患者群内訳：CRC 487例 (median age 58、range 22-91)、polyp 99例、control 64例、COAD 254 / READ 63 / unspecified 170、MSI-H 89 / MSS 268 / unknown 130、TNM stage I 41 / II 121 / III 171 / IV 110。primary tumor 674 / liver metastasis 175 / blood 103 / lymph node 24+18。bioinformatics tool stack: scVI / Harmony統合 (batch correction)、Scanpy v1.9.x、Leiden clustering、CellTypist自動アノテーション、UMAP次元削減、DESeq2 (bulk RNA-seq)、Cell Ranger v6.x → STAR aligner、MSigDB hallmark gene set、QuPath (デジタル病理画像)、Fiji (画像解析)、LIANA+ (細胞間コミュニケーション)、CellChat。Cohort I (n=12 CRC patients)：BD Rhapsody platformで low-mRNA-content cell捕捉に最適化したprotocolで好中球を保護的にsequence、Innsbruck cohortとして統合。spatial transcriptomics: 3.7 million cellsを10x Visium・MERFISH・Stereo-seq・spatial protein profiling (IMC、CO-Detection by indEXing [CODEX]) で 715,000 cells解析。Bulk RNA-seq cohorts: TCGA-CRC (n=573)、AC-ICAM (n=348)、CMC-BSN (n=208) でCIBERSORTx / Shears法 deconvolutionにより好中球浸潤量を推定。Orthotopic mouse model: C57BL/6J wild-typeマウスにMC38 (KRASG12D) 結腸癌細胞を盲腸壁注入、bone marrow / 末梢血 / 脾臓 / 腫瘍を分離しscRNA-seqで「tumor-bone marrow axis」を時系列解析。機能実験: KRAS-mutant patient-derived organoid (PDO) condition media (CM) で好中球を曝露しLOX-1 (OLR1) / CD181 / OLR1 / IL1A mRNA・T cell coculture suppression (CD8/CD4 CD69 expression) を測定。HLA-DR+ TAN spatial定量: cohort II (n=210) のtissue microarray (TMA) でHLA-DR+ TAN密度を腫瘍コア・浸潤端で計数し生存と相関。統計：Cox回帰、Kaplan-Meier、log-rank、Spearman correlation、Fisher exact、ANOVA、Welch t検定をR (v4.x) で実施、Benjamini-Hochberg FDR補正、p<0.05を有意とした。