Clonal architectures and driver mutations in metastatic melanomas

Targeted DNA Sequencing of Cutaneous Melanoma Identifies Prognostic and Predictive Alterations

BACKGROUND

Cutaneous melanoma (CM) can be molecularly classified into four groups: BRAF mutant, NRAS mutant, NF1 mutant and triple wild-type (TWT) tumors lacking any of these three alterations. In the era of immune checkpoint inhibition (ICI) and targeted molecular therapy, the clinical significance of these groups remains unclear. Here, we integrate targeted DNA sequencing with comprehensive clinical follow-up in CM patients.

METHODS

This was a retrospective cohort study that assessed clinical and molecular features from patients with localized or metastatic CM who underwent targeted next-generation sequencing as part of routine clinical care. A total of 254 patients with CM who had a CLIA-certified targeted sequencing assay performed on their tumor tissue were included.

RESULTS

Of the 254 patients with cutaneous melanoma, 77 were BRAF mutant (30.3%), 77 were NRAS mutant (30.3%), 47 were NF1 mutant (18.5%), 33 were TWT (13.0%) and the remaining 20 (7.9%) carried mutations in multiple driver genes (BRAF/NRAS/NF1 co-mutated). The majority of this co-mutation group carried mutations in NF1 (n = 19 or 90%) with co-occurring mutations in BRAF or NRAS, often with a weaker oncogenic variant. Consistently, NF1 mutant tumors harbored numerous significantly co-altered genes compared to BRAF or NRAS mutant tumors. The majority of TWT tumors (n = 29, 87.9%) harbor a pathogenic mutation within a known Ras/MAPK signaling pathway component. Of the 154 cases with available TMB data, the median TMB was 20 (range 0.7-266 mutations/Mb). A total of 14 cases (9.1%) were classified as having a low TMB (≤5 mutations/Mb), 64 of 154 (41.6%) had an intermediate TMB (>5 and ≤20 mutations/Mb), 40 of 154 (26.0%) had a high TMB (>20 and ≤50 mutations/Mb) and 36 of 154 (23.4%) were classified as having a very high TMB (>50 mutations/Mb). NRAS mutant melanoma demonstrated significantly decreased overall survival on multivariable analysis (HR for death 2.95, 95% CI 1.13-7.69, p = 0.027, log-rank test) compared with other TCGA molecular subgroups. Of the 116 patients in our cohort with available treatment data, 36 received a combination of dual ICI with anti-CTLA4 and anti-PD1 inhibition as first-line therapy. Elevated TMB was associated with significantly longer progression-free survival following dual-agent ICI (HR 0.26, 95% CI 0.07-0.90, p = 0.033, log-rank test).

CONCLUSIONS

NRAS mutation in CMs correlated with significantly worse overall survival. Elevated TMB was associated with increased progression-free survival for patients treated with a combination of dual ICI, supporting the potential utility of TMB as a predictive biomarker for ICI response in melanoma.

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.

Crosstalk in skin: Loss of desmoglein 1 in keratinocytes inhibits BRAFV600E-induced cellular senescence in human melanocytes

Melanoma arises from transformation of melanocytes in the basal layer of the epidermis where they are surrounded by keratinocytes, with which they interact through cell contact and paracrine communication. Considerable effort has been devoted to determining how the accumulation of oncogene and tumor suppressor gene mutations in melanocytes drive melanoma development. However, the extent to which alterations in keratinocytes that occur in the developing tumor niche serve as extrinsic drivers of melanoma initiation and progression is poorly understood. We recently identified the keratinocyte-specific cadherin, desmoglein 1 (Dsg1), as an important mediator of keratinocyte:melanoma cell crosstalk, demonstrating that its chronic loss, which can occur through melanoma cell-dependent paracrine signaling, promotes behaviors that mimic a malignant phenotype. Here we address the extent to which Dsg1 loss affects early steps in melanomagenesis. RNA-Seq analysis revealed that paracrine signals from Dsg1-deficient keratinocytes mediate a transcriptional switch from a differentiated to undifferentiated cell state in melanocytes expressing BRAFV600E, a driver mutation commonly present in both melanoma and benign nevi and reported to cause growth arrest and oncogene-induced senescence (OIS). Of ~220 differentially expressed genes in BRAFV600E cells treated with Dsg1-deficient conditioned media (CM), the laminin superfamily member NTN4/Netrin-4, which inhibits senescence in endothelial cells, stood out. Indeed, while BRAFV600E melanocytes treated with Dsg1-deficient CM showed signs of senescence bypass as assessed by increased senescence-associated β-galactosidase activity and decreased p16, knockdown of NTN4 reversed these effects. These results suggest that Dsg1 loss in keratinocytes provides an extrinsic signal to push melanocytes towards oncogenic transformation once an initial mutation has been introduced.

Oncogenic alterations reveal key strategies for precision oncology in melanoma treatment

Background

Molecular profiling with next-generation sequencing (NGS) has been applied in multiple solid tumors, including melanomas, to identify potential drug targets. However, the association between clinical outcomes and the molecular alterations has not yet been fully clarified.

Methods

A total of 108 patients with melanoma were included in this study, 95 of whom had both sequencing data and clinical outcomes were collected. We analyzed the genetic alterations of 108 malignant melanoma patients using the OncoCare panel, which covers 559 genes.

Results

A model was also established to predict side effects through a combination analysis of clinical data and somatic variants, yielding an area under the receiver operating characteristic curve (AUROC) score of 0.8. We also identified epidermal growth factor receptor (EGFR) mutation was excellent predictor for progression-free survival (PFS) for patient who received immunotherapy (log-rank P=0.01), while tumor mutation burden (TMB) was found to not be significantly associated with PFS (log-rank P=0.87). Combining clinical features with genetic analysis, we found that patients carrying both DNA POLD1/ALOX12B or POLD1/PTPRT mutations had a significantly lower survival rate.

Conclusions

Overall, these results demonstrate the benefits of applying NGS clinical panels and shed light on future directions of personalized therapeutics for the treatment of melanoma.

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.

Nanoparticle-Based Combination Therapy for Melanoma

Melanoma is a cutaneous carcinoma, and its incidence is rapidly increasing with every year. The treatment options for melanoma have been comprehensively studied. Conventional treatment methods (e.g., radiotherapy, chemotherapy and photodynamic therapy) with surgical removal inevitably cause serious complications; moreover, resistance is common. Nanoparticles (NPs) combined with conventional methods are new and promising options to treat melanoma, and many combinations have been achieving good success. Due to their physical and biological features, NPs can help target intended melanoma cells more efficiently with less damage. This creates new hope for a better treatment strategy for melanoma with minimum damage and maximum efficacy.

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.

SWI/SNF Chromatin Remodeling Enzymes in Melanoma

Melanoma is an aggressive malignancy that arises from the transformation of melanocytes on the skin, mucosal membranes, and uvea of the eye. SWI/SNF chromatin remodeling enzymes are multi-subunit complexes that play important roles in the development of the melanocyte lineage and in the response to ultraviolet radiation, a key environmental risk factor for developing cutaneous melanoma. Exome sequencing has revealed frequent loss of function mutations in genes encoding SWI/SNF subunits in melanoma. However, some SWI/SNF subunits have also been demonstrated to have pro-tumorigenic roles in melanoma and to affect sensitivity to therapeutics. This review summarizes studies that have implicated SWI/SNF components in melanomagenesis and have evaluated how SWI/SNF subunits modulate the response to current therapeutics.

