Evolution of late-stage metastatic melanoma is dominated by aneuploidy and whole genome doubling

The Molecular Evolution of Melanoma Distant Metastases

The evolution of primary melanoma to lymph node and distant metastasis is incompletely understood. We examined the genomic diversity in melanoma progression in matched primary melanomas and lymph node and distant metastases from 17 patients. FISH analysis revealed cancer cell fractions with monotonic copy number alterations, including PHIP gain and PTEN loss, in the metastatic cascade. By contrast, the cancer cell fraction with copy number alterations for BPTF and MITF was reduced in lymph node metastases but increased in distant metastases. Separately, the cancer cell fraction with NCOA3 copy number alteration was comparable between primary tumors and lymph node metastases yet increased in distant metastases. These results suggest enrichment of the phosphoinositide 3-kinase and MITF pathways in the transition through the metastatic cascade. By contrast, next-generation sequencing analysis did not identify a consistent pattern of changes in variant allele frequency while revealing several intriguing findings, including decreased variant allele frequency in distant metastases and distinct drivers in lymph node versus distant metastases. These results provide evidence that distant melanoma metastasis does not always emanate from lymph node metastasis. These results enhance our understanding of clonal patterns of melanoma metastasis, with possible implications for targeted therapy and metastasis competency.

The intricate dance of tumor evolution: Exploring immune escape, tumor migration, drug resistance, and treatment strategies

Recent research has unveiled fascinating insights into the intricate mechanisms governing tumor evolution. These studies have illuminated how tumors adapt and proliferate by exploiting various factors, including immune evasion, resistance to therapeutic drugs, genetic mutations, and their ability to adapt to different environments. Furthermore, investigations into tumor heterogeneity and chromosomal aberrations have revealed the profound complexity that underlies the evolution of cancer. Emerging findings have also underscored the role of viral influences in the development and progression of cancer, introducing an additional layer of complexity to the field of oncology. Tumor evolution is a dynamic and complex process influenced by various factors, including immune evasion, drug resistance, tumor heterogeneity, and viral influences. Understanding these elements is indispensable for developing more effective treatments and advancing cancer therapies. A holistic approach to studying and addressing tumor evolution is crucial in the ongoing battle against cancer. The main goal of this comprehensive review is to explore the intricate relationship between tumor evolution and critical aspects of cancer biology. By delving into this complex interplay, we aim to provide a profound understanding of how tumors evolve, adapt, and respond to treatment strategies. This review underscores the pivotal importance of comprehending tumor evolution in shaping effective approaches to cancer treatment.

Comparison of Immunotherapy versus Targeted Therapy Effectiveness in BRAF-Mutant Melanoma Patients and Use of cGAS Expression and Aneuploidy as Potential Prognostic Biomarkers

BRAF-mutant melanoma patients can be treated with targeted therapy or immunotherapies, and it is not clear which should be provided first. Targeted treatments do not work in up to one-third of cases, while immunotherapies may only be effective in up to 60% and come with a high risk of immune-related side effects. Determining which treatment to provide first is thus of critical importance. Recent studies suggest that chromosomal instability and aneuploidy and cyclic GMP-AMP synthase (cGAS) can act as biomarkers for cancer severity and patient outcome. Neither potential biomarker has been extensively studied in melanoma. We examined 20 BRAF-mutant melanomas treated with immunotherapy or targeted therapy and measured chromosomal aneuploidy and cGAS expression levels. Treatment type, aneuploidy, and cGAS expression were correlated with progression-free survival (PFS) in these patients. Those treated with immunotherapy first had significantly better outcomes than those treated with targeted therapy, suggesting immunotherapy should be strongly considered as the first-line therapy for patients bearing BRAF-mutant melanoma. We found that there was no correlation of aneuploidy with outcome while there was some positive correlation of cGAS levels with PFS. Further studies are needed to confirm these findings and to test other potential biomarkers.

Aneuploid embryonic stem cells drive teratoma metastasis

Aneuploidy, a deviation of the chromosome number from euploidy, is one of the hallmarks of cancer. High levels of aneuploidy are generally correlated with metastasis and poor prognosis in cancer patients. However, the causality of aneuploidy in cancer metastasis remains to be explored. Here we demonstrate that teratomas derived from aneuploid murine embryonic stem cells (ESCs), but not from isogenic diploid ESCs, disseminated to multiple organs, for which no additional copy number variations were required. Notably, no cancer driver gene mutations were identified in any metastases. Aneuploid circulating teratoma cells were successfully isolated from peripheral blood and showed high capacities for migration and organ colonization. Single-cell RNA sequencing of aneuploid primary teratomas and metastases identified a unique cell population with high stemness that was absent in diploid ESCs-derived teratomas. Further investigation revealed that aneuploid cells displayed decreased proteasome activity and overactivated endoplasmic reticulum (ER) stress during differentiation, thereby restricting the degradation of proteins produced from extra chromosomes in the ESC state and causing differentiation deficiencies. Noticeably, both proteasome activator Oleuropein and ER stress inhibitor 4-PBA can effectively inhibit aneuploid teratoma metastasis.

Multistep tumor genetic evolution and changes in immunogenicity trigger immune-mediated disease eradication in stage IV melanoma: lessons from a single case

