Skip to main content
  • AACR Publications
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

AACR logo

  • Register
  • Log in
  • My Cart
Advertisement

Main menu

  • Home
  • About
    • The Journal
    • AACR Journals
    • Subscriptions
    • Permissions and Reprints
    • Reviewing
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • Meeting Abstracts
    • Collections
      • COVID-19 & Cancer Resource Center
      • Focus on Computer Resources
      • Highly Cited Collection
      • Editors' Picks
      • "Best of" Collection
  • For Authors
    • Information for Authors
    • Author Services
    • Early Career Award
    • Best of: Author Profiles
    • Submit
  • Alerts
    • Table of Contents
    • Editors' Picks
    • OnlineFirst
    • Citations
    • Author/Keyword
    • RSS Feeds
    • My Alert Summary & Preferences
  • News
    • Cancer Discovery News
  • COVID-19
  • Webinars
  • Search More

    Advanced Search

  • AACR Publications
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

User menu

  • Register
  • Log in
  • My Cart

Search

  • Advanced search
Cancer Research
Cancer Research
  • Home
  • About
    • The Journal
    • AACR Journals
    • Subscriptions
    • Permissions and Reprints
    • Reviewing
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • Meeting Abstracts
    • Collections
      • COVID-19 & Cancer Resource Center
      • Focus on Computer Resources
      • Highly Cited Collection
      • Editors' Picks
      • "Best of" Collection
  • For Authors
    • Information for Authors
    • Author Services
    • Early Career Award
    • Best of: Author Profiles
    • Submit
  • Alerts
    • Table of Contents
    • Editors' Picks
    • OnlineFirst
    • Citations
    • Author/Keyword
    • RSS Feeds
    • My Alert Summary & Preferences
  • News
    • Cancer Discovery News
  • COVID-19
  • Webinars
  • Search More

    Advanced Search

Translational Science

Src-Dependent DBL Family Members Drive Resistance to Vemurafenib in Human Melanoma

Charlotte R. Feddersen, Jacob L. Schillo, Afshin Varzavand, Hayley R. Vaughn, Lexy S. Wadsworth, Andrew P. Voigt, Eliot Y. Zhu, Brooke M. Jennings, Sarah A. Mullen, Jeremy Bobera, Jesse D. Riordan, Christopher S. Stipp and Adam J. Dupuy
Charlotte R. Feddersen
1Department of Anatomy and Cell Biology, Carver College of Medicine, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jacob L. Schillo
1Department of Anatomy and Cell Biology, Carver College of Medicine, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Afshin Varzavand
2Department of Biology, College of Liberal Arts and Sciences, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hayley R. Vaughn
1Department of Anatomy and Cell Biology, Carver College of Medicine, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lexy S. Wadsworth
1Department of Anatomy and Cell Biology, Carver College of Medicine, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrew P. Voigt
1Department of Anatomy and Cell Biology, Carver College of Medicine, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eliot Y. Zhu
1Department of Anatomy and Cell Biology, Carver College of Medicine, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brooke M. Jennings
2Department of Biology, College of Liberal Arts and Sciences, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sarah A. Mullen
2Department of Biology, College of Liberal Arts and Sciences, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jeremy Bobera
2Department of Biology, College of Liberal Arts and Sciences, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jesse D. Riordan
1Department of Anatomy and Cell Biology, Carver College of Medicine, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Christopher S. Stipp
2Department of Biology, College of Liberal Arts and Sciences, University of Iowa, Iowa City, Iowa.
3Holden Comprehensive Cancer Center, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: adam-dupuy@uiowa.edu christopher-stipp@uiowa.edu
Adam J. Dupuy
1Department of Anatomy and Cell Biology, Carver College of Medicine, University of Iowa, Iowa City, Iowa.
3Holden Comprehensive Cancer Center, University of Iowa, Iowa City, Iowa.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: adam-dupuy@uiowa.edu christopher-stipp@uiowa.edu
DOI: 10.1158/0008-5472.CAN-19-0244 Published October 2019
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

The use of selective BRAF inhibitors (BRAFi) has produced remarkable outcomes for patients with advanced cutaneous melanoma harboring a BRAFV600E mutation. Unfortunately, the majority of patients eventually develop drug-resistant disease. We employed a genetic screening approach to identify gain-of-function mechanisms of BRAFi resistance in two independent melanoma cell lines. Our screens identified both known and unappreciated drivers of BRAFi resistance, including multiple members of the DBL family. Mechanistic studies identified a DBL/RAC1/PAK signaling axis capable of driving resistance to both current and next-generation BRAFis. However, we show that the SRC inhibitor, saracatinib, can block the DBL-driven resistance. Our work highlights the utility of our straightforward genetic screening method in identifying new drug combinations to combat acquired BRAFi resistance.

Significance: A simple, rapid, and flexible genetic screening approach identifies genes that drive resistance to MAPK inhibitors when overexpressed in human melanoma cells.

Introduction

Melanoma is the deadliest form of skin cancer, with around 90,000 diagnoses of invasive disease and approximately 10,000 deaths per year (1). Patients had few treatment options until the development of vemurafenib, a highly selective kinase inhibitor that specifically targets the BRAFV600E-mutant protein present in approximately 50% of all melanoma cases (2). Initially, vemurafenib provided complete or partial response in over 50% of patients and increased progression-free survival (3). Unfortunately, most patients relapse once tumors acquire resistance to vemurafenib.

Genetic analysis of progression samples has identified resistance mechanisms, including amplification of BRAFV600E, expression of truncated BRAFV600E, and RAS mutation (4–6). However, these mechanisms explain only approximately 60% of cases of BRAF inhibitor (BRAFi) resistance (5, 7, 8). Drug resistance can be delayed by combining vemurafenib with cobimetinib, an MEK inhibitor (MEKi), but most patients eventually develop progressive disease via resistance mechanisms that have not been well characterized (7). Thus, unexplained cases of resistance to MAPK inhibition (MAPKi) in human melanoma represent an important unmet clinical need.

Mechanisms of vemurafenib resistance have been studied in BRAFV600E-mutant human melanoma cell lines using genome-wide shRNA and CRISPR loss-of-function screens (9–11). Overall, these studies showed little overlap in candidate mechanisms. Two screens have been reported that attempted to identify drivers of vemurafenib resistance by high-throughput overexpression of genes via lentiviral libraries (12, 13). Importantly, these screens failed to identify known mechanisms of vemurafenib resistance (e.g., BRAFV600E amplification or N-terminal truncation; refs. 4, 6). These observations led us to develop a simple insertional mutagenesis screening approach using the Sleeping Beauty (SB) transposon system to identify novel drivers of vemurafenib resistance in an unbiased forward genetic screen.

The SB system is a well-established tool for developing mouse models of spontaneous cancer in which transposon-induced somatic mutations drive transformation (14). In this context, the SB system consists of two parts: a mutagenic transposon vector and the transposase enzyme. When introduced into the same cell, the transposase excises the transposon from a donor vector and integrates it at a TA dinucleotide site in the host cell genome. In the context of some selective pressure (e.g., proliferation, drug treatment), cells with specific mutations conferring the ability to outcompete neighboring cells undergo clonal expansion, facilitating subsequent identification of these phenotype-driving mutations. This method has been used to select for specific phenotypes using ex vivo cell-based assays (15–17), as well as to drive the development of drug-resistant tumors in engineered mouse models (18, 19). However, previous ex vivo approaches using human cells have been limited by the relative inefficiency of delivering both transposon and transposase vectors to cells. Moreover, prior studies have required the isolation of clonal cell populations to identify insertional mutations associated with the desired phenotype, a process that greatly reduces screen throughput. Collectively, these challenges have limited the broader application of SB mutagenesis in ex vivo screening approaches.

We present here the results of three SB mutagenesis drug resistance screens conducted in human BRAFV600E-mutant melanoma cells to identify novel drivers of resistance to either BRAFi treatment alone or BRAFi/MEKi combination treatment. Importantly, we detected recurrent N-terminal truncations of BRAF as a driver of BRAFi resistance, a mechanism previously associated with BRAFi resistance in human melanoma (6). We also identified MCF2, VAV1, PDGFRB, and N-terminally truncated RAF1 as drivers of BRAFi resistance. We experimentally verified the ability of candidates to drive drug resistance in independent melanoma cell lines, and analysis of transcriptome data from clinical progression samples revealed that overexpression of our candidates is associated with BRAFi resistance in human patients. Finally, we elucidate a mechanism through which the DBL family members MCF2 and VAV1 act to drive drug resistance and show that this mechanism can be blocked with saracatinib, an inhibitor of the Src family. These findings demonstrate the utility of our genetic screening approach to identify clinically relevant drivers of drug resistance, as well as the potential for discovering new therapeutic approaches to reverse or prevent its occurrence.

