A multi-dimensional omics framework identifies GPR35 as a driver of M2 macrophage activation and poor prognosis in colorectal cancer

Abstract

Introduction:

Colorectal cancer (CRC) remains a leading cause of global cancer mortality, with therapeutic outcomes heavily reliant on the tumor microenvironment (TME). While immunotherapy has revolutionized treatment for distinct subsets, the mechanisms driving immune evasion in the majority of patients remain elusive.

Methods:

In this study, we constructed a comprehensive single-cell atlas of the CRC TME by integrating multi-cohort scRNA-seq data.

Results:

Through non-negative matrix factorization (NMF), we identified nine intratumoral heterogeneity meta-programs (MPs), among which MP8 was robustly linked to M2 macrophage activation. High-dimensional WGCNA further pinpointed GPR35 as a master regulator within the MP8-associated gene network. Clinical analysis across four independent cohorts validated GPR35 as a significant predictor of poor prognosis. Functionally, GPR35 knockdown in vitro markedly impaired CRC cell proliferation, migration, and invasion. Mechanistically, high GPR35 expression orchestrated an immune-excluded microenvironment characterized by diminished cytotoxic T cell and NK cell recruitment, yet paradoxically elevated immune checkpoint expression. Furthermore, GPR35 expression was negatively correlated with eight established immunotherapy response signatures and associated with aggressive mutational landscapes.

Discussion:

Collectively, our findings identify GPR35 as a novel cancer cell-intrinsic driver of immune evasion and immunotherapy resistance, positioning it as a promising therapeutic target to sensitize "cold" CRC tumors to immune checkpoint blockade.

Introduction

Colorectal cancer (CRC) remains one of the most prevalent and deadly malignancies worldwide, accounting for approximately 10% of all cancer cases and deaths globally. Its high mortality is largely attributed to late-stage diagnosis, metastatic progression, and therapeutic resistance (1). CRC is a molecularly and clinically heterogeneous disease, characterized by distinct genetic, epigenetic, and microenvironmental features that influence disease behavior and treatment outcomes (2, 3). In recent years, immunotherapy, particularly immune checkpoint blockade (ICB), has revolutionized oncology by offering durable responses in several cancer types. However, in CRC, its efficacy is largely restricted to the minority of patients with microsatellite instability-high (MSI-H) or mismatch repair-deficient (dMMR) tumors (4). The majority of patients with microsatellite-stable (MSS) CRC exhibit primary resistance to immunotherapy, underscoring the urgent need to elucidate the mechanisms of immune evasion and identify novel targets to sensitize these tumors to immune-mediated attack. Emerging evidence suggests that GPR35 acts as a critical link between cancer cell-intrinsic signaling and the immune microenvironment, yet its specific role in driving immune evasion in CRC remains to be fully elucidated.

The tumor microenvironment (TME) is a dynamic and complex ecosystem composed of malignant cells, immune infiltrates, cancer-associated fibroblasts, endothelial cells, and extracellular matrix components (5, 6). This intricate network mediates critical processes, including immune surveillance, tumor progression, and therapeutic response (7, 8). The development of single-cell transcriptomic technologies has provided unprecedented resolution for dissecting cellular heterogeneity, identifying rare cell states, and mapping intercellular communication within the TME (9). One key cellular component that bridges tumor microenvironment (TME) remodeling and anti-tumor immunity is the M2 macrophage (10). M2 macrophages are immunosuppressive cells that secrete anti-inflammatory cytokines such as IL-10 and TGF-β, promote angiogenesis and matrix remodeling, and inhibit T-cell function. Tumors with abundant M2 macrophage infiltration are generally more resistant to immunotherapy; however, reprogramming M2 macrophages toward an immunostimulatory phenotype may help overcome immune escape (11). Therefore, identifying key regulators of M2 macrophage polarization within cancer cells could unveil new strategies to reprogram the TME toward an immunogenic state.

