ABSTRACT
We performed multi-omic profiling of epidermal keratinocytes, precancerous actinic keratoses, and squamous cell carcinomas to understand the molecular transitions during skin carcinogenesis. Single-cell mutational analyses showed that most keratinocytes in normal skin had lower mutation burdens than melanocytes and fibroblasts, however keratinocytes with TP53 or NOTCH1 mutations had substantially higher mutation burdens, suggesting that these mutations prime keratinocytes for transformation by increasing their mutation rate. Mutational profiling and spatial transcriptomics on squamous cell carcinomas adjacent to actinic keratoses revealed TERT promoter and CDKN2A mutations emerging in actinic keratoses, whereas additional mutations inactivating ARID2 and activating the MAPK-pathway delineated the transition to squamous cell carcinomas. Spatial variation in gene expression patterns was common in both tumor and immune cells, with high expression of checkpoint molecules at the invasive front of tumors. In conclusion, this study documents key events during the evolution of cutaneous squamous cell carcinoma.
ABSTRACT
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.
Subject(s)
Melanoma , Mutation , Skin Neoplasms , Telomerase , Humans , Melanoma/genetics , Melanoma/pathology , Skin Neoplasms/genetics , Skin Neoplasms/pathology , Telomerase/genetics , DNA Copy Number Variations , Evolution, Molecular , Male , Exome Sequencing , Female , MAP Kinase Signaling System/geneticsABSTRACT
Psoriasis vulgaris and other chronic inflammatory diseases improve markedly with therapeutic blockade of interleukin-23 (IL-23) signaling, but the genetic mechanisms underlying clinical responses remain poorly understood. Using single-cell transcriptomics, we profiled immune cells isolated from lesional psoriatic skin before and during IL-23 blockade. In clinically responsive patients, a psoriatic transcriptional signature in skin-resident memory T cells was strongly attenuated. In contrast, poorly responsive patients were distinguished by persistent activation of IL-17-producing T (T17) cells, a mechanism distinct from alternative cytokine signaling or resistance isolated to epidermal keratinocytes. Even in IL-23 blockade-responsive patients, we detected a recurring set of recalcitrant, disease-specific transcriptional abnormalities. This irreversible immunological state may necessitate ongoing IL-23 inhibition. Spatial transcriptomic analyses also suggested that successful IL-23 blockade requires dampening of >90% of IL-17-induced response in lymphocyte-adjacent keratinocytes, an unexpectedly high threshold. Collectively, our data establish a patient-level paradigm for dissecting responses to immunomodulatory treatments.
Subject(s)
Interleukin-17 , Psoriasis , Humans , Interleukin-23 , Skin , Psoriasis/drug therapy , KeratinocytesABSTRACT
To better understand intrinsic resistance to immune checkpoint blockade (ICB), we established a comprehensive view of the cellular architecture of the treatment-naive melanoma ecosystem and studied its evolution under ICB. Using single-cell, spatial multi-omics, we showed that the tumor microenvironment promotes the emergence of a complex melanoma transcriptomic landscape. Melanoma cells harboring a mesenchymal-like (MES) state, a population known to confer resistance to targeted therapy, were significantly enriched in early on-treatment biopsies from non-responders to ICB. TCF4 serves as the hub of this landscape by being a master regulator of the MES signature and a suppressor of the melanocytic and antigen presentation transcriptional programs. Targeting TCF4 genetically or pharmacologically, using a bromodomain inhibitor, increased immunogenicity and sensitivity of MES cells to ICB and targeted therapy. We thereby uncovered a TCF4-dependent regulatory network that orchestrates multiple transcriptional programs and contributes to resistance to both targeted therapy and ICB in melanoma.
Subject(s)
Melanoma , Humans , Gene Regulatory Networks , Immunotherapy , Melanocytes , Melanoma/drug therapy , Melanoma/genetics , Transcription Factor 4/genetics , Tumor MicroenvironmentABSTRACT
Spatial transcriptomic technologies, such as the Visium platform, measure gene expression in different regions of tissues. Here, we describe new software, STmut, to visualize somatic point mutations, allelic imbalance, and copy number alterations in Visium data. STmut is tested on fresh-frozen Visium data, formalin-fixed paraffin-embedded (FFPE) Visium data, and tumors with and without matching DNA sequencing data. Copy number is inferred on all conditions, but the chemistry of the FFPE platform does not permit analyses of single nucleotide variants. Taken together, we propose solutions to add the genetic dimension to spatial transcriptomic data and describe the limitations of different datatypes.
Subject(s)
Formaldehyde , Neoplasms , Humans , Transcriptome , Paraffin Embedding , Neoplasms/genetics , Gene Expression Profiling , High-Throughput Nucleotide SequencingABSTRACT
Acral melanoma is an aggressive type of melanoma with unknown origins, arising on the sole, palm, or nail apparatus. It is the most common type of melanoma in individuals with dark skin and is notoriously challenging to treat. Our study examined exome sequencing data from 139 tissue samples, spanning different progression stages, collected from 37 patients. We found that 78.4% of the melanomas displayed one or more clustered copy number transitions with focal amplifications, recurring predominantly on chromosomes 5, 11, 12, and 22. These genomic "hailstorms" were typically shared across all progression stages within individual patients. Genetic alterations known to activate TERT also arose early. By contrast, mutations in the MAP-kinase pathway appeared later during progression, often leading to different tumor areas harboring non-overlapping driver mutations. We conclude that the evolutionary trajectories of acral melanomas substantially diverge from those of melanomas on sun-exposed skin, where MAP-kinase pathway activation initiates the neoplastic cascade followed by immortalization later. The punctuated formation of hailstorms, paired with early TERT activation, suggests a unique mutational mechanism underlying the origins of acral melanoma. Our findings highlight an essential role for telomerase, likely in re-stabilizing tumor genomes after hailstorms have initiated the tumors. The marked genetic heterogeneity, in particular of MAP-kinase pathway drivers, may partly explain the limited success of targeted and other therapies in treating this melanoma subtype.
ABSTRACT
SUMMARY: Traditionally, omic studies have prioritized the number of patients over the number of tumors per patient, but in a reversal of this study design, Spain and colleagues performed the largest intrapatient analysis of melanoma to date. Their work reveals mechanisms of treatment resistance, patterns of metastatic dissemination, and new insights into the evolutionary trajectories of melanoma. See related article by Spain et al., p. 1364 (1).
Subject(s)
Melanoma , Neoplasms, Second Primary , Humans , Melanoma/genetics , Melanoma/pathology , Genomics , Biological EvolutionABSTRACT
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.
Subject(s)
Melanoma , Skin Neoplasms , DNA Copy Number Variations , Genomics , Humans , Melanoma/pathology , Mutation , Skin Neoplasms/pathology , Transcription Factors/genetics , Melanoma, Cutaneous MalignantABSTRACT
Benign melanocytic nevi frequently emerge when an acquired BRAFV600E mutation triggers unchecked proliferation and subsequent arrest in melanocytes. Recent observations have challenged the role of oncogene-induced senescence in melanocytic nevus formation, necessitating investigations into alternative mechanisms for the establishment and maintenance of proliferation arrest in nevi. We compared the transcriptomes of melanocytes from healthy human skin, nevi, and melanomas arising from nevi and identified a set of microRNAs as highly expressed nevus-enriched transcripts. Two of these microRNAs-MIR211-5p and MIR328-3p-induced mitotic failure, genome duplication, and proliferation arrest in human melanocytes through convergent targeting of AURKB. We demonstrate that BRAFV600E induces a similar proliferation arrest in primary human melanocytes that is both reversible and conditional. Specifically, BRAFV600E expression stimulates either arrest or proliferation depending on the differentiation state of the melanocyte. We report genome duplication in human melanocytic nevi, reciprocal expression of AURKB and microRNAs in nevi and melanomas, and rescue of arrested human nevus cells with AURKB expression. Taken together, our data describe an alternative molecular mechanism for melanocytic nevus formation that is congruent with both experimental and clinical observations.
Lots of people have small dark patches on their skin known as moles. Most moles form when individual cells known as melanocytes in the skin acquire a specific genetic mutation in a gene called BRAF. This mutation causes the cells to divide rapidly to form the mole. After a while, most moles stop growing and remain harmless for the rest of a person's life. Melanoma is a type of skin cancer that develops from damaged melanocytes. The same mutation in BRAF that is found in moles is also present in half of all cases of melanoma. Unlike in moles, the melanoma-causing mutation makes the melanocytes divide rapidly to form a tumor that keeps on growing indefinitely. It remains unclear why the same genetic mutation in the BRAF gene has such different consequences in moles and melanomas. To address this question, McNeal et al. used genetic approaches to study melanocytes from moles and melanomas. The experiments identified some molecules known as microRNAs that are present at higher levels in moles than in melanomas. Increasing the levels of two of these microRNAs in melanocytes from human skin stopped the cells from growing and dividing by inhibiting a gene called AURKB. This suggested that these microRNAs are responsible for halting the growth of moles. Introducing the mutated form of BRAF into melanocytes also stopped cells from growing and dividing by inhibiting AURKB. However, changing the environment surrounding the cells reversed this effect and allowed the melanocytes to resume dividing. In this way the mutated form of BRAF acts like a switch that allows melanocytes in skin cancers to start growing again under certain conditions. Further experiments found that a drug called barasertib is able to inhibit the growth of melanoma cells with the mutant form of BRAF. Future work will investigate whether it is possible to use this drug and other tools to stop skin cancer tumors from growing, and possibly even prevent skin tumors from forming in the first place.
Subject(s)
Aurora Kinase B/genetics , Melanocytes/physiology , MicroRNAs/metabolism , Mitosis/genetics , Proto-Oncogene Proteins B-raf/genetics , Aurora Kinase B/metabolism , Humans , Proto-Oncogene Proteins B-raf/metabolism , Signal TransductionABSTRACT
Cutaneous squamous cell carcinoma is a form of skin cancer originating from keratinocytes in the skin. It is the second most common type of cancer and is responsible for an estimated 8000 deaths per year in the United States. Compared to other cancer subtypes with similar incidences and death tolls, our understanding of the somatic mutations driving cutaneous squamous cell carcinoma is limited. The main challenge is that these tumors have high mutation burdens, primarily a consequence of UV-radiation-induced DNA damage from sunlight, making it difficult to distinguish driver mutations from passenger mutations. We overcame this challenge by performing a meta-analysis of publicly available sequencing data covering 105 tumors from 10 different studies. Moreover, we eliminated tumors with issues, such as low neoplastic cell content, and from the tumors that passed quality control, we utilized multiple strategies to reveal genes under selection. In total, we nominated 30 cancer genes. Among the more novel genes, mutations frequently affected EP300, PBRM1, USP28, and CHUK. Collectively, mutations in the NOTCH and p53 pathways were ubiquitous, and to a lesser extent, mutations affected genes in the Hippo pathway, genes in the Ras/MAPK/PI3K pathway, genes critical for cell-cycle checkpoint control, and genes encoding chromatin remodeling factors. Taken together, our study provides a catalog of driver genes in cutaneous squamous cell carcinoma, offering points of therapeutic intervention and insights into the biology of cutaneous squamous cell carcinoma.
ABSTRACT
BRAF-mutant melanomas are more likely than NRAS-mutant melanomas to arise in anatomical locations protected from chronic sun damage. We hypothesized that this discrepancy in tumor location is a consequence of the differential sensitivity of BRAF and NRAS-mutant melanocytes to ultraviolet light (UV)-mediated carcinogenesis. We tested this hypothesis by comparing the mutagenic consequences of a single neonatal, ultraviolet-AI (UVA; 340-400 nm) or ultraviolet-B (UVB; 280-390 nm) exposure in mouse models heterozygous for mutant Braf or homozygous for mutant Nras Tumor onset was accelerated by UVB, but not UVA, and the resulting melanomas contained recurrent mutations affecting the RING domain of MAP3K1 and Actin-binding domain of Filamin A. Melanomas from UVB-irradiated, Braf-mutant mice averaged twice as many single-nucleotide variants and five times as many dipyrimidine variants than tumors from similarly irradiated Nras-mutant mice. A mutational signature discovered in UVB-accelerated tumors mirrored COSMIC signatures associated with human skin cancer and was more prominent in Braf- than Nras-mutant murine melanomas. These data show that a single UVB exposure yields a greater burden of mutations in murine tumors driven by oncogenic Braf.
Subject(s)
Melanoma/etiology , Monomeric GTP-Binding Proteins/genetics , Mutagenesis/radiation effects , Mutation/radiation effects , Proto-Oncogene Proteins B-raf/genetics , Ultraviolet Rays/adverse effects , Animals , Biomarkers, Tumor , Disease Models, Animal , Disease Susceptibility , Genetic Predisposition to Disease , Melanoma/metabolism , Melanoma/pathology , MiceABSTRACT
Every cell in the human body has a unique set of somatic mutations, but it remains difficult to comprehensively genotype an individual cell1. Here we describe ways to overcome this obstacle in the context of normal human skin, thus offering a glimpse into the genomic landscapes of individual melanocytes from human skin. As expected, sun-shielded melanocytes had fewer mutations than sun-exposed melanocytes. However, melanocytes from chronically sun-exposed skin (for example, the face) had a lower mutation burden than melanocytes from intermittently sun-exposed skin (for example, the back). Melanocytes located adjacent to a skin cancer had higher mutation burdens than melanocytes from donors without skin cancer, implying that the mutation burden of normal skin can be used to measure cumulative sun damage and risk of skin cancer. Moreover, melanocytes from healthy skin commonly contained pathogenic mutations, although these mutations tended to be weakly oncogenic, probably explaining why they did not give rise to discernible lesions. Phylogenetic analyses identified groups of related melanocytes, suggesting that melanocytes spread throughout skin as fields of clonally related cells that are invisible to the naked eye. Overall, our results uncover the genomic landscapes of individual melanocytes, providing key insights into the causes and origins of melanoma.
Subject(s)
Genome, Human/genetics , Genomics , Health , Melanocytes/cytology , Melanoma/genetics , Single-Cell Analysis , Skin/cytology , DNA Mutational Analysis , Female , Genotype , Humans , Male , Melanocytes/metabolism , Melanocytes/pathology , Melanoma/pathology , Mutation , Skin/pathology , WorkflowABSTRACT
An amendment to this paper has been published and can be accessed via a link at the top of the paper.
ABSTRACT
The use of microRNAs as biomarkers has been proposed for many diseases, including the diagnosis of melanoma. Although hundreds of microRNAs have been identified as differentially expressed in melanomas as compared to benign melanocytic lesions, a limited consensus has been achieved across studies, constraining the effective use of these potentially useful markers. In this study, we applied a machine learning-based pipeline to a dataset consisting of genetic features, clinical features, and next-generation microRNA sequencing from micro-dissected formalin-fixed paraffin embedded melanomas and their adjacent benign precursor nevi. We identified patient age and tumor cellularity as variables that frequently confound the measured expression of potentially diagnostic microRNAs. By employing the ratios of microRNAs that were either enriched or depleted in melanoma compared to the nevi as a normalization strategy, we developed a model that classified all the available published cohorts with an area under the receiver operating characteristic curve of 0.98. External validation on an independent cohort classified lesions with 81% sensitivity and 88% specificity and was uninfluenced by the tumor content of the sample or patient age.
Subject(s)
Biomarkers, Tumor/genetics , Melanocytes/physiology , Melanoma/diagnosis , MicroRNAs/genetics , Nevus/diagnosis , Skin Neoplasms/diagnosis , Datasets as Topic , Diagnosis, Differential , High-Throughput Nucleotide Sequencing , Humans , Machine Learning , Prognosis , ROC Curve , Sensitivity and Specificity , Sequence Analysis, RNAABSTRACT
To date, over 1000 melanocytic neoplasms, spanning all stages of tumorigenesis, have been sequenced, offering detailed views into their -omic landscapes. This has coincided with advances in genetic engineering technologies that allow molecular biologists to edit the human genome with extreme precision and new mouse models to simulate disease progression. In this review, we describe how these technologies are being harnessed to provide insights into the evolution of melanoma at an unprecedented resolution, revealing that prior models of melanoma evolution, in which pathways are turned 'on' or 'off' in a binary fashion during the run-up to melanoma, are oversimplified.
Subject(s)
Carcinogenesis/genetics , MAP Kinase Signaling System/genetics , Melanocytes/pathology , Melanoma/genetics , Skin Neoplasms/genetics , Animals , Carcinogenesis/pathology , Cell Cycle Checkpoints/genetics , Cell Proliferation/genetics , DNA Mutational Analysis , Disease Models, Animal , Disease Progression , Gene Dosage , Gene Expression Regulation, Neoplastic , Genetic Engineering , Humans , Melanoma/pathology , Mice , Mice, Transgenic , Mutation , Retinoblastoma Protein/metabolism , Skin Neoplasms/pathologyABSTRACT
BRAF fusions are detected in numerous neoplasms, but their clinical management remains unresolved. We identified six melanoma lines harboring BRAF fusions representative of the clinical cases reported in the literature. Their unexpected heterogeneous responses to RAF and MEK inhibitors could be categorized upon specific features of the fusion kinases. Higher expression level correlated with resistance, and fusion partners containing a dimerization domain promoted paradoxical activation of the mitogen-activated protein kinase (MAPK) pathway and hyperproliferation in response to first- and second-generation RAF inhibitors. By contrast, next-generation αC-IN/DFG-OUT RAF inhibitors blunted paradoxical activation across all lines and had their therapeutic efficacy further increased in vitro and in vivo by combination with MEK inhibitors, opening perspectives in the clinical management of tumors harboring BRAF fusions.
Subject(s)
Drug Resistance, Neoplasm/genetics , Melanoma/pathology , Oncogene Proteins, Fusion/metabolism , Proto-Oncogene Proteins B-raf/genetics , Animals , Dimerization , Drug Resistance, Neoplasm/drug effects , Female , Gene Expression Regulation, Neoplastic/drug effects , Humans , Intracellular Signaling Peptides and Proteins/genetics , Melanoma/genetics , Mice , Mice, Nude , Mitogen-Activated Protein Kinases/metabolism , Oncogene Proteins, Fusion/antagonists & inhibitors , Oncogene Proteins, Fusion/genetics , Protein Isoforms/genetics , Protein Isoforms/metabolism , Protein Kinase Inhibitors/pharmacology , RNA Interference , RNA, Small Interfering/metabolism , Signal Transduction , Vemurafenib/pharmacology , ras Proteins/genetics , ras Proteins/metabolismABSTRACT
Uveal melanoma is a clinically distinct and particularly lethal subtype of melanoma originating from melanocytes in the eye. Here, we performed multi-region DNA sequencing of primary uveal melanomas and their matched metastases from 35 patients. We observed previously unknown driver mutations and established the order in which these and known driver mutations undergo selection. Metastases had genomic alterations distinct from their primary tumors; metastatic dissemination sometimes occurred early during the development of the primary tumor. Our study offers new insights into the genetics and evolution of this melanoma subtype, providing potential biomarkers for progression and therapy.