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.

SUMMARY Immune checkpoint inhibitors (ICIs) produce durable responses in some melanoma patients, but many patients derive no clinical benefit, and the molecular underpinnings of such resistance remain elusive. Here, we leveraged single-cell RNA sequencing (scRNA-seq) from 33 melanoma tumors and computational analyses to interrogate malignant cell states that promote immune evasion. We identified a resistance program expressed by malignant cells that is associated with T cell exclusion and immune evasion. The program is expressed prior to immunotherapy, characterizes cold niches in situ, and predicts clinical responses to anti-PD-1 therapy in an independent cohort of 112 melanoma patients. CDK4/6-inhibition represses this program in individual malignant cells, induces senescence, and reduces melanoma tumor outgrowth in mouse models in vivo when given in combination with immunotherapy. Our study provides a high-resolution landscape of ICI-resistant cell states, identifies clinically predictive signatures, and suggests new therapeutic strategies to overcome immunotherapy resistance.

systems that incorporate features of the tumor microenvironment and model the dynamic response to immune checkpoint blockade (ICB) may facilitate efforts in precision immuno-oncology and the development of effective combination therapies. Here, we demonstrate the ability to interrogate response to ICB using murine- and patient-derived organotypic tumor spheroids (MDOTS/PDOTS). MDOTS/PDOTS isolated from mouse and human tumors retain autologous lymphoid and myeloid cell populations and respond to ICB in short-term three-dimensional microfluidic culture. Response and resistance to ICB was recapitulated using MDOTS derived from established immunocompetent mouse tumor models. MDOTS profiling demonstrated that TBK1/IKKε inhibition enhanced response to PD-1 blockade, which effectively predicted tumor response Systematic profiling of secreted cytokines in PDOTS captured key features associated with response and resistance to PD-1 blockade. Thus, MDOTS/PDOTS profiling represents a novel platform to evaluate ICB using established murine models as well as clinically relevant patient specimens. Resistance to PD-1 blockade remains a challenge for many patients, and biomarkers to guide treatment are lacking. Here, we demonstrate feasibility of profiling of PD-1 blockade to interrogate the tumor immune microenvironment, develop therapeutic combinations, and facilitate precision immuno-oncology efforts..

The architecture of normal and diseased tissues strongly influences the development and progression of disease as well as responsiveness and resistance to therapy. We describe a tissue-based cyclic immunofluorescence (t-CyCIF) method for highly multiplexed immuno-fluorescence imaging of formalin-fixed, paraffin-embedded (FFPE) specimens mounted on glass slides, the most widely used specimens for histopathological diagnosis of cancer and other diseases. t-CyCIF generates up to 60-plex images using an iterative process (a cycle) in which conventional low-plex fluorescence images are repeatedly collected from the same sample and then assembled into a high-dimensional representation. t-CyCIF requires no specialized instruments or reagents and is compatible with super-resolution imaging; we demonstrate its application to quantifying signal transduction cascades, tumor antigens and immune markers in diverse tissues and tumors. The simplicity and adaptability of t-CyCIF makes it an effective method for pre-clinical and clinical research and a natural complement to single-cell genomics.