Durable remissions are observed in 10%-20% of treated patients with advanced metastatic melanoma but the factors associated with long-term complete clinical responses are largely unknown. Here, we report the molecular characteristics of tumor evolution during disease progression along a 9-year clinical course in a patient with advanced disseminated melanoma who received different treatments, including trametinib, ipilimumab, radiation, vemurafenib, surgical tumor debulking and a second ipilimumab course, ultimately achieving complete long-term disease remission.Longitudinal analyses of therapies-resistant metastatic tumors revealed the effects of different treatments on tumor's microenvironment and immunogenicity, ultimately creating a milieu favorable to immunotherapy response. Monitoring of the temporal dynamics of T cells by analysis of the T cell receptor (TCR) repertoire in the tumor and peripheral blood during disease evolution indicated that T-cell clones with common TCR rearrangements, present at low levels at baseline, were maintained and expanded after immunotherapy, and that TCR diversity increased. Analysis of genetic, molecular, and cellular components of the tumor depicted a multistep process in which treatment with kinase inhibitors strongly conditioned the immune microenvironment creating an inflamed milieu converting cold into hot tumors, while ipilimumab impacted and increased the TCR repertoire, a requirement for tumor rejection.Since the optimal sequencing of treatment with antibodies targeting immune checkpoints and kinase inhibitors for advanced melanoma is still clinically debated, this case indicates that immunotherapy success is possible even after progression on targeted therapy.

Whole-genome doubling in tissues and tumors

The overwhelming majority of proliferating somatic human cells are diploid, and this genomic state is typically maintained across successive cell divisions. However, failures in cell division can induce a whole-genome doubling (WGD) event, in which diploid cells transition to a tetraploid state. While some WGDs are developmentally programmed to produce nonproliferative tetraploid cells with specific cellular functions, unscheduled WGDs can be catastrophic: erroneously arising tetraploid cells are ill-equipped to cope with their doubled cellular and chromosomal content and quickly become genomically unstable and tumorigenic. Deciphering the genetics that underlie the genesis, physiology, and evolution of whole-genome doubled (WGD+) cells may therefore reveal therapeutic avenues to selectively eliminate pathological WGD+ cells.

Genetic Concordance in Primary Cutaneous Melanoma and Matched Metastasis: A Systematic Review and Meta-Analysis

Studying primary melanoma and its corresponding metastasis has twofold benefits. Firstly, to better understand tumor biology, and secondly, to determine which sample should be examined in assessing drug targets. This study systematically analyzed all the literature on primary melanoma and its matched metastasis. Following PRISMA guidelines, we searched multiple medical databases for relevant publications from January 2000 to December 2022, assessed the quality of the primary-level studies using the QUIPS tool, and summarized the concordance rate of the most reported genes using the random-effects model. Finally, we evaluated the inter-study heterogeneity using the subgroup analysis. Thirty-one studies investigated the concordance of BRAF and NRAS in 1220 and 629 patients, respectively. The pooled concordance rate was 89.4% [95% CI: 84.5; 93.5] for BRAF and 97.8% [95% CI: 95.8; 99.4] for NRAS. When high-quality studies were considered, only BRAF mutation status consistency increased. Five studies reported the concordance status of c-KIT (93%, 44 patients) and TERT promoter (64%, 53 patients). Lastly, three studies analyzed the concordance of cancer genes involved in the signaling pathways, apoptosis, and proliferation, such as CDKN2A (25%, four patients), TP53 (44%, nine patients), and PIK3CA (20%, five patients). Our study found that the concordance of known drug targets (mainly BRAF) during melanoma progression is higher than in previous meta-analyses, likely due to advances in molecular techniques. Furthermore, significant heterogeneity exists in the genes involved in the melanoma genetic makeup; although our results are based on small patient samples, more research is necessary for validation.

Validation of a microRNA liquid biopsy assay for diagnosis and risk stratification of invasive cutaneous melanoma

BACKGROUND

Noninvasive molecular biomarkers are needed to improve the early, accurate and precise diagnosis of invasive cutaneous melanoma.

OBJECTIVES

To independently validate a previously identified circulating microRNA signature of melanoma (MEL38), and, secondly, to develop a complementary microRNA signature, optimized for prognostication.

PATIENTS AND METHODS

MicroRNA expression profiling was performed on plasma samples from a multicentre observational case-control study, involving patients with primary or metastatic melanoma, melanoma in situ, nonmelanoma skin cancer, or benign naevi. MicroRNA profiles from patients with length of survival, treatment and sentinel lymph node biopsy (SLNB) data were used to develop the prognostic signature. The primary outcome of interest for MEL38 was its association with melanoma status, including area under the curve, binary diagnostic sensitivity and specificity, and incidence-adjusted positive and negative predictive values. The prognostic signature was assessed using rates of survival per risk group and relationship to conventional predictors of outcome.

RESULTS

Circulating microRNA profiles of 372 patients with invasive melanoma and 210 control individuals were generated. The average age of all participants was 59 years; 49% were male. A MEL38 score > 5.5 indicated the presence of invasive melanoma. Overall, 551/582 (95%) of patients were correctly diagnosed, with 93% sensitivity and 98% specificity. MEL38 score ranged from 0 to 10 with an area under the curve of 0.98 (95% confidence interval 0.97-0.99, P < 0.001). A novel prognostic 12-microRNA signature (MEL12) developed from 232 patients identified low-, standard- or high-risk groups, with 94%, 78% and 58% rates of 10-year melanoma-specific survival, respectively (log-rank P < 0.001). MEL12 prognostic risk groups were significantly associated with clinical staging (χ2, P < 0.001) and SLNB status (P = 0.027). Patients who were classified as high risk by MEL12 were approximately three times more likely to have melanoma detected in their sentinel lymph nodes compared to low-risk patients.

CONCLUSIONS

The circulating MEL38 signature may assist in diagnosing patients with invasive melanoma vs. other conditions associated with a lower - or negligible - risk of mortality. A complementary and prognostic MEL12 signature is predictive of SLNB status, clinical stage and probability of survival. Plasma microRNA profiling may help to optimize existing diagnostic pathways as well as enable personalized, risk-informed melanoma treatment decisions.

