Human melanocyte development and melanoma dedifferentiation at single-cell resolution

Melanocyte lineage dynamics in development, growth and disease

Melanocytes evolved to produce the melanin that gives colour to our hair, eyes and skin. The melanocyte lineage also gives rise to melanoma, the most lethal form of skin cancer. The melanocyte lineage differentiates from neural crest cells during development, and most melanocytes reside in the skin and hair, where they are replenished by melanocyte stem cells. Because the molecular mechanisms necessary for melanocyte specification, migration, proliferation and differentiation are co-opted during melanoma initiation and progression, studying melanocyte development is directly relevant to human disease. Here, through the lens of advances in cellular omic and genomic technologies, we review the latest findings in melanocyte development and differentiation, and how these developmental pathways become dysregulated in disease.

The genetic evolution of acral melanoma

Acral melanoma is an aggressive type of melanoma with unknown origins. It is the most common type of melanoma in individuals with dark skin and is notoriously challenging to treat. We examine exome sequencing data of 139 tissue samples, spanning different progression stages, from 37 patients. We find that 78.4% of the melanomas display clustered copy number transitions with focal amplifications, recurring predominantly on chromosomes 5, 11, 12, and 22. These complex genomic aberrations are typically shared across all progression stages of individual patients. TERT activating alterations also arise early, whereas MAP-kinase pathway mutations appear later, an inverted order compared to the canonical evolution. The punctuated formation of complex aberrations and early TERT activation suggest a unique mutational mechanism that initiates acral melanoma. The marked intratumoral heterogeneity, especially concerning MAP-kinase pathway mutations, may partly explain the limited success of therapies for this melanoma subtype.

Unravelling tumour cell diversity and prognostic signatures in cutaneous melanoma through machine learning analysis

Melanoma, a highly malignant tumour, presents significant challenges due to its cellular heterogeneity, yet research on this aspect in cutaneous melanoma remains limited. In this study, we utilized single-cell data from 92,521 cells to explore the tumour cell landscape. Through clustering analysis, we identified six distinct cell clusters and investigated their differentiation and metabolic heterogeneity using multi-omics approaches. Notably, cytotrace analysis and pseudotime trajectories revealed distinct stages of tumour cell differentiation, which have implications for patient survival. By leveraging markers from these clusters, we developed a tumour cell-specific machine learning model (TCM). This model not only predicts patient outcomes and responses to immunotherapy, but also distinguishes between genomically stable and unstable tumours and identifies inflamed ('hot') versus non-inflamed ('cold') tumours. Intriguingly, the TCM score showed a strong association with TOMM40, which we experimentally validated as an oncogene promoting tumour proliferation, invasion and migration. Overall, our findings introduce a novel biomarker score that aids in selecting melanoma patients for improved prognoses and targeted immunotherapy, thereby guiding clinical treatment decisions.

Tumor initiation and early tumorigenesis: molecular mechanisms and interventional targets

Tumorigenesis is a multistep process, with oncogenic mutations in a normal cell conferring clonal advantage as the initial event. However, despite pervasive somatic mutations and clonal expansion in normal tissues, their transformation into cancer remains a rare event, indicating the presence of additional driver events for progression to an irreversible, highly heterogeneous, and invasive lesion. Recently, researchers are emphasizing the mechanisms of environmental tumor risk factors and epigenetic alterations that are profoundly influencing early clonal expansion and malignant evolution, independently of inducing mutations. Additionally, clonal evolution in tumorigenesis reflects a multifaceted interplay between cell-intrinsic identities and various cell-extrinsic factors that exert selective pressures to either restrain uncontrolled proliferation or allow specific clones to progress into tumors. However, the mechanisms by which driver events induce both intrinsic cellular competency and remodel environmental stress to facilitate malignant transformation are not fully understood. In this review, we summarize the genetic, epigenetic, and external driver events, and their effects on the co-evolution of the transformed cells and their ecosystem during tumor initiation and early malignant evolution. A deeper understanding of the earliest molecular events holds promise for translational applications, predicting individuals at high-risk of tumor and developing strategies to intercept malignant transformation.

Malignant Melanoma: An Overview, New Perspectives, and Vitamin D Signaling

Melanoma, originating through malignant transformation of melanin-producing melanocytes, is a formidable malignancy, characterized by local invasiveness, recurrence, early metastasis, resistance to therapy, and a high mortality rate. This review discusses etiologic and risk factors for melanoma, diagnostic and prognostic tools, including recent advances in molecular biology, omics, and bioinformatics, and provides an overview of its therapy. Since the incidence of melanoma is rising and mortality remains unacceptably high, we discuss its inherent properties, including melanogenesis, that make this disease resilient to treatment and propose to use AI to solve the above complex and multidimensional problems. We provide an overview on vitamin D and its anticancerogenic properties, and report recent advances in this field that can provide solutions for the prevention and/or therapy of melanoma. Experimental papers and clinicopathological studies on the role of vitamin D status and signaling pathways initiated by its active metabolites in melanoma prognosis and therapy are reviewed. We conclude that vitamin D signaling, defined by specific nuclear receptors and selective activation by specific vitamin D hydroxyderivatives, can provide a benefit for new or existing therapeutic approaches. We propose to target vitamin D signaling with the use of computational biology and AI tools to provide a solution to the melanoma problem.

Receptor tyrosine kinase inhibition leads to regression of acral melanoma by targeting the tumor microenvironment

Acral melanoma (AM) is an aggressive melanoma variant that arises from palmar, plantar, and nail unit melanocytes. Compared to non-acral cutaneous melanoma (CM), AM is biologically distinct, has an equal incidence across genetic ancestries, typically presents in advanced stage disease, is less responsive to therapy, and has an overall worse prognosis. Independent analysis of published genomic and transcriptomic sequencing identified that receptor tyrosine kinase (RTK) ligands and adapter proteins are frequently amplified, translocated, and/or overexpressed in AM. To target these unique genetic changes, a zebrafish acral melanoma model was exposed to a panel of narrow and broad spectrum multi-RTK inhibitors, revealing that dual FGFR/VEGFR inhibitors decrease acral-analogous melanocyte proliferation and migration. The potent pan-FGFR/VEGFR inhibitor, Lenvatinib, uniformly induces tumor regression in AM patient-derived xenograft (PDX) tumors but only slows tumor growth in CM models. Unlike other multi-RTK inhibitors, Lenvatinib is not directly cytotoxic to dissociated AM PDX tumor cells and instead disrupts tumor architecture and vascular networks. Considering the great difficulty in establishing AM cell culture lines, these findings suggest that AM may be more sensitive to microenvironment perturbations than CM. In conclusion, dual FGFR/VEGFR inhibition may be a viable therapeutic strategy that targets the unique biology of AM.

CSPG4-targeting CAR-macrophages inhibit melanoma growth