G protein-coupled receptor 35 (GPR35) is an orphan receptor with emerging roles in inflammation, metabolism, and cancer (12). It is expressed in various immune cells and epithelial tissues. In cancer, GPR35 has been implicated in cell proliferation, migration, and metastatic potential in certain contexts, yet its precise function and mechanistic role in CRC, particularly in sculpting the immunosuppressive TME, remain largely unexplored. Although its expression and function have been studied in immune cells and certain epithelial contexts, its role in CRC, particularly in shaping the immune landscape, remains poorly understood. Preliminary evidence suggests that GPR35 may influence intracellular signaling pathways related to cell survival, migration, and immune modulation, but systematic studies linking GPR35 to CRC progression and immunotherapy resistance are lacking. Given its signaling potential and expression pattern, we hypothesize that GPR35 may serve as a critical nexus linking cancer cell intrinsic programs to extrinsic immune modulation.

In this study, we hypothesize that cancer cell-intrinsic GPR35 expression modulates the TME and contributes to immunotherapy resistance in CRC. By integrating bulk transcriptomics from multiple cohorts, scRNA-seq data, and in vitro functional validation, we aim to systematically dissect the CRC TME (13). Specifically, we first comprehensively map the cellular and transcriptional landscape at single-cell resolution to identify meta-programs driving intratumoral heterogeneity and M2 macrophage activation. Building on this high-resolution map, we aim to screen and validate key regulator genes associated with patient prognosis. Ultimately, this study focuses on characterizing the immunological and therapeutic implications of the top candidate, GPR35, in driving CRC progression and facilitating immune evasion.

Materials and methods

Data acquisition and preprocessing

Transcriptomic data and corresponding clinical information from CRC patients were obtained from five publicly available cohorts, including TCGA-CRC (14), GSE103479 (15), GSE72968 (16), GSE41258 (17), and GSE72970 (16). We obtained scRNA-seq data from CRC samples of AHCA (5738 cells), GutCellAtlasColon (25104 cells), GutCellAtlasIntestine (41863 cells), TabulaSapiens (8228 cells), and TissueImmuneCellAtlas (453 cells). The R package Seurat was used to process the scRNA-seq data (data normalization, scaling, and integration) (18, 19). Quality control was performed before integration. Cells with unique feature counts (nFeature_RNA) less than 200 or greater than 6000, and cells where >20% of counts originated from mitochondrial genes were filtered out to remove low-quality cells and potential doublets. To correct for batch effects arising from different datasets and sequencing platforms, we applied the Harmony algorithm. The integration performance was assessed visually using UMAP plots to ensure the intermingling of cells from different cohorts while preserving biological clusters. The M2 macrophage signature gene list was obtained from the previous study (20, 21).

Computational analysis

The R package InferCNV was used to identify cancer cells among the TME cells. UMAP was used for two-dimensional visualization (22). A reference set of normal epithelial cells from healthy colon samples was used to infer copy number variation. Cells with elevated large-scale chromosomal aberrations were annotated as cancer cells. To dissect transcriptional heterogeneity within the malignant compartment, non-negative matrix factorization (NMF) was performed on the cancer cell expression matrix. The optimal rank (k = 9) was determined based on cophenetic correlation and the residual sum of squares. The meta-programs (MPs) highlighting intratumor heterogeneity were identified as previously described via the R package GeneNMF (23). Gene modules associated with MP8 were identified via the R package hdWGCNA (24). A soft-thresholding power of six was selected to achieve a scale-free topology fit. The connections between GPR35 and the stromal score, immune score, ESTIMATE score, and tumor purity were determined via the ESTIMATE algorithm (25). The association between GPR35 and nine established immunotherapy response signatures, CYT (26), IFNγ (27), Ayers expIS (27), T-cell-inflamed (27), Roh IS (28), Davoli IS (29), Chemokines (30), RIR (31), and ICBnetIS (6) was explored. The Metascape platform and GSEA of GO terms were used for functional annotation of GPR35 (32). The association between GPR35 and immune cells was predicted through the MCPcounter (33), Porpimol’s study via ssGSEA (34), and TIMER (29). The cancer immune cycles related to GPR35 were quantified (6), with step-specific gene signatures scored via ssGSEA. Somatic mutation landscapes related to GPR35 were visualized and compared using maftools (35).

