Search

Purpose:Molecular tumor heterogeneity may have important implications for the efficacy of targeted therapies in metastatic cancers. Inter-metastatic heterogeneity of sensitivity to anticancer agents has not been well explored in colorectal cancer.Experimental Design:We established a platform for ex vivo pharmacogenomic profiling of patient-derived organoids (PDO) from resected colorectal cancer liver metastases. Drug sensitivity testing (n = 40 clinically relevant agents) and gene expression profiling were performed on 39 metastases from 22 patients.Results:Three drug–response clusters were identified among the colorectal cancer metastases, based primarily on sensitivities to EGFR and/or MDM2 inhibition, and corresponding with RAS mutations and TP53 activity. Potentially effective therapies, including off-label use of drugs approved for other cancer types, could be nominated for eighteen patients (82%). Antimetabolites and targeted agents lacking a decisive genomic marker had stronger differential activity than most approved chemotherapies. We found limited intra-patient drug sensitivity heterogeneity between PDOs from multiple (2–5) liver metastases from each of ten patients. This was recapitulated at the gene expression level, with a highly proportional degree of transcriptomic and pharmacological variation. One PDO with a multi-drug resistance profile, including resistance to EGFR inhibition in a RAS-mutant background, showed sensitivity to MEK plus mTOR/AKT inhibition, corresponding with low-level PTEN expression.Conclusions:Intra-patient inter-metastatic pharmacological heterogeneity was not pronounced and ex vivo drug screening may identify novel treatment options for metastatic colorectal cancer. Variation in drug sensitivities was reflected at the transcriptomic level, suggesting potential to develop gene expression–based predictive signatures to guide experimental therapies.

This file contains Supplementary Text, Supplementary Figures S1-S13 and Supplementary Tables S1, S3-S5, S9-S10, S14-S15, S18-S19. Supplementary Fig. S1. Validation of prognostic associations and biological properties of the CMS groups in primary CRCs Supplementary Fig. S2. The original CMS classifier fails to accurately identify the CMS groups in CRC cell lines Supplementary Fig. S3. Principal components analysis of cell lines based on gene expression Supplementary Fig. S4. Assessment of misclassification risk of CRC cell lines Supplementary Fig. S5. Additional template gene filter for classifier genes expressed in the tumor stroma Supplementary Fig. S6. Biological characteristics of the CMS groups are recapitulated in PDX models of CRC Supplementary Fig. S7. CRC cell lines analyzed by drug screening are representative with respect to gene expression and CMS groups Supplementary Fig. S8. Principal components analysis of in vitro drug screen data Supplementary Fig. S9. Validation of strong relative activity of EGFR inhibition in CMS2 Supplementary Fig. S10. Drug response analyses between CMS1/4 and CMS2/3 cell lines Supplementary Fig. S11. In vitro validation of strong relative sensitivity to HSP90 inhibition in CMS1 and CMS4 cell lines in published data Supplementary Fig. S12. Up-regulation of heat shock response in cell lines treated with luminespib Supplementary Fig. S13. Poor response to fluoropyrimidines in CMS4 cell lines Supplementary Table S1. Clinicopathological and molecular characteristics of consecutive patient series (n = 409) Supplementary Table S3. Enrichment of molecular characteristics among CMS groups in primary CRC Supplementary Table S4. Enrichment of clinicopathological characteristics among CMS groups in primary CRC Supplementary Table S5. Multivariable survival analyses Supplementary Table S9. Classification concordance between the original and new cancer cell-adapted CMS classifier in CRC cell lines Supplementary Table S10. Enrichment of molecular characteristics among CMS groups of CRC cell lines Supplementary Table S14. Overview of the molecular target classes for the highthroughput drug screening library Supplementary Table S15. Drugs with strong effect in MSI+ versus MSS cell lines Supplementary Table S18. Drugs with strong relative activity in CMS2 Supplementary Table S19. Differential gene expression in response to luminespib treatment

Purpose: Improved prognostic stratification of patients with stage II and III colorectal cancer is warranted for postoperative clinical decision making. This study was conducted to develop a clinically feasible and robust prognostic classifier for these patients independent of adjuvant treatment.Experimental Design: Global gene expression profiles from altogether 387 stage II and III colorectal cancer tissue samples from three independent patient series were included in the study. ColoGuidePro, a seven-gene prognostic classifier, was developed from a selected Norwegian learning series (n = 95; no adjuvant treatment) using lasso-penalized multivariate survival modeling with cross-validation.Results: The expression signature significantly stratified patients in a consecutive Norwegian test series, in which patients were treated according to current standards [HR, 2.9 (1.1–7.5); P = 0.03; n = 77] and an external validation series [HR, 3.7 (2.0–6.8); P < 0.001; n = 215] according to survival. ColoGuidePro was also an independent predictor of prognosis in multivariate models including tumor stage in both series (HR, ≥3.1; P ≤ 0.03). In the validation series, which consisted of patients from other populations (United States and Australia), 5-year relapse-free survival was significantly predicted for stage III patients only (P < 0.001; n = 107). Here, prognostic stratification was independent of adjuvant treatment (P = 0.001).Conclusions: We present ColoGuidePro, a prognostic classifier developed for patients with stage II and III colorectal cancer. The test is suitable for transfer to clinical use and has best prognostic prediction potential for stage III patients. Clin Cancer Res; 18(21); 6001–10. ©2012 AACR.