Chimeric antigen receptor (CAR) T-cell therapy has revolutionized the treatment of hematological malignancies but has been clinically less effective in solid tumors. Engineering macrophages with CARs has emerged as a promising approach to overcome some of the challenges faced by CAR-T cells due to the macrophage's ability to easily infiltrate tumors, phagocytose their targets, and reprogram the immune response. We engineered CAR-macrophages (CAR-Ms) to target chondroitin sulfate proteoglycan 4 (CSPG4), an antigen expressed in melanoma, and several other solid tumors. CSPG4-targeting CAR-Ms exhibited specific phagocytosis of CSPG4-expressing melanoma cells. Combining CSPG4-targeting CAR-Ms with CD47 blocking antibodies synergistically enhanced CAR-M-mediated phagocytosis and effectively inhibited melanoma spheroid growth in 3D. Furthermore, CSPG4-targeting CAR-Ms inhibited melanoma tumor growth in mouse models. These results suggest that CSPG4-targeting CAR-M immunotherapy is a promising solid tumor immunotherapy approach for treating melanoma.

STATEMENT OF SIGNIFICANCE

We engineered macrophages with CARs as an alternative approach for solid tumor treatment. CAR-macrophages (CAR-Ms) targeting CSPG4, an antigen expressed in melanoma and other solid tumors, phagocytosed melanoma cells and inhibited melanoma growth in vivo . Thus, CSPG4-targeting CAR-Ms may be a promising strategy to treat patients with CSPG4-expressing tumors.

Unraveling the Functional Heterogeneity of Human Skin at Single-Cell Resolution

The skin consists of several cell populations, including epithelial, immune, and stromal cells. Recently, there has been a significant increase in single-cell RNA-sequencing studies, contributing to the development of a consensus Human Skin Cell Atlas. The aim is to understand skin biology better and identify potential therapeutic targets. The present review utilized previously published single-cell RNA-sequencing datasets to explore human skin's cellular and functional heterogeneity. Additionally, it summarizes the functional significance of newly identified cell subpopulations in processes such as wound healing and aging.

Fibroblast Subpopulations in Systemic Sclerosis: Functional Implications of Individual Subpopulations and Correlations with Clinical Features

Fibroblasts constitute a heterogeneous population of cells. In this study, we integrated single-cell RNA-sequencing and bulk RNA-sequencing data as well as clinical information to study the role of individual fibroblast populations in systemic sclerosis (SSc). SSc skin demonstrated an increased abundance of COMP+, COL11A1+, MYOC+, CCL19+, SFRP4/SFRP2+, and PRSS23/SFRP2+ fibroblasts signatures and decreased proportions of CXCL12+ and PI16+ fibroblast signatures in the Prospective Registry of Early Systemic Sclerosis and Genetics versus Environment in Scleroderma Outcome Study cohorts. Numerical differences were confirmed by multicolor immunofluorescence for selected fibroblast populations. COMP+, COL11A1+, SFRP4/SFRP2+, PRSS23/SFRP2+, and PI16+ fibroblasts were similarly altered between normal wound healing and patients with SSc. The proportions of profibrotic COMP+, COL11A1+, SFRP4/SFRP2+, and PRSS23/SFRP2+ and proinflammatory CCL19+ fibroblast signatures were positively correlated with clinical and histopathological parameters of skin fibrosis, whereas signatures of CXCL12+ and PI16+ fibroblasts were inversely correlated. Incorporating the proportions of COMP+, COL11A1+, SFRP4/SFRP2+, and PRSS23/SFRP2+ fibroblast signatures into machine learning models improved the classification of patients with SSc into those with progressive versus stable skin fibrosis. In summary, the profound imbalance of fibroblast subpopulations in SSc may drive the progression of skin fibrosis. Specific targeting of disease-relevant fibroblast populations may offer opportunities for the treatment of SSc and other fibrotic diseases.

Pulsed electric field induces exocytosis and overexpression of MAGE antigens in melanoma

Nanosecond pulsed electric field (nsPEF) has emerged as a promising approach for inducing cell death in melanoma, either as a standalone treatment or in combination with chemotherapeutics. However, to date, there has been a shortage of studies exploring the impact of nsPEF on the expression of cancer-specific molecules. In this investigation, we sought to assess the effects of nsPEF on melanoma-specific MAGE (Melanoma Antigen Gene Protein Family) expression. To achieve this, melanoma cells were exposed to nsPEF with parameters set at 8 kV/cm, 200 ns duration, 100 pulses, and a frequency of 10 kHz. We also aimed to comprehensively describe the consequences of this electric field on melanoma cells' invasion and proliferation potential. Our findings reveal that following exposure to nsPEF, melanoma cells release microvesicles containing MAGE antigens, leading to a simultaneous increase in the expression and mRNA content of membrane-associated antigens such as MAGE-A1. Notably, we observed an unexpected increase in the expression of PD-1 as well. While we did not observe significant differences in the cells' proliferation or invasion potential, a remarkable alteration in the cells' metabolomic and lipidomic profiles towards a less aggressive phenotype was evident. Furthermore, we validated these results using ex vivo tissue cultures and 3D melanoma culture models. Our study demonstrates that nsPEF can elevate the expression of membrane-associated proteins, including melanoma-specific antigens. The mechanism underlying the overexpression of MAGE antigens involves the initial release of microvesicles containing MAGE antigens, followed by a gradual increase in mRNA levels, ultimately resulting in elevated expression of MAGE antigens post-experiment. These findings shed light on a novel method for modulating cancer cells to overexpress cancer-specific molecules, thereby potentially enhancing their sensitivity to targeted anticancer therapy.

Deciphering cutaneous melanoma prognosis through LDL metabolism: Single-cell transcriptomics analysis via 101 machine learning algorithms

Cutaneous melanoma poses a formidable challenge within the field of oncology, marked by its aggressive nature and capacity for metastasis. Despite extensive research uncovering numerous genetic and molecular contributors to cutaneous melanoma development, there remains a critical knowledge gap concerning the role of lipids, notably low-density lipoprotein (LDL), in this lethal skin cancer. This article endeavours to bridge this knowledge gap by delving into the intricate interplay between LDL metabolism and cutaneous melanoma, shedding light on how lipids influence tumour progression, immune responses and potential therapeutic avenues. Genes associated with LDL metabolism were extracted from the GSEA database. We acquired and analysed single-cell sequencing data (GSE215120) and bulk-RNA sequencing data, including the TCGA data set, GSE19234, GSE22153 and GSE65904. Our analysis unveiled the heterogeneity of LDL across various cell types at the single-cell sequencing level. Additionally, we constructed an LDL-related signature (LRS) using machine learning algorithms, incorporating differentially expressed genes and highly correlated genes. The LRS serves as a valuable tool for assessing the prognosis, immunity and mutation status of patients with cutaneous melanoma. Furthermore, we conducted experiments on A375 and WM-115 cells to validate the function of PPP2R1A, a pivotal gene within the LRS. Our comprehensive approach, combining advanced bioinformatics analyses with an extensive review of current literature, presents compelling evidence regarding the significance of LDL within the cutaneous melanoma microenvironment.

Interleukin-4 induced 1-mediated resistance to an immune checkpoint inhibitor through suppression of CD8+ T cell infiltration in melanoma

