Next Article in Journal
Understanding the Role of Innate Immune Cells and Identifying Genes in Breast Cancer Microenvironment
Next Article in Special Issue
Design and Testing of a Custom Melanoma Next Generation Sequencing Panel for Analysis of Circulating Tumor DNA
Previous Article in Journal
Antibody-Drug Conjugates: A Promising Novel Therapy for the Treatment of Ovarian Cancer
Previous Article in Special Issue
Simultaneously Inhibiting BCL2 and MCL1 Is a Therapeutic Option for Patients with Advanced Melanoma
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 12
Issue 8
10.3390/cancers12082224
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

TERT Promoter Mutation as an Independent Prognostic Marker for Poor Prognosis MAPK Inhibitors-Treated Melanoma

by
Pauline Blateau
1,2,
Etienne Coyaud
3,
Estelle Laurent
3,
Benoit Béganton
1,2,
Vincent Ducros
1,
Géraldine Chauchard
1,
Julie A. Vendrell
1 and
Jérôme Solassol
1,2,*
1
Laboratoire de Biologie des Tumeurs Solides, Département de Pathologie et Oncobiologie, Centre Hospitalier Universitaire de Montpellier, 34000 Montpellier, France
2
Institut de Recherche en Cancérologie de Montpellier, INSERM, Université de Montpellier, Institut du Cancer de Montpellier, Université de Montpellier, 34000 Montpellier, France
3
Laboratoire Protéomique Réponse Inflammatoire Spectrométrie de Masse (PRISM), INSERM U1192, Université de Lille, Centre Hospitalier Universitaire Lille, F-59000 Lille, France
*
Author to whom correspondence should be addressed.
Cancers 2020, 12(8), 2224; https://doi.org/10.3390/cancers12082224
Submission received: 2 June 2020 / Revised: 23 July 2020 / Accepted: 3 August 2020 / Published: 9 August 2020
(This article belongs to the Special Issue Advanced Melanoma)
Browse Figures
Versions Notes

Abstract

:
Although the development of mitogen-activated protein kinase (MAPK) inhibitors has greatly improved the prognosis of BRAFV600 cutaneous melanomas, the identification of molecular indicators for mutated patients at risk of early progression remains a major issue. Using an amplicon-based next-generation-sequencing (NGS) assay that targets cancer-related genes, we investigated co-occurring alterations in 89 melanoma samples. We analyzed both their association with clinicopathological variables and clinical significance in terms of progression-free survival (PFS) and overall survival (OS) according to BRAF genotyping. Among co-occurring mutations, TERT promoter was the most frequently mutated gene. Although no significant difference in PFS was observed in the presence or absence of co-occurring alterations to BRAFV600, there was a trend of longer PFS for patients harboring TERT c.-124C>T mutation. Of most interest, this mutation is an independent marker of good prognosis in subgroups of patients with poor prognosis (presence of brain metastasis and elevated level of lactate dehydrogenase, LDH). Moreover, combination of elevated LDH level, presence of brain metastasis, and TERT c.-124C>T mutation was identified as the best fit model for predicting clinical outcome. Our work revealed the potential interest of c.-124C>T status determination in order to refine the prognosis of BRAFV600 melanoma under mitogen-activated protein kinase (MAPK) inhibitors.

1. Introduction

Over the past few years, the molecular characterization of melanomas has greatly improved, with an emphasis on alteration of cell signaling pathways [1,2]. Approximately 40% of patients with melanoma exhibit exon 15 BRAF mutations in cancer cells, resulting in constitutive activation of the mitogen-activated protein kinase (MAPK) cascade. A therapeutic strategy based on dual inhibition of the MAPK pathway through targeting BRAF and MEK proteins with BRAF inhibitors (e.g., dabrafenib or vemurafenib) in combination with MEK inhibitors (e.g., trametinib or cobimetinib) has significantly improved progression-free survival (PFS) and overall survival (OS) in melanoma patients harboring activating BRAF mutations [3]. Concurrently, immune checkpoint inhibitors targeting Programmed Death -1 (PD-1) and cytotoxic T-lymphocyte associated protein 4 (CTLA-4) showed clinically significant improvements in OS in molecularly unselected populations of advanced melanoma patients. Recent data support the hypothesis that these therapies also provide clinical benefit in melanoma patients with activating BRAF mutations [4].
Although these therapies have significantly improved the prognosis of melanoma advanced forms, their effectiveness in practice remains subject to significant interpersonal variation between patients, with some patients showing primary resistance or early progression. Within this group, prognostic factors conventionally useful in distinguishing individuals at risk of poor clinical outcome or progression from others include the following: stage of disease; baseline serum lactate dehydrogenase (LDH) levels; presence of brain metastases, and the Eastern Cooperative Oncology Group’s (ECOG PS) baseline performance status [5]. However, these prognostic features have been validated years before the advent of targeted therapies and use of BRAF and MEK inhibitors. Hence, they appear as poorly suitable for the genotyping status-based stratification of melanoma patients.
With the recent emergence of next-generation-sequencing (NGS) analyses, concomitant somatic genomic alterations have been identified in samples of BRAF mutant melanomas [6,7,8,9], such as CDKN2A, PTEN, RAC1, and TERT promoter [10]. Most of these co-occurring mutations have been studied individually, leading in some cases to the identification of resistance mechanisms against BRAF and MEK inhibition therapy, such as the activation of the MAPK or PI3K/AKT pathway [11,12,13,14]. However, the relationships between BRAFV600 mutation, concomitant molecular alterations, and clinicopathological parameters as prognosis markers have been barely investigated in melanoma.
Thus, there is a critical need for identifying and validating new tissue biomarkers able to refine patient prognosis and classify them into specific subgroups more or less likely to progress early under targeted therapy, especially in patients with distant metastasis. To tackle this important clinical question, we performed NGS analysis on a cohort of 113 cutaneous melanoma patients focusing on 35 clinically relevant cancer hotspot regions. We analyzed the frequency and type of co-occurring mutations and their association with BRAF status. We also evaluated the correlation between concomitant genomic alterations in BRAF mutant melanomas with their clinical and pathological characteristics, as well as their potential synergistic effect on patient outcome.

2. Results

2.1. Patient Characteristics

A total of 113 samples of cutaneous melanoma were collected and exhaustively analyzed by NGS between April 2014 and September 2019 at the Pathology Laboratory of the University Hospital of Montpellier, France, to assess the presence of molecular alterations (Figure 1). Patients eligible for this retrospective study were diagnosed either for primary or recurrent metastatic melanoma. Their clinicopathological features are shown in Table S1.
The dropout (n = 24) was based on poor DNA quality or lost-to-follow-up. Specifically, we observed that 53 samples (59.6%) harbored a BRAFV600 mutation and 36 samples (40.4%) were BRAF wild type (BRAFWT) (Figure 1).

2.2. NGS Analysis and Co-Occurring Genetic Alteration Detection

We next determined the prevalence of co-occurring aberrations in the 89 formalin-fixed paraffin-embedded tumor samples analyzed by NGS.
In BRAFWT specimens, 30 samples (83.3%) analyzed through our panel NGS assay revealed genetic alterations other than BRAFV600 (Figure 2A, upper panel). A single mutation was found in 12 samples (33.3%), whereas multiple co-occurring mutations were found in 18 patients (50.0%), with two, three, and four mutations in eight (22.2%), seven (19.4%), and three (8.3%) specimens, respectively. In descending order of frequency, the most common mutations were found in TERT promoter (n = 22, 61.1%), then in NRAS (n = 16, 44.4%), KIT, CDKN2A, and CTNNB1 genes (n = 3 for each, 8.3%) (Figure 2A,B). The predominant TERT promoter mutation in BRAFWT samples was c.-124C>T, detected in 12 (33.3%) cases, whereas c.-146C>T and c.-138/139CC>TT mutations were present in eight (22.2%) and two (5.6%) cases, respectively. Several other mutations were detected at a lower incidence with two mutations in RAC1; two in ERBB4; and one in MET, KRAS, HRAS, SMAD4, PDGFRA, and FGFR2.
Among BRAFV600 specimens, 7 samples (13.2%) exhibited a unique BRAFV600 mutation, while 46 samples (86.8%) had at least one co-occurring mutation (Figure 2A, lower panel). Twenty-six samples (49%) had a single additional mutation to BRAFV600, whereas extra-aberrations on other genes were present in 20 samples (37.7%), with two and three mutations for 15 samples (28.3%) and 5 samples (9.4%), respectively. The prevalence of co-occurring mutations is shown in Figure 2B. TERT promoter was the most frequent genetic alteration in BRAFV600 samples (n = 39, 73.6%) with a predominance of c.-146C>T genotype (n = 20, 73.6%) compared with c.-124C>T genotype (n = 17, 32.1%). CDKN2A and PTEN were commonly mutated, but at a lower frequency, in nine (17.0%) and eight (15.1%) cases, respectively. The remaining co-occurring mutations were detected in NRAS (n = 2), MAP2K1 (n = 2), RAC1 (n = 2), IDH1 (n = 2), PIK3CA (n = 1), STK11 (n = 1), ERBB4 (n = 1), GNAS (n = 1), and GNA11 (n = 1).
When we compared the frequencies of co-occurring mutations variations between the BRAFV600 and BRAFWT patients, we only found a statistically significant difference between these two groups for NRAS (3.8% in BRAFV600 samples versus 44.4% in BRAFWT samples, p < 0.001) (Figure 2B). Of note, some mutations were exclusively detected in BRAFWT samples (i.e., KIT, CTNNB1, MET, KRAS, HRAS, SMAD4, PDGFRA, and FGFR2) or in BRAFV600 samples (i.e., MAP2K1, PIK3CA, STK11, GNAS, and GNA11). Finally, we compared the percentage of mutated genes in BRAFWT and BRAFV600 samples according to the presence of one or several of co-occurring genetic alterations (Figure 2C).

2.3. Correlation between Clinicopathological Features with BRAFV600 Mutation and TERT Promoter Mutation

Then, we tested whether BRAF and TERT promoter mutations were correlated with the clinicopathological features in the cohort analyzed by NGS (Table 1). BRAFV600 melanomas significantly occurred in younger patients, with a median age of 57 years old versus 73 years, (p < 0.001) (Figure S1). Regarding primary tumor sites, their locations differed between the two groups, with 21 (39.6%) BRAFV600 primary tumors located on trunk and 18 (53%) BRAFWT on limbs, including 5 (14.7%) acral melanomas (p = 0.06). Clark level also tends to be lower in BRAFV600 samples (p = 0.09). Conversely, no statistical differences were found for histologic subtypes, Breslow thickness, AJCC (American Joint Committee on Cancer) stage at the diagnosis, and patients’ sex between BRAFV600 and BRAFWT melanomas. For TERT promoter mutational status, a significant association was observed with histological subtypes (p = 0.03) and with primary tumor localization (p = 0.003), but not with age, sex, Breslow thickness, Clark level, nor AJCC stage.

2.4. Prognostic Factors and PFS in BRAFV600 Samples

The prognostic value of clinical parameters and mutational status was next assessed for PFS under targeted therapy in BRAFV600 samples (Table 2, Table S2). All the BRAFV600 patients were treated with MAPK inhibitors: 6 received a monotherapy of BRAF inhibitor (vemurafenib (n = 4) or dabrafenib (n = 2)), 46 received a bitherapy of BRAF and MEK inhibitor (vemurafenib-cobimetinib (n = 4) or dabrafenib-trametinib (n = 42)), and one patient was treated alternately by bitherapy and monotherapy owing to poor tolerance. There were no differences in terms of PFS between the patients treated by monotherapy or bitherapy (p = 0.95). As expected, an elevated level of LDH in patient plasma was significantly associated with a shorter PFS (p = 0.01), while the presence of brain metastasis was near significant (p = 0.09). Regarding the mutational status of the samples, no significant differences were observed in terms of PFS for patients with or without co-occurring aberrations (p = 0.12), or in presence or absence of mutation on TERT promoter, CDKN2A, or PTEN, the most commonly altered genes in our cohort (p = 0.13, 0.34, and 0.35, respectively, Table 2).
However, when TERT mutations were considered separately, a trend was observed (p = 0.059) with a longer PFS for patients harboring a c.-124C>T mutation (Figure 3). When these parameters were entered in a multivariate Cox model, the prognostic significance of the LDH level, brain metastasis, and TERT c.-124C>T mutational status remained in the model (Table 3), demonstrating that these biomarkers are independent prognostic markers of PFS.
To better evaluate the added prognostic value of the TERT c.-124C>T mutation, we further assessed it in subgroups of patients with different predicted outcomes. Interestingly, the presence of a TERT c.-124C>T mutation had a significant prognostic value in patients with brain metastasis (p = 0.01, Figure 4A), but not in patients without (p = 0.51, Figure 4B). The same tendency was also observed in patients with an elevated level of LDH in serum (p = 0.057) compared with patients with a normal LDH level (p = 0.403, Figure S2), suggesting that the presence of a TERT c.-124C>T mutation can represent a marker of good prognosis in subgroups of patients with poor prognosis (presence of brain metastasis and elevated level of LDH).
These observations prompt us to evaluate signatures that combine the TERT mutation status with other prognostic markers. Thus, patients with brain metastasis in absence of TERT c.-124C>T mutation had a worse outcome than other patients (p = 0.003, Figure 5A). This signature was identified as the best fit for predicting worse clinical outcomes (likelihood = 243.55), relative to the model with the brain metastasis status alone or the TERT mutation status alone (likelihood = 248.01 with p = 0.034 and likelihood = 248.81 with p = 0.022, respectively). The same trend was observed for the LDH level parameter, as patients with an elevated LDH level in the absence of TERT c.-124C>T mutation had a worse outcome than others (p = 0.001, Figure 5B), and the signature combining the TERT c.-124C>T mutation status and LDH level had a better fit for predicting worse clinical outcomes (likelihood = 203.45) than models with the LDH level alone or the TERT mutation status alone (likelihood = 206.34 with p = 0.08 and likelihood = 210.95 with p = 0.006, respectively). Finally, the combination of parameters having poor prognosis (e.g., elevated LDH level and/or presence of brain metastasis and absence of TERT c.-124C>T mutation) had a worse outcome in terms of PFS than other patients (p < 0.0001, Figure 5C). Finally, this signature was the best fit for identifying patients of poor prognosis (likelihood = 201.78) compared with LDH level alone (likelihood = 206.34, p = 0.033), brain metastasis alone (likelihood = 209.81, p = 0.005), or TERT mutation status alone (likelihood = 210.96, p = 0.002).

2.5. Prognostic Factors and OS in BRAFV600 Samples

The prognostic value of clinical parameters and mutational status was also evaluated for OS in BRAFV600 samples (Table 4). We observed that sex and presence of brain metastasis were significantly correlated with OS. As for PFS, OS was not differentially affected by the presence or absence of co-occurring genetic alterations in the samples (p = 0.39), or TERT promoter, CDKN2A, or PTEN mutational status (p = 0.15, 0.45, and 0.39, respectively, Table 4), whereas there was a statistically significant difference between the TERT promoter mutation types (p = 0.002) (Figure S3). Thus, the specific effect of the TERT c.-124C>T mutation was retrieved, with longer OS in the subsets of patients with elevated LDH level (p = 0.006) or brain metastasis (p = 0.04), as well as in patients with normal LDH level (p = 0.05) and without brain metastasis (p = 0.005) (Figure S4). As observed for PFS, patients displaying a combination of parameters of poor prognosis (i.e., elevated level of LDH and/or presence of brain metastasis and absence of TERT c.-124C>T mutation) had a worse outcome in terms of OS than other patients (p < 0.0001, Figure 5D).

3. Discussion

Important steps have been taken over the past decade to better characterize the somatic melanoma mutation landscape, thanks to numerous high-throughput sequencing initiatives such as whole genome or whole exome sequencing [8,12,14,15,16]. Thus, the cutaneous melanoma genome has a very high and UV-specific mutational burden. Despite these great advances, the segregation between driver and passenger mutations, the cooperation of mutations, as well as the precise contribution in oncogenesis of each alteration in a same tumor remain elusive. To investigate these crucial considerations, we assessed the mutational frequency of well-known targetable-associated genes in a cohort of 89 cutaneous melanoma samples. We identified high mutational rates within BRAF, NRAS, and TERT promoter, as well as in CDKN2A, PTEN, ERBB4, and RAC1. Among them, mutation of the TERT promoter appeared as a frequent co-occurring mutation in both BRAFV600 and BRAFWT patients. Of the utmost interest, we observed that specific TERT promoter mutation c.-124C>T displayed a statistically significant correlation with MAPK inhibitor treatment efficacy, specifically in a subset of BRAFV600 melanoma patients that had poorer PFS and OS. Importantly, this independent prognosis feature remains true for overall survival.
We observed a BRAFV600 mutation in 53 out of 89 (59.6%) melanomas, with specific BRAFV600E and BRAFV600K mutations in 45 and 8 patients, respectively. When we assessed the associations between BRAF genotyping with clinicopathological features, we found that BRAFV600 mutation was correlated with patient age (p < 0.001) and primary tumor location (p = 0.06), as previously reported [17,18]. In our cohort, no correlation was observed between BRAFV600 mutation status and increasing thickness (p = 0.23) or AJCC stage at the diagnosis (p = 0.56). Apart from its theragnostic implications, the prognosis value of BRAFV600 in melanoma remain controversial; several studies reported association of BRAF mutations with reduced survival [17,19,20], whereas others did not [18,21,22]. This suggests that BRAFV600 mutation alone is an insufficient predictive marker of biological aggressiveness.
We then compared the genomic profiling in BRAFV600 and BRAFWT patients and found that NRAS was differentially mutated between the two groups (p < 0.001). This was expected in the context of the molecular landscape of melanomas that can be divided into four distinct subclasses, including BRAF subtype and NRAS subtype [23]. However, NRAS and BRAF mutation did co-occur in two samples. These two samples belonged to patients with recurrent metastatic melanomas that have been previously treated with MAPK inhibitors, which is consistent with the notion that NRAS mutation can arise through a resistance mechanism following targeted therapy [24]. Among the three remaining samples obtained after targeted therapy, one presented a mutation known to induce secondary resistance to targeted therapy (i.e., a deletion-insertion of MAP2K1), whereas the two others did not. Of note, PTEN and RAC1 mutations—leading to the activation of PI3K/AKT pathway and to a MAPK inhibitor primary resistance [14,16]—were also observed in eight and two BRAFV600 samples, respectively.
No significant differences in terms of PFS were found in the presence or absence of co-occurring aberrations in BRAFV600 samples (p = 0.12). This can probably be explained by the lack of statistical power owing to the limited number of samples harboring a single BRAF mutation (n = 7, 13.2% of the BRAFV600 samples). However, this subject would need further explorations because the prognostic value of such co-existing genetic alterations has been demonstrated in melanoma [9,25] and in other cancers (e.g., in non-small cell lung carcinoma harboring activating EGFR mutations [26,27]).
TERT promoter mutations are currently detected in ~40% of melanomas, with mutation frequencies varying between melanoma subtypes. These mutations are most frequent in melanomas arising in non-acral skin and less frequent in mucosal and acral melanomas. Here, TERT promoter alterations were the most frequent genomic alteration identified. As our cohort was mainly composed of non-acral cutaneous melanomas, TERT promoter mutations occurred in 61 (68.5%) samples, which fitted well with previous reports [27,28,29].
BRAFV600 have been reported to cooperate with TERT promoter mutations more frequently than with any other mutations, upregulating TERT mutant expression via the FOS/GABP pathway [28,30,31,32,33,34,35]. Thus, Vinagre et al. [36] demonstrated the association between BRAF and TERT promoter mutations in melanoma, and that TERT messenger RNA (mRNA) levels are higher when TERT promoter and BRAF mutations coexist in thyroid cancers. In our cohort, the frequency of TERT promoter alterations tended to be higher in BRAFV600 samples, with 39 mutations in 53 samples (73.6%), than in BRAFWT samples, with 22 mutations in 36 samples (61.1%), although this difference did not reach statistical significance. We did not observe any correlations between TERT promoter mutation rates and age, sex, and AJCC stage at the time of diagnosis. Other authors did report significant correlations between TERT promoter mutations, age, and regional and distant metastases [31,37]. This discrepancy may be owing to the small number of cases in our study.
TERT encodes the reverse transcriptase component of the telomerase complex, which is necessary for chromosomal telomeres length stabilization and promotion of cell survival. c.-124C>T and c.-146C>T TERT promoter mutations have been described to create de novo E twenty-six (ETS) binding motifs, inducing upregulation of TERT mRNA and increased telomerase activity in malignant cells. However, TERT expression seems to be differentially enhanced by TERT promoter mutations, with a greater effect of c.-124C>T mutation [38]. This mechanism could be modulated by the presence or absence of the rs2853669 polymorphism at the −245 bp position, which is described to disrupt a pre-existing ETS2 binding site [39]. TERT promoter mutations have been collectively associated with more aggressive melanomas and poorer outcomes, allowing to propose TERT promoter as a poor prognostic factor [28,32]. However, with the exception of the c.-138/-139CC>TT mutation, which is associated with the worst survival in several studies [40,41], the specific effect of each mutation is still debated. Discordant results have been published regarding the c.-146C>T mutation, with a better disease-free survival in comparison with other TERT mutants in one study [41], and a poorer PFS in others [34,40]. In the very recent study of Del Bianco et al. [34], median PFS of patients with c.-124C>T mutation was found to be 9.5 months, versus 5.4 months for patients with c.-146C>T mutation. In our study, the c.-124C>T mutation was associated with better survival compared with the c.-146C>T mutation, but only in poor prognosis patients harboring brain metastasis or elevated serum LDH levels, confirming the different role of both TERT mutations on the MAPK pathway blockade. Interestingly, in the multivariate analysis, c.-124C>T mutation remained significantly associated with brain metastasis (n = 51, p = 0.04) and elevated LDH (n = 46, p = 0.006). The c.-138/-139CC>TT TERT promoter mutation was present in 2 BRAFV600 samples, which was too low to confirm its worse prognostic value. Nevertheless, one of the patients who harbored this mutation was indeed a low responder, with a PFS of 1.8 months, while the second had a PFS of 6.0 months, which was close to the median PFS of BRAFV600 samples (6.1 months). TERT promoter mutations may also have a different impact on patient outcome according to the presence or not of the rs2853669 polymorphism. Indeed, the effect of TERT promoter mutations on survival seems to be enhanced in melanoma patients that did not carry the polymorphism [33]. However, discordant results have been reported in glioblastoma [42,43,44,45]. As the corresponding position was not included in our panel, the link between TERT promoter mutations and this polymorphism and its effect on patient survival could not be unraveled in our study and would require further investigations.
Together, these data provided insights into the recurrent genomic aberrations associated with clinical relevance and demonstrated how these insights could inform triage of patients for effective precision cancer treatments. We observed that MAPK inhibitors treatment efficacy and patients’ PFS, in the presence of elevated LDH level and brain metastasis, depended on the mutation type of the TERT promoter, highlighting, first, the functional link that exists between TERT biology and the MAPK pathway, and second, a specific behavior of the TERT c.-124C>T mutations in a specific subset of poor prognosis melanomas. Finally, there is no direct evidence that TERT promoter aberrations have a key role in melanoma progression, but this genotyping needs to be confirmed in future research. The challenging, but essential task now is to seek confirmation in an independent and prospective collection of additional datasets. If validated in larger cohorts, TERT promoter mutations may be used in clinical practice for its therapeutic relevance, either in terms of influencing the efficacy of established therapies (e.g., MAPK inhibitors or immunotherapies) or they might even prove to be valuable direct therapeutic targets.

4. Materials and Methods

4.1. Patient and Sample Collection

From 2014 to 2019, 113 patients with a cutaneous melanoma underwent NGS analysis at the University Hospital of Montpellier. All primary samples were obtained prior to any treatment initiation, as well as most of the metastasis samples, with the exception of five recurrent metastatic BRAFV600 melanomas that have been previously treated by targeted therapy. The noninterventional study was conducted in accordance with local ethical guidelines and was reviewed by the Ethical Committee from the Montpellier University Hospital (March 2019). Samples were obtained following research project approval by the Institutional Review Board from CRB-CHUM (BB-0033-00031). An approved informed consent statement was acquired for all patients. All corresponding lesions were excised and submitted for standard pathological examination. The percentage of tumor cells in the series ranged from 30 to 100% (Table S2). Tissue punches using a 1 mm needle or macrodissected 10 μm thick section were performed from tumor paraffin blocks to increase the percentage of tumor cells in the sample. Genomic DNA was extracted using the Maxwell® RSC DNA FFPE Kit (Promega, Madison, WI, USA) according to the manufacturer’s recommendations. Twenty-four patients were excluded from the study owing to lost-to-follow-up, absence of targeted therapy for BRAFV600 patients, and/or poor quality of extracted DNA (Figure 1). Medical records were reviewed to extract clinicopathological data, including sex, age, diagnoses, therapeutic agents, and survival.

4.2. NGS Analysis

Library preparation was performed as previously described [46]. Briefly, extracted DNA was qualified using KAPA Sybr® Fast qPCR (Kapa Biosystems, Boston, MA, USA). A home-made panel targeting specific exons of 35 clinically relevant cancer genes was used for amplification of regions of interest (Table S3). For each sample, dual-strand libraries were prepared using a TruSeq Custom Amplicon protocol, as described by the manufacturer (Illumina, Evry, France). After amplification, PCR products were purified using AMPure XP beads (Beckman Coulter, Brea, CA, USA), quantified, normalized, and pair-end sequenced on a MiSeq instrument (2 × 150 cycles, Illumina). After sequencing, the four FastQ files generated per samples were automatically analysed using a bioinformatic workflow managed by Jflow [47]. Briefly, reads were trimmed with cutadapt (version 1.18) [48] and aligned to the human genome GRCh37 with BWA (version 0.7.17) [49], and the variant calling was performed using VarDict (version 1.6.0). Variants present in both libraries with a variant allele frequency >5% and a depth coverage of 300× or greater were then annotated with Variant Effect Predictor (version 94) [50] and reported. Variants having a frequency of 1% or more in the population (in database Exome Aggregation Consortium (ExAC) Variants, Variants Exome Sequencing Project (ESP), or 1000 Genomes Project) were considered as polymorphisms and were excluded.

4.3. Statistical Analysis

The relationship between the patients’ clinicopathological characteristics and other markers was analyzed using the χ2 test or the Fisher test depending on the number of patients per group. Mann–Whitney analysis was performed using Statgraphics (Statgraphics Centurion, Neuilly sur Seine, France). The significance of differences between survival rates was estimated using the Kaplan–Meier method and ascertained with the log-rank test using SPSS® Software (SPSS Inc., Armonk, NY, USA). When the number of samples analyzed was low (n < 20), a bootstrap based log-rank method using 1000 bootstrap constructed samples was used to ascertain the statistical significance of the results. Candidate prognostic factors for PFS with a 0.1 significance level in univariate analysis were entered in a multivariate Cox model, and a backward selection procedure was used to determine independent prognostic markers. A likelihood ratio test was applied to select the best fit between models. A p-value < 0.05 was considered significant.

5. Conclusions

In our study, TERT c.-124C>T mutation was predictive of a better clinical outcome in a subset of poor prognosis BRAFV600 patients under MAPK inhibitors. Altogether, given that adaptive therapeutic strategies are required to overcome acquired resistance in BRAFV600 melanoma patients, we believe that our findings will have relevant clinical implications for patient management and that TERT promoter mutation detection should be considered to better anticipate patient relapse and introduction of immunotherapy. This first line of evidence calls for further explorations in larger cohorts and should be considered in the design and interpretation of future clinical trials.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/8/2224/s1, Figure S1: Age repartition in BRAFV600 and BRAFWT samples, Figure S2: Effect of TERT c.-124C>T promoter mutation on clinical outcome of BRAFV600 patients with elevated or normal LDH level, Figure S3: Univariate analysis of TERT promoter mutations with regard to OS in BRAFV600 samples, Figure S4: Effect of TERT c.-124C>T mutation on OS of BRAFV600 patients, Table S1: Cohort clinicopathological features, Table S2: Clinical prognostic factors and targeted therapy modalities in BRAFV600 samples, Table S3: Sample characteristics, Table S4: NGS panel.

Author Contributions