Association of PTPRT mutations with immune checkpoint inhibitors response and outcome in melanoma and non-small cell lung cancer

Purpose

Protein tyrosine phosphatase receptor type T (PTPRT), which is a well‐known phosphatase and mutates frequently in melanoma and non‐small cell lung cancer (NSCLC). Our research aims to elucidate its mutation association with immune checkpoint inhibitors (ICI) efficacy.

Methods

We integrated whole‐exome sequencing (WES)‐based somatic mutation profiles and clinical characteristics of 631 melanoma samples received ICI agents from eight studies and 109 NSCLC samples from two studies. For validation, 321 melanoma and 350 NSCLC immunotherapy samples with targeted next‐generation sequencing (NGS) were employed. Besides, an independent NSCLC cohort contained 240 samples was also collected for further corroboration. Distinct immune infiltration was evaluated according to the PTPRT mutational status.

Results

In the WES melanoma cohort, patients with PTPRT mutations harbored a significantly elevated ICI response rate (40.5% vs. 28.6%, p = 0.036) and a prolonged survival outcome (35.3 vs. 24.9 months, p = 0.006). In the WES NSCLC cohort, the favorable response and immunotherapy survival were also observed in PTPRT‐mutated patients (p = 0.036 and 0.019, respectively). For the validation cohorts, the associations of PTRPT mutations with better prognoses were identified in melanoma, NSCLC, and pan‐cancer patients with targeted‐NGS (all p < 0.05). Moreover, immunology analyses showed the higher mutation burden, increased lymphocyte infiltration, decreased‐ activated‐stroma, and immune response pathways were detected in patients with PTPRT mutations.

Conclusion

Our investigation indicates that PTPRT mutations may be considered as a potential indicator for assessing ICI efficacy in melanoma and NSCLC, even across multiple cancers. Further prospective validation cohorts are warranted.

MUG Mel3 Cell Lines Reflect Heterogeneity in Melanoma and Represent a Robust Model for Melanoma in Pregnancy

Melanomas are aggressive tumors with a high metastatic potential and an increasing incidence rate. They are known for their heterogeneity and propensity to easily develop therapy-resistance. Nowadays they are one of the most common cancers diagnosed during pregnancy. Due to the difficulty in balancing maternal needs and foetal safety, melanoma is challenging to treat. The aim of this study was to provide a potential model system for the study of melanoma in pregnancy and to illustrate melanoma heterogeneity. For this purpose, a pigmented and a non-pigmented section of a lymph node metastasis from a pregnant patient were cultured under different conditions and characterized in detail. All four culture conditions exhibited different phenotypic, genotypic as well as tumorigenic properties, and resulted in four newly established melanoma cell lines. To address treatment issues, especially in pregnant patients, the effect of synthetic human lactoferricin-derived peptides was tested successfully. These new BRAF-mutated MUG Mel3 cell lines represent a valuable model in melanoma heterogeneity and melanoma pregnancy research. Furthermore, treatment with anti-tumor peptides offers an alternative to conventionally used therapeutic options—especially during pregnancy.

Evolution of delayed resistance to immunotherapy in a melanoma responder

Despite initial responses1-3, most melanoma patients develop resistance4 to immune checkpoint blockade (ICB). To understand the evolution of resistance, we studied 37 tumor samples over 9 years from a patient with metastatic melanoma with complete clinical response to ICB followed by delayed recurrence and death. Phylogenetic analysis revealed co-evolution of seven lineages with multiple convergent, but independent resistance-associated alterations. All recurrent tumors emerged from a lineage characterized by loss of chromosome 15q, with post-treatment clones acquiring additional genomic driver events. Deconvolution of bulk RNA sequencing and highly multiplexed immunofluorescence (t-CyCIF) revealed differences in immune composition among different lineages. Imaging revealed a vasculogenic mimicry phenotype in NGFRhi tumor cells with high PD-L1 expression in close proximity to immune cells. Rapid autopsy demonstrated two distinct NGFR spatial patterns with high polarity and proximity to immune cells in subcutaneous tumors versus a diffuse spatial pattern in lung tumors, suggesting different roles of this neural-crest-like program in different tumor microenvironments. Broadly, this study establishes a high-resolution map of the evolutionary dynamics of resistance to ICB, characterizes a de-differentiated neural-crest tumor population in melanoma immunotherapy resistance and describes site-specific differences in tumor-immune interactions via longitudinal analysis of a patient with melanoma with an unusual clinical course.

Cutaneous Melanoma Classification: The Importance of High-Throughput Genomic Technologies

Cutaneous melanoma is an aggressive tumor responsible for 90% of mortality related to skin cancer. In the recent years, the discovery of driving mutations in melanoma has led to better treatment approaches. The last decade has seen a genomic revolution in the field of cancer. Such genomic revolution has led to the production of an unprecedented mole of data. High-throughput genomic technologies have facilitated the genomic, transcriptomic and epigenomic profiling of several cancers, including melanoma. Nevertheless, there are a number of newer genomic technologies that have not yet been employed in large studies. In this article we describe the current classification of cutaneous melanoma, we review the current knowledge of the main genetic alterations of cutaneous melanoma and their related impact on targeted therapies, and we describe the most recent high-throughput genomic technologies, highlighting their advantages and disadvantages. We hope that the current review will also help scientists to identify the most suitable technology to address melanoma-related relevant questions. The translation of this knowledge and all actual advancements into the clinical practice will be helpful in better defining the different molecular subsets of melanoma patients and provide new tools to address relevant questions on disease management. Genomic technologies might indeed allow to better predict the biological - and, subsequently, clinical - behavior for each subset of melanoma patients as well as to even identify all molecular changes in tumor cell populations during disease evolution toward a real achievement of a personalized medicine.

Intertumoral Genetic Heterogeneity Generates Distinct Tumor Microenvironments in a Novel Murine Synchronous Melanoma Model

Simple Summary

Metastatic melanoma patients may present with multiple, simultaneous metastases that are genetically different. This intertumoral heterogeneity can cause these tumors to respond differently to the same systemic therapy. Progression of any one tumor, even when others regress, eventually leads to therapy termination. The mechanism underlying these mixed responses remains unknown due to a lack of clinically representative animal models. In a novel murine model of synchronous melanoma that recapitulates human intertumoral heterogeneity, we show that intertumoral genetic heterogeneity leads to the simultaneous generation of distinct tumor immune microenvironments within the same mouse. Furthermore, each tumor can independently regulate local PD-1 (programmed cell death protein 1) and PD-L1 (PD-1 ligand) expressions, an immunosuppressive axis targeted by popular checkpoint immunotherapies. This model is useful for furthering the study of intertumoral heterogeneity and of lesion-specific therapeutic responses.

Abstract

Metastatic melanoma portends a poor prognosis and patients may present with multiple, simultaneous tumors. Despite recent advances in systemic immunotherapy, a majority of patients fail to respond, or exhibit lesion-specific responses wherein some metastases respond as others progress within the same patient. While intertumoral heterogeneity has been clinically associated with these mixed lesion-specific therapeutic responses, no clear mechanism has been identified, largely due to the scarcity of preclinical models. We developed a novel murine synchronous melanoma model that recapitulates this intertumoral genetic and microenvironmental heterogeneity. We show that genetic differences between tumors are sufficient to generate distinct tumor immune microenvironments (TIME) simultaneously in the same mouse. Furthermore, these TIMEs lead to the independent regulation of PD-1/PD-L1 (programmed cell death protein 1/PD-1 ligand), a popular axis targeted by immune checkpoint therapy, in response to ongoing anti-tumor immunity and the presence of interferon-gamma. Currently, therapeutic selection for metastatic melanoma patients is guided by a single biopsy, which may not represent the immune status of all tumors. As a result, patients can display heterogeneous lesion-specific responses. Further investigations into this synchronous melanoma model will provide mechanistic insight into the effects of intertumoral heterogeneity and guide therapeutic selection in this challenging patient population.

Identifying treatment options for BRAFV600 wild-type metastatic melanoma: A SU2C/MRA genomics-enabled clinical trial

Although combination BRAF and MEK inhibitors are highly effective for the 40-50% of cutaneous metastatic melanomas harboring BRAFV600 mutations, targeted agents have been ineffective for BRAFV600wild-type (wt) metastatic melanomas. The SU2C Genomics-Enabled Medicine for Melanoma Trial utilized a Simon two-stage optimal design to assess whether comprehensive genomic profiling improves selection of molecular-based therapies for BRAFV600wt metastatic melanoma patients who had progressed on standard-of-care therapy, which may include immunotherapy. Of the response-evaluable patients, binimetinib was selected for 20 patients randomized to the genomics-enabled arm, and nine were treated on the alternate treatment arm. Response rates for 27 patients treated with targeted recommendations included one (4%) partial response, 18 (67%) with stable disease, and eight (30%) with progressive disease. Post-trial genomic and protein pathway activation mapping identified additional drug classes that may be considered for future studies. Our results highlight the complexity and heterogeneity of metastatic melanomas, as well as how the lack of response in this trial may be associated with limitations including monotherapy drug selection and the dearth of available single and combination molecularly-driven therapies to treat BRAFV600wt metastatic melanomas.

The Genetic Evolution of Treatment-Resistant Cutaneous, Acral, and Uveal Melanomas

PURPOSE

Melanoma is a biologically heterogeneous disease composed of distinct clinicopathologic subtypes that frequently resist treatment. To explore the evolution of treatment resistance and metastasis, we used a combination of temporal and multilesional tumor sampling in conjunction with whole-exome sequencing of 110 tumors collected from 7 patients with cutaneous (n = 3), uveal (n = 2), and acral (n = 2) melanoma subtypes.

EXPERIMENTAL DESIGN

Primary tumors, metastases collected longitudinally, and autopsy tissues were interrogated. All but 1 patient died because of melanoma progression.

RESULTS

For each patient, we generated phylogenies and quantified the extent of genetic diversity among tumors, specifically among putative somatic alterations affecting therapeutic resistance.

CONCLUSIONS

In 4 patients who received immunotherapy, we found 1-3 putative acquired and intrinsic resistance mechanisms coexisting in the same patient, including mechanisms that were shared by all tumors within each patient, suggesting that future therapies directed at overcoming intrinsic resistance mechanisms may be broadly effective.

SATB2 induction of a neural crest mesenchyme-like program drives melanoma invasion and drug resistance

Recent genomic and scRNA-seq analyses of melanoma demonstrated a lack of recurrent genetic drivers of metastasis, while identifying common transcriptional states correlating with invasion or drug resistance. To test whether transcriptional adaptation can drive melanoma progression, we made use of a zebrafish mitfa:BRAFV600E;tp53-/- model, in which malignant progression is characterized by minimal genetic evolution. We undertook an overexpression-screen of 80 epigenetic/transcriptional regulators and found neural crest-mesenchyme developmental regulator SATB2 to accelerate aggressive melanoma development. Its overexpression induces invadopodia formation and invasion in zebrafish tumors and human melanoma cell lines. SATB2 binds and activates neural crest-regulators, including pdgfab and snai2. The transcriptional program induced by SATB2 overlaps with known MITF^lowAXL^high and AQP1⁺NGFR1^high drug-resistant states and functionally drives enhanced tumor propagation and resistance to Vemurafenib in vivo. In summary, we show that melanoma transcriptional rewiring by SATB2 to a neural crest mesenchyme-like program can drive invasion and drug resistance in autochthonous tumors.

A statistical framework for assessing pharmacological responses and biomarkers using uncertainty estimates

High-throughput testing of drugs across molecular-characterised cell lines can identify candidate treatments and discover biomarkers. However, the cells' response to a drug is typically quantified by a summary statistic from a best-fit dose-response curve, whilst neglecting the uncertainty of the curve fit and the potential variability in the raw readouts. Here, we model the experimental variance using Gaussian Processes, and subsequently, leverage uncertainty estimates to identify associated biomarkers with a new Bayesian framework. Applied to in vitro screening data on 265 compounds across 1074 cancer cell lines, our models identified 24 clinically established drug-response biomarkers, and provided evidence for six novel biomarkers by accounting for association with low uncertainty. We validated our uncertainty estimates with an additional drug screen of 26 drugs, 10 cell lines with 8 to 9 replicates. Our method is applicable to any dose-response data without replicates, and improves biomarker discovery for precision medicine.

Characterizations of Gene Alterations in Melanoma Patients from Chinese Population

Melanoma is a human skin malignant tumor with high invasion and poor prognosis. The limited understanding of genomic alterations in melanomas in China impedes the diagnosis and therapeutic strategy selection. We conducted comprehensive genomic profiling of melanomas from 39 primary and metastatic formalin-fixed paraffin-embedded (FFPE) samples from 27 patients in China based on an NGS panel of 223 genes. No significant difference in gene alterations was found between primary and metastasis melanomas. The status of germline mutation, CNV, and somatic mutation in our cohort was quite different from that reported in Western populations. We further delineated the mutation patterns of 4 molecular subgroups (BRAF, RAS, NF1, and Triple-WT) of melanoma in our cohort. BRAF mutations were more frequently identified in melanomas without chromic sun-induced damage (non-CSD), while RAS mutations were more likely observed in acral melanomas. NF1 and Triple-WT subgroups were unbiased between melanomas arising in non-CSD and acral skin. BRAF, RAS, and NF1 mutations were significantly associated with lymph node metastasis or presence of ulceration, implying that these cancer driver genes were independent prognostic factors. In summary, our results suggest that mutational profiles of malignant melanomas in China are significantly different from Western countries, and both gene mutation and amplification play an important role in the development and progression of melanomas.

Dynamic and unpredictable changes in mutant allele fractions of BRAF and NRAS during visceral progression of cutaneous malignant melanoma

Background

Data indicate that primary cutaneous melanomas are characterized by clonal heterogeneity associated with oncogenic drivers. Less data are available on the clonal changes occurring during melanoma progression. We therefore wished to analyse these changes in skin melanomas in common sites of visceral metastases as compared to the primary tumor.

Methods

An autopsy cohort of 50 patients with BRAF- and NRAS-mutant cutaneous metastatic melanomas including 139 visceral metastases was analysed for mutant allele fractions (MAF), determined by pyrosequencing and corrected for tumor/normal ratio. MAF levels were also classified as high (> 40%), medium (15–40%) or low (< 15%).

Results

Contrary to NRAS mutant cases, in BRAF-mutant melanomas MAFs were found to be significantly increased in visceral metastases compared to the primary due to the significantly higher levels in lung-, adrenal gland-, intestinal- and kidney metastases. The incidence of the three MAF variants in BRAF-mutant primaries was similar, whereas the high MAF cases were found to be increased in metastases. On the other hand, medium MAF levels were more common in case of NRAS-mutant tumors. Only 31.3% of BRAF mutant- and 50% of NRAS mutant cases maintained the MAF profile of the primary in metastasis. In the majority of multiple metastatic tumors, (BRAF:71.8%, NRAS:75%) metastases were relatively homogeneous regarding MAF. However, in 6/32(18.7%) of BRAF mutant cases low MAF primaries switched to high MAF in metastases. In heterogeneous BRAF mutant metastatic cases low to high or high to low MAF conversions occurred in a further 4/32(12.5%) cases in individual metastases as compared to the primary tumors. At lower frequency, in NRAS mutant tumor such changes also observed (2/12,16.7%).

