Activation of innate-adaptive immune machinery by poly(I:C) exposes a therapeutic vulnerability to prevent relapse in stroma-rich colon cancer.

OBJECTIVE: Stroma-rich tumours represent a poor prognostic subtype in stage II/III colon cancer (CC), with high relapse rates and limited response to standard adjuvant chemotherapy. DESIGN: To address the lack of efficacious therapeutic options for patients with stroma-rich CC, we stratified our human tumour cohorts according to stromal content, enabling identification of the biology underpinning relapse and potential therapeutic vulnerabilities specifically within stroma-rich tumours that could be exploited clinically. Following human tumour-based discovery and independent clinical validation, we use a series of in vitro and stroma-rich in vivo models to test and validate the therapeutic potential of elevating the biology associated with reduced relapse in human tumours. RESULTS: By performing our analyses specifically within the stroma-rich/high-fibroblast (HiFi) subtype of CC, we identify and validate the clinical value of a HiFi-specific prognostic signature (HPS), which stratifies tumours based on STAT1-related signalling (High-HPS v Low-HPS=HR 0.093, CI 0.019 to 0.466). Using in silico, in vitro and in vivo models, we demonstrate that the HPS is associated with antigen processing and presentation within discrete immune lineages in stroma-rich CC, downstream of double-stranded RNA and viral response signalling. Treatment with the TLR3 agonist poly(I:C) elevated the HPS signalling and antigen processing phenotype across in vitro and in vivo models. In an in vivo model of stroma-rich CC, poly(I:C) treatment significantly increased systemic cytotoxic T cell activity (p<0.05) and reduced liver metastases (p<0.0002). CONCLUSION: This study reveals new biological insight that offers a novel therapeutic option to reduce relapse rates in patients with the worst prognosis CC.

Image-based consensus molecular subtype (imCMS) classification of colorectal cancer using deep learning.

OBJECTIVE: Complex phenotypes captured on histological slides represent the biological processes at play in individual cancers, but the link to underlying molecular classification has not been clarified or systematised. In colorectal cancer (CRC), histological grading is a poor predictor of disease progression, and consensus molecular subtypes (CMSs) cannot be distinguished without gene expression profiling. We hypothesise that image analysis is a cost-effective tool to associate complex features of tissue organisation with molecular and outcome data and to resolve unclassifiable or heterogeneous cases. In this study, we present an image-based approach to predict CRC CMS from standard H&E sections using deep learning. DESIGN: Training and evaluation of a neural network were performed using a total of n=1206 tissue sections with comprehensive multi-omic data from three independent datasets (training on FOCUS trial, n=278 patients; test on rectal cancer biopsies, GRAMPIAN cohort, n=144 patients; and The Cancer Genome Atlas (TCGA), n=430 patients). Ground truth CMS calls were ascertained by matching random forest and single sample predictions from CMS classifier. RESULTS: Image-based CMS (imCMS) accurately classified slides in unseen datasets from TCGA (n=431 slides, AUC)=0.84) and rectal cancer biopsies (n=265 slides, AUC=0.85). imCMS spatially resolved intratumoural heterogeneity and provided secondary calls correlating with bioinformatic prediction from molecular data. imCMS classified samples previously unclassifiable by RNA expression profiling, reproduced the expected correlations with genomic and epigenetic alterations and showed similar prognostic associations as transcriptomic CMS. CONCLUSION: This study shows that a prediction of RNA expression classifiers can be made from H&E images, opening the door to simple, cheap and reliable biological stratification within routine workflows.

Identification and validation of a machine learning model of complete response to radiation in rectal cancer reveals immune infiltrate and TGFß as key predictors.