Cancer cells adopt multiple strategies to escape tumor surveillance by the host immune system and aberrant amino acid metabolism in the tumor microenvironment suppresses the immune system. Among the amino acid-metabolizing enzymes is an L-amino-acid oxidase called interleukin-4 induced 1 (IL4I1), which depletes essential amino acids in immune cells and is associated with a poor prognosis in various cancer types. Although IL4I1 is involved in immune metabolism abnormalities, its effect on the therapeutic efficacy of immune checkpoint inhibitors is unknown. In this study, we established murine melanoma cells overexpressing IL4I1 and investigated their effects on the intratumor immune microenvironment and the antitumor efficacy of anti-programmed death-ligand 1 (PD-L1) antibodies (Abs) in a syngeneic mouse model. As a result, we found that IL4I1-overexpressing B16-F10-derived tumors showed resistance to anti-PD-L1 Ab therapy. Transcriptome analysis revealed that immunosuppressive genes were globally upregulated in the IL4I1-overexpressing tumors. Consistently, we showed that IL4I1-overexpressing tumors exhibited an altered subset of lymphoid cells and particularly significant suppression of cytotoxic T cell infiltration compared to mock-infected B16-F10-derived tumors. After treatment with anti-PD-L1 Abs, we also found a more prominent elevation of tumor-associated macrophage (TAM) marker, CD68, in the IL4I1-overexpressing tumors than in the mock tumors. Consistently, we confirmed an enhanced TAM infiltration in the IL4I1-overexpressing tumors and a functional involvement of TAMs in the tumor growth. These observations indicate that IL4I1 reprograms the tumor microenvironment into an immunosuppressive state and thereby confers resistance to anti-PD-L1 Abs.

ERK hyperactivation serves as a unified mechanism of escape in intrinsic and acquired CDK4/6 inhibitor resistance in acral lentiginous melanoma

Patients with metastatic acral lentiginous melanoma (ALM) suffer worse outcomes relative to patients with other forms of cutaneous melanoma (CM), and do not benefit as well to approved melanoma therapies. Identification of cyclin-dependent kinase 4 and 6 (CDK4/6) pathway gene alterations in >60% of ALMs has led to clinical trials of the CDK4/6 inhibitor (CDK4i/6i) palbociclib for ALM; however, median progression free survival with CDK4i/6i treatment was only 2.2 months, suggesting existence of resistance mechanisms. Therapy resistance in ALM remains poorly understood; here we report hyperactivation of MAPK signaling and elevated cyclin D1 expression serve as a mechanism of intrinsic early/adaptive CDK4i/6i resistance. ALM cells that have acquired CDK4i/6i resistance following chronic treatment exposure also exhibit hyperactivation of the MAPK pathway. MEK and/or ERK inhibition increases CDK4i/6i efficacy against therapy naïve and CDK4i/6i-resistant AM cells in xenograft and patient-derived xenograft (PDX) models and promotes a defective DNA repair, cell cycle arrested and apoptotic program. Notably, gene alterations poorly correlate with protein expression of cell cycle proteins in ALM or efficacy of CDK4i/6i, urging additional strategies when stratifying patients for CDK4i/6i trial inclusion. Concurrent targeting of the MAPK pathway and CDK4/6 represents a new approach for patients with metastatic ALM to improve outcomes.

Cell state dependent effects of Bmal1 on melanoma immunity and tumorigenicity

The circadian clock regulator Bmal1 modulates tumorigenesis, but its reported effects are inconsistent. Here, we show that Bmal1 has a context-dependent role in mouse melanoma tumor growth. Loss of Bmal1 in YUMM2.1 or B16-F10 melanoma cells eliminates clock function and diminishes hypoxic gene expression and tumorigenesis, which could be rescued by ectopic expression of HIF1α in YUMM2.1 cells. By contrast, over-expressed wild-type or a transcriptionally inactive mutant Bmal1 non-canonically sequester myosin heavy chain 9 (Myh9) to increase MRTF-SRF activity and AP-1 transcriptional signature, and shift YUMM2.1 cells from a Sox10high to a Sox9high immune resistant, mesenchymal cell state that is found in human melanomas. Our work describes a link between Bmal1, Myh9, mouse melanoma cell plasticity, and tumor immunity. This connection may underlie cancer therapeutic resistance and underpin the link between the circadian clock, MRTF-SRF and the cytoskeleton.

Differentiation States of Phenotypic Transition of Melanoma Cells Are Revealed by 3D Cell Cultures

Melanoma is characterized by high metastatic potential favored by the epithelial-to-mesenchymal transition (EMT), leading melanoma cells to exhibit a spectrum of typical EMT markers. This study aimed to analyze the expression of EMT markers in A375 and BLM melanoma cell lines cultured in 2D monolayers and 3D spheroids using morphological and molecular methods. The expression of EMT markers was strongly affected by 3D arrangement and revealed a hybrid phenotype for the two cell lines. Indeed, although E-cadherin was almost undetectable in both A375 and BLM cells, cortical actin was detected in A375 2D monolayers and 3D spheroids and was strongly expressed in BLM 3D spheroids. The mesenchymal marker N-cadherin was significantly up-regulated in A375 3D spheroids while undetectable in BLM cells, but vimentin was similarly expressed in both cell lines at the gene and protein levels. This pattern suggests that A375 cells exhibit a more undifferentiated/mesenchymal phenotype, while BLM cells have more melanocytic/differentiated characteristics. Accordingly, the Zeb1 and 2, Slug, Snail and Twist gene expression analyses showed that they were differentially expressed in 2D monolayers compared to 3D spheroids, supporting this view. Furthermore, A375 cells are characterized by a greater invasive potential, strongly influenced by 3D arrangement, compared to the BLM cell line, as evaluated by SDS-zymography and TIMPs gene expression analysis. Finally, TGF-β1, a master controller of EMT, and lysyl oxidase (LOX), involved in melanoma progression, were strongly up-regulated by 3D arrangement in the metastatic BLM cells alone, likely playing a role in the metastatic phases of melanoma progression. Overall, these findings suggest that A375 and BLM cells possess a hybrid/intermediate phenotype in relation to the expression of EMT markers. The former is characterized by a more mesenchymal/undifferentiated phenotype, while the latter shows a more melanocytic/differentiated phenotype. Our results contribute to the characterization of the role of EMT in melanoma cells and confirm that a 3D cell culture model could provide deeper insight into our understanding of the biology of melanoma.

Delineating the early dissemination mechanisms of acral melanoma by integrating single-cell and spatial transcriptomic analyses

Acral melanoma (AM) is a rare subtype of melanoma characterized by a high incidence of lymph node (LN) metastasis, a critical factor in tumor dissemination and therapeutic decision-making. Here, we employ single-cell and spatial transcriptomic analyses to investigate the dynamic evolution of early AM dissemination. Our findings reveal substantial inter- and intra-tumor heterogeneity in AM, alongside a highly immunosuppressive tumor microenvironment and complex intercellular communication networks, particularly in patients with LN metastasis. Notably, we identify a strong association between MYC+ Melanoma (MYC+MEL) and FGFBP2+NKT cells with LN metastasis. Furthermore, we demonstrate that LN metastasis requires a metabolic shift towards fatty acid oxidation (FAO) induced by MITF in MYC+MEL cells. Etomoxir, a clinically approved FAO inhibitor, can effectively suppress MITF-mediated LN metastasis. This comprehensive dataset enhances our understanding of LN metastasis in AM, and provides insights into the potential therapeutic targeting for the management of early AM dissemination.

