To explore the distinct genotypic and phenotypic states of melanoma tumors we applied single-cell RNA-seq to 4,645 single cells isolated from 19 patients, profiling malignant, immune, stromal and endothelial cells. Malignant cells within the same tumor displayed transcriptional heterogeneity associated with the cell cycle, spatial context, and a drug resistance program. In particular, all tumors harbored malignant cells from two distinct transcriptional cell states, such that “MITF-high” tumors also contained “AXL-high” tumor cells. Single-cell analyses suggested distinct tumor micro-environmental patterns, including cell-to-cell interactions. Analysis of tumor-infiltrating T cells revealed exhaustion programs, their connection to T cell activation and to clonal expansion, and their variability across patients. Overall, we begin to unravel the cellular ecosystem of tumors and how single cell genomics offers insights with implications for both targeted and immune therapies.

Most patients with BRAFV600 metastatic melanoma develop resistance to selective RAF kinase inhibitors. The spectrum of clinical genetic resistance mechanisms to RAF inhibitors and options for salvage therapy are incompletely understood. We performed whole exome sequencing on formalin-fixed, paraffin embedded (FFPE) tumors from 45 patients with BRAFV600 metastatic melanoma who received vemurafenib or dabrafenib monotherapy. Genetic alterations in known or putative RAF inhibitor resistance genes were observed in 23 of 45 patients (51%). Besides previously characterized alterations, we discovered a “long tail” of new MAPK pathway alterations (MAP2K2, MITF) that confer RAF inhibitor resistance. In three cases, multiple resistance gene alterations were observed within the same tumor biopsy. Overall, RAF inhibitor therapy leads to diverse clinical genetic resistance mechanisms, mostly involving MAPK pathway reactivation. Novel therapeutic combinations may be needed to achieve durable clinical control of BRAFV600 melanoma. Integrating clinical genomics with preclinical screens may model subsequent resistance studies.