Phosphoproteomics of patient-derived xenografts identifies targets and markers associated with sensitivity and resistance to EGFR blockade in colorectal cancer.

Epidermal growth factor receptor (EGFR) is a well-exploited therapeutic target in metastatic colorectal cancer (mCRC). Unfortunately, not all patients benefit from current EGFR inhibitors. Mass spectrometry-based proteomics and phosphoproteomics were performed on 30 genomically and pharmacologically characterized mCRC patient-derived xenografts (PDXs) to investigate the molecular basis of response to EGFR blockade and identify alternative drug targets to overcome resistance. Both the tyrosine and global phosphoproteome as well as the proteome harbored distinctive response signatures. We found that increased pathway activity related to mitogen-activated protein kinase (MAPK) inhibition and abundant tyrosine phosphorylation of cell junction proteins, such as CXADR and CLDN1/3, in sensitive tumors, whereas epithelial-mesenchymal transition and increased MAPK and AKT signaling were more prevalent in resistant tumors. Furthermore, the ranking of kinase activities in single samples confirmed the driver activity of ERBB2, EGFR, and MET in cetuximab-resistant tumors. This analysis also revealed high kinase activity of several members of the Src and ephrin kinase family in 2 CRC PDX models with genomically unexplained resistance. Inhibition of these hyperactive kinases, alone or in combination with cetuximab, resulted in growth inhibition of ex vivo PDX-derived organoids and in vivo PDXs. Together, these findings highlight the potential value of phosphoproteomics to improve our understanding of anti-EGFR treatment and response prediction in mCRC and bring to the forefront alternative drug targets in cetuximab-resistant tumors.

Synthetic Lethal Interaction with BCL-XL Blockade Deepens Response to Cetuximab in Patient-Derived Models of Metastatic Colorectal Cancer.

PURPOSE: Approximately 20% of patients with RAS wild-type metastatic colorectal cancer (mCRC) experience objective responses to the anti-EGFR antibody cetuximab, but disease eradication is seldom achieved. The extent of tumor shrinkage correlates with long-term outcome. We aimed to find rational combinations that potentiate cetuximab efficacy by disrupting adaptive dependencies on antiapoptotic molecules (BCL2, BCL-XL, MCL1). EXPERIMENTAL DESIGN: Experiments were conducted in patient-derived xenografts (PDX) and organoids (PDXO). Apoptotic priming was analyzed by BH3 profiling. Proapoptotic and antiapoptotic protein complexes were evaluated by co-immunoprecipitation and electroluminescence sandwich assays. The effect of combination therapies was assessed by caspase activation in PDXOs and by monitoring PDX growth. RESULTS: A population trial in 314 PDX cohorts, established from as many patients, identified 46 models (14.6%) with appreciable (>50% tumor shrinkage) but incomplete response to cetuximab. From these models, 14 PDXOs were derived. Cetuximab primed cells for apoptosis, but only concomitant blockade of BCL-XL precipitated cell death. Mechanistically, exposure to cetuximab induced upregulation of the proapoptotic protein BIM and its sequestration by BCL-XL. Inhibition of BCL-XL resulted in displacement of BIM, which was not buffered by MCL1 and thereby became competent to induce apoptosis. In five PDX models, combination of cetuximab and a selective BCL-XL inhibitor triggered apoptosis and led to more pronounced tumor regressions and longer time to relapse after treatment discontinuation than cetuximab alone. CONCLUSIONS: In mCRC tumors that respond to cetuximab, antibody treatment confers a synthetic-lethal dependency on BCL-XL. Targeting this dependency unleashes apoptosis and increases the depth of response to cetuximab.

Towards precision oncology with patient-derived xenografts.

