Comprehensive mutational scanning of EGFR reveals TKI sensitivities of extracellular domain mutants

The epidermal growth factor receptor, EGFR, is frequently activated in lung cancer and glioblastoma by genomic alterations including missense mutations. The different mutation spectra in these diseases are reflected in divergent responses to EGFR inhibition: significant patient benefit in lung cancer, but limited in glioblastoma. Here, we report a comprehensive mutational analysis of EGFR function. We perform saturation mutagenesis of EGFR and assess function of ~22,500 variants in a human EGFR-dependent lung cancer cell line. This approach reveals enrichment of erlotinib-insensitive variants of known and unknown significance in the dimerization, transmembrane, and kinase domains. Multiple EGFR extracellular domain variants, not associated with approved targeted therapies, are sensitive to afatinib and dacomitinib in vitro. Two glioblastoma patients with somatic EGFR G598V dimerization domain mutations show responses to dacomitinib treatment followed by within-pathway resistance mutation in one case. In summary, this comprehensive screen expands the landscape of functional EGFR variants and suggests broader clinical investigation of EGFR inhibition for cancers harboring extracellular domain mutations.

Comprehensive Mutational Analysis of the BRCA1-Associated DNA Helicase and Tumor-Suppressor FANCJ/BACH1/BRIP1

FANCJ (BRIP1/BACH1) is a hereditary breast and ovarian cancer (HBOC) gene encoding a DNA helicase. Similar to HBOC genes, BRCA1 and BRCA2, FANCJ is critical for processing DNA inter-strand crosslinks (ICL) induced by chemotherapeutics, such as cisplatin. Consequently, cells deficient in FANCJ or its catalytic activity are sensitive to ICL-inducing agents. Unfortunately, the majority of FANCJ clinical mutations remain uncharacterized, limiting therapeutic opportunities to effectively use cisplatin to treat tumors with mutated FANCJ. Here, we sought to perform a comprehensive screen to identify FANCJ loss-of-function (LOF) mutations. We developed a FANCJ lentivirus mutation library representing approximately 450 patient-derived FANCJ nonsense and missense mutations to introduce FANCJ mutants into FANCJ knockout (K/O) HeLa cells. We performed a high-throughput screen to identify FANCJ LOF mutants that, as compared with wild-type FANCJ, fail to robustly restore resistance to ICL-inducing agents, cisplatin or mitomycin C (MMC). On the basis of the failure to confer resistance to either cisplatin or MMC, we identified 26 missense and 25 nonsense LOF mutations. Nonsense mutations elucidated a relationship between location of truncation and ICL sensitivity, as the majority of nonsense mutations before amino acid 860 confer ICL sensitivity. Further validation of a subset of LOF mutations confirmed the ability of the screen to identify FANCJ mutations unable to confer ICL resistance. Finally, mapping the location of LOF mutations to a new homology model provides additional functional information. IMPLICATIONS: We identify 51 FANCJ LOF mutations, providing important classification of FANCJ mutations that will afford additional therapeutic strategies for affected patients.

CloneSifter: enrichment of rare clones from heterogeneous cell populations

BACKGROUND

Many biological processes, such as cancer metastasis, organismal development, and acquisition of resistance to cytotoxic therapy, rely on the emergence of rare sub-clones from a larger population. Understanding how the genetic and epigenetic features of diverse clones affect clonal fitness provides insight into molecular mechanisms underlying selective processes. While large-scale barcoding with NGS readout has facilitated cellular fitness assessment at the population level, this approach does not support characterization of clones prior to selection. Single-cell genomics methods provide high biological resolution, but are challenging to scale across large populations to probe rare clones and are destructive, limiting further functional analysis of important clones.

RESULTS

Here, we develop CloneSifter, a methodology for tracking and enriching rare clones throughout their response to selection. CloneSifter utilizes a CRISPR sgRNA-barcode library that facilitates the isolation of viable cells from specific clones within the barcoded population using a sequence-specific retrieval reporter. We demonstrate that CloneSifter can measure clonal fitness of cancer cell models in vitro and retrieve targeted clones at abundance as low as 1 in 1883 in a heterogeneous cell population.

CONCLUSIONS

CloneSifter provides a means to track and access specific and rare clones of interest across dynamic changes in population structure to comprehensively explore the basis of these changes.

Rhabdoid Tumors Are Sensitive to the Protein-Translation Inhibitor Homoharringtonine

PURPOSE

Rhabdoid tumors are devastating pediatric cancers in need of improved therapies. We sought to identify small molecules that exhibit in vitro and in vivo efficacy against preclinical models of rhabdoid tumor.

EXPERIMENTAL DESIGN

We screened eight rhabdoid tumor cell lines with 481 small molecules and compared their sensitivity with that of 879 other cancer cell lines. Genome-scale CRISPR-Cas9 inactivation screens in rhabdoid tumors were analyzed to confirm target vulnerabilities. Gene expression and CRISPR-Cas9 data were queried across cell lines and primary rhabdoid tumors to discover biomarkers of small-molecule sensitivity. Molecular correlates were validated by manipulating gene expression. Subcutaneous rhabdoid tumor xenografts were treated with the most effective drug to confirm in vitro results.

RESULTS

Small-molecule screening identified the protein-translation inhibitor homoharringtonine (HHT), an FDA-approved treatment for chronic myelogenous leukemia (CML), as the sole drug to which all rhabdoid tumor cell lines were selectively sensitive. Validation studies confirmed the sensitivity of rhabdoid tumor to HHT was comparable with that of CML cell lines. Low expression of the antiapoptotic gene BCL2L1, which encodes Bcl-XL, was the strongest predictor of HHT sensitivity, and HHT treatment consistently depleted Mcl-1, the synthetic-lethal antiapoptotic partner of Bcl-XL. Rhabdoid tumor cell lines and primary-tumor samples expressed low BCL2L1, and overexpression of BCL2L1 induced resistance to HHT in rhabdoid tumor cells. Furthermore, HHT treatment inhibited rhabdoid tumor cell line and patient-derived xenograft growth in vivo.

CONCLUSIONS

Rhabdoid tumor cell lines and xenografts are highly sensitive to HHT, at least partially due to their low expression of BCL2L1. HHT may have therapeutic potential against rhabdoid tumors.

DNA methyltransferase inhibition overcomes diphthamide pathway deficiencies underlying CD123-targeted treatment resistance

The interleukin-3 receptor α subunit, CD123, is expressed in many hematologic malignancies including acute myeloid leukemia (AML) and blastic plasmacytoid dendritic cell neoplasm (BPDCN). Tagraxofusp (SL-401) is a CD123-targeted therapy consisting of interleukin-3 fused to a truncated diphtheria toxin payload. Factors influencing response to tagraxofusp other than CD123 expression are largely unknown. We interrogated tagraxofusp resistance in patients and experimental models and found that it was not associated with CD123 loss. Rather, resistant AML and BPDCN cells frequently acquired deficiencies in the diphthamide synthesis pathway, impairing tagraxofusp's ability to ADP-ribosylate cellular targets. Expression of DPH1, encoding a diphthamide pathway enzyme, was reduced by DNA CpG methylation in resistant cells. Treatment with the DNA methyltransferase inhibitor azacitidine restored DPH1 expression and tagraxofusp sensitivity. We also developed a drug-dependent ADP-ribosylation assay in primary cells that correlated with tagraxofusp activity and may represent an additional novel biomarker. As predicted by these results and our observation that resistance also increased mitochondrial apoptotic priming, we found that the combination of tagraxofusp and azacitidine was effective in patient-derived xenografts treated in vivo. These data have important implications for clinical use of tagraxofusp and led to a phase 1 study combining tagraxofusp and azacitidine in myeloid malignancies.

Genotype-Fitness Maps of EGFR-Mutant Lung Adenocarcinoma Chart the Evolutionary Landscape of Resistance for Combination Therapy Optimization

