Innocenti, Federico

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.

198 Background: The Hedgehog (Hh) signaling pathway regulates tissue patterning, growth and differentiation, including within gastrointestinal tract. Mammalian homologs include Sonic (SHH), Indian (IHH) and Desert (DHH) Hedgehog, which signal through PTCH and SMO and the GLI transcription factors. In the digestive tract, both SHH and IHH are expressed, but the functional role of Hh signaling in mCRC remains incompletely defined. Methods: We analyzed 433 mCRC patients treated with either bevacizumab (bev = 226) or cetuximab (cet = 207) plus first line chemotherapy in CALGB/SWOG 80405 (NCT00265850). RNA expression (Illumina HiSeq 2500) of 9 Hh pathway genes (SHH, IHH, DHH, PTCH1, PTCH2, SMO, GLI1, GLI2, GLI3) was evaluated. Median overall survival (mOS) and progression free survival (mPFS) in months (mo) were compared across expression tertiles (low [L], medium [M], high [H]). Likelihood ratio tests, hazard ratios and 95% confidence intervals were computed from multivariable Cox proportional hazards models, adjusting for age, sex, ethnicity, ECOG PS, tumor location, number of metastatic sites, KRAS , microsatellites instability status and treatment. Prognostic effects were modeled as marginal terms; predictive effects were modeled as multiplicative interactions. Results: Overall, IHH and GLI2 demonstrated prognostic value for OS after adjusting for multiple testing (Benjamini-Hochberg False Discovery Rate (FDR) < 0.05). High IHH expression (H group, n = 143) correlated with improved mOS (H: 36.0 mo vs M: 32.8 mo vs L: 23.5 mo; p= 0.031), whereas low GLI2 expression (L group, n = 145) was associated with longer mOS (H: 25.6 mo vs M: 33.9 mo vs L: 33.7 mo; p= 0.0016). Moreover, IHH showed predictive value for cet benefit in both OS (mOS H: 38.1 mo vs M: 34.0 mo vs L: 18.4 mo; p= 0.0032) and PFS (mPFS H: 12.6 mo vs M: 12.9 mo vs L: 8.8 mo; p= 0.0025) with significant treatment interaction for OS and PFS (p= 0.0046 and p= 0.0049, respectively; FDR< 0.05). No significant associations were found in bev treated pts. No other Hh genes showed significant prognostic or predictive effects after FDR correction. Conclusions: These results highlight a dual role of Hh pathway in mCRC. Consistent with prior preclinical evidence, reduced IHH expression appears to mark poor prognosis, likely reflecting the loss of IHH-mediated restraint on Wnt signaling, which drives proliferation and dedifferentiation. Conversely, low GLI2, a downstream transcription factor, was linked with favorable outcomes, underscoring context-specific effects of Hh signaling. The strong predictive association of IHH with cet efficacy supports its potential as a biomarker in RAS wild type mCRC. These results suggest that Hh pathway modulation, particularly IHH signaling, may influence disease biology and treatment response. Validation in independent cohorts and mechanistic studies is warranted to clarify implications.