Conclusion

We provided evidence for the selection of BRAF-mutant melanoma cells during metastatic progression to the lung, intestine, adrenal gland and kidney. Our findings suggest that in visceral metastases of malignant melanoma BRAF- or NRAS-MAFs are rather heterogeneous and cannot be predicted from data of the primary tumor. These data may have clinical significance when using targeted therapies.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5990-9) contains supplementary material, which is available to authorized users.

PTPRB promotes metastasis of colorectal carcinoma via inducing epithelial-mesenchymal transition

Dysregulation of protein tyrosine phosphatase, receptor type B (PTPRB) correlates with the development of a variety of tumors. Here we show that PTPRB promotes metastasis of colorectal cancer (CRC) cells via inducing epithelial-mesenchymal transition (EMT). We find that PTPRB is expressed at significantly higher levels in CRC tissues compared to adjacent nontumor tissues and in CRC cell lines with high invasion. PTPRB knockdown decreased the number of invasive CRC cells in an in vitro wound healing model, and also reduced tumor metastasis in vivo. Conversely, PTPRB overexpression promoted CRC cell invasion in vitro and metastasis in vivo. PTPRB overexpression decreased vimentin expression and promoted E-cadherin expression, consistent with promotion of EMT, while PTPRB knockdown had the opposite effect. Hypoxic conditions induced EMT and promoted invasion in CRC cells, but these effects were eliminated by PTPRB knockdown. EMT blockade via TWIST1 knockdown inhibited the migration and invasiveness of CRC cells, and even increased PTPRB expression could not reverse this effect. Altogether, these data support the conclusion that PTPRB promotes invasion and metastasis of CRC cells via inducing EMT, and that PTPRB would be a novel therapeutic target for the treatment of CRC.

Regeneration Enhances Metastasis: A Novel Role for Neurovascular Signaling in Promoting Melanoma Brain Metastasis

Neural repair after stroke involves initiation of a cellular proliferative program in the form of angiogenesis, neurogenesis, and molecular growth signals in the surrounding tissue elements. This cellular environment constitutes a niche in which regeneration of new blood vessels and new neurons leads to partial tissue repair after stroke. Cancer metastasis has similar proliferative cellular events in the brain and other organs. Do cancer and CNS tissue repair share similar cellular processes? In this study, we identify a novel role of the regenerative neurovascular niche induced by stroke in promoting brain melanoma metastasis through enhancing cellular interactions with surrounding niche components. Repair-mediated neurovascular signaling induces metastatic cells to express genes crucial to metastasis. Mimicking stroke-like conditions in vitro displays an enhancement of metastatic migration potential and allows for the determination of cell-specific signals produced by the regenerative neurovascular niche. Comparative analysis of both in vitro and in vivo expression profiles reveals a major contribution of endothelial cells in mediating melanoma metastasis. These results point to a previously undiscovered role of the regenerative neurovascular niche in shaping the tumor microenvironment and brain metastatic landscape.

High-resolution clonal mapping of multi-organ metastasis in triple negative breast cancer

Most triple negative breast cancers (TNBCs) are aggressively metastatic with a high degree of intra-tumoral heterogeneity (ITH), but how ITH contributes to metastasis is unclear. Here, clonal dynamics during metastasis were studied in vivo using two patient-derived xenograft (PDX) models established from the treatment-naive primary breast tumors of TNBC patients diagnosed with synchronous metastasis. Genomic sequencing and high-complexity barcode-mediated clonal tracking reveal robust alterations in clonal architecture between primary tumors and corresponding metastases. Polyclonal seeding and maintenance of heterogeneous populations of low-abundance subclones is observed in each metastasis. However, lung, liver, and brain metastases are enriched for an identical population of high-abundance subclones, demonstrating that primary tumor clones harbor properties enabling them to seed and thrive in multiple organ sites. Further, clones that dominate multi-organ metastases share a genomic lineage. Thus, intrinsic properties of rare primary tumor subclones enable the seeding and colonization of metastases in secondary organs in these models.

It is unclear how intra-tumoral heterogeneity contributes to metastasis. Here the authors study the clonal dynamics of triple negative breast cancer metastasis using patient derived xenografts and demonstrate that primary tumor clones harbor properties that support seeding and colonization of multiple organs.

PTPN11 Plays Oncogenic Roles and Is a Therapeutic Target for BRAF Wild-Type Melanomas

Melanoma is one of the most highly mutated cancer types. To identify functional drivers of melanoma, we searched for cross-species conserved mutations utilizing a mouse melanoma model driven by loss of PTEN and CDKN2A, and identified mutations in Kras, Erbb3, and Ptpn11. PTPN11 encodes the SHP2 protein tyrosine phosphatase that activates the RAS/RAF/MAPK pathway. Although PTPN11 is an oncogene in leukemia, lung, and breast cancers, its roles in melanoma are not clear. In this study, we found that PTPN11 is frequently activated in human melanoma specimens and cell lines and is required for full RAS/RAF/MAPK signaling activation in BRAF wild-type (either NRAS mutant or wild-type) melanoma cells. PTPN11 played oncogenic roles in melanoma by driving anchorage-independent colony formation and tumor growth. In Pten- and Cdkn2a-null mice, tet-inducible and melanocyte-specific PTPN11^E76K expression significantly enhanced melanoma tumorigenesis. Melanoma cells derived from this mouse model showed doxycycline-dependent tumor growth in nude mice. Silencing PTPN11^E76K expression by doxycycline withdrawal caused regression of established tumors by induction of apoptosis and senescence, and suppression of proliferation. Moreover, the PTPN11 inhibitor (SHP099) also caused regression of NRAS^Q61K -mutant melanoma. Using a quantitative tyrosine phosphoproteomics approach, we identified GSK3α/β as one of the key substrates that were differentially tyrosine-phosphorylated in these experiments modulating PTPN11. This study demonstrates that PTPN11 plays oncogenic roles in melanoma and regulates RAS and GSK3β signaling pathways. IMPLICATIONS: This study identifies PTPN11 as an oncogenic driver and a novel and actionable therapeutic target for BRAF wild-type melanoma.

miR-665 promotes hepatocellular carcinoma cell migration, invasion, and proliferation by decreasing Hippo signaling through targeting PTPRB

Growing evidence suggests that aberrant microRNA (miRNA) expression contributes to hepatocellular carcinoma (HCC) development and progression. However, the potential role and mechanism of miR-665 in the progression of liver cancer remains largely unknown. Our current study showed that miR-665 expression was upregulated in HCC cells and tissues. High expression of miR-665 exhibited more severe tumor size, vascular invasion and Edmondson grading in HCC patients. Gain- or loss-of-function assays demonstrated that miR-665 promoted cell proliferation, migration, invasion, and the epithelial-mesenchymal transition (EMT) of HCC cells in vitro and in vivo. Tyrosine phosphatase receptor type B (PTPRB) was downregulated in HCC tissues, and was negatively correlated with miR-665 expression. Through western blotting and luciferase reporter assay, PTPRB was identified as a direct downstream target of miR-665. Restoration of PTPRB reverses the effects of miR-665 on HCC migration, invasion, and cell proliferation. A mechanistic study showed that PTPTRB mediated the functional role of miR-665 through regulation of the Hippo signaling pathway. In conclusion, our results suggested that miR-665 was a negative regulator of the PTPRB and could promote tumor proliferation and metastasis in HCC through decreasing Hippo signaling pathway activity, which can be a potential target for HCC treatment.

