Comprehensive genomic profiling on metastatic Melanoma: results from a network screening from 7 Italian Cancer Centres.

BACKGROUND: The current therapeutic algorithm for Advanced Stage Melanoma comprises of alternating lines of Targeted and Immuno-therapy, mostly via Immune-Checkpoint blockade. While Comprehensive Genomic Profiling of solid tumours has been approved as a companion diagnostic, still no approved predictive biomarkers are available for Melanoma aside from BRAF mutations and the controversial Tumor Mutational Burden. This study presents the results of a Multi-Centre Observational Clinical Trial of Comprehensive Genomic Profiling on Target and Immuno-therapy treated advanced Melanoma. METHODS: 82 samples, collected from 7 Italian Cancer Centres of FFPE-archived Metastatic Melanoma and matched blood were sequenced via a custom-made 184-gene amplicon-based NGS panel. Sequencing and bioinformatics analysis was performed at a central hub. Primary analysis was carried out via the Ion Reporter framework. Secondary analysis and Machine Learning modelling comprising of uni and multivariate, COX/Lasso combination, and Random Forest, was implemented via custom R/Python scripting. RESULTS: The genomics landscape of the ACC-mela cohort is comparable at the somatic level for Single Nucleotide Variants and INDELs aside a few gene targets. All the clinically relevant targets such as BRAF and NRAS have a comparable distribution thus suggesting the value of larger scale sequencing in melanoma. No comparability is reached at the CNV level due to biotechnological biases and cohort numerosity. Tumour Mutational Burden is slightly higher in median for Complete Responders but fails to achieve statistical significance in Kaplan-Meier survival analysis via several thresholding strategies. Mutations on PDGFRB, NOTCH3 and RET were shown to have a positive effect on Immune-checkpoint treatment Overall and Disease-Free Survival, while variants in NOTCH4 were found to be detrimental for both endpoints. CONCLUSIONS: The results presented in this study show the value and the challenge of a genomics-driven network trial. The data can be also a valuable resource as a validation cohort for Immunotherapy and Target therapy genomic biomarker research.

Molecular Profile of Subungual Melanoma: A MelaNostrum Consortium Study of 68 Cases Reporting BRAF, NRAS, KIT, and TERT Promoter Status.

BACKGROUND: Subungual melanoma (SM) is an unusual type of melanocytic tumor affecting the nail apparatus. The mutational prevalence of the most prominently mutated genes in melanoma has been reported in small cohorts of SM, with unclear conclusions on whether SM is different from the rest of melanomas arising in acral locations or not. Hence, the molecular profile of a large series of SM is yet to be described. OBJECTIVES: The aim of this study was to describe the molecular characteristics of a large series of SM and their association with demographic and histopathological features. METHODS: Patients diagnosed with SM between 2001 and 2021 were identified from six Spanish and Italian healthcare centers. The mutational status for BRAF, NRAS, KIT, and the promoter region of TERT (TERTp) were determined either by Sanger sequencing or next-generation sequencing. Clinical data were retrieved from the hospital databases to elucidate potential associations. RESULTS: A total of 68 SM cases were included. Mutations were most common in BRAF (10.3%) and KIT (10%), followed by NRAS (7.6%), and TERTp (3.8%). Their prevalence was similar to that of non-subungual acral melanoma but higher in SM located on the hand than on the foot. CONCLUSIONS: To date, this study represents the largest cohort of SM patients with data on the known driver gene mutations. The low mutation rate supports a different etiopathogenic mechanism for SM in comparison of non-acral cutaneous melanoma, particularly for SM of the foot.

Molecular Modeling Unveils the Effective Interaction of B-RAF Inhibitors with Rare B-RAF Insertion Variants.

