Next Article in Journal
Study: International Multicentric Minimally Invasive Liver Resection for Colorectal Liver Metastases (SIMMILR-CRLM)
Next Article in Special Issue
Bone Progenitors Pull the Strings on the Early Metabolic Rewiring Occurring in Prostate Cancer Cells
Previous Article in Journal
Assessing the Acceptability of Home-Based HPV Self-Sampling: A Qualitative Study on Cervical Cancer Screening Conducted in Reunion Island Prior to the RESISTE Trial
Previous Article in Special Issue
Dissecting Tumor Growth: The Role of Cancer Stem Cells in Drug Resistance and Recurrence
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MET Exon 14 Splice-Site Mutations Preferentially Activate KRAS Signaling to Drive Tumourigenesis

1
Department of Integrative Oncology, BC Cancer Research Institute, Vancouver, BC V5Z 1L3, Canada
2
Department of Interdisciplinary Oncology, University of British Columbia, Vancouver, BC V5Z 1L3, Canada
3
Department of Pathology & Laboratory Medicine, University of British Columbia, Vancouver, BC V6T 1Z7, Canada
4
Department of Pathology, Memorial Sloan-Kettering Cancer Center, New York, NY 10065, USA
5
Human Oncology and Pathogenesis Program, Memorial Sloan-Kettering Cancer Center, New York, NY 10065, USA
*
Author to whom correspondence should be addressed.
Cancers 2022, 14(6), 1378; https://doi.org/10.3390/cancers14061378
Submission received: 21 December 2021 / Revised: 17 February 2022 / Accepted: 23 February 2022 / Published: 8 March 2022

Abstract

:

Simple Summary

MET exon 14 splice-site mutations occur in ~3–4% of lung adenocarcinoma cases, defining a cohort of patients which might benefit from anti-MET targeted therapy. Such therapies have yielded mixed results, however, pointing to the need for better treatment design. Our study sought to aid this by characterizing key changes in mutant MET signaling behaviour. We first compared the transcriptional profiles of lung tumours with either METΔex14 or wild-type MET-amplification. METΔex14-mutant tumours exhibited increased activation of the Ras-MAPK pathway, consistent with our observations in an isogenic model system. Furthermore, sustained activity of this pathway is necessary for proliferation and maintenance of METΔex14 tumours, while forced reactivation of this pathway is sufficient to restore growth in the absence of MET activity. Our findings suggest that the MAPK pathway represents a main effector of METΔex14-driven cancer, lending credence to the possibility of combined MET-MAPK inhibition to improve therapeutic outcomes.

Abstract

Targeted therapies for MET exon 14-skipping (METΔex14)-driven lung cancers have generated some promising results but response rates remain below that seen for other kinase-driven cancers. One strategy for improving treatment outcomes is to employ rational combination therapies to enhance the suppression of tumour growth and delay or prevent the emergence of resistance. To this end, we profiled the transcriptomes of MET-addicted lung tumours and cell lines and identified the RAS-mitogen-activated protein kinase (MAPK) pathway as a critical effector required for METΔex14-dependent growth. Ectopic expression of MET in an isogenic cell line model showed that overexpression of the mutant MET receptor led to higher levels of MAPK phosphorylation and nuclear import, resulting in increased expression and phosphorylation of nuclear MAPK targets. In comparison, other known MET effectors were unaffected. Inhibition of this pathway by KRAS knockdown in MET-addicted cells in vitro led to decreased viability in only the METΔex14-mutant cells. Conversely, decoupling RAS-MAPK axis, but not other effector pathways, from MET activity via the introduction of constitutively active mutants conferred resistance to MET inhibitors in vitro. Our results suggest that aberrant hyperactivity of the MET receptor caused by the exon 14-skipping mutation does not uniformly upregulate all known downstream effectors, rather gaining a predilection for aberrantly activating and subsequently relying on the RAS-MAPK pathway. These findings provide a rationale for the co-targeting of the RAS-MAPK pathway alongside MET to prolong therapeutic response and circumvent resistance to improve patient survival.

1. Introduction

The advent of improved sequencing technologies over the past decade has led to routine detection of actionable alterations becoming standard-of-care in patients with non-small cell lung cancer (NSCLC). Tumours harbouring mutations in genes such as EGFR [1,2,3] and ERBB2 [4,5], or gene rearrangements involving ALK [6,7], RET [8], or ROS1 [9] are typically dependent on the oncogene for survival, with profound anti-tumour effects observed in response to the appropriate targeted therapy. Frequently found to be overexpressed or amplified in lung cancers, aberrant activation of the Hepatocyte Growth Factor (MET) receptor typically co-occurs with other established oncogenic targets [10]. However, recent findings demonstrating that activating mutations in the MET gene occur in mutual exclusivity with mutations in other oncogenes have led to intensified interest in its potential as a therapeutic target [11,12].
As part of the receptor tyrosine kinase (RTK) family, MET plays a highly versatile role in the cell by functioning to integrate growth, survival, and signaling cues from the extracellular environment. This role necessitates the use of diverse signaling pathways, tightly regulated via dedicated signal transducers, adaptors, and scaffolding proteins to collectively modulate receptor activity and signal transduction [13]. Importantly, somatic mutations in the MET gene leading to loss of exon 14 (METΔex14) is a recurrent event in ~3–4% of NSCLC cases [12]. These mutations result in loss of the juxtamembrane (JM) regulatory domain of the receptor, impairing the cell’s ability to ubiquitylate and shut off MET receptor activity, leading to oncogenesis [14].
MET has long been known to be a critical mediator in processes ranging from foetal organogenesis during development, to tissue regeneration and inflammation reduction during wound healing [15,16,17]. Its oncogenic potential was also noted decades ago with the identification of TPR-MET fusion proteins in chemically transformed carcinoma cell lines [18]. However, it was the recent NGS-based discovery of recurrent splice-site mutations in MET, leading to oncogenic loss of exon14 in LUAD patients, which has generated interest in its potential as a candidate for targeted therapy [12,19]. MET-targeted therapies using repurposed inhibitors have met with some initial success in the clinic; however, unlike cancers driven by other oncogenes such as EGFR, response rates in METΔex14 mutant tumours are lower [20,21,22]. In addition, as in other instances of targeted therapy, long-term patient survival remains poor due to acquired drug resistance [23,24,25]. A viable strategy for enhancing initial and long-term response is to design combinatorial therapies that target both the primary oncoprotein and its key effector(s) as first-line treatment. This strategy of “vertical inhibition” has been deployed with some success in BRAFV600E/K-positive melanoma patients, resulting in some improvements to overall response and progression-free survival [26,27,28]. In order to design similar therapies targeting RTK-driven cancers, identifying the crucial effector pathways responsible for maintaining oncogene addiction becomes a necessary pre-requisite. Given the relatively novel discovery of MET exon 14-skipping mutations in LUAD, there is a lack of insight towards how the METΔex14 oncogene aberrantly engages its parallel downstream signaling pathways, especially compared to its role when amplified in other oncogene-driven cancers. In this study, we show that MET exon 14-skipping mutations specifically enhance signaling through the RAS/MAPK pathway, providing a potential avenue of combinatorial therapy to enhance patient survival.

2. Materials and Methods

2.1. Cell Lines and Culture Conditions