Molecular Pathways in Melanomagenesis: What We Learned from Next-Generation Sequencing Approaches

PURPOSE OF REVIEW

Conventional clinico-pathological features in melanoma patients should be integrated with new molecular diagnostic, predictive, and prognostic factors coming from the expanding genomic profiles. Cutaneous melanoma (CM), even differing in biological behavior according to sun-exposure levels on the skin areas where it arises, is molecularly heterogeneous. The next-generation sequencing (NGS) approaches are providing data on mutation landscapes in driver genes that may account for distinct pathogenetic mechanisms and pathways. The purpose was to group and classify all somatic driver mutations observed in the main NGS-based studies.

RECENT FINDINGS

Whole exome and whole genome sequencing approaches have provided data on spectrum and distribution of genetic and genomic alterations as well as allowed to discover new cancer genes underlying CM pathogenesis. After evaluating the mutational status in a cohort of 686 CM cases from the most representative NGS studies, three molecular CM subtypes were proposed: BRAFmut, RASmut, and non-BRAFmut/non-RASmut.

Patterns of genomic evolution in advanced melanoma

Genomic alterations occurring during melanoma progression and the resulting genomic heterogeneity between metastatic deposits remain incompletely understood. Analyzing 86 metastatic melanoma deposits from 53 patients with whole-exome sequencing (WES), we show a low branch to trunk mutation ratio and little intermetastatic heterogeneity, with driver mutations almost completely shared between lesions. Branch mutations consistent with UV damage indicate that metastases may arise from different subclones in the primary tumor. Selective gain of mutated BRAF alleles occurs as an early event, contrasting whole-genome duplication (WGD) occurring as a late truncal event in about 40% of cases. One patient revealed elevated mutational diversity, probably related to previous chemotherapy and DNA repair defects. In another patient having received radiotherapy toward a lymph node metastasis, we detected a radiotherapy-related mutational signature in two subsequent distant relapses, consistent with secondary metastatic seeding. Our findings add to the understanding of genomic evolution in metastatic melanomas.

As melanoma progresses, it evolves. Here, in advanced melanoma the authors study genomic evolution, highlighting trunk mutations dominated by the ultraviolet damage signature, common late truncal whole-genome duplication events, as well as selective copy number gain of mutant BRAF.

Genomic and Transcriptomic Analysis Reveals Incremental Disruption of Key Signaling Pathways during Melanoma Evolution

We elucidated genomic and transcriptomic changes that accompany the evolution of melanoma from pre-malignant lesions by sequencing DNA and RNA from primary melanomas and their adjacent precursors, as well as matched primary tumors and regional metastases. In total, we analyzed 230 histopathologically distinct areas of melanocytic neoplasia from 82 patients. Somatic alterations sequentially induced mitogen-activated protein kinase (MAPK) pathway activation, upregulation of telomerase, modulation of the chromatin landscape, G1/S checkpoint override, ramp-up of MAPK signaling, disruption of the p53 pathway, and activation of the PI3K pathway; no mutations were specifically associated with metastatic progression, as these pathways were perturbed during the evolution of primary melanomas. UV radiation-induced point mutations steadily increased until melanoma invasion, at which point copy-number alterations also became prevalent.

Genetics of metastasis: melanoma and other cancers

Melanoma is a malignant neoplasm of melanocytes that accounts for the majority of skin cancer deaths despite comprising less than 5% of all cutaneous malignancies. Its incidence has increased faster than that of any other cancer over the past half-century and the annual costs of treatment in the United States alone have risen rapidly. Although the majority of primary melanomas are cured with local excision, metastatic melanoma historically carries a grim prognosis, with a median survival of 9 months and a long-term survival rate of 10%. Given the urgent need to develop treatment strategies for metastatic melanoma and the explosion of genetic technologies over the past 20 years, there has been extensive research into the genetic alterations that cause melanocytes to become malignant. More recently, efforts have focused on the genetic changes that drive melanoma metastasis. This review aims to summarize the current knowledge of the genetics of primary cutaneous and ocular melanoma, the genetic changes associated with metastasis in melanoma and other cancer types, and non-genetic factors that may contribute to metastasis.

Mutation profiles in early-stage lung squamous cell carcinoma with clinical follow-up and correlation with markers of immune function

Background

Lung squamous cell carcinoma (LUSC) accounts for 20–30% of non-small cell lung cancers (NSCLCs). There are limited treatment strategies for LUSC in part due to our inadequate understanding of the molecular underpinnings of the disease. We performed whole-exome sequencing (WES) and comprehensive immune profiling of a unique set of clinically annotated early-stage LUSCs to increase our understanding of the pathobiology of this malignancy.

Methods

Matched pairs of surgically resected stage I-III LUSCs and normal lung tissues (n = 108) were analyzed by WES. Immunohistochemistry and image analysis-based profiling of 10 immune markers were done on a subset of LUSCs (n = 91). Associations among mutations, immune markers and clinicopathological variables were statistically examined using analysis of variance and Fisher’s exact test. Cox proportional hazards regression models were used for statistical analysis of clinical outcome.

Results

This early-stage LUSC cohort displayed an average of 209 exonic mutations per tumor. Fourteen genes exhibited significant enrichment for somatic mutation: TP53, MLL2, PIK3CA, NFE2L2, CDH8, KEAP1, PTEN, ADCY8, PTPRT, CALCR, GRM8, FBXW7, RB1 and CDKN2A. Among mutated genes associated with poor recurrence-free survival, MLL2 mutations predicted poor prognosis in both TP53 mutant and wild-type LUSCs. We also found that in treated patients, FBXW7 and KEAP1 mutations were associated with poor response to adjuvant therapy, particularly in TP53-mutant tumors. Analysis of mutations with immune markers revealed that ADCY8 and PIK3CA mutations were associated with markedly decreased tumoral PD-L1 expression, LUSCs with PIK3CA mutations exhibited elevated CD45ro levels and CDKN2A-mutant tumors displayed an up-regulated immune response.

Conclusion(s)

Our findings pinpoint mutated genes that may impact clinical outcome as well as personalized strategies for targeted immunotherapies in early-stage LUSC.

BRAF Gene Copy Number and Mutant Allele Frequency Correlate with Time to Progression in Metastatic Melanoma Patients Treated with MAPK Inhibitors

Metastatic melanoma is characterized by complex genomic alterations, including a high rate of mutations in driver genes and widespread deletions and amplifications encompassing various chromosome regions. Among them, chromosome 7 is frequently gained in BRAF-mutant melanoma, inducing a mutant allele-specific imbalance. Although BRAF amplification is a known mechanism of acquired resistance to therapy with MAPK inhibitors, it is still unclear if BRAF copy-number variation and BRAF mutant allele imbalance at baseline can be associated with response to treatment. In this study, we used a multimodal approach to assess BRAF copy number and mutant allele frequency in pretreatment melanoma samples from 46 patients who received MAPK inhibitor-based therapy, and we analyzed the association with progression-free survival. We found that 65% patients displayed BRAF gains, often supported by chromosome 7 polysomy. In addition, we observed that 64% patients had a balanced BRAF-mutant/wild-type allele ratio, whereas 14% and 23% patients had low and high BRAF mutant allele frequency, respectively. Notably, a significantly higher risk of progression was observed in patients with a diploid BRAF status versus those with BRAF gains [HR, 2.86; 95% confidence interval (CI), 1.29-6.35; P = 0.01] and in patients with low percentage versus those with a balanced BRAF mutant allele percentage (HR, 4.54; 95% CI, 1.33-15.53; P = 0.016). Our data suggest that quantitative analysis of the BRAF gene could be useful to select the melanoma patients who are most likely to benefit from therapy with MAPK inhibitors. Mol Cancer Ther; 17(6); 1332-40. ©2018 AACR.

