Transposable Element Expression Profiles in Premalignant Pigment Cell Lesions and Melanoma of Xiphophorus.

Transposable elements (TEs) are characterized by their ability to change their genomic position. Through insertion or recombination leading to deletions and other chromosomal aberrations, they can cause genetic instability. The extent to which they thereby exert regulatory influence on cellular functions is unclear. To better characterize TEs in processes such as carcinogenesis, we used the well-established Xiphophorus melanoma model. By transcriptome sequencing, we show that an increasing total number in transposons correlates with progression of malignancy in melanoma samples from Xiphophorus interspecific hybrids. Further, by comparing the presence of TEs in the parental genomes of Xiphophorus maculatus and Xiphophorus hellerii, we could show that even in closely related species, genomic location and spectrum of TEs are considerably different.

Thiol starvation triggers melanoma state switching in an ATF4 and NRF2-dependent manner.

The cystine/glutamate antiporter xCT is an important source of cysteine for cancer cells. Once taken up, cystine is reduced to cysteine and serves as a building block for the synthesis of glutathione, which efficiently protects cells from oxidative damage and prevents ferroptosis. As melanomas are particularly exposed to several sources of oxidative stress, we investigated the biological role of cysteine and glutathione supply by xCT in melanoma. xCT activity was abolished by genetic depletion in the Tyr::CreER; BrafCA; Ptenlox/+ melanoma model and by acute cystine withdrawal in melanoma cell lines. Both interventions profoundly impacted melanoma glutathione levels, but they were surprisingly well tolerated by murine melanomas in vivo and by most human melanoma cell lines in vitro. RNA sequencing of human melanoma cells revealed a strong adaptive upregulation of NRF2 and ATF4 pathways, which orchestrated the compensatory upregulation of genes involved in antioxidant defence and de novo cysteine biosynthesis. In addition, the joint activation of ATF4 and NRF2 triggered a phenotypic switch characterized by a reduction of differentiation genes and induction of pro-invasive features, which was also observed after erastin treatment or the inhibition of glutathione synthesis. NRF2 alone was capable of inducing the phenotypic switch in a transient manner. Together, our data show that cystine or glutathione levels regulate the phenotypic plasticity of melanoma cells by elevating ATF4 and NRF2.

Critical Evaluation of a microRNA-Based Risk Classifier Predicting Cancer-Specific Survival in Renal Cell Carcinoma with Tumor Thrombus of the Inferior Vena Cava.

(1) Background: Clear cell renal cell carcinoma extending into the inferior vena cava (ccRCCIVC) represents a clinical high-risk setting. However, there is substantial heterogeneity within this patient subgroup regarding survival outcomes. Previously, members of our group developed a microRNA(miR)-based risk classifier-containing miR-21-5p, miR-126-3p and miR-221-3p expression-which significantly predicted the cancer-specific survival (CSS) of ccRCCIVC patients. (2) Methods: Examining a single-center cohort of tumor tissue from n = 56 patients with ccRCCIVC, we measured the expression levels of miR-21, miR-126, and miR-221 using qRT-PCR. The prognostic impact of clinicopathological parameters and miR expression were investigated via single-variable and multivariable Cox regression. Referring to the previously established risk classifier, we performed Kaplan-Meier analyses for single miR expression levels and the combined risk classifier. Cut-off values and weights within the risk classifier were taken from the previous study. (3) Results: miR-21 and miR-126 expression were significantly associated with lymphonodal status at the time of surgery, the development of metastasis during follow-up, and cancer-related death. In Kaplan-Meier analyses, miR-21 and miR-126 significantly impacted CSS in our cohort. Moreover, applying the miR-based risk classifier significantly stratified ccRCCIVC according to CSS. (4) Conclusions: In our retrospective analysis, we successfully validated the miR-based risk classifier within an independent ccRCCIVC cohort.

Oncogenic YAP mediates changes in chromatin accessibility and activity that drive cell cycle gene expression and cell migration.

YAP, the key protein effector of the Hippo pathway, is a transcriptional co-activator that controls the expression of cell cycle genes, promotes cell growth and proliferation and regulates organ size. YAP modulates gene transcription by binding to distal enhancers, but the mechanisms of gene regulation by YAP-bound enhancers remain poorly understood. Here we show that constitutive active YAP5SA leads to widespread changes in chromatin accessibility in untransformed MCF10A cells. Newly accessible regions include YAP-bound enhancers that mediate activation of cycle genes regulated by the Myb-MuvB (MMB) complex. By CRISPR-interference we identify a role for YAP-bound enhancers in phosphorylation of Pol II at Ser5 at MMB-regulated promoters, extending previously published studies that suggested YAP primarily regulates the pause-release step and transcriptional elongation. YAP5SA also leads to less accessible 'closed' chromatin regions, which are not directly YAP-bound but which contain binding motifs for the p53 family of transcription factors. Diminished accessibility at these regions is, at least in part, a consequence of reduced expression and chromatin-binding of the p53 family member ΔNp63 resulting in downregulation of ΔNp63-target genes and promoting YAP-mediated cell migration. In summary, our studies uncover changes in chromatin accessibility and activity that contribute to the oncogenic activities of YAP.

Salvage Nodal Radiotherapy as Metastasis-Directed Therapy for Oligorecurrent Prostate Cancer Detected by Positron Emission Tomography Shows Favorable Outcome in Long-Term Follow-Up.

BACKGROUND: The study aimed to access the long-term outcome of salvage nodal radiotherapy (SNRT) in oligorecurrent prostate cancer. METHODS: A total of 95 consecutive patients received SNRT for pelvic and/or extrapelvic nodal recurrence after prostate-specific membrane antigen (PSMA) or choline PET from 2010 to 2021. SNRT was applied as external beam radiotherapy with simultaneous integrated boost up to a median total dose of 62.9 Gy (EQD21.5Gy) to the recurrent lymph node metastases. The outcome was analyzed by cumulative incidence functions with death as the competing risk. Fine-Gray regression analyses were performed to estimate the relative hazards of the outcome parameters. Genitourinary (GU)/gastrointestinal (GI) toxicity evaluation utilized Common Toxicity Criteria for Adverse Events (v5.0). The results are as follows: the median follow-up was 47.1 months. The five-year biochemical progression rate (95% CI) was 50.1% (35.7-62.9%). Concomitant androgen deprivation therapy (ADT) was adminstered in 60.0% of the patients. The five-year biochemical progression rate was 75.0% (42.0-90.9%) without ADT versus 35.3% (19.6-51.4%) with ADT (p = 0.003). The cumulative five-year late grade 3 GU toxicity rate was 2.1%. No late grade 3 GI toxicity occured. CONCLUSIONS: Metastasis-directed therapy through SNRT for PET-staged oligorecurrent prostate cancer demonstrated a favorable long-term oncologic outcome. Omittance of ADT led to an increased biochemical progression.

The DGCR8 E518K mutation found in Wilms tumors leads to a partial miRNA processing defect that alters gene expression patterns and biological processes.

Wilms tumor (WT) is the most common renal tumor in childhood. We and others have previously identified oncogenic driver mutations affecting the microprocessor genes DROSHA and DGCR8 that lead to altered miRNA expression patterns. In the case of DGCR8, a single recurrent hotspot mutation (E518K) was found in the RNA binding domain. To functionally assess this mutation in vitro, we generated mouse Dgcr8-KO embryonic stem cell (mESC) lines with an inducible expression of wild-type or mutant DGCR8, mirroring the hemizygous mutant expression seen in WT. RNA-seq analysis revealed significant differences of miRNA expression profiles in DGCR8-E518K compared with DGCR8-wild-type mESCs. The E518K mutation only led to a partial rescue of the reported miRNA processing defect in Dgcr8-KO, with selectively reduced expression of numerous canonical miRNAs. Nevertheless, DGCR8-E518K retained significant activity given its ability to still process many miRNAs. Subsequent to altered miRNA levels, the expression of mRNA targets was likewise changed. Functional assays showed that DGCR8-E518K cells still have a partial proliferation and differentiation defect but were able to rescue critical biological processes in embryoid body development. The stem cell program could be shut down and all three germ layers were formed. These findings suggest that the E518K mutation leads to a partial reduction of microprocessor activity and altered specificity with selective impairment only in certain developmental contexts, apparently including nephrogenesis.

Differential expression of transposable elements in the medaka melanoma model.

Malignant melanoma incidence is rising worldwide. Its treatment in an advanced state is difficult, and the prognosis of this severe disease is still very poor. One major source of these difficulties is the high rate of metastasis and increased genomic instability leading to a high mutation rate and the development of resistance against therapeutic approaches. Here we investigate as one source of genomic instability the contribution of activation of transposable elements (TEs) within the tumor. We used the well-established medaka melanoma model and RNA-sequencing to investigate the differential expression of TEs in wildtype and transgenic fish carrying melanoma. We constructed a medaka-specific TE sequence library and identified TE sequences that were specifically upregulated in tumors. Validation by qRT- PCR confirmed a specific upregulation of a LINE and an LTR element in malignant melanomas of transgenic fish.

Oncogenic allelic interaction in Xiphophorus highlights hybrid incompatibility.

Mixing genomes of different species by hybridization can disrupt species-specific genetic interactions that were adapted and fixed within each species population. Such disruption can predispose the hybrids to abnormalities and disease that decrease the overall fitness of the hybrids and is therefore named as hybrid incompatibility. Interspecies hybridization between southern platyfish and green swordtails leads to lethal melanocyte tumorigenesis. This occurs in hybrids with tumor incidence following progeny ratio that is consistent with two-locus interaction, suggesting melanoma development is a result of negative epistasis. Such observations make Xiphophorus one of the only two vertebrate hybrid incompatibility examples in which interacting genes have been identified. One of the two interacting loci has been characterized as a mutant epidermal growth factor receptor. However, the other locus has not been identified despite over five decades of active research. Here we report the localization of the melanoma regulatory locus to a single gene, rab3d, which shows all expected features of the long-sought oncogene interacting locus. Our findings provide insights into the role of egfr regulation in regard to cancer etiology. Finally, they provide a molecular explainable example of hybrid incompatibility.

The transcription factor NRF2 enhances melanoma malignancy by blocking differentiation and inducing COX2 expression.

The transcription factor NRF2 is the major mediator of oxidative stress responses and is closely connected to therapy resistance in tumors harboring activating mutations in the NRF2 pathway. In melanoma, such mutations are rare, and it is unclear to what extent melanomas rely on NRF2. Here we show that NRF2 suppresses the activity of the melanocyte lineage marker MITF in melanoma, thereby reducing the expression of pigmentation markers. Intriguingly, we furthermore identified NRF2 as key regulator of immune-modulating genes, linking oxidative stress with the induction of cyclooxygenase 2 (COX2) in an ATF4-dependent manner. COX2 is critical for the secretion of prostaglandin E2 and was strongly induced by H2O2 or TNFα only in presence of NRF2. Induction of MITF and depletion of COX2 and PGE2 were also observed in NRF2-deleted melanoma cells in vivo. Furthermore, genes corresponding to the innate immune response such as RSAD2 and IFIH1 were strongly elevated in absence of NRF2 and coincided with immune evasion parameters in human melanoma datasets. Even in vitro, NRF2 activation or prostaglandin E2 supplementation blunted the induction of the innate immune response in melanoma cells. Transcriptome analyses from lung adenocarcinomas indicate that the observed link between NRF2 and the innate immune response is not restricted to melanoma.

Cxcl9l and Cxcr3.2 regulate recruitment of osteoclast progenitors to bone matrix in a medaka osteoporosis model.

Bone homeostasis requires continuous remodeling of bone matrix to maintain structural integrity. This involves extensive communication between bone-forming osteoblasts and bone-resorbing osteoclasts to orchestrate balanced progenitor cell recruitment and activation. Only a few mediators controlling progenitor activation are known to date and have been targeted for intervention of bone disorders such as osteoporosis. To identify druggable pathways, we generated a medaka (Oryzias latipes) osteoporosis model, where inducible expression of receptor-activator of nuclear factor kappa-Β ligand (Rankl) leads to ectopic formation of osteoclasts and excessive bone resorption, which can be assessed by live imaging. Here we show that upon Rankl induction, osteoblast progenitors up-regulate expression of the chemokine ligand Cxcl9l. Ectopic expression of Cxcl9l recruits mpeg1-positive macrophages to bone matrix and triggers their differentiation into osteoclasts. We also demonstrate that the chemokine receptor Cxcr3.2 is expressed in a distinct subset of macrophages in the aorta-gonad-mesonephros (AGM). Live imaging revealed that upon Rankl induction, Cxcr3.2-positive macrophages get activated, migrate to bone matrix, and differentiate into osteoclasts. Importantly, mutations in cxcr3.2 prevent macrophage recruitment and osteoclast differentiation. Furthermore, Cxcr3.2 inhibition by the chemical antagonists AMG487 and NBI-74330 also reduced osteoclast recruitment and protected bone integrity against osteoporotic insult. Our data identify a mechanism for progenitor recruitment to bone resorption sites and Cxcl9l and Cxcr3.2 as potential druggable regulators of bone homeostasis and osteoporosis.