Conceptualization, J.S.; methodology, P.B., J.A.V., and J.S.; validation, P.B., B.B., and J.A.V.; formal analysis, P.B. and J.A.V.; investigation, V.D., G.C., and P.B.; data curation, P.B.; writing—original draft preparation, P.B.; writing—review and editing, J.A.V., J.S., E.L., and E.C.; supervision, J.S.; project administration, J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Davies, H.; Bignell, G.R.; Cox, C.; Stephens, P.; Edkins, S.; Clegg, S.; Teague, J.; Woffendin, H.; Garnett, M.J.; Bottomley, W.; et al. Mutations of the BRAF gene in human cancer. Nature 2002, 417, 949–954. [Google Scholar] [CrossRef] [PubMed]
  2. Curtin, J.A.; Fridlyand, J.; Kageshita, T.; Patel, H.N.; Busam, K.J.; Kutzner, H.; Cho, K.-H.; Aiba, S.; Bröcker, E.-B.; LeBoit, P.E.; et al. Distinct Sets of Genetic Alterations in Melanoma. Available online: https://www.nejm.org/doi/10.1056/NEJMoa050092?url_ver=Z39.88-2003&rfr_id=ori%3Arid%3Acrossref.org&rfr_dat=cr_pub++0www.ncbi.nlm.nih.gov (accessed on 29 May 2020).
  3. Robert, C.; Grob, J.J.; Stroyakovskiy, D.; Karaszewska, B.; Hauschild, A.; Levchenko, E.; Chiarion Sileni, V.; Schachter, J.; Garbe, C.; Bondarenko, I.; et al. Five-Year Outcomes with Dabrafenib plus Trametinib in Metastatic Melanoma. N. Engl. J. Med. 2019, 381, 626–636. [Google Scholar] [CrossRef] [PubMed]
  4. Atkins, M.B.; Tarhini, A.; Rael, M.; Gupte-Singh, K.; O’Brien, E.; Ritchings, C.; Rao, S.; McDermott, D.F. Comparative efficacy of combination immunotherapy and targeted therapy in the treatment of BRAF-mutant advanced melanoma: A matching-adjusted indirect comparison. Immunotherapy 2019, 11, 617–629. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Long, G.V.; Grob, J.-J.; Nathan, P.; Ribas, A.; Robert, C.; Schadendorf, D.; Lane, S.R.; Mak, C.; Legenne, P.; Flaherty, K.T.; 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, 1743–1754. [Google Scholar] [CrossRef]
  6. Hodis, E.; Watson, I.R.; Kryukov, G.V.; Arold, S.T.; Imielinski, M.; Theurillat, J.-P.; Nickerson, E.; Auclair, D.; Li, L.; Place, C.; et al. A Landscape of Driver Mutations in Melanoma. Cell 2012, 150, 251–263. [Google Scholar] [CrossRef] [Green Version]
  7. Miraflor, A.P.; de Abreu, F.B.; Peterson, J.D.; Turner, S.A.; Amos, C.I.; Tsongalis, G.J.; Yan, S. Somatic mutation analysis in melanoma using targeted next generation sequencing. Exp. Mol. Pathol. 2017, 103, 172–177. [Google Scholar] [CrossRef]
  8. Ticha, I.; Hojny, J.; Michalkova, R.; Kodet, O.; Krkavcova, E.; Hajkova, N.; Nemejcova, K.; Bartu, M.; Jaksa, R.; Dura, M.; et al. A comprehensive evaluation of pathogenic mutations in primary cutaneous melanomas, including the identification of novel loss-of-function variants. Sci. Rep. 2019, 9. [Google Scholar] [CrossRef] [Green Version]
  9. Zhao, X.; Little, P.; Hoyle, A.P.; Pegna, G.J.; Hayward, M.C.; Ivanova, A.; Parker, J.S.; Marron, D.L.; Soloway, M.G.; Jo, H.; et al. The Prognostic Significance of Low-Frequency Somatic Mutations in Metastatic Cutaneous Melanoma. Front. Oncol. 2019, 8. [Google Scholar] [CrossRef]
  10. Zehir, A.; Benayed, R.; Shah, R.H.; Syed, A.; Middha, S.; Kim, H.R.; Srinivasan, P.; Gao, J.; Chakravarty, D.; Devlin, S.M.; et al. Mutational Landscape of Metastatic Cancer Revealed from Prospective Clinical Sequencing of 10,000 Patients. Nat. Med. 2017, 23, 703–713. [Google Scholar] [CrossRef]
  11. Gopal, Y.N.V.; Deng, W.; Woodman, S.E.; Komurov, K.; Ram, P.; Smith, P.D.; Davies, M.A. Basal and treatment-induced activation of AKT mediates resistance to cell death by AZD6244 (ARRY-142886) in Braf-mutant human cutaneous melanoma cells. Cancer Res. 2010, 70, 8736–8747. [Google Scholar] [CrossRef] [Green Version]
  12. Heidorn, S.J.; Milagre, C.; Whittaker, S.; Nourry, A.; Niculescu-Duvas, I.; Dhomen, N.; Hussain, J.; Reis-Filho, J.S.; Springer, C.J.; Pritchard, C.; et al. Kinase-Dead BRAF and Oncogenic RAS Cooperate to Drive Tumor Progression through CRAF. Cell 2010, 140, 209–221. [Google Scholar] [CrossRef] [PubMed]
  13. Wagle, N.; Emery, C.; Berger, M.F.; Davis, M.J.; Sawyer, A.; Pochanard, P.; Kehoe, S.M.; Johannessen, C.M.; MacConaill, L.E.; Hahn, W.C.; et al. Dissecting Therapeutic Resistance to RAF Inhibition in Melanoma by Tumor Genomic Profiling. J. Clin. Oncol. 2011, 29, 3085–3096. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zuo, Q.; Liu, J.; Huang, L.; Qin, Y.; Hawley, T.; Seo, C.; Merlino, G.; Yu, Y. AXL/AKT axis mediated-resistance to BRAF inhibitor depends on PTEN status in melanoma. Oncogene 2018, 37, 3275. [Google Scholar] [CrossRef] [PubMed]
  15. Dutton-Regester, K.; Irwin, D.; Hunt, P.; Aoude, L.G.; Tembe, V.; Pupo, G.M.; Lanagan, C.; Carter, C.D.; O’Connor, L.; O’Rourke, M.; et al. A High-Throughput Panel for Identifying Clinically Relevant Mutation Profiles in Melanoma. Mol. Cancer Ther. 2012, 11, 888–897. [Google Scholar] [CrossRef] [Green Version]
  16. Krauthammer, M.; Kong, Y.; Ha, B.H.; Evans, P.; Bacchiocchi, A.; McCusker, J.P.; Cheng, E.; Davis, M.J.; Goh, G.; Choi, M.; et al. Exome sequencing identifies recurrent somatic RAC1 mutations in melanoma. Nat. Genet. 2012, 44, 1006–1014. [Google Scholar] [CrossRef] [Green Version]
  17. Long, G.V.; Menzies, A.M.; Nagrial, A.M.; Haydu, L.E.; Hamilton, A.L.; Mann, G.J.; Hughes, T.M.; Thompson, J.F.; Scolyer, R.A.; Kefford, R.F. Prognostic and clinicopathologic associations of oncogenic BRAF in metastatic melanoma. J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 2011, 29, 1239–1246. [Google Scholar] [CrossRef]
  18. Ellerhorst, J.A.; Greene, V.R.; Ekmekcioglu, S.; Warneke, C.L.; Johnson, M.M.; Cooke, C.P.; Wang, L.-E.; Prieto, V.G.; Gershenwald, J.E.; Wei, Q.; et al. Clinical Correlates of NRAS and BRAF Mutations in Primary Human Melanoma. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2011, 17, 229–235. [Google Scholar] [CrossRef] [Green Version]
  19. Mann, G.J.; Pupo, G.M.; Campain, A.E.; Carter, C.D.; Schramm, S.-J.; Pianova, S.; Gerega, S.K.; De Silva, C.; Lai, K.; Wilmott, J.S.; et al. BRAF Mutation, NRAS Mutation, and the Absence of an Immune-Related Expressed Gene Profile Predict Poor Outcome in Patients with Stage III Melanoma. J. Invest. Dermatol. 2013, 133, 509–517. [Google Scholar] [CrossRef] [Green Version]
  20. Barbour, A.P.; Tang, Y.H.; Armour, N.; Dutton-Regester, K.; Krause, L.; Loffler, K.A.; Lambie, D.; Burmeister, B.; Thomas, J.; Smithers, B.M.; et al. BRAF mutation status is an independent prognostic factor for resected stage IIIB and IIIC melanoma: Implications for melanoma staging and adjuvant therapy. Eur. J. Cancer 2014, 50, 2668–2676. [Google Scholar] [CrossRef]
  21. Meckbach, D.; Keim, U.; Richter, S.; Leiter, U.; Eigentler, T.K.; Bauer, J.; Pflugfelder, A.; Büttner, P.; Garbe, C.; Weide, B. BRAF-V600 Mutations Have No Prognostic Impact in Stage IV Melanoma Patients Treated with Monochemotherapy. PLoS ONE 2014, 9. [Google Scholar] [CrossRef] [Green Version]
  22. Sperduto, P.W.; Jiang, W.; Brown, P.D.; Braunstein, S.; Sneed, P.; Wattson, D.A.; Shih, H.A.; Bangdiwala, A.; Shanley, R.; Lockney, N.A.; et al. The Prognostic Value of BRAF, C-KIT, and NRAS Mutations in Melanoma Patients with Brain Metastases. Int. J. Radiat. Oncol. Biol. Phys. 2017, 98, 1069–1077. [Google Scholar] [CrossRef] [PubMed]
  23. Akbani, R.; Akdemir, K.C.; Aksoy, B.A.; Albert, M.; Ally, A.; Amin, S.B.; Arachchi, H.; Arora, A.; Auman, J.T.; Ayala, B.; et al. Genomic Classification of Cutaneous Melanoma. Cell 2015, 161, 1681–1696. [Google Scholar] [CrossRef] [PubMed]
  24. Raaijmakers, M.I.G.; Widmer, D.S.; Narechania, A.; Eichhoff, O.; Freiberger, S.N.; Wenzina, J.; Cheng, P.F.; Mihic-Probst, D.; Desalle, R.; Dummer, R.; et al. Co-existence of BRAF and NRAS driver mutations in the same melanoma cells results in heterogeneity of targeted therapy resistance. Oncotarget 2016, 7, 77163–77174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Olbryt, M.; Pigłowski, W.; Rajczykowski, M.; Pfeifer, A.; Student, S.; Fiszer-Kierzkowska, A. Genetic Profiling of Advanced Melanoma: Candidate Mutations for Predicting Sensitivity and Resistance to Targeted Therapy. Target. Oncol. 2020, 15, 101–113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Barnet, M.B.; O’Toole, S.; Horvath, L.G.; Selinger, C.; Yu, B.; Ng, C.C.; Boyer, M.; Cooper, W.A.; Kao, S. EGFR–Co-Mutated Advanced NSCLC and Response to EGFR Tyrosine Kinase Inhibitors. J. Thorac. Oncol. 2017, 12, 585–590. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Deng, L.-L.; Gao, G.; Deng, H.-B.; Wang, F.; Wang, Z.-H.; Yang, Y. Co-occurring genetic alterations predict distant metastasis and poor efficacy of first-line EGFR-TKIs in EGFR-mutant NSCLC. J. Cancer Res. Clin. Oncol. 2019, 145, 2613–2624. [Google Scholar] [CrossRef]
  28. Griewank, K.G.; Murali, R.; Puig-Butille, J.A.; Schilling, B.; Livingstone, E.; Potrony, M.; Carrera, C.; Schimming, T.; Möller, I.; Schwamborn, M.; et al. TERT Promoter Mutation Status as an Independent Prognostic Factor in Cutaneous Melanoma. JNCI J. Natl. Cancer Inst. 2014, 106. [Google Scholar] [CrossRef]
  29. Vazquez, V.d.L.; Vicente, A.L.; Carloni, A.; Berardinelli, G.; Soares, P.; Scapulatempo, C.; Martinho, O.; Reis, R.M. Molecular profiling, including TERT promoter mutations, of acral lentiginous melanomas. Melanoma Res. 2016, 26, 93–99. [Google Scholar] [CrossRef] [Green Version]
  30. Bai, X.; Kong, Y.; Chi, Z.; Sheng, X.; Cui, C.; Wang, X.; Mao, L.; Tang, B.; Li, S.; Lian, B.; et al. MAPK Pathway and TERT Promoter Gene Mutation Pattern and Its Prognostic Value in Melanoma Patients: A Retrospective Study of 2,793 Cases. Clin. Cancer Res. 2017, 23, 6120–6127. [Google Scholar] [CrossRef] [Green Version]
  31. Heidenreich, B.; Nagore, E.; Rachakonda, P.S.; Garcia-Casado, Z.; Requena, C.; Traves, V.; Becker, J.; Soufir, N.; Hemminki, K.; Kumar, R. Telomerase reverse transcriptase promoter mutations in primary cutaneous melanoma. Nat. Commun. 2014, 5, 1–7. [Google Scholar] [CrossRef]
  32. Macerola, E.; Loggini, B.; Giannini, R.; Garavello, G.; Giordano, M.; Proietti, A.; Niccoli, C.; Basolo, F.; Fontanini, G. Coexistence of TERT promoter and BRAF mutations in cutaneous melanoma is associated with more clinicopathological features of aggressiveness. Virchows Arch. 2015, 467, 177–184. [Google Scholar] [CrossRef] [PubMed]
  33. Nagore, E.; Heidenreich, B.; Rachakonda, S.; Garcia-Casado, Z.; Requena, C.; Soriano, V.; Frank, C.; Traves, V.; Quecedo, E.; Sanjuan-Gimenez, J.; et al. TERT promoter mutations in melanoma survival. Int. J. Cancer 2016, 139, 75–84. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Del Bianco, P.; Stagni, C.; Giunco, S.; Fabozzi, A.; Elefanti, L.; Pellegrini, S.; Vecchiato, A.; Pigozzo, J.; Zamuner, C.; De Rossi, A.; et al. TERT Promoter Mutations Differently Correlate with the Clinical Outcome of MAPK Inhibitor-Treated Melanoma Patients. Cancers 2020, 12, 946. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Liu, R.; Zhang, T.; Zhu, G.; Xing, M. Regulation of mutant TERT by BRAF V600E/MAP kinase pathway through FOS/GABP in human cancer. Nat. Commun. 2018, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Vinagre, J.; Almeida, A.; Pópulo, H.; Batista, R.; Lyra, J.; Pinto, V.; Coelho, R.; Celestino, R.; Prazeres, H.; Lima, L.; et al. Frequency of TERT promoter mutations in human cancers. Nat. Commun. 2013, 4, 2185. [Google Scholar] [CrossRef] [Green Version]
  37. Pópulo, H.; Boaventura, P.; Vinagre, J.; Batista, R.; Mendes, A.; Caldas, R.; Pardal, J.; Azevedo, F.; Honavar, M.; Guimarães, I.; et al. TERT Promoter Mutations in Skin Cancer: The Effects of Sun Exposure and X-Irradiation. J. Invest. Dermatol. 2014, 134, 2251–2257. [Google Scholar] [CrossRef] [Green Version]
  38. Heidenreich, B.; Kumar, R. TERT promoter mutations in telomere biology. Mutat. Res. 2017, 771, 15–31. [Google Scholar] [CrossRef]
  39. Hsu, C.P.; Hsu, N.Y.; Lee, L.W.; Ko, J.L. Ets2 binding site single nucleotide polymorphism at the hTERT gene promoter—effect on telomerase expression and telomere length maintenance in non-small cell lung cancer. Eur J Cancer 2006, 42, 1466–1474. [Google Scholar] [CrossRef]
  40. Nagore, E.; Rachakonda, S.; Kumar, R. TERT promoter mutations in melanoma survival. Oncotarget 2019, 10, 1546–1548. [Google Scholar] [CrossRef]
  41. Andrés-Lencina, J.J.; Rachakonda, S.; García-Casado, Z.; Srinivas, N.; Skorokhod, A.; Requena, C.; Soriano, V.; Kumar, R.; Nagore, E. TERT promoter mutation subtypes and survival in stage I and II melanoma patients. Int. J. Cancer 2019, 144, 1027–1036. [Google Scholar] [CrossRef]
  42. Batista, R.; Cruvinel-Carloni, A.; Vinagre, J.; Peixoto, J.; Catarino, T.A.; Campanella, N.C.; Menezes, W.; Becker, A.P.; de Almeida, G.C.; Matsushita, M.M.; et al. The prognostic impact of TERT promoter mutations in glioblastomas is modified by the rs2853669 single nucleotide polymorphism. Int. J. Cancer 2016, 139, 414–423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Simon, M.; Hosen, I.; Gousias, K.; Rachakonda, S.; Heidenreich, B.; Gessi, M.; Schramm, J.; Hemminki, K.; Waha, A.; Kumar, R. TERT promoter mutations: A novel independent prognostic factor in primary glioblastomas. Neuro Oncol. 2015, 17, 45–52. [Google Scholar] [CrossRef] [PubMed]
  44. Spiegl-Kreinecker, S.; Lotsch, D.; Ghanim, B.; Pirker, C.; Mohr, T.; Laaber, M.; Weis, S.; Olschowski, A.; Webersinke, G.; Pichler, J.; et al. Prognostic quality of activating TERT promoter mutations in glioblastoma: Interaction with the rs2853669 polymorphism and patient age at diagnosis. Neuro Oncol. 2015, 17, 1231–1240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Nencha, U.; Rahimian, A.; Giry, M.; Sechi, A.; Mokhtari, K.; Polivka, M.; Schmitt, Y.; Di Stefano, A.L.; Alentorn, A.; Labussiere, M.; et al. TERT promoter mutations and rs2853669 polymorphism: Prognostic impact and interactions with common alterations in glioblastomas. J. Neurooncol. 2016, 126, 441–446. [Google Scholar] [CrossRef]
  46. Vendrell, J.A.; Grand, D.; Rouquette, I.; Costes, V.; Icher, S.; Selves, J.; Larrieux, M.; Barbe, A.; Brousset, P.; Solassol, J. High-throughput detection of clinically targetable alterations using next-generation sequencing. Oncotarget 2017, 8, 40345–40358. [Google Scholar] [CrossRef] [Green Version]
  47. Mariette, J.; Escudié, F.; Bardou, P.; Nabihoudine, I.; Noirot, C.; Trotard, M.-S.; Gaspin, C.; Klopp, C. Jflow: A workflow management system for web applications. Bioinformatics 2016, 32, 456. [Google Scholar] [CrossRef] [Green Version]
  48. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  49. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [Green Version]
  50. McLaren, W.; Gil, L.; Hunt, S.E.; Riat, H.S.; Ritchie, G.R.S.; Thormann, A.; Flicek, P.; Cunningham, F. The Ensembl Variant Effect Predictor. Genome Biol. 2016, 17, 122. [Google Scholar] [CrossRef] [Green Version]