The Food and Drug Administration (FDA) has approved MAPK inhibitors as a treatment for melanoma patients carrying a mutation in codon V600 of the BRAF gene exclusively. However, BRAF mutations outside the V600 codon may occur in a small percentage of melanomas. Although these rare variants may cause B-RAF activation, their predictive response to B-RAF inhibitor treatments is still poorly understood. We exploited an integrated approach for mutation detection, tumor evolution tracking, and assessment of response to treatment in a metastatic melanoma patient carrying the rare p.T599dup B-RAF mutation. He was addressed to Dabrafenib/Trametinib targeted therapy, showing an initial dramatic response. In parallel, in-silico ligand-based homology modeling was set up and performed on this and an additional B-RAF rare variant (p.A598_T599insV) to unveil and justify the success of the B-RAF inhibitory activity of Dabrafenib, showing that it could adeptly bind both these variants in a similar manner to how it binds and inhibits the V600E mutant. These findings open up the possibility of broadening the spectrum of BRAF inhibitor-sensitive mutations beyond mutations at codon V600, suggesting that B-RAF V600 WT melanomas should undergo more specific investigations before ruling out the possibility of targeted therapy.

Per- and polyfluoroalkyl substances (PFAS) exposure in melanoma patients: a retrospective study on prognosis and histological features.

Per- and polyfluoroalkyl substances (PFAS) are endocrine disrupting chemicals which could be associated with cancer development, such as kidney and testicular cancers, pancreatic and hepatocellular carcinoma and thyroid tumor. Available scientific literature offers no information on the role of PFAS in melanoma development/progression. Since 1965, a massive environmental contamination by PFAS has occurred in northeastern Italy. This study compared histopathology and prognosis between melanoma patients exposed (n = 194) and unexposed (n = 488) to PFAS. All patients were diagnosed and/or treated for melanoma at the Veneto Oncological Institute and the University Hospital of Padua (Italy) in 1998-2014. Patients were categorized in exposed or unexposed groups according to their home address and the geographical classification of municipalities affected by PFAS contamination as provided by Veneto Government in 2018. Presence of mitoses was found in 70.5% of exposed patients and 58.7% of unexposed patients (p = 0.005). Median follow-up was 90 months (IQR 59-136). 5-year overall survival was 83.7% in exposed patients and 88.0% in unexposed patients (p = 0.20); 5-year disease-specific survival was 88.0% in exposed patients and 90.9% in unexposed patients (p = 0.50); 5-year disease-free survival was 83.8% in exposed patients and 87.3% in unexposed patients (p = 0.20). Adjusting for imbalanced characteristics at baseline (presence of mitoses), survival was not statistically different between exposed and unexposed patients (overall survival: HR 1.10, 95% CI 0.77 to 1.58, p = 0.57; disease-specific survival: HR 0.99, 95% CI 0.62 to 1.59, p = 0.99; disease-free survival: HR 1.10, 95% CI 0.74 to 1.64, p = 0.62). Although the magnitude of PFAS exposure was not quantifiable, our findings suggested that exposure to PFAS was associated with higher level of mitosis in melanoma patients, but this did not translate into a survival difference. Further studies are required to investigate this relationship and all effects of PFAS on prognosis.

Contribution of Common Genetic Variants to Familial Aggregation of Disease and Implications for Sequencing Studies.

Despite genetics being accepted as the primary cause of familial aggregation for most diseases, it is still unclear whether afflicted families are likely to share a single highly penetrant rare variant, many minimally penetrant common variants, or a combination of the two types of variants. We therefore use recent estimates of SNP heritability and the liability threshold model to estimate the proportion of afflicted families likely to carry a rare, causal variant. We then show that Polygenic Risk Scores (PRS) may be useful for identifying families likely to carry such a rare variant and therefore for prioritizing families to include in sequencing studies with that aim. Specifically, we introduce a new statistic that estimates the proportion of individuals carrying causal rare variants based on the family structure, disease pattern, and PRS of genotyped individuals. Finally, we consider data from the MelaNostrum consortium and show that, despite an estimated PRS heritability of only 0.05 for melanoma, families carrying putative causal variants had a statistically significantly lower PRS, supporting the idea that PRS prioritization may be a useful future tool. However, it will be important to evaluate whether the presence of rare mendelian variants are generally associated with the proposed test statistic or lower PRS in future and larger studies.

Ataxia-Telangiectasia Mutated Loss of Heterozygosity in Melanoma.