PDF file - 292K, The Supplementary Material contains Supplementary Methods, Supplementary Tables S1 and S2, and Supplementary Figures S1-S4 Table S1. Genes in the prognostic expression signature Table S2. Cox proportional hazards analyses for patients in the learning series Figure S1. Cross-validated partial likelihood and number of active predictors in the learning series as a function of the penalty parameter Figure S2. Expression in the learning series of the genes in the 7-gene prognostic signature. Figure S3. Survival curves for stage III patients in the validation series. Figure S4. Survival curves for stage III patients in the validation series aged 75 years or older

This file contains Supplementary Tables S2, S6-S8, S11-S13, and S16-S17. Supplementary Table S2. Overview and CMS groups of CRC cell line Supplementary Table S6. Gene set analysis among CMS groups in primary CRCs Supplementary Table S7. Potential template genes and filtering steps Supplementary Table S8. Template genes for the adapted CMS classifier Supplementary Table S11. Gene set analyses among CMS groups in CRC cell lines Supplementary Table S12. Assessment of misclassification risk in CRC cell lines Supplementary Table S13. CMS classification of PDX models by adapted CMS classifier Supplementary Table S16. Differential drug sensitivity between cell lines of individual CMS groups Supplementary Table S17. Drug response data for 27 selected drugs in 29 individual cell lines

Purpose: Response to standard oncologic treatment is limited in colorectal cancer. The gene expression–based consensus molecular subtypes (CMS) provide a new paradigm for stratified treatment and drug repurposing; however, drug discovery is currently limited by the lack of translation of CMS to preclinical models.Experimental Design: We analyzed CMS in primary colorectal cancers, cell lines, and patient-derived xenografts (PDX). For classification of preclinical models, we developed an optimized classifier enriched for cancer cell–intrinsic gene expression signals, and performed high-throughput in vitro drug screening (n = 459 drugs) to analyze subtype-specific drug sensitivities.Results: The distinct molecular and clinicopathologic characteristics of each CMS group were validated in a single-hospital series of 409 primary colorectal cancers. The new, cancer cell–adapted classifier was found to perform well in primary tumors, and applied to a panel of 148 cell lines and 32 PDXs, these colorectal cancer models were shown to recapitulate the biology of the CMS groups. Drug screening of 33 cell lines demonstrated subtype-dependent response profiles, confirming strong response to EGFR and HER2 inhibitors in the CMS2 epithelial/canonical group, and revealing strong sensitivity to HSP90 inhibitors in cells with the CMS1 microsatellite instability/immune and CMS4 mesenchymal phenotypes. This association was validated in vitro in additional CMS-predicted cell lines. Combination treatment with 5-fluorouracil and luminespib showed potential to alleviate chemoresistance in a CMS4 PDX model, an effect not seen in a chemosensitive CMS2 PDX model.Conclusions: We provide translation of CMS classification to preclinical models and uncover a potential for targeted treatment repurposing in the chemoresistant CMS4 group. Clin Cancer Res; 24(4); 794–806. ©2017 AACR.

Supplementary methods, materials and figures

Background Colorectal cancer is the 2nd leading cause of cancer-related deaths with few patients benefiting from biomarker-guided therapy. Mutation expression is essential for accurate interpretation of mutations as biomarkers, but surprisingly, little has been done to analyze somatic cancer mutations on the expression level. We report a large-scale analysis of allele-specific mutation expression. Methods Whole-exome and total RNA sequencing was performed on 137 samples from 121 microsatellite stable colorectal cancers, including multiregional samples of primary and metastatic tumors from 4 patients. Data were integrated with allele-specific resolution. Results were validated in an independent set of 241 colon cancers. Therapeutic associations were explored by pharmacogenomic profiling of 15 cell lines or patient-derived organoids. Results The median proportion of expressed mutations per tumor was 34%. Cancer-critical mutations had the highest expression frequency (gene-wise mean of 58%), independent of frequent allelic imbalance. Systematic deviation from the general pattern of expression levels according to allelic frequencies was detected, including preferential expression of mutated alleles dependent on the mutation type and target gene. Translational relevance was suggested by correlations of KRAS/NRAS or TP53 mutation expression levels with downstream oncogenic signatures (p < 0.03), overall survival among patients with stage II and III cancer (KRAS/NRAS: hazard ratio 6.1, p = 0.0070), and targeted drug sensitivity. The latter was demonstrated for EGFR and MDM2 inhibition in pre-clinical models. Conclusions Only a subset of mutations in microsatellite stable colorectal cancers were expressed, and the “expressed mutation dose” may provide an opportunity for more fine-tuned biomarker interpretations.

Gene expression-based subtypes of colorectal cancer have clinical relevance, but the representativeness of primary tumors and the consensus molecular subtypes (CMS) for metastatic cancers is not well known. We investigated the metastatic heterogeneity of CMS. The best approach to subtype translation was delineated by comparisons of transcriptomic profiles from 317 primary tumors and 295 liver metastases, including multi-metastatic samples from 45 patients and 14 primary-metastasis sets. Associations were validated in an external data set (n = 618). Projection of metastases onto principal components of primary tumors showed that metastases were depleted of CMS1-immune/CMS3-metabolic signals, enriched for CMS4-mesenchymal/stromal signals, and heavily influenced by the microenvironment. The tailored CMS classifier (available in an updated version of the R package CMScaller) therefore implemented an approach to regress out the liver tissue background. The majority of classified metastases were either CMS2 or CMS4. Nonetheless, subtype switching and inter-metastatic CMS heterogeneity were frequent and increased with sampling intensity. Poor-prognostic value of CMS1/3 metastases was consistent in the context of intra-patient tumor heterogeneity.