Cell lines were obtained from American Type Culture Collection (ATCC, Manassas, VA, USA) and cultured according to ATCC guidelines. Hs746T, H1993, and HEK293T cells were maintained in DMEM (Gibco, cat #12430-062, Waltham, MA, USA) supplemented with 10% (v/v) foetal bovine serum (FBS) (Gibco, cat #26140079) and 1% (v/v) Glutamax (Gibco, cat #250-30081). H596 and HPL1D cells were maintained in RPMI-1640 (Gibco, cat# 11875119) supplemented with 10% FBS and 1% Glutamax. Cultured cell lines were incubated at 37 °C in a humidified incubator infused with 5% CO2. All lines were maintained in 10 cm tissue culture dishes (Sarstedt, cat #83-3902, Nümbrecht, Germany) and propagated when approaching 75–85% confluency using 0.25% Trypsin-EDTA (Gibco, cat #25200114). All cell lines were tested regularly for mycoplasma by polymerase chain reaction as previously described [29]. All cell lines were authenticated by the supplier and once received were subcultured and stocks frozen. A new vial was thawed and used every 20–25 passages.

2.2. Generation of Plasmids, Lentiviruses, and Viral Transduction to Create Stable Expression Cells

pDONR223 PIK3CA_p.E545K (plasmid #82881), pDONR223 STAT3_p.A662C_N664C (plasmid #82184), and pDONR223 KRAS_p.G12D (plasmid #81651) were ordered from Addgene. Human MET_WT, MET_ΔEx14, and MET_p.Y1003F ORFs were subcloned into pCR8/GW/TOPO (Thermo Fisher Scientific, cat #K250020, Waltham, MA, USA) as we previously described [24]. Target genes were transferred by Gateway LR Clonase II enzyme mix (Thermo Scientific, cat #11791-020, Waltham, MA, USA) to pInducer20 (gift from S. Elledge, Harvard, Cambridge, MA, USA). Presence of insert was verified using PCR and restriction enzyme mapping. Lentivirus was generated by co-transfection of psPAX2 (Addgene, plasmid #12260, Watertown, MA, USA) and pMD2.G (Addgene, plasmid #12259) into HEK293T cells (ATCC) using Lipofectamine2000 (Thermo Scientific, cat #11668019). Lentiviral media containing 8 µg/mL hexadimethrine bromide/polybrene (MilliporeSigma, cat #H9268-10G, Burlington, MA, USA) was applied to target cells in Falcon 6-well plates (VWR, cat#353046, Radnor, PA, USA), followed by spinfection at 1200 G for 60 min. Target cells were maintained polyclonally in 800 µg/mL G-418 Sulphate (Gibco, cat #10131027) 48–72 h post-spinfection.

2.3. Doxycycline-Inducible shRNA and Overexpression Experiments

pTRIPZ:shKRAS and pTRIPZ:shScrambled plasmids were gifts from S. Dedhar, BC Cancer. Target cells virally transduced with a doxycycline-inducible expression system underwent antibiotic selection in 4 µg/mL puromycin and were subsequently maintained in their respective media containing 10% (v/v) Tet-system approved FBS (Takara Bio, cat #631101, Kyoto, Japan). Cells were maintained in 100–500 ng/mL doxycycline (Takara Bio, cat #631311) for 24–72 h prior to plating or harvest. Doxycycline-treated cells seeded for further assays were maintained in doxycycline until endpoint.

2.4. Growth Assays

Dose-response experiments were performed in 96-well flat-bottom plates (Sarstedt, cat #83-3924). To establish dose-response curves to various drugs, doxycycline-treated cells were seeded in quadruplicate at a density of 2000 cells per well and treated with the desired concentrations of chemicals. Viability was assayed 96 h following treatment using alamarBlue viability dye (Thermo Fisher Scientific, cat #A50100). Fluorescence of the dye (Ex: 560 nM, Em: 590 nM) following 2–4 h incubation was measured using a Cytation3 microplate reader (BioTek Instruments, Winooski, VT, USA). Curve-fitting of absorbance readings by non-linear regression with fitting by least sum of squares was performed using Graphpad Prism v8.4.3 to obtain IC50 values. Parameters calculated for each cell line were compared to control parameters by Extra sum-of-squares F-test. p-values < 0.05 were considered statistically significant. Long-term growth assays were performed in six-well flat bottom plates (Sarstedt, cat #83-392). Doxycycline-treated cells were seeded in triplicate at a density of 10,000–20,000 cells per well, based on growth rate, in media containing the desired concentrations of chemicals. Media was refreshed twice weekly for 14 days, following which cells were fixed with 100% methanol (Fisher Scientific, cat #A433P-4, Pittsburgh, PA, USA) and stained with 1:1 DPBS: Crystal Violet solution (MilliporeSigma, cat #HT90132-1L, Burlington, MA, USA). Cell growth was assessed using alamarBlue and by destaining CV-stained plates using 10% acetic acid and transferring the resulting solution into 96-well plates. The absorbance readings (Ex: 560 nM, Em: 590 nM) were measured to compare relative cell growth between conditions.

2.5. RNA Extraction and Microarray Analysis

H1993, Hs746T, and H596 cells were seeded in six-well plates (150,000 cells/well) and incubated overnight. Each cell line, in triplicate, was subjected to 6 h incubation under three different conditions: 0.1% DMSO (baseline), 0.1% DMSO + 25 ng/mL HGF (stimulatory), or 10 µM Cabozantinib + 25 ng/mL HGF (inhibitory). Cells were lysed in 300 µL RNA lysis buffer at room temperature and RNA purified using the Zymo Research Quick-RNA Miniprep kit (Cedarlane, cat #R1055, Burlington, ON, Canada). The RNA extraction, DNase I treatment, and column purification protocol was followed according to the manufacturer’s instructions. RNA concentration was measured using the Qubit 2.0 RNA HS Assay kit (Thermo Scientific, cat #32852). Samples were sent to The Centre for Applied Genomics (TCAG, Toronto, Canada), where RNA integrity, sample labelling, array hybridisation, and data acquisition steps were performed, using the GeneChip Human Gene 2.0 ST Array (Thermo Fisher Scientific, cat #902113, Waltham, MA, USA). Cross-array normalisation and background correction of raw expression values was performed using Robust Multiarray Analysis [30] to obtain triplicate expression values (in log2 scale) for each oligonucleotide probe. Annotation of gene names and symbols to each probeset was done on R using the “hugene20sttranscriptcluster.db” package available through Bioconductor [31]. Expression values for unmapped probes were dropped, while the average expression value for probes corresponding to the same gene was taken for each sample. Gene Set Enrichment Analysis (GSEA) [32,33] was used together with the “Oncogenic Signatures” library of curated gene sets from MSigDB, comparing gene expression profiles between MET_WT (H1993) and MET_ΔEx14 (Hs746T and H596) cells in the stimulatory state to identify top enriched gene sets. Gene Set Variation Analysis (GSVA) [34] was used to compare relative enrichment of KRAS-associated gene sets between cell lines and between stimulatory, inhibitory, and baseline conditions. Leading edge analysis was performed as previously described [33] to identify leading edge genes, and the R package ComplexHeatmap [35] was used to make all the heat maps. Individual differentially expressed (DE) genes were identified using the limma R [36] package available through Bioconductor. Benjamini–Hochberg-adjusted p-value > 0.05 and log2 (fold change) > 0 cut-offs were used to identify positive hits for all analyses.

2.6. Tumour RNA-Seq Expression Analyses

RNA-seq data and accompanying clinicopathological information (Supplementary Table S1) for 585 lung adenocarcinomas, of which 230 were profiled by The Cancer Genome Atlas [12], were downloaded from the NIH GDC Data Portal. Gene expression values were log2(x + 1)-normalised prior to analysis. Using the MSKCC cBioPortal, six tumours with wild-type MET amplification and eight tumours that harboured MET exon 14-skipping mutations were identified for expression analysis. GSVA was used to compare the relative enrichment of KRAS signature gene sets between tumours from the two groups.

2.7. Gene Ontology (GO) Analysis

GO [37] was used to determine the function of leading edge upregulated genes selected from gene sets associated with KRAS-driven cancers. The overrepresentation test of Protein Analysis Through Evolutionary Relationships (PANTHER; version 16.0) with Fisher’s Exact Test (FDR < 0.05 multiple test correction) was used to functionally classify the genes of interest, employing the PANTHER GO Biological Process annotation set [38,39]. Redundant GO terms were trimmed using REVIGO [40], and the remaining terms were visualised on R alongside fold change enrichment and FDR values as previously described [41].

2.8. Ras Activity Assay

The Active Ras Detection Kit (CST, cat #8821) was used to assess relative fractions of GTP-bound Ras in cells. The protocol for cell lysis and GST-RBD pulldown was followed according to the manufacturer’s instructions.

2.9. Protein Extraction and Western Blot Analysis

Cells were grown in 6-well, 6-cm, or 10-cm plates to 70–80% confluency. Prior to harvest, cells were washed with ice-cold PBS (Gibco, cat #14190250) and lysed with an appropriate volume of RIPA lysis buffer (G-Biosciences, cat #CA95029-284, St. Louis, MO, USA) containing Halt Protease and Phosphatase Inhibitor Mixture (Thermo Scientific, cat #PI78446). Cells were harvested by scraping and were frozen overnight at −80 °C. All samples were sonicated, then cleared of cell debris via centrifugation at 15,000× g for 15 min at 4 °C. Protein concentration was measured using the Pierce BCA Protein Assay kit (Thermo Scientific, cat #PI23225) according to the manufacturer’s instructions. Absorbance (Em: 562 nM) was measured using a microplate reader and sample protein concentration calculated using a set of BSA standards. In total, 30 µg of protein from each sample was mixed with NuPAGE LDS Sample Buffer (Thermo Scientific, cat #NP0008) and NuPAGE Sample Reducing Agent (Thermo Scientific, cat #NP0009) and denatured by boiling at 95 °C for 5 min. Samples were loaded into NuPAGE Novex 4–12% Bis Tris Gels (Thermo Scientific) alongside Precision Plus Protein Kaleidoscope Standards (Bio-Rad Laboratories, cat #161-0375, Hercules, CA, USA) and electrophoresed at 100 V for 2 h, then transferred to PVDF Immobilon membranes (MilliporeSigma, cat #IPVH00010, Burlington, MA, USA) at 110 V for 70 min. Membranes were blocked at room temperature for 1 h in 5% (w/v) bovine serum albumin (Sigma, cat #A7906-100 G) dissolved in TBS-T [1x Tris-Buffered Saline (Bio-Rad Laboratories, cat #170-6435), 0.1% (v/v) Tween-20 (Thermo Fisher Scientific, cat #BP337-500)]. Blocked membranes were incubated with primary antibodies (listed below) overnight at 4 °C with shaking. Membranes were then washed 3× with TBS-T and incubated with the appropriate horseradish peroxidase-conjugated secondary antibody at 1:10,000 dilution for 1 h at room temperature. SuperSignal West Pico Plus Chemiluminescent Substrate (Thermo Scientific, cat #PI-34580) and SuperSignal West Femto Chemiluminescent Substrate (Thermo Scientific, cat #PI-34096) were used on the ChemiDocMP Gel Imaging System (Bio-Rad Laboratories) to visualise immunoblot bands. Band intensity was quantified by densitometry using ImageJ software and summarized in Supplementary Table S2. Uncropped western blot images are collated as Supplementary Figure S6.

2.10. Antibodies, Inhibitors, and Other Reagents

All antibodies used for Western blot analysis were diluted in 5% v/v BSA in TBS-T prior to use. Primary antibodies were used at 1:1000 dilution, while secondary antibodies were used at 1:10,000 dilution. All antibodies except anti-ERM/ETV5 (Abcam, cat #ab102010, Cambridge, United Kingdom) and anti-DUSP6 (Santa Cruz Biotechnology, cat #sc-37707, Dallas, TX, USA) were purchased from Cell Signaling Technology (Danvers, MA, USA): MET (8198S), p-MET Tyr1234 (3077S), p-MET Tyr1349 (3121S), AKT (4691S)), p-AKT Ser473 (4060L), p44/42 MAPK (137F5), p-p44/42 MAPK Thr202/Tyr204 (9101S), MEK1/2 (9122S), p-MEK1/2 Ser217/221 (9121S), STAT3 (12640S), p-STAT3 Y705 (9145S), SRC (2108S), p-SRC Y416 (59548S), Ras (8955S), cMYC (5605S), cFOS (2250S), p-cFOS S32 (5348S), FRA1 (5281S), p-FRA1 S265 (5841S), Vinculin (13901S), β-Actin HRP-linked (12620S), Histone H3 (4499S), Anti-Rabbit IgG HRP linked (7074S), and Anti-Mouse IgG HRP linked (7076S). MET kinase inhibitors crizotinib (PF-02341066), cabozantinib (BMS-907351), and capmatinib (INCB28060) were purchased from SelleckChem (Houston, TX, USA) and dissolved in DMSO (Thermo Fisher Scientific, cat #BP231-100) prior to use. Hepatocyte Growth Factor recombinant human protein (PHG0324) was purchased from Thermo Scientific and reconstituted according to manufacturer’s instruction prior to use.

2.11. Xenograft Studies

All animal experiments were conducted in accordance with protocols approved by the University of British Columbia Animal Care Committee (UBC ACC). Cell lines used for the generation of tumour xenografts were as follows: H1993 pTRIPZ:shKRAS, H1993 pTRIPZ:shScrambled, Hs746T pTRIPZ:shKRAS, and Hs746T pTRIPZ:shScrambled. Immediately prior to injection, all cells were trypsinised, counted, and resuspended in an appropriate volume of ice-cold PBS. Moreover, 2 × 106 pTRIPZ:shKRAS and pTRIPZ:shScrambled cells were injected subcutaneously into the left and right flanks of 8-week old male NRG mice (n = 16) in alternating fashion, with half the mice (n = 8) receiving H1993 cells and the other half receiving Hs746T cells. All mice were placed on doxycycline-infused feed (Envigo, cat #TD.130141, Indianapolis, IN, USA) starting 4 days prior to injection and continuing for the duration of the study. Measurement of tumour size was performed using callipers on a biweekly basis. Tumour volume was estimated using the equation V = 0.5 L·W2 and final weights were determined at endpoint.

3. Results

3.1. Cancer Cells Expressing METΔex14 Exhibit an Enhanced KRAS Activation Signature

Overexpression or low-level amplification of MET occurs in a significant fraction of LUADs, yet it represents a poor candidate for targeted monotherapy in these cases owing to its co-occurrence with alterations in other oncogenes [42]. In contrast, MET exon 14-skipping mutations, as well as high-level amplification of the wild type MET gene, constitute true driver alterations, which confer sensitivity to anti-MET therapy in some instances. However, while both receptor amplification and mutations leading to exon 14 skipping can result in aberrant hyperactivity of MET, they differ in the mechanism by which they stimulate receptor activity. There is evidence which suggests this may lead to differences in terms of receptor localisation, recycling, and protein interactions [43]. Thus, we wondered whether the downstream consequences on effector pathways might also differ depending on the route of hyperactivation, as this may impact the development of strategies for combination therapy (Figure 1A).
To determine whether such differences exist, we compared the status of known MET effector pathways in METΔex14-expressing cancer lines (Hs746T and H596) to another MET-addicted line with high-level amplification of the wild type MET gene (H1993). Notably, while all three lines exhibited a comparable degree of ERK1/2 (pERK1/2) and AKT (pAKT) phosphorylation at baseline, treatment with HGF to specifically stimulate MET activity resulted in a greater increase in phosphorylation of both RAS effectors in the METΔex14 lines compared to H1993 cells, suggesting that stimulation of METΔex14 may activate these pathways to a greater degree than the wild-type receptor (Figure 1B). In contrast, the phosphorylation status of the RAS-independent effectors STAT3 and SRC remained at comparable levels regardless of treatment condition (Figure 1B). Furthermore, treatment with the MET inhibitor cabozantinib blocked the HGF-induced increase of pERK1/2 and pAKT, with no effects on pSTAT3 or SRC in the METΔex14 cell lines. Our findings are consistent with recently published data. For example, phospho-kinase profiling MET-amplified (EBC-1) and METΔex14 (Hs746T) cells similarly reported higher ERK1/2 phosphorylation levels in Hs746T cells compared to EBC-1 [44]. These results together suggest that, while all three cell lines are reliant on MET for growth and/or survival, METΔex14-expressing cells exhibit higher levels of HGF-induced RAS pathway activity despite similar or lower levels of MET phosphorylation.
To further explore the potential impact of these observed differences on downstream signaling pathways and transcriptional programs, we performed microarray analysis of RNA collected from Hs746T, H596, and H1993 cells following treatment with HGF (MET agonist), DMSO (baseline), or HGF+cabozantinib (MET inhibitor). We compared expression profiles between each cell line and assessed the resulting differentially expressed genes using gene set enrichment analysis (GSEA) to identify transcriptional programs differentially regulated between METΔex14 and MET-amplified lines under conditions of MET stimulation, with a focus on the “Oncogenic Signatures” library of gene sets representing transcriptional signatures of cellular pathways known to be dysregulated in cancer (Figure 1C). This revealed that in comparison to MET-amplified cells, METΔex14 cell lines specifically demonstrated enrichment of multiple KRAS-mediated gene sets corresponding to genes found upregulated in epithelial cancer cells expressing oncogenic KRAS (Figure 1D, bottom right panel: p-value = 0.033; Fisher’s Exact Test). We repeated this analysis using the alternate “Hallmarks” library gene sets, and observed a similar result. These results suggest that, compared to cancers addicted to the wild type MET receptor, METΔex14-driven cells exhibit an expression signature more similar to KRAS-driven cancers. Given the lack of METΔex14 cell lines available, we expanded our analysis to include publicly available RNA-seq data [12] from tumours with either high-grade wildtype MET amplification (n = 6) or METΔex14 mutations (n = 8). In agreement with our cell line findings, GSEA results show a positive enrichment of KRAS-mediated gene sets in the transcriptomes of METΔex14-mutant tumours compared to those with high-grade wildtype MET amplification (Figure 1E). Gene Set Variation Analysis (GSVA), which provide single-sample enrichment scores (ES) instead of calculating ES based on aggregate sample group gene expression Z-scores, showed that a diverse collection of KRAS-mediated gene sets was positively enriched in most METΔex14-expressing lung tumour samples while the same sets were negatively enriched in all wildtype MET-amplified lung tumour samples by comparison (Figure 1F). Furthermore, GSEA using positional gene sets found most genes in the Chr12Q15 region to be overexpressed in METΔex14 tumour samples, suggesting focal amplification. Interestingly, this region encodes several p53 pathway inhibitors, including MDM2 and YEATS4 [45,46] (Supplementary Figure S1). This finding is consistent with previous reports of concurrent TP53 loss/MDM2 amplification in METΔex14 patient tumours [37].

3.2. METex14 Drives KRAS Signaling to a Greater Degree than the Wild-Type Receptor

While GSEA identified enrichment of gene sets associated with mutant KRAS-driven cancers in both METΔex14-expressing tumours and cell lines, it is unknown whether the observed KRAS activation signature is a direct consequence of mutant MET activity, or whether it is a MET-independent phenomenon associated as a potential co-driver with METΔex14-presenting tumours. To investigate this, we sought to determine how the relative enrichment of the top 5 KRAS-associated gene sets (Figure 1D) changed in response to HGF or cabozantinib treatment using GSVA. While these gene sets are significantly more enriched in METΔex14-expressing lines H596 and Hs746T, we observed that they are further upregulated in response to HGF treatment (H596) and downregulated in response to cabozantinib treatment (Hs746T) (Figure 2A). This observation is consistent with the knowledge that, while Hs746T cells are addicted to active MET signaling with or without HGF, H596 cells show MET-driven growth only in the presence of HGF due to an activating mutation in the PIK3CA gene. Comparison of transcript levels for leading edge genes across treatment conditions revealed a positive correlation between their relative expression levels to MET receptor activity in METΔex14 expressing cells. In contrast, we failed to observe any correlation of transcript levels to MET receptor activity in the wildtype MET-amplified cell line (Figure 2B). A parallel comparison using genes from the “HALLMARKS_KRAS_UP” set found that a larger fraction of these genes was significantly upregulated in METΔex14-expressing Hs746T (33/200 hits) and H596 (36/200 hits) cells in response to HGF stimulation, compared to the H1993 MET-amplified cell line (20/200) (Figure 2C). Some of these METΔex14-controlled genes include SERPINE1, MMP1, and NRP-1, known to promote metastatic and invasive behaviour, as well as promoters of growth and survival ELK3, ETV1, and TMEM158. Indeed, further GO analysis of the upregulated leading edge KRAS target genes revealed a critical involvement in functional processes associated with proliferative and invasive behaviour (e.g., regulation of cell migration and regulation of cell population proliferation) (Figure 2D). Interestingly, PLAU and PLAT were also upregulated in response to HGF treatment in Hs746T and H596 cells, respectively (Figure 2B); both genes code for plasminogen activators that are capable of cleaving environmental pro-HGF into its active form, suggesting a potential positive feedback mechanism promoting MET activity. Taken together, our data suggest that the activity of downstream KRAS pathway target genes, in particular those involved in regulating proliferative and metastatic behaviour, are regulated by MET receptor activity in splice mutant tumours in addition to being overexpressed relative to tumours with only MET amplification, a result that could be due to an apparent bias towards RAS-MAPK pathway activation in METΔex14-expressing cells.

3.3. METΔEx14-Expressing Cells Preferentially Hyperactivate the RAS-MAPK Pathway

Known downstream effectors of the MET receptor include the RAS-MAPK, PI3K-AKT, and JAK-STAT pathways (Figure 1A). On the basis of our transcriptomics results, we hypothesised that the RAS-MAPK pathway stands out as the most critical effector for mediating oncogenic METΔex14 activity. To investigate this, we established an isogenic MET overexpression model system to directly compare the effects of wild type and mutant MET signaling on its downstream effectors. Using a doxycycline-inducible vector, we generated HEK293T cell lines expressing wild type MET, its oncogenic mutant isoform METΔex14, or GFP. Additionally, we further included a MET-Y1003F mutant isoform to investigate whether loss of the c-Cbl binding site encoded by exon 14 phenocopied loss of the exon in its entirety. We found that, despite similar levels of total MET protein across all three established cell lines after the addition of doxycycline, the METΔex14-expressing cells exhibited a notable increase in phosphorylation of its kinase domain (Y1234). By comparison, phosphorylation of this residue in both the wild type and Y1003F mutant MET receptors was significantly lower (Figure 3A). Downstream of MET, we observed an increase in the level of active, GTP-bound RAS within the HEK293T cells, compared to the more modest increases observed in cells overexpressing wild type MET or the MET Y1003F mutant (Figure 3B). In agreement with this observation, phosphorylation levels of downstream RAS effectors MEK1/2 and ERK1/2 were higher, indicating increased activation of the RAS-MAPK pathway as a whole (Figure 3C). As MET is capable of recruiting and phosphorylating additional effectors, such as PI3K and STAT3, we assessed by Western blotting their relative degree of activating phosphorylation across our isogenic cell lines. We found that phosphorylation of these second messengers or their direct downstream targets did not vary between MET receptor status (Figure 3A), thereby suggesting that the mutant MET receptor is preferentially activating the RAS-MAPK pathway, rather than indiscriminately upregulating downstream signaling pathways due to increased receptor phosphorylation. In agreement with the findings from our isogenic system, we also observed greater ERK1/2 phosphorylation in MET-addicted cancer cells expressing METΔex14 (Hs746T/H596) vs. wildtype MET-amplified (H1993) cells (Figure 1A). Interestingly, ERK1/2 phosphorylation was largely lost in response to MET TKI treatment across both H1993 and Hs746T cells. However, only in H1993 cells were pERK1/2 levels able to rebound despite sustained MET inhibition (Supplementary Figure S2). This suggests that, while MET plays an active role in MAPK signaling in MET WT-addicted cells, it is not necessary for ERK1/2 phosphorylation. In contrast, Hs746T cells, which fail to show a rebound in pERK1/2 levels, appear to be wholly dependent on METex14 activity for MAPK signaling.
Though ERK1/2 initially undergoes phosphorylation via the MAPK cascade in the cytoplasm, its nuclear translocation is required for the induction of proliferation via phosphorylation of key transcription factors [47]. Thus, we sought to determine whether increased overall ERK1/2 phosphorylation via METΔex14 stimulation might lead to increased rates of nuclear importation, and subsequently, stronger growth effects. We first compared nuclear levels of pERK1/2 in our isogenic lines via fractionation of nuclear and cytoplasmic cell compartments with and without HGF treatment. Isogenic doxycycline-treated lung epithelial (HPL1D) cells overexpressing the MET WT or MET mutant receptor did not exhibit increased ERK phosphorylation levels overall. However, upon treatment with HGF, we observed a marked increase in both cytosolic and nuclear pERK1/2 levels, specifically in the METΔex14-expressing cells (Figure 3D, Supplementary Figure S5). Notably, METΔex14-expressing cells with greater levels of nuclear pERK1/2 also showed elevated levels of ETV5, cFOS, and FRA1 phosphorylation (Figure 3E). All three transcription factors are known nuclear targets of ERK1/2, and play key roles in promoting proliferation and maintenance of alveolar type II cells, which are the precursor for LUAD. These findings point to the possibility that the METΔex14 mutation enhances the ability of the MET receptor to drive growth through the RAS-MAPK pathway. Despite our above findings, however, overexpression of neither the wild type nor mutant receptors elicited an increase in proliferation in HEK293T cells (Supplementary Figure S3). However, a small but noticeable increase conferred by the METΔex14 receptor was observed in HPL1D lung epithelial lines under anchorage-independent conditions, an effect that was lost upon treatment with ERK1/2 inhibitor SCH772984 (Supplementary Figure S4).

3.4. MET-Driven Cell Growth Is Dependent on RAS/MAPK Pathway Activity

The propensity of the METΔex14 receptor to specifically hyperactivate RAS-MAPK signaling led us to evaluate whether METΔex14-addicted cancer cells are more reliant on this pathway for survival. To address this, we transfected pInducer20 constructs expressing KRAS-targeting shRNA (shKRAS I-01) or non-targeting control (shSCRAMBLE) into Hs746T and H1993 cells, enabling us to conditionally knock down KRAS in these cells upon the addition of doxycycline. Induction of shKRAS with doxycycline led to suppression of total RAS protein in H1993 cells (Figure 4A). However, multiple RAS proteins were detected in Hs746T cells, of which only one was successfully suppressed by shKRAS induction, suggesting additional members of the RAS gene family may play a role in these cancers besides KRAS. Nevertheless, shRNA-mediated KRAS knockdown led to complete and partial loss of MEK1/2 and ERK1/2 phosphorylation, respectively, in METΔex14-expressing Hs746T cells (Figure 4B). In contrast, KRAS knockdown in H1993 cells did not lead to a discernible decrease in MEK1/2 and ERK1/2 phosphorylation, suggesting that KRAS is not necessary for the maintenance of MAPK pathway activity in the context of wild type MET signaling. Interestingly, we found KRAS knockdown to have a significant detrimental effect on METΔex14-addicted Hs746T cells resulting in reduced cell proliferation and colony formation in vitro while growth of H1993 cells remained unaffected (Figure 4C). We attempted to validate our findings in vivo via subcutaneous injection of shKRAS- or shSCR-expressing H1993 and Hs746T cells. In contrast to our in vitro observations, we found that doxycycline-induced expression of shKRAS significantly impaired the growth of subcutaneous xenografts for both H1993 and Hs746T cells (Figure 4D–F), suggesting that the integrity of KRAS-mediated signaling activity remains a vital component of MET-driven oncogenesis.

3.5. MET-Independent KRAS Activation Is Sufficient to Rescue Cell Death following MET TKI Treatment in METΔEx14-Addicted Cells

Our results above suggested that RAS-MAPK signaling remained an important effector pathway for growth and survival in the context of MET-driven oncogenesis in general. However, recent clinical screens found that amplification and/or activating mutations in the gene encoding KRAS decoupled the RAS-MAPK pathway activity from the MET receptor, constituting the most common form of resistance against MET-targeted therapies in METΔex14-presenting tumours [23,24,25]. Thus, we wondered whether activation of this pathway was sufficient to account for the driving force behind the tumourigenic potential of METΔex14-addicted cells specifically. To investigate this hypothesis, we generated cancer lines conditionally expressing oncogenic forms of MET effectors KRAS, PI3K, and STAT3 (p20-KRASG12D, p20-PIK3CAE545A, and p20-STAT3A662C_N664C). In H596 lung adenosquamous cells, which exhibit HGF-dependent growth reversible via MET inhibition, we found that doxycycline-induced expression of constitutively active mutant KRAS led to HGF-independent cell growth (Figure 5H) as well as loss of sensitivity to MET inhibition (Figure 5B,E,G). Doxycycline-induced overexpression of a parallel MET effector, STAT3 (Figure 5A,D,G), by contrast failed to elicit any rescue, similar to GFP-expressing controls (Figure 5C,F,G). While these results point to KRAS as a key mediator of MET-driven growth, a key caveat remains in that these cells do not normally exhibit MET-addicted behaviour despite expressing an oncogenic form of MET receptor. As mentioned earlier, this is likely due to a coexisting PIK3CAE545A oncogenic mutant allele. Therefore, while constitutive KRAS activity is sufficient to decouple HGF-driven MET activity from cell growth in H596 cells, it is not possible to assess the importance of PI3K activity in facilitating this decoupling. To address this, we repeated our observations in MET-addicted H1993 and Hs746T lung adenocarcinoma and gastric carcinoma cells, respectively. Similar to our observation in H596 cells, overexpression of KRASG12D partially rescued MAPK phosphorylation (Figure 5M) and decoupled cell growth and survival from MET activity in METΔex14-expressing Hs746T cells (Figure 5J, bottom panel, Figure 5K). Notably, expression of the same mutant KRAS oncogene, despite similarly restoring MAPK phosphorylation from MET TKI treatment (Figure 5L), failed to rescue viability and growth in H1993 cells, which in contrast to Hs746T cells, express and are addicted to the wild type MET receptor (Figure 5J, top panel, Figure 5K). Finally, we observed that overexpression of oncogenic PIK3CAE545A, similar to STAT3A662C_N664C, did not rescue either line from MET TKI treatment (Figure 5I). As a whole, our results point to KRAS being the crucial signaling effector that is responsible for facilitating MET-driven growth in METΔexon14-addicted cells, while also playing a role in maintaining growth and tumourigenesis in MET WT-addicted cancer cells as well.

4. Discussion

The present study is aimed at facilitating the design of rational combinatorial therapies against METΔexon14-driven cancers, where response rates from single-treatment targeted therapies against MET alone remain largely mixed [21,48,49,50]. The response rates to capmatinib and tepotinib, two FDA-approved MET inhibitors, are below that seen with EGFR and RET inhibitors [51,52]. Here, we provide preclinical evidence pointing to the importance of KRAS as a critical mediator of METΔexon14-driven LUAD. Crucially, this observation was made by comparing the transcriptional profiles of MET-addicted cancers expressing either METΔexon14 or exhibiting high-grade amplification of the wild type MET gene, suggesting divergent mechanisms underpinning MET-driven tumourigenesis. This provides some justification for stratifying patients in the clinic based on MET alteration status when designing combined drug therapies, as one may expect differing mechanisms of primary or secondary resistance to first-line MET targeted therapy based on the type of alteration present. This is especially important given the low ORRs in patients with either alteration when treated with MET TKi [53]. Already, clinical trials testing the efficacy of MET inhibitors in METΔexon14-presenting tumours have identified amplification and/or activating mutations in the KRAS or BRAF genes (as well as activating mutations in the PI3K-AKT pathway) to be a common form of acquired resistance [23,24,25,54,55]. This observation is consistent with our findings that RAS is the primary mediator of METΔexon14-driven tumourigenesis, whereupon decoupling this effector from MET activity via mutant KRAS leads to MET TKi resistance in METΔexon14-mutant, but not MET-amplified, cells. By contrast, in vitro studies have identified alterations in MYC, mTOR, and PIK3CA, but notably rarely KRAS alterations, as potential drivers of resistance in MET-addicted cancer cells with MET-amp [56,57,58].
The ERK1/2 serine-threonine protein kinases sit at the end of a multi-tiered phosphorylation cascade, beginning with KRAS-mediated activation of RAF-1 (MAP3K) and ending with their phosphorylation and nuclear translocation [59]. This latter step is an important one, as it is required for ERK1/2 to promote ETS/AP1-mediated transcriptional regulation of cell proliferation and motility [47,59]. Interestingly, while ectopic expression of MET WT and METΔexon14 in human lung HPL1D cells led to similar levels of ERK1/2 phosphorylation (as opposed to higher ERK phosphorylation from ectopic METΔexon14 expression in HEK293T cells, which have relatively little background MAPK pathway activity), we nevertheless observed a markedly higher degree of pERK1/2 in the nuclear fractions of both cell lines. Consistent with these results, we also found higher levels of phosphorylated cFOS and FRA1, which together with FRA2 and c-JUN form the heterodimeric complex of AP1. Additionally, protein levels of ETV5, a key player in the maintenance of Ras-induced lung adenocarcinoma in alveolar type II cells [60], were higher in HEK293T cells expressing METΔexon14 compared to just ectopic MET WT overexpression. The oncogenic nature of MET exon14-skipping mutations would provide the most direct reason why METΔexon14 expression promotes greater increase in RAS-MAPK signaling compared to amplification of its wild type counterpart. However, high-grade amplification of wild type MET is also considered a driver alteration, such as in H1993 cells, and can be treated successfully in the clinic using anti-MET targeted therapies [61,62,63]. Furthermore, increased MET kinase domain phosphorylation in the mutant receptor appeared to lead to increased ERK1/2 phosphorylation without a concomitant increase in AKT or STAT3 phosphorylation. This suggests that, compared to the wild type receptor, the METΔexon14 mutant receptor has a specific predilection for activating the RAS-MAPK pathway.
One possible reason may be due to a critical importance of the juxtamembrane (JM) domain in regulating the balance between receptor recycling and degradation pathways. Rapid receptor internalisation following activation is a key mechanism by which cells rapidly switch off signaling to prevent sustained receptor activity, hence the oncogenic effects of losing this domain [64]. However, MET is known to continue signaling within its endosomal compartments following endocytosis, with some evidence suggesting that this compartmentalisation of internalised receptors plays an important role in spatially modulating signal transduction processes within the cell [65,66]. For instance, loss of the recycling adaptor GGA3 has been shown to enhance MET trafficking into degradative compartments, leading to attenuation of ERK1/2, but not AKT, phosphorylation following HGF stimulation [67,68]. Additionally, the authors showed that pERK1/2 attenuation as a result of defective MET recycling significantly impaired cell motility and migration, a process we found to be significantly enriched for in our GO analysis of GVSA leading edge genes upregulated in METΔexon14-expressing cell lines in response to HGF.
In addition to spatiotemporal regulation of MET signaling, its pleiotropic behaviour is tightly controlled through a plethora of binding partners, including adaptor proteins GRB2, CRK, and SHC, as well as the scaffolding protein GAB1 [13]. Site-specific recruitment of these proteins to the MET docking domain confers specificity to MET downstream signaling [69]. In particular, GRB2 binding is a necessary precursor for RAS activation through the recruitment of guanine nucleotide exchange factor SOS1 [70]. This MET-GRB2 interaction is an important mediator of MET-driven oncogenic behaviour, capable of predicting response to MET targeted therapy [71]. Changes to receptor structure that have the potential to introduce steric alterations may affect recruitment of binding partners, and of GRB2 especially if its reserved MET Y1356 binding site is affected [69]. X-ray crystallography experiments, similar to those performed on other MET mutants in the past, may reveal whether structural changes in the METΔexon14 receptor, if any, play a causative role in driving its oncogenic behaviour [72]. Finally, future studies incorporating the use of phosphoproteomics would help quantify the dynamics of METΔexon14 activation and downstream signaling networks; such techniques have been used to study the EGF receptor in the past, and can be easily deployed to characterise METΔexon14 signaling in the future [73,74].
To our knowledge, there have been no studies examining the direct molecular consequences of losing the JM domain in its entirety, though the MET Y1003 phosphotyrosine site within this domain has been identified to be functionally indispensable for mediating receptor degradation via its recruitment of the E3 ubiquitin ligase c-CBL [68]. On this basis, we were surprised to observe that introducing MET with a mutation at this site failed to phenocopy exon 14 loss in HEK293T cells in terms of downstream signaling and proliferation, despite the clear loss of MET ubiquitination [14]. However, Y1003 point mutations are comparatively rare in the clinic, suggesting that loss of the exon in its entirety confers additional oncogenic advantages beyond the impairment of CBL-induced MET degradation. Possibilities include the potential for structural impacts due to JM loss to affect MET effector recruitment, as well as involvement of other negative regulatory sites in this region, such as the Ser-985 phosphorylation site targeted by PKCδ and ε [75].

5. Conclusions

In summary, our report highlights distinct signaling and transcriptional properties elicited by different oncogenic driver alterations affecting MET. In the case of exon 14-skipping mutations, growth and survival advantages are conferred mainly through the RAS-MAPK pathway. This hallmark represents a potential point of weakness exploitable by vertical MET + MEK or MET + ERK inhibition, a strategy that was successfully deployed in EML4-ALK-positive LUADs. By comparison, high-grade amplification of wild-type MET appears to be reliant on activation of a more diverse complement of effector signaling pathways to drive tumourigenesis, potentially requiring a more personalised treatment approach. Additional studies comparing molecular mechanisms of MET hyperactivity will be required to fully understand the oncogenic nature of MET exon 14-skipping mutations and aid in the design of rational combinatorial therapies to combat drug resistance and improve treatment outcomes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14061378/s1, Figures S1–S5: Supplemental Figures; Table S1: TCGA clinicopathological information; Table S2: Densitometry analyses of all immunoblots; Figure S6: Original immunoblot images.

Author Contributions

Conceptualization, D.L., R.S. and W.W.L.; methodology, D.L. A.N., J.L.C., Y.T.C. and Q.M.; validation, D.L., A.N., J.L.C., Y.T.C. and Q.M.; formal analysis, D.L. and A.N.; investigation, D.L., A.N., J.L.C., Y.T.C. and Q.M.; resources, W.W.L.; data curation, D.L. and A.N.; writing—original draft preparation, D.L. and W.W.L.; writing—review and editing, D.L., R.S. and W.W.L.; visualization, D.L. and A.N.; supervision, W.W.L.; funding acquisition, W.W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC were funded by the Canadian Institutes of Health Research, grant number PJT-148725, and the BC Cancer Foundation.

Institutional Review Board Statement

All housing, care, and experimental conditions were followed in accordance with the guidelines of the University of British Columbia Animal Care Committee and the Canadian Council for Animal Care (protocol code A19-0090, approved 7/9/2021).

Data Availability Statement

Microarray data are uploaded as GEO: GSE198094.

Acknowledgments

We thank member of the Lockwood lab for critical reading of the manuscript.

Conflicts of Interest

R.S. has received research grants from LOXO Oncology, Merus, Helsinn Healthcare and Elevation Oncology Inc., all unrelated to the current study. W.W.L. is a consultant and has received research support from Hyperbio Therapeutics. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Paez, J.G.; Janne, P.A.; Lee, J.C.; Tracy, S.; Greulich, H.; Gabriel, S.; Herman, P.; Kaye, F.J.; Lindeman, N.; Boggon, T.J.; et al. EGFR Mutations in Lung Cancer: Correlation with Clinical Response to Gefitinib Therapy. Science 2004, 304, 1497–1500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Lynch, T.J.; Bell, D.W.; Sordella, R.; Gurubhagavatula, S.; Okimoto, R.A.; Brannigan, B.W.; Harris, P.L.; Haserlat, S.M.; Supko, J.G.; Haluska, F.G.; et al. Activating Mutations in the Epidermal Growth Factor Receptor Underlying Responsiveness of Non–Small-Cell Lung Cancer to Gefitinib. N. Engl. J. Med. 2004, 350, 2129–2139. [Google Scholar] [CrossRef] [PubMed]
  3. Shepherd, F.A.; Rodrigues Pereira, J.; Ciuleanu, T.; Tan, E.H.; Hirsh, V.; Thongprasert, S.; Campos, D.; Maoleekoonpiroj, S.; Smylie, M.; Martins, R.; et al. Erlotinib in Previously Treated Non–Small-Cell Lung Cancer. N. Engl. J. Med. 2005, 353, 123–132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Stephens, P.; Hunter, C.; Bignell, G.; Edkins, S.; Davies, H.; Teague, J.; Stevens, C.; O’Meara, S.; Smith, R.; Parker, A.; et al. Intragenic ERBB2 kinase mutations in tumours. Nature 2004, 431, 525–526. [Google Scholar] [CrossRef]
  5. Chuang, J.C.; Stehr, H.; Liang, Y.; Das, M.; Huang, J.; Diehn, M.; Wakelee, H.A.; Neal, J.W. ERBB2 -Mutated Metastatic Non–Small Cell Lung Cancer: Response and Resistance to Targeted Therapies. J. Thorac. Oncol. 2017, 12, 833–842. [Google Scholar] [CrossRef] [Green Version]
  6. Kwak, E.L.; Camidge, D.R.; Clark, J.; Shapiro, G.I.; Maki, R.G.; Ratain, M.J.; Solomon, B.; Bang, Y.; Ou, S.; Salgia, R. Clinical activity observed in a phase I dose escalation trial of an oral c-met and ALK inhibitor, PF-02341066. J. Clin. Oncol. 2009, 27, 3509. [Google Scholar] [CrossRef]
  7. Koivunen, J.P.; Mermel, C.; Zejnullahu, K.; Murphy, C.; Lifshits, E.; Holmes, A.J.; Choi, H.G.; Kim, J.; Chiang, D.; Thomas, R.; et al. EML4-ALK Fusion Gene and Efficacy of an ALK Kinase Inhibitor in Lung Cancer. Clin. Cancer Res. 2008, 14, 4275–4283. [Google Scholar] [CrossRef] [Green Version]
  8. Takeuchi, K.; Soda, M.; Togashi, Y.; Suzuki, R.; Sakata, S.; Hatano, S.; Asaka, R.; Hamanaka, W.; Ninomiya, H.; Uehara, H.; et al. RET, ROS1 and ALK fusions in lung cancer. Nat. Med. 2012, 18, 378–381. [Google Scholar] [CrossRef]
  9. Shaw, A.T.; Ou, S.H.I.; Bang, Y.J.; Camidge, D.R.; Solomon, B.J.; Salgia, R.; Riely, G.J.; Varella-Garcia, M.; Shapiro, G.I.; Costa, D.B.; et al. Crizotinib in ROS1-Rearranged Non-Small-Cell Lung Cancer. N. Engl. J. Med. 2014, 371, 1963–1971. [Google Scholar] [CrossRef] [Green Version]
  10. Christensen, J.G.; Burrows, J.; Salgia, R. c-Met as a target for human cancer and characterization of inhibitors for therapeutic intervention. Cancer Lett. 2005, 225, 1–26. [Google Scholar] [CrossRef]
  11. Heist, R.S.; Sequist, L.V.; Borger, D.; Gainor, J.F.; Arellano, R.S.; Le, L.P.; Dias-Santagata, D.; Clark, J.W.; Engelman, J.A.; Shaw, A.T.; et al. Acquired Resistance to Crizotinib in NSCLC with MET Exon 14 Skipping. J. Thorac. Oncol. 2016, 11, 1242–1245. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Cancer Genome Atlas Research Network. Comprehensive molecular profiling of lung adenocarcinoma. Nature 2014, 511, 543–550. [Google Scholar] [CrossRef] [PubMed]
  13. Trusolino, L.; Bertotti, A.; Comoglio, P.M. MET signaling: Principles and functions in development, organ regeneration and cancer. Nat. Rev. Mol. Cell Biol. 2010, 11, 834–848. [Google Scholar] [CrossRef] [PubMed]
  14. Kong-Beltran, M.; Seshagiri, S.; Zha, J.; Zhu, W.; Bhawe, K.; Mendoza, N.; Holcomb, T.; Pujara, K.; Stinson, J.; Fu, L.; et al. Somatic Mutations Lead to an Oncogenic Deletion of Met in Lung Cancer. Cancer Res. 2006, 66, 283–289. [Google Scholar] [CrossRef] [Green Version]
  15. Chmielowiec, J.; Borowiak, M.; Morkel, M.; Stradal, T.; Munz, B.; Werner, S.; Wehland, J.; Birchmeier, C.; Birchmeier, W. c-Met is essential for wound healing in the skin. J. Cell Biol. 2007, 177, 151–162. [Google Scholar] [CrossRef] [Green Version]
  16. Nakamura, T.; Sakai, K.; Nakamura, T.; Matsumoto, K. Hepatocyte growth factor twenty years on: Much more than a growth factor. J. Gastroenterol. Hepatol. 2011, 26, 188–202. [Google Scholar] [CrossRef] [Green Version]
  17. Bladt, F.; Riethmacher, D.; Isenmann, S.; Aguzzi, A.; Birchmeier, C. Essential role for the c-met receptor in the migration of myogenic precursor cells into the limb bud. Nature 1995, 376, 768–771. [Google Scholar] [CrossRef]
  18. Cooper, C.S.; Blair, D.G.; Oskarsson, M.K.; Tainsky, M.A.; Eader, L.A.; Vandewoude, G.F. Characterization of Human Transforming Genes from Chemically Transformed, Teratocarcinoma, and Pancreatic-Carcinoma Cell-Lines. Cancer Res. 1984, 44, 1–10. [Google Scholar]
  19. Seo, J.-S.; Ju, Y.S.; Lee, W.-C.; Shin, J.-Y.; Lee, J.K.; Bleazard, T.; Lee, J.; Jung, Y.J.; Kim, J.-O.; Shin, J.-Y.; et al. The transcriptional landscape and mutational profile of lung adenocarcinoma. Genome Res. 2012, 22, 2109–2119. [Google Scholar] [CrossRef] [Green Version]
  20. Paik, P.K.; Drilon, A.; Fan, P.-D.; Yu, H.; Rekhtman, N.; Ginsberg, M.S.; Borsu, L.; Schultz, N.; Berger, M.F.; Rudin, C.M.; et al. Response to MET Inhibitors in Patients with Stage IV Lung Adenocarcinomas Harboring MET Mutations Causing Exon 14 Skipping. Cancer Discov. 2015, 5, 842–849. [Google Scholar] [CrossRef] [Green Version]
  21. Paik, P.K.; Felip, E.; Veillon, R.; Sakai, H.; Cortot, A.B.; Garassino, M.C.; Mazieres, J.; Viteri, S.; Senellart, H.; Van Meerbeeck, J.; et al. Tepotinib in Non–Small-Cell Lung Cancer with MET Exon 14 Skipping Mutations. N. Engl. J. Med. 2020, 383, 931–943. [Google Scholar] [CrossRef] [PubMed]
  22. Drilon, A.E.; Camidge, D.R.; Ou, S.-H.I.; Clark, J.W.; Socinski, M.A.; Weiss, J.; Riely, G.J.; Winter, M.; Wang, S.C.; Monti, K.; et al. Efficacy and safety of crizotinib in patients (pts) with advanced MET exon 14-altered non-small cell lung cancer (NSCLC). J. Clin. Oncol. 2016, 34, 108. [Google Scholar] [CrossRef]
  23. Rotow, J.K.; Gui, P.; Wu, W.; Raymond, V.M.; Lanman, R.B.; Kaye, F.J.; Peled, N.; Fece de la Cruz, F.; Nadres, B.; Corcoran, R.B.; et al. Co-occurring Alterations in the RAS–MAPK Pathway Limit Response to MET Inhibitor Treatment in MET Exon 14 Skipping Mutation-Positive Lung Cancer. Clin. Cancer Res. 2019, 26, 439–449. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Suzawa, K.; Offin, M.; Lu, D.; Kurzatkowski, C.; Vojnic, M.; Smith, R.S.; Sabari, J.K.; Tai, H.; Mattar, M.; Khodos, I.; et al. Activation of KRAS Mediates Resistance to Targeted Therapy in MET Exon 14–mutant Non–small Cell Lung Cancer. Clin. Cancer Res. 2018, 25, 1248–1260. [Google Scholar] [CrossRef] [Green Version]
  25. Bahcall, M.; Awad, M.M.; Sholl, L.M.; Wilson, F.H.; Xu, M.; Wang, S.; Palakurthi, S.; Choi, J.; Ivanova, E.V.; Leonardi, G.C.; et al. Amplification of Wild-type KRAS Imparts Resistance to Crizotinib in MET Exon 14 Mutant Non–Small Cell Lung Cancer. Clin. Cancer Res. 2018, 24, 5963–5976. [Google Scholar] [CrossRef] [Green Version]
  26. Robert, C.; Karaszewska, B.; Schachter, J.; Rutkowski, P.; Mackiewicz, A.; Stroiakovski, D.; Lichinitser, M.; Dummer, R.; Grange, F.; Mortier, L.; et al. Improved Overall Survival in Melanoma with Combined Dabrafenib and Trametinib. N. Engl. J. Med. 2015, 372, 30–39. [Google Scholar] [CrossRef] [Green Version]
  27. Long, G.V.; Stroyakovskiy, D.; Gogas, H.; Levchenko, E.; de Braud, F.; Larkin, J.; Garbe, C.; Jouary, T.; Hauschild, A.; Grob, J.-J.; et al. Dabrafenib and trametinib versus dabrafenib and placebo for Val600 BRAF-mutant melanoma: A multicentre, double-blind, phase 3 randomised controlled trial. Lancet 2015, 386, 444–451. [Google Scholar] [CrossRef]
  28. Ascierto, P.A.; McArthur, G.A.; Dréno, B.; Atkinson, V.; Liszkay, G.; Di Giacomo, A.M.; Mandalà, M.; Demidov, L.; Stroyakovskiy, D.; Thomas, L.; et al. Cobimetinib combined with vemurafenib in advanced BRAFV600-mutant melanoma (coBRIM): Updated efficacy results from a randomised, double-blind, phase 3 trial. Lancet Oncol. 2016, 17, 1248–1260. [Google Scholar] [CrossRef]
  29. Uphoff, C.C.; Drexler, H.G. Detecting mycoplasma contamination in cell cultures by polymerase chain reaction. Methods Mol. Biol. 2011, 731, 93–103. [Google Scholar] [CrossRef]
  30. Irizarry, R.A.; Hobbs, B.; Collin, F.; Beazer-Barclay, Y.D.; Antonellis, K.J.; Scherf, U.; Speed, T.P. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 2003, 4, 249–264. [Google Scholar] [CrossRef] [Green Version]
  31. MacDonald, J.W. hugene20sttranscriptcluster.db: Affymetrix hugene20 Annotation Data (Chip hugene20sttranscriptcluster). 2017. Available online: https://bioconductor.org/packages/release/data/annotation/html/hugene20sttranscriptcluster.db.html (accessed on 20 December 2021).
  32. Mootha, V.K.; Lindgren, C.M.; Eriksson, K.F.; Subramanian, A.; Sihag, S.; Lehar, J.; Puigserver, P.; Carlsson, E.; Ridderstrale, M.; Laurila, E.; et al. PGC-1 alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34, 267–273. [Google Scholar] [CrossRef] [PubMed]
  33. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Hanzelmann, S.; Castelo, R.; Guinney, J. GSVA: Gene set variation analysis for microarray and RNA-Seq data. Bmc Bioinform 2013, 14, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Gu, Z.E.; Eils, R.; Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 2016, 32, 3. [Google Scholar] [CrossRef] [Green Version]
  36. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.F.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43. [Google Scholar] [CrossRef] [PubMed]
  37. The Gene Ontology Consortium. Expansion of the Gene Ontology knowledgebase and resources. Nucleic Acids Res. 2017, 45, D331–D338. [Google Scholar] [CrossRef] [Green Version]
  38. Mi, H.Y.; Ebert, D.; Muruganujan, A.; Mills, C.; Albou, L.P.; Mushayamaha, T.; Thomas, P.D. PANTHER version 16: A revised family classification, tree-based classification tool, enhancer regions and extensive API. Nucleic Acids Res. 2021, 49, D394–D403. [Google Scholar] [CrossRef]
  39. Mi, H.; Muruganujan, A.; Huang, X.; Ebert, D.; Mills, C.; Guo, X.; Thomas, P.D. Protocol Update for large-scale genome and gene function analysis with the PANTHER classification system (v.14.0). Nat. Protoc. 2019, 14, 703–721. [Google Scholar] [CrossRef]
  40. Supek, F.; Bošnjak, M.; Škunca, N.; Šmuc, T. REVIGO Summarizes and Visualizes Long Lists of Gene Ontology Terms. PLoS ONE 2011, 6, e21800. [Google Scholar] [CrossRef] [Green Version]
  41. Bonnot, T.; Gillard, M.; Nagel, D. A Simple Protocol for Informative Visualization of Enriched Gene Ontology Terms. Bio-Protocol 2019, 9. [Google Scholar] [CrossRef]
  42. Drilon, A.; Cappuzzo, F.; Ignatius, S.H.; Camidge, D.R. Targeting MET in Lung Cancer: Will Expectations Finally Be MET? J. Thorac. Oncol. 2017, 12, 15–26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Barrow-McGee, R.; Kermorgant, S. Met endosomal signaling: In the right place, at the right time. Int. J. Biochem. Cell Biol. 2014, 49, 69–74. [Google Scholar] [CrossRef] [PubMed]
  44. Pudelko, L.; Jaehrling, F.; Reusch, C.; Viteri, S.; Stroh, C.; Linde, N.; Sanderson, M.P.; Musch, D.; Lebrun, C.J.; Keil, M.; et al. SHP2 Inhibition Influences Therapeutic Response to Tepotinib in Tumors with MET Alterations. Iscience 2020, 23. [Google Scholar] [CrossRef] [PubMed]
  45. Pikor, L.A.; Lockwood, W.W.; Thu, K.L.; Vucic, E.A.; Chari, R.; Gazdar, A.F.; Lam, S.; Lam, W.L. YEATS4 Is a Novel Oncogene Amplified in Non–Small Cell Lung Cancer That Regulates the p53 Pathway. Cancer Res. 2013, 73, 7301–7312. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Moll, U.M.; Petrenko, O. The MDM2-p53 interaction. Mol. Cancer Res. 2003, 1, 1001–1008. [Google Scholar] [PubMed]
  47. Plotnikov, A.; Flores, K.; Maik-Rachline, G.; Zehorai, E.; Kapri-Pardes, E.; Berti, D.A.; Hanoch, T.; Besser, M.J.; Seger, R. The nuclear translocation of ERK1/2 as an anticancer target. Nat. Commun. 2015, 6, 6685. [Google Scholar] [CrossRef]
  48. Drilon, A.; Clark, J.W.; Weiss, J.; Ou, S.-H.I.; Camidge, D.R.; Solomon, B.J.; Otterson, G.A.; Villaruz, L.C.; Riely, G.J.; Heist, R.S.; et al. Antitumor activity of crizotinib in lung cancers harboring a MET exon 14 alteration. Nat. Med. 2020, 26, 47–51. [Google Scholar] [CrossRef]
  49. Mazieres, J.; Paik, P.K.; Felip, E.; Veillon, R.; Sakai, H.; Cortot, A.B.; Viteri, S.; Garassino, M.C.; Van Meerbeeck, J.P.; Raskin, J.; et al. 1283P Tepotinib in patients (pts) with advanced NSCLC with MET exon 14 (METex14) skipping: Overall efficacy results from VISION cohort A. Ann. Oncol. 2020, 31, S828–S829. [Google Scholar] [CrossRef]
  50. Lu, S.; Fang, J.; Li, X.; Cao, L.; Zhou, J.; Guo, Q.; Liang, Z.; Cheng, Y.; Jiang, L.; Yang, N.; et al. Phase II study of savolitinib in patients (pts) with pulmonary sarcomatoid carcinoma (PSC) and other types of non-small cell lung cancer (NSCLC) harboring MET exon 14 skipping mutations (METex14+). J. Clin. Oncol. 2020, 38, 9519. [Google Scholar] [CrossRef]
  51. Soria, J.C.; Ohe, Y.; Vansteenkiste, J.; Reungwetwattana, T.; Chewaskulyong, B.; Lee, K.H.; Dechaphunkul, A.; Imamura, F.; Nogami, N.; Kurata, T.; et al. Osimertinib in Untreated EGFR-Mutated Advanced Non-Small-Cell Lung Cancer. N. Engl. J. Med. 2018, 378, 113–125. [Google Scholar] [CrossRef]
  52. Drilon, A.; Oxnard, G.R.; Tan, D.S.W.; Loong, H.H.F.; Johnson, M.; Gainor, J.; McCoach, C.E.; Gautschi, O.; Besse, B.; Cho, B.C.; et al. Efficacy of Selpercatinib in RET Fusion-Positive Non-Small-Cell Lung Cancer. N. Engl. J. Med. 2020, 383, 813–824. [Google Scholar] [CrossRef] [PubMed]
  53. Wolf, J.; Seto, T.; Han, J.-Y.; Reguart, N.; Garon, E.B.; Groen, H.J.M.; Tan, D.S.W.; Hida, T.; de Jonge, M.; Orlov, S.V.; et al. Capmatinib inMETExon 14–Mutated orMET-Amplified Non–Small-Cell Lung Cancer. N. Engl. J. Med. 2020, 383, 944–957. [Google Scholar] [CrossRef] [PubMed]
  54. Jamme, P.; Fernandes, M.; Copin, M.-C.; Descarpentries, C.; Escande, F.; Morabito, A.; Grégoire, V.; Jamme, M.; Baldacci, S.; Tulasne, D.; et al. Alterations in the PI3K Pathway Drive Resistance to MET Inhibitors in NSCLC Harboring MET Exon 14 Skipping Mutations. J. Thorac. Oncol. 2020, 15, 741–751. [Google Scholar] [CrossRef] [PubMed]
  55. Recondo, G.; Bahcall, M.; Spurr, L.F.; Che, J.; Ricciuti, B.; Leonardi, G.C.; Lo, Y.-C.; Li, Y.Y.; Lamberti, G.; Nguyen, T.; et al. Molecular Mechanisms of Acquired Resistance to MET Tyrosine Kinase Inhibitors in Patients with MET Exon 14–Mutant NSCLC. Clin. Cancer Res. 2020, 26, 2615–2625. [Google Scholar] [CrossRef] [PubMed]
  56. Kim, S.; Kim, T.M.; Kim, D.-W.; Kim, S.; Kim, M.; Ahn, Y.-O.; Keam, B.; Heo, D.S. Acquired Resistance of MET-Amplified Non-small Cell Lung Cancer Cells to the MET Inhibitor Capmatinib. Cancer Res. Treat 2019, 51, 951–962. [Google Scholar] [CrossRef] [PubMed]
  57. Shen, A.; Wang, L.; Huang, M.; Sun, J.; Chen, Y.; Shen, Y.-Y.; Yang, X.; Wang, X.; Ding, J.; Geng, M. c-Myc Alterations Confer Therapeutic Response and Acquired Resistance to c-Met Inhibitors in MET-Addicted Cancers. Cancer Res. 2015, 75, 4548–4559. [Google Scholar] [CrossRef] [Green Version]
  58. Henry, R.E.; Barry, E.R.; Castriotta, L.; Ladd, B.; Markovets, A.; Beran, G.; Ren, Y.; Zhou, F.; Adam, A.; Zinda, M.; et al. Acquired savolitinib resistance in non-small cell lung cancer arises via multiple mechanisms that converge on MET-independent mTOR and MYC activation. Oncotarget 2016, 7, 57651–57670. [Google Scholar] [CrossRef] [Green Version]
  59. Roskoski, R. ERK1/2 MAP kinases: Structure, function, and regulation. Pharmacol. Res. 2012, 66, 105–143. [Google Scholar] [CrossRef]
  60. Zhang, Z.; Newton, K.; Kummerfeld, S.K.; Webster, J.; Kirkpatrick, D.S.; Phu, L.; Eastham-Anderson, J.; Liu, J.; Lee, W.P.; Wu, J.; et al. Transcription factor Etv5 is essential for the maintenance of alveolar type II cells. Proc. Natl. Acad. Sci. USA 2017, 114, 3903–3908. [Google Scholar] [CrossRef] [Green Version]
  61. Lennerz, J.K.; Kwak, E.L.; Ackerman, A.; Michael, M.; Fox, S.B.; Bergethon, K.; Lauwers, G.Y.; Christensen, J.G.; Wilner, K.D.; Haber, D.A.; et al. MET Amplification Identifies a Small and Aggressive Subgroup of Esophagogastric Adenocarcinoma with Evidence of Responsiveness to Crizotinib. J. Clin. Oncol. 2011, 29, 4803–4810. [Google Scholar] [CrossRef] [Green Version]
  62. Camidge, D.R.; Ou, S.-H.I.; Shapiro, G.; Otterson, G.A.; Villaruz, L.C.; Villalona-Calero, M.A.; Iafrate, A.J.; Varella-Garcia, M.; Dacic, S.; Cardarella, S.; et al. Efficacy and safety of crizotinib in patients with advanced c-MET-amplified non-small cell lung cancer (NSCLC). J. Clin. Oncol. 2014, 32, 8001. [Google Scholar] [CrossRef]
  63. Song, Z.; Wang, H.; Yu, Z.; Lu, P.; Xu, C.; Chen, G.; Zhang, Y. De Novo MET Amplification in Chinese Patients with Non–Small-Cell Lung Cancer and Treatment Efficacy with Crizotinib: A Multicenter Retrospective Study. Clin. Lung Cancer 2019, 20, e171–e176. [Google Scholar] [CrossRef] [PubMed]
  64. Lai, A.Z.; Abella, J.V.; Park, M. Crosstalk in Met receptor oncogenesis. Trends Cell Biol. 2009, 19, 542–551. [Google Scholar] [CrossRef] [PubMed]
  65. Kermorgant, S.; Zicha, D.; Parker, P.J. PKC controls HGF-dependent c-Met traffic, signaling and cell migration. EMBO J. 2004, 23, 3721–3734. [Google Scholar] [CrossRef] [Green Version]
  66. Ménard, L.; Parker, P.J.; Kermorgant, S. Receptor tyrosine kinase c-Met controls the cytoskeleton from different endosomes via different pathways. Nat. Commun. 2014, 5, 3907. [Google Scholar] [CrossRef] [Green Version]
  67. Parachoniak, C.A.; Luo, Y.; Abella, J.V.; Keen, J.H.; Park, M. GGA3 Functions as a Switch to Promote Met Receptor Recycling, Essential for Sustained ERK and Cell Migration. Dev. Cell 2011, 20, 751–763. [Google Scholar] [CrossRef] [Green Version]
  68. Abella, J.V.; Peschard, P.; Naujokas, M.A.; Lin, T.; Saucier, C.; Urbé, S.; Park, M. Met/Hepatocyte Growth Factor Receptor Ubiquitination Suppresses Transformation and Is Required for Hrs Phosphorylation. Mol. Cell. Biol. 2005, 25, 9632–9645. [Google Scholar] [CrossRef] [Green Version]
  69. Schiering, N.; Knapp, S.; Marconi, M.; Flocco, M.M.; Cui, J.; Perego, R.; Rusconi, L.; Cristiani, C. Crystal structure of the tyrosine kinase domain of the hepatocyte growth factor receptor c-Met and its complex with the microbial alkaloid K-252a. Proc. Natl. Acad. Sci. USA 2003, 100, 12654–12659. [Google Scholar] [CrossRef] [Green Version]
  70. Giubellino, A.; Burke, T.R.; Bottaro, D.P. Grb2 signaling in cell motility and cancer. Expert Opin. Ther. Targets 2008, 12, 1021–1033. [Google Scholar] [CrossRef] [Green Version]
  71. Smith, M.A.; Licata, T.; Lakhani, A.; Garcia, M.V.; Schildhaus, H.-U.; Vuaroqueaux, V.; Halmos, B.; Borczuk, A.C.; Chen, Y.A.; Creelan, B.C.; et al. MET–GRB2 Signaling-Associated Complexes Correlate with Oncogenic MET Signaling and Sensitivity to MET Kinase Inhibitors. Clin. Cancer Res. 2017, 23, 7084–7096. [Google Scholar] [CrossRef] [Green Version]
  72. Collie, G.W.; Koh, C.M.; O’Neill, D.J.; Stubbs, C.J.; Khurana, P.; Eddershaw, A.; Snijder, A.; Mauritzson, F.; Barlind, L.; Dale, I.L.; et al. Structural and Molecular Insight into Resistance Mechanisms of First Generation cMET Inhibitors. ACS Med. Chem. Lett. 2019, 10, 1322–1327. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Chen, Y.; Low, T.-Y.; Choong, L.-Y.; Ray, R.S.; Tan, Y.-L.; Toy, W.; Lin, Q.; Ang, B.K.; Wong, C.H.; Lim, S.; et al. Phosphoproteomics identified Endofin, DCBLD2, and KIAA0582 as novel tyrosine phosphorylation targets of EGF signaling and Iressa in human cancer cells. Proteomics 2007, 7, 2384–2397. [Google Scholar] [CrossRef] [PubMed]
  74. Reddy, R.J.; Curran, T.G.; Zhang, Y.; White, F.M. Measurement of Phosphorylated Peptides with Absolute Quantification. Methods Mol. Biol. 2016, 281–292. [Google Scholar] [CrossRef]
  75. Hashigasako, A.; Machide, M.; Nakamura, T.; Matsumoto, K.; Nakamura, T. Bi-directional Regulation of Ser-985 Phosphorylation of c-Met via Protein Kinase C and Protein Phosphatase 2A Involves c-Met Activation and Cellular Responsiveness to Hepatocyte Growth Factor. J. Biol. Chem. 2004, 279, 26445–26452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. METΔex14 variant-expressing cancers exhibit an enhanced KRAS activation signature. (A) When activated by HGF, the MET receptor is capable of triggering multiple downstream pathways in parallel. This is because its multisubstrate docking domain is capable of recruiting a diverse set of kinases and signaling factors. (B) WB analysis WT MET-amp (H1993) and METΔex14-expressing (Hs746T, H596) cells comparing relative phosphorylation levels of MET effector signaling pathways Ras-MAPK, PI3K-Akt, Stat3, and Src at baseline (0.1% DMSO; 30 min), MET-stimulated (25 ng/mL HGF + 0.1% DMSO; 30 min), or MET-inhibited (25 ng/mL HGF + 10 μM Cabozantinib; 30 min) states. (C) Expression profiling workflow: to identify transcriptomic differences in the molecular profiles of cancer cells addicted to WT MET-amp vs. METΔex14 activity, we performed Gene Set Enrichment Analysis (GSEA) comparing the expression profiles of Hs746T and H596 cells to that of H1993 cells (figure created using BioRender.com). (D) GSEA revealed that gene set expression signatures typically associated with KRAS-driven cancers are among the top-most enriched in METΔex14-expressing cancer lines when compared against the wild type MET cancer line, while gene sets downregulated in KRAS-driven cancers were similarly negatively enriched in METΔex14-expressing cancer lines. As multiple KRAS-related gene sets exist in the “Oncogenic_signatures” MSigDB library curated by UC San Diego and Broad Institute, Fisher’s Exact Test was performed to show significant overrepresentation of KRAS-associated gene sets (bottom right panel). (E) GSEA comparing TCGA RNA-seq data of METΔex14-expressing tumours with WT MET-amp tumours revealed similar enrichment of gene sets associated with KRAS-driven cancers. (F) Gene Set Variation Analysis (GVSA) of TCGA RNA-seq data from (E) showing relative enrichment of KRAS-associated gene sets (compared to all other gene sets in the “Oncogenic_signatures library) for WT MET-amp and in METΔex14-expressing tumours, at single-sample resolution. Heatmap shows positive relative enrichment of KRAS-associated gene sets in METΔex14-expressing lung tumours, and negative relative enrichment in WT MET-amp tumours.
Figure 1. METΔex14 variant-expressing cancers exhibit an enhanced KRAS activation signature. (A) When activated by HGF, the MET receptor is capable of triggering multiple downstream pathways in parallel. This is because its multisubstrate docking domain is capable of recruiting a diverse set of kinases and signaling factors. (B) WB analysis WT MET-amp (H1993) and METΔex14-expressing (Hs746T, H596) cells comparing relative phosphorylation levels of MET effector signaling pathways Ras-MAPK, PI3K-Akt, Stat3, and Src at baseline (0.1% DMSO; 30 min), MET-stimulated (25 ng/mL HGF + 0.1% DMSO; 30 min), or MET-inhibited (25 ng/mL HGF + 10 μM Cabozantinib; 30 min) states. (C) Expression profiling workflow: to identify transcriptomic differences in the molecular profiles of cancer cells addicted to WT MET-amp vs. METΔex14 activity, we performed Gene Set Enrichment Analysis (GSEA) comparing the expression profiles of Hs746T and H596 cells to that of H1993 cells (figure created using BioRender.com). (D) GSEA revealed that gene set expression signatures typically associated with KRAS-driven cancers are among the top-most enriched in METΔex14-expressing cancer lines when compared against the wild type MET cancer line, while gene sets downregulated in KRAS-driven cancers were similarly negatively enriched in METΔex14-expressing cancer lines. As multiple KRAS-related gene sets exist in the “Oncogenic_signatures” MSigDB library curated by UC San Diego and Broad Institute, Fisher’s Exact Test was performed to show significant overrepresentation of KRAS-associated gene sets (bottom right panel). (E) GSEA comparing TCGA RNA-seq data of METΔex14-expressing tumours with WT MET-amp tumours revealed similar enrichment of gene sets associated with KRAS-driven cancers. (F) Gene Set Variation Analysis (GVSA) of TCGA RNA-seq data from (E) showing relative enrichment of KRAS-associated gene sets (compared to all other gene sets in the “Oncogenic_signatures library) for WT MET-amp and in METΔex14-expressing tumours, at single-sample resolution. Heatmap shows positive relative enrichment of KRAS-associated gene sets in METΔex14-expressing lung tumours, and negative relative enrichment in WT MET-amp tumours.
Cancers 14 01378 g001
Figure 2. KRAS activation signature observed in exon14-skipping mutants is dependent on MET receptor activity. (A) GSVA comparing relative enrichment of KRAS-associated gene sets to NES of all gene sets from the entire Oncogenic_signatures library. Heatmaps show positive enrichment of Ras-driven genetic signatures in METΔex14-expressing cancer lines, in contrast to negative enrichment observed for the wild type MET-amp driven cell line. (B) Unsupervised clustering identifies a subset of genes within Ras-associated gene sets which correlate with MET receptor activity in METΔex14-expressing cancer lines, but not in the H1993 cell line. (C) Volcano plots comparing relative changes in gene expression of H1993, Hs746T, and H596 cells between HGF stimulated or MET inhibited states. Significantly upregulated hits (log2[FC] > 0; Adj p-value < 0.05) that fell under the “HALLMARKS_KRAS_UP” gene set were tallied for each line, and the significance of their overrepresentation calculated using Fisher’s Exact Test. (D) Gene ontology analysis of leading edge genes shown in (B,C) heatmaps point to their functional involvement in angiogenesis, cell mobility, and proliferation.
Figure 2. KRAS activation signature observed in exon14-skipping mutants is dependent on MET receptor activity. (A) GSVA comparing relative enrichment of KRAS-associated gene sets to NES of all gene sets from the entire Oncogenic_signatures library. Heatmaps show positive enrichment of Ras-driven genetic signatures in METΔex14-expressing cancer lines, in contrast to negative enrichment observed for the wild type MET-amp driven cell line. (B) Unsupervised clustering identifies a subset of genes within Ras-associated gene sets which correlate with MET receptor activity in METΔex14-expressing cancer lines, but not in the H1993 cell line. (C) Volcano plots comparing relative changes in gene expression of H1993, Hs746T, and H596 cells between HGF stimulated or MET inhibited states. Significantly upregulated hits (log2[FC] > 0; Adj p-value < 0.05) that fell under the “HALLMARKS_KRAS_UP” gene set were tallied for each line, and the significance of their overrepresentation calculated using Fisher’s Exact Test. (D) Gene ontology analysis of leading edge genes shown in (B,C) heatmaps point to their functional involvement in angiogenesis, cell mobility, and proliferation.
Cancers 14 01378 g002
Figure 3. METΔex14-expressing cells preferentially hyperactivate the Ras-MAPK pathway. (A) Western blot analysis of major signaling effectors downstream of the MET receptor, including STAT3, Akt, and MAPK, show that all three pathways are upregulated in response to HGF treatment in MET WT or METΔex14-overexpressing cells, but that MAPK specifically exhibits higher phosphorylation levels in the METΔex14-expressing line. (B) GST-Raf-RBD pulldown demonstrates that HGF treatment (10 ng/mL) of METex14-expressing HEK293T cells results leads to a higher fraction of active, GTP-bound RAS compared to cells overexpressing wild type MET or GFP. (C) Western blot of Ras-MAPK pathway effectors in human embryonic kidney cells (HEK293T) ectopically expressing MET WT or METΔex14 under the dox-inducible TetO promoter. Cells expressing METΔex14 exhibit higher levels of MET receptor phosphorylation, as well as phosphorylation and activation of downstream Ras effectors MEK1/2 and ERK1/2. (D) Western blot analysis following separation of cytoplasmic and nuclear fractions show increased nuclear import of phosphorylated ERK1/2 in HGF-treated MET mutant-overexpressing cells, compared to MET WT or GFP-overexpressing cells. (E) Western blot analysis depicts METex14-expressing cells as exhibiting comparatively higher activating phosphorylation of select ERK1/2 nuclear targets known to promote growth, proliferation, or maintenance of alveolar type II cells.
Figure 3. METΔex14-expressing cells preferentially hyperactivate the Ras-MAPK pathway. (A) Western blot analysis of major signaling effectors downstream of the MET receptor, including STAT3, Akt, and MAPK, show that all three pathways are upregulated in response to HGF treatment in MET WT or METΔex14-overexpressing cells, but that MAPK specifically exhibits higher phosphorylation levels in the METΔex14-expressing line. (B) GST-Raf-RBD pulldown demonstrates that HGF treatment (10 ng/mL) of METex14-expressing HEK293T cells results leads to a higher fraction of active, GTP-bound RAS compared to cells overexpressing wild type MET or GFP. (C) Western blot of Ras-MAPK pathway effectors in human embryonic kidney cells (HEK293T) ectopically expressing MET WT or METΔex14 under the dox-inducible TetO promoter. Cells expressing METΔex14 exhibit higher levels of MET receptor phosphorylation, as well as phosphorylation and activation of downstream Ras effectors MEK1/2 and ERK1/2. (D) Western blot analysis following separation of cytoplasmic and nuclear fractions show increased nuclear import of phosphorylated ERK1/2 in HGF-treated MET mutant-overexpressing cells, compared to MET WT or GFP-overexpressing cells. (E) Western blot analysis depicts METex14-expressing cells as exhibiting comparatively higher activating phosphorylation of select ERK1/2 nuclear targets known to promote growth, proliferation, or maintenance of alveolar type II cells.
Cancers 14 01378 g003
Figure 4. MET-driven cell growth is dependent on activation of the Ras/MAPK pathway. (A) Western blot analysis depicting partial or total loss of Ras protein in H1993 and Hs746T carcinoma cells following doxycycline treatment (500 ng/mL Dox for 72 h). (B) Western blot analysis showing the effect of KRAS knockdown on the phosphorylation of downstream effectors MEK1/2 and ERK1/2. (C) Viability assay by 6-well clonogenics demonstrating a negative impact of KRAS loss on cell viability in METΔex14-expressing Hs746T cells compared to MET WT-expressing H1993 cells in vitro, as measured by Alamar Blue and Crystal Violet staining. (D) Tumour growth kinetics at endpoint of H1993 (shKRAS: n = 8; shSCR: n = 8) and Hs746T (shKRAS: n = 8; shSCR: n = 8) cells injected subcutaneously into flanks of NRG mice maintained on doxycycline diet (starting 4 days before injection) demonstrate that loss of KRAS hinders tumour growth in vivo in both MET-addicted cell lines. (E) Tumour weights at endpoint correlate with our results in (D), demonstrating that loss of KRAS hinders tumour growth in vivo in both cell lines. (F) Visual comparison of H1993 vs. Hs746T xenograft tumours with or without KRAS knockdown. n.s, not significant; * p < 0.05; ** p < 0.01 for the indicated comparisons.
Figure 4. MET-driven cell growth is dependent on activation of the Ras/MAPK pathway. (A) Western blot analysis depicting partial or total loss of Ras protein in H1993 and Hs746T carcinoma cells following doxycycline treatment (500 ng/mL Dox for 72 h). (B) Western blot analysis showing the effect of KRAS knockdown on the phosphorylation of downstream effectors MEK1/2 and ERK1/2. (C) Viability assay by 6-well clonogenics demonstrating a negative impact of KRAS loss on cell viability in METΔex14-expressing Hs746T cells compared to MET WT-expressing H1993 cells in vitro, as measured by Alamar Blue and Crystal Violet staining. (D) Tumour growth kinetics at endpoint of H1993 (shKRAS: n = 8; shSCR: n = 8) and Hs746T (shKRAS: n = 8; shSCR: n = 8) cells injected subcutaneously into flanks of NRG mice maintained on doxycycline diet (starting 4 days before injection) demonstrate that loss of KRAS hinders tumour growth in vivo in both MET-addicted cell lines. (E) Tumour weights at endpoint correlate with our results in (D), demonstrating that loss of KRAS hinders tumour growth in vivo in both cell lines. (F) Visual comparison of H1993 vs. Hs746T xenograft tumours with or without KRAS knockdown. n.s, not significant; * p < 0.05; ** p < 0.01 for the indicated comparisons.
Cancers 14 01378 g004
Figure 5. Decoupling MET from the Ras/MAPK pathway is sufficient to rescue cell death following MET TKI treatment in METΔEx14-addicted cells. (AF) Dose response curves of H596 lung adenosquamous cells to Crizotinib and Capmatinib, depicting loss of HGF-dependent sensitivity to MET tyrosine kinase inhibitors in response to KRASG12D expression. (G) Cell viability assays of H596 cells conditionally expressing oncogenic STAT3, KRAS, or GFP in the presence (top) or absence (bottom) of HGF. Values shown are calculated from Alamar Blue absorbance readings of cells grown in MET TKI (Capmatinib, 100 nM) relative to those grown in 0.1% DMSO. (H) Growth assays of H596 cells cultured under indicated conditions to compare the impact of oncogenic STAT3, KRAS, or GFP expression on H596 cell growth in the presence or absence of HGF-driven proliferation. (I) Western blot analysis of MET effector proteins overexpressed by the pInducer20 promoter in the presence of doxycycline (200 ng/mL), along with their downstream signaling targets with or without HGF (10 ng/mL) or Capmatinib (100 nM) in H596 cells. (J,K) Viability assay in 6-well plates, as measured by (J) Alamar Blue and (K) Crystal Violet staining, demonstrating mutant KRAS-mediated rescue cell viability from MET TKI treatment in Hs746T, but not H1993, cells. (L,M) Western blot analysis of MET effector proteins overexpressed by the pInducer20 promoter in the presence of doxycycline (200 ng/mL), along with their downstream signaling targets with or without Capmatinib (100 ng/mL), in (L) H1993 and (M) Hs746T cells. n.s, not significant; * p < 0.05; *** p < 0.001 for the indicated comparisons.
Figure 5. Decoupling MET from the Ras/MAPK pathway is sufficient to rescue cell death following MET TKI treatment in METΔEx14-addicted cells. (AF) Dose response curves of H596 lung adenosquamous cells to Crizotinib and Capmatinib, depicting loss of HGF-dependent sensitivity to MET tyrosine kinase inhibitors in response to KRASG12D expression. (G) Cell viability assays of H596 cells conditionally expressing oncogenic STAT3, KRAS, or GFP in the presence (top) or absence (bottom) of HGF. Values shown are calculated from Alamar Blue absorbance readings of cells grown in MET TKI (Capmatinib, 100 nM) relative to those grown in 0.1% DMSO. (H) Growth assays of H596 cells cultured under indicated conditions to compare the impact of oncogenic STAT3, KRAS, or GFP expression on H596 cell growth in the presence or absence of HGF-driven proliferation. (I) Western blot analysis of MET effector proteins overexpressed by the pInducer20 promoter in the presence of doxycycline (200 ng/mL), along with their downstream signaling targets with or without HGF (10 ng/mL) or Capmatinib (100 nM) in H596 cells. (J,K) Viability assay in 6-well plates, as measured by (J) Alamar Blue and (K) Crystal Violet staining, demonstrating mutant KRAS-mediated rescue cell viability from MET TKI treatment in Hs746T, but not H1993, cells. (L,M) Western blot analysis of MET effector proteins overexpressed by the pInducer20 promoter in the presence of doxycycline (200 ng/mL), along with their downstream signaling targets with or without Capmatinib (100 ng/mL), in (L) H1993 and (M) Hs746T cells. n.s, not significant; * p < 0.05; *** p < 0.001 for the indicated comparisons.
Cancers 14 01378 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lu, D.; Nagelberg, A.; Chow, J.L.; Chen, Y.T.; Michalchuk, Q.; Somwar, R.; Lockwood, W.W. MET Exon 14 Splice-Site Mutations Preferentially Activate KRAS Signaling to Drive Tumourigenesis. Cancers 2022, 14, 1378. https://doi.org/10.3390/cancers14061378

AMA Style

Lu D, Nagelberg A, Chow JL, Chen YT, Michalchuk Q, Somwar R, Lockwood WW. MET Exon 14 Splice-Site Mutations Preferentially Activate KRAS Signaling to Drive Tumourigenesis. Cancers. 2022; 14(6):1378. https://doi.org/10.3390/cancers14061378

Chicago/Turabian Style

Lu, Daniel, Amy Nagelberg, Justine LM Chow, Yankuan T Chen, Quentin Michalchuk, Romel Somwar, and William W. Lockwood. 2022. "MET Exon 14 Splice-Site Mutations Preferentially Activate KRAS Signaling to Drive Tumourigenesis" Cancers 14, no. 6: 1378. https://doi.org/10.3390/cancers14061378

APA Style

Lu, D., Nagelberg, A., Chow, J. L., Chen, Y. T., Michalchuk, Q., Somwar, R., & Lockwood, W. W. (2022). MET Exon 14 Splice-Site Mutations Preferentially Activate KRAS Signaling to Drive Tumourigenesis. Cancers, 14(6), 1378. https://doi.org/10.3390/cancers14061378

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop