Transcriptome-wide gene expression outlier analysis pinpoints therapeutic vulnerabilities in colorectal cancer.

Multiple strategies are continuously being explored to expand the drug target repertoire in solid tumors. We devised a novel computational workflow for transcriptome-wide gene expression outlier analysis that allows the systematic identification of both overexpression and underexpression events in cancer cells. Here, it was applied to expression values obtained through RNA sequencing in 226 colorectal cancer (CRC) cell lines that were also characterized by whole-exome sequencing and microarray-based DNA methylation profiling. We found cell models displaying an abnormally high or low expression level for 3533 and 965 genes, respectively. Gene expression abnormalities that have been previously associated with clinically relevant features of CRC cell lines were confirmed. Moreover, by integrating multi-omics data, we identified both genetic and epigenetic alternations underlying outlier expression values. Importantly, our atlas of CRC gene expression outliers can guide the discovery of novel drug targets and biomarkers. As a proof of concept, we found that CRC cell lines lacking expression of the MTAP gene are sensitive to treatment with a PRMT5-MTA inhibitor (MRTX1719). Finally, other tumor types may also benefit from this approach.

Targeting the DNA Damage Response Pathways and Replication Stress in Colorectal Cancer.

PURPOSE: Genomic instability is a hallmark of cancer and targeting DNA damage response (DDR) is emerging as a promising therapeutic strategy in different solid tumors. The effectiveness of targeting DDR in colorectal cancer has not been extensively explored. EXPERIMENTAL DESIGN: We challenged 112 cell models recapitulating the genomic landscape of metastatic colorectal cancer with ATM, ATR, CHK1, WEE1, and DNA-PK inhibitors, in parallel with chemotherapeutic agents. We focused then on ATR inhibitors (ATRi) and, to identify putative biomarkers of response and resistance, we analyzed at multiple levels colorectal cancer models highly sensitive or resistant to these drugs. RESULTS: We found that around 30% of colorectal cancers, including those carrying KRAS and BRAF mutations and unresponsive to targeted agents, are sensitive to at least one DDR inhibitor. By investigating potential biomarkers of response to ATRi, we found that ATRi-sensitive cells displayed reduced phospho-RPA32 foci at basal level, while ATRi-resistant cells showed increased RAD51 foci formation in response to replication stress. Lack of ATM and RAD51C expression was associated with ATRi sensitivity. Analysis of mutational signatures and HRDetect score identified a subgroup of ATRi-sensitive models. Organoids derived from patients with metastatic colorectal cancer recapitulated findings obtained in cell lines. CONCLUSIONS: In conclusion, a subset of colorectal cancers refractory to current therapies could benefit from inhibitors of DDR pathways and replication stress. A composite biomarker involving phospho-RPA32 and RAD51 foci, lack of ATM and RAD51C expression, as well as analysis of mutational signatures could be used to identify colorectal cancers likely to respond to ATRi.

The heme synthesis-export system regulates the tricarboxylic acid cycle flux and oxidative phosphorylation.

Heme is an iron-containing porphyrin of vital importance for cell energetic metabolism. High rates of heme synthesis are commonly observed in proliferating cells. Moreover, the cell-surface heme exporter feline leukemia virus subgroup C receptor 1a (FLVCR1a) is overexpressed in several tumor types. However, the reasons why heme synthesis and export are enhanced in highly proliferating cells remain unknown. Here, we illustrate a functional axis between heme synthesis and heme export: heme efflux through the plasma membrane sustains heme synthesis, and implementation of the two processes down-modulates the tricarboxylic acid (TCA) cycle flux and oxidative phosphorylation. Conversely, inhibition of heme export reduces heme synthesis and promotes the TCA cycle fueling and flux as well as oxidative phosphorylation. These data indicate that the heme synthesis-export system modulates the TCA cycle and oxidative metabolism and provide a mechanistic basis for the observation that both processes are enhanced in cells with high-energy demand.