Purpose: Tumor heterogeneity compromises the clinical applicability of transcriptomic subtypes of colorectal cancers (CRCs). We performed multiregional tumor sampling to map the intra-tumor heterogeneity of primary CRCs and to evaluate the prognostic relevance of features less vulnerable to heterogeneity. Methods: Gene expression profiling was performed of 708 primary tumor samples from 515 patients treated by surgery for stage I-IV CRC at Oslo University Hospital (Human Transcriptome 2.0 Arrays). This included 2-4 multiregional samples from each of 97 patients (n=290 samples; mean 2.9 per tumor), and single bulk tumor samples from 418 patients. Intra-tumor heterogeneity was evaluated according to consensus molecular subtype (CMS) classification, including computational assessment in the single samples by CMS enrichment scores. Gene-wise heterogeneity scores were calculated as intra-class correlation coefficients in the multiregional sample set, and genes with low heterogeneity were used as input for subtype discovery by non-negative matrix factorization. Results: Heterogeneous CMS classification of multiregional samples was found in 40% of tumors. Heterogeneous tumors were enriched with CMS3 (odds ratio 5.0) and CMS4 (odds ratio 12.5), and the most common combinations were CMS2/4 and CMS1/3. Microsatellite instability and BRAFV600E mutations were enriched among tumors with a major CMS1 component only. KRAS/NRAS mutations were most frequent in CMS3 tumors without CMS1. Intra-tumor heterogeneity was primarily driven by mesenchymal-like and stromal signals, but was independently associated with a poor 5-year relapse-free survival among patients with stage I-III CRC (multivariable hazard ratio 1.5 [1.0-2.2]). Heterogeneity explained a larger proportion of variation in survival (14%) than cancer-associated fibroblasts (6%). Genes with low intra-tumor heterogeneity were enriched in cancer cell intrinsic signaling pathways, including cell cycle progression and MYC targets. De novo subtyping based on these genes recapitulated the intrinsic iCMS groups identified from single-cell sequencing (Joanito et al., Nat Genet 2022;54:963-75) with high accuracy (Cohen’s κ=0.80). Further sub-stratification identified four subtypes with similar biological activity to the original CMS, but lower classification heterogeneity among multiregional samples (25%). The subtypes were called congruent CMS (cCMS) and provided stronger prognostic stratification than CMS in multivariable analysis with clinicopathological and molecular features (cCMS4: hazard ratio 5.9 [2.4-14.7]; cCMS2: 3.4 [1.5-7.6]; cCMS3: 2.0 [0.9-4.3] with cCMS1 as reference). Conclusion: Multiregional sampling of primary CRCs enabled identification of cancer-intrinsic transcriptomic features robust to intra-tumor heterogeneity and with prognostic relevance. Citation Format: Jonas Langerud, Ina A. Eilertsen, Hossein Moosavi, Ingeborg F. Backe, Merete Hektoen, Marine Jeanmougin, Ole Sjo, Arild Nesbakken, Ragnhild A. Lothe, Anita Sveen. Multiregional transcriptomics of colorectal cancers define prognostic features less vulnerable to tumor heterogeneity. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4540.

Regulatory T cells (Tregs) are a heterogeneous cell population that can either suppress or stimulate immune responses. Tumor-infiltrating Tregs are associated with an adverse outcome from most cancer types, but have generally been found to be associated with a good prognosis in colorectal cancer (CRC). We investigated the prognostic heterogeneity of Tregs in CRC by co-expression patterns and spatial analyses with diverse T cell markers, using multiplex fluorescence immunohistochemistry and digital image analysis in two consecutive series of primary CRCs (total n = 1720). Treg infiltration in tumors, scored as FOXP3+ or CD4+/CD25+/FOXP3+ (triple-positive) cells, was strongly correlated to the overall amount of CD3+ and CD8+ T cells, and consequently associated with a favorable 5-year relapse-free survival rate among patients with stage I–III CRC who underwent complete tumor resection. However, high relative expression of the activation marker CD25 in triple-positive Tregs was independently associated with an adverse outcome in a multivariable model incorporating clinicopathological and known molecular prognostic markers (hazard ratio = 1.35, p = 0.028). Furthermore, spatial marker analysis based on Voronoi diagrams and permutation testing of cellular neighborhoods revealed a statistically significant proximity between Tregs and CD8+-cells in 18% of patients, and this was independently associated with a poor survival (multivariable hazard ratio = 1.36, p = 0.017). These results show prognostic heterogeneity of different Treg populations in primary CRC, and highlight the importance of multi-marker and spatial analyses for accurate immunophenotyping of tumors in relation to patient outcome.

Glioblastoma (GBM), the most aggressive form of brain cancer, is characterized by a high level of molecular heterogeneity, and infiltration by various immune and stromal cell populations. Important advances have been made in deciphering the microenvironment of GBMs, but its association with existing molecular subtypes and its potential prognostic role remain elusive. We have investigated the abundance of infiltrating immune and stromal cells in silico, from gene expression profiles. Two cohorts, including in‐house normal brain and glioma samples (n = 70) and a large sample set from TCGA (n = 393), were combined into a single exploratory dataset. A third independent cohort (n = 124) was used for validation. Tumors were clustered based on their microenvironment infiltration profiles, and associations with known GBM molecular subtypes and patient outcome were tested a posteriori in a multivariable setting. We identified a subset of GBM samples with significantly higher abundances of most immune and stromal cell populations. This subset showed increased expression of both immune suppressor and immune effector genes compared to other GBMs and was enriched for the mesenchymal molecular subtype. Survival analyses suggested that tumor microenvironment infiltration pattern was an independent prognostic factor for GBM patients. Among all, patients with the mesenchymal subtype with low immune and stromal infiltration had the poorest survival. By combining molecular subtyping with gene expression measures of tumor infiltration, the present work contributes with improving prognostic models in GBM., This work combined molecular subtyping of glioblastomas (GBMs) with gene expression measures of tumor infiltration and identified a subset of GBMs with high abundances of most immune and stromal cell populations. Survival analyses suggested that tumor microenvironment infiltration was an independent prognostic factor for GBM patients. Among all, patients with the mesenchymal subtype and low infiltration had the poorest survival.