Materials and Methods

SB mutagenesis screen

A375 cells were stably transfected via Effectene (Qiagen) coupled with a piggyBac transposase integration system (20) with Ef1α-SB100 transgene. After puro selection, SB100-expressing cells were transfected with the pT2-Onc3 transposon plasmid (21). Forty-eight hours later, 1 × 106 cells of SB100 + T2/Onc3 were plated on 10 cm plates. Cells were subsequently treated with vemurafenib (5 μm), vemurafenib (5 μm) and cobimetinib (5 nm), or vehicle (DMSO, 0.2%) 24 hours after plating. Drug or vehicle was renewed every 3 to 4 days. Upon confluency (approximately 3 days after plating), the vehicle plates were collected. Cells treated with vemurafenib or vemurafenib with cobimetinib were collected after approximately 18 or 28 days, respectively.

To determine common transposon insertion sites across plates of resistant cells, genomic DNA from each plate was extracted using the GenElute Mammalian Genome DNA miniprep Kit (Sigma). DNA fragments containing transposon/genome junctions were amplified via ligation-mediated PCR and sequenced using the Illumina Hi-Seq 4000 platform as previously described (22).

Colony staining

The 1 × 106 cells of either A375 SB100 + EGFP or SB100 + T2/Onc3 cells were stained with Coomassie Brilliant Blue after ethanol fixation 25 days after drug treatment (see screen details for concentrations). Colonies were counted with the GelCount Colony Counter (Oxford Optronix).

Creation of candidate overexpression constructs

Gene products mimicking the splice form driven by the SB transposon promoter were amplified from A375 cDNA (BRAF, RAF1, RAC1) or from Transomics Technology human cDNA clones (MCF2, VAV1). Cloned cDNAs were inserted into piggyBac expression vectors containing a human Ef1α promoter along with an IRES-puromycin-polyA cassette. Stable cell lines were obtained by cotransfection of each vector with a piggyBac transposase expression vector via Effectene (Qiagen) transfection reagent. Overexpression was assessed via RT-qPCR and immunoblot. See list of primers and antibodies for specifics.

RNA interference

RHOA and RHOC constructs had a pLKO backbone (Sigma-Aldrich). RAC1 constructs had a pZIP-mCMV vector backbone (Transomics Technologies). A nontargeting shRNA in the appropriate vector backbone was included to produce vector control cell lines. Cells were maintained as stably transduced, polyclonal populations. See Supplementary Information for RNAi-targeting sequences.

Cell culture

The following cell lines were obtained from the ATCC: A375 (CRL-1619), SKMEL28 (HTB-72), A101D (CRL-7898), and 451Lu (CRL-2813). Clones of vemurafenib-resistant A375 were verified as A375 using short tandem repeat analysis performed by IDEXX Bioanalytics. Mycoplasma testing is conducted periodically on all cell lines using a PCR-based method to detect the mycoplasma genome.

All cell lines were grown in DMEM (Gibco) supplemented with penicillin/streptomycin (Gibco) and 10% FBS (Gibco). Spontaneous vemurafenib-resistant clones and populations were created after 2–3 and 4–6 weeks cultured in 3 μm vemurafenib, respectively. Clones were isolated via cloning rings. All spontaneously resistant clones and populations were maintained in media with 3 μm vemurafenib.

CellTiter blue viability assay

Cells were plated in triplicate in 96-well plates at a density of 5 × 102 to 5 × 103 cells per well depending on the cell line. The CellTiter-Blue Viability Assay (Promega) was performed serially on pre- and postinhibitor treated cells. Media containing the specific inhibitor used were renewed every 2 to 6 days. Fold change from day 0 was assessed for each well by comparing pre- and postinhibitor treated cells.

Inhibitors

Unless otherwise indicated, inhibitors and the concentration used in these experiments include vemurafenib (5 μmol/L; Selleckchem), cobimetinib (5 nmol/L; Selleckchem), saracatinib (1 or 2 μmol/L; Selleckchem), FRAX-486 (50 nmol/L; Selleckchem), Fasudil (7 or 10 μmol/L; Sigma), LY3009120 (Selleckchem; 1 μmol/L), imatinib (2 μmol/L; Cayman Chemical), ulixertinib (3 μmol/L; Selleckchem), dabrafenib (200 nmol/L; Selleckchem), encorafenib (500 nmol/L; Selleckchem), trametinib (2 nmol/L; Selleckchem), and binimetinib (50 nmol/L; Selleckchem).

Immunoblotting

RHO protein activation assays were performed using subconfluent 10 cm plates that were treated with the specified inhibitors for the indicated time. Cells were lysed with ice-cold 50 mmol/L Tris Cl, 1% Triton X-100, 0.1% SDS, 10 mmol/L MgCl2 with protease inhibitors, and activated RHO protein was recovered by pulldown with GST-PAK-CRIB fusion protein (23), followed by immunoblotting for active and total RAC1 or CDC42. Protocol was adapted from Pellegrin and Mellor (24). Active GTP was calculated by dividing the quantified signal for RAC1-GTP or CDC42-GTP by the quantified signal for total RAC1 or CDC42, respectively. Total levels of RAC1 and CDC42 were first normalized by the β-actin loading control.

