# Identification of programmed cell death-related subtypes reveals immune heterogeneity and therapeutic divergence in colon cancer

**Authors:** Peng Xia, Ying Qu, Quanzhong Liu, Mengyan Zhu, Bin Huang, Wei Wu, Kening Li, Lingxiang Wu, Ruohan Zhang, Yingli Lv, Qianghu Wang

PMC · DOI: 10.7150/thno.126314 · Theranostics · 2026-02-18

## TL;DR

This study identifies three colon cancer subtypes based on programmed cell death pathways, revealing differences in immune response and treatment resistance.

## Contribution

A novel classification system for colon cancer subtypes based on programmed cell death pathways and their impact on the tumor microenvironment.

## Key findings

- Three COAD subtypes were identified with distinct immune and therapeutic profiles.
- PCDS3 tumors show poor prognosis and resistance to both chemotherapy and immunotherapy.
- Sunitinib shows selective potency against PCDS3 tumors via the MDK-SDC2 axis.

## Abstract

Therapy resistance remains a critical challenge in colon adenocarcinoma (COAD). The dysregulation of programmed cell death (PCD) pathways significantly influences therapeutic response, but its integrated role in shaping the tumor microenvironment (TME) and driving clinical heterogeneity in COAD is poorly defined.

We established a Programmed Cell Death-related Subtype (PCDS) classification by integrating 12 PCD pathways across transcriptomic data from 1,140 COAD patients using non-negative matrix factorization (NMF). The subtypes were validated in independent RNA-sequencing cohorts. We characterized the genomic, TME, and therapeutic features of each PCDS using multi-omics data analysis, and computational drug repositioning. Molecular docking and in silico drug sensitivity analyses were employed to evaluate candidate drugs.

We identified three robust subtypes, including PCDS1 (immune-activated), PCDS2 (WNT and TP53 signaling activation), and PCDS3 (mesenchymal and T-cell dysfunction/exclusion). PCDS3, enriched with inflammatory cancer-associated fibroblasts (iCAFs), exhibited the poorest prognosis and dual resistance to both chemotherapy and immunotherapy (>80% non-response). Analysis of single-cell and spatial transcriptomics data revealed the activation of MDK-SDC2 ligand-receptor axis between tumor cells and fibroblasts in PCDS3, spatially associated with T-cell dysfunction and exclusion. Computational drug repositioning identified the sunitinib as having selective potency against PCDS3 tumors, showing significantly lower IC50 values and high-affinity binding to SDC2 in molecular docking.

This study defines a novel molecular subtype for COAD, linking PCD dysregulation to distinct TME remodeling and therapeutic outcomes. Targeting the MDK-SDC2 axis with agents such as sunitinib may offer a promising strategy to overcome stromal-mediated immunotherapy resistance in the most lethal PCDS3 tumors.

## Linked entities

- **Genes:** Wnt (protein Wnt-2) [NCBI Gene 100641115], TP53 (tumor protein p53) [NCBI Gene 7157], MDK (midkine) [NCBI Gene 4192], SDC2 (syndecan 2) [NCBI Gene 6383]
- **Chemicals:** sunitinib (PubChem CID 5329102)
- **Diseases:** colon adenocarcinoma (MONDO:0002271), colon cancer (MONDO:0002032)

## Full-text entities