Under the selective pressure of therapy, tumours dynamically evolve multiple adaptive mechanisms that make static interrogation of genomic alterations insufficient to guide treatment decisions. Clinical research does not enable the assessment of how various regulatory circuits in tumours are affected by therapeutic insults over time and space. Likewise, testing different precision oncology approaches informed by composite and ever-changing molecular information is hard to achieve in patients. Therefore, preclinical models that incorporate the biology and genetics of human cancers, facilitate analyses of complex variables and enable adequate population throughput are needed to pinpoint randomly distributed response predictors. Patient-derived xenograft (PDX) models are dynamic entities in which cancer evolution can be monitored through serial propagation in mice. PDX models can also recapitulate interpatient diversity, thus enabling the identification of response biomarkers and therapeutic targets for molecularly defined tumour subgroups. In this Review, we discuss examples from the past decade of the use of PDX models for precision oncology, from translational research to drug discovery. We elaborate on how and to what extent preclinical observations in PDX models have confirmed and/or anticipated findings in patients. Finally, we illustrate emerging methodological efforts that could broaden the application of PDX models by honing their predictive accuracy or improving their versatility.

Chromatin Velocity reveals epigenetic dynamics by single-cell profiling of heterochromatin and euchromatin.

Recent efforts have succeeded in surveying open chromatin at the single-cell level, but high-throughput, single-cell assessment of heterochromatin and its underlying genomic determinants remains challenging. We engineered a hybrid transposase including the chromodomain (CD) of the heterochromatin protein-1α (HP-1α), which is involved in heterochromatin assembly and maintenance through its binding to trimethylation of the lysine 9 on histone 3 (H3K9me3), and developed a single-cell method, single-cell genome and epigenome by transposases sequencing (scGET-seq), that, unlike single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq), comprehensively probes both open and closed chromatin and concomitantly records the underlying genomic sequences. We tested scGET-seq in cancer-derived organoids and human-derived xenograft (PDX) models and identified genetic events and plasticity-driven mechanisms contributing to cancer drug resistance. Next, building upon the differential enrichment of closed and open chromatin, we devised a method, Chromatin Velocity, that identifies the trajectories of epigenetic modifications at the single-cell level. Chromatin Velocity uncovered paths of epigenetic reorganization during stem cell reprogramming and identified key transcription factors driving these developmental processes. scGET-seq reveals the dynamics of genomic and epigenetic landscapes underlying any cellular processes.

Molecular Subtyping Combined with Biological Pathway Analyses to Study Regorafenib Response in Clinically Relevant Mouse Models of Colorectal Cancer.

PURPOSE: Regorafenib (REG) is approved for the treatment of metastatic colorectal cancer, but has modest survival benefit and associated toxicities. Robust predictive/early response biomarkers to aid patient stratification are outstanding. We have exploited biological pathway analyses in a patient-derived xenograft (PDX) trial to study REG response mechanisms and elucidate putative biomarkers. EXPERIMENTAL DESIGN: Molecularly subtyped PDXs were annotated for REG response. Subtyping was based on gene expression (CMS, consensus molecular subtype) and copy-number alteration (CNA). Baseline tumor vascularization, apoptosis, and proliferation signatures were studied to identify predictive biomarkers within subtypes. Phospho-proteomic analysis was used to identify novel classifiers. Supervised RNA sequencing analysis was performed on PDXs that progressed, or did not progress, following REG treatment. RESULTS: Improved REG response was observed in CMS4, although intra-subtype response was variable. Tumor vascularity did not correlate with outcome. In CMS4 tumors, reduced proliferation and higher sensitivity to apoptosis at baseline correlated with response. Reverse phase protein array (RPPA) analysis revealed 4 phospho-proteomic clusters, one of which was enriched with non-progressor models. A classification decision tree trained on RPPA- and CMS-based assignments discriminated non-progressors from progressors with 92% overall accuracy (97% sensitivity, 67% specificity). Supervised RNA sequencing revealed that higher basal EPHA2 expression is associated with REG resistance. CONCLUSIONS: Subtype classification systems represent canonical "termini a quo" (starting points) to support REG biomarker identification, and provide a platform to identify resistance mechanisms and novel contexts of vulnerability. Incorporating functional characterization of biological systems may optimize the biomarker identification process for multitargeted kinase inhibitors.

Implementing Systems Modelling and Molecular Imaging to Predict the Efficacy of BCL-2 Inhibition in Colorectal Cancer Patient-Derived Xenograft Models.

