MAESTRO-Pool Enables Highly Parallel and Specific Mutation-Enrichment Sequencing for Minimal Residual Disease Detection in Cohort Studies.

BACKGROUND: Tracing patient-specific tumor mutations in cell-free DNA (cfDNA) for minimal residual disease (MRD) detection is promising but challenging. Assaying more mutations and cfDNA stands to improve MRD detection but requires highly accurate, efficient sequencing methods and proper calibration to prevent false detection with bespoke tests. METHODS: MAESTRO (Minor Allele Enriched Sequencing Through Recognition Oligonucleotides) uses mutation-specific oligonucleotide probes to enrich cfDNA libraries for tumor mutations and enable their accurate detection with minimal sequencing. A new approach, MAESTRO-Pool, which entails pooling MAESTRO probes for all patients and applying these to all samples from all patients, was used to screen for 22 333 tumor mutations from 9 melanoma patients in 98 plasma samples. This enabled quantification of MRD detection in patient-matched samples and false detection in unmatched samples from other patients. To detect MRD, a new dynamic MRD caller was used that computes a probability for MRD detection based on the number of mutations and cfDNA molecules sequenced, thereby calibrating for variations in each bespoke test. RESULTS: MAESTRO-Pool enabled sensitive detection of MRD down to 0.78 parts per million (ppm), reflecting a 10- to 100-fold improvement over existing tests. Of the 8 MRD positive samples with ultra-low tumor fractions <10 ppm, 7 were either in upward-trend preceding recurrence or downward-trend aligning with response. Of 784 patient-unmatched tests, only one was found as MRD positive (tumor fraction = 2.7 ppm), suggesting high specificity. CONCLUSIONS: MAESTRO-Pool enables massively parallel, tumor-informed MRD testing with concurrent benchmarking of bespoke MRD tests. Meanwhile, our new MRD caller enables more mutations and cfDNA molecules to be tested without compromising specificity. These improve the ability for detecting traces of MRD from blood.

Targeting TBK1 to overcome resistance to cancer immunotherapy.

Despite the success of PD-1 blockade in melanoma and other cancers, effective treatment strategies to overcome resistance to cancer immunotherapy are lacking1,2. Here we identify the innate immune kinase TANK-binding kinase 1 (TBK1)3 as a candidate immune-evasion gene in a pooled genetic screen4. Using a suite of genetic and pharmacological tools across multiple experimental model systems, we confirm a role for TBK1 as an immune-evasion gene. Targeting TBK1 enhances responses to PD-1 blockade by decreasing the cytotoxicity threshold to effector cytokines (TNF and IFNγ). TBK1 inhibition in combination with PD-1 blockade also demonstrated efficacy using patient-derived tumour models, with concordant findings in matched patient-derived organotypic tumour spheroids and matched patient-derived organoids. Tumour cells lacking TBK1 are primed to undergo RIPK- and caspase-dependent cell death in response to TNF and IFNγ in a JAK-STAT-dependent manner. Taken together, our results demonstrate that targeting TBK1 is an effective strategy to overcome resistance to cancer immunotherapy.

Fucosylation of HLA-DRB1 regulates CD4+ T cell-mediated anti-melanoma immunity and enhances immunotherapy efficacy.

Immunotherapy efficacy is limited in melanoma, and combinations of immunotherapies with other modalities have yielded limited improvements but also adverse events requiring cessation of treatment. In addition to ineffective patient stratification, efficacy is impaired by paucity of intratumoral immune cells (itICs); thus, effective strategies to safely increase itICs are needed. We report that dietary administration of L-fucose induces fucosylation and cell surface enrichment of the major histocompatibility complex (MHC)-II protein HLA-DRB1 in melanoma cells, triggering CD4+ T cell-mediated increases in itICs and anti-tumor immunity, enhancing immune checkpoint blockade responses. Melanoma fucosylation and fucosylated HLA-DRB1 associate with intratumoral T cell abundance and anti-programmed cell death protein 1 (PD1) responder status in patient melanoma specimens, suggesting the potential use of melanoma fucosylation as a strategy for stratifying patients for immunotherapies. Our findings demonstrate that fucosylation is a key mediator of anti-tumor immunity and, importantly, suggest that L-fucose is a powerful agent for safely increasing itICs and immunotherapy efficacy in melanoma.

