Yang, Yan

Abstract Endocytosis regulates receptor trafficking and signaling, yet its role in colorectal cancer (CRC) remains unclear. We analyzed transcriptomic data from 15,025 CRC tumors to evaluate gene expression in clathrin-mediated endocytosis (CME) and the endosomal sorting complexes required for transport (ESCRT) pathway. Pathway-level signatures were compared across consensus molecular subtypes (CMS) and examined in relation to oncogenic signaling programs, while individual endocytosis-related genes were evaluated for associations with clinical outcomes in an independent randomized trial cohort. Expression of CME and ESCRT pathways varied by CMS, with the highest expression in CMS4 and the lowest in CMS3. Endocytosis signatures correlated with multiple oncogenic signaling pathways, including MAPK, TGF-β, WNT, PI3K, NOTCH, and angiogenesis, suggesting broad links between endocytic activity and tumor biology. In the randomized phase III CALGB/SWOG 80405 trial, low expression of AP2M1 —a core adaptor in the CME complex—was significantly associated with longer overall survival and progression-free survival in patients treated with anti-EGFR agents, but not in those receiving anti-VEGF agents. These findings suggest that endocytic trafficking may be associated with oncogenic signaling and differential therapeutic benefit in patients with metastatic CRC. The endocytosis pathway reflects biologically relevant heterogeneity in CRC and may help guide therapeutic stratification.

3593 Background: CUB domain-containing protein 1 (CDCP1) is a transmembrane protein shown to be over-expressed in colorectal cancer (CRC), with an important role in tumor initiation, growth and metastasis. Through interactions with EGFR, CDCP1 promotes nuclear translocation of the key regulators of Wnt signaling, β-catenin and E-cadherin. High CDCP1 expression has been associated with poor survival outcomes in multiple solid tumors but data in CRC remains limited. We investigated associations between CDCP1 level and survival outcomes in patients with metastatic CRC (mCRC) in the phase III clinical trial CALGB/SWOG 80405 (Alliance). Methods: Data from 433 mCRC patients in the CALGB/SWOG 80405 (Alliance) trial were analyzed. Patients received bevacizumab (n = 226) or cetuximab (n = 207) plus chemotherapy as 1st-line treatment. Tumor RNA was extracted from FFPE samples and processed through the HiSeq 2500 (Illumina) platform. Overall survival (OS) and progression free survival (PFS) curves were evaluated in all treatment subgroups and stratified according to high or low CDCP1 expression. Likelihood ratio tests, hazard ratios (HRs), and 95% confidence intervals (CIs) were calculated using Cox proportional hazards multivariate models, adjusting for age, sex, ECOG performance status, tumor location, number of metastatic sites, KRAS status, consensus molecular subtypes, and treatment arm. Results: Overall, low CDCP1 expression was associated with better OS compared to high expression (HR = 1.38; 95% CI 1.07-1.77; 35.8 vs 25.0 months, p = 0.045). Subgroup analysis re-demonstrated this pattern in males (HR = 1.64; 95% CI 1.64-2.28; 35.9 vs 23.6 months, p = 0.01), patients with MSS tumors (HR = 1.56; 95% CI 1.17-2.08; 36.7 vs 25.1 months, p = 0.0094), and patients treated with cetuximab+FOLFOX (HR = 1.79; 95% CI 1.16-2.77; 40.8 vs 24.0 months, p = 0.03). The effect of CDCP1 expression on PFS was less pronounced. Relative to high expression, low expression was associated with better PFS in males (HR = 1.49; 95% CI 1.10-2.03; 13.1 vs 9.8 months, p = 0.027) and left-sided tumors (HR = 1.50; 95% CI 1.10-2.04; 14.4 vs 10.1 months, p = 0.034). Conclusions: Our findings represent real-world data of the predictive role of CDCP1 in patients with mCRC undergoing anti-EGFR or anti-VEGF treatment. Consistent with studies of other cancers, we show that low CDCP1 expression is predictive of better survival outcomes, with subgroup analysis revealing associations in males, MSS tumors, and left-sided tumors. Among treatment groups, low CDCP1 expression is predictive of better OS in cetuximab+FOLFOX-treated patients, a key finding given CDCP1’s close mechanistic ties to EGFR. Further studies should explore whether these associations are seen in other treatment modalities and whether CDCP1 blockade could represent a viable strategy in patients with high expression.