The sterlet sturgeon genome sequence and the mechanisms of segmental rediploidization.

Sturgeons seem to be frozen in time. The archaic characteristics of this ancient fish lineage place it in a key phylogenetic position at the base of the ~30,000 modern teleost fish species. Moreover, sturgeons are notoriously polyploid, providing unique opportunities to investigate the evolution of polyploid genomes. We assembled a high-quality chromosome-level reference genome for the sterlet, Acipenser ruthenus. Our analysis revealed a very low protein evolution rate that is at least as slow as in other deep branches of the vertebrate tree, such as that of the coelacanth. We uncovered a whole-genome duplication that occurred in the Jurassic, early in the evolution of the entire sturgeon lineage. Following this polyploidization, the rediploidization of the genome included the loss of whole chromosomes in a segmental deduplication process. While known adaptive processes helped conserve a high degree of structural and functional tetraploidy over more than 180 million years, the reduction of redundancy of the polyploid genome seems to have been remarkably random.

miR-221-3p Regulates VEGFR2 Expression in High-Risk Prostate Cancer and Represents an Escape Mechanism from Sunitinib In Vitro.

Downregulation of miR-221-3p expression in prostate cancer (PCa) predicted overall and cancer-specific survival of high-risk PCa patients. Apart from PCa, miR-221-3p expression levels predicted a response to tyrosine kinase inhibitors (TKI) in clear cell renal cell carcinoma (ccRCC) patients. Since this role of miR-221-3p was explained with a specific targeting of VEGFR2, we examined whether miR-221-3p regulated VEGFR2 in PCa. First, we confirmed VEGFR2/KDR as a target gene of miR-221-3p in PCa cells by applying Luciferase reporter assays and Western blotting experiments. Although VEGFR2 was mainly downregulated in the PCa cohort of the TCGA (The Cancer Genome Atlas) database, VEGFR2 was upregulated in our high-risk PCa cohort (n = 142) and predicted clinical progression. In vitro miR-221-3p acted as an escape mechanism from TKI in PC3 cells, as displayed by proliferation and apoptosis assays. Moreover, we confirmed that Sunitinib induced an interferon-related gene signature in PC3 cells by analyzing external microarray data and by demonstrating a significant upregulation of miR-221-3p/miR-222-3p after Sunitinib exposure. Our findings bear a clinical perspective for high-risk PCa patients with low miR-221-3p levels since this could predict a favorable TKI response. Apart from this therapeutic niche, we identified a partially oncogenic function of miR-221-3p as an escape mechanism from VEGFR2 inhibition.

miR-221 Augments TRAIL-Mediated Apoptosis in Prostate Cancer Cells by Inducing Endogenous TRAIL Expression and Targeting the Functional Repressors SOCS3 and PIK3R1.

miR-221 is regarded as an oncogene in many malignancies, and miR-221-mediated resistance towards TRAIL was one of the first oncogenic roles shown for this small noncoding RNA. In contrast, miR-221 is downregulated in prostate cancer (PCa), thereby implying a tumour suppressive function. By using proliferation and apoptosis assays, we show a novel feature of miR-221 in PCa cells: instead of inducing TRAIL resistance, miR-221 sensitized cells towards TRAIL-induced proliferation inhibition and apoptosis induction. Partially responsible for this effect was the interferon-mediated gene signature, which among other things contained an endogenous overexpression of the TRAIL encoding gene TNFSF10. This TRAIL-friendly environment was provoked by downregulation of the established miR-221 target gene SOCS3. Moreover, we introduced PIK3R1 as a target gene of miR-221 in PCa cells. Proliferation assays showed that siRNA-mediated downregulation of SOCS3 and PIK3R1 mimicked the effect of miR-221 on TRAIL sensitivity. Finally, Western blotting experiments confirmed lower amounts of phospho-Akt after siRNA-mediated downregulation of PIK3R1 in PC3 cells. Our results further support the tumour suppressing role of miR-221 in PCa, since it sensitises PCa cells towards TRAIL by regulating the expression of the oncogenes SOCS3 and PIK3R1. Given the TRAIL-inhibiting effect of miR-221 in various cancer entities, our results suggest that the influence of miR-221 on TRAIL-mediated apoptosis is highly context- and entity-dependent.