Myeloid-derived itaconate suppresses cytotoxic CD8+ T cells and promotes tumour growth.

The tumour microenvironment possesses mechanisms that suppress anti-tumour immunity. Itaconate is a metabolite produced from the Krebs cycle intermediate cis-aconitate by the activity of immune-responsive gene 1 (IRG1). While it is known to be immune modulatory, the role of itaconate in anti-tumour immunity is unclear. Here, we demonstrate that myeloid-derived suppressor cells (MDSCs) secrete itaconate that can be taken up by CD8+ T cells and suppress their proliferation, cytokine production and cytolytic activity. Metabolite profiling, stable-isotope tracing and metabolite supplementation studies indicated that itaconate suppressed the biosynthesis of aspartate and serine/glycine in CD8+ T cells to attenuate their proliferation and function. Host deletion of Irg1 in female mice bearing allografted tumours resulted in decreased tumour growth, inhibited the immune-suppressive activities of MDSCs, promoted anti-tumour immunity of CD8+ T cells and enhanced the anti-tumour activity of anti-PD-1 antibody treatment. Furthermore, we found a significant negative correlation between IRG1 expression and response to PD-1 immune checkpoint blockade in patients with melanoma. Our findings not only reveal a previously unknown role of itaconate as an immune checkpoint metabolite secreted from MDSCs to suppress CD8+ T cells, but also establish IRG1 as a myeloid-selective target in immunometabolism whose inhibition promotes anti-tumour immunity and enhances the efficacy of immune checkpoint protein blockade.

In vivo CRISPR screens reveal the landscape of immune evasion pathways across cancer.

The immune system can eliminate tumors, but checkpoints enable immune escape. Here, we identify immune evasion mechanisms using genome-scale in vivo CRISPR screens across cancer models treated with immune checkpoint blockade (ICB). We identify immune evasion genes and important immune inhibitory checkpoints conserved across cancers, including the non-classical major histocompatibility complex class I (MHC class I) molecule Qa-1b/HLA-E. Surprisingly, loss of tumor interferon-γ (IFNγ) signaling sensitizes many models to immunity. The immune inhibitory effects of tumor IFN sensing are mediated through two mechanisms. First, tumor upregulation of classical MHC class I inhibits natural killer cells. Second, IFN-induced expression of Qa-1b inhibits CD8+ T cells via the NKG2A/CD94 receptor, which is induced by ICB. Finally, we show that strong IFN signatures are associated with poor response to ICB in individuals with renal cell carcinoma or melanoma. This study reveals that IFN-mediated upregulation of classical and non-classical MHC class I inhibitory checkpoints can facilitate immune escape.

Proteasomal down-regulation of the proapoptotic MST2 pathway contributes to BRAF inhibitor resistance in melanoma.

The RAS-RAF-MEK-ERK pathway is hyperactivated in most malignant melanomas, and mutations in BRAF or NRAS account for most of these cases. BRAF inhibitors (BRAFi) are highly efficient for treating patients with BRAFV600E mutations, but tumours frequently acquire resistance within a few months. Multiple resistance mechanisms have been identified, due to mutations or network adaptations that revive ERK signalling. We have previously shown that RAF proteins inhibit the MST2 proapoptotic pathway in a kinase-independent fashion. Here, we have investigated the role of the MST2 pathway in mediating resistance to BRAFi. We show that the BRAFV600E mutant protein, but not the wild-type BRAF protein, binds to MST2 inhibiting its proapoptotic signalling. Down-regulation of MST2 reduces BRAFi-induced apoptosis. In BRAFi-resistant cell lines, MST2 pathway proteins are down-regulated by ubiquitination and subsequent proteasomal degradation rendering cells refractory to MST2 pathway-induced apoptosis. Restoration of apoptosis can be achieved by increasing MST2 pathway protein expression using proteasome inhibitors. In summary, we show that the MST2 pathway plays a role in the acquisition of BRAFi resistance in melanoma.

HRS phosphorylation drives immunosuppressive exosome secretion and restricts CD8+ T-cell infiltration into tumors.