ATM germline pathogenic variants were recently found enriched in high-risk melanoma patients. However, ATM loss of heterozygosity (LOH) has never been investigated in melanoma and, therefore, a causal association with melanoma development has not been established yet. The purpose of this study was to functionally characterize 13 germline ATM variants found in high-risk melanoma patients-and classified by in silico tools as pathogenic, uncertain significance, or benign-using multiple assays evaluating ATM/pATM expression and/or LOH in melanoma tissues and cell lines. We assessed ATM status by Immunohistochemistry (IHC), Western Blot, Whole-Exome Sequencing/Copy Number Variation analysis, and RNA sequencing, supported by Sanger sequencing and microsatellite analyses. For most variants, IHC results matched those obtained with in silico classification and LOH analysis. Two pathogenic variants (p.Ser1135_Lys1192del and p.Ser1993ArgfsTer23) showed LOH and complete loss of ATM activation in melanoma. Two variants of unknown significance (p.Asn358Ile and p.Asn796His) showed reduced expression and LOH, suggestive of a deleterious effect. This study, showing a classic two-hit scenario in a well-known tumor suppressor gene, supports the inclusion of melanoma in the ATM-related cancer spectrum.

Gender-specific associations between polymorphisms of the circadian gene RORA and cutaneous melanoma susceptibility.

BACKGROUND: Melanoma is the deadliest of skin cancers and has an increasing annual incidence worldwide. It is a multi-factorial disease most likely arising from both genetic predisposition and environmental exposure to ultraviolet light. Genetic variability of the components of the biological circadian clock is recognized to be a risk factor for different type of cancers. Moreover, two variants of a clock gene, RORA, have been associated with melanoma patient's prognosis. Our aim is to test the hypothesis that specific single nucleotide polymorphisms (SNPs) of the circadian clock genes may significantly influence the predisposition to develop cutaneous melanoma or the outcome of melanoma patients. METHODS: We genotyped 1239 subjects, 629 cases of melanoma and 610 healthy controls in 14 known SNPs of seven selected clock genes: AANAT, CLOCK, NPAS2, PER1, PER2, RORA, and TIMELESS. Genotyping was conducted by q-PCR. Multivariate logistic regression was employed for susceptibility of melanoma assessment, modeled additively. Subgroup analysis was performed by gender. For the female subgroup, a further discrimination was performed by age. For prognosis of melanoma assessment, multivariate Cox proportional hazard regression was employed. The Benjamini-Hochberg method was utilized as adjustment for multiple comparisons. RESULTS: We identified two RORA SNPs statistically significant with respect to the association with melanoma susceptibility. Considering the putative role of RORA as a nuclear steroid hormone receptor, we conducted a subgroup analysis by gender. Interestingly, the RORA rs339972 C allele was associated with a decreased predisposition to develop melanoma only in the female subgroup (OR 0.67; 95% CI 0.51-0.88; P = 0.003) while RORA rs10519097 T allele was associated with a decreased predisposition to develop melanoma only in the male subgroup (OR 0.62; 95% CI 0.44-0.87; P = 0.005). Moreover, the RORA rs339972 C allele had a decreased susceptibility to develop melanoma only in females aged over 50 years old (OR 0.67; 95% CI 0.54-0.83; P = 0.0002). None of the studied SNPs were significantly associated with the prognosis. CONCLUSIONS: Overall, we cannot ascertain that circadian pathway genetic variation is involved in melanoma susceptibility or prognosis. Nevertheless, we identified an interesting relationship between melanoma susceptibility and RORA polymorphisms acting in sex-specific manner and which is worth further future investigation.

The Molecular Landscape of Primary Acral Melanoma: A Multicenter Study of the Italian Melanoma Intergroup (IMI).

Acral melanoma (AM) is a rare and aggressive subtype of melanoma affecting the palms, soles, and nail apparatus with similar incidence among different ethnicities. AM is unrelated to ultraviolet radiation and has a low mutation burden but frequent chromosomal rearrangements and gene amplifications. Next generation sequencing of 33 genes and somatic copy number variation (CNV) analysis with genome-wide single nucleotide polymorphism arrays were performed in order to molecularly characterize 48 primary AMs of Italian patients in association with clinicopathological and prognostic features. BRAF was the most commonly mutated gene, followed by NRAS and TP53, whereas TERT promoter, KIT, and ARID1A were less frequently mutated. Gains and losses were recurrently found in the 1q, 6p, 7, 8q, 20 and 22 chromosomes involving PREX2, RAC1, KMT2C, BRAF, CCND1, TERT, and AKT3 genes, and in the 6q, 9, 10, 11q and 16q chromosomes including CDKN2A, PTEN, and ADAMTS18 genes, respectively. This study confirmed the variety of gene mutations and the high load of CNV in primary AM. Some genomic alterations were associated with histologic prognostic features. BRAF mutations, found with a higher rate than previously reported, correlated with a low Breslow thickness, low mitotic count, low CNV of the AMs, and with early-stage of disease.