Antibodies used for immunoblotting were as follows: p44/42 MAPK (ERK1/2) (#9102, Cell Signaling Technology), phospho-p44/42 MAPK (T202/Y204) (#9101, Cell Signaling Technology), MEK1 (#2352, Cell Signaling Technology), RAC1 (#610651, BD Transduction), phospho-MEK (Ser 217/221) (#9154, Cell Signaling Technology), phospho-MEK1 (Ser298) (#98195, Cell Signaling Technology), β-tubulin (T8328, Sigma-Aldrich), CDC42 (#610928, BD Transduction), β-actin (622102, BioLegend), BRAF (#MA5-24661, Invitrogen), STAT3 (#4904, Cell Signaling Technology), phospho-STAT (#9145, Cell Signaling Technology), AKT (#9272, Cell Signaling Technology), phospho-AKT (S473) (#4060, Cell Signaling Technology), α-tubulin (12G10, DSHB), VAV1 (HPA001864, Sigma-Aldrich), and MCF2 (#2089, Cell Signaling Technology). Secondary antibodies were IR antibodies from LiCOR or Rockland Inc. Immunoblots were imaged on a LiCOR Odyssey blot imager.

Primer sequences

The following primers were used to clone the indicated coding sequences: MCF2ΔN (5′-GCTTAATTAAATGAAGACTATACAACTCAAGC-3′,5′-TGGCGCGCCTCAATATAGGAG AGCCATCT-3′), MCF2-myc (5′-TTAATTAACCATGCAAGACATCGCCTTCT, 5′-GGCGC GCCTTACAGATCTTCTTCAGAAATAAGTTTTTGTTCCAGATCTTCTTCAGAAATAAGTTTTTGTTCCAGATCTTCTTCAGAAATAAGTTTTTGTTCATATAGGAGAGCCATCTCCGACACAG-3′), RAF1-myc full-length (5′-GCTTAATTAAATGGAGCACATACAG GGAGC-3′, 5′-TTAATTAAATGGAACAAAAGTTAATCTCTGAAGAGGATCTGATGGTCAGCACC ACCCTGC-3′, 5′-ATGGCGCGCCCTAGAAGACAGGCAGCCTCG-3′), BRAFV600EΔN-myc (5′-TTAATTAAATGGAACAAAAGTTAATCTCTGAAGAGGATCTGATGAAAACACTTG GTAGAC-3′, 5′-TGGCGCGCCTCAGTGGACAGGAAACGCAC-3′), VAV1 full-length (5′-CTTAATTAAATGATTGTGCTACCACTGTA-3′, 5′-TGGCGCGCCTCAGCAGTATTCAGA ATAAT-3′), RAC1 full-length (5′-TTAATTAACTGATGCAGGCCATCAAGT-3′, 5′-GGCGC GCCGGGGCTGAGACATTTACAACA-3′), and RAC1P29S (5′-AACCAATG CATTTTCTGGAG-3′, 5′-TATATTCTCCAGAAAATGCATTGGT-3′).

The following primers were used to perform RT-qPCR on the indicated targets: SB100X (5′-AATGGGTCTTCCAACACGAC-3′, 5′-GTGATGGCCACTCCAATACC-3′), Puro (5′-GCAA CCTCCCCTTCTACGAG-3′, 5′-ATCGATATCATGGGGTCGTG-3′), VAV1 (5′-AAACTACA TCGCAGGGCTCA-3′, 5′-CAGGGTTGAGCCGTAGAAAG-3′), and TBP (5′-TTCGGAGAGTTC TGGGATTG-3′, 5′-CTCATGATTACCGCAGCAAA-3′).

The following primers were used to amplify a 5′ portion of the BRAF transcript: 5′-CGGCTCT CGGTTATAAGATGG-3′ and 5′-CTCATTTGTTTCAGTGGACAGG-3′.

Statistical analysis

Statistical methods for some experiments are described within the text. The identification of transposon-induced driver mutations was carried out using a gene-centric common insertion site method previously described (25).

Results

SB mutagenesis drives drug resistance in human melanoma cells

We first engineered human melanoma cells (A375, SKMEL28) to express the hyperactive SB100x transposase (26). Cells stably expressing SB100x (Supplementary Fig. S1A) were subsequently transfected with either the mutagenic pT2-Onc3 transposon vector (21) or a control EGFP expression plasmid (Fig. 1A). Cells were then grown for 48 hours in standard culture conditions to allow the SB100x enzyme to integrate the mutagenic transposons into the genome of the transfected cells (Fig. 1B). Independent plates of mutagenized or control cells were pooled, and 1 × 106 cells were seeded onto 10 cm plates. Cells were placed under drug selection 12 to 24 hours after plating, and genomic DNA was collected for analysis once drug-resistant colonies emerged.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

SB mutagenesis drives drug resistance in A375 through known and novel mechanisms. A, A375 SB100x+ cells were transfected with the pT2-Onc3 transposon before undergoing drug treatment with either vemurafenib (5 μm; n = 75), vemurafenib plus cobimetinib (5 μm and 5 nm, respectively; n = 20), or DMSO control (n = 15). B, The T2-Onc3 transposon can cause gain-of-function or loss-of-function mutations based on orientation of the transposon and location of insertion in the gene locus. C, SB mutagenesis significantly increases the frequency of drug-resistant colonies in A375 cells following long-term exposure to vemurafenib or vemurafenib plus cobimetinib. D, Drug resistance drivers were identified by profiling sites of transposon insertion in resistant cells to find genes that were recurrently overexpressed by transposon insertions.

We generated a variety of mutagenized cell populations for our screen: vemurafenib-treated (5 μmol/L), vemurafenib (5 μmol/L) combined with cobimetinib (5 nmol/L), and control cells treated only with DMSO (i.e., vehicle). Vehicle-treated cells expanded rapidly and were collected approximately 3 days after plating. Drug-resistant colonies emerged with varying kinetics on independent plates but generally appeared in approximately 14 days in vemurafenib alone or approximately 21 days in vemurafenib with cobimetinib. In both conditions, drug-resistant colonies appeared faster and in greater numbers than in nonmutagenized control cells (Fig. 1C), suggesting that transposon-induced mutations drive accelerated drug resistance.

Next, we determined if recurrent transposon-induced mutations could be identified in the resistant colonies. As picking individual colonies would limit the throughput of the approach, we instead harvested independent plates of cells as pooled populations, consistent with genome-wide shRNA and CRISPR protocols. We collected cell populations treated with vemurafenib [A375 (n = 75) and SKMEL28 (n = 16)], vemurafenib plus cobimetinib [A375 (n = 20)], or vehicle [A375 (n = 15) and SKMEL28 (n = 5)]. DNA fragments containing the transposon/genome junctions were amplified via ligation-mediated PCR (LM-PCR) and sequenced using the Illumina HiSeq platform as we have previously reported (22). Raw sequence reads were trimmed, mapped to the human reference genome (GRCh38), and filtered to remove nonclonal insertion events.

The final dataset for each treatment group was analyzed using a modified gene-centric common insertion site analysis to identify genes with a higher rate of transposon insertion than predicted based on a random integration pattern (Supplementary Tables S1–S3; ref. 25). The screens identified a set of four genes for which overexpression is significantly associated with resistance to both vemurafenib and vemurafenib plus cobimetinib (Fig. 1D). Overexpression of four additional genes was associated with resistance to vemurafenib alone. However, some of the differences between the two screens in A375 cell could be attributed to the differences in screen depth between the two drug conditions. Finally, differences in resistance mechanisms were also observed between the two cell lines. Although overexpressions of BRAF and RAF1 were common mechanisms in both A375 and SKMEL28, the DBL family of GEFs (VAV1, MCF2) were unique to A375 while overexpression of PDGFRB was seen only in SKMEL28 cells (Fig. 1D; Supplementary Fig. S2).

Closer inspection of the results revealed differences in the mechanism of resistance for individual genes. Transposon insertion within the promoter or first intron likely drives overexpression of full-length proteins (Supplementary Fig. S2). The T2-Onc3 transposon is also capable of expressing truncated proteins (Fig. 1B). The pattern of insertions in the MCF2 locus suggests that overexpression of either full-length or truncated protein is associated with drug resistance (Supplementary Fig. S2). Transposon insertion appears to preferentially drive expression of N-terminal truncations of BRAF and RAF1 (Supplementary Fig. S2). We verified expression of a truncated BRAFV600E protein in an expanded clone of vemurafenib-resistant A375 cells that have a transposon insertion in the BRAF locus (Supplementary Fig. S3). Importantly, a truncation of the BRAFV600E transcript in human melanoma has previously been shown create a similar N-terminally truncated protein associated with vemurafenib resistance in melanoma patients (6). These mechanistic insights highlight the strength of transposon-based genetic screens to identify more complex mechanisms aside from simple knockdown or overexpression.

Validation of candidate drug resistance drivers

We created vectors to mimic the transposon-induced expression of each gene. Each transgene was then stably expressed in A375 (Supplementary Fig. S4). Drug resistance was evaluated by plating cells at low density in a 96-well culture format and serially measuring the relative viable cell number in each well over multiple days of culture in the presence of drug or vehicle (see Materials and Methods). Using this approach, we confirmed that overexpression of VAV1, MCF2, truncated MCF2 (MCF2ΔN), truncated BRAF (BRAFV600EΔN), or truncated RAF1 (RAF1ΔN) significantly increases resistance to vemurafenib and vemurafenib plus cobimetinib (Fig. 2A), as well as multiple RAF inhibitors (Supplementary Fig. S5). BRAFV600EΔN and RAF1ΔN were the strongest drivers of resistance, with little growth inhibition upon drug treatment. Interestingly, full-length RAF1 did not confer resistance, illustrating the necessity of the N-terminal truncation (Fig. 2A). We also overexpressed RAC1P29S, a hotspot mutation known to drive BRAFi resistance (27), to assess the relative strength of our novel resistance drivers. RAC1P29S performed similarly to VAV1 and MCF2 in our assay (Fig. 2A). Notably, none of the candidates identified in our screen (VAV1, MCF2, MCF2ΔN, BRAFV600EΔN, RAF1ΔN) increased proliferation of A375 cells in the absence of drug (Supplementary Fig. S6A). However, overexpression of RAC1, RAC1P29S, and full-length RAF1 did increase the growth rate significantly (Supplementary Fig. S6B).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Validation of candidate drug resistance drivers in a panel of human melanoma cell lines. Growth of various engineered cell lines was assessed via CellTiterBlue (Promega). The ability of candidates to increase resistance to vemurafenib and vemurafenib plus cobimetinib varied in A375 (A), 451Lu (B), A101D (C), and SKMEL28 (D). *, corrected P ≤ 0.05; fold-change is relative to vector cells treated with 5 μm vemurafenib ± 5 nmol/L cobimetinib for each independent assay. E, A brief 12-hour vemurafenib treatment shows distinct patterns in MEK phosphorylation between the DBL- and RAF-driven mechanisms of resistance. F and G, Overexpression of either MCF2ΔN or VAV1 increases the level of CDC42-GTP (F) and RAC1-GTP (G).

Performance of novel MAPKi resistance drivers in independent melanoma cell lines

Demonstration that our candidates confer drug resistance in A375 cells supports the validity of our genetic screen approach. To determine if the same mechanisms provide resistance in other MAPKi-sensitive BRAFV600E-mutant human melanoma cell lines, we generated populations of 451Lu, A101D, and SKMEL28 stably expressing each candidate gene and subjected them to the same 96-well growth assay used for A375 cells. We determined that 451Lu had a similar resistance profile, with all of the candidates providing resistance to both vemurafenib and vemurafenib plus cobimetinib (Fig. 2B). All candidates were able to drive MAPKi resistance in A101D, with the exception of VAV1 (Fig. 2C). Consistent with the findings from our genetic screen, only truncated BRAF and RAF1 were able to drive MAPKi resistance in SKMEL28, although the degree of resistance was much weaker than in A375 or 451Lu (Fig. 2D).

DBL family members VAV1 and MCF2 drive MAPKi resistance through RAC1

Although BRAF truncation is a well-established mechanism of MAPKi resistance, the mechanism of VAV1- and MCF2-driven resistance has not previously been investigated. Both VAV1 and MCF2 are members of the DBL family of guanine exchange factors (GEF) and function by activating members of the RHO family of small GTPases (28). We evaluated AKT, STAT3, and MAPK signaling in the presence and absence of vemurafenib because MAPK reactivation is frequently observed in vemurafenib-resistant tumors (29, 30). Overexpression of VAV1 or MCF2 did not alter activation of AKT and STAT3 relative to the vector control (Supplementary Fig. S7). Vemurafenib abolished phosphorylated ERK in vector control A375 cells (Fig. 2E). In contrast, phosphorylated ERK was partially restored in cells expressing MCF2, VAV1, and BRAFΔN, and completely restored in RAF1ΔN. Next, we examined MEK phosphorylation upstream of ERK in cells expressing each candidate. As expected, expression of either truncated RAF1 or BRAFV600E partially or completely restored phosphorylation at the RAF-controlled serine 217 on MEK in cells exposed to vemurafenib (Fig. 2E). In contrast, overexpression of VAV1 or MCF2 promoted maintenance of phosphorylation of serine 298 on MEK, an indication of elevated PAK activity in these cells (Fig. 2E). This observation is consistent with the finding that increased PAK signaling can drive acquired resistance to MAPK inhibitors in melanoma (31). Taken together, the dichotomous activation of MEK suggests a mechanistic distinction between resistance driven by RAF truncation versus DBL GEF overexpression.

We next investigated how VAV1 and MCF2 work to reestablish ERK signaling in the presence of MAPKi. Both VAV1 and MCF2 were originally identified as oncogenes that act as GEFs for the RHO family (32, 33). Although the precise RHO family targets for these proteins are not entirely clear, it is generally accepted that VAV1 and MCF2 have activity for RHO, RAC1, and CDC42 (34). We verified that overexpression of VAV1 and MCF2ΔN increases the levels of GTP-bound CDC42 (Fig. 2F) and RAC1 (Fig. 2G), particularly in the presence of vemurafenib.

Overexpression of VAV1 and MCF2 could drive vemurafenib resistance via a RAC1/CDC42/PAK pathway and/or by signaling via RHO through RHO-associated kinase (ROCK; Fig. 3A). To distinguish between these mechanisms, we assessed if inhibition of either ROCK or PAK alters vemurafenib resistance driven by MCF2 or VAV1. First, we identified drug concentrations that had a minimal impact on cell growth in the absence of vemurafenib (Supplementary Fig. S6C and S6D). Vemurafenib resistance was consistently reduced by the addition of the PAK inhibitor FRAX-486 but not by the ROCK inhibitor Fasudil (Fig. 3B). Inhibition of PAK also prevented the emergence of spontaneous vemurafenib resistance in long-term A375 cultures (Fig. 3C; Supplementary Fig. S6E), whereas ROCK inhibition did not have a substantial effect on long-term resistance (Fig. 3D; Supplementary Fig. S6F). Finally, knockdown of RAC1 in parental A375 showed synergistic cell killing in long-term culture with vemurafenib (Fig. 3E) but not vehicle (i.e., DMSO; Supplementary Fig. S6G). Conversely, knockdown of CDC42 did not show the same effect with long-term vemurafenib treatment (Supplementary Fig. S6H). Finally, knockdown of RHOA or RHOC did show any significant effects in parental A375 under standard culture conditions (Fig. 3F; Supplementary Fig. S6I).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Evaluation of DBL GEF downstream signaling mechanism in A375 cells. A, DBL family members have exchange activity for members of the RHO and RAC1/CDC42 family, each having distinct signaling mechanisms. B, The addition of the PAK inhibitor FRAX-486 (50 nmol/L) is able to reduce vemurafenib resistance driven by MCF2 and VAV1. However, the ROCK inhibitor Fasudil did not show a consistent effect *, adjusted P < 0.001 relative to vemurafenib alone. C, Long-term treatment of A375 cells with a PAK inhibitor (FRAX-486) is capable of converting vemurafenib-induced cytostasis to cell killing, whereas long-term treatment with a ROCK inhibitor (Fasudil; D) shows a trend of enhanced growth in vemurafenib. E and F, Knockdown of RAC1 using three independent shRNAs converts vemurafenib-induced cytostasis to cell killing in long-term cultures of A375 (E), whereas knockdown of either RHOA or RHOC does not affect cell growth with long-term vemurafenib treatment (F).

The role of novel vemurafenib resistance drivers in spontaneous melanoma progression

Although N-terminal truncations of BRAF have previously been associated with vemurafenib resistance in patients (6), the other candidate resistance drivers we identified have not. We sought to determine if alterations in these candidates would arise spontaneously to provide vemurafenib resistance in the absence of experimentally induced mutagenesis. First, we generated A375 cells with spontaneously acquired resistance after long-term culture in 3 μmol/L vemurafenib. Two phenotypically distinct populations of resistant cells were derived from this process (Supplementary Fig. S8A–S8D). The first vemurafenib-resistant cells that grew out of these cultures were derived from colonies showing a compact morphology. We isolated nine such colonies from eight different populations of vemurafenib-treated A375 cells. We also isolated three independent populations of the second class of vemurafenib-resistant A375 cells, which grew more diffusely and slowly. Importantly, cell populations of this second class do not generate the more rapidly growing colonies of vemurafenib-resistant cells. We hypothesized that the two distinct cell morphologies represent distinct drug resistance mechanisms.

We first looked for evidence of BRAFV600E alterations in the nine colonies and three populations of vemurafenib-resistant A375 cells by performing Western blotting using an antibody specific to BRAFV600E isoform (Supplementary Fig. S8E). No alterations were detected in the vemurafenib-resistant populations (VRP1-3), but six of the nine vemurafenib-resistant clones showed evidence of BRAFV600E protein alteration (Supplementary Fig. S8E). Two of the six clones express a truncated form of BRAFV600E that is approximately 40 kD, consistent with truncated forms of BRAFV600E known to drive vemurafenib resistance in patients (6). We confirmed that this truncation is caused by an aberrant splice event between BRAF exons 1 and 11 (Supplementary Fig. S8F). An additional four vemurafenib-resistant clones appear to express a BRAFV600E fusion that is approximately 120 kD in size (Supplementary Fig. S8E). Consistent with the presence of BRAF alterations, the level of pMEK-S217 (RAF phosphorylation site) was elevated in all but one of the vemurafenib-resistant clones, even in the presence of vemurafenib (Supplementary Fig. S8G). A similar approach did not provide evidence of RAF1 alterations in vemurafenib-resistant clones and populations (Supplementary Fig. S8H).

Next, we performed Western blotting to determine if either MCF2 or VAV1 expression is altered in the A375 cells with acquired spontaneous resistance to vemurafenib. Although no changes in MCF2 expression were detected, VAV1 protein expression is significantly higher in all three vemurafenib-resistant cell populations (VRP1-3; Supplementary Fig. S8F). Levels of pMEK-S298 were also elevated in these populations, consistent with our findings that VAV1 drives vemurafenib resistance via RAC1/PAK signaling (Supplementary Fig. S8F). The increased levels of VAV1 protein were not accompanied by increased mRNA expression (Supplementary Table S4), suggesting that the increased protein level is achieved through posttranscriptional regulation (e.g., increased protein stability). These results provide corroborating evidence that VAV1 is involved in spontaneous vemurafenib resistance.

Finally, we used several independent approaches to evaluate the novel candidate vemurafenib resistance drivers (RAF1, MCF2, VAV1) in human melanoma samples. First, we analyzed 159 BRAFV600E-mutant melanoma samples present in The Cancer Genome Atlas cutaneous melanoma project (TCGA-SKCM). Although BRAFi response data are not available for these samples, we evaluated the expression levels of candidates identified in our vemurafenib resistance screen that are predicted to drive resistance when overexpressed (Supplementary Table S1). Expression of VAV1 initially appeared to be elevated in a portion of melanomas in the TCGA-SKCM panel. However, it has been shown that VAV1 is highly expressed in lymphocytes. Indeed, further analysis showed that VAV1 expression correlated strongly with T-cell markers, suggesting that the majority of VAV1 expression in the TCGA-SKCM biopsies was contributed by infiltrating lymphocytes. Therefore, we were unable to evaluate VAV1 expression using transcriptome sequencing data derived from bulk analysis of melanoma biopsies. However, each of the remaining genes showed overexpression (z-score ≥ 2) in a subset of TCGA-SKCM samples. Overall, 56 samples (∼37%) exhibit overexpression of one or more genes from the set (Fig. 4A).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Role of candidate drug resistance drivers in human melanoma patients. A, Evaluation of the top vemurafenib resistance drivers in primary cutaneous melanoma samples from the TCGA-SKCM project. Over 37% of patient samples show increased expression of at least one of the candidate resistance drivers. B, Expression of the candidate resistance drivers in an independent RNA-seq data set obtained from a cohort of 18 patients, each with a diagnostic sample (i.e., pretreatment) matched with a second sample taken after progression on BRAFi therapy (35). Top, six of the 18 patients had elevated expression of at least one resistance driver in the pretreatment sample. Additional 8 patients showed overexpression of at least one candidate resistance driver in the progression samples.

Unfortunately, the majority of samples in the TCGA-SKCM panel lack sufficient patient treatment histories to allow correlation of gene expression with response to MAPKi treatment. Therefore, we evaluated an RNA-seq data set obtained from a collection of primary melanomas with matched progression samples (35). Hugo and colleagues obtained RNA-seq data on pretreatment melanoma samples taken from eighteen patients along with matched biopsies obtained after patients had progressed during treatment with a BRAFi alone or a BRAFi/MEKi combination.

We evaluated expression of the vemurafenib resistance gene set using this RNA-seq data to identify patterns in expression that correlate with treatment response. Pretreatment tumor samples showed an expression pattern similar to that observed in the TCGA-SKCM panel with 6 of 18 (∼33%) patient samples showing overexpression of at least one of the resistance driver genes (Fig. 4B). We then examined expression of the gene set in progression samples (n = 33) taken from patients whose primary tumor sample did not initially exhibit overexpression. Interestingly, 8 of 12 patients (∼66%) appeared to acquire overexpression of at least one resistance driver (Fig. 4B).

Src inhibition blocks MAPKi resistance driven by Rac signaling

Our forward genetic screen identified multiple members of the DBL family of guanine nucleotide exchange factors (MCF2, VAV1) as novel drivers of vemurafenib resistance (Fig. 1D). Unfortunately, there are no drugs available that can inhibit the activity of GEFs or their small GTPase targets. However, the activity of the DBL family is regulated by a variety of upstream signals. For instance, both VAV1 and MCF2 can be activated by Src-mediated phosphorylation (36, 37). Interestingly, several prior publications have shown that the SRC inhibitor saracatinib exhibits synergism with vemurafenib, although this synergism was not attributed to SRC's role in regulating VAV1 or MCF2 (38, 39).

We hypothesized that inhibiting SRC using saracatinib would block vemurafenib resistance driven by MCF2 and VAV1 overexpression. We performed a short-term growth assay in A375 cells engineered to express each resistance driver. As we had seen previously, expression of MCF2, VAV1, RAC1P29S, RAF1ΔN, and BRAFΔN all drove growth of A375 cells in vemurafenib and vemurafenib combined with cobimetinib (Fig. 5A). As predicted, the addition of saracatinib to vemurafenib inhibited the growth of cells overexpressing MCF2 and VAV1. Surprisingly, saracatinib was also able to block vemurafenib resistance driven by RAC1P29S expression, suggesting that the P29S activation mechanism may still depend on SRC-dependent GEF activity. In all of these cases, the combination of vemurafenib and saracatinib not only blocked growth but caused cell death over the course of the 12-day assay (Fig. 5A). Although vemurafenib resistance driven by truncated BRAFV600E and RAF1 was modestly reduced by the addition of saracatinib, these cells were still able to grow in the presence of both drugs (Fig. 5A).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

An SRC inhibitor, saracatinib, shows synergistic cytotoxicity in combination with a next-generation Raf inhibitor. A, The resistance driven by MCF2 and VAV1 can be overcome with the addition of saracatinib. However, saracatinib cannot block resistance driven by the expression of truncated BRAFV600E or RAF1. *, adjusted P < 0.05 relative to vemurafenib alone). B, The addition of saracatinib does not uniformly increase vemurafenib response in a panel of BRAFV600E-mutant human melanoma cell lines. *, adjusted P < 0.0001. C, Saracatinib, but not imatinib, is able to induce cytotoxicity with vemurafenib in A375 cells. D, The addition of saracatinib is able to induce cytotoxicity in cells that have acquired spontaneous resistance to vemurafenib. E and F, Saracatinib treatment reverses the increase in RAC1-GTP (E) and CDC42-GTP (F) observed in A375 cells expressing MCF2 and VAV1. G, The resistance candidates perform differently in response to treatment with LY3009120, a next-generation Raf inhibitor. Importantly, the combination of LY3009120 and saracatinib is cytotoxic to all cell populations (note independent y-axes).