Inactivation of RASA1 promotes melanoma tumorigenesis via R-Ras activation

Inactivation of Ras GTPase activating proteins (RasGAPs) can activate Ras, increasing the risk for tumor development. Utilizing a melanoma whole genome sequencing (WGS) data from 13 patients, we identified two novel, clustered somatic missense mutations (Y472H and L481F) in RASA1 (RAS p21 protein activator 1, also called p120RasGAP). We have shown that wild type RASA1, but not identified mutants, suppresses soft agar colony formation and tumor growth of BRAF mutated melanoma cell lines via its RasGAP activity toward R-Ras (related RAS viral (r-ras) oncogene homolog) isoform. Moreover, R-Ras increased and RASA1 suppressed Ral-A activation among Ras downstream effectors. In addition to mutations, loss of RASA1 expression was frequently observed in metastatic melanoma samples on melanoma tissue microarray (TMA) and a low level of RASA1 mRNA expression was associated with decreased overall survival in melanoma patients with BRAF mutations. Thus, these data support that RASA1 is inactivated by mutation or by suppressed expression in melanoma and that RASA1 plays a tumor suppressive role by inhibiting R-Ras, a previously less appreciated member of the Ras small GTPases.

Intratumor and Intertumor Heterogeneity in Melanoma

Melanoma is a cancer that exhibits one of the most aggressive and heterogeneous features. The incidence rate escalates. A high number of clones harboring various mutations contribute to an exceptional level of intratumor heterogeneity of melanoma. It also refers to metastases which may originate from different subclones of primary lesion. Such component of the neoplasm biology is termed intertumor and intratumor heterogeneity. These levels of tumor heterogeneity hinder accurate diagnosis and effective treatment. The increasing number of research on the topic reflects the need for understanding limitation or failure of contemporary therapies. Majority of analyses concentrate on mutations in cancer-related genes. Novel high-throughput techniques reveal even higher degree of variations within a lesion. Consolidation of theories and researches indicates new routes for treatment options such as targets for immunotherapy. The demand for personalized approach in melanoma treatment requires extensive knowledge on intratumor and intertumor heterogeneity on the level of genome, transcriptome/proteome, and epigenome. Thus, achievements in exploration of melanoma variety are described in details. Particularly, the issue of tumor heterogeneity or homogeneity given BRAF mutations is discussed.

The receptor protein tyrosine phosphatase PTPRB negatively regulates FGF2-dependent branching morphogenesis

PTPRB is a transmembrane protein tyrosine phosphatase known to regulate blood vessel remodelling and angiogenesis. Here, we demonstrate that PTPRB negatively regulates branching morphogenesis in the mouse mammary epithelium. We show that Ptprb is highly expressed in adult mammary stem cells and also, although at lower levels, in oestrogen receptor-positive luminal cells. During mammary development, Ptprb expression is downregulated during puberty, a period of extensive ductal outgrowth and branching. In vivo shRNA knockdown of Ptprb in the cleared mammary fat pad transplant assay resulted in smaller epithelial outgrowths with an increased branching density and also increased branching in an in vitro organoid assay. Organoid branching was dependent on stimulation by FGF2, and Ptprb knockdown in mammary epithelial cells resulted in a higher level of fibroblast growth factor receptor (FGFR) activation and ERK1/2 phosphorylation, both at baseline and following FGF2 stimulation. Therefore, PTPRB regulates branching morphogenesis in the mammary epithelium by modulating the response of the FGFR signalling pathway to FGF stimulation. Considering the importance of branching morphogenesis in multiple taxa, our findings have general importance outside mammary developmental biology.

Sequencing of New and Old Therapies for Metastatic Melanoma

OPINION STATEMENT

Identification of BRAF driver mutations and agents that block their activity combined with development of immune checkpoint inhibitor therapies have dramatically changed survival and quality of life for patients with metastatic melanoma. Approximately half of patients with metastatic melanoma do not harbor mutations in the BRAF gene and therefore cannot benefit from currently available agents that target this mutation. Additionally, few patients with metastatic melanoma achieve durable disease control with these targeted therapies alone. Conversely, immune-based therapies have the potential to treat melanomas with or without mutations and produce durable responses following discontinuation of therapy, but responses can be delayed. Defining the goals of therapy (rapid response vs durable disease control), establishing the presence of targetable mutations, and considering the toxicities associated with each therapy can inform a treatment strategy. Incorporating both recent therapeutic modalities and older treatment options can provide the greatest potential for durable response. Overall, we recommend using immunotherapies (anti-CTLA4, anti-PD-1, combined anti-CTLA4/anti-PD-1, or interleukin-2) as the backbone of treatment for metastatic melanoma due to their potential for durable response. The targeted therapies and cytotoxic therapies can then be used intermittently to rescue patients from symptomatic disease progression. Of course, available clinical trials should always be considered, whenever possible.

Discrepancy in BRAF status among patients with metastatic malignant melanoma: A meta-analysis

The incidence of malignant melanoma is growing rapidly. Approximately half of the cases are BRAF mutated, making treatment with kinase inhibitors a (MEK and BRAF inhibitors) preferred choice in the advanced setting. The vast majority of these patients will benefit from the treatment. It is therefore of vital importance that the BRAF analysis is reliable and reflects the true nature of the tumour. Intraindividual tumour BRAF heterogeneity may exist, and changes of BRAF status over time might occur. We reviewed the literature by searching the PubMed database and 630 potentially relevant studies were identified. Thereafter, studies that investigated intralesional heterogeneity only, studies with ≤10 patients and studies that did not include adequate data to calculate discrepancy rates were excluded. Twenty-two studies met our inclusion criteria and were included in the meta-analysis. The pooled discrepancy rate between primary and metastatic lesions was 13.4% (95% confidence interval [CI]: 9.2-18.2%) while it was 7.3% (95% CI: 3.3-12.6) between two metastatic lesions. The number of patients whose tumoural BRAF status was changed from mutation to wild type and from wild type to mutation, respectively, was comparable. We conclude that a clinically meaningful discrepancy rate in BRAF status both between primary-metastatic and metastatic-metastatic melanoma lesions exists. Our results support the polyclonal model of melanomas in which subclones with different BRAF status co-exist in the same melanoma lesion. In addition, the results indicate a need for biopsy of a metastatic lesion for subsequent BRAF analysis when treatment with kinase inhibitors is considered.

From fish bowl to bedside: The power of zebrafish to unravel melanoma pathogenesis and discover new therapeutics

Melanoma is the most aggressive and deadliest form of skin cancer. A detailed knowledge of the cellular, molecular, and genetic events underlying melanoma progression is highly relevant to diagnosis, prognosis and risk stratification, and the development of new therapies. In the last decade, zebrafish have emerged as a valuable model system for the study of melanoma. Pathway conservation, coupled with the availability of robust genetic, transgenic, and chemical tools, has made the zebrafish a powerful model for identifying novel disease genes, visualizing cancer initiation, interrogating tumor-microenvironment interactions, and discovering new therapeutics that regulate melanocyte and melanoma development. In this review, we will give an overview of these studies, and highlight recent advancements that will help unravel melanoma pathogenesis and impact human disease.

