- Open Access
Identification of coexistence of BRAF V600E mutation and EZH2 gain specifically in melanoma as a promising target for combination therapy
Journal of Translational Medicine volume 15, Article number: 243 (2017)
Coexistence of enhancer of zeste homolog 2 (EZH2) and BRAF gene aberrations has been described in many cancer types. In this study, we aim to explore the coexistence status of BRAF V600E mutation and the copy number variation of EZH2 and explore the potential of this combination as a therapeutic target.
A total of 138 cases of melanoma samples harboring BRAF V600E mutation were included, and EZH2 copy numbers were examined by QuantiGenePlex DNA Assays. Clinical pathological distinction between patient groups with or without EZH2 amplification (hereafter referred to as EZH2 gain) was statistically analyzed. The sensitivity of melanoma cell lines and patient-derived xenograft (PDX) models containing BRAF V600E mutation with or without EZH2 gain to vemurafenib (BRAF inhibitor), GSK2816126 (EZH2 inhibitor) and a combination of both agents was evaluated.
In our cohort, the coexistence rate of BRAF V600E mutation and EZH2 gain was up to 29.0%, and significant differences in overall survival and disease-free survival were found between no EZH2 copy number gain and gain groups (P = 0.038, P = 0.030), gain and high EZH2 copy number gain groups (P = 0.006, P = 0.010). Combination with BRAF and EZH2 inhibition showed better inhibitory efficacy in melanoma prevention compared with vemurafenib monotherapy. More importantly, this improved therapeutic effect was observed especially in melanoma cell lines and PDX models containing concurrently BRAF V600E mutation and EZH2 gain.
Coexistence of BRAF V600E mutation and EZH2 gain is rather prevalent in melanoma. Our findings provided evidence for the feasibility of combination therapy with EZH2 and BRAF inhibitors in melanoma with concurrent BRAF V600E mutation and EZH2 gain.
The incidence of melanoma, one of the most malignant cancer, is increasing worldwide . The aggressiveness of melanoma is dependent on the high metastatic potential of melanoma cells, which can still not be effectively targeted despite recent progresses in targeted therapy and immunotherapy .
The mutation rates of BRAF V600E in Caucasians and Asians are approximately 50 and 25%, respectively . Vemurafenib, a BRAF inhibitor, has been shown to improve outcomes in the majority of melanoma patients harbouring BRAF V600E mutation, with a median overall survival (OS) of approximately 16 months . However, most patients treated with vemurafenib show disease progression within 6–8 months due to invariable drug resistance [4,5,6,7,8,9,10,11,12]. Recently, combined therapy has significantly improved response rates, along with progression-free and overall survival compared with single agent, such as vemurafenib plus trametinib, dabrafenib (BRAF inhibitor) plus trametinib, vemurafenib plus pembrolizumab and nivolumab plus ipilimumab [13,14,15]. However, despite rapid early response and high response rate to these combination therapeutic regimens, progression of disease occurs at a median of 11 months, with few patients remaining progression-free beyond 15 months , thus novel combination targets are urgently needed to be found.
The enhancer of zeste homolog 2 (EZH2) gene, located on chromosome 7q36.1, is adjacent to the BRAF gene which located on chromosome 7q34. Abnormalities in these two genes often coexist in various types of cancer, including papillary thyroid carcinoma . EZH2 is core component of the polycomb repressive complex 2, which catalyzes trimethylation of lysine 27 in histone 3 (H3K27me3), inducing chromatin compaction and preventing the transcription of target genes which are mostly tumor suppressor genes . Dysregulation of the EZH2 gene has been observed in several types of cancers, including lung, breast, and prostate cancer [17, 19, 20]. Growing evidence demonstrates that EZH2 is imperative for cancer initiation, development, progression, metastasis, and drug resistance. Therefore, EZH2 is currently considered a promising drug target, and multiple inhibitors of EZH2 have been developed, some of which are in clinical trials . Moreover, recent studies have shown that EZH2 also plays a critical role in the proliferation and survival of melanoma [19, 22,23,24,25,26,27,28].
EZH2 gain-of-function mutations often occur concurrently with the BRAF V600E mutation in melanoma . Knockdown of BRAF directly downregulates EZH2 gene expression in melanoma cells  and prostate cancer . Therefore, there may be a close relationship between aberrant EZH2 gene expression and the BRAF V600E mutation, thus highlighting EZH2 as a promising combination therapeutic target with vemurafenib for targeted therapy.
Herein, to reveal the potential of this combination as a therapeutic target, we will explore the coexistence status of BRAF V600E mutation and EZH2 gain in large-scale melanoma samples and show whether or not a combination of BRAF and EZH2 inhibition is effective for inhibition of melanoma cell growth in vitro and in vivo.
Patients and tumor tissue samples
This study involved samples from primary lesions of 138 patients with melanoma, who were hospitalized from January 2007 to January 2015 at Beijing Cancer Hospital & Institute. The diagnosis of melanoma was confirmed by hematoxylin and eosin staining and immunohistochemistry. BRAF V600E mutation status was verified by Sanger sequencing. This study was approved by the Medical Ethics Committee of the Beijing Cancer Hospital & Institute and was conducted according to the principles of the Declaration of Helsinki.
DNA preparation and BRAF V600E mutation screening
Genomic DNA was extracted from formalin-fixed, paraffin-embedded sections using a QIAamp DNA FFPE Tissue Kit (Qiagen, Valencia, CA, USA). To detect the BRAF V600E mutation, we amplified exon 15 of the BRAF gene by polymerase chain reaction (PCR) in at least two separate preparations of genomic DNA. The primer sequences were as follows: 15F1, 5′-AAACTCTTCATAATGCTTGCTC-3′; 15R1, 5′-TAATCAGTGGAAAAATAGCCTC-3′; 15F2, 5′-CTTGCTCTGATAGGAAAATG-3′; and 15R2, 5′-AGCCTCAATTCTTACCATCC-3′. All BRAF sequencing reactions were performed from forward and reverse directions and in duplicate. We purified PCR products with QIAquick (Qiagen) and directly sequenced them using Big Dye Terminator sequencing chemistry on an ABI 3130 automated sequencer (Applied Biosystems, Foster City, CA, USA). All mutations were confirmed by repeat bidirectional sequencing on the ABI sequencer.
QuantiGene Plex DNA assay
Tissue homogenates were prepared using the QuantiGene sample processing kit for formalin-fixed paraffin-embedded tissues (FFPE; Panomics of Affymetrics, Santa Clara, CA, USA) according to the manufacturer’s instructions. Briefly, 5–8 pieces of deparaffinized sections (4–10 μm) were incubated with 150 μl homogenizing solution supplemented with 1.5 μl of proteinase K (50 μg/μl) at 65 °C for 6 h. The tissue homogenate was separated from debris by brief centrifugation and transferred to a new tube.
The branched DNA (bDNA) assay was performed using the QuantiGene Plex DNA kit (Panomics) according to the manufacturer’s instructions. Briefly, the homogenate DNA was sheared using Covaris S2 (Covaris, Woburn, MA, USA) with the following settings: duty cycle 5%, intensity 3, cycles/burst 200, 80 s. For each assay well, 40 μl homogenate was denatured with 2.5 M NaOH (final concentration 0.18 M) in the presence of a DNA probe. The neutralized tissue homogenate was transferred to each well of hybridization plates containing the working bead mix. All samples were in processed in duplicates. The hybridization plate was sealed and incubated at 54 ± 1 °C in shaking incubator (600 rpm) for 18–22 h. The unbound samples were washed away using the Bio-plex Pro II wash station (Bio-Rad, Hercules, CA, USA). Then, the beads were sequentially hybridized with DNA pre-amplifier, DNA amplifier, labeled probe, and SAPE (streptavidin-conjugated R-phycoerythrin). Fluorescence intensities were measured using the Bio-plex 100 system (Bio-Rad).
The mean fluorescence intensities of the duplicates were calculated for all genes. Background values were subtracted from each probe set signal. Values of tested genes were normalized to the geometric means of Rpph1, Rpp30, and Rplp0. For each test gene, the normalized signal was divided by the signal of the reference DNA sample (G1521, Promega, Madison, WI, USA) and the values were multiplied by the known copy number (usually 2 copies) of each gene in the reference genome.
EZH2 gene copy number variation analysis by real-time PCR
To validate the results of QuantiGenePlex DNA Assays, the copy numbers of EZH2 were further quantified by TaqMan Copy Number Assays (Applied Biosystems, ThermoFisher, Waltham, MA, USA). A TaqMan probe targeting the Rnasep gene was used as a control. The EZH2 Tagman probe (4400291) was purchased from Invitrogen (ThermoFisher). Quantitative real-time PCR was performed using an ABI 7500 FAST real-time PCR system (Applied Biosystems). Copy numbers were then determined using CopyCaller v2.0 software (Applied Biosystems) with the comparative Ct (ΔΔCt) method. The EZH2 copy numbers for melanoma samples were calculated by dividing the sample values by the control values (benign nevus). No gain was considered for copy number ≤ 2, whereas gain was considered for copy number ≥ 2 and high gain was considered for copy number ≥ 4.
SK-MEL-5 (cat. no. HTB-70), WM-1664 (cat. no. CRL-1676), WM-115 (cat. no. CRL-1675), A2058 (cat. no. CRL-11147), and A375 cells were obtained from American Type Culture Collection (ATCC) and cultured at 37 °C in Dulbecco’s modified Eagle’s medium (DMEM; Invitrogen, ThermoFisher) supplemented with penicillin, streptomycin (Invitrogen), and 10% fetal bovine serum (HyClone, GE Healthcare, Logan, UT, USA).
All cells and tumor tissues of PDX models were lysed using the PhosphoSafe extraction reagent (Millipore). Supernatants were collected by centrifugation. Western blotting analysis of protein complexes were performed with the following antibodies: EZH2 (1:1000; Cell Signaling Technology, Danvers, MA, USA), H3K27me3 (1:1000; Cell Signaling Technology, Danvers, MA, USA), GAPDH (1:10,000; abcam, Cambridge, UK).
IHC analysis of EZH2, P-ERK and P-AKT protein was performed using anti-EZH2 polyclonal rabbit antibodies, anti-P-AKT polyclonal rabbit antibodies and anti-P-ERK polyclonal rabbit antibodies (Cell Signaling Technology, Danvers, MA, USA) at a dilution of 1:1000, followed by a standard avidin–biotin detection protocol using diaminobenzidine. Hematoxylin-counterstained slides were mounted and examined to determine the intensity of staining. The staining intensity was scored as −, +, ++, or +++ (with − indicating negative staining, and +++ indicating the strongest staining) independently by three pathologists.
Cell cycle and apoptosis assays
Cell cycle and apoptosis assays were performed using A2058 cells at 24 or 48 h after exposure to vemurafenib, GSK126, or both vemurafenib and GSK126 using a BD Pharmingen (BD Biosciences, Sparks, MD, USA) and Annexin V FITC Apoptosis Detection Kit (DOJINDO, Shanghai, China) according to the manufacturers’ protocols. Analysis was performed on a FACS Calibur flow cytometer (BD Biosciences).
Effects of EZH2 and BRAF inhibitors on proliferation
GSK126 and vemurafenib were purchased from Selleck Chemicals (Houston, TX, USA). All inhibitors were dissolved at 10 mM in dimethylsulfoxide (DMSO) as stock solutions. After treatment with various concentrations of inhibitors or DMSO for 24 h, cell proliferation was evaluated by CellTiter-Glo Luminescent Cell Viability Assays (Promega, Madison, WI, USA) according to the manufacturer’s instructions.
Calculation of the combination index (CI)
The combined activity of vemurafenib and GSK126 was determined by calculating the CI for both compounds in A2058 and WM115 cells with the following equation in Compusyn software: CI = (D1/Dx1) + (D2/Dx2) + α [(D1 × D2)/(Dx1 × Dx2)], where Dx1 is the concentration of drug 1 required to produce 50% cell death alone, D1 is the concentration of drug 1 required to produce 50% cell death in combination with D2, Dx2 is the concentration of drug 2 required to produce 50% cell death alone, D2 is the concentration of drug 2 required to produce 50% cell death in combination with D1, and α is 0 for mutually exclusive or 1 for mutually nonexclusive modes of drug action. The results are interpreted as: 0 < CI < 1 indicates synergism; CI = 1 indicates an addictive effect; and CI > 1 indicates antagonism.
Patient-derived xenograft (PDX) model and treatment
Fragments of patient-derived melanoma tissues bearing BRAF V600E mutation were cut into fragments and then subcutaneously inoculated into a 5 week-old NOD/SCID (non-obese diabetic and severe combined immunodeficiency) female mouse (4–6 week-old, 18–22 g-weight) to establish the PDX model. When the tumor size reached approximately 1 cm3, the mice were sacrificed, and tumor tissues were separated and re-inoculated into new mice. 2 PDX models (detailed information are listed in Additional file 1: Table S1) were finally established: PDX 001 (concurrently containing BRAF V600E mutation and EZH2 gain), and PDX 002 (containing BRAF V600E mutation alone).
When the tumor size reached approximately 600 mm3, each type of PDX models were randomized (treatment arm versus control arm; n = 5) into four groups: group 1 received oral vemurafenib (30 mg/kg) twice a day; group 2 received GSK126 (50 mg/kg) intraperitoneally once a day; group 3 received a combination of the two agents with doses and administration methods similar to those introduced in single treatment; and group 4 served as a negative control and received vehicle once a day by intraperitoneal injection of 100 µl corn oil and oral gavage of 100 µl PBS. Body weights and the lengths and widths of tumors were measured every 3 days using scales and calipers, and the tumor volume was calculated as ([length × width2] × 0.5). After 14 days, the mice were sacrificed, and the tumors were fixed in 10% formalin for histological and immunohistological analysis. All animal care and experimental procedures were carried out in accordance with the Animal Care Ethics guidelines approved by the Medical Ethics Committee of Beijing Cancer Hospital & Institute.
Statistical analyses were performed using SPSS 22.0 software. Continuous data, such as age and lesion thickness, were described using means ± standard deviations (SDs) for normally distributed data. The correlations between aberration status and clinical parameters were evaluated by Chi square test or Fisher’s exact test. Kaplan–Meier estimates of time-to-event OS were calculated. Log-rank tests were used to estimate the statistical significance between the time-dependent outcomes of OS. Cox hazard proportion models were used to estimate the hazard ratios (HRs) and corresponding 95% interval confidences (CIs). All statistical analyses were two-sided, and differences with P values of less than 0.05 were considered statistically significant.
EZH2 gain in BRAF V600E-mutated melanoma
Among the 138 patients with BRAF V600E-mutated melanoma, 40 cases (29.0%) showed EZH2 gain (Table 1). Additionally, we detected the EZH2 gain frequency in melanoma subtypes and found that the frequencies of EZH2 gain in acral melanoma, mucosal melanoma, melanoma on skin with chronic sun-induced damage (CSD) and melanoma on skin without chronic sun-induced damage (non-CSD) were 21.9, 38.9, 25.0 and 33.3%, respectively (Table 1).
Altogether, our data suggested that the coexistence of EZH2 gain and the BRAF V600E mutation was prevalent in melanoma, especially mucosal subtype.
Correlation of EZH2 gain with clinicopathological features in patients with BRAF V600E-mutated melanoma
A summary of the correlations between EZH2 gain and clinical characteristics are shown in Table 2. The mean age (P = 0.312) and gender (P = 0.553) were not significantly different between melanoma patients harbouring BRAF V600E mutation with or without EZH2 gain.
Thickness, ulceration, and stage are important clinical characteristics that can be used to predict the prognosis in patients with melanoma. In our cohort, the data for stage (138 cases), thickness (59 cases), and ulceration (122 cases) were available for further analysis (Fig. 1a).
The average thickness of the melanoma lesions in 59 patients with the BRAF V600E mutation was 4.8 mm, consistent with the observation that thickness in Chinese patients is thicker than that in Caucasians . Additionally, the average thickness of lesions in melanomas harboring EZH2 gain was 5.525 ± 3.475 mm, whereas that in lesions without EZH2 gain was 4.65 ± 3.35 mm; there was no significant difference between these two groups (P = 0.907).
Ulceration is another adverse prognostic factor in melanoma. In our cohort, the overall ulceration rate for the 122 available cases was 59.5% (72/121), and the ulceration rates in acral, mucosal, CSD, and non-CSD melanomas were 61.0% (25/41), 44.4% (8/18), 30.6% (11/36), 41.0% (16/39), respectively. The ulceration rate in melanomas harboring EZH2 gain was 35.3%, whereas in melanomas without EZH2 gain was 69% (P = 0.001). Thus, EZH2 gain did not appear to be correlated with the formation of ulceration, further studies with larger sample sizes are needed to confirm this result.
Among melanoma patients containing concurrently BRAF V600E mutation and EZH2 gain, the percentages of patients with stages I, II, III and IV were 10% (four cases), 12.5% (five cases), 42.5% (17 cases), and 35.0% (14 cases), respectively (P = 0.071); these data were not significantly different from those in patients without EZH2 gain. However, the percentage of patients with stages III and IV in the EZH2 gain group was 77.5%, which was substantially higher than that in the EZH2 no gain group (67.3%); thus, further studies with larger sample sizes are needed to confirm the differences in disease stage between the two groups.
Moreover, correlation of EZH2 gain with clinicopathological features among three groups (no gain, gain, high gain) were also performed, tables were provided in Additional file 2.
Prognostic significance of EZH2 gain for OS and DFS in patients with BRAF V600E mutated melanoma patients
The thickness and ulceration of primary melanoma lesions and the stage of the disease are adverse prognostic factors. Consistent with previous reports, we found that the overall survival (OS) of patients with thicker lesions and more advanced-stage disease was significantly shorter than those with thinner lesions (P = 0.023) and early-stage disease (P < 0.001), respectively. In addition, in our cohort, the OS of patients with mucosal melanoma and unknown primary melanoma was significantly shorter than that of patients with other subtypes (P < 0.001; data not shown).
Next, we analyzed the prognostic significance of EZH2 gain for overall survival (OS) and disease-free survival (DFS) in melanoma patients harbouring BRAF V600E mutation. Significant differences in overall survival and disease-free survival were found between no EZH2 copy number gain and gain groups (P = 0.038, HR 0.50, 95% CI 0.02–0.61; P = 0.030, HR 0.49, 95% CI 0.26–0.94), gain and high EZH2 copy number gain groups (P = 0.006, P = 0.010) (OS median value: no gain group: 78.1 months; gain group: 40.8 months; high gain group: 31.3 months; DFS median value: no gain group: 39.5 months; gain group: 22.3 months; high gain group: 11.1 months) (Fig. 1b).
In univariate Cox analysis, the clinicopathologic factors, such as TNM stage and EZH2 amplification, may be of prognostic significance for melanoma patients; For multivariate Cox regression assay, TNM stage is independent prognostic factors for OS (Additional file 3).
Taken together, our findings indicated that EZH2 gain may be an adverse prognostic biomarker in BRAF V600E-mutated melanoma.
EZH2 copy number levels and corresponding EZH2 protein levels in four BRAF V600E-mutated cell lines, PDX models and patient tissues
EZH2 copy number levels, and corresponding western blotting results for four cell lines and two PDX models are shown in Fig. 2a–d. The copy number decreased in the order of A2058 > SK-MEL-5 > A375 > WM115 in cell lines and PDX-001 > PDX-002 in mouse models, corresponding protein levels were roughly consistent with copy number levels.
We also examined EZH2 protein expression and H3K27me3 levels by IHC (Fig. 2e) which supported our previous detected EZH2 copy number aberrations. H3K27me3 levels reflected the active status of EZH2 protein.
These findings indicated that EZH2 protein expression was positively correlated with copy number alterations.
Inhibitory effect of combination on BRAF V600E-mutated cell lines
A2058, WM115, A375, and SK-MEL-5, all four BRAF V600E-mutated cell lines were evaluated to determine the efficacy of GSK126 and vemurafenib by analysis of cell proliferation rates in vitro. The concentration–response curve for each cell line is shown in Fig. 3a, b. For vemurafenib, WM115 cells (half-maximal inhibitory concentration [IC50] = 1.7 μM) were the most sensitive, and A2058 cells (IC50 = 4 μM) were the least sensitive. For GSK126, A2058 cells (IC50 = 4 μM) were the most sensitive, whereas WM115 cells (IC50 = 10 μM) were the least sensitive. The IC50 values were then used to generate fixed ratios for subsequent combination treatment regimens and calculation of CIs.
Combination therapy with these two agents was also performed in all four cell lines. The concentration–response curve is shown in Fig. 3c–f. For evaluation of the combination of vemurafenib and GSK126, the cells were exposed to medium containing both drugs concurrently at a fixed ratio approximating their individual IC50s for 48 h to ensure that both agents contributed equally to the observed antiproliferative effects. As shown in Fig. 3c, for A2058 cells, the drug combination showed increased inhibitory efficacy compared either agent used alone; for SK-MEL-5 cells, the drug combination showed advantages only under high doses condition (vemurafenib: 0.8 μM/l, GSK126: 5.6 μM/l); for A375 cells, combination therapy showed advantages only under low doses condition (vemurafenib: 0.8 μM/l, GSK126: 9.6 μM/l). However, in WM115 cells, combination therapy showed no obvious advantage. To further elucidate the interaction patterns of these two drugs, CIs and fractional cell growth inhibition (Fa) values were calculated (Fig. 3g–j). CIs of significantly less than 1 were obtained in A2058 cells (mean CI at Fa 0.5 = 0.885), SK-MEL-5 cells (mean CI at Fa 0.5 = 1.260) under high doses condition and A375 cells (mean CI at Fa 0.5 = 0.972) under low doses condition, demonstrating that the two drugs interacted synergistically to inhibit cell growth totally or partially. In contrast, values more than 1 were obtained in WM115 cells, indicating that the two drugs interacted antagonistically (mean CI at Fa 0.5 = 1.370) (Compusyn reports can be seen in Additional files 4, 5, 6, 7 and 8). After treatment, P-AKT protein levels were detected in GSK126 monotherapy groups of 4 cell lines, we found that EZH2 inhibition was associated with downregulation of P-AKT levels (Additional file 9: Figure S3).
This result indicated that combination therapy could achieve better anticancer efficacy only in melanoma cell lines with concurrent BRAF V600E mutation and EZH2 gain, and the effects of GSK126 treatment may involve inhibition of PI3K signaling pathway.
Variations in cell cycle distributions and apoptosis rates in cells after combination therapy
We next examined whether the suppression of cell proliferation by combination therapy was caused by cell cycle arrest or apoptosis. A2058 cells were used for further study.
The cell cycle phase distribution of cells treated with vemurafenib (1 μM), GSK126 (1 μM), and combined treatment (1 μM each) for 48 h was detected by flow cytometry. We found that compared with the monotherapy group, cells treated with combined treatment showed a significant increase in the proportion of cells at G1 phase and a reduction in the proportion of cells in the S phase (P < 0.05; Fig. 4a, b). The percentages of cells in the G1 phase following combination therapy increased by approximately 24.9, 11.14, and 14.89% compared with the control, vemurafenib, and GSK126 groups, respectively. In addition, Sub-G0 population detected by flow cytometry showed that combination therapy caused more obvious cell death than monotherapy groups (Additional file 9: Figure S1).
The apoptosis rates in A2058 cells treated with vemurafenib (1.2 μM), GSK126 (1.2 μM), or both drugs (1.2 μM each) for 24 h were detected by flow cytometry with Annexin V-FITC/propidium iodide (PI) double staining. We found that combined treatment caused obvious apoptosis effects, with apoptosis rates of 53.2, 39.7, and 40.3% more than those of the control, vemurafenib, and GSK126 groups, respectively. Results of western blotting for PARP (apoptosis biomarker) further confirmed this conclusion (Fig. 4e). Then, we probed this result in other three cell lines, we found that combination therapy caused more apparently cell apoptosis especially in cells harboring EZH2 gain. (A2058, SK-MEL-5 and A375) (Additional file 9).
Taken together, these results indicated that vemurafenib and GSK126 could function synergistically to induce cell cycle arrest and apoptosis, especially cells concurrently harboring BRAF V600E mutation and EZH2 gain.
Inhibitory effect of combination on PDX models
To verify the results in cell lines, in vivo assays were performed using PDX mouse models. We found that the combined treatment with vemurafenib and GSK126 significantly inhibited the growth of xenograft tumors in PDX models, especially PDX 001 (concurrently containing BRAF V600E and EZH2 gain) compared with treatment with single drug alone (Fig. 5a, c; P < 0.05). In terms of side effects, we did not observe a significant decrease in body weight (Fig. 5e, f). In contrast, in PDX 002 (without EZH2 gain), the inhibitory efficacy of combined treatment was not significantly different from that of vemurafenib alone but was stronger than single treatment with GSK126 (Fig. 5b, d). Moreover, vemurafenib alone showed stronger inhibitory efficacy than GSK126 alone in two types of PDX mouse models, whereas PDX 001 were more sensitive to GSK126 than PDX 002 models. As further evidence of the inhibitory effects of the three therapeutic regimens on tumor growth, we examined the proliferation of tumor cells in all mice by IHC analysis of Ki-67 (Fig. 6a, b). At last, tumor growth inhibition rate of three therapeutic regimens in same PDX model and each regimen in two types of PDX models was calculated (Figs. 5g, h, 6c). The findings were consistent with the results of changes in tumor volumes. In addition, IHC analysis of P-AKT and P-ERK (Fig. 6a) was also performed, the P-AKT protein levels obviously downregulated in GSK126 monotherapy and combination therapy groups whereas P-ERK protein levels obviously downregulated in vemurafenib monotherapy and combination therapy groups.
These results showed that the combined treatment with vemurafenib and GSK126 inhibited tumor growth of PDX models more efficiently than a single treatment with vemurafenib, especially in those concurrently containing BRAF V600E mutation and EZH2 gain with no obvious toxicity.
Novel therapeutic strategies for the treatment of advanced melanoma is the urgent priority all long, especially BRAF V600E mutated melanoma . Even though BRAF inhibitor-vemurafenib represented remarkable clinical efficacy, most patients show disease progression within 6–8 months .
Recently, progress in combination therapies of melanoma treatment have significantly addressed the limitation of vemurafenib alone, the most successful combination therapeutic regimen was dabrafenib (the second BRAF inhibitor approved by FDA) plus trametinib, which induced complete and partial responses in almost 70% of patients and prolonged median progression-free survival in several phase 3 clinical trials, approximately 3–5 months longer than single vemurafenib alone . Therefore, this combination regimen was identified by US NCCN guideline as first-line targeted therapy of advanced melanoma patients containing BRAF V600E mutation. However, in many patients with dabrafenib plus trametinib treatment, the disease rapidly progresses and leading to death, such as patients with extremely high LDH concentration  and brain metastatic lesions , lacking of effective targeted therapies in clinical practice. Therefore, novel combination targets are urgently needed to be found.
Given the high coexistence rate (29%) of EZH2 gain and BRAF V600E mutation in our cohort, we for the first time evaluated the efficacy of combination therapy with EZH2 and BRAF inhibitors in melanoma cell lines and PDX mouse models. We demonstrated that combination therapy was more efficacious than single agent alone in vitro and in vivo, whereas only in melanoma cell lines and PDX models concurrently containing BRAF V600E mutation and EZH2 gain. Therefore, we provided a promising novel therapeutic strategy for melanoma. In patients with BRAF V600E-mutated melanoma, the individual EZH2 gain status should be tested first, and cotreatment with BRAF and EZH2 inhibitors can be applied to patients containing EZH2 gain and BRAF V600E mutation simultaneously to enhance the anticancer effects.
Notably, in our cohort, the frequencies of EZH2 gain in mucosal (38.9%) and unknown primary melanoma (50%) are markedly higher than acral (21.9%), CSD (25.0%) and non-CSD melanoma (33.3%) subtypes. Therefore, combination therapy with BRAF and EZH2 inhibitors may be particularly beneficial to patients with mucosal and unknown primary melanoma. Moreover, it has been reported that increased activation of the phosphatidylinositol 3-kinase-protein kinase B (PI3K-AKT) signaling pathway has been associated with resistance to BRAF and MEK inhibitors . Maria E and colleagues found EZH2 overexpression is sufficient for activation of the PI3K-Akt pathway . According to our western blot and IHC results for P-AKT protein, EZH2 inhibition was significantly associated with downregulation of P-AKT levels. Thus, combination with EZH2 inhibitors may help to down-regulate the activity of PI3K-AKT pathway and deal with resistance to BRAF inhibitor alone and BRAF plus MEK inhibitors.
We acknowledge several limitations exist in our study. Firstly, the detailed mechanisms involved in the molecular interactions of BRAF and EZH2 gene are still unknown. We speculate that intersection exists between BRAF and EZH2 signaling pathway. Evidence shows that BRAF signaling has a role in regulating DNA methylation status, and EZH2 is a methylation target of BRAF . Given the central role of gene methylation in tumorigenesis, EZH2 and BRAF may cooperate to promote the initiation and progression of melanoma . In addition, BRAF is the upstream component of the mitogen-activated protein kinase (MAPK) pathway that regulates cell proliferation and survival . EZH2 is the downstream component of CDK4 pathway that drives cell-cycle progression and regulates cell proliferation . Given the obvious cell-cycle arrest and apoptosis induction was observed after combined treatment, BRAF and EZH2 may cooperate in accelerating cell-cycle progression and inhibiting apoptosis to promote melanoma cell growth. These hypotheses need further study. Secondly, sensitivity of cell lines and PDX models with different EZH2 gain levels to combined treatment should have been tested to verify whether efficacy of combination therapy is positively correlated with EZH2 gain levels. In vitro results of our cell lines basically confirmed this idea, however, more cell lines and PDX models with representative, different EZH2 gain levels should be established to further verify this hypothesis.
Thirdly, although no obvious side effects were observed in our PDX models of combination treatment, further large-scale clinical trials are needed to test the toxicities, potential drug–drug interactions, and optimize concentrations of these two agents in patients with melanoma. Lastly, the effects of other EZH2 inhibitors, such as EPZ6438 and GSK343, still need to be clarified. All above will be the focus of our subsequent studies.
In this study, we indicated that coexistence of BRAF V600E mutation and EZH2 gain is rather prevalent in melanoma. Combination with BRAF and EZH2 inhibition showed better inhibitory efficacy compared with vemurafenib monotherapy in vitro and in vivo, especially in melanoma containing concurrently BRAF V600E mutation and EZH2 gain. In conclusion, our study identified EZH2 gain as a biomarker for selecting patients with BRAF V600E-mutated melanoma who may benefit from combination therapy with BRAF and EZH2 inhibitors, thus providing insights into a promising novel therapeutic strategy for the treatment of melanoma.
Kong Y, Si L, Zhu Y, et al. Large-scale analysis of KIT aberrations in Chinese patients with melanoma. Clin Cancer Res. 2011;17(7):1684–91.
Zingg D, Debbache J, Schaefer SM, et al. The epigenetic modifier EZH2 controls melanoma growth and metastasis through silencing of distinct tumour suppressors. Nat Commun. 2015;6:6051.
Si L, Kong Y, Xu X, et al. Prevalence of BRAF V600E mutation in Chinese melanoma patients: large scale analysis of BRAF and NRAS mutations in a 432-case cohort. Eur J Cancer. 2012;48(1):94–100.
Young K, Minchom A, Larkin J. BRIM-1, -2 and -3 trials: improved survival with vemurafenib in metastatic melanoma patients with a BRAF(V600E) mutation. Future Oncol. 2012;8(5):499–507.
Sosman JA, Kim KB, Schuchter L, et al. Survival in BRAF V600-mutant advanced melanoma treated with vemurafenib. N Engl J Med. 2012;366(8):707–14.
Shaw AT, Engelman JA. ALK in lung cancer: past, present, and future. J Clin Oncol. 2013;31(8):1105–11.
Villanueva J, Vultur A, Lee JT, et al. Acquired resistance to BRAF inhibitors mediated by a RAF kinase switch in melanoma can be overcome by co-targeting MEK and IGF-1R/PI3K. Cancer Cell. 2010;18(6):683–95.
Nazarian R, Shi H, Wang Q, et al. Melanomas acquire resistance to B-RAF(V600E) inhibition by RTK or N-RAS upregulation. Nature. 2010;468(7326):973–7.
Poulikakos PI, Persaud Y, Janakiraman M, et al. RAF inhibitor resistance is mediated by dimerization of aberrantly spliced BRAF(V600E). Nature. 2011;480(7377):387–90.
Wagle N, Emery C, Berger MF, et al. Dissecting therapeutic resistance to RAF inhibition in melanoma by tumor genomic profiling. J Clin Oncol. 2011;29(22):3085–96.
Lito P, Pratilas CA, Joseph EW, et al. Relief of profound feedback inhibition of mitogenic signaling by RAF inhibitors attenuates their activity in BRAFV600E melanomas. Cancer Cell. 2012;22(5):668–82.
Liu S, Zhang B, Zhao Y, et al. Association of BRAFV600E mutation with clinicopathological features of papillary thyroid carcinoma: a study on a Chinese population. Int J Clin Exp Pathol. 2014;7(10):6922–8.
Robert C, Karaszewska B, Schachter J, et al. Improved overall survival in melanoma with combined dabrafenib and trametinib. N Engl J Med. 2015;372(1):30–9.
Long GV, Grob JJ, Nathan P, et al. Factors predictive of response, disease progression, and overall survival after dabrafenib and trametinib combination treatment: a pooled analysis of individual patient data from randomised trials. Lancet Oncol. 2016;17(12):1743–54.
Khushalani NI, Sondak VK. Are we there yet? Prolonged MAPK inhibition in BRAF(V600)-mutant melanoma. Lancet Oncol. 2016;17(9):1178–9.
Paulson KG, Thompson JA. A second chance for success with BRAF and MEK inhibitors in melanoma. Lancet Oncol. 2017;18(4):418–9.
Barsotti AM, Ryskin M, Zhong W, et al. Epigenetic reprogramming by tumor-derived EZH2 gain-of-function mutations promotes aggressive 3D cell morphologies and enhances melanoma tumor growth. Oncotarget. 2015;6(5):2928–38.
Kim NY, Pyo JS. Clinicopathological significance and prognostic role of EZH2 expression in non-small cell lung cancer. Pathol Res Pract. 2017;213(7):778–82.
Bachmann IM, Halvorsen OJ, Collett K, et al. EZH2 expression is associated with high proliferation rate and aggressive tumor subgroups in cutaneous melanoma and cancers of the endometrium, prostate, and breast. J Clin Oncol. 2006;24(2):268–73.
Liu C, Liu R, Zhang D, et al. MicroRNA-141 suppresses prostate cancer stem cells and metastasis by targeting a cohort of pro-metastasis genes. Nat Commun. 2017;8:14270.
Yamaguchi H, Hung MC. Regulation and role of EZH2 in cancer. Cancer Res Treat. 2014;46(3):209–22.
Mahmoud F, Shields B, Makhoul I, et al. Role of EZH2 histone methyltrasferase in melanoma progression and metastasis. Cancer Biol Ther. 2016;17(6):579–91.
Fane ME, Chhabra Y, Hollingsworth DE, et al. NFIB mediates BRN2 driven melanoma cell migration and invasion through regulation of EZH2 and MITF. EBioMedicine. 2017;16:63–75.
Montagnani V, Benelli M, Apollo A, et al. Thin and thick primary cutaneous melanomas reveal distinct patterns of somatic copy number alterations. Oncotarget. 2016;7(21):30365–78.
Gallagher SJ, Tiffen JC, Hersey P. Histone modifications, modifiers and readers in melanoma resistance to targeted and immune therapy. Cancers. 2015;7(4):1959–82.
McHugh JB, Fullen DR, Ma L, Kleer CG, Su LD. Expression of polycomb group protein EZH2 in nevi and melanoma. J Cutan Pathol. 2007;34(8):597–600.
Fan T, Jiang S, Chung N, et al. EZH2-dependent suppression of a cellular senescence phenotype in melanoma cells by inhibition of p21/CDKN1A expression. Mol Cancer Res. 2011;9(4):418–29.
Chen X, Wang J, Shen H, et al. Epigenetics, microRNAs, and carcinogenesis: functional role of microRNA-137 in uveal melanoma. Investig Ophthalmol Vis Sci. 2011;52(3):1193–9.
Chen H, Tu SW, Hsieh JT. Down-regulation of human DAB2IP gene expression mediated by polycomb Ezh2 complex and histone deacetylase in prostate cancer. J Biol Chem. 2005;280(23):22437–44.
Hou P, Liu D, Dong J, Xing M. The BRAF(V600E) causes widespread alterations in gene methylation in the genome of melanoma cells. Cell Cycle. 2012;11(2):286–95.
Li L, Zheng H, et al. DNA methylation signatures and coagulation factors in the peripheral blood leucocytes of epithelial ovarian cancer. Carcinogenesis. 2017;38(8):797–805.
Davies MA, Saiag P, Robert C, et al. Dabrafenib plus trametinib in patients with BRAFV600-mutant melanoma brain metastases (COMBI-MB): a multicentre, multicohort, open-label, phase 2 trial. Lancet Oncol. 2017;7(18):836–7.
Gonzalez ME, DuPrie ML, Krueger H, et al. Histone methyltransferase EZH2 induces Akt-dependent genomic instability and BRCA1 inhibition in breast cancer. Cancer Res. 2011;71(6):2360–70.
Obenauf AC, Zou Y, Ji AL, et al. Therapy-induced tumour secretomes promote resistance and tumour progression. Nature. 2015;520(7547):368–72.
Perico L, Mandalà M, Schieppati A, et al. BRAF signaling pathway inhibition, podocyte injury, and nephrotic syndrome. Am J Kidney Dis. 2017;70(1):145–50.
Santos M, Martínez-Fernández M, Dueñas M, et al. In vivo disruption of an Rb-E2F-Ezh2 signaling loop causes bladder cancer. Cancer Res. 2014;74(22):6565–77.
YH, KY designed the study. YH, XLW performed the experiments. YH, YJY analyzed the data. YH, KY wrote the manuscript. All authors read and approved the final manuscript.
We thank Editage company for English language editing.
The authors declare that they have no competing interests.
Availability of data and materials
All the data and materials supporting the conclusions were included in the main paper.
Consent for publication
Ethics approval and consent to participate
This study was approved by the Medical Ethics Committee of the Beijing Cancer Hospital & Institute and was conducted according to the principles of the Declaration of Helsinki. All animal care and experimental procedures were carried out in accordance with the Animal Care Ethics guidelines approved by the Medical Ethics Committee of Beijing Cancer Hospital & Institute.
This work was supported by grants from the Major State Basic Research Development Program of China (2013CB911004), National Natural Science Foundation of China (81672696), Beijing Municipal Natural Science Foundation (7152033), Beijing Baiqianwan Talents Project, Beijing Municipal Administration of Hospitals Clinical medicine Development of special funding support (ZYLX201603), and Beijing Municipal Science & Technology Commission (Z151100003915074).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 2. EZH2 gain in BRAF V600E mutated melanoma. Correlation of EZH2 gain to clinicopathologic features of BRAF V600E mutated melanomas.
Additional file 3. Univariate and multivariate analysis of risk factors associated with overall survival.
Additional file 5. Correlation of EZH2 amplification to clinicopathologic features of BRAFV600E mutated mucosal melanomas.
Additional file 9: Figure S1. Sub-G0 cells detected in 2058 cells after combination therapy. Figure S2. Variations in apoptosis rate in all four BRAF V600E mutated cell lines after combination therapy. Figure S3. The levels of P-AKT at baseline and after treatment with GSK drug in all cell lines.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Yu, H., Ma, M., Yan, J. et al. Identification of coexistence of BRAF V600E mutation and EZH2 gain specifically in melanoma as a promising target for combination therapy. J Transl Med 15, 243 (2017). https://doi.org/10.1186/s12967-017-1344-z
- EZH2 gain
- BRAF V600E mutation
- Combination therapy