BACKGROUND: It is uncertain which biological features underpin the response of rectal cancer (RC) to radiotherapy. No biomarker is currently in clinical use to select patients for treatment modifications. METHODS: We identified two cohorts of patients (total N = 249) with RC treated with neoadjuvant radiotherapy (45Gy/25) plus fluoropyrimidine. This discovery set included 57 cases with pathological complete response (pCR) to chemoradiotherapy (23%). Pre-treatment cancer biopsies were assessed using transcriptome-wide mRNA expression and targeted DNA sequencing for copy number and driver mutations. Biological candidate and machine learning (ML) approaches were used to identify predictors of pCR to radiotherapy independent of tumour stage. Findings were assessed in 107 cases from an independent validation set (GSE87211). FINDINGS: Three gene expression sets showed significant independent associations with pCR: Fibroblast-TGFß Response Signature (F-TBRS) with radioresistance; and cytotoxic lymphocyte (CL) expression signature and consensus molecular subtype CMS1 with radiosensitivity. These associations were replicated in the validation cohort. In parallel, a gradient boosting machine model comprising the expression of 33 genes generated in the discovery cohort showed high performance in GSE87211 with 90% sensitivity, 86% specificity. Biological and ML signatures indicated similar mechanisms underlying radiation response, and showed better AUC and p-values than published transcriptomic signatures of radiation response in RC. INTERPRETATION: RCs responding completely to chemoradiotherapy (CRT) have biological characteristics of immune response and absence of immune inhibitory TGFß signalling. These tumours may be identified with a potential biomarker based on a 33 gene expression signature. This could help select patients likely to respond to treatment with a primary radiotherapy approach as for anal cancer. Conversely, those with predicted radioresistance may be candidates for clinical trials evaluating addition of immune-oncology agents and stromal TGFß signalling inhibition. FUNDING: The Stratification in Colorectal Cancer Consortium (S:CORT) was funded by the Medical Research Council and Cancer Research UK (MR/M016587/1).

Tumor-associated mesenchymal stromal cells modulate macrophage phagocytosis in stromal-rich colorectal cancer via PD-1 signaling.

CMS4 colorectal cancer (CRC), based on the consensus molecular subtype (CMS), stratifies patients with the poorest disease-free survival rates. It is characterized by a strong mesenchymal stromal cell (MSC) signature, wound healing-like inflammation and therapy resistance. We utilized 2D and 3D in vitro, in vivo, and ex vivo models to assess the impact of inflammation and stromal cells on immunosuppression in CMS4 CRC. RNA sequencing data from untreated stage II/III CRC patients showed enriched TNF-α signatures in CMS1 and CMS4 tumors. Secretome from TNF-α treated cancer cells induced an immunomodulatory and chemotactic phenotype in MSC and cancer-associated fibroblasts (CAFs). Macrophages in CRC tumours migrate and preferentially localise in stromal compartment. Inflammatory CRC secretome enhances expression of PD-L1 and CD47 on both human and murine stromal cells. We demonstrate that TNF-α-induced inflammation in CRC suppresses macrophage phagocytosis via stromal cells. We show that stromal cell-mediated suppression of macrophage phagocytosis is mediated in part through PD-1 signaling. These data suggest that re-stratification of CRC by CMS may reveal patient subsets with microsatellite stable tumors, particularly CMS4-like tumors, that may respond to immunotherapies.

Stratification to Neoadjuvant Radiotherapy in Rectal Cancer by Regimen and Transcriptional Signatures.

Response to neoadjuvant radiotherapy (RT) in rectal cancer has been associated with immune and stromal features that are captured by transcriptional signatures. However, how such associations perform across different chemoradiotherapy regimens and within individual consensus molecular subtypes (CMS) and how they affect survival remain unclear. In this study, gene expression and clinical data of pretreatment biopsies from nine cohorts of primary rectal tumors were combined (N = 826). Exploratory analyses were done with transcriptomic signatures for the endpoint of pathologic complete response (pCR), considering treatment regimen or CMS subtype. Relevant findings were tested for overall survival and recurrence-free survival. Immune and stromal signatures were strongly associated with pCR and lack of pCR, respectively, in RT and capecitabine (Cap)/5-fluorouracil (5FU)-treated patients (N = 387), in which the radiosensitivity signature (RSS) showed the strongest association. Upon addition of oxaliplatin (Ox; N = 123), stromal signatures switched direction and showed higher chances to achieve pCR than without Ox (p for interaction 0.02). Among Cap/5FU patients, most signatures performed similarly across CMS subtypes, except cytotoxic lymphocytes that were associated with pCR in CMS1 and CMS4 cases compared with other CMS subtypes (p for interaction 0.04). The only variables associated with survival were pCR and RSS. Although the frequency of pCR across different chemoradiation regimens is relatively similar, our data suggest that response rates may differ depending on the biological landscape of rectal cancer. Response to neoadjuvant RT in stroma-rich tumors may potentially be improved by the addition of Ox. RSS in preoperative biopsies provides predictive information for response specifically to neoadjuvant RT with 5FU. SIGNIFICANCE: Rectal cancers with stromal features may respond better to RT and 5FU/Cap with the addition of Ox. Within patients not treated with Ox, high levels of cytotoxic lymphocytes associate with response only in immune and stromal tumors. Our analyses provide biological insights about the outcome by different radiotherapy regimens in rectal cancer.

Pathway level subtyping identifies a slow-cycling biological phenotype associated with poor clinical outcomes in colorectal cancer.

Molecular stratification using gene-level transcriptional data has identified subtypes with distinctive genotypic and phenotypic traits, as exemplified by the consensus molecular subtypes (CMS) in colorectal cancer (CRC). Here, rather than gene-level data, we make use of gene ontology and biological activation state information for initial molecular class discovery. In doing so, we defined three pathway-derived subtypes (PDS) in CRC: PDS1 tumors, which are canonical/LGR5+ stem-rich, highly proliferative and display good prognosis; PDS2 tumors, which are regenerative/ANXA1+ stem-rich, with elevated stromal and immune tumor microenvironmental lineages; and PDS3 tumors, which represent a previously overlooked slow-cycling subset of tumors within CMS2 with reduced stem populations and increased differentiated lineages, particularly enterocytes and enteroendocrine cells, yet display the worst prognosis in locally advanced disease. These PDS3 phenotypic traits are evident across numerous bulk and single-cell datasets, and demark a series of subtle biological states that are currently under-represented in pre-clinical models and are not identified using existing subtyping classifiers.

Targeting stromal cell sialylation reverses T cell-mediated immunosuppression in the tumor microenvironment.

Immunosuppressive tumor microenvironments (TMEs) reduce the effectiveness of immune responses in cancer. Mesenchymal stromal cells (MSCs), precursors to cancer-associated fibroblasts (CAFs), promote tumor progression by enhancing immune cell suppression in colorectal cancer (CRC). Hyper-sialylation of glycans promotes immune evasion in cancer through binding of sialic acids to their receptors, Siglecs, expressed on immune cells, which results in inhibition of effector functions. The role of sialylation in shaping MSC/CAF immunosuppression in the TME is not well characterized. In this study, we show that tumor-conditioned stromal cells have increased sialyltransferase expression, α2,3/6-linked sialic acid, and Siglec ligands. Tumor-conditioned stromal cells and CAFs induce exhausted immunomodulatory CD8+ PD1+ and CD8+ Siglec-7+/Siglec-9+ T cell phenotypes. In vivo, targeting stromal cell sialylation reverses stromal cell-mediated immunosuppression, as shown by infiltration of CD25 and granzyme B-expressing CD8+ T cells in the tumor and draining lymph node. Targeting stromal cell sialylation may overcome immunosuppression in the CRC TME.

Molecular Subtyping Resource: a user-friendly tool for rapid biological discovery from transcriptional data.

Generation of transcriptional data has dramatically increased in the past decade, driving the development of analytical algorithms that enable interrogation of the biology underpinning the profiled samples. However, these resources require users to have expertise in data wrangling and analytics, reducing opportunities for biological discovery by 'wet-lab' users with a limited programming skillset. Although commercial solutions exist, costs for software access can be prohibitive for academic research groups. To address these challenges, we have developed an open source and user-friendly data analysis platform for on-the-fly bioinformatic interrogation of transcriptional data derived from human or mouse tissue, called Molecular Subtyping Resource (MouSR). This internet-accessible analytical tool, https://mousr.qub.ac.uk/, enables users to easily interrogate their data using an intuitive 'point-and-click' interface, which includes a suite of molecular characterisation options including quality control, differential gene expression, gene set enrichment and microenvironmental cell population analyses from RNA sequencing. The MouSR online tool provides a unique freely available option for users to perform rapid transcriptomic analyses and comprehensive interrogation of the signalling underpinning transcriptional datasets, which alleviates a major bottleneck for biological discovery. This article has an associated First Person interview with the first author of the paper.

Biological Misinterpretation of Transcriptional Signatures in Tumor Samples Can Unknowingly Undermine Mechanistic Understanding and Faithful Alignment with Preclinical Data.

PURPOSE: Precise mechanism-based gene expression signatures (GES) have been developed in appropriate in vitro and in vivo model systems, to identify important cancer-related signaling processes. However, some GESs originally developed to represent specific disease processes, primarily with an epithelial cell focus, are being applied to heterogeneous tumor samples where the expression of the genes in the signature may no longer be epithelial-specific. Therefore, unknowingly, even small changes in tumor stroma percentage can directly influence GESs, undermining the intended mechanistic signaling. EXPERIMENTAL DESIGN: Using colorectal cancer as an exemplar, we deployed numerous orthogonal profiling methodologies, including laser capture microdissection, flow cytometry, bulk and multiregional biopsy clinical samples, single-cell RNA sequencing and finally spatial transcriptomics, to perform a comprehensive assessment of the potential for the most widely used GESs to be influenced, or confounded, by stromal content in tumor tissue. To complement this work, we generated a freely-available resource, ConfoundR; https://confoundr.qub.ac.uk/, that enables users to test the extent of stromal influence on an unlimited number of the genes/signatures simultaneously across colorectal, breast, pancreatic, ovarian and prostate cancer datasets. RESULTS: Findings presented here demonstrate the clear potential for misinterpretation of the meaning of GESs, due to widespread stromal influences, which in-turn can undermine faithful alignment between clinical samples and preclinical data/models, particularly cell lines and organoids, or tumor models not fully recapitulating the stromal and immune microenvironment. CONCLUSIONS: Efforts to faithfully align preclinical models of disease using phenotypically-designed GESs must ensure that the signatures themselves remain representative of the same biology when applied to clinical samples.

Epithelial TGFß engages growth-factor signalling to circumvent apoptosis and drive intestinal tumourigenesis with aggressive features.

The pro-tumourigenic role of epithelial TGFß signalling in colorectal cancer (CRC) is controversial. Here, we identify a cohort of born to be bad early-stage (T1) colorectal tumours, with aggressive features and a propensity to disseminate early, that are characterised by high epithelial cell-intrinsic TGFß signalling. In the presence of concurrent Apc and Kras mutations, activation of epithelial TGFß signalling rampantly accelerates tumourigenesis and share transcriptional signatures with those of the born to be bad T1 human tumours and predicts recurrence in stage II CRC. Mechanistically, epithelial TGFß signalling induces a growth-promoting EGFR-signalling module that synergises with mutant APC and KRAS to drive MAPK signalling that re-sensitise tumour cells to MEK and/or EGFR inhibitors. Together, we identify epithelial TGFß signalling both as a determinant of early dissemination and a potential therapeutic vulnerability of CRC's with born to be bad traits.

In-depth Clinical and Biological Exploration of DNA Damage Immune Response as a Biomarker for Oxaliplatin Use in Colorectal Cancer.

PURPOSE: The DNA damage immune response (DDIR) assay was developed in breast cancer based on biology associated with deficiencies in homologous recombination and Fanconi anemia pathways. A positive DDIR call identifies patients likely to respond to platinum-based chemotherapies in breast and esophageal cancers. In colorectal cancer, there is currently no biomarker to predict response to oxaliplatin. We tested the ability of the DDIR assay to predict response to oxaliplatin-based chemotherapy in colorectal cancer and characterized the biology in DDIR-positive colorectal cancer. EXPERIMENTAL DESIGN: Samples and clinical data were assessed according to DDIR status from patients who received either 5-fluorouracil (5-FU) or 5FUFA (bolus and infusion 5-FU with folinic acid) plus oxaliplatin (FOLFOX) within the FOCUS trial (n = 361, stage IV), or neoadjuvant FOLFOX in the FOxTROT trial (n = 97, stage II/III). Whole transcriptome, mutation, and IHC data of these samples were used to interrogate the biology of DDIR in colorectal cancer. RESULTS: Contrary to our hypothesis, DDIR-negative patients displayed a trend toward improved outcome for oxaliplatin-based chemotherapy compared with DDIR-positive patients. DDIR positivity was associated with microsatellite instability (MSI) and colorectal molecular subtype 1. Refinement of the DDIR signature, based on overlapping IFN-related chemokine signaling associated with DDIR positivity across colorectal cancer and breast cancer cohorts, further confirmed that the DDIR assay did not have predictive value for oxaliplatin-based chemotherapy in colorectal cancer. CONCLUSIONS: DDIR positivity does not predict improved response following oxaliplatin treatment in colorectal cancer. However, data presented here suggest the potential of the DDIR assay in identifying immune-rich tumors that may benefit from immune checkpoint blockade, beyond current use of MSI status.

Fibroblast-derived Gremlin1 localises to epithelial cells at the base of the intestinal crypt.

Gremlin1 (GREM1) is a secreted glycoprotein member of the differential screening-selected gene in aberrant neuroblastoma (DAN) family of bone morphogenetic protein (BMP) antagonists, which binds to BMPs preventing their receptor engagement. Previous studies have identified that stage II colorectal cancer (CRC) patients with high levels of GREM1 gene expression in their tumour tissue have a poorer prognosis. Using a series of in silico and in situ methodologies, we demonstrate that GREM1 gene expression is significantly higher (p < 0.0001) in CRC consensus molecular subtype 4 (CMS4), compared to the other CMS subtypes and correlates (p < 0.0001) with levels of cancer-associated fibroblasts (CAFs) within the CRC tumour microenvironment (TME). Our optimised immunohistochemistry protocol identified endogenous GREM1 protein expression in both the muscularis mucosa and adjacent colonic crypt bases in mouse intestine, in contrast to RNA expression which was shown to localise specifically to the muscularis mucosa, as determined by in situ hybridisation. Importantly, we demonstrate that cells with high levels of GREM1 expression display low levels of phospho-Smad1/5, consistent with reduced BMP signalling. Taken together, these data highlight a novel paracrine signalling circuit, which involves uptake of mature GREM1 protein by colonic crypt cells following secretion from neighbouring fibroblasts in the TME.

EphA2 Expression Is a Key Driver of Migration and Invasion and a Poor Prognostic Marker in Colorectal Cancer.

PURPOSE: EphA2, a member of the Eph receptor tyrosine kinases family, is an important regulator of tumor initiation, neovascularization, and metastasis in a wide range of epithelial and mesenchymal cancers; however, its role in colorectal cancer recurrence and progression is unclear. EXPERIMENTAL DESIGN: EphA2 expression was determined by immunohistochemistry in stage II/III colorectal tumors (N = 338), and findings correlated with clinical outcome. The correlation between EphA2 expression and stem cell markers CD44 and Lgr5 was examined. The role of EphA2 in migration/invasion was assessed using a panel of KRAS wild-type (WT) and mutant (MT) parental and invasive colorectal cancer cell line models. RESULTS: Colorectal tumors displayed significantly higher expression levels of EphA2 compared with matched normal tissue, which positively correlated with high CD44 and Lgr5 expression levels. Moreover, high EphA2 mRNA and protein expression were found to be associated with poor overall survival in stage II/III colorectal cancer tissues, in both univariate and multivariate analyses. Preclinically, we found that EphA2 was highly expressed in KRASMT colorectal cancer cells and that EphA2 levels are regulated by the KRAS-driven MAPK and RalGDS-RalA pathways. Moreover, EphA2 levels were elevated in several invasive daughter cell lines, and downregulation of EphA2 using RNAi or recombinant EFNA1 suppressed migration and invasion of KRASMT colorectal cancer cells. CONCLUSIONS: These data show that EpHA2 is a poor prognostic marker in stage II/III colorectal cancer, which may be due to its ability to promote cell migration and invasion, providing support for the further investigation of EphA2 as a novel prognostic biomarker and therapeutic target.