We next determined if the addition of saracatinib increases vemurafenib sensitivity across a panel of BRAFV600E-mutant melanoma cell lines. Saracatinib exhibited synergistic cell killing in short-term growth assays with vemurafenib in three cell lines (A375, 451Lu, and COLO858) that was not observed in two additional cell lines (SKMEL28 and WM2664; Fig. 5B). Prior work has shown that saracatinib can inhibit the ABL kinases, albeit with lower activity (40). We performed long-term growth assays using A375 cells to determine if ABL kinase inhibition using imatinib had a similar effect as saracatinib. Two different concentrations of saracatinib were able to block the emergence of A375 cells with spontaneous vemurafenib resistance, whereas imatinib treatment did not affect the acquisition of spontaneous resistance (Fig. 5C). Moreover, saracatinib was also able to kill A375 cells after spontaneous vemurafenib resistance developed (Fig. 5D). Finally, the addition of saracatinib reverses the activation of RAC1 and CDC42 observed in A375 cells overexpressing MCF2ΔN and VAV1, consistent with SRC acting upstream as an activator of DBL GEF activity (Fig. 5E and F).

As previously mentioned, N-terminal truncations of BRAF have been shown to drive resistance to vemurafenib in patients (6). Subsequent work has shown that BRAF truncations and fusions can drive resistance by functioning as constitutively active dimers, which cannot be blocked by vemurafenib (6, 41, 42). However, next-generation BRAFis (e.g., LY3009120) have been developed that are capable of inhibiting Raf dimers in addition to monomeric BRAFV600E (41, 42). These dimer-blocking compounds are active against RAS mutant cells because they are able to inhibit RAS-dependent BRAF dimers, which cannot be blocked by vemurafenib (41). Nevertheless, cells with intrinsic resistance to LY3009120 have been reported (41), suggesting that not all BRAFi resistance mechanisms act by enforcing BRAF dimerization.