The lack of tumor infiltration by CD8+ T cells is associated with poor patient response to anti-PD-1 therapy. Understanding how tumor infiltration is regulated is key to improving treatment efficacy. Here, we report that phosphorylation of HRS, a pivotal component of the ESCRT complex involved in exosome biogenesis, restricts tumor infiltration of cytolytic CD8+ T cells. Following ERK-mediated phosphorylation, HRS interacts with and mediates the selective loading of PD-L1 to exosomes, which inhibits the migration of CD8+ T cells into tumors. In tissue samples from patients with melanoma, CD8+ T cells are excluded from the regions where tumor cells contain high levels of phosphorylated HRS. In murine tumor models, overexpression of phosphorylated HRS increases resistance to anti-PD-1 treatment, whereas inhibition of HRS phosphorylation enhances treatment efficacy. Our study reveals a mechanism by which phosphorylation of HRS in tumor cells regulates anti-tumor immunity by inducing PD-L1+ immunosuppressive exosomes, and suggests HRS phosphorylation blockade as a potential strategy to improve the efficacy of cancer immunotherapy.

Microenvironmental Landscape of Human Melanoma Brain Metastases in Response to Immune Checkpoint Inhibition.

Melanoma-derived brain metastases (MBM) represent an unmet clinical need because central nervous system progression is frequently an end stage of the disease. Immune checkpoint inhibitors (ICI) provide a clinical opportunity against MBM; however, the MBM tumor microenvironment (TME) has not been fully elucidated in the context of ICI. To dissect unique elements of the MBM TME and correlates of MBM response to ICI, we collected 32 fresh MBM and performed single-cell RNA sequencing of the MBM TME and T-cell receptor clonotyping on T cells from MBM and matched blood and extracranial lesions. We observed myeloid phenotypic heterogeneity in the MBM TME, most notably multiple distinct neutrophil states, including an IL8-expressing population that correlated with malignant cell epithelial-to-mesenchymal transition. In addition, we observed significant relationships between intracranial T-cell phenotypes and the distribution of T-cell clonotypes intracranially and peripherally. We found that the phenotype, clonotype, and overall number of MBM-infiltrating T cells were associated with response to ICI, suggesting that ICI-responsive MBMs interact with peripheral blood in a manner similar to extracranial lesions. These data identify unique features of the MBM TME that may represent potential targets to improve clinical outcomes for patients with MBM.

Landscape of helper and regulatory antitumour CD4+ T cells in melanoma.

Within the tumour microenvironment, CD4+ T cells can promote or suppress antitumour responses through the recognition of antigens presented by human leukocyte antigen (HLA) class II molecules1,2, but how cancers co-opt these physiologic processes to achieve immune evasion remains incompletely understood. Here we performed in-depth analysis of the phenotype and tumour specificity of CD4+ T cells infiltrating human melanoma specimens, finding that exhausted cytotoxic CD4+ T cells could be directly induced by melanoma cells through recognition of HLA class II-restricted neoantigens, and also HLA class I-restricted tumour-associated antigens. CD4+ T regulatory (TReg) cells could be indirectly elicited through presentation of tumour antigens via antigen-presenting cells. Notably, numerous tumour-reactive CD4+ TReg clones were stimulated directly by HLA class II-positive melanoma and demonstrated specificity for melanoma neoantigens. This phenomenon was observed in the presence of an extremely high tumour neoantigen load, which we confirmed to be associated with HLA class II positivity through the analysis of 116 melanoma specimens. Our data reveal the landscape of infiltrating CD4+ T cells in melanoma and point to the presentation of HLA class II-restricted neoantigens and direct engagement of immunosuppressive CD4+ TReg cells as a mechanism of immune evasion that is favoured in HLA class II-positive melanoma.

STAG2 regulates interferon signaling in melanoma via enhancer loop reprogramming.

