Longitudinal molecular trajectories of diffuse glioma in adults.

The evolutionary processes that drive universal therapeutic resistance in adult patients with diffuse glioma remain unclear1,2. Here we analysed temporally separated DNA-sequencing data and matched clinical annotation from 222 adult patients with glioma. By analysing mutations and copy numbers across the three major subtypes of diffuse glioma, we found that driver genes detected at the initial stage of disease were retained at recurrence, whereas there was little evidence of recurrence-specific gene alterations. Treatment with alkylating agents resulted in a hypermutator phenotype at different rates across the glioma subtypes, and hypermutation was not associated with differences in overall survival. Acquired aneuploidy was frequently detected in recurrent gliomas and was characterized by IDH mutation but without co-deletion of chromosome arms 1p/19q, and further converged with acquired alterations in the cell cycle and poor outcomes. The clonal architecture of each tumour remained similar over time, but the presence of subclonal selection was associated with decreased survival. Finally, there were no differences in the levels of immunoediting between initial and recurrent gliomas. Collectively, our results suggest that the strongest selective pressures occur during early glioma development and that current therapies shape this evolution in a largely stochastic manner.

Profiling of Copy Number Alterations Using Low-Coverage Whole-Genome Sequencing Informs Differential Diagnosis and Prognosis in Primary Cutaneous Follicle Center Lymphoma.

Primary cutaneous follicle center lymphoma (PCFCL) has an excellent prognosis using local treatment, whereas nodal follicular lymphoma (nFL), occasionally presenting with cutaneous spread, often requires systemic therapy. Distinction of the 2 diseases based on histopathology alone might be challenging. Copy number alterations (CNAs) have scarcely been explored on a genome-wide scale in PCFCL; however, they might serve as potential biomarkers during differential diagnosis and risk stratification. Low-coverage whole-genome sequencing is a robust, high-throughput method for genome-wide copy number profiling. In this study, we analyzed 28 PCFCL samples from 20 patients and compared the copy number profiles with a cohort of diagnostic samples of 64 nFL patients. Although the copy number profile of PCFCL was similar to that of nFL, PCFCL lacked amplifications of 18q, with the frequency peaking at 18q21.33 in nFL cases involving the BCL2 locus (PCFCL: 5.0% vs nFL: 31.3%, P = .018, Fisher exact test). Development of distant cutaneous spread was significantly associated with higher genomic instability including the proportion of genome altered (0.02 vs 0.13, P = .033) and number of CNAs (2 vs 9 P = .017), as well as the enrichment of 2p22.2-p15 amplification involving REL and XPO1 (6.3% vs 60.0%, P = .005), 3q23-q24 amplification (0.0% vs 50.0%, P = .004), 6q16.1-q23.3 deletion (6.3% vs 50.0%, P = .018), and 9p21.3 deletion covering CDKN2A and CDKN2B loci (0.0% vs 40.0%, P = .014, all Fisher exact test) in PCFCL. Analysis of sequential tumor samples in 2 cases harboring an unfavorable clinical course pointed to the acquisition of 2p amplification in the earliest common progenitor underlining its pivotal role in malignant transformation. By performing genome-wide copy number profiling on the largest patient cohort to date, we identified distinctive CNA alterations conceivably facilitating the differential diagnosis of PCFCL and secondary cutaneous involvement of nFL and potentially aiding the risk stratification of patients with PCFCL in the future.

PersonALL: a genetic scoring guide for personalized risk assessment in pediatric B-cell precursor acute lymphoblastic leukemia.

BACKGROUND: Recurrent genetic lesions provide basis for risk assessment in pediatric acute lymphoblastic leukemia (ALL). However, current prognostic classifiers rely on a limited number of predefined sets of alterations. METHODS: Disease-relevant copy number aberrations (CNAs) were screened genome-wide in 260 children with B-cell precursor ALL. Results were integrated with cytogenetic data to improve risk assessment. RESULTS: CNAs were detected in 93.8% (n = 244) of the patients. First, cytogenetic profiles were combined with IKZF1 status (IKZF1normal, IKZF1del and IKZF1plus) and three prognostic subgroups were distinguished with significantly different 5-year event-free survival (EFS) rates, IKAROS-low (n = 215): 86.3%, IKAROS-medium (n = 27): 57.4% and IKAROS-high (n = 18): 37.5%. Second, contribution of genetic aberrations to the clinical outcome was assessed and an aberration-specific score was assigned to each prognostically relevant alteration. By aggregating the scores of aberrations emerging in individual patients, personalized cumulative values were calculated and used for defining four prognostic subgroups with distinct clinical outcomes. Two favorable subgroups included 60% of patients (n = 157) with a 5-year EFS of 96.3% (excellent risk, n = 105) and 87.2% (good risk, n = 52), respectively; while 40% of patients (n = 103) showed high (n = 74) or ultra-poor (n = 29) risk profile (5-year EFS: 67.4% and 39.0%, respectively). CONCLUSIONS: PersonALL, our conceptually novel prognostic classifier considers all combinations of co-segregating genetic alterations, providing a highly personalized patient stratification.

Parallel testing of liquid biopsy (ctDNA) and tissue biopsy samples reveals a higher frequency of EZH2 mutations in follicular lymphoma.

BACKGROUND: Recent genomic studies revealed enhancer of zeste homolog 2 (EZH2) gain-of-function mutations, representing novel therapeutic targets in follicular lymphoma (FL) in around one quarter of patients. However, these analyses relied on single-site tissue biopsies and did not investigate the spatial heterogeneity and temporal dynamics of these alterations. OBJECTIVES: We aimed to perform a systematic analysis of EZH2 mutations using paired tissue (tumor biopsies [TB]) and liquid biopsies (LB) collected prior to treatment within the framework of a nationwide multicentric study. METHODS: Pretreatment LB and TB samples were collected from 123 patients. Among these, 114 had paired TB and LB, with 39 patients characterized with paired diagnostic and relapse samples available. The EZH2 mutation status and allele burden were assessed using an in-house-designed, highly sensitive multiplex droplet digital PCR assay. RESULTS: EZH2 mutation frequency was found to be 41.5% in the entire cohort. In patients with paired TB and LB samples, EZH2 mutations were identified in 37.8% of the patients with mutations exclusively found in 5.3% and 7.9% of TB and LB samples, respectively. EZH2 mutation status switch was documented in 35.9% of the patients with paired diagnostic and relapse samples. We also found that EZH2 wild-type clones may infiltrate the bone marrow more frequently compared to the EZH2 mutant ones. CONCLUSION: The in-depth spatio-temporal analysis identified EZH2 mutations in a considerably higher proportion of patients than previously reported. This expands the subset of FL patients who most likely would benefit from EZH2 inhibitor therapy.

Genomic Landscape of Normal and Breast Cancer Tissues in a Hungarian Pilot Cohort.

A limited number of studies have focused on the mutational landscape of breast cancer in different ethnic populations within Europe and compared the data with other ethnic groups and databases. We performed whole-genome sequencing of 63 samples from 29 Hungarian breast cancer patients. We validated a subset of the identified variants at the DNA level using the Illumina TruSight Oncology (TSO) 500 assay. Canonical breast-cancer-associated genes with pathogenic germline mutations were CHEK2 and ATM. Nearly all the observed germline mutations were as frequent in the Hungarian breast cancer cohort as in independent European populations. The majority of the detected somatic short variants were single-nucleotide polymorphisms (SNPs), and only 8% and 6% of them were deletions or insertions, respectively. The genes most frequently affected by somatic mutations were KMT2C (31%), MUC4 (34%), PIK3CA (18%), and TP53 (34%). Copy number alterations were most common in the NBN, RAD51C, BRIP1, and CDH1 genes. For many samples, the somatic mutational landscape was dominated by mutational processes associated with homologous recombination deficiency (HRD). Our study, as the first breast tumor/normal sequencing study in Hungary, revealed several aspects of the significantly mutated genes and mutational signatures, and some of the copy number variations and somatic fusion events. Multiple signs of HRD were detected, highlighting the value of the comprehensive genomic characterization of breast cancer patient populations.

Landscape of BCL2 Resistance Mutations in a Real-World Cohort of Patients with Relapsed/Refractory Chronic Lymphocytic Leukemia Treated with Venetoclax.

The oral, highly selective Bcl2 inhibitor venetoclax has substantially improved the therapeutic landscape of chronic lymphocytic leukemia (CLL). Despite the remarkable response rates in patients with relapsed/refractory (R/R) disease, acquired resistance is the leading cause of treatment failure, with somatic BCL2 mutations being the predominant genetic drivers underpinning venetoclax resistance. To assess the correlation between disease progression and the most common BCL2 mutations G101V and D103Y, sensitive (10-4) screening for the most common BCL2 mutations G101V and D103Y was performed in 67 R/R CLL patients during venetoclax single-agent or venetoclax-rituximab combination therapy. With a median follow-up time of 23 months, BCL2 G101V and D103Y were detected in 10.4% (7/67) and 11.9% (8/67) of the cases, respectively, with four patients harboring both resistance mutations. Ten out of eleven patients carrying BCL2 G101V and/or D103Y experienced relapse during the follow-up period, representing 43.5% of the cases (10/23) showing clinical signs of disease progression. All BCL2 G101V or D103Y variants were detected in patients receiving venetoclax as a continuous single-agent treatment while these mutations were not observed during or after fixed-duration venetoclax therapy. Targeted ultra-deep sequencing of BCL2 uncovered three additional variants in four patient samples obtained at relapse, suggesting convergent evolution and implying a cooperating role of BCL2 mutations in driving venetoclax resistance. This cohort is the largest R/R CLL patient population reported to date in which BCL2 resistance mutations were investigated. Our study demonstrates the feasibility and clinical value of sensitive screening for BCL2 resistance mutations in R/R CLL.

Screening and monitoring of the BTKC481S mutation in a real-world cohort of patients with relapsed/refractory chronic lymphocytic leukaemia during ibrutinib therapy.

The Bruton's tyrosine kinase (BTK) inhibitor ibrutinib has revolutionised the therapeutic landscape of chronic lymphocytic leukaemia (CLL). Acquired mutations emerging at position C481 in the BTK tyrosine kinase domain are the predominant genetic alterations associated with secondary ibrutinib resistance. To assess the correlation between disease progression, and the emergence and temporal dynamics of the most common resistance mutation BTKC481S , sensitive (10-4 ) time-resolved screening was performed in 83 relapsed/refractory CLL patients during single-agent ibrutinib treatment. With a median follow-up time of 40 months, BTKC481S was detected in 48·2% (40/83) of the patients, with 80·0% (32/40) of them showing disease progression during the examined period. In these 32 cases, representing 72·7% (32/44) of all patients experiencing relapse, emergence of the BTKC481S mutation preceded the symptoms of clinical relapse with a median of nine months. Subsequent Bcl-2 inhibition therapy applied in 28/32 patients harbouring BTKC481S and progressing on ibrutinib conferred clinical and molecular remission across the patients. Our study demonstrates the clinical value of sensitive BTKC481S monitoring with the largest longitudinally analysed real-world patient cohort reported to date and validates the feasibility of an early prediction of relapse in the majority of ibrutinib-treated relapsed/refractory CLL patients experiencing disease progression.

NGS-Based Application for Routine Non-Invasive Pre-Implantation Genetic Assessment in IVF.

Although non-invasive pre-implantation genetic testing for aneuploidy (NIPGT-A) is potentially appropriate to assess chromosomal ploidy of the embryo, practical application of it in a routine IVF centre have not been started in the absence of a recommendation. Our objective in this study was to provide a comprehensive workflow for a clinically applicable strategy for NIPGT-A based on next-generation sequencing (NGS) technology with the corresponding bioinformatic pipeline. In a retrospective study, we performed NGS on spent blastocyst culture media of Day 3 embryos fertilised with intracytoplasmic sperm injection (ICSI) with quality score on morphology assessment using the blank culture media as background control. Chromosomal abnormalities were identified by an optimised bioinformatics pipeline applying copy number variation (CNV) detecting algorithm. In this study, we demonstrate a comprehensive workflow covering both wet- and dry-lab procedures supporting a clinically applicable strategy for NIPGT-A that can be carried out within 48 h, which is critical for the same-cycle blastocyst transfer. The described integrated approach of non-invasive evaluation of embryonic DNA content of the culture media can potentially supplement existing pre-implantation genetic screening methods.

Dissection of subclonal evolution by temporal mutation profiling in chronic lymphocytic leukemia patients treated with ibrutinib.

The Bruton's tyrosine kinase (BTK) inhibitor ibrutinib is inducing durable responses in chronic lymphocytic leukemia (CLL) patients with refractory/relapsed disease or with TP53 defect, with BTK and phospholipase C gamma 2 (PLCG2) mutations representing the predominant mechanisms conferring secondary ibrutinib resistance. To understand the landscape of genomic changes and the dynamics of subclonal architecture associated with ibrutinib treatment, an ultra-deep next-generation sequencing analysis of 30 recurrently mutated genes was performed on sequential samples of 20 patients, collected before and during single-agent ibrutinib treatment. Mutations in the SF3B1, MGAand BIRC3 genes were enriched during ibrutinib treatment, while aberrations in the BTK, PLCG2, RIPK1, NFKBIE and XPO1 genes were exclusively detected in posttreatment samples. Besides the canonical mutations, four novel BTK mutations and three previously unreported PLCG2 variants were identified. BTK and PLCG2 mutations were backtracked in five patients using digital droplet PCR and were detectable on average 10.5 months before clinical relapse. With a median follow-up time of 36.5 months, 7/9 patients harboring BTK mutations showed disease progression based on clinical and/or laboratory features. In conclusion, subclonal heterogeneity, dynamic clonal selection and various patterns of clonal variegation were identified with novel resistance-associated BTK mutations in individual patients treated with ibrutinib.

Pooled CRISPR screening with single-cell transcriptome readout.

CRISPR-based genetic screens are accelerating biological discovery, but current methods have inherent limitations. Widely used pooled screens are restricted to simple readouts including cell proliferation and sortable marker proteins. Arrayed screens allow for comprehensive molecular readouts such as transcriptome profiling, but at much lower throughput. Here we combine pooled CRISPR screening with single-cell RNA sequencing into a broadly applicable workflow, directly linking guide RNA expression to transcriptome responses in thousands of individual cells. Our method for CRISPR droplet sequencing (CROP-seq) enables pooled CRISPR screens with single-cell transcriptome resolution, which will facilitate high-throughput functional dissection of complex regulatory mechanisms and heterogeneous cell populations.

Comprehensive profiling of disease-relevant copy number aberrations for advanced clinical diagnostics of pediatric acute lymphoblastic leukemia.

Acute lymphoblastic leukemia is the most common pediatric cancer characterized by a heterogeneous genomic landscape with copy number aberrations occurring at various stages of pathogenesis, disease progression, and treatment resistance. In this study, disease-relevant copy number aberrations were profiled in bone marrow samples of 91 children with B- or T-cell precursor acute lymphoblastic leukemia using digital multiplex ligation-dependent probe amplification (digitalMLPATM). Whole chromosome gains and losses, subchromosomal copy number aberrations, as well as unbalanced alterations conferring intrachromosomal gene fusions were simultaneously identified with results available within 36 hours. Aberrations were observed in 96% of diagnostic patient samples, and increased numbers of copy number aberrations were detected at the time of relapse as compared with diagnosis. Comparative scrutiny of 24 matching diagnostic and relapse samples from 11 patients revealed three different patterns of clonal relationships with (i) one patient displaying identical copy number aberration profiles at diagnosis and relapse, (ii) six patients showing clonal evolution with all lesions detected at diagnosis being present at relapse, and (iii) four patients displaying conserved as well as lost or gained copy number aberrations at the time of relapse, suggestive of the presence of a common ancestral cell compartment giving rise to clinically manifest leukemia at different time points during the disease course. A newly introduced risk classifier combining cytogenetic data with digitalMLPATM-based copy number aberration profiles allowed for the determination of four genetic subgroups of B-cell precursor acute lymphoblastic leukemia with distinct event-free survival rates. DigitalMLPATM provides fast, robust, and highly optimized copy number aberration profiling for the genomic characterization of acute lymphoblastic leukemia samples, facilitates the decipherment of the clonal origin of relapse and provides highly relevant information for clinical prognosis assessment.

[Multiplex ligation-dependent probe amplification in oncohematological diagnostics and research]. / Multiplex ligatiofüggo szondaamplifikáció az onkohematológiai kutatásban és diagnosztikában.

Genetic abnormalities associated with the development, progression and treatment resistance of hematological malignancies are extensively characterized. Rapid, reliable and cost-efficient techniques are needed to screen the clinically relevant aberrations in routine diagnostics. Multiplex ligation-dependent probe amplification is an efficient tool to analyze genomic copy number aberrations at 55-60 different genomic loci. The method allows the profiling of prognostic and predictive markers; thus, it can efficiently be combined with karyotyping and fluorescence in situ hybridization, the most commonly used diagnostic techniques to detect cytogenetic lesions. Furthermore, the method can interrogate methylation status and unravel point mutations at specific sites, providing results in 24 hours. Here, we describe the technical background of multiplex ligation-dependent probe amplification, summarize its advantages and limitations as well as discuss its role in oncohematological diagnostics and research. Finally, future outlook is provided, with emphasis on recent technological advances related to next-generation sequencing. Orv Hetil. 2018; 159(15): 583-592.

[Introduction to the molecular diagnostic methods of oncohematology]. / Az onkohematológia molekuláris diagnosztikai vizsgálómódszereinek alapjai.

Owing to our rapidly expanding knowledge on the genetic background of various oncohematologic diseases and the introduction of novel targeted therapies, molecular genetic techniques have been playing an increasingly important role in the diagnostics and follow-up of hematological malignancies. The various DNA- and RNA-based in situ hybridization, polymerase chain reaction and sequencing technologies are of key significance in diagnostics, classification and prognostic assessment of these diseases, as well as in the monitoring of minimal residual disease and selection of the most appropriate targeted therapy. This review provides an overview on the background and applications of the molecular methods most commonly used in oncohematological diagnostics.

[The genetic landscape of chronic lymphocytic leukemia]. / A krónikus limfocitás leukémia genetikai háttere az újgenerációs szekvenálás korszakában.

Chronic lymphocytic leukemia (CLL) is the most frequent mature B-cell non-Hodgkin's lymphoma in the Western countries. The recent next-generation sequencing (NGS) studies lead to an exponential increase in our knowledge of the pathogenesis and progression of CLL. Whole genome and exome sequencing studies revealed a remarkable inter- and intra-patient genetic heterogeneity with a significant therapy-induced clonal evolution in the majority of the patients. Driver mutations were identified in components of various signalling pathways and cellular processes with notable prognostic and therapeutic relevance. Interestingly, these studies revealed only a few genes mutated in at least 15-20% of the patients with a larger number of genes mutated in a smaller proportion of patients. This improved understanding of the genomic landscape of CLL has opened new avenues for a more precise patient stratification and rational application of novel, more effective targeted therapies.