The copy number and mutational landscape of recurrent ovarian high-grade serous carcinoma

The drivers of recurrence and resistance in ovarian high grade serous carcinoma remain unclear. We investigate the acquisition of resistance by collecting tumour biopsies from a cohort of 276 women with relapsed ovarian high grade serous carcinoma in the BriTROC-1 study. Panel sequencing shows close concordance between diagnosis and relapse, with only four discordant cases. There is also very strong concordance in copy number between diagnosis and relapse, with no significant difference in purity, ploidy or focal somatic copy number alterations, even when stratified by platinum sensitivity or prior chemotherapy lines. Copy number signatures are strongly correlated with immune cell infiltration, whilst diagnosis samples from patients with primary platinum resistance have increased rates of CCNE1 and KRAS amplification and copy number signature 1 exposure. Our data show that the ovarian high grade serous carcinoma genome is remarkably stable between diagnosis and relapse and acquired chemotherapy resistance does not select for common copy number drivers.

Pathology and Molecular Biology of Melanoma

Almost every death in young patients with an advanced skin tumor is caused by melanoma. Today, with the help of modern treatments, these patients survive longer or can even achieve a cure. Advanced stage melanoma is frequently related with poor prognosis and physicians still find this disease difficult to manage due to the absence of a lasting response to initial treatment regimens and the lack of randomized clinical trials in post immunotherapy/targeted molecular therapy settings. New therapeutic targets are emerging from preclinical data on the genetic profile of melanocytes and from the identification of molecular factors involved in the pathogenesis of malignant transformation. In the current paper, we present the diagnostic challenges, molecular biology and genetics of malignant melanoma, as well as the current therapeutic options for patients with this diagnosis.

Whole-Genome Doubling as a source of cancer: how, when, where, and why?

Chromosome instability is a well-known hallmark of cancer, leading to increased genetic plasticity of tumoral cells, which favors cancer aggressiveness, and poor prognosis. One of the main sources of chromosomal instability are events that lead to a Whole-Genome Duplication (WGD) and the subsequently generated cell polyploidy. In recent years, several studies showed that WGD occurs at the early stages of cell transformation, which allows cells to later become aneuploid, thus leading to cancer progression. On the other hand, other studies convey that polyploidy plays a tumor suppressor role, by inducing cell cycle arrest, cell senescence, apoptosis, and even prompting cell differentiation, depending on the tissue cell type. There is still a gap in understanding how cells that underwent WGD can overcome the deleterious effect on cell fitness and evolve to become tumoral. Some laboratories in the chromosomal instability field recently explored this paradox, finding biomarkers that modulate polyploid cells to become oncogenic. This review brings a historical view of how WGD and polyploidy impact cell fitness and cancer progression, and bring together the last studies that describe the genes helping cells to adapt to polyploidy.

Late-Stage Metastatic Melanoma Emerges through a Diversity of Evolutionary Pathways

Understanding the evolutionary pathways to metastasis and resistance to immune-checkpoint inhibitors (ICI) in melanoma is critical for improving outcomes. Here, we present the most comprehensive intrapatient metastatic melanoma dataset assembled to date as part of the Posthumous Evaluation of Advanced Cancer Environment (PEACE) research autopsy program, including 222 exome sequencing, 493 panel-sequenced, 161 RNA sequencing, and 22 single-cell whole-genome sequencing samples from 14 ICI-treated patients. We observed frequent whole-genome doubling and widespread loss of heterozygosity, often involving antigen-presentation machinery. We found KIT extrachromosomal DNA may have contributed to the lack of response to KIT inhibitors of a KIT-driven melanoma. At the lesion-level, MYC amplifications were enriched in ICI nonresponders. Single-cell sequencing revealed polyclonal seeding of metastases originating from clones with different ploidy in one patient. Finally, we observed that brain metastases that diverged early in molecular evolution emerge late in disease. Overall, our study illustrates the diverse evolutionary landscape of advanced melanoma.

SIGNIFICANCE

Despite treatment advances, melanoma remains a deadly disease at stage IV. Through research autopsy and dense sampling of metastases combined with extensive multiomic profiling, our study elucidates the many mechanisms that melanomas use to evade treatment and the immune system, whether through mutations, widespread copy-number alterations, or extrachromosomal DNA. See related commentary by Shain, p. 1294. This article is highlighted in the In This Issue feature, p. 1275.

Melanoma Genomics: Shifting Focus from Inter- to Intrapatient Variation

SUMMARY

Traditionally, omic studies have prioritized the number of patients over the number of tumors per patient, but in a reversal of this study design, Spain and colleagues performed the largest intrapatient analysis of melanoma to date. Their work reveals mechanisms of treatment resistance, patterns of metastatic dissemination, and new insights into the evolutionary trajectories of melanoma. See related article by Spain et al., p. 1364 (1).

Insights into the metastatic cascade through research autopsies

Metastasis is a complex process and the leading cause of cancer-related death globally. Recent studies have demonstrated that genomic sequencing data from paired primary and metastatic tumours can be used to trace the evolutionary origins of cells responsible for metastasis. This approach has yielded new insights into the genomic alterations that engender metastatic potential, and the mechanisms by which cancer spreads. Given that the reliability of these approaches is contingent upon how representative the samples are of primary and metastatic tumour heterogeneity, we review insights from studies that have reconstructed the evolution of metastasis within the context of their cohorts and designs. We discuss the role of research autopsies in achieving the comprehensive sampling necessary to advance the current understanding of metastasis.

Targeting homologous recombination deficiency in uterine leiomyosarcoma

BACKGROUND

Uterine leiomyosarcoma (uLMS) is a rare and aggressive gynaecological malignancy, with individuals with advanced uLMS having a five-year survival of < 10%. Mutations in the homologous recombination (HR) DNA repair pathway have been observed in ~ 10% of uLMS cases, with reports of some individuals benefiting from poly (ADP-ribose) polymerase (PARP) inhibitor (PARPi) therapy, which targets this DNA repair defect. In this report, we screened individuals with uLMS, accrued nationally, for mutations in the HR repair pathway and explored new approaches to therapeutic targeting.

METHODS

A cohort of 58 individuals with uLMS were screened for HR Deficiency (HRD) using whole genome sequencing (WGS), whole exome sequencing (WES) or NGS panel testing. Individuals identified to have HRD uLMS were offered PARPi therapy and clinical outcome details collected. Patient-derived xenografts (PDX) were generated for therapeutic targeting.

RESULTS

All 13 uLMS samples analysed by WGS had a dominant COSMIC mutational signature 3; 11 of these had high genome-wide loss of heterozygosity (LOH) (> 0.2) but only two samples had a CHORD score > 50%, one of which had a homozygous pathogenic alteration in an HR gene (deletion in BRCA2). A further three samples harboured homozygous HRD alterations (all deletions in BRCA2), detected by WES or panel sequencing, with 5/58 (9%) individuals having HRD uLMS. All five individuals gained access to PARPi therapy. Two of three individuals with mature clinical follow up achieved a complete response or durable partial response (PR) with the subsequent addition of platinum to PARPi upon minor progression during initial PR on PARPi. Corresponding PDX responses were most rapid, complete and sustained with the PARP1-specific PARPi, AZD5305, compared with either olaparib alone or olaparib plus cisplatin, even in a paired sample of a BRCA2-deleted PDX, derived following PARPi therapy in the patient, which had developed PARPi-resistance mutations in PRKDC, encoding DNA-PKcs.

CONCLUSIONS

Our work demonstrates the value of identifying HRD for therapeutic targeting by PARPi and platinum in individuals with the aggressive rare malignancy, uLMS and suggests that individuals with HRD uLMS should be included in trials of PARP1-specific PARPi.

The Genomic Landscape of Melanoma and Its Therapeutic Implications

Melanoma is one of the most aggressive malignancies of the skin. The genetic composition of melanoma is complex and varies among different subtypes. With the aid of recent technologies such as next generation sequencing and single-cell sequencing, our understanding of the genomic landscape of melanoma and its tumor microenvironment has become increasingly clear. These advances may provide explanation to the heterogenic treatment outcomes of melanoma patients under current therapeutic guidelines and provide further insights to the development of potential new therapeutic targets. Here, we provide a comprehensive review on the genetics related to melanoma tumorigenesis, metastasis, and prognosis. We also review the genetics affecting the melanoma tumor microenvironment and its relation to tumor progression and treatment.

The journey from melanocytes to melanoma

Over the past decade, melanoma has led the field in new cancer treatments, with impressive gains in on-treatment survival but more modest improvements in overall survival. Melanoma presents heterogeneity and transcriptional plasticity that recapitulates distinct melanocyte developmental states and phenotypes, allowing it to adapt to and eventually escape even the most advanced treatments. Despite remarkable advances in our understanding of melanoma biology and genetics, the melanoma cell of origin is still fiercely debated because both melanocyte stem cells and mature melanocytes can be transformed. Animal models and high-throughput single-cell sequencing approaches have opened new opportunities to address this question. Here, we discuss the melanocytic journey from the neural crest, where they emerge as melanoblasts, to the fully mature pigmented melanocytes resident in several tissues. We describe a new understanding of melanocyte biology and the different melanocyte subpopulations and microenvironments they inhabit, and how this provides unique insights into melanoma initiation and progression. We highlight recent findings on melanoma heterogeneity and transcriptional plasticity and their implications for exciting new research areas and treatment opportunities. The lessons from melanocyte biology reveal how cells that are present to protect us from the damaging effects of ultraviolet radiation reach back to their origins to become a potentially deadly cancer.

Whole-Exome Sequencing and cfDNA Analysis Uncover Genetic Determinants of Melanoma Therapy Response in a Real-World Setting

Although several studies have explored the molecular landscape of metastatic melanoma, the genetic determinants of therapy resistance are still largely unknown. Here, we aimed to determine the contribution of whole-exome sequencing and circulating free DNA (cfDNA) analysis in predicting response to therapy in a consecutive real-world cohort of 36 patients, undergoing fresh tissue biopsy and followed during treatment. Although the underpowered sample size limited statistical analysis, samples from non-responders had higher copy number variations and mutations in melanoma driver genes compared to responders in the BRAF V600+ subset. In the BRAF V600- subset, Tumor Mutational Burden (TMB) was twice that in responders vs. non-responders. Genomic layout revealed commonly known and novel potential intrinsic/acquired resistance driver gene variants. Among these, RAC1, FBXW7, GNAQ mutations, and BRAF/PTEN amplification/deletion were present in 42% and 67% of patients, respectively. Both Loss of Heterozygosity (LOH) load and tumor ploidy were inversely associated with TMB. In immunotherapy-treated patients, samples from responders showed higher TMB and lower LOH and were more frequently diploid compared to non-responders. Secondary germline testing and cfDNA analysis proved their efficacy in finding germline predisposing variants carriers (8.3%) and following dynamic changes during treatment as a surrogate of tissue biopsy, respectively.

Genomic Profiling of Metastatic Basal cell Carcinoma Reveals Candidate Drivers of Disease and Therapeutic Targets