The Myb-MuvB Complex Is Required for YAP-Dependent Transcription of Mitotic Genes.

YAP and TAZ, downstream effectors of the Hippo pathway, are important regulators of proliferation. Here, we show that the ability of YAP to activate mitotic gene expression is dependent on the Myb-MuvB (MMB) complex, a master regulator of genes expressed in the G2/M phase of the cell cycle. By carrying out genome-wide expression and binding analyses, we found that YAP promotes binding of the MMB subunit B-MYB to the promoters of mitotic target genes. YAP binds to B-MYB and stimulates B-MYB chromatin association through distal enhancer elements that interact with MMB-regulated promoters through chromatin looping. The cooperation between YAP and B-MYB is critical for YAP-mediated entry into mitosis. Furthermore, the expression of genes coactivated by YAP and B-MYB is associated with poor survival of cancer patients. Our findings provide a molecular mechanism by which YAP and MMB regulate mitotic gene expression and suggest a link between two cancer-relevant signaling pathways.

Expression Signatures of Cisplatin- and Trametinib-Treated Early-Stage Medaka Melanomas.

Small aquarium fish models provide useful systems not only for a better understanding of the molecular basis of many human diseases, but also for first-line screening to identify new drug candidates. For testing new chemical substances, current strategies mostly rely on easy to perform and efficient embryonic screens. Cancer, however, is a disease that develops mainly during juvenile and adult stage. Long-term treatment and the challenge to monitor changes in tumor phenotype make testing of large chemical libraries in juvenile and adult animals cost prohibitive. We hypothesized that changes in the gene expression profile should occur early during anti-tumor treatment, and the disease-associated transcriptional change should provide a reliable readout that can be utilized to evaluate drug-induced effects. For the current study, we used a previously established medaka melanoma model. As proof of principle, we showed that exposure of melanoma developing fish to the drugs cisplatin or trametinib, known cancer therapies, for a period of seven days is sufficient to detect treatment-induced changes in gene expression. By examining whole body transcriptome responses we provide a novel route toward gene panels that recapitulate anti-tumor outcomes thus allowing a screening of thousands of drugs using a whole-body vertebrate model. Our results suggest that using disease-associated transcriptional change to screen therapeutic molecules in small fish model is viable and may be applied to pre-clinical research and development stages in new drug discovery.

Application of the Transcriptional Disease Signature (TDSs) to Screen Melanoma-Effective Compounds in a Small Fish Model.

Cell culture and protein target-based compound screening strategies, though broadly utilized in selecting candidate compounds, often fail to eliminate candidate compounds with non-target effects and/or safety concerns until late in the drug developmental process. Phenotype screening using intact research animals is attractive because it can help identify small molecule candidate compounds that have a high probability of proceeding to clinical use. Most FDA approved, first-in-class small molecules were identified from phenotypic screening. However, phenotypic screening using rodent models is labor intensive, low-throughput, and very expensive. As a novel alternative for small molecule screening, we have been developing gene expression disease profiles, termed the Transcriptional Disease Signature (TDS), as readout of small molecule screens for therapeutic molecules. In this concept, compounds that can reverse, or otherwise affect known disease-associated gene expression patterns in whole animals may be rapidly identified for more detailed downstream direct testing of their efficacy and mode of action. To establish proof of concept for this screening strategy, we employed a transgenic strain of a small aquarium fish, medaka (Oryzias latipes), that overexpresses the malignant melanoma driver gene xmrk, a mutant egfr gene, that is driven by a pigment cell-specific mitf promoter. In this model, melanoma develops with 100% penetrance. Using the transgenic medaka malignant melanoma model, we established a screening system that employs the NanoString nCounter platform to quantify gene expression within custom sets of TDS gene targets that we had previously shown to exhibit differential transcription among xmrk-transgenic and wild-type medaka. Compound-modulated gene expression was identified using an internet-accessible custom-built data processing pipeline. The effect of a given drug on the entire TDS profile was estimated by comparing compound-modulated genes in the TDS using an activation Z-score and Kolmogorov-Smirnov statistics. TDS gene probes were designed that target common signaling pathways that include proliferation, development, toxicity, immune function, metabolism and detoxification. These pathways may be utilized to evaluate candidate compounds for potential favorable, or unfavorable, effects on melanoma-associated gene expression. Here we present the logistics of using medaka to screen compounds, as well as, the development of a user-friendly NanoString data analysis pipeline to support feasibility of this novel TDS drug-screening strategy.