Canonical Wnt and TGF-β/BMP signaling enhance melanocyte regeneration but suppress invasiveness, migration, and proliferation of melanoma cells

Melanoma is the deadliest form of skin cancer and develops from the melanocytes that are responsible for the pigmentation of the skin. The skin is also a highly regenerative organ, harboring a pool of undifferentiated melanocyte stem cells that proliferate and differentiate into mature melanocytes during regenerative processes in the adult. Melanoma and melanocyte regeneration share remarkable cellular features, including activation of cell proliferation and migration. Yet, melanoma considerably differs from the regenerating melanocytes with respect to abnormal proliferation, invasive growth, and metastasis. Thus, it is likely that at the cellular level, melanoma resembles early stages of melanocyte regeneration with increased proliferation but separates from the later melanocyte regeneration stages due to reduced proliferation and enhanced differentiation. Here, by exploiting the zebrafish melanocytes that can efficiently regenerate and be induced to undergo malignant melanoma, we unravel the transcriptome profiles of the regenerating melanocytes during early and late regeneration and the melanocytic nevi and malignant melanoma. Our global comparison of the gene expression profiles of melanocyte regeneration and nevi/melanoma uncovers the opposite regulation of a substantial number of genes related to Wnt signaling and transforming growth factor beta (TGF-β)/(bone morphogenetic protein) BMP signaling pathways between regeneration and cancer. Functional activation of canonical Wnt or TGF-β/BMP pathways during melanocyte regeneration promoted melanocyte regeneration but potently suppressed the invasiveness, migration, and proliferation of human melanoma cells in vitro and in vivo. Therefore, the opposite regulation of signaling mechanisms between melanocyte regeneration and melanoma can be exploited to stop tumor growth and develop new anti-cancer therapies.

Disrupting cellular memory to overcome drug resistance

Gene expression states persist for varying lengths of time at the single-cell level, a phenomenon known as gene expression memory. When cells switch states, losing memory of their prior state, this transition can occur in the absence of genetic changes. However, we lack robust methods to find regulators of memory or track state switching. Here, we develop a lineage tracing-based technique to quantify memory and identify cells that switch states. Applied to melanoma cells without therapy, we quantify long-lived fluctuations in gene expression that are predictive of later resistance to targeted therapy. We also identify the PI3K and TGF-β pathways as state switching modulators. We propose a pretreatment model, first applying a PI3K inhibitor to modulate gene expression states, then applying targeted therapy, which leads to less resistance than targeted therapy alone. Together, we present a method for finding modulators of gene expression memory and their associated cell fates.

GABA Regulates Electrical Activity and Tumor Initiation in Melanoma

Oncogenes can initiate tumors only in certain cellular contexts, which is referred to as oncogenic competence. In melanoma, whether cells in the microenvironment can endow such competence remains unclear. Using a combination of zebrafish transgenesis coupled with human tissues, we demonstrate that GABAergic signaling between keratinocytes and melanocytes promotes melanoma initiation by BRAFV600E. GABA is synthesized in melanoma cells, which then acts on GABA-A receptors in keratinocytes. Electron microscopy demonstrates specialized cell-cell junctions between keratinocytes and melanoma cells, and multielectrode array analysis shows that GABA acts to inhibit electrical activity in melanoma/keratinocyte cocultures. Genetic and pharmacologic perturbation of GABA synthesis abrogates melanoma initiation in vivo. These data suggest that GABAergic signaling across the skin microenvironment regulates the ability of oncogenes to initiate melanoma.

SIGNIFICANCE

This study shows evidence of GABA-mediated regulation of electrical activity between melanoma cells and keratinocytes, providing a new mechanism by which the microenvironment promotes tumor initiation. This provides insights into the role of the skin microenvironment in early melanomas while identifying GABA as a potential therapeutic target in melanoma. See related commentary by Ceol, p. 2128. This article is featured in Selected Articles from This Issue, p. 2109.

Dedifferentiated and undifferentiated melanomas: a practical approach to a challenging diagnosis

Dedifferentiation, defined as the loss of cellular features of terminal differentiation resulting in a more primitive, unspecialized phenotype, is a rare phenomenon in melanoma and more commonly observed in the metastatic setting than in primary melanomas. The diagnosis of dedifferentiated melanoma poses a significant challenge, and the tumors need to be carefully sampled and worked up to identify any residual morphologic, immunohistochemical, or molecular evidence of their melanocytic lineage. This article reviews our current knowledge of the clinical, histopathological, immunohistochemical, and molecular features of these rare tumors and provides a practical diagnostic approach and discussion of the relevant differential diagnoses and associated diagnostic pitfalls.

Robust dimethyl-based multiplex-DIA doubles single-cell proteome depth via a reference channel

Single-cell proteomics aims to characterize biological function and heterogeneity at the level of proteins in an unbiased manner. It is currently limited in proteomic depth, throughput, and robustness, which we address here by a streamlined multiplexed workflow using data-independent acquisition (mDIA). We demonstrate automated and complete dimethyl labeling of bulk or single-cell samples, without losing proteomic depth. Lys-N digestion enables five-plex quantification at MS1 and MS2 level. Because the multiplexed channels are quantitatively isolated from each other, mDIA accommodates a reference channel that does not interfere with the target channels. Our algorithm RefQuant takes advantage of this and confidently quantifies twice as many proteins per single cell compared to our previous work (Brunner et al, PMID 35226415), while our workflow currently allows routine analysis of 80 single cells per day. Finally, we combined mDIA with spatial proteomics to increase the throughput of Deep Visual Proteomics seven-fold for microdissection and four-fold for MS analysis. Applying this to primary cutaneous melanoma, we discovered proteomic signatures of cells within distinct tumor microenvironments, showcasing its potential for precision oncology.

MacroH2A restricts inflammatory gene expression in melanoma cancer-associated fibroblasts by coordinating chromatin looping

MacroH2A has established tumour suppressive functions in melanoma and other cancers, but an unappreciated role in the tumour microenvironment. Using an autochthonous, immunocompetent mouse model of melanoma, we demonstrate that mice devoid of macroH2A variants exhibit increased tumour burden compared with wild-type counterparts. MacroH2A-deficient tumours accumulate immunosuppressive monocytes and are depleted of functional cytotoxic T cells, characteristics consistent with a compromised anti-tumour response. Single cell and spatial transcriptomics identify increased dedifferentiation along the neural crest lineage of the tumour compartment and increased frequency and activation of cancer-associated fibroblasts following macroH2A loss. Mechanistically, macroH2A-deficient cancer-associated fibroblasts display increased myeloid chemoattractant activity as a consequence of hyperinducible expression of inflammatory genes, which is enforced by increased chromatin looping of their promoters to enhancers that gain H3K27ac. In summary, we reveal a tumour suppressive role for macroH2A variants through the regulation of chromatin architecture in the tumour stroma with potential implications for human melanoma.

A Roadmap for a Consensus Human Skin Cell Atlas and Single-Cell Data Standardization