In vitro validation

Human CRC cell lines, including LOVO and RKO, were maintained in DMEM containing 10% FBS and 1% penicillin/streptomycin under standard culture conditions at 37 °C in a 5% CO₂ humidified atmosphere. To investigate the functional role of GPR35, small interfering RNA (siRNA) targeting GPR35 and a non-targeting control (NC) siRNA were transfected into cells according to the manufacturer’s protocol. Cell proliferation capacity was assessed using the EdU incorporation assay. Migratory and invasive capacities were evaluated using Transwell chambers without or with Matrigel coating, respectively. After 48 hours of incubation, cells that migrated/invaded through the membrane were stained and counted.

Statistical analysis

For comparisons between two groups with non-normal distribution, the Wilcoxon rank-sum test was used. For survival analysis, the log-rank test was employed. P-values from multiple comparisons were adjusted using the Benjamini-Hochberg False Discovery Rate (FDR) method. An adjusted p-value (FDR) < 0.05 was considered statistically significant.

Results

Single-cell deconvolution defines the CRC tumor microenvironment

Unified analysis of scRNA-seq data from multiple public cohorts resolved the major cellular compartments within the CRC TME, including B cells, T cells, Epithelial cells, Stromal cells, Endothelial cells, Innate Lymphoid Cells (ILC), and Myeloid cells (Figure 1A). InferCNV analysis distinguished malignant epithelial cells from the surrounding stroma (Figures 1B, D), while further clustering revealed distinct minor TME subsets, validated by specific marker gene expression (Figure 1C, 2).

Figure 1

Figure 2

NMF identifies intratumoral heterogeneity programs linking MP8 to M2 macrophage activation

Non-negative matrix factorization of cancer cell transcriptomes revealed nine robust MPs representing core biological processes (Figure 3A). Among these, MP8 demonstrated the strongest association with a validated M2 macrophage gene signature (Figures 3B, C), suggesting its role in regulating a potential immunogenic state within tumor cells.

Figure 3

hdWGCNA identifies the MP8-associated green module and nominates GPR35 as a key regulator

To elucidate the gene network underlying MP8, we performed hdWGCNA. Network construction employed a soft power threshold of 6, optimized for scale-free topology (Figure 3D). This yielded multiple co-expression modules (Figure 3E), among which the green module showed the highest positive correlation with the MP8 signature (r = 0.73) (Figure 3F). This module contained 133 genes. Intersection of these genes with 1,206 transcripts significantly upregulated in cancer cells versus all other TME cells yielded 60 overlapping candidates (Figure 4A). Univariate Cox regression analysis prioritized GPR35 as the most significant prognostic hazard gene from this intersection (Figure 4B).

Figure 4

GPR35 expression is a robust prognostic marker across independent CRC cohorts

Kaplan-Meier survival analysis consistently demonstrated that high tumor expression of GPR35 was associated with significantly shorter overall survival across four independent validation cohorts: GSE103479 (Figure 4C), GSE41258 (Figure 4D), GSE72968 (Figure 4E), and GSE72970 (Figure 4F). This establishes GPR35 as a reproducible biomarker of poor prognosis in CRC. To further validate the prognostic value of GPR35, we performed multivariate Cox regression analysis adjusting for age, gender, and tumor stage. GPR35 remained an independent prognostic factor (HR > 1, P < 0.05), indicating its robustness beyond standard clinical variables.

In vitro functional validation confirms the pro-tumorigenic role of GPR35

siRNA-mediated knockdown of GPR35 in LOVO and RKO CRC cell lines was confirmed at the mRNA level (Figures 5A, B). Functional assays revealed that GPR35 depletion significantly impaired cellular proliferation, as measured by EdU incorporation (Figures 5C–F). Furthermore, GPR35 knockdown markedly reduced the migratory (Figures 5G–J) and invasive (Figures 5K–N) capacities of both cell lines in Transwell assays, confirming its direct role in promoting aggressive CRC phenotypes.

Figure 5