Cancers 12 02224 g001 550
Figure 1. Analytical flowchart of the study. BRAFWT: BRAF wild type, NGS: next generation sequencing.
Figure 1. Analytical flowchart of the study. BRAFWT: BRAF wild type, NGS: next generation sequencing.
Cancers 12 02224 g001
Cancers 12 02224 g002 550
Figure 2. Mutational landscape profiling in BRAFWT and BRAFV600 samples. (A) Co-occurring genetic alterations detected in BRAFWT (n = 36, upper panel) and BRAFV600 samples (n = 53, lower panel). Alteration types are specified (substitution, stop, frameshift, deletion/insertion, or splice variant), except for TERT c.-146C>T, c.-124C>T, or c.-138/139CC>TT mutations. The total number of mutations is shown for each mutated gene in the histogram at the right side of the figure. (B) Frequency of mutated genes in BRAFWT and BRAFV600 samples. *** p < 0.001; # p-values close to significance (p = 0.06). (C) Percentage of mutated genes in BRAFWT (upper panel) and BRAFV600 (lower panel) samples according to the number of co-occurring genetic alterations. Left pie charts show the percentage of mutated genes in samples harboring a single genetic alteration in BRAFWT (n = 12) and BRAFV600 (n = 26). Right pie charts show the percentage of mutated genes in samples harboring several genetic alterations in BRAFWT (n = 18) and BRAFV600 (n = 20).
Figure 2. Mutational landscape profiling in BRAFWT and BRAFV600 samples. (A) Co-occurring genetic alterations detected in BRAFWT (n = 36, upper panel) and BRAFV600 samples (n = 53, lower panel). Alteration types are specified (substitution, stop, frameshift, deletion/insertion, or splice variant), except for TERT c.-146C>T, c.-124C>T, or c.-138/139CC>TT mutations. The total number of mutations is shown for each mutated gene in the histogram at the right side of the figure. (B) Frequency of mutated genes in BRAFWT and BRAFV600 samples. *** p < 0.001; # p-values close to significance (p = 0.06). (C) Percentage of mutated genes in BRAFWT (upper panel) and BRAFV600 (lower panel) samples according to the number of co-occurring genetic alterations. Left pie charts show the percentage of mutated genes in samples harboring a single genetic alteration in BRAFWT (n = 12) and BRAFV600 (n = 26). Right pie charts show the percentage of mutated genes in samples harboring several genetic alterations in BRAFWT (n = 18) and BRAFV600 (n = 20).
Cancers 12 02224 g002
Cancers 12 02224 g003 550
Figure 3. Univariate analysis of TERT promoter mutations with regard to PFS in BRAFV600 samples. Forest plot showing the hazard ratio for PFS associated to TERT promoter mutational status. TERT wild type (WT) samples were taken as reference (ref).
Figure 3. Univariate analysis of TERT promoter mutations with regard to PFS in BRAFV600 samples. Forest plot showing the hazard ratio for PFS associated to TERT promoter mutational status. TERT wild type (WT) samples were taken as reference (ref).
Cancers 12 02224 g003
Cancers 12 02224 g004 550
Figure 4. Effect of TERT c.-124C>T mutation on clinical outcome of BRAFV600 patients with or without brain metastasis. (A) Kaplan–Meier analyses of PFS in patients with brain metastasis in function of the TERT c.-124C>T mutation status. (B) Kaplan–Meier analyses of PFS in patients without brain metastasis in function of the TERT c.-124C>T mutation status.
Figure 4. Effect of TERT c.-124C>T mutation on clinical outcome of BRAFV600 patients with or without brain metastasis. (A) Kaplan–Meier analyses of PFS in patients with brain metastasis in function of the TERT c.-124C>T mutation status. (B) Kaplan–Meier analyses of PFS in patients without brain metastasis in function of the TERT c.-124C>T mutation status.
Cancers 12 02224 g004
Cancers 12 02224 g005 550
Figure 5. Prognostic value of a signature combining parameters of poor prognosis: elevated LDH level, brain metastasis, and absence of TERT c.-124C>T promoter mutation. (A) Kaplan–Meier analyses of PFS in patients combining brain metastasis and absence of TERT c.-124C>T promoter mutation (group A, green curve) and in the remaining patients (group B, blue curve). (B) Kaplan–Meier analyses of PFS in patients combining elevated LDH level and absence of the TERT c.-124C>T promoter mutation (group A, green curve) and in the remaining patients (group B, blue curve). (C) Kaplan–Meier analyses of PFS in patients combining parameters of poor prognosis (group A, green curve) and in the remaining patients (group B, blue curve). (D) Kaplan–Meier analyses of overall survival (OS) in patients combining parameters of poor prognosis (group A, green curve) and in the remaining patients (group B, blue curve).
Figure 5. Prognostic value of a signature combining parameters of poor prognosis: elevated LDH level, brain metastasis, and absence of TERT c.-124C>T promoter mutation. (A) Kaplan–Meier analyses of PFS in patients combining brain metastasis and absence of TERT c.-124C>T promoter mutation (group A, green curve) and in the remaining patients (group B, blue curve). (B) Kaplan–Meier analyses of PFS in patients combining elevated LDH level and absence of the TERT c.-124C>T promoter mutation (group A, green curve) and in the remaining patients (group B, blue curve). (C) Kaplan–Meier analyses of PFS in patients combining parameters of poor prognosis (group A, green curve) and in the remaining patients (group B, blue curve). (D) Kaplan–Meier analyses of overall survival (OS) in patients combining parameters of poor prognosis (group A, green curve) and in the remaining patients (group B, blue curve).
Cancers 12 02224 g005
Table
Table 1. Correlation between clinicopathological features and BRAF and TERT promoter mutational status.
Table 1. Correlation between clinicopathological features and BRAF and TERT promoter mutational status.
Clinicopathological FeaturesBRAF StatusTERT Promoter Status
BRAFWT (n = 36)BRAFV600
(n = 53)		p-ValueWT (n = 28)Mutated (n = 61)p-Value
Age≤604 (11.1)28 (52.8)<0.0019 (32.1)23 (37.7)0.61
>6032 (88.9)25 (47.2)19 (67.9)38 (62.3)
SexMale22 (61.1)28 (52.8)0.4415 (53.6)35 (57.4)0.74
Female14 (38.9)25 (47.2)13 (46.4)26 (42.6)
Histological typeNM8 (25.8)10 (20.8)0.185 (20)13 (24.1)0.03
SSM15 (48.4)24 (50)10 (40)29 (53.7)
MUP3 (9.7)9 (18.8)3 (12)9 (16.7)
Unclassified0 (0)3 (6.3)1 (4)2 (3.7)
ALM5 (16.1)2 (4.2)6 (24)1 (1.9)
Missing data5537
Primary tumor siteHead/neck7 (20.6)11 (20.8)0.061 (3.7)17 (28.3)0.003
Upper limbs4 (11.8)2 (3.8)1 (3.7)5 (8.3)
Trunk6 (17.6)21 (39.6)8 (29.6)19 (31.7)
Lower limbs9 (26.5)8 (15.1)8 (29.6)9 (15.0)
MUP3 (8.8)9 (17)3 (11.1)9 (15.0)
Acral5 (14.7)2 (3.8)6 (22.2)1 (1.7)
Missing data2011
Breslow thickness *<1 mm1 (3.3)4 (10.2)0.252 (8.7)3 (6.5)0.66
1–1.99 mm4 (13.3)11 (28.2)3 (13.0)12 (26.1)
2–3.99 mm8 (26.6)9 (23.1)6 (26.1)11 (23.9)
≥4 mm17 (56.6)15 (38.4)12 (52.2)20 (43.5)
Missing data3526
Clark level *II0 (0)2 (5.4)0.090 (0)2 (4.6)0.66
III2 (6.9)10 (27.0)3 (13.6)9 (20.4)
IV21 (72.4)20 (54.1)15 (68.2)26 (59.1)
V6 (20.7)5 (13.5)4 (18.2)7 (15.9)
Missing data4738
AJCCI3 (8.8)7 (14.0)0.563 (11.1)7 (12.3)0.97
II16 (47.1)17 (34.0)10 (37.0)23 (40.4)
III8 (23.5)11 (22.0)6 (22.2)13 (22.8)
IV7 (20.6)15 (30.0)8 (29.6)14 (24.6)
Missing data2314
WT: wild-type, NM: nodular melanoma, SSM: superficial spreading melanoma, MUP: melanoma of unknown primary, ALM: acral lentiginous melanoma, AJCC: American Joint Committee on Cancer. * MUP are not taken into account for Breslow thickness and Clark level.
Table
Table 2. Univariate analysis of clinical parameters and sample mutational status with regard to progression-free survival (PFS) in BRAFV600 patients.
Table 2. Univariate analysis of clinical parameters and sample mutational status with regard to progression-free survival (PFS) in BRAFV600 patients.
Clinical ParametersSamples (n)HR95% CIp-Value
Sex (Male; Female)511.160.61–2.20NS (0.65)
Age (≤60; >60-year-old)511.660.88–3.14NS (0.12)
LDH level (normal; high)462.381.19–4.750.01
Brain metastases (absence; presence)511.720.90–3.29NS (0.09)
Histological type (NM; SSM; MUP; ALM; unclassified)451.350.93–1.95NS (0.11)
Breslow thickness (<1 mm; 1–2 mm; 2–4 mm; >4)371.100.75–1.62NS (0.61)
Treatment modalities (monotherapy; bitherapy)501.030.40–2.65NS (0.95)
Co-occurring mutation (absence; presence)510.490.20–1.20NS (0.12)
TERT promoter mutation (absence; presence)510.580.29–1.17NS (0.13)
CDKN2A mutation (absence; presence)511.500.65–3.44NS (0.34)
PTEN mutation (absence; presence)511.780.65–3.40NS (0.35)
HR: hazard ratio, 95% CI: 95% confidence interval, NM: nodular melanoma, SSM: superficial spreading melanoma, MUP: melanoma of unknown primary, ALM: acral lentiginous melanoma, LDH: lactate dehydrogenase, NS: not significant.
Table
Table 3. Multivariate analysis of clinical parameters and TERT c.-124C>T mutational status with regard to PFS in BRAFV600 patients.
Table 3. Multivariate analysis of clinical parameters and TERT c.-124C>T mutational status with regard to PFS in BRAFV600 patients.
Clinical ParametersHR95% CIp-Value
LDH level (normal; high)2.861.36–6.020.006
Brain metastasis (absence; presence)2.341.03–5.300.04
TERT c.-124C>T mutation (absence; presence)1.371.08–1.730.009
HR: hazard ratio, 95%CI: 95% confidence interval, LDH: lactate dehydrogenase.
Table
Table 4. Univariate analysis of clinical parameters and sample mutational status with regard to overall survival (OS) in BRAFV600 patients.
Table 4. Univariate analysis of clinical parameters and sample mutational status with regard to overall survival (OS) in BRAFV600 patients.
Clinical ParametersSamples (n)HR95% CIp-Value
Sex (Male; Female)512.201.10–4.420.03
Age (≤60; >60-year-old)511.630.83–3.20NS (0.15)
LDH level (normal; high)461.710.83–3.52NS (0.15)
Brain metastases (absence; presence)512.211.08–4.500.03
Histological type (NM; SSM; MUP; ALM; unclassified)451.270.89–1.82NS (0.19)
Breslow thickness (<1 mm; 1–2 mm; 2–4 mm; >4)371.020.67–1.55NS (0.93)
Treatment modalities (monotherapy; bitherapy)500.610.25–1.49NS (0.27)
Co-occurring mutation (absence; presence)510.680.28–1.65NS (0.39)
TERT promoter mutation (absence; presence)510.580.28–1.21NS (0.15)
CDKN2A mutation (absence; presence)511.410.58–3.45NS (0.45)
PTEN mutation (absence; presence)511.440.63–3.31NS (0.39)
HR: hazard ratio, 95% CI: 95% confidence interval, NM: nodular melanoma, SSM: superficial spreading melanoma, MUP: melanoma of unknown primary, ALM: acral lentiginous melanoma.