Combining common genetic variants and non-genetic risk factors to predict risk of cutaneous melanoma.

Melanoma heritability is among the highest for cancer and single nucleotide polymorphisms (SNPs) contribute to it. To date, only SNPs that reached statistical significance in genome-wide association studies or few candidate SNPs have been included in melanoma risk prediction models. We compared four approaches for building polygenic risk scores (PRS) using 12 874 melanoma cases and 23 203 controls from Melanoma Meta-Analysis Consortium as a training set, and newly genotyped 3102 cases and 2301 controls from the MelaNostrum consortium for validation. We estimated adjusted odds ratios (ORs) for melanoma risk using traditional melanoma risk factors and the PRS with the largest area under the receiver operator characteristics curve (AUC). We estimated absolute risks combining the PRS and other risk factors, with age- and sex-specific melanoma incidence and competing mortality rates from Italy as an example. The best PRS, including 204 SNPs (AUC = 64.4%; 95% confidence interval (CI) = 63-65.8%), developed using winner's curse estimate corrections, had a per-quintile OR = 1.35 (95% CI = 1.30-1.41), corresponding to a 3.33-fold increase comparing the 5th to the 1st PRS quintile. The AUC improvement by adding the PRS was up to 7%, depending on adjusted factors and country. The 20-year absolute risk estimates based on the PRS, nevus count and pigmentation characteristics for a 60-year-old Italian man ranged from 0.5 to 11.8% (relative risk  = 26.34), indicating good separation.

MelaNostrum: a consensus questionnaire of standardized epidemiologic and clinical variables for melanoma risk assessment by the melanostrum consortium.

BACKGROUND: Many melanoma observational studies have been carried out across different countries and geographic areas using heterogeneous assessments of epidemiologic risk factors and clinical variables. AIM: To develop a consensus questionnaire to standardize epidemiologic and clinical data collection for melanoma risk assessment. METHODS: We used a stepwise strategy that included: compilation of variables from case-control datasets collected at various centres of the MelaNostrum Consortium; integration of variables from published case-control studies; consensus discussion of the collected items by MelaNostrum members; revision by independent experts; addition of online tools and image-based charts; questionnaire testing across centres and generation of a final draft. RESULTS: We developed a core consensus questionnaire (MelanoQ) that includes four separate sections: A. general and demographic data; B. phenotypic and ultraviolet radiation exposure risk factors and lifestyle habits; C. clinical examination, medical and family history; and D. diagnostic data on melanoma (cases only). Accompanying online tools, informative tables, and image-based charts aid standardization. Different subsections of the questionnaire are designed for self-administration, patient interviews performed by a physician or study nurse, and data collection from medical records. CONCLUSIONS: The MelanoQ questionnaire is a useful tool for the collection and standardization of epidemiologic and clinical data across different studies, centres, cultures and languages. This will expedite ongoing efforts to compile high-quality data for pooled analyses or meta-analyses and offer a solid base for the design of clinical, epidemiologic and translational studies on melanoma.

Performance of in silico tools for the evaluation of p16INK4a (CDKN2A) variants in CAGI.

Correct phenotypic interpretation of variants of unknown significance for cancer-associated genes is a diagnostic challenge as genetic screenings gain in popularity in the next-generation sequencing era. The Critical Assessment of Genome Interpretation (CAGI) experiment aims to test and define the state of the art of genotype-phenotype interpretation. Here, we present the assessment of the CAGI p16INK4a challenge. Participants were asked to predict the effect on cellular proliferation of 10 variants for the p16INK4a tumor suppressor, a cyclin-dependent kinase inhibitor encoded by the CDKN2A gene. Twenty-two pathogenicity predictors were assessed with a variety of accuracy measures for reliability in a medical context. Different assessment measures were combined in an overall ranking to provide more robust results. The R scripts used for assessment are publicly available from a GitHub repository for future use in similar assessment exercises. Despite a limited test-set size, our findings show a variety of results, with some methods performing significantly better. Methods combining different strategies frequently outperform simpler approaches. The best predictor, Yang&Zhou lab, uses a machine learning method combining an empirical energy function measuring protein stability with an evolutionary conservation term. The p16INK4a challenge highlights how subtle structural effects can neutralize otherwise deleterious variants.