We tested the ability of the dimer-blocking drug LY3009120 to inhibit growth of A375 cells expressing each of the resistance drivers identified by our screen. As expected, LY3009120 significantly reduced the growth of cells expressing either truncated BRAFV600E or RAF1 (Fig. 5G). However, expression of MCF2ΔN and VAV1 were still able to drive proliferation in the presence of LY3009120, suggesting that the Rac-driven resistance mechanism may not rely on RAF dimerization. Nevertheless, the combination of LY3009120 and saracatinib was able to induce cytotoxicity in all cell populations, suggesting that this drug combination can thwart both resistance mechanisms identified by our genetic screen.

Each drug resistance driver was then tested in A375 cells using a panel of MAPK inhibitors including four RAF inhibitors, three MEKis, and an ERK inhibitor (Fig. 6). These experiments demonstrated several trends. First, none of the currently approved RAF (vemurafenib, dabrafenib, or encorafenib) or MEKis (cobimetinib, trametinib, or binimetinib) were able to control the growth driven by the resistance drivers as mono- or combination therapies (i.e., RAFi + MEKi). The ERK inhibitor ulixertinib was also unable to suppress the growth of cells overexpressing MCF2 and VAV1. However, although saracatinib (SRCi) alone had little effect on cell growth, the combination of RAFi with SRCi was much more effective. The combination of LY3009120 with saracatinib was the most effective drug combination (Fig. 6). As previously observed, RAFi with SRCi was effective in suppressing the growth of 451Lu (Fig. 6) but not A101D or SKMEL28 (Fig. 6).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Performance of candidate resistance drivers across a panel of MAPK inhibitors. The ability of candidates to increase resistance to the indicated drugs is shown as a relative colorimetric endpoint (see legend). Each candidate was tested in the indicated drug combination in A375, 451Lu, A101D, and SKMEL28. Rows represent independent populations engineered to express the indicated candidate drug resistance driver shown on the left.