GPR35 is linked to an immunosuppressive transcriptional landscape

GSEA revealed that high GPR35 expression was negatively associated with key immune-activating pathways, including “Immune response,” “T cell activation,” and “T cell receptor signaling pathway” (Figure 6A). Metascape pathway annotation corroborated these findings, highlighting the suppression of immune system processes (Figure 6B).

Figure 6

GPR35 expression correlates with an immune-excluded microenvironment and altered cancer-immunity cycle

Analysis of microenvironment scores revealed that high GPR35 expression was inversely correlated with stromal, immune, and ESTIMATE scores but positively correlated with tumor purity (Figure 7A). This suggested an immunologically “cold” tumor niche. Consistent with this, bioinformatics analyses using MCPcounter, ssGSEA, and TIMER algorithms consistently demonstrated that elevated GPR35 levels were associated with significantly reduced inferred infiltration of multiple anti-tumor immune cell types, including CD4 T cells, CD8 T cells, dendritic cells, and neutrophils (Figure 7A). To functionally dissect this immune exclusion, we quantified key steps of the cancer-immunity cycle. This revealed that tumors with high GPR35 expression exhibited significant defects in the recruitment of both T cells and natural killer (NK) cells (Figure 7B). Paradoxically, despite the overall paucity of immune infiltrates, GPR35 expression showed a significant positive correlation with the transcriptional levels of critical immune checkpoint molecules, such as PD-L1 and CTLA-4 (Figure 7C), indicating a co-suppressive microenvironment.

Figure 7

GPR35 expression predicts diminished response to immunotherapy

Evaluation using nine established transcriptional signatures predictive of immunotherapy response demonstrated that high GPR35 expression was associated with significantly lower activity scores for eight of these signatures, including CYT, IFNγ, and T-cell-inflamed (Figure 8). This pattern suggests that GPR35-high tumors exhibit a broadly suppressed adaptive immune microenvironment.

Figure 8

Distinct somatic mutation landscapes associate with GPR35 expression

Analysis of somatic mutations revealed that tumors with high GPR35 expression were characterized by a distinct mutation profile. APC and TP53 were the most frequently mutated genes in this group (Supplementary Figure S1A). A comparative analysis identified several genes with significantly different mutation frequencies between the high- and low-GPR35 groups, including the epigenetic regulator KMT2D, the netrin receptor UNC5B, and the potassium channel KCNH8 (Supplementary Figure S1B, S2).

Discussion

Our integrative multi-omics study provides a comprehensive single-cell atlas of the CRC TME and systematically identifies GPR35 as a novel cancer cell-intrinsic driver of poor prognosis and immune evasion. The robust association between high GPR35 expression and adverse survival across multiple independent cohorts solidifies its clinical relevance as a promising prognostic biomarker.

The functional validation in vitro unambiguously establishes the pro-oncogenic role of GPR35, demonstrating its necessity for the proliferation, migration, and invasion of CRC cells. This positions GPR35 as a potential direct therapeutic target. More critically, our bioinformatic analyses uncover a profound and consistent link between GPR35 and a profoundly immunosuppressive TME. The strong negative correlations with T cell activation pathways, cytotoxic immune cell infiltration (CD4/CD8 T cells, DCs), and essential steps in the cancer-immunity cycle (e.g., T/NK cell recruitment) collectively describe an “immune-excluded” phenotype (36). This exclusion is paradoxically accompanied by increased expression of immune checkpoint molecules, such as PD-L1, suggesting that an adaptive immune resistance mechanism may be concurrently activated, further complicating immune-mediated tumor elimination (37).

The predictive power of GPR35 for immunotherapy resistance is a particularly significant translational finding of our study. The uniformly low activity of eight key immunotherapy response signatures in GPR35-high tumors indicates that GPR35 may orchestrate a broad-spectrum resistance program, potentially independent of microsatellite status. This is crucial because most CRC cases are MSS and currently derive minimal benefit from immunotherapy. GPR35 could thus serve as a biomarker to identify MSS patients unlikely to respond to immune checkpoint inhibitors, thereby sparing them from ineffective and potentially toxic treatments. The isolated elevation of the chemokine in GPR35-high tumors is intriguing and warrants mechanistic investigation. It may represent a maladaptive or immunosuppressive chemokine milieu that fails to productively recruit cytotoxic lymphocytes, or recruits immunosuppressive myeloid populations such as MDSCs or M2-like macrophages. Future studies should profile the chemokine and cytokine secretion profiles of GPR35-overexpressing tumor cells and determine their functional impact on immune cell migration and polarization.