DNA Mutations May Not Be the Cause of Cancer

Cancer is the most challenging disease of our time with increasing numbers of new cases each year, worldwide. Great achievements have been reached in cancer research through deep sequencing which helped define druggable targets. However, the still-evolving targeted therapy suffers resistance suggesting that DNA mutations considered as drivers may not have a role in tumor initiation. The present work discusses the role of DNA mutations as drivers and passengers in cancer initiation and development. First, it is important to discern the role of these DNA mutations as initiating events causing cancer or as contributors crucial for the development of a tumor once it has initiated. Second, breast cancer shown here illustrates how identification of DNA mutations in cancerous cells has influenced our approach for anti-cancer drug design. The cancer trilogy we have reached and described as: initial drug; resistance/recurrence; drug/treatment combinations, calls for a paradigm shift. To design more effective cancer drugs with durable and positive outcome, future cancer research needs to move beyond the sequencing era and explore changes which are taking place in cancer cells at levels other than the DNA. Evolutionary constraints may be acting as a barrier to preserve the human species from being transformed and, for that matter, all multi-cellular species which can incur cancer. Furthermore, mutations in the DNA do occur and for a multitude of reasons but without necessarily causing cancer. New directions will draw themselves when more focus is given to the event responsible for the switch of a cell from normalcy to malignancy. Until then, targeted therapy will certainly continue to improve the outcome of patients; however, it is unlikely to eradicate breast cancer depicted here.

Whole-genome landscapes of major melanoma subtypes

Melanoma of the skin is a common cancer only in Europeans, whereas it arises in internal body surfaces (mucosal sites) and on the hands and feet (acral sites) in people throughout the world. Here we report analysis of whole-genome sequences from cutaneous, acral and mucosal subtypes of melanoma. The heavily mutated landscape of coding and non-coding mutations in cutaneous melanoma resolved novel signatures of mutagenesis attributable to ultraviolet radiation. However, acral and mucosal melanomas were dominated by structural changes and mutation signatures of unknown aetiology, not previously identified in melanoma. The number of genes affected by recurrent mutations disrupting non-coding sequences was similar to that affected by recurrent mutations to coding sequences. Significantly mutated genes included BRAF, CDKN2A, NRAS and TP53 in cutaneous melanoma, BRAF, NRAS and NF1 in acral melanoma and SF3B1 in mucosal melanoma. Mutations affecting the TERT promoter were the most frequent of all; however, neither they nor ATRX mutations, which correlate with alternative telomere lengthening, were associated with greater telomere length. Most melanomas had potentially actionable mutations, most in components of the mitogen-activated protein kinase and phosphoinositol kinase pathways. The whole-genome mutation landscape of melanoma reveals diverse carcinogenic processes across its subtypes, some unrelated to sun exposure, and extends potential involvement of the non-coding genome in its pathogenesis.

Genomic and immune heterogeneity are associated with differential responses to therapy in melanoma

Appreciation for genomic and immune heterogeneity in cancer has grown though the relationship of these factors to treatment response has not been thoroughly elucidated. To better understand this, we studied a large cohort of melanoma patients treated with targeted therapy or immune checkpoint blockade (n = 60). Heterogeneity in therapeutic responses via radiologic assessment was observed in the majority of patients. Synchronous melanoma metastases were analyzed via deep genomic and immune profiling, and revealed substantial genomic and immune heterogeneity in all patients studied, with considerable diversity in T cell frequency, and few shared T cell clones (<8% on average) across the cohort. Variables related to treatment response were identified via these approaches and through novel radiomic assessment. These data yield insight into differential therapeutic responses to targeted therapy and immune checkpoint blockade in melanoma, and have key translational implications in the age of precision medicine.

Patients with metastatic melanoma display molecular and immune differences across tumor sites associated with differential drug responses. A team led by Jennifer Wargo from the University of Texas MD Anderson Cancer Center, Houston, USA, studied the radiological responses of 60 patients with metastatic melanoma, half of whom received targeted drug therapy and half of whom received an immune checkpoint inhibitor. The majority (83%) showed differences in responses across metastases. The group then profiled tumors in a subset, and found molecular and immune heterogeneity in different tumors within the same patient. Heterogeneity in mutational and immune profiles within tumors from individual patients could explain differences in treatment response. Knowing this, the authors emphasize the importance of acquiring biopsies from more than one tumor site in order to best tailor therapies to the features of metastatic cancer.

Recurrent epigenetic silencing of the PTPRD tumor suppressor in laryngeal squamous cell carcinoma

Cellular processes like differentiation, mitotic cycle, and cell growth are regulated by tyrosine kinases with known oncogenic potential and tyrosine phosphatases that downmodulate the first. Therefore, tyrosine phosphatases are recurrent targets of gene alterations in human carcinomas. We and others suggested recently a tumor suppressor function of the PTPRD tyrosine phosphatase and reported homozygous deletions of the PTPRD locus in laryngeal squamous cell carcinoma. In this study, we investigated other gene-inactivating mechanisms potentially targeting PTPRD, including loss-of-function mutations and also epigenetic alterations like promoter DNA hypermethylation. We sequenced the PTPRD gene in eight laryngeal squamous cell carcinoma cell lines but did not identify any inactivating mutations. In contrast, by bisulfite pyrosequencing of the gene promoter region, we identified significantly higher levels of methylation (p = 0.001 and p = 0.0002, respectively) in 9/14 (64%) laryngeal squamous cell carcinoma cell lines and 37/79 (47%) of primary laryngeal squamous cell carcinoma tumors as compared to normal epithelium of the upper aerodigestive tract. There was also a strong correlation (p = 0.0001) between methylation and transcriptional silencing for the PTPRD gene observed in a cohort of 497 head and neck tumors from The Cancer Genome Atlas dataset suggesting that DNA methylation is the main mechanism of PTPRD silencing in these tumors. In summary, our data provide further evidence of the high incidence of PTPRD inactivation in laryngeal squamous cell carcinoma. We suggest that deletions and loss-of-function mutations are responsible for PTPRD loss only in a fraction of cases, whereas DNA methylation is the dominating mechanism of PTPRD inactivation.

Genetic progression of malignant melanoma

Malignant melanoma of the skin is the most aggressive human cancer given that a primary tumor a few millimeters in diameter frequently has full metastatic competence. In view of that, revealing the genetic background of this potential may also help to better understand tumor dissemination in general. Genomic analyses have established the molecular classification of melanoma based on the most frequent driver oncogenic mutations (BRAF, NRAS, KIT) and have also revealed a long list of rare events, including mutations and amplifications as well as genetic microheterogeneity. At the moment, it is unclear whether any of these rare events have role in the metastasis initiation process since the major drivers do not have such a role. During lymphatic and hematogenous dissemination, the clonal selection process is evidently reflected by differences in oncogenic drivers in the metastases versus the primary tumor. Clonal selection is also evident during lymphatic progression, though the genetic background of this immunoselection is less clear. Genomic analyses of metastases identified further genetic alterations, some of which may correspond to metastasis maintenance genes. The natural genetic progression of melanoma can be modified by targeted (BRAF or MEK inhibitor) or immunotherapies. Some of the rare events in primary tumors may result in primary resistance, while further new genetic lesions develop during the acquired resistance to both targeted and immunotherapies. Only a few genetic lesions of the primary tumor are constant during natural or therapy-modulated progression. EGFR4 and NMDAR2 mutations, MITF and MET amplifications and PTEN loss can be considered as metastasis drivers. Furthermore, BRAF and MITF amplifications as well as PTEN loss are also responsible for resistance to targeted therapies, whereas NRAS mutation is the only founder genetic lesion showing any association with sensitivity to immunotherapies. Unfortunately, there are hardly any data on the possible organ-specific metastatic drivers in melanoma. These observations suggest that clinical management of melanoma patients must rely on the genetic analysis of the metastatic lesions to be able to monitor progression-associated changes and to personalize therapies.