Discussion

We report here the results of a forward genetic screen designed to identify gain-of-function mutations that drive resistance to targeted kinase inhibitors (Fig. 1A). Our approach utilized a hyperactive form of the SB transposase to increase the mutagenesis efficiency in cultured cells, allowing us to establish a simple screening method that can be easily replicated and applied in other contexts. We used this approach to identify novel drug resistance drivers in a variety of conditions. The screen results suggest that although drug resistance can be driven by common mechanisms (e.g., BRAF/RAF1 truncation), there are diverse mechanisms that are unique to specific drug combinations or cell lines (Fig. 1D). However, additional work is needed to validate the broader findings from our screens.

Prior work has already shown that N-terminal truncation of BRAF is associated with vemurafenib resistance in patients (6). We observed transposon clustering in the middle of the BRAF locus predicted to drive expression of a similar truncated isoform of BRAF (Supplementary Fig. S2C), strongly supporting the relevance of our screening method. Furthermore, we observed a similar pattern of insertions in the RAF1 locus (Supplementary Fig. S2D), and we subsequently showed that expression of an N-terminal truncated isoform of RAF1 can drive drug resistance to a level comparable with that of truncated BRAF in a collection of vemurafenib-sensitive human melanoma cells (Fig. 2). This finding suggests that RAF1 N-terminal truncation through intragenic deletion or gene fusion events could account for resistance in some human melanoma patients.

The other major finding from our genetic screens is that overexpression of the DBL-family guanine nucleotide exchange factors MCF2 and VAV1 can drive resistance to MAPK inhibition (Fig. 6). Although the DBL family of GEFs act on both RHO and Rac proteins (28), we show that the DBL-driven resistance likely signals through RAC1/PAK rather than RHO/ROCK (Fig. 3). A prior study showed that gain-of-function mutations in RAC1 are associated with vemurafenib resistance (27), and overexpression of the DBL GEFs is another mechanism through which RAC1 can be activated to drive resistance. To our knowledge, this is the first time the DBL family has been implicated in vemurafenib resistance.

All of the candidate resistance drivers we chose for further study significantly increased drug resistance in at least two independent sensitive human melanoma cell lines (Fig. 2). To further establish the relevancy of our candidates, we also evaluated the candidate drivers for a role in vemurafenib resistance by studying 12 independent clonal populations of A375 cells that had spontaneously acquired vemurafenib resistance, finding evidence of BRAFV600E alterations (Supplementary Fig. S8E) and increased VAV1 protein expression (Supplementary Fig. S8F). These two mechanisms could account for resistance in 75% of the spontaneously resistant clones and populations (9 of 12). The remaining three clones (c2.2.1, c2.2.2, and c7.1) all show increased phosphorylation of MEK1 on serine 298, indicating increased PAK activity in these cells (Supplementary Fig. S8F). This observation suggests that these populations have active RAC1/CDC42 signaling, consistent with the DBL GEF–driven mechanism we identified.

We also evaluated the relevancy of our screen results with sequencing data from melanoma patients who had progressed on MAPKi treatment (35). Analysis of this previously published dataset of 18 patients (GSE65186) revealed that over 75% of patients (14 of 18) had overexpression of at least one resistance gene in the initial diagnostic biopsy or had acquired overexpression in at least one progression sample (Fig. 4). This suggests that melanomas with overexpression of resistance drivers at the time of treatment are less likely to respond to MAPKi and that acquired overexpression of our candidates during treatment may drive progression. However, analysis of a larger cohort is needed to determine if these trends are significant.

The identification of DBL GEFs as drivers of vemurafenib resistance also provides a direct mechanistic link between SRC and vemurafenib resistance. Prior studies have implicated SRC in mediating vemurafenib resistance (39, 43–45), but none have tied SRC mechanistically to a pathway known to drive vemurafenib resistance in human melanoma. We have shown here that vemurafenib resistance driven by MCF2 and VAV1 can be blocked using the SRC family kinase inhibitor saracatinib (Fig. 5A). Thus, one mechanism by which SRC can drive vemurafenib resistance is through the activation of RAC1/PAK through DBL GEFs such as MCF2 and VAV1. Importantly, we show that saracatinib blocks vemurafenib resistance driven by both overexpression of DBL GEFs and by overexpression of RAC1P29S, a previously identified mutation associated with vemurafenib resistance in cutaneous melanoma (27). Furthermore, we show for the first time that saracatinib can prevent the emergence of spontaneous vemurafenib resistance in long-term cultures of A375 cells (Fig. 5D).

Prior work has shown that N-terminal truncation of either BRAFV600E or RAF1 promotes the formation of Ras-independent Raf dimers that cannot be inhibited with vemurafenib or dabrafenib (6, 41, 42). However, next-generation Raf inhibitors have been developed that can block both monomeric and dimeric Raf (42). Importantly, one of these compounds, LY3009120, has recently been tested in a phase I clinical trial (NCT02014116). As predicted, we show that LY3009120 is more effective than vemurafenib at inhibiting the proliferation of melanoma cells expressing truncated BRAFV600E or RAF1 (Fig. 5G). It is important to note, however, that LY3009120 was unable to completely inhibit the growth of these cells at the concentration used in our experiment. Nevertheless, the addition of saracatinib to LY3009120 was able to block growth of cells driven by all mechanisms we validated from our forward genetic screen (Figs. 5G and 6).

We also evaluated a panel of MAPK inhibitors for their ability to control the growth driven by the various resistance drivers (Fig. 6). As expected, each driver was able to support cell growth in the presence of either RAF or MEK inhibitors alone or in combination. Interestingly, an ERK inhibitor was able to control growth driven by truncated BRAF or RAF1 but could not completely inhibit the growth driven by the DBL GEFs. This suggests that the DBL-driven drug resistance mechanism may involve a mechanism independent of MAPK. Nevertheless, DBL-driven resistance could be controlled with the addition of saracatinib.

Based on our findings, we propose a model of BRAFi resistance that involves two distinct mechanisms: one utilizing RAF/BRAFV600E truncation (Fig. 7A and B) and a second involving RAC1 activation via the DBL family members MCF2 and VAV1 (Fig. 7C and D). Some aspects of the model will require additional studies. For instance, we do not yet understand the connection between BRAFi treatment and SRC signaling. It is possible that MAPK inhibition leads to changes in SRC activity through posttranslational modification of the SRC kinases and/or through transcriptional feedback mechanisms. It is also important to acknowledge that not all melanoma cell lines are responsive to saracatinib treatment, suggesting that not all melanomas show this BRAFi-induced SRC alteration. Consistent with this idea, the DBL GEFs are unable to drive BRAFi resistance in some melanoma cell lines.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Model of drug resistance mechanisms. We have shown that N-terminal truncation of either RAF1 or BRAFV600E leads to vemurafenib resistance (A) that can be overcome with a pan-RAF inhibitor (e.g., LY3009120; B). Our genetic screen identified a DBL-RAC1-PAK resistance mechanism that can drive proliferation in the presence of either BRAF inhibitor (e.g., vemurafenib) or a pan-RAF inhibitor (C). However, the combination of an SRC inhibitor (e.g., saracatinib) and a pan-RAF inhibitor can block both mechanisms identified by our screen (D).

Many of the questions raised by our experiments can be addressed by employing our forward genetic screening method in independent melanoma cell lines to elucidate shared and unique drivers of MAPKi resistance. Beyond resistance to vemurafenib, we can easily utilize our approach to identify novel mechanisms of resistance to next-generation inhibitors, such as LY3009120 (41, 46). Lastly, our screening approach could inform the effective utilization of other targeted agents in cancer.

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.

Authors' Contributions

Conception and design: C.R. Feddersen, J.L. Schillo, C.S. Stipp, A.J. Dupuy

Development of methodology: C.R. Feddersen, J.L. Schillo, H.R. Vaughn, L.S. Wadsworth, A.J. Dupuy

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): C.R. Feddersen, J.L. Schillo, L.S. Wadsworth, A.P. Voigt, E.Y. Zhu, B.M. Jennings, S.A. Mullen, J. Bobera, J.D. Riordan, C.S. Stipp

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): C.R. Feddersen, J.L. Schillo, A. Varzavand, B.M. Jennings, J.D. Riordan, C.S. Stipp, A.J. Dupuy