The distinct mutational landscape associated with high GPR35 expression provides additional context. The enrichment of canonical APC (38) and TP53 (39) mutations aligns this subtype with more aggressive, chromosomally unstable CRC pathways. The differential mutation of genes like KMT2D (40), a histone methyltransferase involved in chromatin remodeling, and UNC5B (41), a dependence receptor implicated in apoptosis and angiogenesis, points to potential cooperative genetic events that could synergize with GPR35 signaling to shape both tumor evolution and the immune contexture. For instance, loss-of-function mutations in KMT2D can lead to global alterations in histone methylation, potentially silencing tumor suppressor genes and immune-stimulatory genes, thereby reinforcing the immunosuppressive program driven by GPR35. The co-occurrence of these mutations suggests a possible convergent evolutionary pathway toward immune evasion and therapy resistance.

While our study establishes strong correlative and functional links, several key mechanistic questions remain. The precise molecular mechanism by which GPR35 signaling in cancer cells extrinsically suppresses immune cell infiltration and function is unknown. As an orphan receptor, the relevant endogenous or tumor-derived ligands for GPR35 in the CRC TME need to be identified. Candidates include kynurenic acid, a tryptophan metabolite known to activate GPR35 and accumulate in immunosuppressive environments, and 5-hydroxyindoleacetic acid (5-HIAA), a serotonin metabolite recently shown to promote neutrophil recruitment via GPR35 in inflammatory contexts (23, 36). It is plausible that tumor-derived ligands activate GPR35 in an autocrine manner, leading to the secretion of immunosuppressive factors (e.g., IL-10, TGF-β, VEGF) or the downregulation of chemokines required for T cell homing (e.g., CXCL9, CXCL10). Downstream effector pathways must also be delineated—whether GPR35 exerts its effects primarily through Gαi/o protein-mediated inhibition of cAMP, through β-arrestin-mediated scaffolding and receptor internalization, or through cross-talk with other oncogenic pathways, such as the Wnt/β-catenin pathway, which is frequently dysregulated in APC-mutant CRC.

Future work must validate these interactions in vivo using syngeneic or genetically engineered immunocompetent mouse models of CRC. Such models would enable the dissection of cell-type-specific functions of GPR35, whether its tumor-promoting and immunosuppressive effects are primarily driven by expression in cancer cells or by expression in stromal or immune cells (e.g., myeloid cells). Furthermore, the clinical predictive value of GPR35 should be prospectively tested in cohorts of CRC patients treated with immunotherapy. Assessing GPR35 expression (via RNA sequencing or immunohistochemistry) in pre-treatment tumor biopsies from clinical trials of anti-PD-1/PD-L1 agents would determine its utility as a companion diagnostic.

A primary limitation of this study is the lack of in vivo validation using immunocompetent animal models (42). Future work must employ syngeneic or genetically engineered mouse models of CRC to confirm the causal role of tumor cell-intrinsic GPR35 in shaping an immunosuppressive TME and driving immunotherapy resistance in vivo. Furthermore, the precise downstream signaling pathways (e.g., cAMP, β-arrestin, or cross-talk with Wnt/β-catenin) through which GPR35 exerts its pro-tumorigenic and immunomodulatory effects remain unknown. Systematic approaches such as phosphoproteomics or RNA-seq following GPR35 modulation in relevant models are needed to map these networks. Additionally, the cohorts analyzed were primarily from public databases with limited diversity in race and ethnicity. Multi-center prospective studies incorporating patient samples from varied geographic and ethnic backgrounds, as well as the use of patient-derived organoids (43, 44), are essential to establish the robustness and generalizability of GPR35 as a biomarker. Furthermore, potential selection bias exists as the cohorts analyzed were retrospective and sourced from public databases. To mitigate this, future multi-center prospective studies incorporating patient samples from varied geographic and ethnic backgrounds are essential to generalize our findings.