Background. MicroRNAs are small, single-stranded non-coding RNA molecules that play major role in gene regulation. miRNAs have been reported as both tumor suppressive and oncogenic and are promising candidates for cancer prognosis and treatment. miR-494 is known to promote cell migration and invasion in CRC. Our study aimed to further explore the role of miR-494 modulation (inhibition and overexpression) in CRC. Methods. We evaluated the impact of tumor gene expression of MIR494 on outcomes (overall survival: OS; progression-free survival: PFS) in data for 433 metastatic CRC (mCRC) patients (pts) (bevacizumab [bev] (N = 226) vs cetuximab [cet] (N = 207)) enrolled in randomized phase III CALGB/SWOG 80405 trial obtained from Alliance for Clinical Trials in Oncology. RNA extracted from FFPE tumor samples was sequenced using the HiSeq 2500 (Illumina). Gene expression was divided into high, medium and low. Associations were assessed by log-rank test and multivariable Cox proportional hazards regression model. Treatment-gene interactions were analyzed using the likelihood ratio test. CRC cell lines HCT116 and HT29 were transfected with miR-494 inhibitors and mimics for 48 hours and total RNA was collected, followed by RT-qPCR (miR-494 and target genes BMAL1, PAK1, CDH1, CDH2, SNAI1, TWIST) and clonogenic assays. Results. High miR-494 was significantly associated with worse PFS in overall (median survival: high vs low: 9.3 vs 11.3 months, P<0.025) and cet-treated (median survival high vs low: 9.2 vs 9.6 months, P<0.025) mCRC cohorts, while no significance was observed in bev-treated mCRC. Multivariable Cox model showed significant association between high miR-494 and worse PFS and OS in the cet-treated cohort; however, treatment-gene interaction was not statistically significant. Expression analysis revealed miR-494 was significantly reduced and increased following inhibitor and mimic transfection, respectively. Differential gene expression of the target genes analyzed was observed between the miR-494 inhibitor and mimic-treated HCT116 and HT29 samples. In HCT116 cells, inhibition of miR-494 significantly (P < 0.01) reduced expression of PAK1, while increased the expression of CDH1. The number of colonies formed relatively increased after transfection with miR-494 inhibitor compared to miR-494 mimic. Conclusion. Our findings showed that high miR-494 expression correlated with worse PFS overall and in cet-treated mCRC pts. It would be interesting to check for potential associations with EGFR signaling (direct or indirect). Further studies are warranted for efficient targeting and to gain better understanding of the associated signaling involved in miR-494 modulation in CRC. Support: U10CA180821, U10CA180882; U10CA180888 (SWOG): Genentech. https://acknowledgments.alliancefound.org; P30CA014089, UG1CA180830 [HJL]. Clinicaltrials.gov Identifier: NCT00265850 Citation Format: Pooja Mittal, Joshua Millstein, Shivani Soni, Francesca Battaglin, Karam Ashouri, Sandra Algaze, Priya Jayachandran, Yan Yang, Lesly T. Gonzalez, Unnati H. Shah, Wu Zhang, Fang-Shu Ou, Alan P. Venook, Indrakant K. Singh, Heinz J. Lenz. Effect of miR-494 modulation in colorectal cancer (CRC): Clinical outcome, target expression and cell characteristics [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 2787.