Writing, review, and/or revision of the manuscript: C.R. Feddersen, J.L. Schillo, A.P. Voigt, E.Y. Zhu, J.D. Riordan, C.S. Stipp, A.J. Dupuy

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): C.R. Feddersen, A. Varzavand, A.J. Dupuy

Study supervision: C.S. Stipp, A.J. Dupuy

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Footnotes

  • Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org/).

  • Cancer Res 2019;79:5074–87

  • Received January 24, 2019.
  • Revision received June 5, 2019.
  • Accepted August 6, 2019.
  • Published first August 15, 2019.
  • ©2019 American Association for Cancer Research.

References

  1. 1.↵
    1. Cronin KA,
    2. Lake AJ,
    3. Scott S,
    4. Sherman RL,
    5. Noone AM,
    6. Howlader N,
    7. et al.
    Annual report to the nation on the status of cancer, part I: national cancer statistics. Cancer 2018;124:2785–800.
    OpenUrl
  2. 2.↵
    1. Davies H,
    2. Bignell GR,
    3. Cox C,
    4. Stephens P,
    5. Edkins S,
    6. Clegg S,
    7. et al.
    Mutations of the BRAF gene in human cancer. Nature 2002;417:949–54.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Chapman PB,
    2. Hauschild A,
    3. Robert C,
    4. Haanen JB,
    5. Ascierto P,
    6. Larkin J,
    7. et al.
    Improved survival with vemurafenib in melanoma with BRAF V600E mutation. N Engl J Med 2011;364:2507–16.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Corcoran RB,
    2. Dias-Santagata D,
    3. Bergethon K,
    4. Iafrate AJ,
    5. Settleman J,
    6. Engelman JA
    . BRAF gene amplification can promote acquired resistance to MEK inhibitors in cancer cells harboring the BRAF V600E mutation. Sci Signal 2010;3:ra84.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Long GV,
    2. Fung C,
    3. Menzies AM,
    4. Pupo GM,
    5. Carlino MS,
    6. Hyman J,
    7. et al.
    Increased MAPK reactivation in early resistance to dabrafenib/trametinib combination therapy of BRAF-mutant metastatic melanoma. Nat Commun 2014;5:5694.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Poulikakos PI,
    2. Persaud Y,
    3. Janakiraman M,
    4. Kong X,
    5. Ng C,
    6. Moriceau G,
    7. et al.
    RAF inhibitor resistance is mediated by dimerization of aberrantly spliced BRAF(V600E). Nature 2011;480:387–90.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Ascierto PA,
    2. McArthur GA,
    3. Dreno B,
    4. Atkinson V,
    5. Liszkay G,
    6. Di Giacomo AM,
    7. et al.
    Cobimetinib combined with vemurafenib in advanced BRAF(V600)-mutant melanoma (coBRIM): updated efficacy results from a randomised, double-blind, phase 3 trial. Lancet Oncol 2016;17:1248–60.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Long GV,
    2. Stroyakovskiy D,
    3. Gogas H,
    4. Levchenko E,
    5. de Braud F,
    6. Larkin J,
    7. et al.
    Combined BRAF and MEK inhibition versus BRAF inhibition alone in melanoma. N Engl J Med 2014;371:1877–88.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Konermann S,
    2. Brigham MD,
    3. Trevino AE,
    4. Joung J,
    5. Abudayyeh OO,
    6. Barcena C,
    7. et al.
    Genome-scale transcriptional activation by an engineered CRISPR-Cas9 complex. Nature 2015;517:583–8.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Shalem O,
    2. Sanjana NE,
    3. Hartenian E,
    4. Shi X,
    5. Scott DA,
    6. Mikkelson T,
    7. et al.
    Genome-scale CRISPR-Cas9 knockout screening in human cells. Science 2014;343:84–7.
    OpenUrlAbstract/FREE Full Text
  11. 11.↵
    1. Whittaker SR,
    2. Theurillat JP,
    3. Van Allen E,
    4. Wagle N,
    5. Hsiao J,
    6. Cowley GS,
    7. et al.
    A genome-scale RNA interference screen implicates NF1 loss in resistance to RAF inhibition. Cancer Discov 2013;3:350–62.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Johannessen CM,
    2. Boehm JS,
    3. Kim SY,
    4. Thomas SR,
    5. Wardwell L,
    6. Johnson LA,
    7. et al.
    COT drives resistance to RAF inhibition through MAP kinase pathway reactivation. Nature 2010;468:968–72.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Johannessen CM,
    2. Johnson LA,
    3. Piccioni F,
    4. Townes A,
    5. Frederick DT,
    6. Donahue MK,
    7. et al.
    A melanocyte lineage program confers resistance to MAP kinase pathway inhibition. Nature 2013;504:138–42.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. O'Donnell KA
    . Advances in functional genetic screening with transposons and CRISPR/Cas9 to illuminate cancer biology. Curr Opin Genet Dev 2018;49:85–94.
    OpenUrl
  15. 15.↵
    1. Chen HJ,
    2. Wei Z,
    3. Sun J,
    4. Bhattacharya A,
    5. Savage DJ,
    6. Serda R,
    7. et al.
    A recellularized human colon model identifies cancer driver genes. Nat Biotechnol 2016;34:845–51.
    OpenUrlCrossRef
  16. 16.↵
    1. Guo Y,
    2. Updegraff BL,
    3. Park S,
    4. Durakoglugil D,
    5. Cruz VH,
    6. Maddux S,
    7. et al.
    Comprehensive ex vivo transposon mutagenesis identifies genes that promote growth factor independence and leukemogenesis. Cancer Res 2016;76:773–86.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Kodama T,
    2. Newberg JY,
    3. Kodama M,
    4. Rangel R,
    5. Yoshihara K,
    6. Tien JC,
    7. et al.
    Transposon mutagenesis identifies genes and cellular processes driving epithelial-mesenchymal transition in hepatocellular carcinoma. Proc Natl Acad Sci U S A 2016;113:E3384–93.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Perna D,
    2. Karreth FA,
    3. Rust AG,
    4. Perez-Mancera PA,
    5. Rashid M,
    6. Iorio F,
    7. et al.
    BRAF inhibitor resistance mediated by the AKT pathway in an oncogenic BRAF mouse melanoma model. Proc Natl Acad Sci U S A 2015;112:E536–45.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. Kas SM,
    2. de Ruiter J,
    3. Schipper K,
    4. Schut E,
    5. Bombardelli L,
    6. Wientjens E,
    7. et al.
    Transcriptomics and transposon mutagenesis identify multiple mechanisms of resistance to the FGFR inhibitor AZD4547. Cancer Res 2018;78:5668–79.
    OpenUrlCrossRef
  20. 20.↵
    1. Yusa K,
    2. Zhou L,
    3. Li MA,
    4. Bradley A,
    5. Craig NL
    . A hyperactive piggyBac transposase for mammalian applications. Proc Natl Acad Sci U S A 2011;108:1531–6.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    1. Dupuy AJ,
    2. Rogers LM,
    3. Kim J,
    4. Nannapaneni K,
    5. Starr TK,
    6. Liu P,
    7. et al.
    A modified sleeping beauty transposon system that can be used to model a wide variety of human cancers in mice. Cancer Res 2009;69:8150–6.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Riordan JD,
    2. Drury LJ,
    3. Smith RP,
    4. Brett BT,
    5. Rogers LM,
    6. Scheetz TE,
    7. et al.
    Sequencing methods and datasets to improve functional interpretation of sleeping beauty mutagenesis screens. BMC Genomics 2014;15:1150.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Benard V,
    2. Bohl BP,
    3. Bokoch GM
    . Characterization of rac and cdc42 activation in chemoattractant-stimulated human neutrophils using a novel assay for active GTPases. J Biol Chem 1999;274:13198–204.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Pellegrin S,
    2. Mellor H
    . Rho GTPase activation assays. Curr Protoc Cell Biol 2008;Chapter 14:Unit 14 8.
  25. 25.↵
    1. Brett BT,
    2. Berquam-Vrieze KE,
    3. Nannapaneni K,
    4. Huang J,
    5. Scheetz TE,
    6. Dupuy AJ
    . Novel molecular and computational methods improve the accuracy of insertion site analysis in Sleeping Beauty-induced tumors. PLoS One 2011;6:e24668.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Mates L,
    2. Chuah MK,
    3. Belay E,
    4. Jerchow B,
    5. Manoj N,
    6. Acosta-Sanchez A,
    7. et al.
    Molecular evolution of a novel hyperactive Sleeping Beauty transposase enables robust stable gene transfer in vertebrates. Nat Genet 2009;41:753–61.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Watson IR,
    2. Li L,
    3. Cabeceiras PK,
    4. Mahdavi M,
    5. Gutschner T,
    6. Genovese G,
    7. et al.
    The RAC1 P29S hotspot mutation in melanoma confers resistance to pharmacological inhibition of RAF. Cancer Res 2014;74:4845–52.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Zheng Y
    . Dbl family guanine nucleotide exchange factors. Trends Biochem Sci 2001;26:724–32.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Lim SY,
    2. Menzies AM,
    3. Rizos H
    . Mechanisms and strategies to overcome resistance to molecularly targeted therapy for melanoma. Cancer 2017;123:2118–29.
    OpenUrl
  30. 30.↵
    1. Luebker SA,
    2. Koepsell SA
    . Diverse mechanisms of BRAF inhibitor resistance in melanoma identified in clinical and preclinical studies. Front Oncol 2019;9:268.
    OpenUrl
  31. 31.↵
    1. Lu H,
    2. Liu S,
    3. Zhang G,
    4. Bin Wu,
    5. Zhu Y,
    6. Frederick DT,
    7. et al.
    PAK signalling drives acquired drug resistance to MAPK inhibitors in BRAF-mutant melanomas. Nature 2017;550:133–6.
    OpenUrlCrossRef
  32. 32.↵
    1. Eva A,
    2. Aaronson SA
    . Isolation of a new human oncogene from a diffuse B-cell lymphoma. Nature 1985;316:273–5.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Katzav S,
    2. Martin-Zanca D,
    3. Barbacid M
    . vav, a novel human oncogene derived from a locus ubiquitously expressed in hematopoietic cells. EMBO J 1989;8:2283–90.
    OpenUrlPubMed
  34. 34.↵
    1. Jaiswal M,
    2. Dvorsky R,
    3. Ahmadian MR
    . Deciphering the molecular and functional basis of Dbl family proteins: a novel systematic approach toward classification of selective activation of the Rho family proteins. J Biol Chem 2013;288:4486–500.
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    1. Hugo W,
    2. Shi H,
    3. Sun L,
    4. Piva M,
    5. Song C,
    6. Kong X,
    7. et al.
    Non-genomic and immune evolution of melanoma acquiring MAPKi resistance. Cell 2015;162:1271–85.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Bustelo XR
    . Vav family exchange factors: an integrated regulatory and functional view. Small GTPases 2014;5:9.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Gupta M,
    2. Qi X,
    3. Thakur V,
    4. Manor D
    . Tyrosine phosphorylation of Dbl regulates GTPase signaling. J Biol Chem 2014;289:17195–202.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. Fallahi-Sichani M,
    2. Becker V,
    3. Izar B,
    4. Baker GJ,
    5. Lin JR,
    6. Boswell SA,
    7. et al.
    Adaptive resistance of melanoma cells to RAF inhibition via reversible induction of a slowly dividing de-differentiated state. Mol Syst Biol 2017;13:905.
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    1. Girotti MR,
    2. Pedersen M,
    3. Sanchez-Laorden B,
    4. Viros A,
    5. Turajlic S,
    6. Niculescu-Duvaz D,
    7. et al.
    Inhibiting EGF receptor or SRC family kinase signaling overcomes BRAF inhibitor resistance in melanoma. Cancer Discov 2013;3:158–67.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    1. Nam HJ,
    2. Im SA,
    3. Oh DY,
    4. Elvin P,
    5. Kim HP,
    6. Yoon YK,
    7. et al.
    Antitumor activity of saracatinib (AZD0530), a c-Src/Abl kinase inhibitor, alone or in combination with chemotherapeutic agents in gastric cancer. Mol Cancer Ther 2013;12:16–26.
    OpenUrlAbstract/FREE Full Text
  41. 41.↵
    1. Peng SB,
    2. Henry JR,
    3. Kaufman MD,
    4. Lu WP,
    5. Smith BD,
    6. Vogeti S,
    7. et al.
    Inhibition of RAF isoforms and active dimers by LY3009120 leads to anti-tumor activities in RAS or BRAF mutant cancers. Cancer Cell 2015;28:384–98.
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Yao Z,
    2. Torres NM,
    3. Tao A,
    4. Gao Y,
    5. Luo L,
    6. Li Q,
    7. et al.
    BRAF mutants evade ERK-dependent feedback by different mechanisms that determine their sensitivity to pharmacologic inhibition. Cancer Cell 2015;28:370–83.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Girotti MR,
    2. Lopes F,
    3. Preece N,
    4. Niculescu-Duvaz D,
    5. Zambon A,
    6. Davies L,
    7. et al.
    Paradox-breaking RAF inhibitors that also target SRC are effective in drug-resistant BRAF mutant melanoma. Cancer Cell 2015;27:85–96.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Vergani E,
    2. Vallacchi V,
    3. Frigerio S,
    4. Deho P,
    5. Mondellini P,
    6. Perego P,
    7. et al.
    Identification of MET and SRC activation in melanoma cell lines showing primary resistance to PLX4032. Neoplasia 2011;13:1132–42.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Yadav V,
    2. Zhang X,
    3. Liu J,
    4. Estrem S,
    5. Li S,
    6. Gong XQ,
    7. et al.
    Reactivation of mitogen-activated protein kinase (MAPK) pathway by FGF receptor 3 (FGFR3)/Ras mediates resistance to vemurafenib in human B-RAF V600E mutant melanoma. J Biol Chem 2012;287:28087–98.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Dummer R,
    2. Ascierto PA,
    3. Gogas HJ,
    4. Arance A,
    5. Mandala M,
    6. Liszkay G,
    7. et al.
    Encorafenib plus binimetinib versus vemurafenib or encorafenib in patients with BRAF-mutant melanoma (COLUMBUS): a multicentre, open-label, randomised phase 3 trial. Lancet Oncol 2018;19:603–15.
    OpenUrlCrossRef