Cancer evolution poses a central obstacle to cure, as resistant clones expand under therapeutic selection pressures. Genome sequencing of relapsed disease can nominate genomic alterations conferring resistance but sample collection lags behind, limiting therapeutic innovation. Genome-wide screens offer a complementary approach to chart the compendium of escape genotypes, anticipating clinical resistance. We report genome-wide open reading frame (ORF) resistance screens for first- and third-generation epidermal growth factor receptor (EGFR) inhibitors and a MEK inhibitor. Using serial sampling, dose gradients, and mathematical modeling, we generate genotype-fitness maps across therapeutic contexts and identify alterations that escape therapy. Our data expose varying dose-fitness relationship across genotypes, ranging from complete dose invariance to paradoxical dose dependency where fitness increases in higher doses. We predict fitness with combination therapy and compare these estimates to genome-wide fitness maps of drug combinations, identifying genotypes where combination therapy results in unexpected inferior effectiveness. These data are applied to nominate combination optimization strategies to forestall resistant disease.

Abstract AP14: POOLED GENOMIC SCREENS IDENTIFY ANTI-APOPTOTIC GENES AS MEDIATORS OF CHEMOTHERAPY RESISTANCE IN OVARIAN CANCER

*Co-senior authors

Primary high-grade serous ovarian cancer (HGSOC) is often sensitive to platinum and taxane combination chemotherapy, but most patients relapse with chemotherapy-resistant disease. Although alterations in DNA repair function, gene expression, apoptosis, and other pathways have been described that can mediate chemotherapy resistance in HGSOC, the full landscape of HGSOC drug resistance mechanisms and the optimal strategies to eliminate resistant disease have not been fully elucidated. We performed systematic, unbiased near-genome-scale pooled overexpression and CRISPR/Cas9 knockout screens in two BRCA2-mutant HGSOC cell lines to identify genes promoting survival following cisplatin, paclitaxel, or cisplatin/paclitaxel treatment. Anti-apoptotic genes including BCL2L1 (BCL-XL), and BCL2L2 (BCL-W) were among the top hits mediating chemotherapy resistance in the overexpression screen. In the CRISPR/Cas9 screen, loss of pro-apoptotic genes (caspases, APAF1) conferred resistance, and knockout of BCL2L1 sensitized to platinum. A secondary overexpression screen of ~400 genes confirmed anti-apoptotic proteins BCL-XL, BCL-W and BCL-2 as top resistance genes, and validated numerous other candidates. Of note, anti-apoptotic genes BCL2L1 and MCL1 are focally amplified and overexpressed in patients with primary HGSOC. In HGSOC cell lines, overexpression of BCL-XL or BCL-W, and to a lesser extent BCL-2 or MCL1, conferred platinum and taxane resistance and decreased chemotherapy-induced apoptosis in HGSOC cell lines. We systematically tested small molecule inhibitors of BCL-2, BCL-XL, MCL1, or BCL2/BCL-XL as single agents or combined with chemotherapy in HGSOC cell lines. Inhibiting BCL-XL, MCL1, or BCL2/BCL-XL, but not BCL-2, significantly increased cell death when combined with cisplatin or paclitaxel. BCL-XL, MCL1, or BCL2/BCL-XL inhibitors also synergized with olaparib, a poly- ADP-ribose inhibitor. Concomitant overexpression of BCL-XL, BCL-W, or MCL1 abrogated the sensitizing effect of the anti-apoptotic protein inhibitors, depending upon the specific inhibitor. Taken together, unbiased near-genome-scale overexpression screens and patient genomic data highlight the role of the intrinsic pathway of apoptosis in HGSOC chemotherapy resistance. Our studies validate that anti-apoptotic proteins mediate resistance to several clinically relevant drugs in HGSOC, and support that BCL-XL and MCL1 may be therapeutic targets in HGSOC, particularly in combination with DNA-damaging agents.

Citation Format: Elizabeth H. Stover, Maria B. Baco, Ofir Cohen, Yvonne Li, Elizabeth Christie, Mukta Bagul, Amy Goodale, Yenarae Lee, Sasha Pantel, Matthew Rees, Guo Wei, Adam Presser, Ioannis Zervantonakis, Patrick Bhola, Jeremy Ryan, Jennifer Guerriero, Felice Liang, Andrew Cherniack, Federica Piccioni, Ursula A. Matulonis, David D. L. Bowtell, Anthony Letai, Kris Sarosiek, Levi Garraway, Cory M. Johannessen, Matthew Meyerson. POOLED GENOMIC SCREENS IDENTIFY ANTI-APOPTOTIC GENES AS MEDIATORS OF CHEMOTHERAPY RESISTANCE IN OVARIAN CANCER [abstract]. In: Proceedings of the 12th Biennial Ovarian Cancer Research Symposium; Sep 13-15, 2018; Seattle, WA. Philadelphia (PA): AACR; Clin Cancer Res 2019;25(22 Suppl):Abstract nr AP14.

Mitochondrial Reprogramming Underlies Resistance to BCL-2 Inhibition in Lymphoid Malignancies

Mitochondrial apoptosis can be effectively targeted in lymphoid malignancies with the FDA-approved B cell lymphoma 2 (BCL-2) inhibitor venetoclax, but resistance to this agent is emerging. We show that venetoclax resistance in chronic lymphocytic leukemia is associated with complex clonal shifts. To identify determinants of resistance, we conducted parallel genome-scale screens of the BCL-2-driven OCI-Ly1 lymphoma cell line after venetoclax exposure along with integrated expression profiling and functional characterization of drug-resistant and engineered cell lines. We identified regulators of lymphoid transcription and cellular energy metabolism as drivers of venetoclax resistance in addition to the known involvement by BCL-2 family members, which were confirmed in patient samples. Our data support the implementation of combinatorial therapy with metabolic modulators to address venetoclax resistance.

Pooled Genomic Screens Identify Anti-apoptotic Genes as Targetable Mediators of Chemotherapy Resistance in Ovarian Cancer

High-grade serous ovarian cancer (HGSOC) is often sensitive to initial treatment with platinum and taxane combination chemotherapy, but most patients relapse with chemotherapy-resistant disease. To systematically identify genes modulating chemotherapy response, we performed pooled functional genomic screens in HGSOC cell lines treated with cisplatin, paclitaxel, or cisplatin plus paclitaxel. Genes in the intrinsic pathway of apoptosis were among the top candidate resistance genes in both gain-of-function and loss-of-function screens. In an open reading frame overexpression screen, followed by a mini-pool secondary screen, anti-apoptotic genes including BCL2L1 (BCL-XL) and BCL2L2 (BCL-W) were associated with chemotherapy resistance. In a CRISPR-Cas9 knockout screen, loss of BCL2L1 decreased cell survival whereas loss of proapoptotic genes promoted resistance. To dissect the role of individual anti-apoptotic proteins in HGSOC chemotherapy response, we evaluated overexpression or inhibition of BCL-2, BCL-XL, BCL-W, and MCL1 in HGSOC cell lines. Overexpression of anti-apoptotic proteins decreased apoptosis and modestly increased cell viability upon cisplatin or paclitaxel treatment. Conversely, specific inhibitors of BCL-XL, MCL1, or BCL-XL/BCL-2, but not BCL-2 alone, enhanced cell death when combined with cisplatin or paclitaxel. Anti-apoptotic protein inhibitors also sensitized HGSOC cells to the poly (ADP-ribose) polymerase inhibitor olaparib. These unbiased screens highlight anti-apoptotic proteins as mediators of chemotherapy resistance in HGSOC, and support inhibition of BCL-XL and MCL1, alone or combined with chemotherapy or targeted agents, in treatment of primary and recurrent HGSOC. IMPLICATIONS: Anti-apoptotic proteins modulate drug resistance in ovarian cancer, and inhibitors of BCL-XL or MCL1 promote cell death in combination with chemotherapy.