224 Background: Metabolic pathways are reprogrammed during CRC development, leading to increased glycolysis, glutaminolysis, and fatty acid synthesis. Understanding the impact of metabolic reprogramming on treatment outcomes is paramount. Hence, we investigated whether the tumor gene expression of 4 selected major metabolic genes ( PKM, SLC2A1 , SLC16A1 , and CAV1 ) could affect treatment response in pts enrolled in the CALGB/SWOG 80405 trial. Methods: 433 mCRC pts treated with either bevacizumab (bev, n = 226) or cetuximab (cet, n = 207) in combination with first-line chemotherapy were analyzed. RNA was isolated from FFPE tumor samples and sequenced on the HiSeq 2500 (Illumina). Overall survival (OS) and progression-free survival (PFS) were compared between groups of pts categorized by tertiles of gene expression (high [H], medium [M] and low [L]). Logrank P -values provide a non-parametric unadjusted assessment of differences. Likelihood ratio tests (LRT) were computed from multivariate Cox proportional hazards models, adjusting for age, sex, ECOG, tumor side, number of metastatic sites, KRAS , CMS, and treatment. Results: PKM gene expression was associated with cet treatment outcomes with PKM- low tumors showing significantly longer PFS and OS (median PFS: 14.3 vs 9.8 vs 8.1 months, L vs M vs H, P < 0.0001, LRT P = 0.011; median OS: 45.9 vs 31.2 vs 20.9 months, P < 0.0001, LRT P = 0.12). Low SLC2A1 was also associated with longer OS in cet-treated pts (32.4 vs 35.8 vs 25.2 months, P = 0.02, LRT P < 0.0001). Similar results were observed in bev-treated pts where low SLC2A1 expressing tumors had longer OS (37.4 vs 26.1 vs 29 months, P = 0.037, LRT P = 0.069) and a non-significant trend for longer PFS (13.1 vs 11 vs 9.5 months, P = 0.076). CAV1 low gene expression was associated with longer OS in cet-treated pts (37.4 vs 34 vs 21.5 months, P = 0.0048, LRT P = 0.013), while no significant associations were found in bev-treated (interaction LRT P = 0.0009). No significant results were observed for SLC16A1 . Conclusions: The intratumoral expression of metabolic genes was prognostic and predictive in mCRC pts treated with first-line therapy. GLUT1 (encoded by SLC2A1 ) and PKM2 (encoded by PKM ) promote aerobic glycolysis of cancer cells, known as Warburg effect, supporting rapid cell proliferation and survival. CAV1 has been shown to be involved in mitochondrial bioenergetics and fatty acid metabolism and to play either a tumor suppressor or oncogenic role depending on the cancer type and stage. Lower expression of SLC2A1 and PKM was associated with improved survival and increased benefit from cet treatment in our cohort. Low CAV1 expression was also associated with increase survival suggesting a cancer-promoting role in mCRC. Our results suggest that targeting cancer cell metabolism through novel inhibitors of the aforementioned pathways may be a promising therapeutic strategy.