Resistance to chemotherapy often results from dysfunctional apoptosis, however multiple proteins with overlapping functions regulate this pathway. We sought to determine whether an extensively validated, deterministic apoptosis systems model, 'DR_MOMP', could be used as a stratification tool for the apoptosis sensitiser and BCL-2 antagonist, ABT-199 in patient-derived xenograft (PDX) models of colorectal cancer (CRC). Through quantitative profiling of BCL-2 family proteins, we identified two PDX models which were predicted by DR_MOMP to be sufficiently sensitive to 5-fluorouracil (5-FU)-based chemotherapy (CRC0344), or less responsive to chemotherapy but sensitised by ABT-199 (CRC0076). Treatment with ABT-199 significantly improved responses of CRC0076 PDXs to 5-FU-based chemotherapy, but showed no sensitisation in CRC0344 PDXs, as predicted from systems modelling. 18F-Fluorodeoxyglucose positron emission tomography/computed tomography (18F-FDG-PET/CT) scans were performed to investigate possible early biomarkers of response. In CRC0076, a significant post-treatment decrease in mean standard uptake value was indeed evident only in the combination treatment group. Radiomic CT feature analysis of pre-treatment images in CRC0076 and CRC0344 PDXs identified features which could phenotypically discriminate between models, but were not predictive of treatment responses. Collectively our data indicate that systems modelling may identify metastatic (m)CRC patients benefitting from ABT-199, and that 18F-FDG-PET could independently support such predictions.

Colorectal cancer residual disease at maximal response to EGFR blockade displays a druggable Paneth cell-like phenotype.

Blockade of epidermal growth factor receptor (EGFR) causes tumor regression in some patients with metastatic colorectal cancer (mCRC). However, residual disease reservoirs typically remain even after maximal response to therapy, leading to relapse. Using patient-derived xenografts (PDXs), we observed that mCRC cells surviving EGFR inhibition exhibited gene expression patterns similar to those of a quiescent subpopulation of normal intestinal secretory precursors with Paneth cell characteristics. Compared with untreated tumors, these pseudodifferentiated tumor remnants had reduced expression of genes encoding EGFR-activating ligands, enhanced activity of human epidermal growth factor receptor 2 (HER2) and HER3, and persistent signaling along the phosphatidylinositol 3-kinase (PI3K) pathway. Clinically, properties of residual disease cells from the PDX models were detected in lingering tumors of responsive patients and in tumors of individuals who had experienced early recurrence. Mechanistically, residual tumor reprogramming after EGFR neutralization was mediated by inactivation of Yes-associated protein (YAP), a master regulator of intestinal epithelium recovery from injury. In preclinical trials, Pan-HER antibodies minimized residual disease, blunted PI3K signaling, and induced long-term tumor control after treatment discontinuation. We found that tolerance to EGFR inhibition is characterized by inactivation of an intrinsic lineage program that drives both regenerative signaling during intestinal repair and EGFR-dependent tumorigenesis. Thus, our results shed light on CRC lineage plasticity as an adaptive escape mechanism from EGFR-targeted therapy and suggest opportunities to preemptively target residual disease.

Systems analysis of protein signatures predicting cetuximab responses in KRAS, NRAS, BRAF and PIK3CA wild-type patient-derived xenograft models of metastatic colorectal cancer.

Antibodies targeting the human epidermal growth factor receptor (EGFR) are used for the treatment of RAS wild-type metastatic colorectal cancer. A significant proportion of patients remains unresponsive to this therapy. Here, we performed a reverse-phase protein array-based (phospho)protein analysis of 63 KRAS, NRAS, BRAF and PIK3CA wild-type metastatic CRC tumours. Responses of tumours to anti-EGFR therapy with cetuximab were recorded in patient-derived xenograft (PDX) models. Unsupervised hierarchical clustering of pretreatment tumour tissue identified three clusters, of which Cluster C3 was exclusively composed of responders. Clusters C1 and C2 exhibited mixed responses. None of the three protein clusters exhibited a significant correlation with transcriptome-based subtypes. Analysis of protein signatures across all PDXs identified 14 markers that discriminated cetuximab-sensitive and cetuximab-resistant tumours: PDK1 (S241), caspase-8, Shc (Y317), Stat3 (Y705), p27, GSK-3ß (S9), HER3, PKC-α (S657), EGFR (Y1068), Akt (S473), S6 ribosomal protein (S240/244), HER3 (Y1289), NF-κB-p65 (S536) and Gab-1 (Y627). Least absolute shrinkage and selection operator and binominal logistic regression analysis delivered refined protein signatures for predicting response to cetuximab. (Phospo-)protein analysis of matched pretreated and posttreated models furthermore showed significant reduction of Gab-1 (Y627) and GSK-3ß (S9) exclusively in responding models, suggesting novel targets for treatment.

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.

