Next Article in Journal
Molecular Classification and Tumor Microenvironment Characterization of Gallbladder Cancer by Comprehensive Genomic and Transcriptomic Analysis
Next Article in Special Issue
Brain Metastases from Colorectal Cancer: A Systematic Review of the Literature and Meta-Analysis to Establish a Guideline for Daily Treatment
Previous Article in Journal
Blocking Migration of Polymorphonuclear Myeloid-Derived Suppressor Cells Inhibits Mouse Melanoma Progression
Previous Article in Special Issue
Brain Metastases from Uterine Cervical and Endometrial Cancer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 4
10.3390/cancers13040731
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

Integrative Genomic Analyses of Patient-Matched Intracranial and Extracranial Metastases Reveal a Novel Brain-Specific Landscape of Genetic Variants in Driver Genes of Malignant Melanoma

by
Renáta Váraljai
1,2,†,
Susanne Horn
1,2,3,†,
Antje Sucker
1,2,
Daniela Piercianek
2,4,
Verena Schmitt
5,
Alexander Carpinteiro
6,
Katrin Anne Becker
6,
Julia Reifenberger
7,
Alexander Roesch
1,2,
Jörg Felsberg
8,
Guido Reifenberger
2,8,
Ulrich Sure
2,4,
Dirk Schadendorf
1,2 and
Iris Helfrich
1,2,*
1
Skin Cancer Unit of the Dermatology Department, Medical Faculty, West German Cancer Center, University Duisburg-Essen, 45147 Essen, Germany
2
German Cancer Consortium (DKTK), partner site Essen/Düsseldorf, 45147 Essen, Germany
3
Faculty Rudolf-Schönheimer-Institute for Biochemistry, University of Leipzig, 04103 Leipzig, Germany
4
Department of Neurosurgery, Medical Faculty, West German Cancer Center, University Duisburg-Essen, 45147 Essen, Germany
5
Institute of Anatomy, Medical Faculty, University Duisburg-Essen, 45147 Essen, Germany
6
Department of Molecular Biology, Medical Faculty, University Duisburg-Essen, 45147 Essen, Germany
7
Department of Dermatology, Medical Faculty, Heinrich Heine University, 40225 Düsseldorf, Germany
8
Institute of Neuropathology, Heinrich Heine University, 40225 Düsseldorf, Germany
*
Author to whom correspondence should be addressed.
†
Equally contributing authors.
Cancers 2021, 13(4), 731; https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13040731
Submission received: 5 January 2021 / Revised: 28 January 2021 / Accepted: 8 February 2021 / Published: 10 February 2021
(This article belongs to the Special Issue Brain Metastases (Secondary Brain Tumor))
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Melanoma is the third most common cause of brain metastasis with a reported incidence of up to 80% leading to patients’ early mortality. Clinical activity at intracranial sites is often less and unsatisfactory when compared to extracranial metastases by using novel targeted or immune therapies. Thus, the identification of genetic alterations may provide new insights into the pathogenesis of brain metastases and this will facilitate the improvement of precision oncology. Therefore, the aim of our study was to address site-specific oncogenic alterations in intracranial metastases of 29 recurrently mutated driver genes in melanoma by next generation sequencing. In line with the branched evolution model of metastasis, we identified in our cohort of intracranial and corresponding patient-matched extracranial melanoma metastases novel genetic variants and site-specific nucleotide modifications. Therapeutic targeting of the new-identified genetic variants could help to facilitate novel, more effective therapies for prevention and/or treatment of melanoma brain metastases.

Abstract

Background: Development of brain metastases in advanced melanoma patients is a frequent event that limits patients’ quality of life and survival. Despite recent insights into melanoma genetics, systematic analyses of genetic alterations in melanoma brain metastasis formation are lacking. Moreover, whether brain metastases harbor distinct genetic alterations beyond those observed at different anatomic sites of the same patient remains unknown. Experimental Design and Results: In our study, 54 intracranial and 18 corresponding extracranial melanoma metastases were analyzed for mutations using targeted next generation sequencing of 29 recurrently mutated driver genes in melanoma. In 11 of 16 paired samples, we detected nucleotide modifications in brain metastases that were absent in matched metastases at extracranial sites. Moreover, we identified novel genetic variants in ARID1A, ARID2, SMARCA4 and BAP1, genes that have not been linked to brain metastases before; albeit most frequent mutations were found in ARID1A, ARID2 and BRAF. Conclusion: Our data provide new insights into the genetic landscape of intracranial melanoma metastases supporting a branched evolution model of metastasis formation.

1. Introduction

Metastasis of malignant melanoma to the brain is a clinically challenging issue that may develop in up to 40% of patients with advanced disease [1]. The incidence of brain metastases (BM) is rising partly due to improved diagnostic techniques and advances in systemic treatment approaches directing prolonged survival of cancer patients. However, therapeutic strategies specifically targeting the metastatic cascade demonstrating survival benefit are still lacking. In consequence, metastatic spread is responsible for about 90% of cancer-related deaths across all entities [2]. Different cohort studies demonstrated cutaneous melanoma as the third most common cause of BM development [3]. BM in malignant melanoma is frequent during melanoma progression, dominating prognosis and quality of life of affected patients [4,5,6]. The incidence of BM at first presentation is about 20%, in patients suffering from advanced melanoma around 50%, and even higher as autopsy studies reported frequencies of 55 up to 80% [7,8]. Prognosis of these patients is particularly poor, resulting in median overall survival of only 17 to 22 weeks [9,10]. In consequence, the presence of BM was recently incorporated into the American Joint Committee on Cancer (AJCC) staging system as an independent prognostic factor in patients with malignant melanoma [11]. Treatment options targeting established metastases in the central nervous system (CNS) are rather limited, mainly caused by inefficient drug penetration across the blood-brain barrier (BBB). Moreover, patients with BM are commonly excluded from clinical trials, including those investigating novel targeted therapies, as the limited survival associated with BM prevents reaching study endpoints.
Large-scale sequencing studies have revealed the widespread genetic diversity of melanoma, directing therapy efficacy and patient’s survival [12,13,14]. Thus, in the current era of molecularly targeted therapies and personalized medicine, therapeutic decisions are increasingly being tailored according to the individual genetic profile of a tumor. Primary melanomas and extracranial metastases have been extensively studied, but largely due to lack of available tissue, the biology of intracranial melanoma metastases remains poorly understood. However, genomic sequencing of primary tumors and intracranial metastases from different types of solid cancers revealed remarkable genetic heterogeneity but also further genetic modifications from primary tumors to metastases at distant sites, or BM specific mutations, implicating that intracerebral metastases develop from “branched genetic evolution” [15,16,17,18].
Several large-scale sequencing studies in melanoma identified significantly mutated melanoma genes, such as NF), ARID2, TERT or RAC beside the well-established oncogenes including NRAS, BRAF or KIT [13,19,20,21]. About half of the melanoma patients harbor activating mutation in oncogenes involved in the mitogen-activated protein kinase (MAPK) pathway, e.g., BRAF or NRAS and/or constitutive activation of the cell cycle modulating phosphoinositide 3-kinase (PI3K) signaling pathway. Combining BRAF with MEK inhibitors shows meaningful clinical efficacy in terms of objective response rates in extracranial melanoma [22] and quite recently, also in intracranial activity with response rates of 44 to 59% in patients with melanoma BM. The median intracranial duration of response was 6.5 month in a prospective study in melanoma patients [23]. However, integrative genetic analyses focusing on oncogenic alterations of the “hotspot” cancer genes of anatomically and temporally distinct melanoma metastases from the same patient are still limited. These analyses are of pivotal clinical importance for the understanding of the metastatic seeding of cancer cells to the brain, but also to guide the development of next generation targeted therapies.
The aim of this study was to identify BM-associated driver mutations in melanoma with potential for translation into future therapeutic strategies. Therefore, we performed targeted next-generation sequencing of 29 melanoma-associated genes in intracerebral melanoma metastases and “matched” melanoma metastases from extracranial sites.