Our study nominates GPR35 as both a prognostic biomarker and a potential therapeutic target to overcome immunotherapy resistance in MSS CRC. A direct translational application would be to evaluate GPR35 expression in pretreatment biopsies from patients enrolled in ICB trials as a predictive biomarker. From a therapeutic perspective, developing selective GPR35 antagonists or leveraging siRNA/nanoparticle delivery systems to silence GPR35 in tumors could be promising strategies. Combining GPR35 inhibition with anti-PD-1/PD-L1 therapy might synergistically reverse the immune-excluded phenotype. This dual-target approach could potentially convert “cold” tumors into “hot” ones, thereby expanding the reach of immunotherapy in colorectal cancer.

References

• 1

  WangNLinJLiWLyuYJiangYNiZet al. Deep multimodal state-space fusion of endoscopic-radiomic and clinical data for survival prediction in colorectal cancer. NPJ Digit Med. (2025) 8:801. doi: 10.1038/s41746-025-02236-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  AhnHMLeeTGShinHRLeeJYangIJSuhJWet al. Oncologic impact of technical difficulties during the early experience with laparoscopic surgery for colorectal cancer: long-term follow-up results of a prospective cohort study. Curr Probl Surg. (2025) 63:101694. doi: 10.1016/j.cpsurg.2024.101694

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  BillerLHSchragD. Diagnosis and treatment of metastatic colorectal cancer: A review. JAMA. (2021) 325:669–85. doi: 10.1001/jama.2021.0106

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  GaneshKStadlerZKCercekAMendelsohnRBShiaJSegalNHet al. Immunotherapy in colorectal cancer: rationale, challenges and potential. Nat Rev Gastroenterol Hepatol. (2019) 16:361–75. doi: 10.1038/s41575-019-0126-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  WangRZhuangJZhangQWuWYuXZhangHet al. Decoding the metabolic dialogue in the tumor microenvironment: from immune suppression to precision cancer therapies. Exp Hematol Oncol. (2025) 14:99. doi: 10.1186/s40164-025-00689-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  ZhangNZhangHLiSWuWLuoPLiuZet al. Uncovering the predictive and immunomodulatory potential of transient receptor potential melastatin family-related CCNE1 in pan-cancer. Mol Cancer. (2024) 23:258. doi: 10.1186/s12943-024-02169-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  WuYSunRRenSZenginGLiM-Y. Neuronal reshaping of the tumor microenvironment in tumorigenesis and metastasis: bench to clinic. Med Adv. (2025) 3:364–71. doi: 10.1002/med4.70044

  • CrossRef
  • Google Scholar