Basal cell carcinomas (BCCs) are human beings' most common malignant tumors. Most are easily managed by surgery or topical therapies, and metastasis is rare. Although BCCs can become locally advanced, metastatic BCCs are very uncommon and may be biologically distinct. We assessed the clinicopathologic characteristics of 17 patients with metastatic BCC and pursued whole-exome sequencing of tumor and germline DNA from 8 patients. Genomic profiling revealed aberrant activation of Hedgehog signaling and alterations in GLI transcriptional regulators and Notch and Hippo signaling. Matched local recurrences of primary BCCs and metastases from 3 patients provided evidence of a clonal origin in all cases. Mutations associated with YAP inhibition were found exclusively in 2 hematogenously-spread lung metastases, and metastatic BCCs were enriched for mutations in the YAP/TAZ-binding domain of TEAD genes. Accordingly, YAP/TAZ nuclear localization was associated with metastatic types and Hippo mutations, suggesting an enhanced oncogenic role in hematogenously-spread metastases. Mutations in RET, HGF, and phosphatidylinositol 3‑kinase (PI3K)/protein kinase B (AKT) signaling were enriched compared with a cohort of low clinical-risk BCCs. Our results implicate Hippo and PI3K/AKT dysregulation in metastatic progression of BCCs, making these potential therapeutic targets in metastatic disease. The common clonal origin of matched recurrent and metastatic BCCs suggests that molecular profiling can assist in determining the nature/origin of poorly differentiated metastatic tumors of uncertain type. Genes and pathways enriched for mutations in this cohort are candidate drivers of metastasis and can be used to identify patients at high risk of metastasis who may benefit from aggressive local treatment and careful clinical follow-up.

Comparative Genomics Provides Etiologic and Biological Insight into Melanoma Subtypes

Melanoma is a cancer of melanocytes, with multiple subtypes based on body site location. Cutaneous melanoma is associated with skin exposed to ultraviolet radiation; uveal melanoma occurs in the eyes; mucosal melanoma occurs in internal mucous membranes; and acral melanoma occurs on the palms, soles, and nail beds. Here, we present the largest whole-genome sequencing study of melanoma to date, with 570 tumors profiled, as well as methylation and RNA sequencing for subsets of tumors. Uveal melanoma is genomically distinct from other melanoma subtypes, harboring the lowest tumor mutation burden and with significantly mutated genes in the G-protein signaling pathway. Most cutaneous, acral, and mucosal melanomas share alterations in components of the MAPK, PI3K, p53, p16, and telomere pathways. However, the mechanism by which these pathways are activated or inactivated varies between melanoma subtypes. Additionally, we identify potential novel germline predisposition genes for some of the less common melanoma subtypes.

SIGNIFICANCE

This is the largest whole-genome analysis of melanoma to date, comprehensively comparing the genomics of the four major melanoma subtypes. This study highlights both similarities and differences between the subtypes, providing insights into the etiology and biology of melanoma. This article is highlighted in the In This Issue feature, p. 2711.

From Samples to Germline and Somatic Sequence Variation: A Focus on Next-Generation Sequencing in Melanoma Research

Next-generation sequencing (NGS) applications have flourished in the last decade, permitting the identification of cancer driver genes and profoundly expanding the possibilities of genomic studies of cancer, including melanoma. Here we aimed to present a technical review across many of the methodological approaches brought by the use of NGS applications with a focus on assessing germline and somatic sequence variation. We provide cautionary notes and discuss key technical details involved in library preparation, the most common problems with the samples, and guidance to circumvent them. We also provide an overview of the sequence-based methods for cancer genomics, exposing the pros and cons of targeted sequencing vs. exome or whole-genome sequencing (WGS), the fundamentals of the most common commercial platforms, and a comparison of throughputs and key applications. Details of the steps and the main software involved in the bioinformatics processing of the sequencing results, from preprocessing to variant prioritization and filtering, are also provided in the context of the full spectrum of genetic variation (SNVs, indels, CNVs, structural variation, and gene fusions). Finally, we put the emphasis on selected bioinformatic pipelines behind (a) short-read WGS identification of small germline and somatic variants, (b) detection of gene fusions from transcriptomes, and (c) de novo assembly of genomes from long-read WGS data. Overall, we provide comprehensive guidance across the main methodological procedures involved in obtaining sequencing results for the most common short- and long-read NGS platforms, highlighting key applications in melanoma research.

Molecular Mechanisms Driving the Formation of Brain Metastases

Simple Summary

Brain metastases are the most common brain tumor in adults and are associated with poor prognosis. The propensity of different solid tumors to metastasize varies greatly, with lung, breast, and melanoma primary tumors commonly leading to brain metastases, while other primaries such as prostate rarely metastasize to the brain. The molecular mechanisms that predispose and facilitate brain metastasis development are poorly understood. In this review, we present the current data on the genomic landscape of brain metastases that arise from various primary cancers and also outline potential molecular mechanisms that drive the formation of distant metastases in the brain.

Abstract

Targeted therapies for cancers have improved primary tumor response rates, but concomitantly, brain metastases (BM) have become the most common brain tumors in adults and are associated with a dismal prognosis of generally less than 6 months, irrespective of the primary cancer type. They most commonly occur in patients with primary breast, lung, or melanoma histologies; however, they also appear in patients with other primary cancers including, but not limited to, prostate cancer, colorectal cancer, and renal cell carcinoma. Historically, molecular biomarkers have normally been identified from primary tumor resections. However, clinically informative genomic alterations can occur during BM development and these potentially actionable alterations are not always detected in the primary tumor leading to missed opportunities for effective targeted therapy. The molecular mechanisms that facilitate and drive metastasis to the brain are poorly understood. Identifying the differences between the brain and other extracranial sties of metastasis, and between primary tumors and BM, is essential to improving our understanding of BM development and ultimately patient management and survival. In this review, we present the current data on the genomic landscape of BM from various primary cancers which metastasize to the brain and outline potential mechanisms which may play a role in promoting the formation of the distant metastases in the brain.