The cohesin complex participates in the organization of 3D genome through generating and maintaining DNA loops. Stromal antigen 2 (STAG2), a core subunit of the cohesin complex, is frequently mutated in various cancers. However, the impact of STAG2 inactivation on 3D genome organization, especially the long-range enhancer-promoter contacts and subsequent gene expression control in cancer, remains poorly understood. Here we show that depletion of STAG2 in melanoma cells leads to expansion of topologically associating domains (TADs) and enhances the formation of acetylated histone H3 lysine 27 (H3K27ac)-associated DNA loops at sites where binding of STAG2 is switched to its paralog STAG1. We further identify Interferon Regulatory Factor 9 (IRF9) as a major direct target of STAG2 in melanoma cells via integrated RNA-seq, STAG2 ChIP-seq and H3K27ac HiChIP analyses. We demonstrate that loss of STAG2 activates IRF9 through modulating the 3D genome organization, which in turn enhances type I interferon signaling and increases the expression of PD-L1. Our findings not only establish a previously unknown role of the STAG2 to STAG1 switch in 3D genome organization, but also reveal a functional link between STAG2 and interferon signaling in cancer cells, which may enhance the immune evasion potential in STAG2-mutant cancer.

Combined tumor and immune signals from genomes or transcriptomes predict outcomes of checkpoint inhibition in melanoma.

Immune checkpoint blockade (CPB) improves melanoma outcomes, but many patients still do not respond. Tumor mutational burden (TMB) and tumor-infiltrating T cells are associated with response, and integrative models improve survival prediction. However, integrating immune/tumor-intrinsic features using data from a single assay (DNA/RNA) remains underexplored. Here, we analyze whole-exome and bulk RNA sequencing of tumors from new and published cohorts of 189 and 178 patients with melanoma receiving CPB, respectively. Using DNA, we calculate T cell and B cell burdens (TCB/BCB) from rearranged TCR/Ig sequences and find that patients with TMBhigh and TCBhigh or BCBhigh have improved outcomes compared to other patients. By combining pairs of immune- and tumor-expressed genes, we identify three gene pairs associated with response and survival, which validate in independent cohorts. The top model includes lymphocyte-expressed MAP4K1 and tumor-expressed TBX3. Overall, RNA or DNA-based models combining immune and tumor measures improve predictions of melanoma CPB outcomes.

Benefit and toxicity of programmed death-1 blockade vary by ethnicity in patients with advanced melanoma: an international multicentre observational study.

BACKGROUND: Programmed cell death receptor-1 (PD-1) monotherapy is a standard treatment for advanced cutaneous melanoma, but its efficacy and toxicity are defined in white populations and remain poorly characterized in other ethnic groups, such as East Asian, Hispanic and African. OBJECTIVES: To determine the efficacy and toxicity of PD-1 monotherapy in different ethnic groups. METHODS: Clinical data for patients with unresectable or advanced melanoma treated with anti-PD-1 monotherapy between 2009 and 2019 were collected retrospectively from five independent institutions in the USA, Australia and China. Tumour response, survival and immune-related adverse events (irAEs) were compared by ethnicity (white vs. East Asian/Hispanic/African) across different melanoma subtypes: nonacral cutaneous (NAC)/unknown primary (UP) and acral/mucosal/uveal. RESULTS: In total, 1135 patients were included. White patients had significantly higher objective response rate (ORR) [54%, 95% confidence interval (CI) 50-57% vs. 20%, 95% CI 13-28%; adjusted P < 0·001] and longer progression-free survival (14·2 months, 95% CI 10·7-20·3 vs. 5·4 months, 95% CI 4·5-7·0; adjusted P < 0·001) than East Asian, Hispanic and African patients in the NAC and UP subtypes. White ethnicity remained independently associated with a higher ORR (odds ratio 4·10, 95% CI 2·48-6·81; adjusted P < 0·001) and longer PFS (hazard ratio 0·58, 95% CI 0·46-0·74; adjusted P < 0·001) in multivariate analyses after adjustment for age, sex, primary anatomical location, metastasis stage, baseline lactate dehydrogenase level, mutational status and prior systemic treatment. White and East Asian/Hispanic/African patients shared similar ORR and progression-free survival in acral/mucosal/uveal melanomas. Similar melanoma-subtype-specific ethnic discrepancies were observed in complete response rate and overall survival. White patients had higher rates of gastrointestinal irAEs but lower rates of endocrine, liver and other rare types of irAEs. These differences in irAEs by ethnicity were not attributable to varying melanoma subtypes. CONCLUSIONS: Ethnic discrepancy in clinical benefit is specific to melanoma subtype, and East Asian, Hispanic and African patients with NAC and UP melanomas have poorer clinical benefits than previously recognized. The ethnic discrepancy in toxicity observed across different melanoma subtypes warrants an ethnicity-based irAE surveillance strategy. More research is needed to elucidate the molecular and immunological determinants of these differences. What is already known about this topic? There is a great difference in response to immunotherapy between different subtypes of melanoma (cutaneous, mucosal, acral and uveal) in patients with advanced disease. What does this study add? Our data show for the first time that there are differences between different ethnic groups in terms of both response and toxicity to immunotherapy beyond the well-appreciated discrepancies due to melanoma subtype.