• 8

  YangJMZhangNLuoTYangMShenWKTanZLet al. TCellSI: A novel method for T cell state assessment and its applications in immune environment prediction. Imeta. (2024) 3:e231. doi: 10.1002/imt2.231

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  HuXHuZZhangHZhangNFengHJiaXet al. Deciphering the tumor-suppressive role of PSMB9 in melanoma through multi-omics and single-cell transcriptome analyses. Cancer Lett. (2024) 581:216466. doi: 10.1016/j.canlet.2023.216466

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  GuoQLiLHouSYuanZLiCZhangWet al. The role of iron in cancer progression. Front Oncol. (2021) 11:778492. doi: 10.3389/fonc.2021.778492

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  NianZDouYShenYLiuJDuXJiangYet al. Interleukin-34-orchestrated tumor-associated macrophage reprogramming is required for tumor immune escape driven by p53 inactivation. Immunity. (2024) 57:2344–2361 e7. doi: 10.1016/j.immuni.2024.08.015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  TakkarSSharmaGKaushalJBAbdullahKMBatraSKSiddiquiJA. From orphan to oncogene: The role of GPR35 in cancer and immune modulation. Cytokine Growth Factor Rev. (2024) 77:56–66. doi: 10.1016/j.cytogfr.2024.03.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  LiYDuanYLiYGuYZhouLXiaoZet al. Cascade loop of ferroptosis induction and immunotherapy based on metal-phenolic networks for combined therapy of colorectal cancer. Explor (Beijing). (2025) 5:20230117. doi: 10.1002/EXP.20230117

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  BarbieDATamayoPBoehmJSKimSYMoodySEDunnIFet al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. (2009) 462:108–12. doi: 10.1038/nature08460

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  AllenWLDunnePDMcDadeSScanlonELoughreyMColemanHet al. Transcriptional subtyping and CD8 immunohistochemistry identifies poor prognosis stage II/III colorectal cancer patients who benefit from adjuvant chemotherapy. JCO Precis Oncol. (2018) 2018:PO.17.00241. doi: 10.1200/PO.17.00241

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  Del RioMMolleviCBibeauFVieNSelvesJEmileJFet al. Molecular subtypes of metastatic colorectal cancer are associated with patient response to irinotecan-based therapies. Eur J Cancer. (2017) 76:68–75. doi: 10.1016/j.ejca.2017.02.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  ShefferMBacolodMDZukOGiardinaSFPincasHBaranyFet al. Association of survival and disease progression with chromosomal instability: a genomic exploration of colorectal cancer. Proc Natl Acad Sci U.S.A. (2009) 106:7131–6. doi: 10.1073/pnas.0902232106

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  SatijaRFarrellJAGennertDSchierAFRegevA. Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. (2015) 33:495–502. doi: 10.1038/nbt.3192

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  NiZZhuPLiuLWangHZhouJWangXet al. Single-cell transcriptomic analysis reveals epithelial-mesenchymal transition and key gene AGRN as a universal programme in gastrointestinal tumours by an artificial intelligence-derived prognostic index. Med Res. (2025). doi: 10.1002/mdr2.70038

  • CrossRef
  • Google Scholar