Causes, consequences and clinical significance of aneuploidy across melanoma subtypes

Aneuploidy, the state of the cell in which the number of whole chromosomes or chromosome arms becomes imbalanced, has been recognized as playing a pivotal role in tumor evolution for over 100 years. In melanoma, the extent of aneuploidy, as well as the chromosomal regions that are affected differ across subtypes, indicative of distinct drivers of disease. Multiple studies have suggested a role for aneuploidy in diagnosis and prognosis of melanomas, as well as in the context of immunotherapy response. A number of key constituents of the cell cycle have been implicated in aneuploidy acquisition in melanoma, including several driver mutations. Here, we review the state of the art on aneuploidy in different melanoma subtypes, discuss the potential drivers, mechanisms underlying aneuploidy acquisition as well as its value in patient diagnosis, prognosis and response to immunotherapy treatment.

Association of Circular RNA and Long Non-Coding RNA Dysregulation with the Clinical Response to Immune Checkpoint Blockade in Cutaneous Metastatic Melanoma

Cutaneous melanoma (CM) is the most lethal form of skin cancer if it becomes metastatic, where treatment options and survival chances decrease dramatically. Immunotherapy treatments based on the immunologic checkpoint inhibitors programmed death cell protein 1 (PD-1) and cytotoxic T-lymphocyte antigen 4 (CTLA-4) constituted a main breakthrough in the treatment of metastatic CM, particularly for the achievement of long-term benefits. Even though it is a very promising therapy, resistance to primary immune checkpoint blockade (ICB) arises in about 70% of CM patients treated with a CTLA-4 inhibitor, and 40–65% of CM patients administered with a PD-1-targeting treatment. Some long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) are implicated in triggering pro- and anti-tumorigenic responses to various cancer treatments. The relationship between lncRNAs, circRNAs and ICB immunotherapy has not been explored in cutaneous metastatic melanoma (CMM). The aim of this pilot study is to evaluate the potential role of circRNA and lncRNA expression as pre-treatment predictors of the clinical response to immunotherapy in CMM patients. RNA-seq from 12 formalin-fixed paraffin-embedded (FFPE) samples from the metastatic biopsies of CMM patients treated with nivolumab was used to identify response-associated transcripts. Our findings indicate that specific lncRNAs and circRNAs, probably acting as competitive endogenous RNAs (ceRNAs), are involved in the regulatory networks of the immune response against metastatic melanoma that these patients have under treatment with nivolumab. Moreover, we established a risk score that yields predictions of the overall survival (OS) and progression-free survival (PFS) of CMM patients with high accuracy. This proof-of-principle work provides a possible insight into the function of ceRNAs, contributing to efforts to decipher the complex molecular mechanisms of ICB cancer treatment response.

Gain of Chromosome 1q Perturbs a Competitive Endogenous RNA Network to Promote Melanoma Metastasis

Somatic copy-number alterations (CNA) promote cancer, but the underlying driver genes may not be comprehensively identified if only the functions of the encoded proteins are considered. mRNAs can act as competitive endogenous RNAs (ceRNA), which sponge miRNAs to posttranscriptionally regulate gene expression in a protein coding-independent manner. We investigated the contribution of ceRNAs to the oncogenic effects of CNAs. Chromosome 1q gains promoted melanoma progression and metastasis at least in part through overexpression of three mRNAs with ceRNA activity: CEP170, NUCKS1, and ZC3H11A. These ceRNAs enhanced melanoma metastasis by sequestering tumor suppressor miRNAs. Orthogonal genetic assays with miRNA inhibitors and target site blockers, along with rescue experiments, demonstrated that miRNA sequestration is critical for the oncogenic effects of CEP170, NUCKS1, and ZC3H11A mRNAs. Furthermore, chromosome 1q ceRNA-mediated miRNA sequestration alleviated the repression of several prometastatic target genes. This regulatory RNA network was evident in other cancer types, suggesting chromosome 1q ceRNA deregulation as a common driver of cancer progression. Taken together, this work demonstrates that ceRNAs mediate the oncogenicity of somatic CNAs.

SIGNIFICANCE

The function of CEP170, NUCKS1, and ZC3H11A mRNAs as competitive endogenous RNAs that sequester tumor suppressor microRNAs underlies the oncogenic activity of chromosome 1q gains.

UV-Induced Somatic Mutations Driving Clonal Evolution in Healthy Skin, Nevus, and Cutaneous Melanoma

Introduction: Due to its aggressiveness, cutaneous melanoma (CM) is responsible for most skin cancer-related deaths worldwide. The origin of CM is closely linked to the appearance of UV-induced somatic mutations in melanocytes present in normal skin or in CM precursor lesions (nevi or dysplastic nevi). In recent years, new NGS studies performed on CM tissue have increased the understanding of the genetic somatic changes underlying melanomagenesis and CM tumor progression. Methods: We reviewed the literature using all important scientific databases. All articles related to genomic mutations in CM as well as normal skin and nevi were included, in particular those related to somatic mutations produced by UV radiation. Conclusions: CM development and progression are strongly associated with exposure to UV radiation, although each melanoma subtype has different characteristic genetic alterations and evolutionary trajectories. While BRAF and NRAS mutations are common in the early stages of tumor development for most CM subtypes, changes in CDKN2A, TP53 and PTEN, together with TERT promoter mutations, are especially common in advanced stages. Additionally, large genome duplications, loss of heterozygosity, and copy number variations are hallmarks of metastatic disease. Finally, the mutations driving melanoma targeted-therapy drug resistance are also summarized. The complete sequential stages of clonal evolution leading to CM onset from normal skin or nevi are still unknown, so further studies are needed in this field to shed light on the molecular pathways involved in CM malignant transformation and in melanoma acquired drug resistance.