2. Results

2.1. Patient Characteristics

We included 47 patients, 20 females and 27 males, diagnosed with malignant cutaneous melanoma. The mean age at the time of first diagnosis of malignant melanoma was 60 ± 14 and 55 ± 14 (mean ± SD), respectively (incomplete clinical data from 14 patients). The time of first intracranial tumor presentation occurred at 56 ± 14 (mean ± SD) years in the whole cohort. From 16 of the 47 patients matched biopsies were available from intracerebral and extracranial metastases allowing for comparative genetic analyses. In total, we included 54 melanoma metastases with cerebellar (5 cases), cerebral (48 cases) or spinal (1 case) localization. The set of 18 extracranial metastases included 16 cutaneous, 1 lymph node and 1 adrenal gland melanoma metastases. Relevant clinical data of these patients are listed in Table 1.

2.2. Comparison of Brain Metastases and Metastases of Other Sites

We called single nucleotide variants including known mutations and rare single-nucleotide-polymorphisms (SNPs) in 47 of 54 intracerebral and in 15 of 18 extracranial metastases (Figure 1). To evaluate co-occurrence and mutual exclusivity of the detected variants, we tested sample pairs between the 29 genes as shown on the OncoPrint. We evaluated genes that were altered in at least 10% of the samples (n ≥ 7). SMARCA4 mutations significantly co-occurred with ARID2 (p = 0.001, log2 odds ratio [OR] >3, n = 9) and ARID1A (p < 0.0001, log2 OR > 3, n = 7) mutations. Similarly, TERT mutations significantly co-occurred with NRAS and MAP2K2 (both: p = 0.046, log2 OR = 2.373, n = 9), RAC1 (p = 0.004, log2 OR > 3, n = 8), and NF1 (p = 0.035, log2 OR > 3, n = 7) mutations. ARID2 was mutated in samples that also carried ARID1A (p < 0.0001, log2 OR > 3, n = 8) or NF1 mutations (p < 0.0001, log2 OR > 3, n = 8). Finally, TP53 mutations were highly likely present in NRAS mutated samples (p = 0.005, log2 OR > 3, n = 8).
The 16 matched biopsies allowed for the direct comparison of anatomically and temporally distinct melanoma metastases. In 11 of the 16 pairs, we detected genetic variants in BM that were absent in the corresponding metastases obtained from other sites of the body. These included mutations in 12 well-known cancer genes associated with cutaneous and/or uveal melanomas: BRAF, NRAS, PIK3CA, TERT, IDH1, EZH2, NF1, PTEN, TP53, MAP2K1, GNAQ and GNA11. Interestingly, we identified variants in five new candidate genes, i.e., ARID1A, ARID2, SMARCA4, PIK3R1 and BAP1, which neither have been previously identified in melanoma BM nor reported in CNS metastases from other cancers (Figure 2A). We assessed the progression patterns of matched extracranial and BM samples to determine which gene variants are significantly associated with BM within our matched sample dataset. Variants in 6 genes showed a statistically significant occurrence specific to BM. Mutations in ARID1A (p = 0.0205), ARID2 (p = 0.0207), BRAF (p = 0.023) SMARCA4 (p = 0.015), TERT (p = 0.0041), and IDH1 (p = 0.0041) were significantly exclusive to BM while the same variants remained undetected in the corresponding extracranial tumor tissues.
Next, we analyzed the mutation frequency in the mutational driver genes. Overall, the most frequently mutated genes in BM included ARID1A, ARID2, BRAF, and SMARCA4 (Figure 2B, Table 2). ARID1A variants (P650S, P1568S, L2119S, P1562L, G831A, G423E4) and ARID2 variants (D239N, S1489L, P1022S, S297F) were detected in 4 intracerebral melanoma tissues. BRAF variants (P336L, G73E and V207E) and SMARCA4 variants (P262L, P919S, E1113K, H884Y) were detected in 3 melanoma BM. We also detected nucleotide transitions resulting in gene activation in NRAS (Q61K, Q61R), in the promoter region of TERT (chr5:1295250 G > A), and in GNA11 (R183C), the latter being frequent in uveal melanoma [24] and primary meningeal melanocytic tumors [25]. Moreover, recurrence was seen for three brain-specific variants, in BRAF (V207E), and TP53 (E171G) in 2 of 16 tested pairs. These genetic variants include a typical cancer driver mutations in BRAF. Brain-specific loss-of-function mutations were detected in PIK3CA (Y1157*) and in PTEN (Q219*, A328fs). Previously not described variants in PIK3R1 (S43R, F41C) were detected private to BM. The detailed description of these variants are reported in Table 2. Furthermore, we evaluated whether some mutations tend to co-occurred in BM. SMARCA4 significantly co-occurred with IDH1 and PIK3CA mutations (both: p = 0.042, log2 OR > 3, n = 2). As seen before, SMARCA4 tended to co-occurred with ARID1A and ARID2 mutations, however, here significance was not reached (p = 0.083 and p = 0.583, respectively). IDH1 mutations were highly likely present in samples that also carried PIK3CA mutations (p = 0.042, log2 OR > 3, n = 2). We also searched in 180 different cancer sequencing studies (including over 47000 samples using the cBioPortal site) for evidence of BM associated variants, and while all of the detected variants in the 17 genes were found in different cancer types, only 4 variants (NRAS Q61, TP53 E171G, PTEN A328*, ARID2 S297F) were detected corresponding to CNS cancers (n = 10 tumor tissues). The presence of these mutations in primary CNS tumors further support their significance in brain metastases.

2.3. UV Mutation Signature

Metastases located in sun-exposed areas of the skin might accumulate more mutations compared to metastases located of sun-protected areas of the body, such as metastases of the CNS. Indeed, using our multi-gene amplicon panel the number of called SNVs in BM was significantly lower compared to metastases at other sites (p = 0.003, Figure 3).
In addition, large-scale genome sequencing indicated that considerable elevated baseline mutation rates in cutaneous melanoma are primarily represented by increased abundance of cytidine to thymidine (C > T or CC > TT) transitions that are characteristic for an UV-light-induced mutational signature [13,26,27]. As most of the metastases at extracranial sites were indeed cutaneous metastases (16/18), we analyzed the frequency of C > T transitions to address confounding impact of high mutational load due to UV mutagenesis. As expected, BM showed less frequent (C > T) variants than extracranial metastases (p = 0.038, Figure 4A). Next, we analyzed C > T transitions in four mutational subtypes of melanoma, BRAF-, NRAS-, NF1-mutant samples and triple-wildtype samples [28]. As reported, NF1-mutant samples showed a significantly higher abundance of C > T transitions when compared to samples from the other three subgroups (p = 0.004, Figure 4B).

3. Discussion

Melanoma is the third common cause of intracerebral metastases after lung and breast carcinoma [29]. In recent years, treatment of melanoma has made remarkable progress in the era of precision medicine and targeted oncology. However, treatment of BM in melanoma patients remains one of the major clinical challenges. Although recent data show therapeutic activity against melanoma BM by targeting the BRAF pathway or by immune checkpoint inhibition using inhibitors against the cytotoxic T-lymphocyte-associated Protein 4 (CTLA-4) or the Programmed cell death protein 1 (PD-1), response rates are rather limited and, overall, patients still have a poor prognosis [30,31]. Intracranial response rates for BRAF inhibitors in patients harboring mutant BRAF (V600E) have been reported in the range of 30–50%, with an increased response in asymptomatic patients and patients without previous local BM therapy [10,23,32]. However, while stable disease was achieved in the majority of patients with partial or complete response, CNS progression has been reported in almost 75% of patients undergoing BRAF-targeted therapy [33,34].
Furthermore, recent data indicate activity of anti-CTLA-4 [35] and anti-PD-1 monotherapy [36,37], or a combination of both, resulting in intracranial response rates of up to 47 %; however, the level of intracranial response was not translated into improved overall survival [38]. Thus, there is a time-critical clinical need to understand the biology of BM and their molecular signatures to guide the development of novel therapeutic approaches against metastases of the CNS.
The divergence of genetic alterations detected in multiple clinical sequencing projects including primary tumors and “patient-matched” BM raised the notion, that each distinct lesion might harbor a unique set of oncogenic alterations. Collecting patient-matched tissue sets is a major challenge mainly due to the fact that resection of BM is performed by neurosurgeons while primary melanomas and cutaneous melanoma metastases are resected by dermatologists, often in private practice. Nevertheless, in this study we were able to screen 54 melanoma brain metastases of 47 patients and 16 matched extracranial melanoma metastases collected from 16 corresponding patients to focus on genetic variants in anatomically and temporally distinct metastases. We focused on 29 genes, identified for their cancer driving capability in cutaneous or uveal melanoma. In line with the literature, we found a reduced mutation load in tumor samples from intracranial localization, including less frequent UV-light-induced mutational signatures when compared to extracranial, mostly cutaneous metastases. In a recently published study by Fischer et al., [7] the relevance of mutations were also evaluated in patient matched extracranial and BM samples from 17 patients across 74 genes. While only 16 genes overlapped between our gene panel and the gene panel by Fischer et al., in 5 of those genes (BRAF, GNAQ, NF1, BAP1, PIK3CA, PIK3R1) mutations private to BM were confirmed. Although that the exact nucleotide transitions within those genes could not be directly compared to the variants disclosed in our study, the clinical significance of distinct genetic BM alterations is further supported. Notably, 56% of our “matched-pair” metastases presented genetic variants private to BM. Moreover, we detected single nucleotide transitions in ARID1A, ARID2, SMARCA4 and BAP1, all genes which have not been associated before with CNS metastases of melanomas. However, beside modifications in ARID1A, genetic variants in BRAF were detected most frequently compared to alterations in other genes. Point mutations in BRAF, mostly resulting by the substitution of valine for glutamin acid at codon 600 (VAL600 Glu or V600E) but also NRAS are the most frequent activating somatic events in melanoma and have been extensively studied in large cohorts of clinical trials for their impact on disease progression and therapeutic intervention [13,23,38,39,40,41]. In the metastatic situation, mutations in both genes were implicated in worse survival outcomes and increased risk for development of BM [9;10]. In line with previous studies of extracranial melanomas, none of our 16 sample pairs presented coexistence of BRAF and NRAS mutations. However, the clinical correlates of BRAF and NRAS mutations in melanoma BM are limited. Fang and colleagues recently described the relationship between tumor BRAF/NRAS mutation status, clinical characteristics and response to conventional therapy in patients with melanoma BM [42]. The authors resumed that either an activating BRAF mutation or an NRAS mutation significantly increases the local failure rate compared to patients who carried tumors that were wildtype for both genes.
ARID1A encodes the ATPase subunit interacting component BAF250, which, together with a set of core subunits, represents the switch/sucrose non-fermentable (SWI/SNF) chromatin remodeling machinery, providing access of proteins to DNA [43]. ARID1A mutations were first reported in ovarian clear cell carcinoma and subsequently described in endometrium-derived carcinomas [44,45]. However, growing evidence indicated that ARID1A may have a widespread suppressive role in various cancer entities [46,47]. In line with this finding, it has been reported that ARID1A negatively regulates TERT transcriptional regulation and activity via binding to TERT regulatory elements. In consequence, cells showed a repressive histone mode and a survival advantage via telomere maintenance [48]. Moreover, recent data reported by Shen and colleagues correlate loss-of-function mutations in ARID1A with the efficacy of anti-PD1 therapy, as a result of enhanced intratumoral lymphocyte recruitment and up-regulation of PD-L1 in preclinical models of ovarian cancer [49]. However, the clinical significance of such differential expressions and the function of the ARID1A protein remain undefined due to the lack of studies using fresh human tumor samples. In addition to ARID1A mutations, our data also indicated variants of unknown significance in other components of the SWI/SNF complex, such as ARID2 and SMARCA4. Single nucleotide modifications in ARID2 and SMARCA4, all of unknown function but together with mutated splice site and frameshifts events often summarized and reviewed as “loss-of-function mutations”, have already been reported in extracranial melanomas [13,50,51,52]. In summary, 11/54 (20%, Figure 2A) and 6 of the 11 patient-matched brain pairs (55%, Figure 2B) of our discovery samples harbored genetic modifications in a component of the SWI/SNF machinery, implicating a potential role for dysregulation of chromatin remodeling in promoting the formation of melanoma CNS metastases. Interestingly, SMARCA4, first identified to be highly mutated in NSCLC patients [53], can act as a corepressor of ZEB1, promoting epithelial-to-mesenchymal transition (EMT) [54]. EMT is a major signature of highly tumorigenic cancer cells, inducing E-cadherin mediated disruption of cell-cell contacts and, consequently, metastatic dissemination of cancer cells to distant organs [55]. In line with this finding, RNASeq data of spontaneously growing brain metastases of our novel preclinical melanoma models but also human BM biopsies clearly indicate strong EMT signatures when compared to RNASeq profiles of extracerebral melanoma metastases).
In 2 of 16 paired patient samples, we found mutations of unknown relevance in BAP1 which were absent in the corresponding extracranial metastases. In both cases, we detected single nucleotide changes A > T or C > T in position 52436829 (protein change V295M and L650H). Germline pathogenic variants and somatic mutations in the BAP1 gene have been described in ocular and cutaneous melanoma and paralleled with a highly aggressive ocular melanoma phenotype [56,57] but also detected in several other cancer subtypes [58].
In summary, our data support the concept that melanoma metastasis formation is frequently driven by branched-evolution with unique gene modifications restricted to intracranial metastases and absent in corresponding metastases of distant sites. However, larger cohorts of matched primary and metastatic tumors at different sites as well as large-scale genomic sequencing are needed to comprehensively elucidate the landscape of genetic alterations in melanoma brain metastases and uncover their clinical impact.

4. Materials and Methods

4.1. Patient Material

Patients suffering from advanced skin melanoma were included in this study. Samples of intracranial melanoma metastases were provided by the tissue banks at the Department of Neurosurgery, University Hospital Essen, as well as the Institute of Neuropathology, Heinrich Heine University Düsseldorf, Germany. Extracranial metastases of corresponding patients (“matched-pair” samples), including cutaneous, lymph node and adrenal gland melanoma metastases, were retrieved from the Skin Cancer Biobank (SCABIO) of the Department of Dermatology, University Hospital Essen, or the Department of Dermatology, Heinrich Heine University Düsseldorf, Germany. Intracranial and extracranial melanoma metastases were histopathologically diagnosed by our in-house neuropathologists (for intracranial metastases) and dermatopathologists (for extracranial metastases). Informed patient consent and the appropriate IRB approval was obtained for all patients. The study was performed with approval by the ethics committee of the Medical Faculty, University Duisburg-Essen (ethics approvals no. 11-4715 and no. 15-6723-BO), and the ethics committee of the Medical Faculty, Heinrich Heine University Düsseldorf (ethics approval no. 5246).

4.2. DNA Isolation

Sections of 10 µm thickness (3–4 sections per sample) were cut from fresh frozen samples of CNS metastases or cutaneous, lymph node and adrenal gland metastases of patients diagnosed with malignant melanoma. Prior DNA extraction standard hematoxylin and eosin (H&E) staining was performed for visualization of the tissue morphology. After histopathological confirmation, the tumor area was marked as “Region of Interest (ROI)” by our inhouse pathologists and manually dissected from the slide. Samples from intracranial and extracranial origin presenting a tumor content of >80% in the ROI have been included in the analysis. Genomic DNA was isolated by using the QIAmp DNA Mini Kit (Qiagen, Hilden, Germany), including RNase A treatment, according to the manufacturer’s instructions.

4.3. Targeted Sequencing

A custom designed amplicon-based sequencing panel covering 29 genes known as being recurrently mutated in cutaneous and uveal melanomas was applied [59]. Sequencing libraries were prepared applying the GeneRead Library Prep Kit from Qiagen® according to the manufacturer’s instructions. For adapter ligation and barcoding samples, we applied the NEBNext Ultra DNA Library Prep Mastermix Set and NEBNext Multiplex Oligos for Illumina from New England Biolabs. Twelve to 24 samples were sequenced in parallel on an Illumina MiSeq next generation sequencer.
Sequence analysis was performed with CLC Cancer Research Workbench (Qiagen®). In brief, the following steps were applied. The CLC workflow included adapter trimming and merging paired reads before mapping to the human reference genome (hg19). Insertions and deletions as well as single nucleotide variants were detected, local realignment and primer trimming followed. Additional information was then obtained regarding potential mutation type, single nucleotide polymorphisms with a minor allele frequency of 1 percent and higher, and conservation scores by cross-referencing databases (COSMIC, ClinVar, dbSNP, 1000 Genomes Project, HAPMAP and PhastCons-Conservation_scores_hg19). After CLC processing, csv files were analyzed manually. Mutations affecting protein coding regions of genes were regarded if predicted to result in non-synonymous amino acid changes. Questionable low frequency background mutations calls, not uncommon in FFPE amplicon sequencing approaches were excluded [60] by applying the following thresholds: an overall coverage of the mutation site ≥ 30 reads, ≥15 reads reporting the mutated variant and ≥5% frequency of mutated reads. Genetic variants private to BM were reported only when the target region was covered sufficiently as described above in the corresponding metastasis of extracranial origin.

4.4. Data Analysis

To confirm validity of the analyzed genomic variants, we checked for pathogenic score in ClinVar [61] and in COSMIC databases [62]. Genomic variants were furthermore searched in 180 different cancer sequencing studies including over 47000 samples using the cBioPortal website [63]. OncoPrint maps were generated using the OncoPrinter tool [64,65]. Sankey diagram was created using the SankeyMATIC tool. Two-sided unpaired t-test was used for the determination of statistical significance between experimental groups. Mutational co-occurrence or mutual exclusivity was tested by two-sided Fisher’s exact test with multiple test correction. Mutational association with organ site (extracranial vs. BM vs. both) was analyzed by two-sided Chi-square test. Chi-square p value was reported only for those genes where a valid comparison could be made i.e., all groups (extracranial vs. BM vs. both) contained values > 0. p values ≤ 0.05 were considered as statistically significant.

5. Conclusions

Melanoma brain metastases demonstrate site-specific genetic variants (62%). Mutations in ARID1A, ARID2, BRAF, SMARCA4, TERT and IDH1 were significantly exclusive to intracranial metastases while the same variants remained undetected in the corresponding extracranial tumor tissues. We also detected previously not described genetic variants in the “hot spot” oncogenic driver genes private to melanoma brain metastases. Analyses of large cancer sequencing data sets of primary CNS tumors point the significance of our novel variants for tumor establishment in the brain.

Author Contributions

I.H., S.H. and A.S. conceptualized, designed and performed the MiSeq experiments, R.V. and A.R. analyzed the data and generated statistical analyses. D.P., V.S., J.R., J.F., G.R., U.S. and D.S. provided clinical material and designed the selection of the patient cohort. V.S., A.C., and K.A.B. performed experiments and generated data. D.P., A.R. and A.S. collected and provided clinical data. I.H. conceptualized, coordinated and directed the project. I.H., R.V., A.R. and S.H. wrote the manuscript. Manuscript review was performed by all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the German Cancer Aid (DKH, no.70112507 to I.H., D.S. and G.R., and no.111402 to S.H.), Deutsche Forschungsgemeinschaft (DFG KFO 337 to I.H.: HE 5294/2-1; SH: HO 6389/2-1; AR: RO 3577/7-1; DS: SCHA 422/17-1), Hiege Stiftung gegen Hautkrebs, Monika Kutzner Stiftung (P#18) and Brigitte & Konstanze Wegener Stiftung (all to I.H.).

Institutional Review Board Statement

Informed patient consent and the appropriate IRB approval was obtained for all patients. The study was performed with approval by the ethics committee of the Medical Faculty, University Duisburg-Essen (ethics approvals no. 11-4715 and no. 15-6723-BO), and the ethics committee of the Medical Faculty, Heinrich Heine University Düsseldorf (ethics ap-proval no. 5246).

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be shared upon request.

Acknowledgments

We thank M. Eilebrecht and I.-V. Westedt (Department of Dermatology, Medical Faculty Essen) for technical assistance. We gratefully acknowledge Kathy Keyvani and her team in Institute of Neuropathology, University Hospital Essen for providing the histological diagnoses of the cases in Essen. We also acknowledge support by the Open Access Publication Fund of the University of Duisburg-Essen.

Conflicts of Interest

G. Reifenberger has received research grants from Roche and Merck, as well.

Abbreviations

AJCC: American Joint Committee on Cancer; ARID1: AT-rich interactive domain-containing protein 1; ARID2: AT-rich interactive domain-containing protein 2; BBB: blood-brain barrier; BM: brain metastases; BRAF: v-Raf murine sarcoma viral oncogene homolog B; CNS: Central Nervous System; CTLA-4: cytotoxic T-lymphocyte-associated Protein 4; EMT: epithelial-to-mesencymal transition; EZH2: Enhancer of zeste homolog 2; GNA11: G Protein Subunit Alpha 11; GNAQ: Guanine nucleotide-binding protein G(q) subunit alpha; IDH1: Isocitrate dehydrogenase 1; KIT: cluster of differentiation 117; MAPK: mitogen-activated protein kinase; MEK: Mitogen-activated protein kinase kinase; NF1: Neurofibromin 1; NRAS: Neuroblastoma RAS viral oncogene homolog; NSCLC: Non-small cell lung cancer; PI3K: phosphoinisitide 3-kinase; PD-1: Programmed cell death protein 1; PIK3CA: Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Alpha; PTEN: Phosphatase and tensin homolog deleted in from chromosome ten; RAC: Ras-related C3 botulinum toxin substrate; SD: Standard deviation; SWI/SNF: switch/sucrose non-fermentable; TERT: Telomerase reverse transcriptase, ZEB1: Zinc Finger E-Box Binding Homeobox 1.

References

  1. Cagney, D.N.; Martin, A.M.; Catalano, P.J.; Redig, A.J.; Lin, N.U.; Lee, E.Q.; Wen, P.Y.; Dunn, I.F.; Bi, W.L.; Weiss, E.S.; et al. Incidence and prognosis of patients with brain metastases at diagnosis of systemic malignancy: A population-based study. Neuro-Oncology 2017, 19, 1511–1521. [Google Scholar] [CrossRef] [Green Version]
  2. Chaffer, C.L.; Weinberg, R.A. A Perspective on Cancer Cell Metastasis. Science 2011, 331, 1559–1564. [Google Scholar] [CrossRef]
  3. Fabi, A.; Felici, A.; Metro, G.; Mirri, A.; Bria, E.; Telera, S.; Moscetti, L.; Russillo, M.; Lanzetta, G.; Mansueto, G.; et al. Brain metastases from solid tumors: Disease outcome according to type of treatment and therapeutic resources of the treating center. J. Exp. Clin. Cancer Res. 2011, 30, 10. [Google Scholar] [CrossRef] [Green Version]
  4. Balch, C.M.; Balch, G.C.; Sharma, R.R. Identifying Early Melanomas at Higher Risk for Metastases. J. Clin. Oncol. 2012, 30, 1406–1407. [Google Scholar] [CrossRef] [PubMed]
  5. Gramsch, C.; Göricke, S.L.; Behrens, F.; Zimmer, L.; Schadendorf, D.; Krasny, A.; Forsting, M.; Schlamann, M.U. Isolated cerebral susceptibility artefacts in patients with malignant melanoma: Metastasis or not? Eur. Radiol. 2013, 23, 2622–2627. [Google Scholar] [CrossRef]
  6. Miller, D.; Zappala, V.; El Hindy, N.; Livingstone, E.; Schadendorf, D.; Sure, U.; Sandalcioglu, I.E. Intracerebral metastases of malignant melanoma and their recurrences—A clinical analysis. Clin. Neurol. Neurosurg. 2013, 115, 1721–1728. [Google Scholar] [CrossRef]
  7. Fischer, G.M.; Jalali, A.; Kircher, D.A.; Lee, W.-C.; McQuade, J.L.; Haydu, L.E.; Joon, A.Y.; Reuben, A.; De Macedo, M.P.; Carapeto, F.C.L.; et al. Molecular Profiling Reveals Unique Immune and Metabolic Features of Melanoma Brain Metastases. Cancer Discov. 2019, 9, 628–645. [Google Scholar] [CrossRef] [Green Version]
  8. Gorantla, V.; Kirkwood, J.M.; Tawbi, H.A. Melanoma Brain Metastases: An Unmet Challenge in the Era of Active Therapy. Curr. Oncol. Rep. 2013, 15, 483–491. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Falchook, G.S.; Long, G.V.; Kurzrock, R.; Kim, K.B.; Arkenau, T.H.; Brown, M.P.; Hamid, O.; Infante, J.R.; Millward, M.; Pavlick, A.C.; et al. Dabrafenib in patients with melanoma, untreated brain metastases, and other solid tumours: A phase 1 dose-escalation trial. Lancet 2012, 379, 1893–1901. [Google Scholar] [CrossRef] [Green Version]
  10. Long, G.V.; Trefzer, U.; Davies, M.A.; Kefford, R.F.; Ascierto, P.A.; Chapman, P.B.; Puzanov, I.; Hauschild, A.; Robert, C.; Algazi, A.; et al. Dabrafenib in patients with Val600Glu or Val600Lys BRAF-mutant melanoma metastatic to the brain (BREAK-MB): A multicentre, open-label, phase 2 trial. Lancet Oncol. 2012, 13, 1087–1095. [Google Scholar] [CrossRef]
  11. Gershenwald, J.E.; Scolyer, R.A.; Hess, K.R.; Sondak, V.K.; Long, G.V.; Ross, M.I.; Lazar, A.J.; Faries, M.B.; Kirkwood, J.M.; McArthur, G.A.; et al. Melanoma staging: Evidence-based changes in the American Joint Committee on Cancer eighth edition cancer staging manual. CA Cancer J. Clin. 2017, 67, 472–492. [Google Scholar] [CrossRef] [Green Version]
  12. Hayward, N.K.; Wilmott, J.S.; Waddell, N.; Johansson, P.A.; Field, M.A.; Nones, K.; Patch, A.-M.; Kakavand, H.; Alexandrov, L.B.; Burke, H.; et al. Whole-genome landscapes of major melanoma subtypes. Nat. Cell Biol. 2017, 545, 175–180. [Google Scholar] [CrossRef]
  13. 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]
  14. Krauthammer, M.; Kong, Y.; Bacchiocchi, A.; Evans, P.; Pornputtapong, N.; Wu, C.; McCusker, J.P.; Ma, S.; Cheng, E.; Straub, R.; et al. Exome sequencing identifies recurrent mutations in NF1 and RASopathy genes in sun-exposed melanomas. Nat. Genet. 2015, 47, 996–1002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Brastianos, P.K.; Carter, S.L.; Santagata, S.; Cahill, D.P.; Taylor-Weiner, A.; Jones, R.T.; Van Allen, E.M.; Lawrence, M.S.; Horowitz, P.M.; Cibulskis, K.; et al. Genomic Characterization of Brain Metastases Reveals Branched Evolution and Potential Therapeutic Targets. Cancer Discov. 2015, 5, 1164–1177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Chen, G.; Chakravarti, N.; Aardalen, K.; Lazar, A.J.; Tetzlaff, M.T.; Wubbenhorst, B.; Kim, S.-B.; Kopetz, S.; LeDoux, A.A.; Gopal, Y.V.; et al. Molecular Profiling of Patient-Matched Brain and Extracranial Melanoma Metastases Implicates the PI3K Pathway as a Therapeutic Target. Clin. Cancer Res. 2014, 20, 5537–5546. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Liu, W.; Laitinen, S.; Khan, S.; Vihinen, M.; Kowalski, J.; Yu, G.; Chen, L.; Ewing, C.M.; Eisenberger, M.A.; Carducci, M.A.; et al. Copy number analysis indicates monoclonal origin of lethal metastatic prostate cancer. Nat. Med. 2009, 15, 559–565. [Google Scholar] [CrossRef] [Green Version]
  18. Lu, Y.-W.; Zhang, H.-F.; Liang, R.; Xie, Z.-R.; Luo, H.-Y.; Zeng, Y.-J.; Xu, Y.; Wang, L.-M.; Kong, X.-Y.; Wang, K.-H. Colorectal Cancer Genetic Heterogeneity Delineated by Multi-Region Sequencing. PLoS ONE 2016, 11, e0152673. [Google Scholar] [CrossRef] [Green Version]
  19. Horn, S.; Figl, A.; Rachakonda, P.S.; Fischer, C.; Sucker, A.; Gast, A.; Kadel, S.; Moll, I.; Nagore, E.; Hemminki, K.; et al. TERT Promoter Mutations in Familial and Sporadic Melanoma. Science 2013, 339, 959–961. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. 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]
  21. Van Allen, E.M.; Wagle, N.; Sucker, A.; Treacy, D.J.; Johannessen, C.M.; Goetz, E.M.; Place, C.S.; Taylor-Weiner, A.; Whittaker, S.; Kryukov, G.V.; et al. The genetic landscape of clinical resistance to RAF inhibition in metastatic melanoma. Cancer Discov. 2014, 4, 94–109. [Google Scholar] [CrossRef] [Green Version]
  22. Robert, C.; Karaszewska, B.; Schachter, J.; Rutkowski, P.; Mackiewicz, A.; Stroiakovski, D.; Lichinitser, M.; Dummer, R.; Grange, F.; Mortier, L.; et al. Improved Overall Survival in Melanoma with Combined Dabrafenib and Trametinib. N. Engl. J. Med. 2015, 372, 30–39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Davies, M.A.; Saiag, P.; Robert, C.; Grob, J.J.; Flaherty, K.T.; Arance, A.; Chiarion-Sileni, V.; Thomas, L.; Lesimple, T.; Mortier, L.; et al. Dabrafenib plus trametinib in patients with BRAF(V600)-mutant melanoma brain metastases (COMBI-MB): A multicentre, multicohort, open-label, phase 2 trial. Lancet Oncol. 2017, 18, 863–873. [Google Scholar] [CrossRef]
  24. Van Raamsdonk, C.D.; Griewank, K.G.; Crosby, M.B.; Garrido, M.C.; Vemula, S.; Wiesner, T.; Obenauf, A.C.; Wackernagel, W.; Green, G.; Bouvier, N.; et al. Mutations in GNA11 in uveal melanoma. N. Engl. J. Med. 2010, 363, 2191–2199. [Google Scholar] [CrossRef] [Green Version]
  25. Van de Nes, J.; Wrede, K.; Ringelstein, A.; Stiller, M.; Horn, S.; Sucker, A.; Möller, I.; Scholz, S.L.; Murali, R.; Gessi, M.; et al. Diagnosing a Primary Leptomeningeal Melanoma by Gene Mutation Signature. J. Investig. Dermatol. 2016, 136, 1526–1528. [Google Scholar] [CrossRef] [Green Version]
  26. Berger, M.F.; Hodis, E.; Heffernan, T.P.; Deribe, Y.L.; Lawrence, M.S.; Protopopov, A.; Ivanova, E.; Watson, I.R.; Nickerson, E.; Ghosh, P.; et al. Melanoma genome sequencing reveals frequent PREX2 mutations. Nat. Cell Biol. 2012, 485, 502–506. [Google Scholar] [CrossRef] [PubMed]
  27. Pleasance, E.D.; Cheetham, R.K.; Stephens, P.J.; McBride, D.J.; Humphray, S.J.; Greenman, C.; Varela, I.; Lin, M.-L.; Ordóñez, G.R.; Bignell, G.R.; et al. A comprehensive catalogue of somatic mutations from a human cancer genome. Nat. Cell Biol. 2009, 463, 191–196. [Google Scholar] [CrossRef]
  28. Shain, A.H.; Garrido, M.C.; Botton, T.; Talevich, E.; Yeh, I.; Sanborn, J.Z.; Chung, J.; Wang, N.J.; Kakavand, H.; Mann, G.J.; et al. Exome sequencing of desmoplastic melanoma identifies recurrent NFKBIE promoter mutations and diverse activating mutations in the MAPK pathway. Nat. Genet. 2015, 47, 1194–1199. [Google Scholar] [CrossRef] [PubMed]
  29. Brastianos, P.K.; Curry, W.T.; Oh, K.S. Clinical Discussion and Review of the Management of Brain Metastases. J. Natl. Compr. Cancer Netw. 2013, 11, 1153–1164. [Google Scholar] [CrossRef] [Green Version]
  30. Milsch, L.; Gesierich, A.; Kreft, S.; Livingstone, E.; Zimmer, L.; Goebeler, M.; Schadendorf, D.; Schilling, B. Patterns of disease control and survival in patients with melanoma brain metastases undergoing immune-checkpoint blockade. Eur. J. Cancer 2018, 99, 58–65. [Google Scholar] [CrossRef]
  31. Sloot, S.; Chen, Y.A.; Zhao, X.; Weber, J.S.; Benedict, J.J.; Mulé, J.J.; Smalley, K.S.; Weber, J.S.; Zager, J.S.; Forsyth, P.A.; et al. Improved survival of patients with melanoma brain metastases in the era of targeted BRAF and immune checkpoint therapies. Cancer 2018, 124, 297–305. [Google Scholar] [CrossRef] [PubMed]
  32. Rastrelli, M.; Pigozzo, J.; Di Maggio, A.; Tosi, A.L.; Sileni, V.C.; Rossi, C.R. Neoadjuvant treatment with dabrafenib of unresectable localizations from occult melanoma. Melanoma Res. 2014, 24, 413–414. [Google Scholar] [CrossRef] [PubMed]
  33. Peuvrel, L.; Saint-Jean, M.; Quéreux, G.; Brocard, A.; Khammari, A.; Knol, A.C.; Dréno, B. Incidence and characteristics of melanoma brain metastases developing during treatment with vemurafenib. J. Neuro-Oncology 2014, 120, 147–154. [Google Scholar] [CrossRef]
  34. Puzanov, I.; Amaravadi, R.K.; McArthur, G.A.; Flaherty, K.T.; Chapman, P.B.; Sosman, J.A.; Ribas, A.; Shackleton, M.; Hwu, P.; Chmielowski, B.; et al. Long-term outcome in BRAF(V600E) melanoma patients treated with vemurafenib: Patterns of disease progression and clinical management of limited progression. Eur. J. Cancer 2015, 51, 1435–1443. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Margolin, K.A.; Ernstoff, M.S.; Hamid, O.; Lawrence, D.; McDermott, D.; Puzanov, I.; Wolchok, J.D.; Clark, J.I.; Sznol, M.; Logan, T.F.; et al. Ipilimumab in patients with melanoma and brain metastases: An open-label, phase 2 trial. Lancet Oncol. 2012, 13, 459–465. [Google Scholar] [CrossRef]
  36. Goldberg, S.B.; Gettinger, S.N.; Mahajan, A.; Chiang, A.C.; Herbst, R.S.; Sznol, M.; Tsiouris, A.J.; Cohen, J.; Vortmeyer, A.; Jilaveanu, L.; et al. Pembrolizumab for patients with melanoma or non-small-cell lung cancer and untreated brain metastases: Early analysis of a non-randomised, open-label, phase 2 trial. Lancet Oncol. 2016, 17, 976–983. [Google Scholar] [CrossRef] [Green Version]
  37. Nishino, M.; Giobbie-Hurder, A.; Manos, M.P.; Bailey, N.D.; Buchbinder, E.I.; Ott, P.A.; Ramaiya, N.H.; Hodi, F.S. Immune-Related Tumor Response Dynamics in Melanoma Patients Treated with Pembrolizumab: Identifying Markers for Clinical Outcome and Treatment Decisions. Clin. Cancer Res. 2017, 23, 4671–4679. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Long, G.V.; Atkinson, V.; Lo, S.; Sandhu, S.; Guminski, A.D.; Brown, M.P.; Wilmott, J.S.; Edwards, J.; Gonzalez, M.A.; Scolyer, R.; et al. Combination nivolumab and ipilimumab or nivolumab alone in melanoma brain metastases: A multicentre randomised phase 2 study. Lancet Oncol. 2018, 19, 672–681. [Google Scholar] [CrossRef]
  39. Chapman, P.; Hauschild, A.; Robert, C.; Haanen, J.B.; Ascierto, P.; Larkin, J.; Dummer, R.; Garbe, C.; Testori, A.; Maio, M.; et al. Improved Survival with Vemurafenib in Melanoma with BRAF V600E Mutation. N. Engl. J. Med. 2011, 364, 2507–2516. [Google Scholar] [CrossRef] [Green Version]
  40. Flaherty, K.T.; Robert, C.; Hersey, P.; Nathan, P.; Garbe, C.; Milhem, M.; Demidov, L.V.; Hassel, J.C.; Rutkowski, P.; Mohr, P.; et al. Improved Survival with MEK Inhibition in BRAF-Mutated Melanoma. N. Engl. J. Med. 2012, 367, 107–114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Goyal, S.; Silk, A.W.; Tian, S.; Mehnert, J.; Danish, S.; Ranjan, S.; Kaufman, H.L. Clinical Management of Multiple Melanoma Brain Metastases: A Systematic Review. JAMA Oncol. 2015, 1, 668–676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Fang, P.; Boehling, N.S.; Koay, E.J.; Bucheit, A.D.; Jakob, J.A.; Settle, S.H.; Brown, P.D.; Davies, M.A.; Sulman, E.P. Melanoma brain metastases harboring BRAF V600K or NRAS mutations are associated with an increased local failure rate following conventional therapy. J. Neuro-Oncology 2018, 137, 67–75. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, X.; Nagl, N.G.; Jr Flowers, S.; Zweitzig, D.; Dallas, P.B.; Moran, E. Expression of p270 (ARID1A), a component of human SWI/SNF complexes, in human tumors. Int. J. Cancer 2004, 112, 636. [Google Scholar] [CrossRef] [PubMed]
  44. Jones, S.; Wang, T.-L.; Shih, I.-M.; Mao, T.-L.; Nakayama, K.; Roden, R.; Glas, R.; Slamon, D.; Diaz, L.A.; Vogelstein, B.; et al. Frequent Mutations of Chromatin Remodeling Gene ARID1A in Ovarian Clear Cell Carcinoma. Science 2010, 330, 228–231. [Google Scholar] [CrossRef] [Green Version]
  45. Wiegand, K.C.; Sy, K.; Kalloger, S.E.; Li-Chang, H.; Woods, R.; Kumar, A.; Streutker, C.J.; Hafezi-Bakhtiari, S.; Zhou, C.; Lim, H.J.; et al. ARID1A/BAF250a as a prognostic marker for gastric carcinoma: A study of 2 cohorts. Hum. Pathol. 2014, 45, 1258–1268. [Google Scholar] [CrossRef] [PubMed]
  46. Wilson, B.G.; Roberts, C.W.M. SWI/SNF nucleosome remodellers and cancer. Nat. Rev. Cancer 2011, 11, 481–492. [Google Scholar] [CrossRef]
  47. Wu, R.; Wang, T.-L.; Shih, I.-M. The emerging roles of ARID1A in tumor suppression. Cancer Biol. Ther. 2014, 15, 655–664. [Google Scholar] [CrossRef] [Green Version]
  48. Rahmanto, Y.S.; Jung, J.-G.; Wu, R.-C.; Kobayashi, Y.; Heaphy, C.M.; Meeker, A.K.; Wang, T.-L.; Shih, I.-M. Inactivating ARID1A Tumor Suppressor Enhances TERT Transcription and Maintains Telomere Length in Cancer Cells. J. Biol. Chem. 2016, 291, 9690–9699. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Shen, J.; Ju, Z.; Zhao, W.; Wang, L.; Peng, Y.; Ge, Z.; Nagel, Z.D.; Zou, J.; Wang, C.; Kapoor, P.; et al. ARID1A deficiency promotes mutability and potentiates therapeutic antitumor immunity unleashed by immune checkpoint blockade. Nat. Med. 2018, 24, 556–562. [Google Scholar] [CrossRef]
  50. Nikolaev, S.I.; Rimoldi, D.; Iseli, C.; Valsesia, A.; Robyr, D.; Gehrig, C.; Harshman, K.; Guipponi, M.; Bukach, O.; Zoete, V.; et al. Exome sequencing identifies recurrent somatic MAP2K1 and MAP2K2 mutations in melanoma. Nat. Genet. 2011, 44, 133–139. [Google Scholar] [CrossRef]
  51. Stark, M.S.; Woods, S.L.; Gartside, M.G.; Bonazzi, V.F.; Dutton-Regester, K.; Aoude, L.G.; Chow, D.; Sereduk, C.; Niemi, N.M.; Tang, N.; et al. Frequent somatic mutations in MAP3K5 and MAP3K9 in metastatic melanoma identified by exome sequencing. Nat. Genet. 2011, 44, 165–169. [Google Scholar] [CrossRef]
  52. Wei, X.; Program, N.C.S.; Walia, V.; Lin, J.C.; Teer, J.K.; Prickett, T.D.; Gartner, J.; Davis, S.; Stemke-Hale, K.A.; Davies, M.; et al. Exome sequencing identifies GRIN2A as frequently mutated in melanoma. Nat. Genet. 2011, 43, 442–446. [Google Scholar] [CrossRef] [Green Version]
  53. Medina, P.P.; Carretero, J.; Fraga, M.F.; Esteller, M.; Sidransky, D.; Sanchez-Cespedes, M. Genetic and Epigenetic screening for gene alterations of the chromatin-remodeling factor, SMARCA4/BRG1, in lung tumors. Genes Chromosom. Cancer 2004, 41, 170–177. [Google Scholar] [CrossRef]
  54. Sanchez-Tillo, E.; Lazaro, A.; Torrent, R.; Cuatrecasas, M.; Vaquero, E.C.; Castells, A.; Engel, P.; Postigo, A. ZEB1 represses E-cadherin and induces an EMT by recruiting the SWI/SNF chromatin-remodeling protein BRG1. Oncogene 2010, 29, 3490–3500. [Google Scholar] [CrossRef] [Green Version]
  55. Brabletz, T.; Kalluri, R.; Nieto, M.A.; Weinberg, R.A. EMT in cancer. Nat. Rev. Cancer 2018, 18, 128–134. [Google Scholar] [CrossRef]
  56. Njauw, C.-N.J.; Kim, I.; Piris, A.; Gabree, M.; Taylor, M.; Lane, A.M.; DeAngelis, M.M.; Gragoudas, E.; Duncan, L.M.; Tsao, H. Germline BAP1 Inactivation Is Preferentially Associated with Metastatic Ocular Melanoma and Cutaneous-Ocular Melanoma Families. PLoS ONE 2012, 7, e35295. [Google Scholar] [CrossRef] [Green Version]
  57. Wiesner, T.; Obenauf, A.C.; Murali, R.; Fried, I.; Griewank, K.G.; Ulz, P.; Windpassinger, C.; Wackernagel, W.; Loy, S.; Wolf, I.; et al. Germline mutations in BAP1 predispose to melanocytic tumors. Nat. Genet. 2011, 43, 1018–1021. [Google Scholar] [CrossRef] [Green Version]
  58. Tesch, M.E.; Pater, J.A.; Vandekerkhove, G.; Wang, G.; Binnington, K.; So, A.I.; Wyatt, A.W.; Eigl, B.J. Concurrent germline and somatic pathogenic BAP1 variants in a patient with metastatic bladder cancer. NPJ Genom. Med. 2020, 5, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Van de Nes, J.; Gessi, M.; Sucker, A.; Moller, I.; Stiller, M.; Horn, S.; Scholz, S.L.; Pischler, C.; Stadtler, N.; Schilling, B.; et al. Targeted next generation sequencing reveals unique mutation profile of primary melanocytic tumors of the central nervous system. J. Neurooncol. 2016, 127, 435–444. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Murali, R.; Chandramohan, R.; Möller, I.; Scholz, S.L.; Berger, M.; Huberman, K.; Viale, A.; Pirun, M.; Socci, N.D.; Bouvier, N.; et al. Targeted massively parallel sequencing of angiosarcomas reveals frequent activation of the mitogen activated protein kinase pathway. Oncotarget 2015, 6, 36041–36052. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Landrum, M.J.; Lee, J.M.; Benson, M.; Brown, G.R.; Chao, C.; Chitipiralla, S.; Gu, B.; Hart, J.; Hoffman, D.; Jang, W.; et al. ClinVar: Improving access to variant interpretations and supporting evidence. Nucleic Acids Res. 2018, 46, D1062–D1067. [Google Scholar] [CrossRef] [Green Version]
  62. Tate, J.G.; Bamford, S.; Jubb, H.C.; Sondka, Z.; Beare, D.M.; Bindal, N.; Boutselakis, H.; Cole, C.G.; Creatore, C.; Dawson, E.; et al. COSMIC: The Catalogue of Somatic Mutations in Cancer. Nucleic Acids Res. 2018, 47, D941–D947. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. cBioPortal. Available online: http://www.cbioportal.org/ (accessed on 5 October 2020).
  64. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data: Figure 1. Cancer Discov. 2012, 2, 401–404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal. Sci. Signal. 2013, 6, pl1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cancers 13 00731 g001 550
Figure 1. OncoPrint chart displaying the type of detected genomic alterations in the 29 gene-panel across 72 sequenced melanoma metastases. Patient identification (ID) numbers are given. Each column represents a sample, sample type is indicated at the top in a color-coded manner. The percentage of samples with detectable alteration is listed on the left for each gene.
Figure 1. OncoPrint chart displaying the type of detected genomic alterations in the 29 gene-panel across 72 sequenced melanoma metastases. Patient identification (ID) numbers are given. Each column represents a sample, sample type is indicated at the top in a color-coded manner. The percentage of samples with detectable alteration is listed on the left for each gene.
Cancers 13 00731 g001
Cancers 13 00731 g002 550
Figure 2. (A) Sankey diagram illustrating the shared and unique variants in 20 genes across 11 matched samples (left: extracranial, right: brain tumor tissues). Variants that were detected in extracranial, but not in brain metastases are not shown. Asterisk indicate p ≤ 0.05 significance from Chi-square test, where the proportion of samples with shared and unique variants were tested (extracranial vs. BM vs. both). (B) OncoPrint chart displaying the frequency of detectable alterations in 17 genes that were specific to brain tissues.
Figure 2. (A) Sankey diagram illustrating the shared and unique variants in 20 genes across 11 matched samples (left: extracranial, right: brain tumor tissues). Variants that were detected in extracranial, but not in brain metastases are not shown. Asterisk indicate p ≤ 0.05 significance from Chi-square test, where the proportion of samples with shared and unique variants were tested (extracranial vs. BM vs. both). (B) OncoPrint chart displaying the frequency of detectable alterations in 17 genes that were specific to brain tissues.
Cancers 13 00731 g002
Cancers 13 00731 g003 550
Figure 3. Box and whiskers plot showing the sum of called variants in extracranial and brain meta-static tissues including unpaired samples. The number of detected variants were summed and plotted in extracranial metastases (mean: 17, median: 4, range: 1–83) and in brain metastases (mean: 4, median: 2, range: 1–12). p-value indicates two-sided unpaired t-test.
Figure 3. Box and whiskers plot showing the sum of called variants in extracranial and brain meta-static tissues including unpaired samples. The number of detected variants were summed and plotted in extracranial metastases (mean: 17, median: 4, range: 1–83) and in brain metastases (mean: 4, median: 2, range: 1–12). p-value indicates two-sided unpaired t-test.
Cancers 13 00731 g003
Cancers 13 00731 g004 550
Figure 4. (A) Box and whiskers plot showing the frequency of C > T/G > A UV mutation signature in extracranial and brain metastatic tissues including unpaired samples. p-value indicates two-sided unpaired t-test. (B) Box and whiskers plot showing the frequency of C > T/G > A UV mutation signature across mutational subtypes of brain metastases. p-value from two-sided unpaired t-test.
Figure 4. (A) Box and whiskers plot showing the frequency of C > T/G > A UV mutation signature in extracranial and brain metastatic tissues including unpaired samples. p-value indicates two-sided unpaired t-test. (B) Box and whiskers plot showing the frequency of C > T/G > A UV mutation signature across mutational subtypes of brain metastases. p-value from two-sided unpaired t-test.
Cancers 13 00731 g004
Table
Table 1. Patient characteristics.
Table 1. Patient characteristics.
CharacteristicsNumber of CasesAge at Diagnosis [Years ± Range]Age at Surgery [Years ± Range]
Total patients47
Sex
Male2755 ± 1459 ± 14
Female2060 ± 1460 ± 11
Metastases72
Brain5456 ± 1457 ± 14
Skin1653 ± 2055 ± 1
Lymph node14445
Adrenal gland14747
Localization of BM
Cerebellar5--
Cerebral48--
Spinal1--
Therapy
Mono-CT2
Mono-RT4
Mono-IMT2
CT + RT1
CT + IMT3
RT + IMT5
Unknown30
Status --
Alive0--
Dead45--
Lost to follow up2--
BM: brain metastasis, CT: chemotherapy, RT: radiotherapy, IMT: immunotherapy.
Table
Table 2. Variants specific to brain-metastases of melanoma.
Table 2. Variants specific to brain-metastases of melanoma.
PairGeneChrRegionCovFreq (%)SubsProtein ChangeRelevance
9EZH271485235878214.6C > TC289YMissense-US
9ARID1A12708737481335.6C > TP650SMissense-US
9ARID1A1271014202665.6C > TP1568SMissense-US
9BAP13524398299032.2C > TV295MMissense-US
9BRAF714049424124465.2G > AP336LMissense-US
8SMARCA4191109760537839.2C > TP262LMissense-US
8MAP2K11566727441111848.3T > AF53IMissense-US
8ARID1A1271067451012357.1T > CL2119SMissense-US
8IDH12209106800..20910680122341.3GG > AAA212VMissense-US
8PIK3CA317895190033043.3C > AY1157*inactivating
8GNA11193115011..311501213040.8CC > TTR183Cactivating
2TERT512554741485.4G > AP1029SMissense-US
7IDH122091081642295.2C > TE229KMissense-US
7BRAF714048141179612.8C > TG73EMissense-US
7ARID1A1271014035136.0C > TP1562LMissense-US
7ARID1A1270895366866.1G > AG831AMissense-US
7SMARCA4191113253936310.2C > TP919SMissense-US
7PIK3CA31789274231116.3C > TL396FMissense-US
7NF117295600404856.4C > TL1173FMissense-US
7NF117296634034125.3C > TT1999IMissense-US
7SMARCA419111385819512.6G > AE1113KMissense-US
7ARID21246230381918.8G > AD239NMissense-US
7IDH122091132462165.1C > GK87NMissense-US
7PTEN10897176301156.1C > TQ219*inactivating
5BAP13524368291395.8A > TL650HMissense-US
5ARID2124624637237739.0C > TS1489LMissense-US
6ARID21246244970179525.5C > TP1022SMissense-US
10NRAS111525652937843.4T > CQ61Ractivating
10TP5317757838944038.2G > AR142CMissense-US
10TERT promoter5129525032433.0G > AC250Tactivating
11TP5317757841844493.5T > CE171GMissense-US
11BRAF714045313612454.0A > TV207Eactivating
11TP53177578418233228.1T > CE171GMissense-US
11GNAQ980537098..805370991716.4TG > GAP100LMissense-US
11PTEN108972083322429.0A > -A328fsinactivating
11BRAF714045313690623.2A > TV207Eactivating
1GNAQ9805371134209.3G > CD95EMissense-US
1GNAQ980537098..8053709921817.4TG > GAP100LMissense-US
3NRAS 111525653037041.9G > TQ61Kactivating
3SMARCA41911132434156328.9C > TH884YMissense-US
4TERT 5125459429434.4C > TA1062TMissense-US
4ARID1A 12705627275835.0G > AG423EMissense-US
4PIK3R1 5675226321817.2T > GS43RMissense-US
4PIK3R1 5675226281827.7T > GF41CMissense-US
4ARID2 124623064174924.0C > TS297FMissense-US
Brain specific variants discovered in 11 matched sample pairs. Chr: Chromosome, Cov: Coverage; Freq: Frequency of variant in sequencing reads, Subs: Substitution; AA: Amino acid, NA: Not available; US: unknown significance.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Váraljai, R.; Horn, S.; Sucker, A.; Piercianek, D.; Schmitt, V.; Carpinteiro, A.; Becker, K.A.; Reifenberger, J.; Roesch, A.; Felsberg, J.; et al. Integrative Genomic Analyses of Patient-Matched Intracranial and Extracranial Metastases Reveal a Novel Brain-Specific Landscape of Genetic Variants in Driver Genes of Malignant Melanoma. Cancers 2021, 13, 731. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13040731

AMA Style

Váraljai R, Horn S, Sucker A, Piercianek D, Schmitt V, Carpinteiro A, Becker KA, Reifenberger J, Roesch A, Felsberg J, et al. Integrative Genomic Analyses of Patient-Matched Intracranial and Extracranial Metastases Reveal a Novel Brain-Specific Landscape of Genetic Variants in Driver Genes of Malignant Melanoma. Cancers. 2021; 13(4):731. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13040731

Chicago/Turabian Style

Váraljai, Renáta, Susanne Horn, Antje Sucker, Daniela Piercianek, Verena Schmitt, Alexander Carpinteiro, Katrin Anne Becker, Julia Reifenberger, Alexander Roesch, Jörg Felsberg, and et al. 2021. "Integrative Genomic Analyses of Patient-Matched Intracranial and Extracranial Metastases Reveal a Novel Brain-Specific Landscape of Genetic Variants in Driver Genes of Malignant Melanoma" Cancers 13, no. 4: 731. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13040731

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