The genomic landscape of response to EGFR blockade in colorectal cancer.

Colorectal cancer is the third most common cancer worldwide, with 1.2 million patients diagnosed annually. In late-stage colorectal cancer, the most commonly used targeted therapies are the monoclonal antibodies cetuximab and panitumumab, which prevent epidermal growth factor receptor (EGFR) activation. Recent studies have identified alterations in KRAS and other genes as likely mechanisms of primary and secondary resistance to anti-EGFR antibody therapy. Despite these efforts, additional mechanisms of resistance to EGFR blockade are thought to be present in colorectal cancer and little is known about determinants of sensitivity to this therapy. To examine the effect of somatic genetic changes in colorectal cancer on response to anti-EGFR antibody therapy, here we perform complete exome sequence and copy number analyses of 129 patient-derived tumour grafts and targeted genomic analyses of 55 patient tumours, all of which were KRAS wild-type. We analysed the response of tumours to anti-EGFR antibody blockade in tumour graft models and in clinical settings and functionally linked therapeutic responses to mutational data. In addition to previously identified genes, we detected mutations in ERBB2, EGFR, FGFR1, PDGFRA, and MAP2K1 as potential mechanisms of primary resistance to this therapy. Novel alterations in the ectodomain of EGFR were identified in patients with acquired resistance to EGFR blockade. Amplifications and sequence changes in the tyrosine kinase receptor adaptor gene IRS2 were identified in tumours with increased sensitivity to anti-EGFR therapy. Therapeutic resistance to EGFR blockade could be overcome in tumour graft models through combinatorial therapies targeting actionable genes. These analyses provide a systematic approach to evaluating response to targeted therapies in human cancer, highlight new mechanisms of responsiveness to anti-EGFR therapies, and delineate new avenues for intervention in managing colorectal cancer.

Sustained Inhibition of HER3 and EGFR Is Necessary to Induce Regression of HER2-Amplified Gastrointestinal Carcinomas.

PURPOSE: Preclinical studies in HER2-amplified gastrointestinal cancer models have shown that cotargeting HER2 with a monoclonal antibody and a small molecule is superior to monotherapy with either inhibitor, but the underlying cooperative mechanisms remain unexplored. We investigated the molecular underpinnings of this synergy to identify key vulnerabilities susceptible to alternative therapeutic opportunities. EXPERIMENTAL DESIGN: The phosphorylation/activation of HER2, HER3, EGFR (HER receptors), and downstream transducers was evaluated in HER2-overexpressing colorectal and gastric cancer cell lines by Western blotting and/or multiplex phosphoproteomics. The in vivo outcome of antibody-mediated HER2 blockade by trastuzumab, reversible HER2 inhibition by lapatinib, and irreversible HER2 inhibition by afatinib was assessed in patient-derived tumorgrafts and cell-line xenografts by monitoring tumor growth curves and by using antibody-based proximity assays. RESULTS: Trastuzumab monotherapy reduced HER3 phosphorylation, with minor consequences on downstream transducers. Lapatinib alone acutely inhibited all HER receptors and effectors but led to delayed rephosphorylation of HER3 and EGFR and partial restoration of ERK and AKT activity. When combined with lapatinib, trastuzumab prevented HER3/EGFR reactivation and caused prolonged inhibition of ERK/AKT. Afatinib alone was also very effective in counteracting the reinstatement of HER3, EGFR, and downstream signaling activation. In vivo, the combination of trastuzumab and lapatinib-or, importantly, monotherapy with afatinib-resulted in overt tumor shrinkage. CONCLUSIONS: Only prolonged inhibition of HER3 and EGFR, achievable by dual blockade with trastuzumab and lapatinib or irreversible HER2 inhibition by single-agent afatinib, led to regression of HER2-amplified gastrointestinal carcinomas. Clin Cancer Res; 21(24); 5519-31. ©2015 AACR.