269 Background: Sirtuins (SIRT1-7) are class III histone deacetylases with different roles in CRC as tumor suppressive and oncogenic, mainly influencing proliferation, migration, invasion, DNA damage and chemoresistance. In the present study, we explored the association between tumor gene expression of sirtuins and first-line treatment outcomes in mCRC pts. Methods: We evaluated the impact of the tumor gene expression of sirtuin genes ( SIRT1 -7) on treatment outcomes (overall survival: OS; progression free survival: PFS) in 433 pts enrolled in the randomized, phase III CALGB/SWOG 80405 (NCT00265850) trial. RNA extracted from FFPE tumor samples was sequenced using the HiSeq 2500 (Illumina). Subgroup analyses were based on treatment (bevacizumab [bev] (N = 226) vs. cetuximab [cet] (N = 207)). For each gene, expression was divided into tertiles of high, medium and low. The associations between tumor gene expression and outcomes (PFS and OS) were assessed using log-rank test in univariate analysis. Multivariable Cox proportional hazards regression model was used to compute hazard ratios (HR) and 95% confidence intervals while adjusting for baseline characteristics. Treatment-gene interactions were analyzed using the likelihood ratio test. Results: Of the seven genes analyzed, SIRT3 high was associated with both longer PFS (HR: 0.59 [95% CI: 0.46-0.74], P < 0.001) and OS (HR:0.64 [95% CI: 0.50-0.83], P = 0.002) overall. Treatment-gene interaction testing found significant interactions for SIRT1 and SIRT5 with PFS ( P values: 0.021 and 0.016, respectively) and OS ( P values: 0.038 and 0.014, respectively) with targeted therapies (bev vs. cet). Multivariate analysis showed SIRT1 high associated with longer PFS (HR: 0.62 [95% CI: 0.42-0.89], P = 0.014) and OS (HR: 0.59 [95% CI: 0.40-0.89], P = 0.036) in cet-treated pts; in contrast with shorter observed PFS (HR: 1.38 [95% CI: 0.96-1.99], P = 0.22) and OS (HR: 1.31 [95% CI: 0.90-1.93], P = 0.23) in bev-treated pts. SIRT5 high showed similar results and associated with longer PFS (HR: 0.64 [95% CI: 0.45-0.91], P = 0.039) and OS (HR: 0.57 [95% CI: 0.39-0.83], P = 0.013) in cet-treated pts, with no similar associations in bev-treated pts. Conclusions: Our results highlight a potential role for SIRT3 as prognostic and, SIRT1 and SIRT5 as predictive biomarkers for first-line treatment in mCRC, with differential effects based on targeted agents. Significant interaction favored cet treatment in SIRT1 high and SIRT5 high tumors. However, previous reports suggest an opposite role of SIRT5 that has been linked to chemotherapy and/or cet resistance in CRC via the inactivation of succinate dehydrogenase complex subunit A. Therefore, further studies are warranted to dissect the mechanism behind the interaction of SIRT5 with cetuximab activity and establish the potential of targeting SIRT1/5 in CRC.

254 Background: SLFN11 was recently identified as a dominant determinant of response to DNA-damaging anticancer agents. Several studies showed SLFN11 expression led to better sensitivity to topoisomerase inhibitors and PARP inhibitors. However, clinical evidence in mCRC has been lacking. Methods: Two independent cohorts (CALGB/SWOG 80405 trial and a real-world patient [Caris CODEai] cohort) were included in this study. Patients were stratified into SLFN11 -high and -low groups based on the median SLFN11 gene expression levels. Survival outcomes were compared between SLFN11 -high and SLFN11 -low in irinotecan-treated and oxaliplatin-treated patients. Median time was estimated using Kaplan-Meier method and Cox models were used for covariate adjustment. Endpoints of interest were progression-free survival (PFS) and overall survival (OS). Results: In the CALGB/SWOG 80405 trial cohort, 433 patients were included (110 treated with FOLFIRI-based chemotherapy and 323 treated with FOLFOX-based chemotherapy). In patients treated with FOLFIRI-based chemotherapy, SLFN11 -low patients exhibited numerically better PFS (median PFS, 13.4 vs 11.2 months, log-rank p = 0.03, adjusted hazard ratio [HR] = 0.86, 95% confidence interval [CI] 0.55-1.39, p = 0.50) and OS (median OS, 39.6 vs 30.3 months, log-rank p = 0.02, adjusted HR = 0.73, 95% CI 0.46-1.16, p = 0.19) compared to SLFN11 -high patients. In patients treated with FOLFOX-based chemotherapy, SLFN11 -low patients exhibited numerically better PFS (median PFS 11.2 vs 9.8 months, log-rank p = 0.03, adjusted HR = 0.83, 95% CI 0.65-1.05, p = 0.12) and OS (31.1 vs 26.1 months, log-rank p = 0.07, adjusted HR = 0.84, 95% CI 0.65-1.08, p = 0.16) compared to SLFN11 -high patients. No significant interaction was observed between treatment and SLFN11 expression in terms of PFS ( p = 0.69) and OS ( p = 0.44). In the Caris CODEai cohort, a better OS in SLFN11 -low patients than in SLFN11 -high patients was observed in irinotecan-treated patients ( N = 2906) (median OS, 25.5 vs 21.7 months, HR = 0.90, 95% CI 0.82-1.00, p = 0.048). However, no significant difference in OS between SLFN11 -low and SLFN11-high was observed in oxaliplatin-treated patients ( N = 4428) (median OS, 39.5 vs 35.8 months, HR = 0.98, 95% CI 0.90-1.06, p = 0.56). Conclusions: Low SLFN11 expression showed a potential for better prognosis in mCRC patients receiving irinotecan- and oxaliplatin-containing chemotherapies. These findings are not consistently aligned with preclinical data and clinical evidence in other cancer types previously reported. Further investigation is warranted to better understand the biological role and clinical implication of SLFN11 in mCRC.