179 Background: Neuroligins are a family of postsynaptic cell adhesion molecules that play a crucial role in synapse formation and function. Recent studies suggest Neuroligin family genes (NLGN1,2,3,4X) may promote cancer cell growth and metastasis in various cancers, including colorectal cancer. Here, we investigated the association between neuroligin family gene expression and first-line treatment outcomes in mCRC. Methods: We evaluated the expression of 4 genes in the neuroligin family (NLGN1, NLGN2, NLGN3 and NLGN4X) on treatment outcomes of 433 pts with mCRC treated with either bevacizumab (Bev, n = 226) or cetuximab (Cet, n = 207) in combination with first-line chemotherapy within the CALGB/SWOG 80405 trial. mRNA 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 low expression (Q1) (≤25 percentile), medium expression (Q2-3) (25th-75th percentile) and high expression (Q4)(>75th percentile) of NLGN family genes. The associations between gene expression and clinical outcomes were assessed via the multivariable Cox model and plotted with Kaplan-Meier curve. Treatment-gene interactions were evaluated using the likelihood ratio test. False discovery rate (FDR) was controlled through the Benjamini-Hochberg procedure. Results: High NLGN2 expression was associated with shorter PFS (median 10.0 mos, 95% CI (8.3, 11.3)) and OS (median 25.2mos, 95%CI (20.9, 30.0)) compared to low NLGN2 expression (PFS: median 11.2mos, 95% CI (9.8,14.3))(HR 1.41, 95% CI (1.06, 1.89)), P = 0.029). (OS: (median 31.8mos, 95%CI (26.8, 37.5)( HR 1.49, 95% CI (1.10, 2.01) P =0.0095) in multivariate analysis. Patients treated in the cetuximab arm(n=207), high HLGN2 expression significantly associated with shorter PFS (median 8.0 vs 13.4, HR 2.02, 95% CI (1.35, 3.03), P = 0.0021) and OS (median19.6 vs 36.3, HR 2.30, 95% CI (1.50, 3.53), P =8.3e-05). High NLGN3 expression was associated with shorter OS (23.9 vs 36.0 mos, P =0.0083) in cetuximab treatment arms, no significance found in PFS. We found significant treatment interaction for OS with NLGN2 and NLGN3 favoring cetuximab. In pts treated with cetuximab, High NLGN2 expression was significantly associated with shorter OS (median 29.1vs 36.3) (p=0.024) compared to pts treated with bevacizumab. For NLGN3, pts treated with cetuximab, high NLGN3 expression was significantly associated with shorter OS (median 25.0 vs 36.0) (p=0.018) compared to pts treated with bevacizumab. No significant association was found in NLGN1 and NLGN4X. Conclusions: Our results suggest NLGN2/3 gene expression may predict mCRC treatment outcomes. Potential therapeutic drug target NLGN in combination with anti-EGFR treatment in mCRC is warranted.

3145 Background: SnoRNAs are non-coding RNAs that primarily guide the chemical modification of ribosomal RNA. Emerging evidence suggests snoRNAs play critical roles in cancer, including CRC, by regulating cell proliferation, apoptosis and tumor progression. Aberrant snoRNA expression has been linked to CRC development and poor prognosis, offering potential as diagnostic biomarkers and therapeutic targets. We investigated whether the tumor expression levels of 3 types of snoRNAs (SCARNA, SNORA, SNORD) affect treatment response in pts enrolled in CALGB/SWOG 80405 (NCT00265850). Methods: 433 mCRC pts treated with 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). 422 snoRNAs were evaluated (23 SCARNA, 140 SNORA, 259 SNORD). Overall survival (OS) and progression-free survival (PFS) were compared across tertiles of gene expression (high [H], medium [M], low [L]) by multivariable Cox proportional hazards models adjusted for age, sex, ECOG performance status, tumor side, number of metastatic sites, KRAS , CMS subtypes, and treatment. Interaction tests for the predictive effect (bev vs cet) were performed. P -values were corrected for multiple testing using the Benjamini-Hochberg approach ( q < 0.05). Results: Only SCARNA21 achieved statistical significance for OS after false discovery rate (FDR) adjustment. Tumors with H levels of SCARNA21 had shorter survival compared to tumors with M or L expression (median OS 24.4 vs 32.4 vs 33.9 months, respectively; P = 0.0015, q = 0.033), independent of treatment. No snoRNAs achieved FDR significance for PFS. However, several snoRNAs showed significant treatment interactions with biologic agents. SCARNA6-H and SCARNA5-H tumors had longer PFS and OS when treated with cet, but shorter PFS and OS when treated with bev, compared to the M and L expression groups (PFS interaction q = 0.0067 and 0.045, respectively; OS interaction q = 0.022 and 0.034, respectively). SCARNA7 also showed significant treatment interaction for OS, favoring cet in the H expression group ( q = 0.018); the opposite was observed for SNORA63B, SNORA63D, SNORA35B, and SNORA36C ( q = 0.027, 0.027, 0.027 and 0.047, respectively). No significant results were observed for any of the tested SNORDs. Conclusions: SnoRNAs dysregulation affects key pathways such as cell cycle control and immune evasion, making them promising players in CRC biology. Our study highlights the prognostic and predictive potential of specific snoRNAs in mCRC. Notably, high SCARNA21 expression was linked to shorter OS, while SCARNA5 and 6 showed predictive value for treatment response, indicating their potential for guiding treatment decisions. Further validation is needed to confirm these findings, and mechanistic studies are warranted.

3128 Background: LncRNAs have emerged as key regulators of cancer progression and therapeutic responses. In CRC, several lncRNAs have been implicated in modulating tumor growth, metastasis and treatment resistance by interacting with critical tumorigenic pathways. Here, we investigate the potential prognostic and predictive value of lncRNA expression in patients (pts) with mCRC enrolled in CALGB/SWOG 80405 (NCT00265850) trial. Methods: We analyzed 433 mCRC pts treated with either bevacizumab (bev, n = 226) or cetuximab (cet, n = 207) plus first-line chemotherapy. Tumor RNA expression (Illumina HiSeq 2500) of 13 candidate lncRNAs (SNHG11, HOTAIR, FGF14-AS2, H19, YWHAE, NEAT1, MIR100HG, UCA1, LINC00973, SLCO4A1-AS1, POU5F1P4, MALAT1, HCG18) was explored. Median overall survival (mOS) and progression-free survival (mPFS) in months (mo) were compared between pts grouped by tertiles (low [L] vs medium [M] vs high [H]) of gene expression. Likelihood ratio tests, hazard ratios and 95% confidence intervals were computed from multivariable Cox proportional hazards models, adjusting for age, sex, ECOG PS, tumor location, number of metastatic sites, KRAS , Consensus Molecular Subtypes (CMS), and treatment. Results: Overall, SNHG11 was strongly associated with OS and PFS after adjusting for multiple tests (Benjamini-Hochberg False Discovery Rate < 0.05). High SNHG11 expression (H group, n = 144) was associated with improved mPFS (H: 14.3 vs M: 11.2 vs L: 8.3 mo; p = 0.038) and mOS (H: 39.6 vs M: 31.1 vs L: 20.5 mo; p = 0.033) in the combined treatment analysis. Among cet-treated pts, SNHG11-H showed a numerically longer mPFS (H: 14.2 vs M: 11.1 vs L: 7.6 mo; p = 0.19) and significantly longer mOS (H: 41.1 vs M: 32.4 vs L: 14.3 mo; p = 0.012). In contrast, no statistically significant OS or PFS differences were observed in bev-treated pts. SNHG11-H tumors had a significant OS benefit from cet compared to bev (mOS 41.1 vs 36.5 mo, respectively; p = 0.016), with a nominally significant treatment interaction observed for OS (p = 0.030). No significant differences were observed in the L or M expression groups. Additional analyses showed that SNHG11 expression was high in the CMS2 (canonical) subtype and substantially lower in CMS1 (immune). Conclusions: LncRNA SNHG11 plays a significant role in CRC progression and metastasis via tumorigenic pathways, including c-Myc and HIF-1α. Moreover, elevated circulating SNHG11 levels show promise as a non-invasive biomarker for early CRC detection. In CALGB/SWOG 80405, high SNHG11 expression correlated with improved PFS and OS, particularly in cet-treated pts, supporting its role as a prognostic and predictive biomarker. Its strong association with CMS2 aligns with its reported involvement in c-Myc-driven pathways. Further validation is needed to confirm the clinical utility of this biomarker and elucidate underlying mechanisms.

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.

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.