The majority of prostate cancer patients are diagnosed with multiple primary malignant foci. The distinct foci are exceptionally heterogeneous with regard to DNA mutations, but whether this is recapitulated at the transcriptome level remains unknown. In this study, inter- and intrafocal heterogeneity has been assessed by whole-transcriptome sequencing of 87 tissue samples from 23 patients with localized prostate cancer treated with radical prostatectomy. From each patient, multiple samples were taken from one or more malignant foci, in addition to one sample from benign prostate tissue. Transcriptomic profiles of different malignant foci from the same patient showed a similar level of heterogeneity as tumors from different patients. This applies to expression of genes, fusion genes, and somatic mutations. Within-patient pair-wise analyses identified expression patterns linked to ETS status and extraprostatic extension. A set of 62 genes were found with low intrapatient heterogeneity and high interpatient heterogeneity, retaining stable expression profiles across foci within the same patient. Among these, 16 genes are associated with biochemical recurrence in a separately published study and are therefore nominated as biomarkers with prognostic value regardless of which malignant focus is sampled. In conclusion, an extensive heterogeneity in multifocal prostate cancer is confirmed at the gene expression level. Diagnostic biomarkers were identified for ETS positive samples and samples from extraprostatic extensions. Finally, prognostic biomarkers independent of multifocal heterogeneity were found.

Most patients with colorectal cancer (CRC) do not have an “actionable” target based on molecular tumor profiling. Ex vivo drug sensitivity testing of tumor organoids holds promise as a complementary approach in precision oncology. We have established a pre-clinical pharmacogenomics platform for CRC at Oslo University Hospital. A living biobank of 208 organoids of distinct liver lesions from 100 patients treated by hepatic resection for metastatic CRC was generated in the observational phase of the project. Sensitivity testing of custom drug libraries (n=40-47 drugs) and molecular profiling (targeted DNA sequencing, RNA sequencing, multiplex immunohistochemistry) indicated that the organoids recapitulate well-known pharmacogenomic associations, phenotypic tumor features, and clinical treatment responses (n=35 patients published (ref 1,2)). Evaluation of multiple lesions (n=2-6) per patient suggested diversity among patients in the level of inter-metastatic heterogeneity of drug sensitivities, supporting the need for multi-tumor evaluations in a clinical setting. A collection of CRC cell lines (n=107) was screened for sensitivity to more than 500 drugs and this suggested a molecular basis for response to selected targeted agents, including to PARP inhibition in a small subset of the cell lines (ref 3). Combination screens suggested possibility to overcome acquired resistance to dual BRAF and EGFR blockade in BRAFV600E cells (unpublished data). The living biobank serves as a reference for drug nominations in an ongoing phase II umbrella trial of pharmacogenomics-guided treatment in patients with metastatic CRC (EVIDENT; EU Clinical Trials Register no. 2020-003395-41). Several ad hoc patient cases have been screened as a pilot for the trial, and data of clinical benefit from this functional precision oncology strategy will be presented. References: 1. Bruun J. et al. Patient-derived organoids from multiple colorectal cancer liver metastases reveal moderate intra-patient pharmacotranscriptomic heterogeneity. Clin Cancer Res 2020;26:4107-19. 2. Kryeziu K. et al. Increased sensitivity to SMAC mimetic LCL161 identified by longitudinal ex vivo pharmacogenomics of recurrent, KRAS mutated rectal cancer liver metastases. J. Transl Med 2021;19:384. 3. Smeby J. et al. Molecular correlates of sensitivity to PARP inhibition beyond homologous recombination deficiency in pre-clinical models of colorectal cancer point to wild-type TP53 activity. EBioMed 2020;59:102923. Citation Format: Anita Sveen, Kushtrim Kryeziu, Solveig M.K. Klokkerud, Kaja C.G. Berg, Max Z. Totland, Christian H. Bergsland, Barbara Niederdorfer, Seyed H. Moosavi, Eva Hofsli, Morten Brændengen, Kristoffer Lassen, Arild Nesbakken, Sheraz Yaqub, Tormod K. Guren, Ragnhild A. Lothe. Tumor organoids of multi-metastatic colorectal cancer: From research tools to treatment decision tools [abstract]. In: Proceedings of the AACR Special Conference on Colorectal Cancer; 2022 Oct 1-4; Portland, OR. Philadelphia (PA): AACR; Cancer Res 2022;82(23 Suppl_1):Abstract nr IA009.