Sequential and hierarchical chromosomal changes and chromosome instability are distinct features of high hyperdiploid pediatric acute lymphoblastic leukemia.

BACKGROUND: Pathogenesis of the non-random accumulation of extra chromosomes in the low and high hyperdiploid (HeL, HeH) pre-B pediatric acute lymphoblastic leukemia (B-pALL) is largely unknown, and has been clarified with respect only to tetrasomic chromosomes. We analyzed the hierarchy of changes in chromosome number and chromosomal instability, as well as clonal heterogeneity and evolution, in the untreated bone marrow cell samples from 214 B-pALL patients. PROCEDURE: Applying relocation, 2 × 4 color interphase fluorescence in situ hybridization was used to detect copy number alterations (CNAs) of the most commonly involved chromosomes, 4, 6, 10, 14, 17, 18, 21, and X. This approach allowed us to acquire a dataset correlated for all eight parameters. RESULTS: Based on chromosome number, an average of 6.9 and 10.2, whereas according to unique constellation 15.3 and 26.7 subclones could be identified in the HeL and HeH subgroups, respectively. Cluster analysis revealed the order of CNAs to chromosomes was highly conserved, and network analysis indicated changes in chromosome number were sequential for 80-90% of all numerical aberrations. Significant chromosome instability was revealed in both subgroups of leukemia. CONCLUSIONS: Data generated using this new approach indicate that chromosomal instability, which causes heterogeneity in the subclonal landscape, and the sequential changes to chromosome numbers, are both determining factors in the pathomechanism of the hyperdiploid B-pALL. These new observations could prompt research into the mitotic machinery of leukemic cells to identify new therapeutic targets for treating this disease.

Multiplex ligation-dependent probe amplification and fluorescence in situ hybridization are complementary techniques to detect cytogenetic abnormalities in multiple myeloma.

Multiple myeloma (MM) is a genetically heterogeneous disease with diverse clinical outcomes. Interphase fluorescence in situ hybridization (i-FISH) is the most commonly used approach to detect recurrent cytogenetic abnormalities in this malignancy. We aimed to assess the performance of multiplex ligation-dependent probe amplification (MLPA) to reveal copy number abnormalities (CNAs) in MM. Diagnostic bone marrow samples from 81 patients were analyzed using 42 MLPA probes for the following regions: 1p32-31, 1p21, 1q21.3, 1q23.3, 5q31.3, 12p13.31, 13q14, 16q12, 16q23, and 17p13. All samples were also screened by i-FISH for the presence of hyperdiploidy, deletion/monosomy of chromosome 13, deletion of TP53, disruption of the immunoglobulin heavy-chain gene, t(4;14), t(11;14), t(14;16), t(8;14), gain of 5q and abnormalities of chromosome 1. A total of 245 alterations were detected in 79 cases (98%). Investigating the same aberrations, the two methods showed a congruency of higher than 90%. A low proportion of cells with the relevant abnormality, focal CNAs and unmatched probes were responsible for the discrepancies. MLPA revealed 95 CNAs not detected by i-FISH providing additional information in 53 cases (65%). Scrutiny of CNAs on chromosome 1, using more than 20 probes, revealed significant heterogeneity in size and location, and variable intra-chromosomal and intra-clonal rates of loss or gain. Our results suggest that MLPA is a reliable high-throughput technique to detect CNAs in MM. Since balanced aberrations are key to prognostic classification of this disease, MLPA and i-FISH should be applied as complementary techniques in diagnostic pathology.

Co-targeting SOS1 enhances the antitumor effects of KRASG12C inhibitors by addressing intrinsic and acquired resistance.