Multiple primary melanomas (MPMs) and criteria for genetic assessment: MultiMEL, a multicenter study of the Italian Melanoma Intergroup.

BACKGROUND: Multiple primary melanoma (MPM), in concert with a positive family history, is a predictor of cyclin-dependent kinase (CDK) inhibitor 2A (CDKN2A) germline mutations. A rule regarding the presence of either 2 or 3 or more cancer events (melanoma and pancreatic cancer) in low or high melanoma incidence populations, respectively, has been established to select patients for genetic referral. OBJECTIVE: We sought to determine the CDKN2A/CDK4/microphthalmia-associated transcription factor mutation rate among Italian patients with MPM to appropriately direct genetic counseling regardless of family history. METHODS: In all, 587 patients with MPM and an equal number with single primary melanomas and control subjects were consecutively enrolled at the participating centers and tested for CDKN2A, CDK4, and microphthalmia-associated transcription factor. RESULTS: CDKN2A germline mutations were found in 19% of patients with MPM versus 4.4% of patients with single primary melanoma. In familial MPM cases the mutation rate varied from 36.6% to 58.8%, whereas in sporadic MPM cases it varied from 8.2% to 17.6% in patients with 2 and 3 or more melanomas, respectively. The microphthalmia-associated transcription factor E318K mutation accounted for 3% of MPM cases altogether. LIMITATIONS: The study was hospital based, not population based. Rare novel susceptibility genes were not tested. CONCLUSION: Italian patients who developed 2 melanomas, even in situ, should be referred for genetic counseling even in the absence of family history.

CDKN2A unclassified variants in familial malignant melanoma: combining functional and computational approaches for their assessment.

CDKN2A codes for two oncosuppressors by alternative splicing of two first exons: p16INK4a and p14ARF. Germline mutations are found in about 40% of melanoma-prone families, and most of them are missense mutations mainly affecting p16INK4a. A growing number of p16INK4a variants of uncertain significance (VUS) are being identified but, unless their pathogenic role can be demonstrated, they cannot be used for identification of carriers at risk. Predicting the effect of these VUS by either a "standard" in silico approach, or functional tests alone, is rather difficult. Here, we report a protocol for the assessment of any p16INK4a VUS, which combines experimental and computational tools in an integrated approach. We analyzed p16INK4a VUS from melanoma patients as well as variants derived through permutation of conserved p16INK4a amino acids. Variants were expressed in a p16INK4a-null cell line (U2-OS) and tested for their ability to block proliferation. In parallel, these VUS underwent in silico prediction analysis and molecular dynamics simulations. Evaluation of in silico and functional data disclosed a high agreement for 15/16 missense mutations, suggesting that this approach could represent a pilot study for the definition of a protocol applicable to VUS in general, involved in other diseases, as well.

A multiparameter liquid biopsy approach allows to track melanoma dynamics and identify early treatment resistance.

Melanoma heterogeneity is a hurdle in metastatic disease management. Although the advent of targeted therapy has significantly improved patient outcomes, the occurrence of resistance makes monitoring of the tumor genetic landscape mandatory. Liquid biopsy could represent an important biomarker for the real-time tracing of disease evolution. Thus, we aimed to correlate liquid biopsy dynamics with treatment response and progression by devising a multiplatform approach applied to longitudinal melanoma patient monitoring. We conceived an approach that exploits Next Generation Sequencing (NGS) and droplet digital PCR, as well as the FDA-cleared platform CellSearch, to analyze circulating tumor DNA (ctDNA) trend and circulating melanoma cell (CMC) count, together with their customized genetic and copy number variation analysis. The approach was applied to 17 stage IV melanoma patients treated with BRAF/MEK inhibitors, followed for up to 28 months. BRAF mutations were detected in the plasma of 82% of patients. Single nucleotide variants known or suspected to confer resistance were identified in 70% of patients. Moreover, the amount of ctDNA, both at baseline and during response, correlated with the type and duration of the response itself, and the CMC count was confirmed to be a prognostic biomarker. This work provides proof of principle of the power of this approach and paves the way for a validation study aimed at evaluating early ctDNA-guided treatment decisions in stage IV melanoma. The NGS-based molecular profile complemented the analysis of ctDNA trend and, together with CMC analysis, revealed to be useful in capturing tumor evolution.