Share and Cite

MDPI and ACS Style

Blateau, P.; Coyaud, E.; Laurent, E.; Béganton, B.; Ducros, V.; Chauchard, G.; Vendrell, J.A.; Solassol, J. TERT Promoter Mutation as an Independent Prognostic Marker for Poor Prognosis MAPK Inhibitors-Treated Melanoma. Cancers 2020, 12, 2224. https://doi.org/10.3390/cancers12082224

AMA Style

Blateau P, Coyaud E, Laurent E, Béganton B, Ducros V, Chauchard G, Vendrell JA, Solassol J. TERT Promoter Mutation as an Independent Prognostic Marker for Poor Prognosis MAPK Inhibitors-Treated Melanoma. Cancers. 2020; 12(8):2224. https://doi.org/10.3390/cancers12082224

Chicago/Turabian Style

Blateau, Pauline, Etienne Coyaud, Estelle Laurent, Benoit Béganton, Vincent Ducros, Géraldine Chauchard, Julie A. Vendrell, and Jérôme Solassol. 2020. "TERT Promoter Mutation as an Independent Prognostic Marker for Poor Prognosis MAPK Inhibitors-Treated Melanoma" Cancers 12, no. 8: 2224. https://doi.org/10.3390/cancers12082224

APA Style

Blateau, P., Coyaud, E., Laurent, E., Béganton, B., Ducros, V., Chauchard, G., Vendrell, J. A., & Solassol, J. (2020). TERT Promoter Mutation as an Independent Prognostic Marker for Poor Prognosis MAPK Inhibitors-Treated Melanoma. Cancers, 12(8), 2224. https://doi.org/10.3390/cancers12082224

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