IGF2 is an actionable target that identifies a distinct subpopulation of colorectal cancer patients with marginal response to anti-EGFR therapies.

Among patients with colorectal cancer who benefit from therapy targeted to the epidermal growth factor receptor (EGFR), stable disease (SD) occurs more frequently than massive regressions. Exploring the mechanisms of this incomplete sensitivity to devise more efficacious treatments will likely improve patients' outcomes. We tested therapies tailored around hypothesis-generating molecular features in patient-derived xenografts ("xenopatients"), which originated from 125 independent samples that did not harbor established resistance-conferring mutations. Samples from xenopatients that responded to cetuximab, an anti-EGFR agent, with disease stabilization displayed high levels of EGFR family ligands and receptors, indicating high EGFR pathway activity. Five of 21 SD models (23.8%) characterized by particularly high expression of EGFR and EGFR family members regressed after intensified EGFR blockade by cetuximab and a small-molecule inhibitor. In addition, a subset of cases in which enhanced EGFR inhibition was unproductive (6 of 16, 37.5%) exhibited marked overexpression of insulin-like growth factor 2 (IGF2). Enrichment of IGF2 overexpressors among cases with SD was demonstrated in the entire xenopatient collection and was confirmed in patients by mining clinical gene expression data sets. In functional studies, IGF2 overproduction attenuated the efficacy of cetuximab. Conversely, interception of IGF2-dependent signaling in IGF2-overexpressing xenopatients potentiated the effects of cetuximab. The clinical implementation of IGF inhibitors awaits reliable predictors of response, but the results of this study suggest rational combination therapies for colorectal cancer and provide evidence for IGF2 as a biomarker of reduced tumor sensitivity to anti-EGFR therapy and a determinant of response to combined IGF2/EGFR targeting.

Inhibition of MEK and PI3K/mTOR suppresses tumor growth but does not cause tumor regression in patient-derived xenografts of RAS-mutant colorectal carcinomas.

PURPOSE: Gene mutations along the Ras pathway (KRAS, NRAS, BRAF, PIK3CA) occur in approximately 50% of colorectal cancers (CRC) and correlate with poor response to anti-EGF receptor (EGFR) therapies. We assessed the effects of mitogen-activated protein (MAP)/extracellular signal-regulated kinase (ERK) kinase (MEK) and phosphoinositide 3-kinase (PI3K)/mTOR inhibitors, which neutralize the major Ras effectors, in patient-derived xenografts from RAS/RAF/PIK3CA-mutant metastatic CRCs (mCRC). EXPERIMENTAL DESIGN: Forty mCRC specimens harboring KRAS, NRAS, BRAF, and/or PIK3CA mutations were implanted in nonobese diabetic/severe combined immunodeficient (NOD/SCID) mice. Each xenograft was expanded into four treatment arms: placebo, the MEK inhibitor AZD6244, the PI3K/mTOR inhibitor, BEZ235, or AZD6244 + BEZ235. Cases initially treated with placebo crossed over to AZD6244, BEZ235, and the anti-EGFR monoclonal antibody cetuximab. RESULTS: At the 3-week evaluation time point, cotreatment of established tumors with AZD6244 + BEZ235 induced disease stabilization in the majority of cases (70%) but did not lead to overt tumor regression. Monotherapy was less effective, with BEZ235 displaying higher activity than AZD6244 (disease control rates, DCRs: AZD6244, 27.5%; BEZ235, 42.5%). Triple therapy with cetuximab provided further advantage (DCR, 88%). The extent of disease control declined at the 6-week evaluation time point (DCRs: AZD6244, 13.9%; BEZ235, 16.2%; AZD6244 + BEZ235, 34%). Cross-analysis of mice harboring xenografts from the same original tumor and treated with each of the different modalities revealed subgroups with preferential sensitivity to AZD6244 (12.5%), BEZ235 (35%), or AZD6244 + BEZ235 (42.5%); another subgroup (10%) showed equivalent response to any treatment. CONCLUSIONS: The prevalent growth-suppressive effects produced by MEK and PI3K/mTOR inhibition suggest that this strategy may retard disease progression in patients. However, data offer cautionary evidence against the occurrence of durable responses.