Pathway signatures derived from on-treatment tumor specimens predict response to anti-PD1 blockade in metastatic melanoma.

Both genomic and transcriptomic signatures have been developed to predict responses of metastatic melanoma to immune checkpoint blockade (ICB) therapies; however, most of these signatures are derived from pre-treatment biopsy samples. Here, we build pathway-based super signatures in pre-treatment (PASS-PRE) and on-treatment (PASS-ON) tumor specimens based on transcriptomic data and clinical information from a large dataset of metastatic melanoma treated with anti-PD1-based therapies as the training set. Both PASS-PRE and PASS-ON signatures are validated in three independent datasets of metastatic melanoma as the validation set, achieving area under the curve (AUC) values of 0.45-0.69 and 0.85-0.89, respectively. We also combine all test samples and obtain AUCs of 0.65 and 0.88 for PASS-PRE and PASS-ON signatures, respectively. When compared with existing signatures, the PASS-ON signature demonstrates more robust and superior predictive performance across all four datasets. Overall, we provide a framework for building pathway-based signatures that is highly and accurately predictive of response to anti-PD1 therapies based on on-treatment tumor specimens. This work would provide a rationale for applying pathway-based signatures derived from on-treatment tumor samples to predict patients' therapeutic response to ICB therapies.

Absolute quantification of tumor antigens using embedded MHC-I isotopologue calibrants.

Absolute quantification measurements (copies per cell) of peptide major histocompatibility complex (pMHC) antigens are necessary to inform targeted immunotherapy drug design; however, existing methods for absolute quantification have critical limitations. Here, we present a platform termed SureQuant-IsoMHC, utilizing a series of pMHC isotopologues and internal standard-triggered targeted mass spectrometry to generate an embedded multipoint calibration curve to determine endogenous pMHC concentrations for a panel of 18 tumor antigens. We apply SureQuant-IsoMHC to measure changes in expression of our target panel in a melanoma cell line treated with a MEK inhibitor and translate this approach to estimate antigen concentrations in melanoma tumor biopsies.

Early Use of High-Dose Glucocorticoid for the Management of irAE Is Associated with Poorer Survival in Patients with Advanced Melanoma Treated with Anti-PD-1 Monotherapy.

PURPOSE: Programmed cell death receptor-1 (PD-1) inhibitors are frontline therapy in advanced melanoma. Severe immune-related adverse effects (irAEs) often require immunosuppressive treatment with glucocorticoids (GCCs), but GCC use and its correlation with patient survival outcomes during anti-PD-1 monotherapy remains unclear. EXPERIMENTAL DESIGN: In this multicenter retrospective analysis, patients treated with anti-PD-1 monotherapy between 2009 and 2019 and detailed GCC use, data were identified from five independent cohorts, with median follow-up time of 206 weeks. IrAEs were tracked from the initiation of anti-PD-1 until disease progression, initiation of a new therapy, or last follow-up. Correlations between irAEs, GCC use, and survival outcomes were analyzed. RESULTS: Of the entire cohort of 947 patients, 509 (54%) developed irAEs. In the MGH cohort [irAE(+) n = 90], early-onset irAE (within 8 weeks of anti-PD-1 initiation) with high-dose GCC use (≥60-mg prednisone equivalent once a day) was independently associated with poorer post-irAE PFS/OS (progression-free survival/overall survival) [post-irAE PFS: HR, 5.37; 95% confidence interval (CI), 2.10-13.70; P < 0.001; post-irAE OS: HR, 5.95; 95% CI, 2.20-16.09; P < 0.001] compared with irAEs without early high-dose GCC use. These findings were validated in the combined validation cohort [irAE(+) n = 419, post-irAE PFS: HR, 1.69; 95% CI, 1.04-2.76; P = 0.04; post-irAE OS: HR, 1.97; 95% CI, 1.15-3.39; P = 0.01]. Similar findings were also observed in the 26-week landmark analysis for post-irAE-PFS but not for post-irAE-OS. A sensitivity analysis using accumulated GCC exposure as the measurement achieved similar results. CONCLUSIONS: Early high-dose GCC use was associated with poorer PFS and OS after irAE onset. Judicious use of GCC early during anti-PD-1 monotherapy should be considered. Further prospective randomized control clinical trials designed to explore alternative irAE management options are warranted.