3588 Background: The brain-gut axis, a bidirectional network between the central and enteric nervous systems, plays a critical role in modulating the gastrointestinal tract function. Increasing evidence suggests that neuronal signaling molecules can promote gastrointestinal cancers; however, the mechanisms remain unclear. Gamma-aminobutyric acid (GABA) is a major inhibitory neurotransmitter and dysregulated GABA signaling has been found in CRC, affecting both tumor and tumor microenvironment. Here, we investigated the association between tumor GABA pathway gene expression and first-line treatment outcomes in mCRC. Methods: We evaluated the impact of the expression of 25 genes in the GABA pathway ( GABBR1, GABBR2, GABRA1, GABRA2, GABRA4, GABRA5, GABRA6, GABRB1, GABRB3, GABRD, GABRE, GABRG1, GABRG2, GABRG3, GABRQ, GABRR1, GABRR2, ABAT, ALDH5A1, GAD1, GAD2, SLC6A1, SLC6A13, SLC6A11, SLC6A12) on treatment outcomes (overall survival: OS; progression free survival: PFS) of 433 pts enrolled in the randomized, phase III CALGB/SWOG 80405 trial. RNA extracted from FFPE tumor samples was sequenced using the HiSeq 2500 (Illumina). The associations between gene expression (continuous, per standard deviation increase) and outcomes were assessed via the multivariable Cox proportional framework, with hazard ratios (HR) and 95% confidence intervals (CI) computed. The likelihood ratio test was used for interaction tests. False discovery rate (FDR) was controlled at 5% with the Benjamini-Hochberg method. Results: Overall, low GAD1 (encoding for glutamate decarboxylase 1, the enzyme responsible for catalyzing the production of GABA from L-glutamic acid) was associated with longer PFS (HR: 1.26, 95% CI [1.13, 1.41], P = 3.2E-05) and OS (HR 1.27 [1.13, 1.43], P = 5.9E-05) across all treatment arms. High ABAT (encoding for the aminotransferase responsible for GABA catabolism) was also associated with longer OS (HR 0.81 [0.72, 0.91], P = 6E-04). To evaluate the predictive value, interactions with targeted therapies bevacizumab (n = 226) vs cetuximab (n = 207) were tested. We found significant interactions for PFS with ABAT ( P= 0.0017) and OS with SLC6A13 ( P= 0.0014 ).In pts treated with cetuximab, high ABAT expression was significantly associated with longer PFS (HR 0.73 [0.61, 0.87], P = 2.8E-04), and high SLC6A13 was marginally associated with worse OS (HR 1.18 [1.00, 1.39], P = 0.052). We did not observe similar associations in pts treated with bevacizumab. Conclusions: Our results highlight a potential role for GABA signaling in affecting mCRC treatment outcomes, suggesting that GABA pathway genes may be promising novel predictive and prognostic biomarkers and supporting the rationale to investigate GABA as a target for drug development. ClinicalTrials.gov Identifier: NCT00265850