Single-cell technologies have become essential to driving discovery in both basic and translational investigative dermatology. Despite the multitude of available datasets, a central reference atlas of normal human skin, which can serve as a reference resource for skin cell types, cell states, and their molecular signatures, is still lacking. For any such atlas to receive broad acceptance, participation by many investigators during atlas construction is an essential prerequisite. As part of the Human Cell Atlas project, we have assembled a Skin Biological Network to build a consensus Human Skin Cell Atlas and outline a roadmap toward that goal. We define the drivers of skin diversity to be considered when selecting sequencing datasets for the atlas and list practical hurdles during skin sampling that can result in data gaps and impede comprehensive representation and technical considerations for tissue processing and computational analysis, the accounting for which should minimize biases in cell type enrichments and exclusions and decrease batch effects. By outlining our goals for Atlas 1.0, we discuss how it will uncover new aspects of skin biology.

Atlas-scale single-cell multi-sample multi-condition data integration using scMerge2

The recent emergence of multi-sample multi-condition single-cell multi-cohort studies allows researchers to investigate different cell states. The effective integration of multiple large-cohort studies promises biological insights into cells under different conditions that individual studies cannot provide. Here, we present scMerge2, a scalable algorithm that allows data integration of atlas-scale multi-sample multi-condition single-cell studies. We have generalized scMerge2 to enable the merging of millions of cells from single-cell studies generated by various single-cell technologies. Using a large COVID-19 data collection with over five million cells from 1000+ individuals, we demonstrate that scMerge2 enables multi-sample multi-condition scRNA-seq data integration from multiple cohorts and reveals signatures derived from cell-type expression that are more accurate in discriminating disease progression. Further, we demonstrate that scMerge2 can remove dataset variability in CyTOF, imaging mass cytometry and CITE-seq experiments, demonstrating its applicability to a broad spectrum of single-cell profiling technologies.

A single-cell level comparison of human inner ear organoids with the human cochlea and vestibular organs

Inner ear disorders are among the most common congenital abnormalities; however, current tissue culture models lack the cell type diversity to study these disorders and normal otic development. Here, we demonstrate the robustness of human pluripotent stem cell-derived inner ear organoids (IEOs) and evaluate cell type heterogeneity by single-cell transcriptomics. To validate our findings, we construct a single-cell atlas of human fetal and adult inner ear tissue. Our study identifies various cell types in the IEOs including periotic mesenchyme, type I and type II vestibular hair cells, and developing vestibular and cochlear epithelium. Many genes linked to congenital inner ear dysfunction are confirmed to be expressed in these cell types. Additional cell-cell communication analysis within IEOs and fetal tissue highlights the role of endothelial cells on the developing sensory epithelium. These findings provide insights into this organoid model and its potential applications in studying inner ear development and disorders.

Lipid droplets are a metabolic vulnerability in melanoma

Melanoma exhibits numerous transcriptional cell states including neural crest-like cells as well as pigmented melanocytic cells. How these different cell states relate to distinct tumorigenic phenotypes remains unclear. Here, we use a zebrafish melanoma model to identify a transcriptional program linking the melanocytic cell state to a dependence on lipid droplets, the specialized organelle responsible for lipid storage. Single-cell RNA-sequencing of these tumors show a concordance between genes regulating pigmentation and those involved in lipid and oxidative metabolism. This state is conserved across human melanoma cell lines and patient tumors. This melanocytic state demonstrates increased fatty acid uptake, an increased number of lipid droplets, and dependence upon fatty acid oxidative metabolism. Genetic and pharmacologic suppression of lipid droplet production is sufficient to disrupt cell cycle progression and slow melanoma growth in vivo. Because the melanocytic cell state is linked to poor outcomes in patients, these data indicate a metabolic vulnerability in melanoma that depends on the lipid droplet organelle.

Transcriptomic Analysis of the Ocular Posterior Segment Completes a Cell Atlas of the Human Eye

Although the visual system extends through the brain, most vision loss originates from defects in the eye. Its central element is the neural retina, which senses light, processes visual signals, and transmits them to the rest of the brain through the optic nerve (ON). Surrounding the retina are numerous other structures, conventionally divided into anterior and posterior segments. Here we used high-throughput single nucleus RNA sequencing (snRNA-seq) to classify and characterize cells in the extraretinal components of the posterior segment: ON, optic nerve head (ONH), peripheral sclera, peripapillary sclera (PPS), choroid, and retinal pigment epithelium (RPE). Defects in each of these tissues are associated with blinding diseases - for example, glaucoma (ONH and PPS), optic neuritis (ON), retinitis pigmentosa (RPE), and age-related macular degeneration (RPE and choroid). From ∼151,000 single nuclei, we identified 37 transcriptomically distinct cell types, including multiple types of astrocytes, oligodendrocytes, fibroblasts, and vascular endothelial cells. Our analyses revealed a differential distribution of many cell types among distinct structures. Together with our previous analyses of the anterior segment and retina, the new data complete a "Version 1" cell atlas of the human eye. We used this atlas to map the expression of >180 genes associated with the risk of developing glaucoma, which is known to involve ocular tissues in both anterior and posterior segments as well as neural retina. Similar methods can be used to investigate numerous additional ocular diseases, many of which are currently untreatable.

Critical Considerations for Investigating MicroRNAs during Tumorigenesis: A Case Study in Conceptual and Contextual Nuances of miR-211-5p in Melanoma

MicroRNAs are non-coding RNAs fundamental to metazoan development and disease. Although the aberrant regulation of microRNAs during mammalian tumorigenesis is well established, investigations into the contributions of individual microRNAs are wrought with conflicting observations. The underlying cause of these inconsistencies is often attributed to context-specific functions of microRNAs. We propose that consideration of both context-specific factors, as well as underappreciated fundamental concepts of microRNA biology, will permit a more harmonious interpretation of ostensibly diverging data. We discuss the theory that the biological function of microRNAs is to confer robustness to specific cell states. Through this lens, we then consider the role of miR-211-5p in melanoma progression. Using literature review and meta-analyses, we demonstrate how a deep understating of domain-specific contexts is critical for moving toward a concordant understanding of miR-211-5p and other microRNAs in cancer biology.

The journey from melanocytes to melanoma

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

ERK Hyperactivation Serves as a Unified Mechanism of Escape in Intrinsic and Acquired CDK4/6 Inhibitor Resistance in Acral Lentiginous Melanoma

Patients with metastatic acral lentiginous melanoma (ALM) suffer worse outcomes relative to patients with other forms of cutaneous melanoma (CM), and do not benefit as well to approved melanoma therapies. Identification of cyclin-dependent kinase 4 and 6 (CDK4/6) pathway gene alterations in > 60% of ALMs has led to clinical trials of the CDK4/6 inhibitor (CDK4i/6i) palbociclib for ALM; however, median progression free survival with CDK4i/6i treatment was only 2.2 months, suggesting existence of resistance mechanisms. Therapy resistance in ALM remains poorly understood; here we report hyperactivation of MAPK signaling and elevated cyclin D1 expression are a unified mechanism of both intrinsic and acquired CDK4i/6i resistance. MEK and/or ERK inhibition increases CDK4i/6i efficacy in a patient-derived xenograft (PDX) model of ALM and promotes a defective DNA repair, cell cycle arrested and apoptotic program. Notably, gene alterations poorly correlate with protein expression of cell cycle proteins in ALM or efficacy of CDK4i/6i, urging additional strategies when stratifying patients for CDK4i/6i trial inclusion. Concurrent targeting of the MAPK pathway and CDK4/6 represents a new approach to improve outcomes for patients with advanced ALM.