Altered tumor formation and evolutionary selection of genetic variants in the human MDM4 oncogene.

A large body of evidence strongly suggests that the p53 tumor suppressor pathway is central in reducing cancer frequency in vertebrates. The protein product of the haploinsufficient mouse double minute 2 (MDM2) oncogene binds to and inhibits the p53 protein. Recent studies of human genetic variants in p53 and MDM2 have shown that single nucleotide polymorphisms (SNPs) can affect p53 signaling, confer cancer risk, and suggest that the pathway is under evolutionary selective pressure (1-4). In this report, we analyze the haplotype structure of MDM4, a structural homolog of MDM2, in several different human populations. Unusual patterns of linkage disequilibrium (LD) in the haplotype distribution of MDM4 indicate the presence of candidate SNPs that may also modify the efficacy of the p53 pathway. Association studies in 5 different patient populations reveal that these SNPs in MDM4 confer an increased risk for, or early onset of, human breast and ovarian cancers in Ashkenazi Jewish and European cohorts, respectively. This report not only implicates MDM4 as a key regulator of tumorigenesis in the human breast and ovary, but also exploits for the first time evolutionary driven linkage disequilibrium as a means to select SNPs of p53 pathway genes that might be clinically relevant.

Altitude Effect on Cutaneous Melanoma Epidemiology in the Veneto Region (Northern Italy): A Pilot Study.

The incidence of cutaneous melanoma has been increasing in the last decades among the fair-skinned population. Despite its complex and multifactorial etiology, the exposure to ultraviolet radiation (UVR) is the most consistent modifiable risk factor for melanoma. Several factors influence the amount of UVR reaching the Earth's surface. Our study aimed to explore the relationship between melanoma and altitude in an area with mixed geographic morphology, such as the Veneto region (Italy). We included 2752 melanoma patients who were referred to our centers between 1998 and 2014. Demographics, histological and clinical data, and survival information were extracted from a prospectively maintained local database. Head/neck and acral melanoma were more common in patients from the hills and the mountains, while limb and trunk melanoma were more common in patients living in plain and coastal areas. Breslow thickness, ulceration and mitotic rate impaired with increased altitude. However, the geographical area of origin was not associated with overall or disease-free survival. The geographical area of origin of melanoma patients and the "coast-plain-hill gradient" could help to estimate the influence of different sun exposure and to explain the importance of vitamin D levels in skin-cancer control.

Diagnostics of BAP1-Tumor Predisposition Syndrome by a Multitesting Approach: A Ten-Year-Long Experience.

Germline mutations in the tumor suppressor gene BRCA1-associated protein-1 (BAP1) lead to BAP1 tumor predisposition syndrome (BAP1-TPDS), characterized by high susceptibility to several tumor types, chiefly melanoma, mesothelioma, renal cell carcinoma, and basal cell carcinoma. Here, we present the results of our ten-year experience in the molecular diagnosis of BAP1-TPDS, along with a clinical update and cascade genetic testing of previously reported BAP1-TPDS patients and their relatives. Specifically, we sequenced germline DNA samples from 101 individuals with suspected BAP1-TPDS and validated pathogenic variants (PVs) by assessing BAP1 somatic loss in matching tumor specimens. Overall, we identified seven patients (7/101, 6.9%) carrying six different germline BAP1 PVs, including one novel variant. Consistently, cascade testing revealed a total of seven BAP1 PV carriers. In addition, we explored the mutational burden of BAP1-TPDS tumors by targeted next-generation sequencing. Lastly, we found that certain tumors present in PV carriers retain a wild-type BAP1 allele, suggesting a sporadic origin of these tumors or a functional role of heterozygous BAP1 in neoplastic development. Altogether, our findings have important clinical implications for therapeutic response of BAP1-TPDS patients.