Neuronal differentiation and cell-cycle programs mediate response to BET-bromodomain inhibition in MYC-driven medulloblastoma

BET-bromodomain inhibition (BETi) has shown pre-clinical promise for MYC-amplified medulloblastoma. However, the mechanisms for its action, and ultimately for resistance, have not been fully defined. Here, using a combination of expression profiling, genome-scale CRISPR/Cas9-mediated loss of function and ORF/cDNA driven rescue screens, and cell-based models of spontaneous resistance, we identify bHLH/homeobox transcription factors and cell-cycle regulators as key genes mediating BETi's response and resistance. Cells that acquire drug tolerance exhibit a more neuronally differentiated cell-state and expression of lineage-specific bHLH/homeobox transcription factors. However, they do not terminally differentiate, maintain expression of CCND2, and continue to cycle through S-phase. Moreover, CDK4/CDK6 inhibition delays acquisition of resistance. Therefore, our data provide insights about the mechanisms underlying BETi effects and the appearance of resistance and support the therapeutic use of combined cell-cycle inhibitors with BETi in MYC-amplified medulloblastoma.

A Functional Landscape of Resistance to MEK1/2 and CDK4/6 Inhibition in NRAS-Mutant Melanoma

Combinatorial inhibition of MEK1/2 and CDK4/6 is currently undergoing clinical investigation in NRAS-mutant melanoma. To prospectively map the landscape of resistance to this investigational regimen, we utilized a series of gain- and loss-of-function forward genetic screens to identify modulators of resistance to clinical inhibitors of MEK1/2 and CDK4/6 alone and in combination. First, we identified NRAS-mutant melanoma cell lines that were dependent on NRAS for proliferation and sensitive to MEK1/2 and CDK4/6 combination treatment. We then used a genome-scale ORF overexpression screen and a CRISPR knockout screen to identify modulators of resistance to each inhibitor alone or in combination. These orthogonal screening approaches revealed concordant means of achieving resistance to this therapeutic modality, including tyrosine kinases, RAF, RAS, AKT, and PI3K signaling. Activated KRAS was sufficient to cause resistance to combined MEK/CDK inhibition and to replace genetic depletion of oncogenic NRAS. In summary, our comprehensive functional genetic screening approach revealed modulation of resistance to the inhibition of MEK1/2, CDK4/6, or their combination in NRAS-mutant melanoma. SIGNIFICANCE: These findings reveal that NRAS-mutant melanomas can acquire resistance to genetic ablation of NRAS or combination MEK1/2 and CDK4/6 inhibition by upregulating activity of the RTK-RAS-RAF and RTK-PI3K-AKT signaling cascade.

Mutational processes shape the landscape of TP53 mutations in human cancer

Unlike most tumor suppressor genes, the most common genetic alterations in TP53 are missense mutations^1,2. Mutant p53 protein is often abundantly expressed in cancers, and specific allelic variants exhibit dominant-negative or gain-of-function activities in experimental models^3–8. To gain a systematic view of p53 function, we interrogated loss-of-function screens conducted in hundreds of human cancer cell lines and performed TP53 saturation mutagenesis screens in an isogenic pair of TP53-wild-type and -null cell lines. We found that loss or dominant-negative inhibition of p53 function reliably enhanced cellular fitness. By integrating these data with the COSMIC mutational signatures database^9,10, we developed a statistical model that describes the TP53 mutational spectrum as a function of the baseline probability of acquiring each mutation and the fitness advantage conferred by attenuation of p53 activity. Collectively, these observations show that widely-acting and tissue-specific mutational processes combine with phenotypic selection to dictate the frequencies of recurrent TP53 mutations.

Abstract 1815: Massively parallel identification of conserved drug resistant mutations in kinases

Drug resistant mutations that arise in therapeutic targets often limit clinical responses. However, the discovery of such mutations has historically been performed one gene or mutation at a time, often over decades of experimental and clinical testing, limiting our understanding of conserved mechanisms of drug resistance.

We hypothesized that deep mutational scanning of canonical kinases may expedite this process and identify novel conserved elements that cause drug resistance when mutated (similar to the well-studied “Gatekeeper” residue). To test this, we generated cDNA-expression libraries containing all possible amino acid substitutions in CDK6, CDK4, ERK2, and EGFR. We screened each library against clinically utilized, ATP-competitive small molecule inhibitors. We then mapped the phenotypic data for over 40,000 missense mutations onto the aligned crystal structures of each protein and searched for shared structural attributes associated with drug resistance.

This analysis revealed 4 equivalent amino acid sites whose mutation conferred drug resistance to ATP-competitive inhibitors in all of our screens: the Gatekeeper residue, as well as three uncharacterized residues. One of these sites, which we have termed the “Keymaster”, was additionally found to cause resistance in published data sets of sub-saturation BRAF, HER2, BCR-ABL, and MEK1 mutagenesis screens against their respective inhibitors. We confirmed that drug resistant phenotypes are caused by these alterations utilizing growth assays and protein target phosphorylation detection assays. Mechanistically, we show preliminary evidence that Keymaster-mutant proteins are competent for drug binding, but may display elevated basal activity. Consistent with our findings, we additionally identified mutations at Keymaster residues in reported patient tumors in a number of oncogene kinases, suggesting that Keymaster mutations could be drivers of tumorigenesis, as well as drug resistance. These efforts may prove useful for characterizing somatic kinase mutations of unknown function, designing next-generation therapeutics and deepening our understanding of kinase regulation.