Laser-assisted drug delivery of synthetic alpha melanocyte stimulating hormone and L-tyrosine leads to increased pigmentation area and expression of melanogenesis genes in a porcine hypertrophic scar model

OBJECTIVES

One symptom of hypertrophic scar (HTS) that can develop after burn injury is dyschromia with hyper- and hypopigmentation. There are limited treatments for these conditions. Previously, we showed there is no expression of alpha melanocyte stimulating hormone (α-MSH) in hypopigmented scars, and if these melanocytes are treated with synthetic α-MSH in vitro, they respond by repigmenting. The current study tested the same hypothesis in the in vivo environment using laser-assisted drug delivery (LADD).

METHODS

HTSs were created in red Duroc pigs. At Day 77 (pre), they were treated with CO₂ fractional ablative laser (FLSR). Synthetic α-MSH was delivered as a topical solution dissolved in  l-tyrosine (n = 6, treated). Control scars received LADD of  l-tyrosine only (n = 2, control). Scars were treated and examined weekly through Week 4. Digital images and punch biopsies of hyper, hypo-, and normally pigmented scar and skin were collected. Digital pictures were analyzed with ImageJ by tracing the area of hyperpigmentation. Epidermal sheets were obtained from punch biopsies through dispase separation and RNA was isolated. qRT-PCR was run for melanogenesis-related genes: tyrosinase (TYR), tyrosinase-related protein-1 (TYRP1), and dopachrome tautomerase (DCT). Two-way ANOVA with multiple comparisons and Dunnett's correction compared the groups.

RESULTS

The areas of hyperpigmentation were variable before treatment. Therefore, data is represented as fold-change where each scar was normalized to its own pre value. Within the LADD of NDP α-MSH + l-tyrosine group, hyperpigmented areas gradually increased each week, reaching 1.3-fold over pre by Week 4. At each timepoint, area of hyperpigmentation was greater in the treated versus the control (1.04 ± 0.05 vs. 0.89 ± 0.08, 1.21 ± 0.07 vs. 0.98 ± 0.24, 1.21 ± 0.08 vs. 1.04 ± 0.11, 1.28 ± 0.11 vs. 0.94 ± 0.25; fold-change from pre-). Within the treatment group, pretreatment, levels of TYR were decreased -17.76 ± 4.52 below the level of normal skin in hypopigmented scars. After 1 treatment, potentially due to laser fractionation, the levels decreased to -43.49 ± 5.52. After 2, 3, and 4 treatments, there was ever increasing levels of TYR to almost the level of normally pigmented skin (-35.74 ± 15.72, -23.25 ± 6.80, -5.52 ± 2.22 [p < 0.01, Week 4]). This pattern was also observed for TYRP1 (pre = -12.94 ± 1.82, Week 1 = -48.85 ± 13.25 [p < 0.01], Weeks 2, 3, and 4 = -34.45 ± 14.64, -28.19 ± 4.98, -6.93 ± 3.05 [p < 0.01, Week 4]) and DCT (pre = -214.95 ± 89.42, Week 1 = -487.93 ± 126.32 [p < 0.05], Weeks 2, 3, and 4 = -219.06 ± 79.33, -72.91 ± 20.45 [p < 0.001], -76.00 ± 24.26 [p < 0.001]). Similar patterns were observed for scars treated with LADD of  l-tyrosine alone without NDP α-MSH. For each gene, in hyperpigmented scar, levels increased at Week 4 of treatment compared to Week 1 (p < 0.01).

CONCLUSIONS

A clinically-relevant FLSR treatment method can be combined with topical delivery of synthetic α-MSH and l-tyrosine to increase the area of pigmentation and expression of melanogenesis genes in hypopigmented HTS. LADD of  l-tyrosine alone leads to increased expression of melanogenesis genes. Future studies will aim to optimize drug delivery, timing, and dosing.

Loss of MC1R signaling implicates TBX3 in pheomelanogenesis and melanoma predisposition

The human Red Hair Color (RHC) trait is caused by increased pheomelanin (red-yellow) and reduced eumelanin (black-brown) pigment in skin and hair due to diminished melanocortin 1 receptor (MC1R) function. In addition, individuals harboring the RHC trait are predisposed to melanoma development. While MC1R variants have been established as causative of RHC and are a well-defined risk factor for melanoma, it remains unclear mechanistically why decreased MC1R signaling alters pigmentation and increases melanoma susceptibility. Here, we use single-cell RNA-sequencing (scRNA-seq) of melanocytes isolated from RHC mouse models to reveal a Pheomelanin Gene Signature (PGS) comprising genes implicated in melanogenesis and oncogenic transformation. We show that TBX3, a well-known anti-senescence transcription factor implicated in melanoma progression, is part of the PGS and binds both E-box and T-box elements to regulate genes associated with melanogenesis and senescence bypass. Our results provide key insights into mechanisms by which MC1R signaling regulates pigmentation and how individuals with the RHC phenotype are predisposed to melanoma.

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.

Differential cell composition and split epidermal differentiation in human palm, sole, and hip skin

Palmoplantar skin is structurally and functionally unique, but the transcriptional programs driving this specialization are unclear. Here, we use bulk and single-cell RNA sequencing of human palm, sole, and hip skin to describe the distinguishing characteristics of palmoplantar and non-palmoplantar skin while also uncovering differences between palmar and plantar sites. Our approach reveals an altered immune environment in palmoplantar skin, with downregulation of diverse immunological processes and decreased immune cell populations. Further, we identify specific fibroblast populations that appear to orchestrate key differences in cell-cell communication in palm, sole, and hip. Dedicated keratinocyte analysis highlights major differences in basal cell fraction among the three sites and demonstrates the existence of two spinous keratinocyte populations constituting parallel, site-selective epidermal differentiation trajectories. In summary, this deep characterization of highly adapted palmoplantar skin contributes key insights into the fundamental biology of human skin and provides a valuable data resource for further investigation.

Single-cell chromatin landscapes of mouse skin development

The coat of mammals is produced by hair follicles, and hair follicle is an important and complex accessory organ of skin. As a complex physiological regulation process, hair follicle morphogenesis is regulated by a series of signal pathway factors, involves the interaction of multiple cell types and begins in the early embryonic stage. However, its transcriptional regulatory mechanism is unclear. We have therefore utilized single-cell ATAC sequencing to obtain the chromatin accessibility landscapes of 6,928, 6,961 and 7,374 high-quality cells from the dorsal skins of E13.5, E16.5 and P0 mice (Mus musculus), respectively. Based on marker gene activity clustering, we defined 6, 8 and 5 distinct cell types in E13.5, E16.5 and P0 stages, respectively. Furtherly, we integrated the fibroblasts and keratinocytes clusters, performed further analysis and re-clustered. The single cell map of the chromatin open area was drawn from each cell type and the mechanism of cell transcription regulation was explored. Collectively, our data provide a reference for deeply exploring the epigenetic regulation mechanism of mouse hair follicles development.