• 20

  BiKHeMXBakounyZKanodiaANapolitanoSWuJet al. Tumor and immune reprogramming during immunotherapy in advanced renal cell carcinoma. Cancer Cell. (2021) 39:649–661 e5. doi: 10.1016/j.ccell.2021.02.015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  LiuCYuHHuangRLeiTLiXLiuMet al. Radioimmunotherapy-induced intratumoral changes in cervical squamous cell carcinoma at single-cell resolution. Cancer Commun (Lond). (2022) 42:1407–11. doi: 10.1002/cac2.12342

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  LiSZhangNZhangHYangZChengQWeiKet al. Deciphering the role of LGALS2: insights into tertiary lymphoid structure-associated dendritic cell activation and immunotherapeutic potential in breast cancer patients. Mol Cancer. (2024) 23:216. doi: 10.1186/s12943-024-02126-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  GavishATylerMGreenwaldACHoefflinRSimkinDTschernichovskyRet al. Hallmarks of transcriptional intratumour heterogeneity across a thousand tumours. Nature. (2023) 618:598–606. doi: 10.1038/s41586-023-06130-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  MorabitoSReeseFRahimzadehNMiyoshiESwarupV. hdWGCNA identifies co-expression networks in high-dimensional transcriptomics data. Cell Rep Methods. (2023) 3:100498. doi: 10.1016/j.crmeth.2023.100498

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  YoshiharaKShahmoradgoliMMartínezEVegesnaRKimHTorres-GarciaWet al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. (2013) 4:2612. doi: 10.1038/ncomms3612

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  RooneyMSShuklaSAWuCJGetzGHacohenN. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. (2015) 160:48–61. doi: 10.1016/j.cell.2014.12.033

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  AyersMLuncefordJNebozhynMMurphyELobodaAKaufmanDRet al. IFN-gamma-related mRNA profile predicts clinical response to PD-1 blockade. J Clin Invest. (2017) 127:2930–40. doi: 10.1172/JCI91190

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  RohWChenPLReubenASpencerCNPrietoPAMillerJPet al. Integrated molecular analysis of tumor biopsies on sequential CTLA-4 and PD-1 blockade reveals markers of response and resistance. Sci Transl Med. (2017) 9(379):eaah3560. doi: 10.1126/scitranslmed.aah3560

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  DavoliTUnoHWootenECElledgeSJ. Tumor aneuploidy correlates with markers of immune evasion and with reduced response to immunotherapy. Science. (2017) 355(6322):eaaf8399. doi: 10.1126/science.aaf8399

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  MessinaJLFenstermacherDAEschrichSQuXBerglundAELloydMCet al. 12-Chemokine gene signature identifies lymph node-like structures in melanoma: potential for patient selection for immunotherapy? Sci Rep. (2012) 2:765. doi: 10.1038/srep00765

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  Jerby-ArnonLShahPCuocoMSRodmanCSuMJMelmsJCet al. A cancer cell program promotes T cell exclusion and resistance to checkpoint blockade. Cell. (2018) 175:984–997 e24. doi: 10.1016/j.cell.2018.09.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  ZhouYZhouBPacheLChangMKhodabakhshiAHTanaseichukOet al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. (2019) 10:1523. doi: 10.1038/s41467-019-09234-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  BechtEGiraldoNALacroixLButtardBElarouciNPetitprezFet al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. (2016) 17:218. doi: 10.1186/s13059-016-1070-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  CharoentongPFinotelloFAngelovaMMayerCEfremovaMRiederDet al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. (2017) 18:248–62. doi: 10.1016/j.celrep.2016.12.019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  MayakondaALinDCAssenovYPlassCKoefflerHP. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res. (2018) 28:1747–56. doi: 10.1101/gr.239244.118

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  De GiovanniMTamHValetCXuYLooneyMRCysterJG. GPR35 promotes neutrophil recruitment in response to serotonin metabolite 5-HIAA. Cell. (2022) 185:815–830 e19. doi: 10.1016/j.cell.2022.01.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  LeeDChoMKimESeoYChaJH. PD-L1: From cancer immunotherapy to therapeutic implications in multiple disorders. Mol Ther. (2024) 32:4235–55. doi: 10.1016/j.ymthe.2024.09.026

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  ZhangLShayJW. Multiple roles of APC and its therapeutic implications in colorectal cancer. J Natl Cancer Inst. (2017) 109(8):djw332. doi: 10.1093/jnci/djw332

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  VeschiVVeronaFDi BellaSTurdoAGaggianesiMDi FrancoSet al. C1Q(+) TPP1(+) macrophages promote colon cancer progression through SETD8-driven p53 methylation. Mol Cancer. (2025) 24:102. doi: 10.1186/s12943-025-02293-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  PanYHanHHuHWangHSongYHaoYet al. KMT2D deficiency drives lung squamous cell carcinoma and hypersensitivity to RTK-RAS inhibition. Cancer Cell. (2023) 41:88–105 e8. doi: 10.1016/j.ccell.2022.11.015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  PradellaDDeflorianGPezzottaADi MatteoABelloniECampolungoDet al. A ligand-insensitive UNC5B splicing isoform regulates angiogenesis by promoting apoptosis. Nat Commun. (2021) 12:4872. doi: 10.1038/s41467-021-24998-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  LiSYZhangNZhangHWangNDuYYLiHNet al. Deciphering the TCF19/miR-199a-5p/SP1/LOXL2 pathway: Implications for breast cancer metastasis and epithelial-mesenchymal transition. Cancer Lett. (2024) 597:216995. doi: 10.1016/j.canlet.2024.216995

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  GronholmMFeodoroffMAntignaniGMartinsBHamdanFCerulloV. Patient-derived organoids for precision cancer immunotherapy. Cancer Res. (2021) 81:3149–55. doi: 10.1158/0008-5472.CAN-20-4026

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  WuYSunRZenginGRenSLiM-Y. Tumor organoids: breakthroughs in clinical decision making, drug development, and translational advances beyond conventional models. Med Res. (2025) 1:1–15. doi: 10.1002/mdr2.70048

  • CrossRef
  • Google Scholar