Background Squamous cell carcinoma of the anus (SCCA) is a rare malignancy with rising incidence, associated with human papilloma virus (HPV). Chemoradiotherapy (CRT) is the preferred treatment. The purpose was to investigate treatment failure, survival and prognostic factors after CRT. Material and methods In this prospective observational study from a large regional centre, 141 patients were included from 2013 to 2017, and 132 were eligible for analysis. The main inclusion criteria were SCCA, planned radiotherapy, and performance status (ECOG) ≤2. Patient characteristics, disease stage, treatment, and treatment response were prospectively registered. Disease-free survival (DFS), overall survival (OS), and locoregional treatment failure after CRT were analysed. Hazard ratios (HRs) were estimated with Cox`s proportional hazards model. Results Median follow-up was 54 (range 6–71) months. Eighteen patients (14%) had treatment failures after CRT; of these 10 (8%) had residual tumour, and 8 (6%) relapse as first failure. The first treatment failure was locoregional (11 patients), distant (5 patients), and both (2 patients). Salvage abdomino-perineal resection was performed in 10 patients, 2 had resections of metastases, and 3 both. DFS was 85% at 3 years and 78% at 5 years. OS was 93% at 3 years and 86% at 5 years. In analyses adjusted for age and gender, HPV negative tumours (HR 2.5, p = 0.024), N3 disease (HR 2.6, p = 0.024), and tumour size ≥4 cm (HR 2.4, p = 0.038) were negative prognostic factors for DFS. Conclusion State-of-the-art chemoradiotherapy for SCCA resulted in excellent outcomes, and improved survival compared with previous national data, with 80%.

About 80% of colorectal cancers (CRCs) have chromosomal instability, which is an integral part of aggressive malignancy development, but the importance of specific copy number aberrations (CNAs) in modulating gene expression, particularly within the framework of clinically relevant molecular subtypes, remains mostly elusive. We performed DNA copy number profiling of 257 stage I-IV primary CRCs and integrative gene expression analysis in 151 microsatellite stable (MSS) tumors, focusing on high-level amplifications and the effect of CNAs on the characteristics of the gene expression-based consensus molecular subtypes (CMS). The results were validated in 323 MSS tumors from TCGA. Novel recurrent high-level amplifications (≥15 additional copies) with a major impact on gene expression were found for TOX3 (16q) at 1.5% frequency, as well as for CCND2 (12p) and ANXA11 (10q) at 1% frequency, in addition to the well-known targets ERBB2 (17q) and MYC (8q). Focal amplifications with ≥15 or ≥5 additional copies of at least one of these regions were associated with a poor overall survival among patients with stage I-III MSS CRCs (multivariable hazard ratio ≥3.2, p ≤ 0.01). All high-level amplifications were focal and had a more consistent relationship with gene expression than lower amplitude and/or broad-range amplifications, suggesting specific targeting during carcinogenesis. Genome-wide, copy number driven gene expression was enriched for pathways characteristic of the CMS2-epithelial/canonical subtype, including DNA repair and cell cycle progression. Furthermore, 50% of upregulated genes in CMS2-epithelial/canonical MSS CRCs were driven by CNAs, an enrichment compared with the other CMS groups, and associated with the stronger correspondence between CNAs and gene expression in malignant epithelial cells than in the cells of the tumor microenvironment (fibroblasts, endothelial cells, leukocytes). In conclusion, we identify novel recurrent amplifications with impact on gene expression in CRC and provide the first evidence that CMS2 may have a stronger copy-number related genetic basis than subtypes more heavily influenced by gene expression signals from the tumor microenvironment. publishedVersion

Elevated miR‐31 expression is associated with poor outcome in colorectal cancer (CRC). Whether the prognostic information is independent of known molecular subgroups and gene expression‐based consensus molecular subtypes (CMS) is currently unknown. To investigate this, we analyzed nearly 2000 CRC biopsies and preclinical models. The expression of miR‐31‐5p and its host transcript, long noncoding RNA MIR31HG, was strongly correlated (Spearman's ρ > 0.80). MIR31HG outlier expression was observed in 158/1265 (12%) of pCRCs and was associated with depletion of CMS2‐canonical subgroup (odds ratio = 0.21 [0.11–0.35]) and shorter relapse‐free survival (RFS) in multivariable analysis (adjusted hazard ratio = 2.2 [1.6–3.0]). For stage II disease, 5‐year RFS for patients with MIR31HG outlier status was 49% compared to 77% for those with normal‐like expression. MIR31HG outlier status was associated with inferior outcome also within clinical high risk groups and within the poor prognostic CMS4‐mesenchymal gene expression subtype specifically. Preclinical models with MIR31HG outlier expression were characterized by reduced expression of MYC targets as well as elevated epithelial‐mesenchymal transition, TNF‐α/NFκB, TGF‐β, and IFN‐α/γ gene expression signatures, indicating cancer cell‐intrinsic properties resembling the CMS4 subgroup—associations which were recapitulated in patient biopsies. Moreover, the prognostic value of MIR31HG outlier status was independent of cytotoxic T lymphocyte and fibroblast infiltration. We here present evidence that MIR31HG expression provides clinical stratification beyond major gene expression phenotypes and tumor immune and stromal cell infiltration and propose a model where increased expression is an indicator of a cellular state conferring intrinsic invasive and/or immuno‐evasive capabilities., What's new? Expression of miR‐31 is associated with poor prognosis in colorectal cancer, but it's not known whether the prognostic value of the microRNA depends on consensus molecular subtypes (CMS). These authors looked at expression of MIR31HG, a long non‐coding RNA derived from the same primary transcript as miR‐31. They found that outlier expression of MIR31HG was associated with poor prognosis independent of CMS. Cells expressing high levels of MIR31HG also had strong enrichment for EMT and immune‐related gene sets. The authors propose that MIR31HG outlier expression is a biomarker for pro‐invasive and/or immunosuppressive characteristics, and poor prognosis.

KRAS mutation is a well-known marker for poor response to targeted treatment and patient prognosis in microsatellite stable (MSS) colorectal cancer (CRC). However, variation in clinical outcomes among patients wild-type for KRAS underlines that this is not a homogeneous population. Here, we evaluated the prognostic impact of KRAS alternative splicing in relation to mutation status in a single-hospital series of primary MSS CRCs (N = 258). Using splicing-sensitive microarrays and RNA sequencing, the relative expression of KRAS-4A versus KRAS-4B transcript variants was confirmed to be down-regulated in CRC compared to normal colonic mucosa (N = 41; p ≤ 0.001). This was independent of mutation status, however, gene set enrichment analysis revealed that the effect of splicing on KRAS signaling was specific to the KRAS wild-type subgroup, in which low relative KRAS-4A expression was associated with a higher level of KRAS signaling (p = 0.005). In concordance, the prognostic value of KRAS splicing was also dependent on mutation status, and for patients with Stage I-III KRAS wild-type MSS CRC, low relative KRAS-4A expression was associated with inferior overall survival (HR: 2.36, 95% CI: 1.07-5.18, p = 0.033), a result not found in mutant cases (pinteraction = 0.026). The prognostic association in the wild-type subgroup was independent of clinicopathological factors, including cancer stage in multivariable analysis (HR: 2.68, 95% CI: 1.18-6.09, p = 0.018). This suggests that KRAS has prognostic value beyond mutation status in MSS CRC, and highlights the importance of molecular heterogeneity in the clinically relevant KRAS wild-type subgroup.

Background Gene expression-based subtyping has the potential to form a new paradigm for stratified treatment of colorectal cancer. However, current frameworks are based on the transcriptomic profiles of primary tumors, and metastatic heterogeneity is a challenge. Here we aimed to develop a de novo metastasis-oriented framework. Methods In total, 829 transcriptomic profiles from patients with colorectal cancer were analyzed, including primary tumors, liver metastases, and non-malignant liver samples. High-resolution microarray gene expression profiling was performed of 283 liver metastases from 171 patients treated by hepatic resection, including multiregional and/or multi-metastatic samples from each of 47 patients. A single randomly selected liver metastasis sample from each patient was used for unsupervised subtype discovery by nonnegative matrix factorization, and a random forest prediction model was trained to classify multi-metastatic samples, as well as liver metastases from two independent series of 308 additional patients. Results Initial comparisons with non-malignant liver samples and primary colorectal tumors showed a highly variable degree of influence from the liver microenvironment in metastases, which contributed to inter-metastatic transcriptomic heterogeneity, but did not define subtype distinctions. The de novo liver metastasis subtype (LMS) framework recapitulated the main distinction between epithelial-like and mesenchymal-like tumors, with a strong immune and stromal component only in the latter. We also identified biologically distinct epithelial-like subtypes originating from different progenitor cell types. LMS1 metastases had several transcriptomic features of cancer aggressiveness, including secretory progenitor cell origin, oncogenic addictions, and microsatellite instability in a microsatellite stable background, as well as frequent RAS/TP53 co-mutations. The poor-prognostic association of LMS1 metastases was independent of mutation status, clinicopathological variables, and current subtyping frameworks (consensus molecular subtypes and colorectal cancer intrinsic subtypes). LMS1 was also the least heterogeneous subtype in comparisons of multiple metastases per patient, and tumor heterogeneity did not confound the prognostic value of LMS1. Conclusions We report the first large study of multi-metastatic gene expression profiling of colorectal cancer. The new metastasis-oriented subtyping framework showed potential for clinically relevant transcriptomic classification in the context of metastatic heterogeneity, and an LMS1 mini-classifier was constructed to facilitate prognostic stratification and further clinical testing.

Purpose:Molecular tumor heterogeneity may have important implications for the efficacy of targeted therapies in metastatic cancers. Inter-metastatic heterogeneity of sensitivity to anticancer agents has not been well explored in colorectal cancer.Experimental Design:We established a platform for ex vivo pharmacogenomic profiling of patient-derived organoids (PDO) from resected colorectal cancer liver metastases. Drug sensitivity testing (n = 40 clinically relevant agents) and gene expression profiling were performed on 39 metastases from 22 patients.Results:Three drug–response clusters were identified among the colorectal cancer metastases, based primarily on sensitivities to EGFR and/or MDM2 inhibition, and corresponding with RAS mutations and TP53 activity. Potentially effective therapies, including off-label use of drugs approved for other cancer types, could be nominated for eighteen patients (82%). Antimetabolites and targeted agents lacking a decisive genomic marker had stronger differential activity than most approved chemotherapies. We found limited intra-patient drug sensitivity heterogeneity between PDOs from multiple (2–5) liver metastases from each of ten patients. This was recapitulated at the gene expression level, with a highly proportional degree of transcriptomic and pharmacological variation. One PDO with a multi-drug resistance profile, including resistance to EGFR inhibition in a RAS-mutant background, showed sensitivity to MEK plus mTOR/AKT inhibition, corresponding with low-level PTEN expression.Conclusions:Intra-patient inter-metastatic pharmacological heterogeneity was not pronounced and ex vivo drug screening may identify novel treatment options for metastatic colorectal cancer. Variation in drug sensitivities was reflected at the transcriptomic level, suggesting potential to develop gene expression–based predictive signatures to guide experimental therapies.

Background: There is a need for improved selection of patients for adjuvant chemotherapy after resection of non-metastatic colorectal cancer (CRC). Regulator of chromosome condensation 2 (RCC2) is a potential prognostic biomarker. We report on the establishment of a robust protocol for RCC2 expression analysis and prognostic tumour biomarker evaluation in patients who did and did not receive adjuvant chemotherapy. Materials and Methods: RCC2 was analysed in 2916 primary CRCs from the QUASAR2 randomised trial and two single-hospital Norwegian series. A new protocol using fluorescent antibody staining and digital image analysis was optimised. Biomarker value for 5-year relapse-free survival was analysed in relation to tumour stage, adjuvant chemotherapy and the molecular markers microsatellite instability, KRAS/BRAFV600E/TP53 mutations and CDX2 expression. Results: Low RCC2 expression was scored in 41% of 2696 evaluable samples. Among patients with stage I–III CRC who had not received adjuvant chemotherapy, low RCC2 expression was an independent marker of inferior 5-year relapse-free survival in multivariable Cox models including clinicopathological factors and molecular markers (HR 1.45, 95% CI 1.09 to 1.94, p=0.012, N=521). RCC2 was not prognostic in patients who had received adjuvant chemotherapy, neither in QUASAR2 nor the pooled Norwegian series. The interaction between RCC2 and adjuvant chemotherapy for prediction of patient outcome was significant in stage III, and strongest among patients with microsatellite stable tumours (pinteraction=0.028). Conclusions: Low expression of RCC2 is a biomarker for poor prognosis in patients with stage I–III CRC and seems to be a predictive biomarker for effect of adjuvant chemotherapy.

Colorectal cancer commonly metastasizes into multiple liver foci, but intermetastatic heterogeneity remains poorly described. Here, we demonstrate that mutations of RAS/BRAF/TP53 are homogeneous within patients, while DNA copy number aberrations can vary greatly. RAS/BRAF/TP53 co‐mutations conferred a dismal prognosis, and co‐mutations combined with a high burden and heterogeneity of copy number aberrations identified patients with the poorest outcome., Hepatic resection is potentially curative for patients with colorectal liver metastases, but the treatment benefit varies. KRAS/NRAS (RAS)/TP53 co‐mutations are associated with a poor prognosis after resection, but there is large variation in patient outcome within the mutation groups, and genetic testing is currently not used to evaluate benefit from surgery. We have investigated the potential for improved prognostic stratification by combined biomarker analysis with DNA copy number aberrations (CNAs), and taking tumor heterogeneity into account. We determined the mutation status of RAS, BRAF V600, and TP53 in 441 liver lesions from 171 patients treated by partial hepatectomy for metastatic colorectal cancer. CNAs were profiled in 232 tumors from 67 of the patients. Mutations and high‐level amplifications of cancer‐critical genes, the latter including ERBB2 and EGFR, were predominantly homogeneous within patients. RAS/BRAF V600E and TP53 co‐mutations were associated with a poor patient outcome (hazard ratio, HR, 3.9, 95% confidence interval, CI, 1.3–11.1, P = 0.012) in multivariable analyses with clinicopathological variables. The genome‐wide CNA burden and intrapatient intermetastatic CNA heterogeneity varied within the mutation groups, and the CNA burden had prognostic associations in univariable analysis. Combined prognostic analyses of RAS/BRAF V600E/TP53 mutations and CNAs, either as a high CNA burden or high intermetastatic CNA heterogeneity, identified patients with a particularly poor outcome (co‐mutation/high CNA burden: HR 2.7, 95% CI 1.2–5.9, P = 0.013; co‐mutation/high CNA heterogeneity: HR 2.5, 95% CI 1.1–5.6, P = 0.022). In conclusion, DNA copy number profiling identified genomic and prognostic heterogeneity among patients with resectable colorectal liver metastases with co‐mutated RAS/BRAF V600E/TP53.

Background PARP inhibitors are active in various tumour types beyond BRCA-mutant cancers, but their activity and molecular correlates in colorectal cancer (CRC) are not well studied. Methods Mutations and genome-wide mutational patterns associated with homologous recombination deficiency (HRD) were investigated in 255 primary CRCs with whole-exome sequencing and/or DNA copy number data. Efficacy of five PARP inhibitors and their molecular correlates were evaluated in 93 CRC cell lines partly annotated with mutational-, DNA copy number-, and/or gene expression profiles. Post-treatment gene expression profiling and specific protein expression analyses were performed in two pairs of PARP inhibitor sensitive and resistant cell lines. Findings A subset of microsatellite stable (MSS) CRCs had truncating mutations in homologous recombination-related genes, but these were not associated with genomic signatures of HRD. Eight CRC cell lines (9%) were sensitive to PARP inhibition, but sensitivity was not predicted by HRD-related genomic and transcriptomic signatures. In contrast, drug sensitivity in MSS cell lines was strongly associated with TP53 wild-type status (odds ratio 15.7, p = 0.023) and TP53-related expression signatures. Increased downstream TP53 activity was among the primary response mechanisms, and TP53 inhibition antagonized the effect of PARP inhibitors. Wild-type TP53-mediated suppression of RAD51 was identified as a possible mechanism of action for sensitivity to PARP inhibition. Interpretation PARP inhibitors are active in a subset of CRC cell lines and preserved TP53 function may increase the likelihood of response.

Intratumor heterogeneity of colorectal cancers (CRCs) is manifested both at the genomic and epigenomic levels. Early genetic aberrations in carcinogenesis are clonal and present throughout the tumors, but less is known about the heterogeneity of the epigenetic CpG island methylator phenotype (CIMP). CIMP characterizes a subgroup of CRCs thought to originate from specific precursor lesions, and it is defined by widespread DNA methylation within promoter regions. In this work, we investigated CIMP in two to four multiregional samples from 30 primary tumors (n = 86 samples) using the consensus Weisenberger gene panel (CACNA1G, IGF2, NEUROG1, RUNX3 and SOCS1). Twenty‐nine of 30 tumors (97%) showed concordant CIMP status in all samples, and percent methylated reference (PMR) values of all five markers had higher intertumor than intratumor variation (P value = 1.5e−09). However, a third of the CIMP+ tumors exhibited discrepancies in methylation status in at least one of the five gene markers. To conclude, CIMP status was consistent within primary CRCs, and it is likely a clonal phenotype. However, spatial discordances of the individual genes suggest that large‐scale analysis of multiregional samples could be of interest for identifying CIMP markers that are robust to intratumor heterogeneity., What's new? Colorectal cancers (CRCs) exhibit significant intratumoral genetic and epigenetic heterogeneity. A subgroup of CRCs is characterized in particular by the epigenetic CpG island methylator phenotype (CIMP), though the extent to which CIMP contributes to intratumoral heterogeneity in these tumors is unknown. Here, investigation of CIMP in multiregional samples from primary CRCs shows that CIMP status is highly homogenous within tumors. In one‐third of CIMP‐positive primary CRCs, methylation status differed in at least one of five gene markers investigated. The findings suggest that inclusion of multiregional CRC samples could aid the development of more robust marker panels for CIMP assessment.

Key Points Diagnostic and relapse diffuse large B-cell lymphoma (DLBCL) biopsies reveal increased mutational burden/loss of heterozygosity in HLA-A. Serially sampled tumor biopsies provide insight into therapeutic targets and evolutionary divergence in relapsed/refractory DLBCL.

Colorectal cancers (CRCs) can be divided into four gene expression-based biologically distinct consensus molecular subtypes (CMS). This classification provides a potential framework for stratified treatment, but to identify novel CMS-drug associations, translation of the subtypes to pre-clinical models is essential. The currently available classifier is dependent on gene expression signals from the immune and stromal compartments of tumors and fails to identify the poor-prognostic CMS4-mesenchymal group in immortalized cell lines, patient-derived organoids and xenografts. To address this, we present a novel CMS classifier based on a filtered set of cancer cell-intrinsic, subtype-enriched gene expression markers. This new classifier, referred to as CMScaller, recapitulated the subtypes in both in vitro and in vivo models (551 in total). Importantly, by analyzing public drug response data from patient-derived xenografts and cell lines, we show that the subtypes are predictive of response to standard CRC drugs. CMScaller is available as an R package.

The number of molecularly stratified treatment options available to patients with colorectal cancer (CRC) is increasing, with a parallel rise in the use of biomarkers to guide prognostication and treatment decision-making. The increase in both the number of biomarkers and their use has resulted in a progressively complex situation, evident both from the extensive interactions between biomarkers and from their sometimes complex associations with patient prognosis and treatment benefit. Current and emerging biomarkers also reflect the genomic complexity of CRC, and include a wide range of aberrations such as point mutations, amplifications, fusions and hypermutator phenotypes, in addition to global gene expression subtypes. In this Review, we provide an overview of current and emerging clinically relevant biomarkers and their role in the management of patients with CRC, illustrating the intricacies of biomarker interactions and the growing treatment opportunities created by the availability of comprehensive molecular profiling. With the development of novel targeted therapies for patients with colorectal cancer, comes a wealth of new, and increasingly complex information on biomarkers. In this Review, the authors describe this increased complexity, with a focus on interactions between more than one biomarker and the implications of these interactions for patient management.

Background Accumulating evidence suggests immunomodulatory and context-dependent effects of TP53 mutations in cancer. We performed an exploratory analysis of the transcriptional, immunobiological and prognostic associations of TP53 mutations within the gene expression-based consensus molecular subtypes (CMSs) of colorectal cancer (CRC). Materials and methods In a single-hospital series of 401 stage I–IV primary CRCs, we sequenced the whole coding region of TP53 and analysed CMS-dependent transcriptional consequences of the mutations by gene expression profiling. Immunomodulatory associations were validated by multiplex, fluorescence-based immunohistochemistry of immune cell markers. Prognostic associations of TP53 mutations were analysed in an aggregated series of 635 patients classified according to CMS, including publicly available data from a French multicentre cohort (GSE39582). Results TP53 mutations were found in 60% of the CRCs. However, gene set enrichment analyses indicated that their transcriptional consequences varied among the CMSs and were most pronounced in CMS1-immune and CMS4-mesenchymal. Subtype specificity was primarily seen as an upregulation of gene sets reflecting cell cycle progression in CMS4 and a downregulation of T cell activity in CMS1. The subtype-dependent immunomodulatory associations were reinforced by significant depletion of several immune cell populations in mutated tumours compared with wild-type (wt) tumours exclusively in CMS1, including cytotoxic lymphocytes (adjusted p value in CMS1=0.002 and CMS2−4>0.9, Microenvironment Cell Populations (MCP)-counter algorithm). This was validated by immunohistochemistry-based quantification of tumour infiltrating CD8+ cells. Within CMS1, the immunomodulatory association of TP53 mutations was strongest among microsatellite stable (MSS) tumours, and this translated into a propensity for metastatic disease and poor prognostic value of the mutations specifically in the CMS1/MSS subtype (both series overall survival: TP53 mutation vs wt: HR 5.52, p=0.028). Conclusions Integration of TP53 mutation status with the CMS framework in primary CRC suggested subtype-dependent immunobiological associations with prognostic and potentially immunotherapeutic implications, warranting independent validation.