Citation Format: Nicole S. Persky, Desiree Hernandez, Jonathon Cordova, Amanda Walker, Lisa Brenan, Federica Piccioni, Sasha Pantel, Yenarae Lee, Amy Goodale, Xiaoping Yang, Yoichiro Mitsuishi, Mariana Do Carmo, Cong Zhu, Aleksandr Andreev, David E. Root, Cory M. Johannessen. Massively parallel identification of conserved drug resistant mutations in kinases [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 1815.

Abstract 954: Predicting synergistic drug combinations and resistance mechanisms from genomic features and single-agent response profiles

Drug combinations promise to improve clinical responses and/or forestall drug resistance. To capitalize on this promise, we need to know which drugs to combine, and which patients to give them to based on the genetic or pathological features of their disease. However, progress towards this goal has been hindered by the infeasibility of performing comprehensive drug-combination studies across thousands of cellular contexts.

We hypothesized that the basal gene-transcription state of cancer cell lines, in concert with the cell-viability profiles of single-agent small molecules, might be leveraged to nominate specific synergistic drug combinations and identify mechanisms of drug resistance, eliminating the need to test all possible drug/drug combinations across cellular models. Specifically, we predicted that inhibiting the protein product of transcripts associated with drug resistance to a given small molecule might induce drug synergy.

To test this notion, we analyzed nearly 400,000 drug-sensitivity profiles in >800 cancer cell lines to identify candidate compound-gene pairs. We identified over 100 examples where outlier expression of a single transcript was correlated with resistance to a small molecule. Of these gene/drug pairs, 9 genes represented imminently druggable targets, including established clinically-relevant relationships between the alkylating agent temozolomide and MGMT expression, and between a subset of chemotherapeutics including paclitaxel and the efflux pump ABCB1.

Inhibition of candidate “co-targets”, which included 3 previously characterized relationships and 6 novel relationships, resulted in cell-line-specific synergistic cell killing across multiple cell-line models. For validated compound-gene pairs, exogenous expression of the “co-target” was sufficient to confer resistance. For example, we found that high expression of MGLL, encoding monoglyceride lipase, was uniquely associated with lack of response to the histone lysine demethylase inhibitor GSK-J4. Endogenous or exogenous MGLL expression conferred resistance to GSK-J4, while MGLL-proficient cell lines could be sensitized to GSK-J4 up to 50-fold by co-treatment with an irreversible MGLL inhibitor.

These initial studies highlight the potential of integrating basal gene expression features with small-molecule response to nominate rational candidates for drug combinations. As public repositories of single agent response data from diverse cellular contexts continue to expand, so too will our repertoire of therapeutic combinations. Moreover, this approach permits the parallel identification of genomic features that indicate which patient populations are most likely to benefit from such combinations.

Citation Format: Matthew G. Rees, Lisa Brenan, Patrick Duggan, Cory M. Johannessen. Predicting synergistic drug combinations and resistance mechanisms from genomic features and single-agent response profiles [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 954.

A Convergence-Based Framework for Cancer Drug Resistance

Despite advances in cancer biology and therapeutics, drug resistance remains problematic. Resistance is often multifactorial, heterogeneous, and prone to undersampling. Nonetheless, many individual mechanisms of targeted therapy resistance may coalesce into a smaller number of convergences, including pathway reactivation (downstream re-engagement of original effectors), pathway bypass (recruitment of a parallel pathway converging on the same downstream output), and pathway indifference (development of a cellular state independent of the initial therapeutic target). Similar convergences may also underpin immunotherapy resistance. Such parsimonious, convergence-based frameworks may help explain resistance across tumor types and therapeutic categories and may also suggest strategies to overcome it.

Phenotypic Characterization of a Comprehensive Set of MAPK1/ERK2 Missense Mutants

Tumor-specific genomic information has the potential to guide therapeutic strategies and revolutionize patient treatment. Currently, this approach is limited by an abundance of disease-associated mutants whose biological functions and impacts on therapeutic response are uncharacterized. To begin to address this limitation, we functionally characterized nearly all (99.84%) missense mutants of MAPK1/ERK2, an essential effector of oncogenic RAS and RAF. Using this approach, we discovered rare gain- and loss-of-function ERK2 mutants found in human tumors, revealing that, in the context of this assay, mutational frequency alone cannot identify all functionally impactful mutants. Gain-of-function ERK2 mutants induced variable responses to RAF-, MEK-, and ERK-directed therapies, providing a reference for future treatment decisions. Tumor-associated mutations spatially clustered in two ERK2 effector-recruitment domains yet produced mutants with opposite phenotypes. This approach articulates an allele-characterization framework that can be scaled to meet the goals of genome-guided oncology.

Abstract A18: Predicting synergistic drug combinations from genomic features and single-agent response profiles

Drug combinations promise to improve clinical responses and/or forestall drug resistance. To capitalize on this promise, we need to know which drugs to combine, and whom to give them to based on the genetic or pathological features of their disease. However, accomplishing this goal has been precluded by the infeasibility of performing comprehensive drug-combination studies across thousands of cellular contexts. We hypothesized that the basal gene-transcription state of cancer cell lines, in concert with the response profiles of hundreds of single-agent small molecules, might be leveraged to nominate synergistic drug combinations, eliminating the need to test all possible drug/drug combinations across cellular models. Specifically, we predicted that inhibiting the protein product of transcripts associated with drug resistance to a given small molecule might induce drug synergy. To test this notion, we analyzed public cell-line drug-sensitivity data to identify candidate compound-gene pairs. We identified 7 examples in which outlier expression of a druggable candidate protein was associated with lack of single-agent response. Inhibition of 6/7 candidate co-targets resulted in cell-line-specific synergistic cell killing across multiple cell line models, validating the overall approach. For example, consistent with clinical findings, we found that high expression of the MGMT gene, encoding O-6-methylguanine-DNA methyltransferase, was uniquely associated with response to alkylating agents such as temozolomide, and that combination of the MGMT inhibitor O-6-benzylguanine with temozolomide resulted in synergistic killing. These initial studies highlight the potential of integrating basal gene expression features with small-molecule response to nominate rational candidates for drug combinations. As public repositories of single-agent response data from diverse cellular contexts continue to expand, so too will our repertoire of therapeutic combinations. Moreover, this approach permits the parallel identification of genomic features that indicate which patient populations are most likely to benefit from such combinations.

Citation Format: Matthew G. Rees, Lisa Brenan, Amanda Walker, Cory M. Johannessen. Predicting synergistic drug combinations from genomic features and single-agent response profiles [abstract]. In: Proceedings of the AACR Precision Medicine Series: Opportunities and Challenges of Exploiting Synthetic Lethality in Cancer; Jan 4-7, 2017; San Diego, CA. Philadelphia (PA): AACR; Mol Cancer Ther 2017;16(10 Suppl):Abstract nr A18.

Abstract 394: Saturation mutagenesis of KRAS reveals the functional landscape of missense variants

KRAS is the most commonly mutated oncogenes and is a major driver of tumor initiation and progression. Understanding the functional consequences of cancer-associated KRAS variants may have important clinical implications. For example, KRAS mutation status defines those that are likely to respond to EGFR-directed therapy in KRAS-mutant metastatic colorectal cancer. A compendium of all possible oncogenic KRAS alleles would serve as a roadmap for future therapeutic strategies directed at KRAS itself or downstream signaling effectors. Comprehensive mutagenesis of KRAS may also elucidate structure-function relationships that reveal novel biochemical properties that may be exploited for therapeutic gain. We performed saturation mutagenesis of both a wild-type (WT) and a G12D mutant form of KRAS cDNA and generated lentiviral expression libraries of 3,553 and 3,534 single amino acid substitution mutants of each backbone. We utilized these WT and G12D mutagenesis libraries for functional genetic screening to identify gain- and loss-of-function missense variants that alter critical oncogenic properties of KRAS. First, we sought to comprehensively identify all possible oncogenic missense mutations in KRAS that mediate oncogenic transformation. We stably transduced the WT library into immortalized human epithelial cells and evaluated growth in low attachment (GILA), an assay that is highly correlated with in vivo tumor formation. We identified all previously known hotspot oncogenic alleles of KRAS as well as many functionally relevant alleles that are also discovered at lower frequency in human tumors. Moreover, we also discovered a group of transforming KRAS variants that have not been well described in human tumors, thus revealing potentially novel activating mechanisms for oncogenic KRAS. In parallel, we utilized the G12D mutagenesis library to perform second-site suppressor screening to identify loss-of-function single amino acid changes that abrogate the transforming ability of oncogenic KRAS. We performed positive-selection screening in primary cell lines for variants that enable bypass of oncogene-induced senescence. Additionally, we conducted a negative-selection screen with the G12D library in a KRAS-dependent cancer cell line with inducible suppression of endogenous KRAS, thus identifying all possible second-site mutations that abolish KRAS-driven signaling necessary for maintenance of cellular proliferation and viability. Structure-function analysis of these data may reveal novel patterns of amino-acid changes that result in inactivation of oncogenic KRAS. In summary, this comprehensive dictionary of gain- and loss-of-function KRAS mutants will facilitate understanding of clinically important mutations and also yield novel insights into structure-function relationships that may improve our understanding of the KRAS oncogene.

Citation Format: Eejung Kim, Seav Huong Ly, Nicole S. Persky, Belinda Wang, Xiaoping Yang, Federica Piccioni, Katherine Labella, Mihir Doshi, Robert E. Lintner, Cong Zhu, Scott Steelman, David E. Root, Cory M. Johannessen, Alex B. Burgin, Laura E. MacConaill, William C. Hahn, Andrew J. Aguirre. Saturation mutagenesis of KRAS reveals the functional landscape of missense variants [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 394. doi:10.1158/1538-7445.AM2017-394

Creation of Novel Protein Variants with CRISPR/Cas9-Mediated Mutagenesis: Turning a Screening By-Product into a Discovery Tool

CRISPR/Cas9 screening has proven to be a versatile tool for genomics research. Based on unexpected results from a genome-wide screen, we developed a CRISPR/Cas9-mediated approach to mutagenesis, exploiting the allelic diversity generated by error-prone non-homologous end-joining (NHEJ) to identify novel gain-of-function and drug resistant alleles of the MAPK signaling pathway genes MEK1 and BRAF. We define the parameters of a scalable technique to easily generate cell populations containing thousands of endogenous allelic variants to map gene functions. Further, these results highlight an unexpected but important phenomenon, that Cas9-induced gain-of-function alleles are an inherent by-product of normal Cas9 loss-of-function screens and should be investigated during analysis of data from large-scale positive selection screens.

Characterizing genomic alterations in cancer by complementary functional associations

Systematic efforts to sequence the cancer genome have identified large numbers of relevant mutations and copy number alterations in human cancers; however, elucidating their functional consequences, and their interactions to drive or maintain oncogenic states, is still a significant challenge. Here we introduce REVEALER, a computational method that identifies combinations of mutually exclusive genomic alterations correlated with functional phenotypes, such as the activation or gene-dependency of oncogenic pathways or the sensitivity to a drug treatment. We use REVEALER to uncover complementary genomic alterations associated with the transcriptional activation of β-catenin and NRF2, MEK-inhibitor sensitivity, and KRAS dependency. REVEALER successfully identified both known and new associations demonstrating the power of combining functional profiles with extensive characterization of genomic alterations in cancer genomes.

Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq

To explore the distinct genotypic and phenotypic states of melanoma tumors, we applied single-cell RNA sequencing (RNA-seq) to 4645 single cells isolated from 19 patients, profiling malignant, immune, stromal, and endothelial cells. Malignant cells within the same tumor displayed transcriptional heterogeneity associated with the cell cycle, spatial context, and a drug-resistance program. In particular, all tumors harbored malignant cells from two distinct transcriptional cell states, such that tumors characterized by high levels of the MITF transcription factor also contained cells with low MITF and elevated levels of the AXL kinase. Single-cell analyses suggested distinct tumor microenvironmental patterns, including cell-to-cell interactions. Analysis of tumor-infiltrating T cells revealed exhaustion programs, their connection to T cell activation and clonal expansion, and their variability across patients. Overall, we begin to unravel the cellular ecosystem of tumors and how single-cell genomics offers insights with implications for both targeted and immune therapies.

A functional landscape of resistance to ALK inhibition in lung cancer

We conducted a large-scale functional genetic study to characterize mechanisms of resistance to ALK inhibition in ALK-dependent lung cancer cells. We identify members of known resistance pathways and additional putative resistance drivers. Among the latter were members of the P2Y purinergic receptor family of G-protein-coupled receptors (P2Y1, P2Y2, and P2Y6). P2Y receptors mediated resistance in part through a protein-kinase-C (PKC)-dependent mechanism. Moreover, PKC activation alone was sufficient to confer resistance to ALK inhibitors, whereas combined ALK and PKC inhibition restored sensitivity. We observed enrichment of gene signatures associated with several resistance drivers (including P2Y receptors) in crizotinib-resistant ALK-rearranged lung tumors compared to treatment-naive controls, supporting a role for these identified mechanisms in clinical ALK inhibitor resistance.

Abstract B33: Integration of genome-wide datasets identifies SOX10 as a lineage-specific genetic dependency in melanoma

Melanoma accounts for less than 5% of skin cancers, but a majority of skin cancer deaths. Genomic studies have firmly established the role of recurrent somatic alterations in the pathogenesis of melanoma, including mutations in BRAF, NRAS, TP53, PTEN, and CDKN2A. While the discovery of activating mutations in BRAF in over 60% of melanomas has led to the development of FDA-approved BRAF and MEK inhibitors that induce dramatic responses in BRAF V600-mutant melanomas, the inevitable resistance to BRAF/MEK inhibitors necessitates a greater understanding of melanoma biology and genetics. By integrating genome-wide functional genomic and expression data, we have identified SOX10 as a lineage-specific genetic dependency in melanoma. SOX10 dependency is highly correlated with SOX10 gene expression in melanoma. Subsequent functional in vitro experiments have confirmed the role of SOX10 in the cellular proliferation and growth of melanoma cell lines. The combination of SOX10 depletion with current FDA-approved targeted therapies for melanoma hold promise for the delay or prevention of resistance to these therapeutics and may lead to more durable control of melanoma in human patients.

Citation Format: Terence C. Wong, Cory M. Johannessen, Levi A. Garraway. Integration of genome-wide datasets identifies SOX10 as a lineage-specific genetic dependency in melanoma. [abstract]. In: Proceedings of the AACR Precision Medicine Series: Drug Sensitivity and Resistance: Improving Cancer Therapy; Jun 18-21, 2014; Orlando, FL. Philadelphia (PA): AACR; Clin Cancer Res 2015;21(4 Suppl): Abstract nr B33.

Abstract PR10: A cyclic AMP-regulated melanocyte lineage program confers resistance to MAP kinase pathway inhibition

BRAFV600E-mutant malignant melanomas depend on RAF/MEK/ERK (MAP kinase) signaling for tumor cell growth. RAF and MEK inhibitors show remarkable clinical efficacy in BRAFV600E melanoma; however, resistance to these agents remains a formidable challenge. Systematic characterization of resistance mechanisms may inform the development of more effective therapeutic combinations. To chart the landscape of potential resistance mechanisms, we performed systematic gain-of-function resistance studies by expressing >15,500 genes individually in a BRAFV600E melanoma cell line treated with RAF, MEK, ERK, or combined RAF/MEK inhibitors. These studies revealed a cyclic AMP-dependent melanocytic signaling network not previously associated with drug resistance, including G-protein coupled receptors, adenyl cyclase, protein kinase A and cAMP response element binding protein (CREB). In biopsies from BRAFV600E melanoma patients, phosphorylated (active) CREB was suppressed by RAF/MEK-inhibition but restored in relapsing tumors. Expression of transcription factors activated downstream of MAP kinase and cAMP pathways also conferred resistance, including c-FOS, NR4A1, NR4A2 and MITF. While direct small-molecule inhibition of transcription factors is a current challenge, MITF expression can be suppressed with histone deacetylase (HDAC) inhibitors. Combined treatment with MAP kinase pathway and a panel of HDAC inhibitors suppressed MITF expression and cAMP-mediated resistance. Collectively, these data suggest that oncogenic dysregulation of a melanocyte lineage dependency can cause resistance to RAF/MEK/ERK inhibition, which may be overcome by combining signaling- and chromatin-directed therapeutics.

This abstract is also presented as Poster B23.

Citation Format: Cory M. Johannessen, Laura A. Johnson, Federica Piccioni, Aisha Townes, Dennie T. Frederick, Keith T. Flaherty, Jennifer A. Wargo, David E. Root, Levi A. Garraway. A cyclic AMP-regulated melanocyte lineage program confers resistance to MAP kinase pathway inhibition. [abstract]. In: Proceedings of the AACR Precision Medicine Series: Synthetic Lethal Approaches to Cancer Vulnerabilities; May 17-20, 2013; Bellevue, WA. Philadelphia (PA): AACR; Mol Cancer Ther 2013;12(5 Suppl):Abstract nr PR10.

A melanoma cell state distinction influences sensitivity to MAPK pathway inhibitors

Most melanomas harbor oncogenic BRAF(V600) mutations, which constitutively activate the MAPK pathway. Although MAPK pathway inhibitors show clinical benefit in BRAF(V600)-mutant melanoma, it remains incompletely understood why 10% to 20% of patients fail to respond. Here, we show that RAF inhibitor-sensitive and inhibitor-resistant BRAF(V600)-mutant melanomas display distinct transcriptional profiles. Whereas most drug-sensitive cell lines and patient biopsies showed high expression and activity of the melanocytic lineage transcription factor MITF, intrinsically resistant cell lines and biopsies displayed low MITF expression but higher levels of NF-κB signaling and the receptor tyrosine kinase AXL. In vitro, these MITF-low/NF-κB-high melanomas were resistant to inhibition of RAF and MEK, singly or in combination, and ERK. Moreover, in cell lines, NF-κB activation antagonized MITF expression and induced both resistance marker genes and drug resistance. Thus, distinct cell states characterized by MITF or NF-κB activity may influence intrinsic resistance to MAPK pathway inhibitors in BRAF(V600)-mutant melanoma.

SIGNIFICANCE

Although most BRAF(V600)-mutant melanomas are sensitive to RAF and/or MEK inhibitors, a subset fails to respond to such treatment. This study characterizes a transcriptional cell state distinction linked to MITF and NF-κB that may modulate intrinsic sensitivity of melanomas to MAPK pathway inhibitors.

MAP kinase pathway alterations in BRAF-mutant melanoma patients with acquired resistance to combined RAF/MEK inhibition

Treatment of BRAF-mutant melanoma with combined dabrafenib and trametinib, which target RAF and the downstream MAP-ERK kinase (MEK)1 and MEK2 kinases, respectively, improves progression-free survival and response rates compared with dabrafenib monotherapy. Mechanisms of clinical resistance to combined RAF/MEK inhibition are unknown. We performed whole-exome sequencing (WES) and whole-transcriptome sequencing (RNA-seq) on pretreatment and drug-resistant tumors from five patients with acquired resistance to dabrafenib/trametinib. In three of these patients, we identified additional mitogen-activated protein kinase (MAPK) pathway alterations in the resistant tumor that were not detected in the pretreatment tumor, including a novel activating mutation in MEK2 (MEK2(Q60P)). MEK2(Q60P) conferred resistance to combined RAF/MEK inhibition in vitro, but remained sensitive to inhibition of the downstream kinase extracellular signal-regulated kinase (ERK). The continued MAPK signaling-based resistance identified in these patients suggests that alternative dosing of current agents, more potent RAF/MEK inhibitors, and/or inhibition of the downstream kinase ERK may be needed for durable control of BRAF-mutant melanoma.

The genetic landscape of clinical resistance to RAF inhibition in metastatic melanoma

Most patients with BRAF(V600)-mutant metastatic melanoma develop resistance to selective RAF kinase inhibitors. The spectrum of clinical genetic resistance mechanisms to RAF inhibitors and options for salvage therapy are incompletely understood. We performed whole-exome sequencing on formalin-fixed, paraffin-embedded tumors from 45 patients with BRAF(V600)-mutant metastatic melanoma who received vemurafenib or dabrafenib monotherapy. Genetic alterations in known or putative RAF inhibitor resistance genes were observed in 23 of 45 patients (51%). Besides previously characterized alterations, we discovered a "long tail" of new mitogen-activated protein kinase (MAPK) pathway alterations (MAP2K2, MITF) that confer RAF inhibitor resistance. In three cases, multiple resistance gene alterations were observed within the same tumor biopsy. Overall, RAF inhibitor therapy leads to diverse clinical genetic resistance mechanisms, mostly involving MAPK pathway reactivation. Novel therapeutic combinations may be needed to achieve durable clinical control of BRAF(V600)-mutant melanoma. Integrating clinical genomics with preclinical screens may model subsequent resistance studies.

Whole exome and whole transcriptome sequencing in melanoma patients to identify mechanisms of resistance to combined RAF/MEK inhibition

9015

Background: The RAF inhibitors vemurafenib and dabrafenib (D) and the MEK inhibitor trametinib (T) improve survival as monotherapies in BRAF-mutant melanoma. Since clinical mechanisms of resistance (MoR) result in MAPK pathway reactivation, recent efforts have focused on combined targeting of RAF and MEK. The combination of D and T (D/T) increased progression-free survival and response rate compared with D alone (Flaherty et al, NEJM, 2012). The MoR to this combination remain unknown. Methods: To look for clinical MoR to combined RAF/MEK inhibition, we performed whole exome (WES) and whole transcriptome sequencing (RNASeq) on tumors from 4 patients (pts) with acquired resistance and 1 pt with intrinsic resistance to D/T. Pre-treatment and post-resistance tumors from all pts were analyzed for point mutations, insertions/deletions, copy number alterations, alternatively spliced transcripts, rearrangements, and expression changes. Results: In 2 of 4 pts with acquired resistance, WES identified mutations in MEK1 and MEK2 that were undetectable in the pre-treatment tumors. In the 3rd pt, RNASeq identified an alternatively spliced isoform of BRAF lacking exons 2-10, also undetectable in the pre-treatment tumor. In the 4th pt, no obvious MoR were seen, though multiple alterations were enriched in the post-resistance tumor. The pt with intrinsic resistance had several alterations in genes that conferred resistance to RAF/MEK inhibition when overexpressed in BRAF-mutant cell lines. Integration of WES and RNASeq data also identified several co-existing alterations that may synergize to increase resistance. Conclusions: Analysis of combined WES and RNASeq data from pt samples provides a more complete picture of clinical MoR to MAPK-targeted therapy. Post-resistance tumors from 3 of 4 pts with acquired resistance to D/T had alterations in MAPK genes not detectable in the pre-treatment tumors, suggesting that resistance involves reactivation of the MAPK pathway despite combined RAF/MEK inhibition. Alternative dosing of current agents, more potent RAF/MEK inhibitors, and/or inhibition of the downstream kinase ERK may be needed for durable control of BRAF-mutant melanoma.

MicroSCALE screening reveals genetic modifiers of therapeutic response in melanoma

Cell microarrays are a promising tool for performing large-scale functional genomic screening in mammalian cells at reasonable cost, but owing to technical limitations they have been restricted for use with a narrow range of cell lines and short-term assays. Here, we describe MicroSCALE (Microarrays of Spatially Confined Adhesive Lentiviral Features), a cell microarray-based platform that enables application of this technology to a wide range of cell types and longer-term assays. We used MicroSCALE to uncover kinases that when overexpressed partially desensitized B-RAFV600E-mutant melanoma cells to inhibitors of the mitogen-activated protein kinase kinase kinase (MAPKKK) RAF, the MAPKKs MEK1 and 2 (MEK1/2, mitogen-activated protein kinase kinase 1 and 2), mTOR (mammalian target of rapamycin), or PI3K (phosphatidylinositol 3-kinase). These screens indicated that cells treated with inhibitors acting through common mechanisms were affected by a similar profile of overexpressed proteins. In contrast, screens involving inhibitors acting through distinct mechanisms yielded unique profiles, a finding that has potential relevance for small-molecule target identification and combination drugging studies. Further, by integrating large-scale functional screening results with cancer cell line gene expression and pharmacological sensitivity data, we validated the nuclear factor κB pathway as a potential mediator of resistance to MAPK pathway inhibitors. The MicroSCALE platform described here may enable new classes of large-scale, resource-efficient screens that were not previously feasible, including those involving combinations of cell lines, perturbations, and assay outputs or those involving limited numbers of cells and limited or expensive reagents.

Abstract PR5: Identifying mechanisms of drug resistance to MAPK pathway inhibition via genome-scale rescue screens

The recent clinical success of RAF inhibitors in BRAF(V600E)-mutant malignant melanoma highlights the promise of targeted therapeutics as well as the inevitability of acquired drug resistance: responses are sustained for less than a year and all patients ultimately become drug resistant. To identify causative drug resistance genes in BRAF(V600E)-mutant melanoma, we have individually expressed over 14, 000 unique genes and queried their ability to affect sensitivity to a panel of mitogen activated protein kinase (MAPK) pathway inhibitors, including a RAF, MEK and ERK inhibitor, both as single agents or in combination (RAFi+MEKi). Out of this genome-scale collection, less than 1% of genes are sufficient to induce resistance to the MAPK pathway inhibitors tested. Here, we present the results of this rescue screen and extensive validation/characterization efforts that reveal a number of novel insights into mechanisms of drug resistance, melanoma biology and the MAPK signaling pathway. Additionally, the scale of this first-of-kind functional rescue screen permits its integration with existing genome characterization efforts, including melanoma sequencing studies, gene expression data and parallel functional assays to identify the most relevant genetic modifiers of resistance to MAPK pathway inhibition.

This proffered talk is also presented as Poster B27.

Dissecting therapeutic resistance to RAF inhibition in melanoma by tumor genomic profiling

A detailed understanding of the mechanisms by which tumors acquire resistance to targeted anticancer agents should speed the development of treatment strategies with lasting clinical efficacy. RAF inhibition in BRAF-mutant melanoma exemplifies the promise and challenge of many targeted drugs; although response rates are high, resistance invariably develops. Here, we articulate overarching principles of resistance to kinase inhibitors, as well as a translational approach to characterize resistance in the clinical setting through tumor mutation profiling. As a proof of principle, we performed targeted, massively parallel sequencing of 138 cancer genes in a tumor obtained from a patient with melanoma who developed resistance to PLX4032 after an initial dramatic response. The resulting profile identified an activating mutation at codon 121 in the downstream kinase MEK1 that was absent in the corresponding pretreatment tumor. The MEK1(C121S) mutation was shown to increase kinase activity and confer robust resistance to both RAF and MEK inhibition in vitro. Thus, MEK1(C121S) or functionally similar mutations are predicted to confer resistance to combined MEK/RAF inhibition. These results provide an instructive framework for assessing mechanisms of acquired resistance to kinase inhibition and illustrate the use of emerging technologies in a manner that may accelerate personalized cancer medicine.

A public genome-scale lentiviral expression library of human ORFs

Functional characterization of the human genome requires tools for systematically modulating gene expression in both loss-of-function and gain-of-function experiments. We describe the production of a sequence-confirmed, clonal collection of over 16,100 human open-reading frames (ORFs) encoded in a versatile Gateway vector system. Using this ORFeome resource, we created a genome-scale expression collection in a lentiviral vector, thereby enabling both targeted experiments and high-throughput screens in diverse cell types.

Abstract 3933: Dissecting therapeutic resistance to RAF inhibition in melanoma by tumor genomic profiling

Although targeted anticancer agents often result in dramatic responses, tumors invariably become resistant to these agents. A detailed understanding of the mechanisms by which tumors acquire resistance to such inhibitors should speed the development of more durable treatment strategies. The promise of RAF inhibition in metastatic melanoma is exemplary in this regard. Mutations in the BRAF oncogene have been detected in more than 50% of metastatic melanomas, and inhibition of mutant BRAF has shown high response rates in early clinical trials of patients with melanoma. However, resistance to RAF inhibition invariably develops. Mechanisms of de novo and acquired resistance to RAF inhibition remain poorly understood. Moreover, the clinical application of genomic approaches that might be capable of diagnosing salient resistance mechanisms remains underdeveloped.

Here, we describe an approach to characterize genetic mechanisms of resistance through systematic tumor mutation profiling. We performed massively parallel sequencing of 138 cancer genes in a tumor specimen from a melanoma patient who developed resistance to PLX4032 after a dramatic initial response. The resulting profile identified a novel mutation in the downstream kinase MEK1 that was absent in the corresponding pre-treatment tumor. This MEK1 mutation was shown to increase kinase activity and confer robust resistance to both RAF and MEK inhibition in vitro. Our results provide an instructive framework both for assessing mechanisms of acquired resistance to kinase inhibition and deploying new technologies for elaboration of resistance mechanisms in a manner that may accelerate personalized cancer medicine.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 3933. doi:10.1158/1538-7445.AM2011-3933

COT drives resistance to RAF inhibition through MAP kinase pathway reactivation

Oncogenic mutations in the serine/threonine kinase B-RAF (also known as BRAF) are found in 50-70% of malignant melanomas. Pre-clinical studies have demonstrated that the B-RAF(V600E) mutation predicts a dependency on the mitogen-activated protein kinase (MAPK) signalling cascade in melanoma-an observation that has been validated by the success of RAF and MEK inhibitors in clinical trials. However, clinical responses to targeted anticancer therapeutics are frequently confounded by de novo or acquired resistance. Identification of resistance mechanisms in a manner that elucidates alternative 'druggable' targets may inform effective long-term treatment strategies. Here we expressed ∼600 kinase and kinase-related open reading frames (ORFs) in parallel to interrogate resistance to a selective RAF kinase inhibitor. We identified MAP3K8 (the gene encoding COT/Tpl2) as a MAPK pathway agonist that drives resistance to RAF inhibition in B-RAF(V600E) cell lines. COT activates ERK primarily through MEK-dependent mechanisms that do not require RAF signalling. Moreover, COT expression is associated with de novo resistance in B-RAF(V600E) cultured cell lines and acquired resistance in melanoma cells and tissue obtained from relapsing patients following treatment with MEK or RAF inhibitors. We further identify combinatorial MAPK pathway inhibition or targeting of COT kinase activity as possible therapeutic strategies for reducing MAPK pathway activation in this setting. Together, these results provide new insights into resistance mechanisms involving the MAPK pathway and articulate an integrative approach through which high-throughput functional screens may inform the development of novel therapeutic strategies.

Abstract 5769: A novel, high-throughput approach to discovering drug-resistance genes identifies CRAF, CrkL, and TPL2/COT as drivers of resistance to BRAFV600E inhibition

Activating mutations in the BRAF kinase are found in 50-70% of melanomas and have nominated BRAF, along with downstream components of the MAP kinase pathway, as therapeutic targets. Until recently, MAPK pathway inhibitors have shown poor efficacy in clinical trials, leading many to doubt the therapeutic potential of BRAF/MAPK inhibition. However, the recent clinical success of the BRAFV600E-specific kinase inhibitor PLX4032 has validated MAPK pathway inhibition as a viable treatment strategy. Despite initial success, most targeted therapeutics will eventually succumb to acquired resistance. It is therefore imperative that we identify drivers of targeted drug resistance in a rapid, systematic manner that identifies parallel drug targets, ultimately facilitating an effective long-term treatment strategy. Accordingly, we report here a novel, high-throughput approach to systematically identify kinases capable of driving resistance to clinically efficacious, BRAFV600E-specific kinase inhibitors. We screened an open reading frame (ORF) expression library consisting of ∼600 kinases and kinase-related ORFs for their ability to drive resistance to PLX4720 in BRAFV600E-mutant cell lines. Parallel analytical approaches and subsequent validation identified RAF1 (CRAF), CRKL (CrkL) and MAP3K8 (TPL2/COT) as the strongest candidate resistance genes. CRAF and TPL2/COT are both MAP3 kinases that directly phosphorylate and activate MEK, whereas CrkL is an adapter protein that is amplified and phosphorylated in melanoma. CRAF, CrkL and TPL2/COT expressing cells proliferate and show constitutive MAPK pathway activation at concentrations of PLX4720 that are 10-fold over IG50. While CRAF-expressing cells show a dose-dependent increase in MAPK activation when treated with PLX4720, MAPK pathway activation is weakly dampened in CrkL and TPL2/COT-expressing cells upon PLX4720 treatment, suggesting that CrkL and TPL2/COT drive resistance largely, but not exclusively, through MAPK pathway activation. Accordingly, CRAF, CrkL and TPL2/COT segregate into two categories when queried for their ability to drive resistance to a panel of MAPK pathway inhibitors: CRAF mediates resistance only to BRAFV600E-specific inhibitors, while CrkL and TPL2/COT drive resistance to BRAFV600E-specific inhibitors and modify sensitivity to allosteric inhibitors of MEK kinases. Combinatorial treatment with BRAFV600E inhibitors and MEK inhibitors effectively inhibits the proliferation of cells expressing CRAF, CrkL and TPL2/COT, suggesting that combinatorial strategies may address acquired resistance in the clinic. These data highlight the utility and broad applicability of rapid, ORF-based methods of uncovering drivers of targeted resistance and identify CRAF, CrkL and TPL2/COT as putative resistance alleles to BRAFV600E-specific inhibition.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 5769.

TORC1 is essential for NF1-associated malignancies

Inactivating mutations in NF1 underlie the prevalent familial cancer syndrome neurofibromatosis type 1 [1]. The NF1-encoded protein is a Ras GTPase-activating protein (RasGAP) [2]. Accordingly, Ras is aberrantly activated in NF1-deficient tumors; however, it is unknown which effector pathways critically function in tumor development. Here we provide in vivo evidence that TORC1/mTOR activity is essential for tumorigenesis. Specifically, we show that the mTOR inhibitor rapamycin potently suppresses the growth of aggressive NF1-associated malignancies in a genetically engineered murine model. However, in these tumors rapamycin does not function via mechanisms generally assumed to mediate tumor suppression, including inhibition of HIF-1alpha and indirect suppression of AKT, but does suppress the mTOR target Cyclin D1 [3]. These results demonstrate that mTOR inhibitors may be an effective targeted therapy for this commonly untreatable malignancy. Moreover, they indicate that mTOR inhibitors do not suppress all tumor types via the same mechanism, suggesting that current biomarkers that rely on HIF-1alpha suppression may not be informative for all cancers. Finally, our results reveal important differences between the effects of mTOR inhibition on the microvasculature in genetically engineered versus xenograft models and indicate that the former may be required for effective preclinical screening with this class of inhibitors.

Tumor-specific p73 up-regulation mediates p63 dependence in squamous cell carcinoma

p63 is essential for normal epithelial development and is overexpressed in the vast majority of squamous cell carcinomas (SCC). Recent work had shown that DeltaNp63alpha is essential for survival of SCC cells, raising the possibility that the p63 pathway may be an attractive therapeutic target in these tumors. Nevertheless, it is unknown whether a therapeutic window exists for inhibiting p63 in tumor cells versus normal epithelia. Here, we show that SCC cells are uniquely dependent on DeltaNp63alpha for survival, unlike normal p63-expressing epithelial cells, and that dependence is mediated through tumor-specific up-regulation of the related protein p73. In normal primary human keratinocytes, we find that inhibition of endogenous p63 by RNA interference (RNAi) induces p21(CIP1) expression, inhibits cell cycle progression, and ultimately promotes cellular senescence. In contrast, p63 inhibition in SCC cells induces proapoptotic bcl-2 family members and rapidly triggers apoptosis. Expression of p73 is low in uncultured basal keratinocytes but is markedly up-regulated in both SCC cell lines and primary tumors in vivo. Whereas p21(CIP1) induction following loss of p63 in normal cells is independent of p53 and p73, both proapoptotic gene induction and cell death following p63 RNAi in tumor cells are p73 dependent. Finally, ectopic p73 expression in primary keratinocytes does not affect baseline cell proliferation but is sufficient to trigger cell death following loss of p63. Together, these findings define a specific molecular mechanism of p63 dependence through p73 up-regulation, and they provide a rationale for targeting the p63 pathway as a therapeutic strategy in SCCs.

A negative feedback signaling network underlies oncogene-induced senescence

Oncogene-induced senescence functions to limit tumor development. However, a complete understanding of the signals that trigger this type of senescence is currently lacking. We found that mutations affecting NF1, Raf, and Ras induce a global negative feedback response that potently suppresses Ras and/or its effectors. Moreover, these signals promote senescence by inhibiting the Ras/PI3K pathway, which can impact the senescence machinery through HDM2 and FOXO. This negative feedback program is regulated in part by RasGEFs, Sprouty proteins, RasGAPs, and MKPs. Moreover, these signals function in vivo in benign human tumors. Thus, the ultimate response to the aberrant activation of the Ras pathway is a multifaceted negative feedback signaling network that terminates the oncogenic signal and participates in the senescence response.

Regulation of mTOR and cell growth in response to energy stress by REDD1

The tuberous sclerosis tumor suppressors TSC1 and TSC2 regulate the mTOR pathway to control translation and cell growth in response to nutrient and growth factor stimuli. We have recently identified the stress response REDD1 gene as a mediator of tuberous sclerosis complex (TSC)-dependent mTOR regulation by hypoxia. Here, we demonstrate that REDD1 inhibits mTOR function to control cell growth in response to energy stress. Endogenous REDD1 is induced following energy stress, and REDD1-/- cells are highly defective in dephosphorylation of the key mTOR substrates S6K and 4E-BP1 following either ATP depletion or direct activation of the AMP-activated protein kinase (AMPK). REDD1 likely acts on the TSC1/2 complex, as regulation of mTOR substrate phosphorylation by REDD1 requires TSC2 and is blocked by overexpression of the TSC1/2 downstream target Rheb but is not blocked by inhibition of AMPK. Tetracycline-inducible expression of REDD1 triggers rapid dephosphorylation of S6K and 4E-BP1 and significantly decreases cellular size. Conversely, inhibition of endogenous REDD1 by short interfering RNA increases cell size in a rapamycin-sensitive manner, and REDD1-/- cells are defective in cell growth regulation following ATP depletion. These results define REDD1 as a critical transducer of the cellular response to energy depletion through the TSC-mTOR pathway.

The NF1 tumor suppressor critically regulates TSC2 and mTOR

Loss-of-function mutations in the NF1 tumor suppressor gene underlie the familial cancer syndrome neurofibromatosis type I (NF1). The NF1-encoded protein, neurofibromin, functions as a Ras-GTPase activating protein (RasGAP). Accordingly, deregulation of Ras is thought to contribute to NF1 development. However, the critical effector pathways involved in disease pathogenesis are still unknown. We show here that the mTOR pathway is tightly regulated by neurofibromin. mTOR is constitutively activated in both NF1-deficient primary cells and human tumors in the absence of growth factors. This aberrant activation depends on Ras and PI3 kinase, and is mediated by the phosphorylation and inactivation of the TSC2-encoded protein tuberin by AKT. Importantly, tumor cell lines derived from NF1 patients, and a genetically engineered cell system that requires Nf1-deficiency for transformation, are highly sensitive to the mTOR inhibitor rapamycin. Furthermore, while we show that the activation of endogenous Ras leads to constitutive mTOR signaling in this disease state, we also demonstrate that in normal cells Ras is differentially required for mTOR signaling in response to various growth factors. Thus, these findings identify the NF1 tumor suppressor as an indispensable regulator of TSC2 and mTOR. Furthermore, our results also demonstrate that Ras plays a critical role in the activation of mTOR in both normal and tumorigenic settings. Finally, these data suggest that rapamycin, or its derivatives, may represent a viable therapy for NF1.

REDD1, a developmentally regulated transcriptional target of p63 and p53, links p63 to regulation of reactive oxygen species

We identified REDD1 as a novel transcriptional target of p53 induced following DNA damage. During embryogenesis, REDD1 expression mirrors the tissue-specific pattern of the p53 family member p63, and TP63 null embryos show virtually no expression of REDD1, which is restored in mouse embryo fibroblasts following p63 expression. In differentiating primary keratinocytes, TP63 and REDD1 expression are coordinately downregulated, and ectopic expression of either gene inhibits in vitro differentiation. REDD1 appears to function in the regulation of reactive oxygen species (ROS); we show that TP63 null fibroblasts have decreased ROS levels and reduced sensitivity to oxidative stress, which are both increased following ectopic expression of either TP63 or REDD1. Thus, REDD1 encodes a shared transcriptional target that implicates ROS in the p53-dependent DNA damage response and in p63-mediated regulation of epithelial differentiation.