3144 Background: Longitudinal brain structure changes across the lifespan are associated with neurodegenerative disorders. We previously showed that genetic variants in dopamine and Parkinson’s disease genes are predictive and prognostic in mCRC. Indeed, the brain-gut axis is emerging as a critical player in CRC biology. Hence, we investigated whether the tumor gene expression of 3 metabolic genes involved in brain development and aging could affect treatment response in patients enrolled in the CALGB/SWOG 80405 trial. Methods: 433 mCRC pts treated with either bevacizumab (bev, n = 226) or cetuximab (cet, n = 207) in combination with first-line chemotherapy were included in the analysis. RNA was isolated from FFPE tumor samples and sequenced on the HiSeq 2500 (Illumina). Overall survival (OS) and progression-free survival (PFS) were compared between groups of pts categorized by tertiles of expression of 3 genes: APOE, DACH1, and GPR139. Likelihood ratio tests, hazard ratios and 95% confidence intervals were computed from multivariate Cox proportional hazards models, adjusting for age, sex, ethnicity, ECOG PS, tumor location, number of metastatic sites, KRAS, MSI status, and treatment. Logrank P-values describe differences without adjustment for patient characteristics. Results: Both APOE and DACH1 (but not GPR139) met the 0.05 false discovery rate threshold for association with both PFS ( P = .00085, P = .015) and OS ( P = .021, P = .021, respectively) in the full Cox models. In cet-treated pts, APOE-high (H) showed significantly shorter PFS (median 9.2 vs 11.8 vs 13.0 months (mo), high vs low (as reference) logrank P = .0021) and OS (median 21.9 vs 35.8 vs 34.6 mo; logrank P = .0017) compared to APOE-median and -low, respectively. Opposite results were observed for DACH1, where DACH1-H tumors had longer PFS (median 14.5 vs 10.6 vs 8.07 mo; logrank P = .0011) and OS (median 39.4 vs 30.9 vs 21.5 mo; logrank P = .0021) when treated with cet. In bev-treated pts, similar results were observed with APOE-H showing shorter OS (median 25.8 vs 26.3 vs 37.5 mo; logrank P = .0016) and DACH1-H longer OS (median 36.5 vs 28.3 vs 22.4 mo; logrank P = .0017). No significant gene-treatment interaction was observed for either gene for PFS and OS. Conclusions: The tumor expression of metabolic genes associated with changes in the brain structure was prognostic in mCRC pts treated with first-line therapy. Notably, APOE and DACH1 have been recently reported to play a role in CRC where APOE overexpression has been associated with aggressive biological behavior whereas DACH1 has been shown to suppress CRC cell growth and invasion. Our results support this evidence and suggest that further investigation into the role of APOE and DACH1 in CRC biology and into their potential as targets for drug development is warranted. Trial Identifier: NCT00265850.