Applications of single-cell RNA sequencing in atopic dermatitis and psoriasis

Single-cell RNA sequencing (scRNA-seq) is a novel technology that characterizes molecular heterogeneity at the single-cell level. With the development of more automated, sensitive, and cost-effective single-cell isolation methods, the sensitivity and efficiency of scRNA-seq have improved. Technological advances in single-cell analysis provide a deeper understanding of the biological diversity of cells present in tissues, including inflamed skin. New subsets of cells have been discovered among common inflammatory skin diseases, such as atopic dermatitis (AD) and psoriasis. ScRNA-seq technology has also been used to analyze immune cell distribution and cell-cell communication, shedding new light on the complex interplay of components involved in disease responses. Moreover, scRNA-seq may be a promising tool in precision medicine because of its ability to define cell subsets with potential treatment targets and to characterize cell-specific responses to drugs or other stimuli. In this review, we briefly summarize the progress in the development of scRNA-seq technologies and discuss the latest scRNA-seq-related findings and future trends in AD and psoriasis. We also discuss the limitations and technical problems associated with current scRNA-seq technology.

Tumor-intrinsic SIRPA promotes sensitivity to checkpoint inhibition immunotherapy in melanoma

Checkpoint inhibition immunotherapy has revolutionized cancer treatment, but many patients show resistance. Here we perform integrative transcriptomic and proteomic analyses on emerging immuno-oncology targets across multiple clinical cohorts of melanoma under anti-PD-1 treatment, on both bulk and single-cell levels. We reveal a surprising role of tumor-intrinsic SIRPA in enhancing antitumor immunity, in contrast to its well-established role as a major inhibitory immune modulator in macrophages. The loss of SIRPA expression is a marker of melanoma dedifferentiation, a key phenotype linked to immunotherapy efficacy. Inhibition of SIRPA in melanoma cells abrogates tumor killing by activated CD8+ T cells in a co-culture system. Mice bearing SIRPA-deficient melanoma tumors show no response to anti-PD-L1 treatment, whereas melanoma-specific SIRPA overexpression significantly enhances immunotherapy response. Mechanistically, SIRPA is regulated by its pseudogene, SIRPAP1. Our results suggest a complicated role of SIRPA in the tumor ecosystem, highlighting cell-type-dependent antagonistic effects of the same target on immunotherapy.

SRC-RAC1 signaling drives drug resistance to BRAF inhibition in de-differentiated cutaneous melanomas

Rare gain-of-function mutations in RAC1 drive drug resistance to targeted BRAF inhibition in cutaneous melanoma. Here, we show that wildtype RAC1 is a critical driver of growth and drug resistance, but only in a subset of melanomas with elevated markers of de-differentiation. Similarly, SRC inhibition also selectively sensitized de-differentiated melanomas to BRAF inhibition. One possible mechanism may be the suppression of the de-differentiated state, as SRC and RAC1 maintained markers of de-differentiation in human melanoma cells. The functional differences between melanoma subtypes suggest that the clinical management of cutaneous melanoma can be enhanced by the knowledge of differentiation status. To simplify the task of classification, we developed a binary classification strategy based on a small set of ten genes. Using this gene set, we reliably determined the differentiation status previously defined by hundreds of genes. Overall, our study informs strategies that enhance the precision of BRAFi by discovering unique vulnerabilities of the de-differentiated cutaneous melanoma subtype and creating a practical method to resolve differentiation status.

AP-1 transcription factor network explains diverse patterns of cellular plasticity in melanoma cells

Cellular plasticity associated with fluctuations in transcriptional programs allows individual cells in a tumor to adopt heterogeneous differentiation states and switch phenotype during their adaptive responses to therapies. Despite increasing knowledge of such transcriptional programs, the molecular basis of cellular plasticity remains poorly understood. Here, we combine multiplexed transcriptional and protein measurements at population and single-cell levels with multivariate statistical modeling to show that the state of AP-1 transcription factor network plays a unifying role in explaining diverse patterns of plasticity in melanoma. We find that a regulated balance among AP-1 factors cJUN, JUND, FRA2, FRA1, and cFOS determines the intrinsic diversity of differentiation states and adaptive responses to MAPK inhibitors in melanoma cells. Perturbing this balance through genetic depletion of specific AP-1 proteins, or by MAPK inhibitors, shifts cellular heterogeneity in a predictable fashion. Thus, AP-1 may serve as a critical node for manipulating cellular plasticity with potential therapeutic implications.

A robust experimental and computational analysis framework at multiple resolutions, modalities and coverages

The ability to study cancer-immune cell communication across the whole tumor section without tissue dissociation is needed, especially for cancer immunotherapy development, which requires understanding of molecular mechanisms and discovery of more druggable targets. In this work, we assembled and evaluated an integrated experimental framework and analytical process to enable genome-wide scale discovery of ligand-receptors potentially used for cellular crosstalks, followed by targeted validation. We assessed the complementarity of four different technologies: single-cell RNA sequencing and Spatial transcriptomic (measuring over >20,000 genes), RNA In Situ Hybridization (RNAscope, measuring 4-12 genes) and Opal Polaris multiplex protein staining (4-9 proteins). To utilize the multimodal data, we implemented existing methods and also developed STRISH (Spatial TRanscriptomic In Situ Hybridization), a computational method that can automatically scan across the whole tissue section for local expression of gene (e.g. RNAscope data) and/or protein markers (e.g. Polaris data) to recapitulate an interaction landscape across the whole tissue. We evaluated the approach to discover and validate cell-cell interaction in situ through in-depth analysis of two types of cancer, basal cell carcinoma and squamous cell carcinoma, which account for over 70% of cancer cases. We showed that inference of cell-cell interactions using scRNA-seq data can misdetect or detect false positive interactions. Spatial transcriptomics still suffers from misdetecting lowly expressed ligand-receptor interactions, but reduces false discovery. RNAscope and Polaris are sensitive methods for defining the location of potential ligand receptor interactions, and the STRISH program can determine the probability that local gene co-expression reflects true cell-cell interaction. We expect that the approach described here will be widely applied to discover and validate ligand receptor interaction in different types of solid cancer tumors.

Analysis of Melanoma Gene Expression Signatures at the Single-Cell Level Uncovers 45-Gene Signature Related to Prognosis