Loss or oncogenic mutation of DROSHA impairs kidney development and function, but is not sufficient for Wilms tumor formation.

Wilms tumor (WT) is the most common kidney cancer in childhood. Mutations in the microprocessor genes DROSHA and DGCR8 have been identified as putative oncogenic drivers, indicating a critical role of aberrant miRNA processing in WT formation. To characterize the in vivo role of DROSHA mutations during kidney development and their oncogenic potential, we analyzed mouse lines with either a targeted deletion of Drosha or an inducible expression of human DROSHA carrying a tumor-specific E1147K mutation that acts in a dominant negative manner. Both types of mutation induce striking changes in miRNA patterns. Six2-cre mediated deletion of Drosha in nephron progenitors led to perinatal lethality with apoptotic loss of progenitor cells and early termination of nephrogenesis. Mosaic deletions via Wt1-creERT2 resulted in a milder phenotype with viable offspring that developed proteinuria after 2-4 weeks, but no evidence of tumor formation. Activation of the DROSHA-E1147K transgene via Six2-cre, on the other hand, induced a more severe phenotype with apoptosis of progenitor cells, proteinuria and glomerular sclerosis. The severely growth retarded mice died within the first 2 months of life, confirming the predicted dominant-negative effect of DROSHA-E1147K in vivo. While our data underscores the importance of a viable self-renewing progenitor pool for kidney development, there was no evidence of tumor formation through impaired DROSHA function. This suggests that either additional alterations in mitogenic or antiapoptotic pathways are needed for malignant transformation, or premature loss of a susceptible target cell population and early lethality prevent WT formation.

Analysis of the putative tumor suppressor gene cdkn2ab in pigment cells and melanoma of Xiphophorus and medaka.

In humans, the CDKN2A locus encodes two transcripts, INK4A and ARF. Inactivation of either one by mutations or epigenetic changes is a frequent signature of malignant melanoma and one of the most relevant entry points for melanomagenesis. To analyze whether cdkn2ab, the fish ortholog of CDKN2A, has a similar function as its human counterpart, we studied its action in fish models for human melanoma. Overexpression of cdkn2ab in a Xiphophorus melanoma cell line led to decreased proliferation and induction of a senescence-like phenotype, indicating a melanoma-suppressive function analogous to mammals. Coexpression of Xiphophorus cdkn2ab in medaka transgenic for the mitfa:xmrk melanoma-inducing gene resulted in full suppression of melanoma development, whereas CRISPR/Cas9 knockout of cdkn2ab resulted in strongly enhanced tumor growth. In summary, this provides the first functional evidence that cdkn2ab acts as a potent tumor suppressor gene in fish melanoma models.

RNA-seq analysis identifies different transcriptomic types and developmental trajectories of primary melanomas.

Recent studies revealed trajectories of mutational events in early melanomagenesis, but the accompanying changes in gene expression are far less understood. Therefore, we performed a comprehensive RNA-seq analysis of laser-microdissected melanocytic nevi (n = 23) and primary melanoma samples (n = 57) and characterized the molecular mechanisms of early melanoma development. Using self-organizing maps, unsupervised clustering, and analysis of pseudotime (PT) dynamics to identify evolutionary trajectories, we describe here two transcriptomic types of melanocytic nevi (N1 and N2) and primary melanomas (M1 and M2). N1/M1 lesions are characterized by pigmentation-type and MITF gene signatures, and a high prevalence of NRAS mutations in M1 melanomas. N2/M2 lesions are characterized by inflammatory-type and AXL gene signatures with an equal distribution of wild-type and mutated BRAF and low prevalence of NRAS mutations in M2 melanomas. Interestingly, N1 nevi and M1 melanomas and N2 nevi and M2 melanomas, respectively, cluster together, but there is no clustering in a stage-dependent manner. Transcriptional signatures of M1 melanomas harbor signatures of BRAF/MEK inhibitor resistance and M2 melanomas harbor signatures of anti-PD-1 antibody treatment resistance. Pseudotime dynamics of nevus and melanoma samples are suggestive for a switch-like immune-escape mechanism in melanoma development with downregulation of immune genes paralleled by an increasing expression of a cell cycle signature in late-stage melanomas. Taken together, the transcriptome analysis identifies gene signatures and mechanisms underlying development of melanoma in early and late stages with relevance for diagnostics and therapy.