Combination approaches are needed to strengthen and extend the clinical response to KRASG12C inhibitors (KRASG12Ci). Here, we assessed the antitumor responses of KRASG12C mutant lung and colorectal cancer models to combination treatment with a SOS1 inhibitor (SOS1i), BI-3406, plus the KRASG12C inhibitor, adagrasib. We found that responses to BI-3406 plus adagrasib were stronger than to adagrasib alone, comparable to adagrasib with SHP2 (SHP2i) or EGFR inhibitors and correlated with stronger suppression of RAS-MAPK signaling. BI-3406 plus adagrasib treatment also delayed the emergence of acquired resistance and elicited antitumor responses from adagrasib-resistant models. Resistance to KRASG12Ci seemed to be driven by upregulation of MRAS activity, which both SOS1i and SHP2i were found to potently inhibit. Knockdown of SHOC2, a MRAS complex partner, partially restored response to KRASG12Ci treatment. These results suggest KRASG12C plus SOS1i to be a promising strategy for treating both KRASG12Ci naive and relapsed KRASG12C-mutant tumors.

Low-burden TP53 mutations represent frequent genetic events in CLL with an increased risk for treatment initiation.

TP53 aberrations predict chemoresistance and represent a contraindication for the use of standard chemoimmunotherapy in chronic lymphocytic leukaemia (CLL). Recent next-generation sequencing (NGS)-based studies have identified frequent low-burden TP53 mutations with variant allele frequencies below 10%, but the clinical impact of these low-burden TP53 mutations is still a matter of debate. In this study, we aimed to scrutinise the subclonal architecture and clinical impact of TP53 mutations using a sensitive, NGS-based mutation analysis in a 'real-world' cohort of 901 patients with CLL. In total, 225 TP53 mutations were identified in 17.5% (158/901) of the patients; 48% of these alterations represented high-burden mutations, while 52% were low-burden TP53 mutations. Low-burden mutations as sole alterations were identified in 39% (62/158) of all mutated cases with 82% (51/62) of these being represented by a single low-burden TP53 mutation. Patients harbouring low-burden TP53 mutations had significantly lower time to first treatment compared to patients with wild-type TP53. Our study has expanded the knowledge on the frequency, clonal architecture, and clinical impact of low-burden TP53 mutations. By demonstrating that patients with sole low-burden TP53 variants represent more than one-third of patients with TP53 mutations and have an increased risk for treatment initiation, our findings strengthen the need to redefine the threshold of TP53 variant reporting to below 10% in the routine diagnostic setting.

The Driverless Triple-Wild-Type (BRAF, RAS, KIT) Cutaneous Melanoma: Whole Genome Sequencing Discoveries.

The genetic makeup of the triple-wild-type melanoma (BRAF, NRAS and NF1) has been known for some time, but those studies grouped together rare histopathological versions with common ones, as well as mucosal and even uveal ones. Here we used whole genome sequencing to genetically characterize the triple-wild-type melanoma (TWM), termed here as BRAF, RAS and KIT wild type (the most frequent oncogenic drivers of skin melanoma), using the most common histological forms and excluding rare ones. All these tumors except one were clearly induced by UV based on the mutational signature. The tumor mutational burden was low in TWM, except in the NF1 mutant forms, and a relatively high frequency of elevated LOH scores suggested frequent homologue recombination deficiency, but this was only confirmed by the mutation signature in one case. Furthermore, all these TWMs were microsatellite-stabile. In this driverless setting, we revealed rare oncogenic drivers known from melanoma or other cancer types and identified rare actionable tyrosine kinase mutations in NTRK1, RET and VEGFR1. Mutations of TWM identified genes involved in antitumor immunity (negative and positive predictors of immunotherapy), Ca++ and BMP signaling. The two regressed melanomas of this cohort shared a 17-gene mutation signature, containing genes involved in antitumor immunity and several cell surface receptors. Even with this comprehensive genomic approach, a few cases remained driverless, suggesting that unrecognized drivers are hiding among passenger mutations.