- **Genes:** COL3A1 (collagen type III alpha 1 chain) [NCBI Gene 1281] {aka EDS4A, EDSVASC, PMGEDSV}, SDC2 (syndecan 2) [NCBI Gene 6383] {aka CD362, HSPG, HSPG1, SYND2}, TIGIT (T cell immunoreceptor with Ig and ITIM domains) [NCBI Gene 201633] {aka VSIG9, VSTM3, WUCAM}, CD79B (CD79b molecule) [NCBI Gene 974] {aka AGM6, B29, IGB, Igbeta}, CCL2 (C-C motif chemokine ligand 2) [NCBI Gene 6347] {aka GDCF-2, HC11, HSMCR30, MCAF, MCP-1, MCP1}, KRT8 (keratin 8) [NCBI Gene 3856] {aka CARD2, CK-8, CK8, CYK8, K2C8, K8}, KIT (KIT proto-oncogene, receptor tyrosine kinase) [NCBI Gene 3815] {aka C-Kit, CD117, MASTC, PBT, SCFR}, CPA3 (carboxypeptidase A3) [NCBI Gene 1359] {aka MC-CPA}, RPL17 (ribosomal protein L17) [NCBI Gene 6139] {aka DBA22, L17, PD-1, RPL23, uL22}, APC (APC regulator of Wnt signaling pathway) [NCBI Gene 324] {aka BTPS2, DESMD, DP2, DP2.5, DP3, GS}, CD14 (CD14 molecule) [NCBI Gene 929], TAGLN (transgelin) [NCBI Gene 6876] {aka SM22, SM22-alpha, SMCC, TAGLN1, TGLN, WS3-10}, GNLY (granulysin) [NCBI Gene 10578] {aka D2S69E, LAG-2, LAG2, NKG5, TLA519}, LAG3 (lymphocyte activating 3) [NCBI Gene 3902] {aka CD223}, KRAS (KRAS proto-oncogene, GTPase) [NCBI Gene 3845] {aka 'C-K-RAS, C-K-RAS, CFC2, K-RAS2A, K-RAS2B, K-RAS4A}, HAVCR2 (hepatitis A virus cellular receptor 2) [NCBI Gene 84868] {aka CD366, HAVcr-2, KIM-3, SPTCL, TIM3, TIMD-3}, KRT18 (keratin 18) [NCBI Gene 3875] {aka CK-18, CYK18, K18}, CD79A (CD79a molecule) [NCBI Gene 973] {aka IGA, IGAlpha, MB-1, MB1}, CTNNB1 (catenin beta 1) [NCBI Gene 1499] {aka CTNNB, EVR7, MRD19, NEDSDV, armadillo}, FCGR3A (Fc gamma receptor IIIa) [NCBI Gene 2214] {aka CD16-II, CD16A, FCG3, FCGR3, FCRIIIA, FcGRIIIA}, RBBP4 (RB binding protein 4, chromatin remodeling factor) [NCBI Gene 5928] {aka NURF55, RBAP48, lin-53}, EZR (ezrin) [NCBI Gene 7430] {aka CVIL, CVL, HEL-S-105, VIL2}, COL1A1 (collagen type I alpha 1 chain) [NCBI Gene 1277] {aka CAFYD, EDSARTH1, EDSC, OI1, OI2, OI3}, DCN (decorin) [NCBI Gene 1634] {aka CSCD, DSPG2, PG40, PGII, PGS2, SLRR1B}, PECAM1 (platelet and endothelial cell adhesion molecule 1) [NCBI Gene 5175] {aka CD31, CD31/EndoCAM, GPIIA', PECA1, PECAM-1, endoCAM}, PDCD1LG2 (programmed cell death 1 ligand 2) [NCBI Gene 80380] {aka B7DC, Btdc, CD273, PD-L2, PDCD1L2, PDL2}, TNFRSF17 (TNF receptor superfamily member 17) [NCBI Gene 608] {aka BCM, BCMA, CD269, TNFRSF13A}, STAT3 (signal transducer and activator of transcription 3) [NCBI Gene 6774] {aka ADMIO, ADMIO1, APRF, HIES}, SYNE1 (spectrin repeat containing nuclear envelope protein 1) [NCBI Gene 23345] {aka 8B, AMC3, AMCM, ARCA1, C6orf98, CPG2}, NKG7 (natural killer cell granule protein 7) [NCBI Gene 4818] {aka GIG1, GMP-17, p15-TIA-1}, MZB1 (marginal zone B and B1 cell specific protein) [NCBI Gene 51237] {aka MEDA-7, PACAP, pERp1}, TTN (titin) [NCBI Gene 7273] {aka CMD1G, CMH9, CMPD4, CMYO5, CMYP5, EOMFC}, IL6 (interleukin 6) [NCBI Gene 3569] {aka BSF-2, BSF2, CDF, HGF, HSF, IFN-beta-2}, PRF1 (perforin 1) [NCBI Gene 5551] {aka HPLH2, P1, PFP}, CTLA4 (cytotoxic T-lymphocyte associated protein 4) [NCBI Gene 1493] {aka ALPS5, CD, CD152, CELIAC3, CTLA-4, GRD4}, CD24 (CD24 molecule) [NCBI Gene 100133941] {aka CD24A}, MDK (midkine) [NCBI Gene 4192] {aka ARAP, MK, NEGF2}, CD2 (CD2 molecule) [NCBI Gene 914] {aka LFA-2, SRBC, T11}, S100A8 (S100 calcium binding protein A8) [NCBI Gene 6279] {aka 60B8AG, CAGA, CFAG, CGLA, CP-10, L1Ag}, CEACAM5 (CEA cell adhesion molecule 5) [NCBI Gene 1048] {aka CD66e, CEA}, CXCL14 (C-X-C motif chemokine ligand 14) [NCBI Gene 9547] {aka BMAC, BRAK, KEC, KS1, MIP-2g, MIP2G}, LYZ (lysozyme) [NCBI Gene 4069] {aka AMYLD5, LYZF1, LZM}, ADA2 (adenosine deaminase 2) [NCBI Gene 51816] {aka ADGF, CECR1, IDGFL, PAN, SNEDS, VAIHS}, SMAD4 (SMAD family member 4) [NCBI Gene 4089] {aka DPC4, JIP, MADH4, MYHRS}, LRP1 (LDL receptor related protein 1) [NCBI Gene 4035] {aka A2MR, APOER, APR, CD91, DDH3, IGFBP-3R}, BRAF (B-Raf proto-oncogene, serine/threonine kinase) [NCBI Gene 673] {aka B-RAF1, B-raf, BRAF-1, BRAF1, NS7, RAFB1}, EPCAM (epithelial cell adhesion molecule) [NCBI Gene 4072] {aka Ber-Ep4, BerEp4, DIAR5, EGP-2, EGP314, EGP40}, TPSAB1 (tryptase alpha/beta 1) [NCBI Gene 7177] {aka TPS1, TPS2, TPSB1, TPSB2, Tryptase-2}, PTEN (phosphatase and tensin homolog) [NCBI Gene 5728] {aka 10q23del, BZS, CWS1, DEC, GLM2, MHAM}, GZMB (granzyme B) [NCBI Gene 3002] {aka C11, CCPI, CGL-1, CGL1, CSP-B, CSPB}, TP53 (tumor protein p53) [NCBI Gene 7157] {aka BCC7, BMFS5, LFS1, P53, TRP53}, RET (ret proto-oncogene) [NCBI Gene 5979] {aka CDHF12, CDHR16, HSCR1, MEN2A, MEN2B, MTC1}, ARID1A (AT-rich interaction domain 1A) [NCBI Gene 8289] {aka B120, BAF250, BAF250a, BM029, C1orf4, CSS2}, CLDN5 (claudin 5) [NCBI Gene 7122] {aka AWAL, BEC1, CPETRL1, TMDVCF, TMVCF}, ENG (endoglin) [NCBI Gene 2022] {aka END, HHT1, ORW1}, CD3D (CD3 delta subunit of T-cell receptor complex) [NCBI Gene 915] {aka CD3-DELTA, CD3DELTA, IMD19, T3D}, MIF (macrophage migration inhibitory factor) [NCBI Gene 4282] {aka GIF, GLIF, MMIF}, SDC1 (syndecan 1) [NCBI Gene 6382] {aka CD138, SDC, SYND1, syndecan}, MS4A1 (membrane spanning 4-domains A1) [NCBI Gene 931] {aka B1, Bp35, CD20, CVID5, FMC7, LEU-16}, CD3G (CD3 gamma subunit of T-cell receptor complex) [NCBI Gene 917] {aka CD3-GAMMA, CD3GAMMA, IMD17, T3G}
- **Diseases:** Tumor Immune Dysfunction (MESH:D007154), Cancer (MESH:D009369), MSI-H (MESH:D053842), cytotoxicity (MESH:D064420), T cell dysfunction (MESH:C536780), Inflammatory (MESH:D007249), NMF (MESH:C538347), GIST (MESH:D046152), PCD (MESH:D003643), immunodeficient (MESH:D007153), colon cancer (MESH:D015179), UMAP (MESH:C567162), clear cell renal cell carcinoma (MESH:D002292), COAD (MESH:D003110), CMS (MESH:C535673)
- **Chemicals:** NTP (-), 4,6-Diamidino-2-phenylindole (MESH:C007293), Sunitinib (MESH:D000077210), imatinib (MESH:D000068877)
- **Species:** Mus musculus (house mouse, species) [taxon 10090], Homo sapiens (human, species) [taxon 9606]
- **Cell lines:** PDX — Mus musculus (Mouse), Hybridoma (CVCL_KI95), hu — Homo sapiens (Human), Finite cell line (CVCL_B0BH), PCDS3 — Mus musculus (Mouse), Hybridoma (CVCL_C6V6)

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12964230/full.md

## Figures

8 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12964230/full.md

## References

92 references — full list in the complete paper: https://tomesphere.com/paper/PMC12964230/full.md

---
Source: https://tomesphere.com/paper/PMC12964230