Overcoming Resistance to Targeted Therapies in Cancer.

The recent discovery of oncogenic drivers and subsequent development of novel targeted strategies has significantly added to the therapeutic armamentarium of anti-cancer therapies. Targeting BCR-ABL in chronic myeloid leukemia (CML) or HER2 in breast cancer has led to practice-changing clinical benefits, while promising therapeutic responses have been achieved by precision medicine approaches in EGFR mutant lung cancer, colorectal cancer and BRAF mutant melanoma. However, although initial therapeutic responses to targeted therapies can be substantial, many patients will develop disease progression within 6-12 months. An increasing application of powerful omics-based approaches and improving preclinical models have enabled the rapid identification of secondary resistance mechanisms. Herein, we discuss how this knowledge has translated into rational, novel treatment strategies for relapsed patients in genomically selected cancer populations.

AXL is a key regulator of inherent and chemotherapy-induced invasion and predicts a poor clinical outcome in early-stage colon cancer.

PURPOSE: Despite the use of 5-fluorouracil (5-FU)-based adjuvant treatments, a large proportion of patients with high-risk stage II/III colorectal cancer will relapse. Thus, novel therapeutic strategies are needed for early-stage colorectal cancer. Residual micrometastatic disease from the primary tumor is a major cause of patient relapse. EXPERIMENTAL DESIGN: To model colorectal cancer tumor cell invasion/metastasis, we have generated invasive (KRASMT/KRASWT/+chr3/p53-null) colorectal cancer cell subpopulations. Receptor tyrosine kinase (RTK) screens were used to identify novel proteins that underpin the migratory/invasive phenotype. Migration/invasion was assessed using the XCELLigence system. Tumors from patients with early-stage colorectal cancer (N = 336) were examined for AXL expression. RESULTS: Invasive colorectal cancer cell subpopulations showed a transition from an epithelial-to-mesenchymal like phenotype with significant increases in migration, invasion, colony-forming ability, and an attenuation of EGF receptor (EGFR)/HER2 autocrine signaling. RTK arrays showed significant increases in AXL levels in all invasive sublines. Importantly, 5-FU treatment resulted in significantly increased migration and invasion, and targeting AXL using pharmacologic inhibition or RNA interference (RNAi) approaches suppressed basal and 5-FU-induced migration and invasion. Significantly, high AXL mRNA and protein expression were found to be associated with poor overall survival in early-stage colorectal cancer tissues. CONCLUSIONS: We have identified AXL as a poor prognostic marker and important mediator of cell migration/invasiveness in colorectal cancer. These findings provide support for the further investigation of AXL as a novel prognostic biomarker and therapeutic target in colorectal cancer, in particular in the adjuvant disease in which EGFR/VEGF-targeted therapies have failed.

ADAM17-dependent c-MET-STAT3 signaling mediates resistance to MEK inhibitors in KRAS mutant colorectal cancer.

There are currently no approved targeted therapies for advanced KRAS mutant (KRASMT) colorectal cancer (CRC). Using a unique systems biology approach, we identified JAK1/2-dependent activation of STAT3 as the key mediator of resistance to MEK inhibitors in KRASMT CRC in vitro and in vivo. Further analyses identified acute increases in c-MET activity following treatment with MEK inhibitors in KRASMT CRC models, which was demonstrated to promote JAK1/2-STAT3-mediated resistance. Furthermore, activation of c-MET following MEK inhibition was found to be due to inhibition of the ERK-dependent metalloprotease ADAM17, which normally inhibits c-MET signaling by promoting shedding of its endogenous antagonist, soluble "decoy" MET. Most importantly, pharmacological blockade of this resistance pathway with either c-MET or JAK1/2 inhibitors synergistically increased MEK-inhibitor-induced apoptosis and growth inhibition in vitro and in vivo in KRASMT models, providing clear rationales for the clinical assessment of these combinations in KRASMT CRC patients.