Oncogenic BRAF induces whole-genome doubling through suppression of cytokinesis

Melanomas and other solid tumors commonly have increased ploidy, with near-tetraploid karyotypes being most frequently observed. Such karyotypes have been shown to arise through whole-genome doubling events that occur during early stages of tumor progression. The generation of tetraploid cells via whole-genome doubling is proposed to allow nascent tumor cells the ability to sample various pro-tumorigenic genomic configurations while avoiding the negative consequences that chromosomal gains or losses have in diploid cells. Whereas a high prevalence of whole-genome doubling events has been established, the means by which whole-genome doubling arises is unclear. Here, we find that BRAFV600E, the most common mutation in melanomas, can induce whole-genome doubling via cytokinesis failure in vitro and in a zebrafish melanoma model. Mechanistically, BRAFV600E causes decreased activation and localization of RhoA, a critical cytokinesis regulator. BRAFV600E activity during G1/S phases of the cell cycle is required to suppress cytokinesis. During G1/S, BRAFV600E activity causes inappropriate centriole amplification, which is linked in part to inhibition of RhoA and suppression of cytokinesis. Together these data suggest that common abnormalities of melanomas linked to tumorigenesis - amplified centrosomes and whole-genome doubling events - can be induced by oncogenic BRAF and other mutations that increase RAS/MAPK pathway activity.

Inactivation of the Hippo tumor suppressor pathway promotes melanoma

Melanoma is commonly driven by activating mutations in the MAP kinase BRAF; however, oncogenic BRAF alone is insufficient to promote melanomagenesis. Instead, its expression induces a transient proliferative burst that ultimately ceases with the development of benign nevi comprised of growth-arrested melanocytes. The tumor suppressive mechanisms that restrain nevus melanocyte proliferation remain poorly understood. Here we utilize cell and murine models to demonstrate that oncogenic BRAF leads to activation of the Hippo tumor suppressor pathway, both in melanocytes in vitro and nevus melanocytes in vivo. Mechanistically, we show that oncogenic BRAF promotes both ERK-dependent alterations in the actin cytoskeleton and whole-genome doubling events, which independently reduce RhoA activity to promote Hippo activation. We also demonstrate that functional impairment of the Hippo pathway enables oncogenic BRAF-expressing melanocytes to bypass nevus formation and rapidly form melanomas. Our data reveal that the Hippo pathway enforces the stable arrest of nevus melanocytes and represents a critical barrier to melanoma development.

Activating mutations of BRAF alone are inadequate to drive melanoma formation. Here the authors show that activation of Hippo signalling by oncogenic BRAF represents an additional safeguard to limit BRAF-dependent human melanocyte growth and melanoma formation.

Heterogeneity in Melanoma

There is growing evidence that tumour heterogeneity has an imperative role in cancer development, evolution and resistance to therapy. Continuing advancements in biomedical research enable tumour heterogeneity to be observed and studied more critically. As one of the most heterogeneous human cancers, melanoma displays a high level of biological complexity during disease progression. However, much is still unknown regarding melanoma tumour heterogeneity, as well as the role it plays in disease progression and treatment response. This review aims to provide a concise summary of the importance of tumour heterogeneity in melanoma.

The Genetic Evolution of Metastasis

Cancer is an evolutionary process that is characterized by the emergence of multiple genetically distinct populations or clones within the primary tumor. Intratumor heterogeneity provides a substrate for the selection of adaptive clones, such as those that lead to metastasis. Comparative molecular studies of primary tumors and metastases have identified distinct genomic features associated with the development of metastases. In this review, we discuss how these insights could inform clinical decision-making and uncover rational antimetastasis treatment strategies.

Unravelling the tumour genome: the evolutionary and clinical impacts of structural variants in tumourigenesis

Structural variants (SVs) represent a major source of aberration in tumour genomes. Given the diversity in the size and type of SVs present in tumours, the accurate detection and interpretation of SVs in tumours is challenging. New classes of complex structural events in tumours are discovered frequently, and the definitions of the genomic consequences of complex events are constantly being refined. Detailed analyses of short‐read whole‐genome sequencing (WGS) data from large tumour cohorts facilitate the interrogation of SVs at orders of magnitude greater scale and depth. However, the inherent technical limitations of short‐read WGS prevent us from accurately detecting and investigating the impact of all the SVs present in tumours. The expanded use of long‐read WGS will be critical for improving the accuracy of SV detection, and in fully resolving complex SV events, both of which are crucial for determining the impact of SVs on tumour progression and clinical outcome. Despite the present limitations, we demonstrate that SVs play an important role in tumourigenesis. In particular, SVs contribute significantly to late‐stage tumour development and to intratumoural heterogeneity. The evolutionary trajectories of SVs represent a window into the clonal dynamics in tumours, a comprehensive understanding of which will be vital for influencing patient outcomes in the future. Recent findings have highlighted many clinical applications of SVs in cancer, from early detection to biomarkers for treatment response and prognosis. As the methods to detect and interpret SVs improve, elucidating the full breadth of the complex SV landscape and determining how these events modulate tumour evolution will improve our understanding of cancer biology and our ability to capitalise on the utility of SVs in the clinical management of cancer patients. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Molecular characterization of fast-growing melanomas

BACKGROUND

The rate of growth of primary melanoma is a robust predictor of aggressiveness, but the mutational profile of fast-growing melanomas (FGMM) and the potential to stratify patients at high risk of death has not been comprehensively studied.

OBJECTIVE

To investigate the epidemiologic, clinical, and mutational profile of primary cutaneous melanomas with a thickness ≥ 1 mm, stratified by rate of growth.

METHODS