3523 Background: Evidence suggests neurotransmitter signaling may play a crucial role in tumorigenesis of gastrointestinal (GI) cancers, including migration, proliferation, and apoptosis, as well as angiogenesis and the microenvironment. We developed and applied a strategy to identify upstream regulators of neurotransmitters and their associations with patient outcomes. Methods: Prior literature yielded 85 neurotransmitter signaling genes mapped to RNA-seq gene expression in tumor tissue from 545 patients with mCRC enrolled in the phase III CALGB/SWOG 80405 trial. These samples were also sequenced by FoundationOne CDx for DNA damage. Using Cox proportional hazards regression, we identified a subset of CDx genes associated with overall survival (OS) or progression-free survival (PFS). We extended the Causal Inference Test (CIT) for time-to-event data and identified a subset of the 85 signaling genes consistent with mediation of DNA mutation effects on OS and PFS. Another CIT analysis identified potential upstream neurotransmitter regulators, gene expression mediators of mutation effects on the signaling genes. Top genes were then assessed in whole-transcriptome sequencing (WTS) expression data with Caris MI profile collected from 398 patients with advanced CRC enrolled in the MONSTAR-SCREEN-2 trial. Results: At the 0.05 alpha level, we identified 20 signaling genes with gene expression associated with PFS and/or OS, and 68 genes with evidence of DNA mutations associated with PFS and/or OS. False discovery rate (FDR) < 0.05 was applied to filter the 56,674 expression features on PFS and OS associations, yielding 173 unique genes as possible mediators. The CIT yielded 53 putative mediators at FDR < 0.05, suggesting upstream regulation of neurotransmitter signaling. Five genes were identified as targeting at least 8 signaling genes each, PIGU, ANO9, LY6G6D, and GPR160, and FSCN1, where the first 4 are favorable and only FSCN1 is a poorer prognostic. Following up on these 5 genes in MONSTAR-SCREEN-2 by dichotomizing at the median value of WTS expression, 3 were nominally associated with OS, PIGU ( P = 0.050; HR, 0.59; 95%CI, 0.34-1.00), ANO9( P = 0.029; HR, 0.55; 95%CI, 0.32-0.95), LY6G6D( P = 0.0024; HR, 0.44; 95%CI, 0.25-0.76), and FSCN1 was suggestive ( P = 0.079; HR, 1.62; 95%CI, 0.94-2.77), with direction of effects in agreement with 80405 results for all four. Conclusions: Our novel approach identified previously unknown upstream regulators of neurotransmitter signaling associated with treatment outcomes in mCRC. PIGU, FSCN1, and LY6G6D have been suggested in prior literature as targets for CRC and ANO9 has been suggested as a target for GI cancer. These results bolster the strategy of using mediation analysis with large-scale molecular data to identify new therapeutic targets in CRC. NCT00265850.

Background The C-C motif chemokine receptor 5 (CCR5)/C-C motif chemokine ligand 5 (CCL5) axis plays a major role in colorectal cancer (CRC). We aimed to characterize the molecular features associated with CCR5 / CCL5 expression in CRC and to determine whether CCR5 / CCL5 levels could impact treatment outcomes. Methods 7604 CRCs tested with NextGen Sequencing on DNA and RNA were analyzed. Molecular features were evaluated according to CCR5 and CCL5 tumor gene expression quartiles. The impact on treatment outcomes was assessed in two cohorts, including 6341 real-world patients and 429 patients from the Cancer and Leukemia Group B (CALGB)/SWOG 80405 trial. Results CCR5 / CCL5 expression was higher in right-sided versus left-sided tumors, and positively associated with consensus molecular subtypes 1 and 4. Higher CCR5 / CCL5 expression was associated with higher tumor mutational burden, deficiency in mismatch repair and programmed cell death ligand 1 (PD-L1) levels. Additionally, high CCR5 / CCL5 were associated with higher immune cell infiltration in the tumor microenvironment (TME) of MMR proficient tumors. Ingenuity pathway analysis revealed upregulation of the programmed cell death protein 1 (PD-1)/PD-L1 cancer immunotherapy pathway, phosphatase and tensin homolog (PTEN) and peroxisome proliferator-activated receptors (PPAR) signaling, and cytotoxic T-lymphocyte antigen 4 (CTLA-4) signaling in cytotoxic T lymphocytes, whereas several inflammation-related pathways were downregulated. Low CCR5 / CCL5 expression was associated with increased benefit from cetuximab-FOLFOX treatment in the CALGB/SWOG 80405 trial, where significant treatment interaction was observed with biologic agents and chemotherapy backbone. Conclusions Our data show a strong association between CCR5/CCL5 gene expression and distinct molecular features, gene expression profiles, TME cell infiltration, and treatment benefit in CRC. Targeting the CCR5/CCL5 axis may have clinical applications in selected CRC subgroups and may play a key role in developing and deploying strategies to modulate the immune TME for CRC treatment.