View Abstract
PreviousNext
Back to top
Cancer Research: 79 (19)
October 2019
Volume 79, Issue 19
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Editorial Board (PDF)

Sign up for alerts

View this article with LENS

Open full page PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Cancer Research article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Src-Dependent DBL Family Members Drive Resistance to Vemurafenib in Human Melanoma
(Your Name) has forwarded a page to you from Cancer Research
(Your Name) thought you would be interested in this article in Cancer Research.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Src-Dependent DBL Family Members Drive Resistance to Vemurafenib in Human Melanoma
Charlotte R. Feddersen, Jacob L. Schillo, Afshin Varzavand, Hayley R. Vaughn, Lexy S. Wadsworth, Andrew P. Voigt, Eliot Y. Zhu, Brooke M. Jennings, Sarah A. Mullen, Jeremy Bobera, Jesse D. Riordan, Christopher S. Stipp and Adam J. Dupuy
Cancer Res October 1 2019 (79) (19) 5074-5087; DOI: 10.1158/0008-5472.CAN-19-0244

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Src-Dependent DBL Family Members Drive Resistance to Vemurafenib in Human Melanoma
Charlotte R. Feddersen, Jacob L. Schillo, Afshin Varzavand, Hayley R. Vaughn, Lexy S. Wadsworth, Andrew P. Voigt, Eliot Y. Zhu, Brooke M. Jennings, Sarah A. Mullen, Jeremy Bobera, Jesse D. Riordan, Christopher S. Stipp and Adam J. Dupuy
Cancer Res October 1 2019 (79) (19) 5074-5087; DOI: 10.1158/0008-5472.CAN-19-0244
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Disclosure of Potential Conflicts of Interest
    • Authors' Contributions
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF
Advertisement

Related Articles

Cited By...

More in this TOC Section

  • ELF3 Regulates Cell Motility and Tumor Immunity
  • Unbiased IGF1R Downregulation through GRK2 Inhibition
  • Semaphorin 7a in ER+ Breast Cancer
Show more Translational Science
  • Home
  • Alerts
  • Feedback
  • Privacy Policy
Facebook  Twitter  LinkedIn  YouTube  RSS

Articles

  • Online First
  • Current Issue
  • Past Issues
  • Meeting Abstracts

Info for

  • Authors
  • Subscribers
  • Advertisers
  • Librarians

About Cancer Research

  • About the Journal
  • Editorial Board
  • Permissions
  • Submit a Manuscript
AACR logo

Copyright © 2021 by the American Association for Cancer Research.

Cancer Research Online ISSN: 1538-7445
Cancer Research Print ISSN: 0008-5472
Journal of Cancer Research ISSN: 0099-7013
American Journal of Cancer ISSN: 0099-7374

Advertisement