Observational prospective study. Deep-targeted sequencing of 40 melanoma driver genes on formalin fixed, paraffin-embedded primary melanoma samples. Comparison of FGMM (rate of growth > 0.5 mm/month) and nonFGMM (rate of growth ≤ 0.5 mm/month).

RESULTS

Two hundred patients were enrolled, among wom 70 had FGMM. The relapse-free survival was lower in the FGMM group (P = .014). FGMM had a higher number of predicted deleterious mutations within the 40 genes than nonFGMM (P = .033). Ulceration (P = .032), thickness (P = .006), lower sun exposure (P = .049), and fibroblast growth factor receptor 2 (FGFR2) mutations (P = .037) were significantly associated with fast growth.

LIMITATIONS

Single-center study, cohort size, potential memory bias, number of investigated genes.

CONCLUSION

Fast growth is linked to specific tumor biology and environmental factors. Ulceration, thickness, and FGFR2 mutations are associated with fast growth. Screening for FGFR2 mutations might provide an additional tool to better identify FGMM, which are probably good candidates for adjuvant therapies.

Genomic and Transcriptomic Underpinnings of Melanoma Genesis, Progression, and Metastasis

Melanoma is a deadly skin cancer with rapidly increasing incidence worldwide. The discovery of the genetic drivers of melanomagenesis in the last decade has led the World Health Organization to reclassify melanoma subtypes by their molecular pathways rather than traditional clinical and histopathologic features. Despite this significant advance, the genomic and transcriptomic drivers of metastatic progression are less well characterized. This review describes the known molecular pathways of cutaneous and uveal melanoma progression, highlights recently identified pathways and mediators of metastasis, and touches on the influence of the tumor microenvironment on metastatic progression and treatment resistance. While targeted therapies and immune checkpoint blockade have significantly aided in the treatment of advanced disease, acquired drug resistance remains an unfortunately common problem, and there is still a great need to identify potential prognostic markers and novel therapeutic targets to aid in such cases.

BRAFV600E induces reversible mitotic arrest in human melanocytes via microrna-mediated suppression of AURKB

Benign melanocytic nevi frequently emerge when an acquired BRAFV600E mutation triggers unchecked proliferation and subsequent arrest in melanocytes. Recent observations have challenged the role of oncogene-induced senescence in melanocytic nevus formation, necessitating investigations into alternative mechanisms for the establishment and maintenance of proliferation arrest in nevi. We compared the transcriptomes of melanocytes from healthy human skin, nevi, and melanomas arising from nevi and identified a set of microRNAs as highly expressed nevus-enriched transcripts. Two of these microRNAs-MIR211-5p and MIR328-3p-induced mitotic failure, genome duplication, and proliferation arrest in human melanocytes through convergent targeting of AURKB. We demonstrate that BRAFV600E induces a similar proliferation arrest in primary human melanocytes that is both reversible and conditional. Specifically, BRAFV600E expression stimulates either arrest or proliferation depending on the differentiation state of the melanocyte. We report genome duplication in human melanocytic nevi, reciprocal expression of AURKB and microRNAs in nevi and melanomas, and rescue of arrested human nevus cells with AURKB expression. Taken together, our data describe an alternative molecular mechanism for melanocytic nevus formation that is congruent with both experimental and clinical observations.

High-Dimensional Single-Cell Transcriptomics in Melanoma and Cancer Immunotherapy

Recent advances in single-cell transcriptomics have greatly improved knowledge of complex transcriptional programs, rapidly expanding our knowledge of cellular phenotypes and functions within the tumour microenvironment and immune system. Several new single-cell technologies have been developed over recent years that have enabled expanded understanding of the mechanistic cells and biological pathways targeted by immunotherapies such as immune checkpoint inhibitors, which are now routinely used in patient management with high-risk early-stage or advanced melanoma. These technologies have method-specific strengths, weaknesses and capabilities which need to be considered when utilising them to answer translational research questions. Here, we provide guidance for the implementation of single-cell transcriptomic analysis platforms by reviewing the currently available experimental and analysis workflows. We then highlight the use of these technologies to dissect the tumour microenvironment in the context of cancer patients treated with immunotherapy. The strategic use of single-cell analytics in clinical settings are discussed and potential future opportunities are explored with a focus on their use to rationalise the design of novel immunotherapeutic drug therapies that will ultimately lead to improved cancer patient outcomes.

Squaring the circle: circRNAs in melanoma

Non-coding RNAs are emerging as critical molecules in the genesis, progression, and therapy resistance of cutaneous melanoma. This includes circular RNAs (circRNAs), a class of non-coding RNAs with distinct characteristics that forms through non-canonical back-splicing. In this review, we summarize the features and functions of circRNAs and introduce the current knowledge of the roles of circRNAs in melanoma. We also highlight the various mechanisms of action of the well-studied circRNA CDR1as and describe how it acts as a melanoma tumor suppressor. We further discuss the utility of circRNAs as biomarkers, therapeutic targets, and therapeutic agents in melanoma and outline challenges that must be overcome to comprehensively characterize circRNA functions.

Custom workflows to improve joint variant calling from multiple related tumour samples: FreeBayesSomatic and Strelka2Pass

SUMMARY

This work describes two novel workflows for variant calling that extend the widely used algorithms of Strelka2 and FreeBayes to call somatic mutations from multiple related tumour samples and one matched normal sample. We show that these workflows offer higher precision and recall than their single tumour-normal pair equivalents in both simulated and clinical sequencing data.

AVAILABILITY AND IMPLEMENTATION

Source code freely available at the following link: https://atlassian.petermac.org.au/bitbucket/projects/DAW/repos/multisamplevariantcalling and executable through Janis (https://github.com/PMCC-BioinformaticsCore/janis) under the GPLv3 licence.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.