The Many Roles of BAF (mSWI/SNF) and PBAF Complexes in Cancer

During the last decade, a host of epigenetic mechanisms were found to contribute to cancer and other human diseases. Several genomic studies have revealed that ∼20% of malignancies have alterations of the subunits of polymorphic BRG-/BRM-associated factor (BAF) and Polybromo-associated BAF (PBAF) complexes, making them among the most frequently mutated complexes in cancer. Recurrent mutations arise in genes encoding several BAF/PBAF subunits, including ARID1A, ARID2, PBRM1, SMARCA4, and SMARCB1 These subunits share some degree of conservation with subunits from related adenosine triphosphate (ATP)-dependent chromatin remodeling complexes in model organisms, in which a large body of work provides insight into their roles in cancer. Here, we review the roles of BAF- and PBAF-like complexes in these organisms, and relate these findings to recent discoveries in cancer epigenomics. We review several roles of BAF and PBAF complexes in cancer, including transcriptional regulation, DNA repair, and regulation of chromatin architecture and topology. More recent results highlight the need for new techniques to study these complexes.

From melanocytes to melanomas

Melanomas on sun-exposed skin are heterogeneous tumours, which can be subtyped on the basis of their cumulative levels of exposure to ultraviolet (UV) radiation. A melanocytic neoplasm can also be staged by how far it has progressed, ranging from a benign neoplasm, such as a naevus, to a malignant neoplasm, such as a metastatic melanoma. Each subtype of melanoma can evolve through distinct evolutionary trajectories, passing through (or sometimes skipping over) various stages of transformation. This Review delineates several of the more common progression trajectories that occur in the patient setting and proposes models for tumour evolution that integrate genetic, histopathological, clinical and biological insights from the melanoma literature.

Single cell genotyping of exome sequencing-identified mutations to characterize the clonal composition and evolution of inv(16) AML in a CBL mutated clonal hematopoiesis

We recently described the development of an inv(16) acute myeloid leukemia (AML) in a CBL mutated clonal hematopoiesis. Here, we further characterized the clonal composition and evolution of the AML based on the genetic information from the bulk specimen and analyses of individual bone marrow cells for mutations in CAND1, PTPRT, and DOCK6. To control for allele dropout, heterozygous polymorphisms located close to the respective mutation loci were assessed in parallel. The clonal composition concluded from exome sequencing suggested a proliferation advantage associated with the acquisition of mutations in CAND1, PTPRT, and DOCK6. Out of 102 single cell sequencing reactions on these mutations and the respective polymorphisms, analyses yielded conclusive results for at least 2 mutation sites in 12 cells. The single cell genotyping not only confirmed the co-occurrence of the PTPRT, CAND1 and DOCK6 mutations in the same AML clone but also revealed a clonal hierarchy, as the PTPRT mutation was likely acquired after the CAND1 and DOCK6 mutations. This insight had not been possible based solely on the exome sequencing data and suggests that the mutation in PTPRT, which encodes a STAT3-inhibiting protein tyrosine phosphatase, contributed to the AML development at a later stage by enhancing proliferation.

Integrated clinical, whole-genome, and transcriptome analysis of multisampled lethal metastatic prostate cancer

We report the first combined analysis of whole-genome sequence, detailed clinical history, and transcriptome sequence of multiple prostate cancer metastases in a single patient (A21). Whole-genome and transcriptome sequence was obtained from nine anatomically separate metastases, and targeted DNA sequencing was performed in cancerous and noncancerous foci within the primary tumor specimen removed 5 yr before death. Transcriptome analysis revealed increased expression of androgen receptor (AR)-regulated genes in liver metastases that harbored an AR p.L702H mutation, suggesting a dominant effect by the mutation despite being present in only one of an estimated 16 copies per cell. The metastases harbored several alterations to the PI3K/AKT pathway, including a clonal truncal mutation in PIK3CG and present in all metastatic sites studied. The list of truncal genomic alterations shared by all metastases included homozygous deletion of TP53, hemizygous deletion of RB1 and CHD1, and amplification of FGFR1. If the patient were treated today, given this knowledge, the use of second-generation androgen-directed therapies, cessation of glucocorticoid administration, and therapeutic inhibition of the PI3K/AKT pathway or FGFR1 receptor could provide personalized benefit. Three previously unreported truncal clonal missense mutations (ABCC4 p.R891L, ALDH9A1 p.W89R, and ASNA1 p.P75R) were expressed at the RNA level and assessed as druggable. The truncal status of mutations may be critical for effective actionability and merit further study. Our findings suggest that a large set of deeply analyzed cases could serve as a powerful guide to more effective prostate cancer basic science and personalized cancer medicine clinical trials.

The genomic landscape of cutaneous melanoma

Somatic mutation analysis of melanoma has been performed at the single gene level extensively over the past several decades. This has provided considerable insight into the critical pathways controlling melanoma initiation and progression. During the last 5 yr, next-generation sequencing (NGS) has enabled even more comprehensive mutational screening at the level of multigene panels, exomes and genomes. These studies have uncovered many new and unexpected players in melanoma development. The recent landmark study from The Cancer Genome Atlas (TCGA) consortium describing the genomic architecture of 333 cutaneous melanomas provides the largest and broadest analysis to date on the somatic aberrations underlying melanoma genesis. It thus seems timely to review the mutational landscape of melanoma and highlight the key genes and cellular pathways that appear to drive this cancer.

densityCut: an efficient and versatile topological approach for automatic clustering of biological data

Motivation: Many biological data processing problems can be formalized as clustering problems to partition data points into sensible and biologically interpretable groups.

Results: This article introduces densityCut, a novel density-based clustering algorithm, which is both time- and space-efficient and proceeds as follows: densityCut first roughly estimates the densities of data points from a K-nearest neighbour graph and then refines the densities via a random walk. A cluster consists of points falling into the basin of attraction of an estimated mode of the underlining density function. A post-processing step merges clusters and generates a hierarchical cluster tree. The number of clusters is selected from the most stable clustering in the hierarchical cluster tree. Experimental results on ten synthetic benchmark datasets and two microarray gene expression datasets demonstrate that densityCut performs better than state-of-the-art algorithms for clustering biological datasets. For applications, we focus on the recent cancer mutation clustering and single cell data analyses, namely to cluster variant allele frequencies of somatic mutations to reveal clonal architectures of individual tumours, to cluster single-cell gene expression data to uncover cell population compositions, and to cluster single-cell mass cytometry data to detect communities of cells of the same functional states or types. densityCut performs better than competing algorithms and is scalable to large datasets.

Availability and Implementation: Data and the densityCut R package is available from https://bitbucket.org/jerry00/densitycut_dev.

Contact: condon@cs.ubc.ca or sshah@bccrc.ca or jiaruid@cs.ubc.ca

Supplementary information:Supplementary data are available at Bioinformatics online.