Werner Helicase Is a Synthetic-Lethal Vulnerability in Mismatch Repair-Deficient Colorectal Cancer Refractory to Targeted Therapies, Chemotherapy, and Immunotherapy.

Targeted therapies, chemotherapy, and immunotherapy are used to treat patients with mismatch repair-deficient (dMMR)/microsatellite instability-high (MSI-H) colorectal cancer. The clinical effectiveness of targeted therapy and chemotherapy is limited by resistance and drug toxicities, and about half of patients receiving immunotherapy have disease that is refractory to immune checkpoint inhibitors. Loss of Werner syndrome ATP-dependent helicase (WRN) is a synthetic lethality in dMMR/MSI-H cells. To inform the development of WRN as a therapeutic target, we performed WRN knockout or knockdown in 60 heterogeneous dMMR colorectal cancer preclinical models, demonstrating that WRN dependency is an almost universal feature and a robust marker for patient selection. Furthermore, models of resistance to clinically relevant targeted therapy, chemotherapy, and immunotherapy retain WRN dependency. These data show the potential of therapeutically targeting WRN in patients with dMMR/MSI-H colorectal cancer and support WRN as a therapeutic option for patients with dMMR/MSI-H cancers refractory to current treatment strategies. SIGNIFICANCE: We found that a large, diverse set of dMMR/MSI-H colorectal cancer preclinical models, including models of treatment-refractory disease, are WRN-dependent. Our results support WRN as a promising synthetic-lethal target in dMMR/MSI-H colorectal cancer tumors as a monotherapy or in combination with targeted agents, chemotherapy, or immunotherapy.This article is highlighted in the In This Issue feature, p. 1861.

EGFR Blockade Reverts Resistance to KRASG12C Inhibition in Colorectal Cancer.

Most patients with KRAS G12C-mutant non-small cell lung cancer (NSCLC) experience clinical benefit from selective KRASG12C inhibition, whereas patients with colorectal cancer bearing the same mutation rarely respond. To investigate the cause of the limited efficacy of KRASG12C inhibitors in colorectal cancer, we examined the effects of AMG510 in KRAS G12C colorectal cancer cell lines. Unlike NSCLC cell lines, KRAS G12C colorectal cancer models have high basal receptor tyrosine kinase (RTK) activation and are responsive to growth factor stimulation. In colorectal cancer lines, KRASG12C inhibition induces higher phospho-ERK rebound than in NSCLC cells. Although upstream activation of several RTKs interferes with KRASG12C blockade, we identify EGFR signaling as the dominant mechanism of colorectal cancer resistance to KRASG12C inhibitors. The combinatorial targeting of EGFR and KRASG12C is highly effective in colorectal cancer cells and patient-derived organoids and xenografts, suggesting a novel therapeutic strategy to treat patients with KRAS G12C colorectal cancer. SIGNIFICANCE: The efficacy of KRASG12C inhibitors in NSCLC and colorectal cancer is lineage-specific. RTK dependency and signaling rebound kinetics are responsible for sensitivity or resistance to KRASG12C inhibition in colorectal cancer. EGFR and KRASG12C should be concomitantly inhibited to overcome resistance to KRASG12C blockade in colorectal tumors.See related commentary by Koleilat and Kwong, p. 1094.This article is highlighted in the In This Issue feature, p. 1079.

A Subset of Colorectal Cancers with Cross-Sensitivity to Olaparib and Oxaliplatin.

PURPOSE: Defects in the homologous recombination (HR) repair pathway are of clinical interest due to sensitivity of HR-deficient cells to PARP inhibitors. We were interested in defining PARP vulnerability in patients with metastatic colorectal cancer (mCRC) carrying KRAS and BRAF mutations who display poor prognosis, have limited therapeutic options, and represent an unmet clinical need. EXPERIMENTAL DESIGN: We tested colorectal cancer cell lines, patient-derived organoids (PDO), and patient-derived xenografts (PDX) enriched for KRAS and BRAF mutations for sensitivity to the PARP inhibitor olaparib, and the chemotherapeutic agents oxaliplatin and 5-fluorouracil (5-FU). Genomic profiles and DNA repair proficiency of colorectal cancer models were compared with pharmacologic response. RESULTS: Thirteen of 99 (around 13%) colorectal cancer cell lines were highly sensitive to clinically active concentrations of olaparib and displayed functional deficiency in HR. Response to PARP blockade was positively correlated with sensitivity to oxaliplatin in colorectal cancer cell lines as well as patient-derived organoids. Treatment of PDXs with olaparib impaired tumor growth and maintenance therapy with PARP blockade after initial oxaliplatin response delayed disease progression in mice. CONCLUSIONS: These results indicate that a colorectal cancer subset characterized by poor prognosis and limited therapeutic options is vulnerable to PARP inhibition and suggest that PDO-based drug-screening assays can be used to identify patients with colorectal cancer likely to benefit from olaparib. As patients with mCRC almost invariably receive therapies based on oxaliplatin, "maintenance" treatment with PARP inhibitors warrants further clinical investigation.

Patient-Derived Xenografts and Matched Cell Lines Identify Pharmacogenomic Vulnerabilities in Colorectal Cancer.

PURPOSE: Patient-derived xenograft (PDX) models accurately recapitulate the tumor of origin in terms of histopathology, genomic landscape, and therapeutic response, but some limitations due to costs associated with their maintenance and restricted amenability for large-scale screenings still exist. To overcome these issues, we established a platform of 2D cell lines (xeno-cell lines, XL), derived from PDXs of colorectal cancer with matched patient germline gDNA available. EXPERIMENTAL DESIGN: Whole-exome and transcriptome sequencing analyses were performed. Biomarkers of response and resistance to anti-HER therapy were annotated. Dependency on the WRN helicase gene was assessed in MSS, MSI-H, and MSI-like XLs using a reverse genetics functional approach. RESULTS: XLs recapitulated the entire spectrum of colorectal cancer transcriptional subtypes. Exome and RNA-seq analyses delineated several molecular biomarkers of response and resistance to EGFR and HER2 blockade. Genotype-driven responses observed in vitro in XLs were confirmed in vivo in the matched PDXs. MSI-H models were dependent upon WRN gene expression, while loss of WRN did not affect MSS XLs growth. Interestingly, one MSS XL with transcriptional MSI-like traits was sensitive to WRN depletion. CONCLUSIONS: The XL platform represents a preclinical tool for functional gene validation and proof-of-concept studies to identify novel druggable vulnerabilities in colorectal cancer.

Evolving neoantigen profiles in colorectal cancers with DNA repair defects.

BACKGROUND: Neoantigens that arise as a consequence of tumor-specific mutations can be recognized by T lymphocytes leading to effective immune surveillance. In colorectal cancer (CRC) and other tumor types, a high number of neoantigens is associated with patient response to immune therapies. The molecular processes governing the generation of neoantigens and their turnover in cancer cells are poorly understood. We exploited CRC as a model system to understand how alterations in DNA repair pathways modulate neoantigen profiles over time. METHODS: We performed whole exome sequencing (WES) and RNA sequencing (RNAseq) in CRC cell lines, in vitro and in vivo, and in CRC patient-derived xenografts (PDXs) to track longitudinally genomic profiles, clonal evolution, mutational signatures, and predicted neoantigens. RESULTS: The majority of CRC models showed remarkably stable mutational and neoantigen profiles; however, those carrying defects in DNA repair genes continuously diversified. Rapidly evolving and evolutionary stable CRCs displayed characteristic genomic signatures and transcriptional profiles. Downregulation of molecules implicated in antigen presentation occurred selectively in highly mutated and rapidly evolving CRC. CONCLUSIONS: These results indicate that CRCs carrying alterations in DNA repair pathways display dynamic neoantigen patterns that fluctuate over time. We define CRC subsets characterized by slow and fast evolvability and link this phenotype to downregulation of antigen-presenting cellular mechanisms. Longitudinal monitoring of the neoantigen landscape could be relevant in the context of precision medicine.

Discovery of methylated circulating DNA biomarkers for comprehensive non-invasive monitoring of treatment response in metastatic colorectal cancer.

OBJECTIVE: Mutations in cell-free circulating DNA (cfDNA) have been studied for tracking disease relapse in colorectal cancer (CRC). This approach requires personalised assay design due to the lack of universally mutated genes. In contrast, early methylation alterations are restricted to defined genomic loci allowing comprehensive assay design for population studies. Our objective was to identify cancer-specific methylated biomarkers which could be measured longitudinally in cfDNA (liquid biopsy) to monitor therapeutic outcome in patients with metastatic CRC (mCRC). DESIGN: Genome-wide methylation microarrays of CRC cell lines (n=149) identified five cancer-specific methylated loci (EYA4, GRIA4, ITGA4, MAP3K14-AS1, MSC). Digital PCR assays were employed to measure methylation of these genes in tumour tissue DNA (n=82) and cfDNA from patients with mCRC (n=182). Plasma longitudinal assessment was performed in a patient subset treated with chemotherapy or targeted therapy. RESULTS: Methylation in at least one marker was detected in all tumour tissue samples and in 156 mCRC patient cfDNA samples (85.7%). Plasma marker prevalence was 71.4% for EYA4, 68.5% for GRIA4, 69.7% for ITGA4, 69.1% for MAP3K14-AS1% and 65.1% for MSC. Dynamics of methylation markers was not affected by treatment type and correlated with objective tumour response and progression-free survival. CONCLUSION: This five-gene methylation panel can be used to circumvent the absence of patient-specific mutations for monitoring tumour burden dynamics in liquid biopsy under different therapeutic regimens. This method might be proposed for assessing pharmacodynamics in clinical trials or when conventional imaging has limitations.

Molecular Landscape of Acquired Resistance to Targeted Therapy Combinations in BRAF-Mutant Colorectal Cancer.

Although recent clinical trials of BRAF inhibitor combinations have demonstrated improved efficacy in BRAF-mutant colorectal cancer, emergence of acquired resistance limits clinical benefit. Here, we undertook a comprehensive effort to define mechanisms underlying drug resistance with the goal of guiding development of therapeutic strategies to overcome this limitation. We generated a broad panel of BRAF-mutant resistant cell line models across seven different clinically relevant drug combinations. Combinatorial drug treatments were able to abrogate ERK1/2 phosphorylation in parental-sensitive cells, but not in their resistant counterparts, indicating that resistant cells escaped drug treatments through one or more mechanisms leading to biochemical reactivation of the MAPK signaling pathway. Genotyping of resistant cells identified gene amplification of EGFR, KRAS, and mutant BRAF, as well as acquired mutations in KRAS, EGFR, and MAP2K1 These mechanisms were clinically relevant, as we identified emergence of a KRAS G12C mutation and increase of mutant BRAF V600E allele frequency in the circulating tumor DNA of a patient at relapse from combined treatment with BRAF and MEK inhibitors. To identify therapeutic combinations capable of overcoming drug resistance, we performed a systematic assessment of candidate therapies across the panel of resistant cell lines. Independent of the molecular alteration acquired upon drug pressure, most resistant cells retained sensitivity to vertical MAPK pathway suppression when combinations of ERK, BRAF, and EGFR inhibitors were applied. These therapeutic combinations represent promising strategies for future clinical trials in BRAF-mutant colorectal cancer. Cancer Res; 76(15); 4504-15. ©2016 AACR.

The EGFR-specific antibody cetuximab combined with chemotherapy triggers immunogenic cell death.

Cetuximab is a monoclonal antibody that is effective in the treatment of metastatic colorectal cancer (mCRC). Cetuximab blocks epidermal growth factor receptor (EGFR)-ligand interaction and inhibits downstream RAS-ERK activation. However, only some activating mutations in RAS affect cetuximab efficacy, and it is not clear what else mediates treatment success. Here we hypothesized that cetuximab induces immunogenic cell death (ICD) that activates a potent antitumor response. We found that cetuximab, in combination with chemotherapy, fostered ICD in CRC cells, which we measured via the endoplasmic reticulum (ER) stress response and an increase in phagocytosis by dendritic cells. ICD induction depended on the mutational status of the EGFR signaling pathway and on the inhibition of the splicing of X-box binding protein 1 (XBP1), an unfolded protein response (UPR) mediator. We confirmed the enhanced immunogenicity elicited by cetuximab in a mouse model of human EGFR-expressing CRC. Overall, we demonstrate a new, immune-related mechanism of action of cetuximab that may help to tailor personalized medicine.

Oncogenic KRAS sensitizes premalignant, but not malignant cells, to Noxa-dependent apoptosis through the activation of the MEK/ERK pathway.

KRAS is mutated in about 20-25% of all human cancers and especially in pancreatic, lung and colorectal tumors. Oncogenic KRAS stimulates several pro-survival pathways, but it also triggers the trans-activation of pro-apoptotic genes. In our work, we show that G13D mutations of KRAS activate the MAPK pathway, and ERK2, but not ERK1, up-regulates Noxa basal levels. Accordingly, premalignant epithelial cells are sensitized to various cytotoxic compounds in a Noxa-dependent manner. In contrast to these findings, colorectal cancer cell sensitivity to treatment is independent of KRAS status and Noxa levels are not up-regulated in the presence of mutated KRAS despite the fact that ERK2 still promotes Noxa expression. We therefore speculated that other survival pathways are counteracting the pro-apoptotic effect of mutated KRAS and found that the inhibition of AKT restores sensitivity to treatment, especially in presence of oncogenic KRAS. In conclusion, our work suggests that the pharmacological inhibition of the pathways triggered by mutated KRAS could also switch off its oncogene-activated pro-apoptotic stimulation. On the contrary, the combination of chemotherapy to inhibitors of specific pro-survival pathways, such as the one controlled by AKT, could enhance treatment efficacy by exploiting the pro-death stimulation derived by oncogene activation.

The molecular landscape of colorectal cancer cell lines unveils clinically actionable kinase targets.

The development of molecularly targeted anticancer agents relies on large panels of tumour-specific preclinical models closely recapitulating the molecular heterogeneity observed in patients. Here we describe the mutational and gene expression analyses of 151 colorectal cancer (CRC) cell lines. We find that the whole spectrum of CRC molecular and transcriptional subtypes, previously defined in patients, is represented in this cell line compendium. Transcriptional outlier analysis identifies RAS/BRAF wild-type cells, resistant to EGFR blockade, functionally and pharmacologically addicted to kinase genes including ALK, FGFR2, NTRK1/2 and RET. The same genes are present as expression outliers in CRC patient samples. Genomic rearrangements (translocations) involving the ALK and NTRK1 genes are associated with the overexpression of the corresponding proteins in CRC specimens. The approach described here can be used to pinpoint CRCs with exquisite dependencies to individual kinases for which clinically approved drugs are already available.

RAF suppression synergizes with MEK inhibition in KRAS mutant cancer cells.

KRAS is the most frequently mutated oncogene in human cancer, yet no therapies are available to treat KRAS mutant cancers. We used two independent reverse genetic approaches to identify components of the RAS-signaling pathways required for growth of KRAS mutant tumors. Small interfering RNA (siRNA) screening of 37 KRAS mutant colorectal cancer cell lines showed that RAF1 suppression was synthetic lethal with MEK inhibition. An unbiased kinome short hairpin RNA (shRNA)-based screen confirmed this synthetic lethal interaction in colorectal as well as in lung cancer cells bearing KRAS mutations. Compounds targeting RAF kinases can reverse resistance to the MEK inhibitor selumetinib. MEK inhibition induces RAS activation and BRAF-RAF1 dimerization and sustains MEK-ERK signaling, which is responsible for intrinsic resistance to selumetinib. Prolonged dual blockade of RAF and MEK leads to persistent ERK suppression and efficiently induces apoptosis. Our data underlie the relevance of developing combinatorial regimens of drugs targeting the RAF-MEK pathway in KRAS mutant tumors.

BRAF V600E is a determinant of sensitivity to proteasome inhibitors.

A critical step toward defining tailored therapy in patients with cancer is the identification of genetic interactions that may impair-or boost-the efficacy of selected therapeutic approaches. Cell models able to recapitulate combinations of genetic aberrations are important to find drug-genotype interactions poorly affected by the heterogeneous genetics of human tumors. In order to identify novel pharmacogenomic relationships, we employed an isogenic cell panel that reconstructs cancer genetic scenarios. We screened a library of 43 compounds in human hTERT-HME1 epithelial cells in which PTEN or RB1 were silenced in combination with the targeted knockin of cancer-associated mutations in EGFR, KRAS, BRAF, or PIK3CA oncogenes. Statistical analysis and clustering algorithms were applied to display similar drug response profiles and mutation-specific patterns of activity. From the screen, we discovered that proteasome inhibitors show selectivity toward BRAF V600E-mutant cells, irrespective of PTEN or RB1 expression. Preferential targeting of BRAF-mutant cells by proteasome inhibitors was corroborated in a second BRAF V600E isogenic model, as well as in a panel of colorectal cancer cell lines by the use of the proteasome inhibitor carfilzomib. Notably, carfilzomib also showed striking in vivo activity in a BRAF-mutant human colorectal cancer xenograft model. Vulnerability to proteasome inhibitors is dependent on persistent BRAF signaling, because BRAF V600E blockade by PLX4720 reversed sensitivity to carfilzomib in BRAF-mutant colorectal cancer cells. Our findings indicated that proteasome inhibition might represent a valuable targeting strategy in BRAF V600E-mutant colorectal tumors.

Deregulation of the PI3K and KRAS signaling pathways in human cancer cells determines their response to everolimus.

Personalized cancer medicine is based on the concept that targeted therapies are effective on subsets of patients whose tumors carry specific molecular alterations. Several mammalian target of rapamycin (mTOR) inhibitors are in preclinical or clinical trials for cancers, but the molecular basis of sensitivity or resistance to these inhibitors among patients is largely unknown. Here we have identified oncogenic variants of phosphoinositide-3-kinase, catalytic, alpha polypeptide (PIK3CA) and KRAS as determinants of response to the mTOR inhibitor everolimus. Human cancer cells carrying alterations in the PI3K pathway were responsive to everolimus, both in vitro and in vivo, except when KRAS mutations occurred concomitantly or were exogenously introduced. In human cancer cells with mutations in both PIK3CA and KRAS, genetic ablation of mutant KRAS reinstated response to the drug. Consistent with these data, PIK3CA mutant cells, but not KRAS mutant cells, displayed everolimus-sensitive translation. Importantly, in a cohort of metastatic cancer patients, the presence of oncogenic KRAS mutations was associated with lack of benefit after everolimus therapy. Thus, our results demonstrate that alterations in the KRAS and PIK3CA genes may represent biomarkers to optimize treatment of patients with mTOR inhibitors.

Replacement of normal with mutant alleles in the genome of normal human cells unveils mutation-specific drug responses.

Mutations in oncogenes and tumor suppressor genes are responsible for tumorigenesis and represent favored therapeutic targets in oncology. We exploited homologous recombination to knock-in individual cancer mutations in the genome of nontransformed human cells. Sequential introduction of multiple mutations was also achieved, demonstrating the potential of this strategy to construct tumor progression models. Knock-in cells displayed allele-specific activation of signaling pathways and mutation-specific phenotypes different from those obtainable by ectopic oncogene expression. Profiling of a library of pharmacological agents on the mutated cells showed striking sensitivity or resistance phenotypes to pathway-targeted drugs, often matching those of tumor cells carrying equivalent cancer mutations. Thus, knock-in of single or multiple cancer alleles provides a pharmacogenomic platform for the rational design of targeted therapies.