Spatially organized multicellular immune hubs in human colorectal cancer.

Immune responses to cancer are highly variable, with mismatch repair-deficient (MMRd) tumors exhibiting more anti-tumor immunity than mismatch repair-proficient (MMRp) tumors. To understand the rules governing these varied responses, we transcriptionally profiled 371,223 cells from colorectal tumors and adjacent normal tissues of 28 MMRp and 34 MMRd individuals. Analysis of 88 cell subsets and their 204 associated gene expression programs revealed extensive transcriptional and spatial remodeling across tumors. To discover hubs of interacting malignant and immune cells, we identified expression programs in different cell types that co-varied across tumors from affected individuals and used spatial profiling to localize coordinated programs. We discovered a myeloid cell-attracting hub at the tumor-luminal interface associated with tissue damage and an MMRd-enriched immune hub within the tumor, with activated T cells together with malignant and myeloid cells expressing T cell-attracting chemokines. By identifying interacting cellular programs, we reveal the logic underlying spatially organized immune-malignant cell networks.

Neural Crest-Like Stem Cell Transcriptome Analysis Identifies LPAR1 in Melanoma Progression and Therapy Resistance.

Metastatic melanoma is challenging to clinically address. Although standard-of-care targeted therapy has high response rates in patients with BRAF-mutant melanoma, therapy relapse occurs in most cases. Intrinsically resistant melanoma cells drive therapy resistance and display molecular and biologic properties akin to neural crest-like stem cells (NCLSC) including high invasiveness, plasticity, and self-renewal capacity. The shared transcriptional programs and vulnerabilities between NCLSCs and cancer cells remains poorly understood. Here, we identify a developmental LPAR1-axis critical for NCLSC viability and melanoma cell survival. LPAR1 activity increased during progression and following acquisition of therapeutic resistance. Notably, genetic inhibition of LPAR1 potentiated BRAFi ± MEKi efficacy and ablated melanoma migration and invasion. Our data define LPAR1 as a new therapeutic target in melanoma and highlights the promise of dissecting stem cell-like pathways hijacked by tumor cells. SIGNIFICANCE: This study identifies an LPAR1-axis critical for melanoma invasion and intrinsic/acquired therapy resistance.

Phenotype, specificity and avidity of antitumour CD8+ T cells in melanoma.

Interactions between T cell receptors (TCRs) and their cognate tumour antigens are central to antitumour immune responses1-3; however, the relationship between phenotypic characteristics and TCR properties is not well elucidated. Here we show, by linking the antigenic specificity of TCRs and the cellular phenotype of melanoma-infiltrating lymphocytes at single-cell resolution, that tumour specificity shapes the expression state of intratumoural CD8+ T cells. Non-tumour-reactive T cells were enriched for viral specificities and exhibited a non-exhausted memory phenotype, whereas melanoma-reactive lymphocytes predominantly displayed an exhausted state that encompassed diverse levels of differentiation but rarely acquired memory properties. These exhausted phenotypes were observed both among clonotypes specific for public overexpressed melanoma antigens (shared across different tumours) or personal neoantigens (specific for each tumour). The recognition of such tumour antigens was provided by TCRs with avidities inversely related to the abundance of cognate targets in melanoma cells and proportional to the binding affinity of peptide-human leukocyte antigen (HLA) complexes. The persistence of TCR clonotypes in peripheral blood was negatively affected by the level of intratumoural exhaustion, and increased in patients with a poor response to immune checkpoint blockade, consistent with chronic stimulation mediated by residual tumour antigens. By revealing how the quality and quantity of tumour antigens drive the features of T cell responses within the tumour microenvironment, we gain insights into the properties of the anti-melanoma TCR repertoire.

Evolution of delayed resistance to immunotherapy in a melanoma responder.

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