Since the current melanoma clinicopathological staging system remains restricted to predicting survival outcomes, establishing precise prognostic targets is needed. Here, we used gene expression signature (GES) classification and Cox regression analyses to biologically characterize melanoma cells at the single-cell level and construct a prognosis-related gene signature for melanoma. By analyzing publicly available scRNA-seq data, we identified six distinct GESs (named: “Anti-apoptosis”, “Immune cell interactions”, “Melanogenesis”, “Ribosomal biogenesis”, “Extracellular structure organization”, and “Epithelial-Mesenchymal Transition (EMT)”). We verified these GESs in the bulk RNA-seq data of patients with skin cutaneous melanoma (SKCM) from The Cancer Genome Atlas (TCGA). Four GESs (“Immune cell interactions”, “Melanogenesis”, “Ribosomal biogenesis”, and “Extracellular structure organization”) were significantly correlated with prognosis (p = 1.08 × 10−5, p = 0.042, p = 0.001, and p = 0.031, respectively). We identified a prognostic signature of melanoma composed of 45 genes (MPS_45). MPS_45 was validated in TCGA-SKCM (HR = 1.82, p = 9.08 × 10−6) and three other melanoma datasets (GSE65904: HR = 1.73, p = 0.006; GSE19234: HR = 3.83, p = 0.002; and GSE53118: HR = 1.85, p = 0.037). MPS_45 was independently associated with survival (p = 0.002) and was proved to have a high potential for predicting prognosis in melanoma patients.

Integrated genomic analyses of acral and mucosal melanomas nominate novel driver genes

BACKGROUND

Acral and mucosal melanomas are aggressive subtypes of melanoma, which have a significantly lower burden of somatic mutations than cutaneous melanomas, but more frequent copy number variations, focused gene amplifications, and structural alterations. The landscapes of their genomic alterations remain to be fully characterized.

METHODS

We compiled sequencing data of 240 human acral and mucosal melanoma samples from 11 previously published studies and applied a uniform pipeline to call tumor cell content, ploidy, somatic and germline mutations, as well as CNVs, LOH, and SVs. We identified genes that are significantly mutated or recurrently affected by CNVs and implicated in oncogenesis. We further examined the difference in the frequency of recurrent pathogenic alterations between the two melanoma subtypes, correlation between pathogenic alterations, and their association with clinical features.

RESULTS

We nominated PTPRJ, mutated and homozygously deleted in 3.8% (9/240) and 0.8% (2/240) of samples, respectively, as a probable tumor suppressor gene, and FER and SKP2, amplified in 3.8% and 11.7% of samples, respectively, as probable oncogenes. We further identified a long tail of infrequent pathogenic alterations, involving genes such as CIC and LZTR1. Pathogenic germline mutations were observed on MITF, PTEN, ATM, and PRKN. We found BRAF V600E mutations in acral melanomas with fewer structural variations, suggesting that they are distinct and related to cutaneous melanomas. Amplifications of PAK1 and GAB2 were more commonly observed in acral melanomas, whereas SF3B1 R625 codon mutations were unique to mucosal melanomas (12.9%). Amplifications at 11q13-14 were frequently accompanied by fusion to a region on chromosome 6q12, revealing a recurrent novel structural rearrangement whose role remains to be elucidated.

CONCLUSIONS

Our meta-analysis expands the catalog of driver mutations in acral and mucosal melanomas, sheds new light on their pathogenesis and broadens the catalog of therapeutic targets for these difficult-to-treat cancers.

Stepwise-edited, human melanoma models reveal mutations' effect on tumor and microenvironment

Establishing causal relationships between genetic alterations of human cancers and specific phenotypes of malignancy remains a challenge. We sequentially introduced mutations into healthy human melanocytes in up to five genes spanning six commonly disrupted melanoma pathways, forming nine genetically distinct cellular models of melanoma. We connected mutant melanocyte genotypes to malignant cell expression programs in vitro and in vivo, replicative immortality, malignancy, rapid tumor growth, pigmentation, metastasis, and histopathology. Mutations in malignant cells also affected tumor microenvironment composition and cell states. Our melanoma models shared genotype-associated expression programs with patient melanomas, and a deep learning model showed that these models partially recapitulated genotype-associated histopathological features as well. Thus, a progressive series of genome-edited human cancer models can causally connect genotypes carrying multiple mutations to phenotype.

The Genetic Evolution of Metastasis

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

Anatomic position determines oncogenic specificity in melanoma

Oncogenic alterations to DNA are not transforming in all cellular contexts1,2. This may be due to pre-existing transcriptional programmes in the cell of origin. Here we define anatomic position as a major determinant of why cells respond to specific oncogenes. Cutaneous melanoma arises throughout the body, whereas the acral subtype arises on the palms of the hands, soles of the feet or under the nails3. We sequenced the DNA of cutaneous and acral melanomas from a large cohort of human patients and found a specific enrichment for BRAF mutations in cutaneous melanoma and enrichment for CRKL amplifications in acral melanoma. We modelled these changes in transgenic zebrafish models and found that CRKL-driven tumours formed predominantly in the fins of the fish. The fins are the evolutionary precursors to tetrapod limbs, indicating that melanocytes in these acral locations may be uniquely susceptible to CRKL. RNA profiling of these fin and limb melanocytes, when compared with body melanocytes, revealed a positional identity gene programme typified by posterior HOX13 genes. This positional gene programme synergized with CRKL to amplify insulin-like growth factor (IGF) signalling and drive tumours at acral sites. Abrogation of this CRKL-driven programme eliminated the anatomic specificity of acral melanoma. These data suggest that the anatomic position of the cell of origin endows it with a unique transcriptional state that makes it susceptible to only certain oncogenic insults.

Single-cell transcriptomics in human skin research: available technologies, technical considerations, and disease applications

Single‐cell technologies have revolutionized research in the last decade, including for skin biology. Single‐cell RNA sequencing has emerged as a powerful tool allowing the dissection of human disease pathophysiology at unprecedented resolution by assessing cell‐to‐cell variation, facilitating identification of rare cell populations and elucidating cellular heterogeneity. In dermatology, this technology has been widely applied to inflammatory skin disorders, fibrotic skin diseases, wound healing complications and cutaneous neoplasms. Here, we discuss the available technologies and technical considerations of single‐cell RNA sequencing and describe its applications to a broad spectrum of dermatological diseases.

Parallel single cell multi-omics analysis of neonatal skin reveals transitional fibroblast states that restricts differentiation into distinct fates

One of the keys to achieving skin regeneration lies within understanding the heterogeneity of neonatal fibroblasts, which support skin regeneration. However, the molecular underpinnings regulating the cellular states and fates of these cells are not fully understood. To investigate this, we performed a parallel multi-omics analysis by processing neonatal murine skin for single-cell ATAC-sequencing (scATAC-seq) and single-cell RNA-sequencing (scRNA-seq) separately. Our approach revealed that fibroblast clusters could be sorted into papillary and reticular lineages based on transcriptome profiling, as previously published. However, scATAC-seq analysis of neonatal fibroblast lineage markers, such as, Dpp4/CD26, Corin, and Dlk1 along with markers of myofibroblasts, revealed accessible chromatin in all fibroblast populations despite their lineage-specific transcriptome profiles. These results suggests that accessible chromatin does not always translate to gene expression and that many fibroblast lineage markers reflect a fibroblast state, which includes neonatal papillary, reticular, and myofibroblasts. This analysis also provides a possible explanation as to why these marker genes can be promiscuously expressed in different fibroblast populations under different conditions. Our scATAC-seq analysis also revealed that the functional lineage restriction between dermal papilla and adipocyte fates are regulated by distinct chromatin landscapes. Finally, we have developed a webtool for our multi-omics analysis: https://skinregeneration.org/scatacseq-and-scrnaseq-data-from-thompson-et-al-2021-2/.