Targeted BRAF inhibition (BRAFi) and combined BRAF and MEK inhibition (BRAFi and MEKi) therapies have markedly improved the clinical outcomes of patients with metastatic melanoma. Unfortunately, the efficacy of these treatments is often countered by the acquisition of drug resistance. Here we investigated the molecular mechanisms that underlie acquired resistance to BRAFi and to the combined therapy. Consistent with previous studies, we show that resistance to BRAFi is mediated by ERK pathway reactivation. Resistance to the combined therapy, however, is mediated by mechanisms independent of reactivation of ERK in many resistant cell lines and clinical samples. p21-activated kinases (PAKs) become activated in cells with acquired drug resistance and have a pivotal role in mediating resistance. Our screening, using a reverse-phase protein array, revealed distinct mechanisms by which PAKs mediate resistance to BRAFi and the combined therapy. In BRAFi-resistant cells, PAKs phosphorylate CRAF and MEK to reactivate ERK. In cells that are resistant to the combined therapy, PAKs regulate JNK and β-catenin phosphorylation and mTOR pathway activation, and inhibit apoptosis, thereby bypassing ERK. Together, our results provide insights into the molecular mechanisms underlying acquired drug resistance to current targeted therapies, and may help to direct novel drug development efforts to overcome acquired drug resistance.

Diverse pathways drive resistance to BRAF/MEK inhibitors in BRAF-mutant melanoma, suggesting that durable control of resistance will be a challenge. By combining statistical modeling of genomic data from matched pre-treatment and post-relapse patient tumors with functional interrogation of >20 in vitro and in vivo resistance models, we discovered that major pathways of resistance converge to activate the transcription factor, c-MYC (MYC). MYC expression and pathway gene signatures were suppressed following drug treatment, and then rebounded during progression. Critically, MYC activation was necessary and sufficient for resistance, and suppression of MYC activity using genetic approaches or BET bromodomain inhibition was sufficient to resensitize cells and delay BRAFi resistance. Finally, MYC-driven, BRAFi-resistant cells are hypersensitive to the inhibition of MYC synthetic lethal partners, including SRC family and c-KIT tyrosine kinases, as well as glucose, glutamine, and serine metabolic pathways. These insights enable the design of combination therapies that select against resistance evolution.

Treatment of cancer has been revolutionized by immune checkpoint blockade therapies. Despite the high rate of response in advanced melanoma, the majority of patients succumb to disease. To identify factors associated with success or failure of checkpoint therapy, we profiled transcriptomes of 16,291 individual immune cells from 48 tumor samples of melanoma patients treated with checkpoint inhibitors. Two distinct states of CD8+ T cells were defined by clustering and associated with patient tumor regression or progression. A single transcription factor, TCF7, was visualized within CD8+ T cells in fixed tumor samples and predicted positive clinical outcome in an independent cohort of checkpoint-treated patients. We delineated the epigenetic landscape and clonality of these T cell states and demonstrated enhanced antitumor immunity by targeting novel combinations of factors in exhausted cells. Our study of immune cell transcriptomes from tumors demonstrates a strategy for identifying predictors, mechanisms, and targets for enhancing checkpoint immunotherapy.

Natural mitochondrial DNA (mtDNA) mutations enable the inference of clonal relationships among cells. mtDNA can be profiled along with measures of cell state, but has not yet been combined with the massively parallel approaches needed to tackle the complexity of human tissue. Here, we introduce a high-throughput, droplet-based mitochondrial single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq), a method that combines high-confidence mtDNA mutation calling in thousands of single cells with their concomitant high-quality accessible chromatin profile. This enables the inference of mtDNA heteroplasmy, clonal relationships, cell state and accessible chromatin variation in individual cells. We reveal single-cell variation in heteroplasmy of a pathologic mtDNA variant, which we associate with intra-individual chromatin variability and clonal evolution. We clonally trace thousands of cells from cancers, linking epigenomic variability to subclonal evolution, and infer cellular dynamics of differentiating hematopoietic cells in vitro and in vivo. Taken together, our approach enables the study of cellular population dynamics and clonal properties in vivo.

Melanomas can switch to a dedifferentiated cell state upon exposure to cytotoxic T cells. However, it is unclear whether such tumor cells pre-exist in patients and whether they can be resensitized to immunotherapy. Here, we chronically expose (patient-derived) melanoma cell lines to differentiation antigen-specific cytotoxic T cells and observe strong enrichment of a pre-existing NGFRhi population. These fractions are refractory also to T cells recognizing non-differentiation antigens, as well as to BRAF + MEK inhibitors. NGFRhi cells induce the neurotrophic factor BDNF, which contributes to T cell resistance, as does NGFR. In melanoma patients, a tumor-intrinsic NGFR signature predicts anti-PD-1 therapy resistance, and NGFRhi tumor fractions are associated with immune exclusion. Lastly, pharmacologic NGFR inhibition restores tumor sensitivity to T cell attack in vitro and in melanoma xenografts. These findings demonstrate the existence of a stable and pre-existing NGFRhi multitherapy-refractory melanoma subpopulation, which ought to be eliminated to revert intrinsic resistance to immunotherapeutic intervention.

Immune checkpoint inhibitors (ICIs) show promise, but most patients do not respond. We identify and validate biomarkers from extracellular vesicles (EVs), allowing non-invasive monitoring of tumor- intrinsic and host immune status, as well as a prediction of ICI response. We undertook transcriptomic profiling of plasma-derived EVs and tumors from 50 patients with metastatic melanoma receiving ICI, and validated with an independent EV-only cohort of 30 patients. Plasma-derived EV and tumor transcriptomes correlate. EV profiles reveal drivers of ICI resistance and melanoma progression, exhibit differentially expressed genes/pathways, and correlate with clinical response to ICI. We created a Bayesian probabilistic deconvolution model to estimate contributions from tumor and non-tumor sources, enabling interpretation of differentially expressed genes/pathways. EV RNA-seq mutations also segregated ICI response. EVs serve as a non-invasive biomarker to jointly probe tumor-intrinsic and immune changes to ICI, function as predictive markers of ICI responsiveness, and monitor tumor persistence and immune activation.

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.

Immune checkpoint blockade (ICB) has demonstrated efficacy in patients with melanoma, but many exhibit poor responses. Using single cell RNA sequencing of melanoma patient-derived circulating tumor cells (CTCs) and functional characterization using mouse melanoma models, we show that the KEAP1/NRF2 pathway modulates sensitivity to ICB, independently of tumorigenesis. The NRF2 negative regulator, KEAP1, shows intrinsic variation in expression, leading to tumor heterogeneity and subclonal resistance.

The poor efficacy of chimeric antigen receptor T-cell therapy (CAR T) for solid tumors is due to insufficient CAR T cell tumor infiltration, in vivo expansion, persistence, and effector function, as well as exhaustion, intrinsic target antigen heterogeneity or antigen loss of target cancer cells, and immunosuppressive tumor microenvironment (TME). Here we describe a broadly applicable nongenetic approach that simultaneously addresses the multiple challenges of CAR T as a therapy for solid tumors. The approach reprograms CAR T cells by exposing them to stressed target cancer cells which have been exposed to the cell stress inducer disulfiram (DSF) and copper (Cu)(DSF/Cu) plus ionizing irradiation (IR). The reprogrammed CAR T cells acquire early memory-like characteristics, potent cytotoxicity, enhanced in vivo expansion, persistence, and decreased exhaustion. Tumors stressed by DSF/Cu and IR also reprogram and reverse the immunosuppressive TME in humanized mice. The reprogrammed CAR T cells, derived from peripheral blood mononuclear cells of healthy donors or metastatic female breast cancer patients, induce robust, sustained memory and curative anti-solid tumor responses in multiple xenograft mouse models, establishing proof of concept for empowering CAR T by stressing tumor as a promising therapy for solid tumors.

Lineage tracing provides key insights into the fate of individual cells in complex organisms. Although effective genetic labeling approaches are available in model systems, in humans, most approaches require detection of nuclear somatic mutations, which have high error rates, limited scale, and do not capture cell state information. Here, we show that somatic mutations in mtDNA can be tracked by single-cell RNA or assay for transposase accessible chromatin (ATAC) sequencing. We leverage somatic mtDNA mutations as natural genetic barcodes and demonstrate their utility as highly accurate clonal markers to infer cellular relationships. We track native human cells both in vitro and in vivo and relate clonal dynamics to gene expression and chromatin accessibility. Our approach should allow clonal tracking at a 1,000-fold greater scale than with nuclear genome sequencing, with simultaneous information on cell state, opening the way to chart cellular dynamics in human health and disease.

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.

BRAF is a serine/threonine kinase that harbors activating mutations in ∼7% of human malignancies and ∼60% of melanomas. Despite initial clinical responses to BRAF inhibitors, patients frequently develop drug resistance. To identify candidate therapeutic targets for BRAF inhibitor resistant melanoma, we conduct CRISPR screens in melanoma cells harboring an activating BRAF mutation that had also acquired resistance to BRAF inhibitors. To investigate the mechanisms and pathways enabling resistance to BRAF inhibitors in melanomas, we integrate expression, ATAC-seq, and CRISPR screen data. We identify the JUN family transcription factors and the ETS family transcription factor ETV5 as key regulators of CDK6, which together enable resistance to BRAF inhibitors in melanoma cells. Our findings reveal genes contributing to resistance to a selective BRAF inhibitor PLX4720, providing new insights into gene regulation in BRAF inhibitor resistant melanoma cells.

We determined the frequency of recurrent hotspot mutations in 46 cancer-related genes across tumor histologies in patients with advanced cancer.

Checkpoint blockade with antibodies specific for the PD-1 and CTLA-4 inhibitory receptors can induce durable responses in a wide range of human cancers. However, the immunological mechanisms responsible for severe inflammatory side effects remain poorly understood. Here we report a comprehensive single-cell analysis of immune cell populations in colitis, a common and severe side effect of checkpoint blockade. We observed a striking accumulation of CD8 T cells with highly cytotoxic and proliferative states and no evidence of regulatory T cell depletion. T cell receptor (TCR) sequence analysis demonstrated that a substantial fraction of colitis-associated CD8 T cells originated from tissue-resident populations, explaining the frequently early onset of colitis symptoms following treatment initiation. Our analysis also identified cytokines, chemokines, and surface receptors that could serve as therapeutic targets for colitis and potentially other inflammatory side effects of checkpoint blockade.

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.

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.

Multiple mechanisms have been described that confer BRAF inhibitor resistance to melanomas, yet the basis of this resistance remains undefined in a sizable portion of patient samples. Here, we characterized samples from a set of patients with melanoma that included individuals at baseline diagnosis, on BRAF inhibitor treatment, and with resistant tumors at both the protein and RNA levels. Using RNA and DNA sequencing, we identified known resistance-conferring mutations in 50% (6 of 12) of the resistant samples. In parallel, targeted proteomic analysis by protein array categorized the resistant samples into 3 stable groups, 2 of which were characterized by reactivation of MAPK signaling to different levels and 1 that was MAPK independent. The molecular relevance of these classifications identified in patients was supported by both mutation data and the similarity of resistance patterns that emerged during a co-clinical trial in a genetically engineered mouse (GEM) model of melanoma that recapitulates the development of BRAF inhibitor resistance. Additionally, we defined candidate biomarkers in pre- and early-treatment patient samples that have potential for predicting clinical responses. On the basis of these observations, we suggest that BRAF inhibitor-resistant melanomas can be actionably classified using protein expression patterns, even without identification of the underlying genetic alteration.

Treatment with immune checkpoint blockade (CPB) therapies often leads to prolonged responses in patients with metastatic melanoma, but the common mechanisms of primary and acquired resistance to these agents remain incompletely characterized and have yet to be validated in large cohorts. By analyzing longitudinal tumor biopsies from 17 metastatic melanoma patients treated with CPB therapies, we observed point mutations, deletions or loss of heterozygosity (LOH) in beta-2-microglobulin (B2M), an essential component of MHC class I antigen presentation, in 29.4% of patients with progressing disease. In two independent cohorts of melanoma patients treated with anti-CTLA4 and anti-PD1, respectively, we find that B2M LOH is enriched threefold in non-responders (~30%) compared to responders (~10%) and associated with poorer overall survival. Loss of both copies of B2M is found only in non-responders. B2M loss is likely a common mechanism of resistance to therapies targeting CTLA4 or PD1.

Despite increasing amounts of experimental evidence depicting the involvement of non-coding RNAs in cancer, the study of BRAFV600E-regulated genes has thus far focused mainly on protein-coding ones. Here, we identify and study the microRNAs that BRAFV600E regulates through the ERK pathway.By performing small RNA sequencing on A375 melanoma cells and a vemurafenib-resistant clone that was taken as negative control, we discover miR-204 and miR-211 as the miRNAs most induced by vemurafenib. We also demonstrate that, although belonging to the same family, these two miRNAs have distinctive features. miR-204 is under the control of STAT3 and its expression is induced in amelanotic melanoma cells, where it acts as an effector of vemurafenib's anti-motility activity by targeting AP1S2. Conversely, miR-211, a known transcriptional target of MITF, is induced in melanotic melanoma cells, where it targets EDEM1 and consequently impairs the degradation of TYROSINASE (TYR) through the ER-associated degradation (ERAD) pathway. In doing so, miR-211 serves as an effector of vemurafenib's pro-pigmentation activity. We also show that such an increase in pigmentation in turn represents an adaptive response that needs to be overcome using appropriate inhibitors in order to increase the efficacy of vemurafenib.In summary, we unveil the distinct and context-dependent activities exerted by miR-204 family members in melanoma cells. Our work challenges the widely accepted "same miRNA family = same function" rule and provides a rationale for a novel treatment strategy for melanotic melanomas that is based on the combination of ERK pathway inhibitors with pigmentation inhibitors.

Despite compelling rates of durable clinical responses to programmed cell death-1 (PD-1) blockade, advances are needed to extend these benefits to resistant tumors. We found that tumor-bearing mice deficient in the chemokine receptor CXCR3 responded poorly to anti-PD-1 treatment. CXCR3 and its ligand CXCL9 were critical for a productive CD8+ T cell response in tumor-bearing mice treated with anti-PD-1 but were not required for the infiltration of CD8+ T cells into tumors. The anti-PD-1-induced anti-tumor response was facilitated by CXCL9 production from intratumoral CD103+ dendritic cells, suggesting that CXCR3 facilitates dendritic cell-T cell interactions within the tumor microenvironment. CXCR3 ligands in murine tumors and in plasma of melanoma patients were an indicator of clinical response to anti-PD-1, and their induction in non-responsive murine tumors promoted responsiveness to anti-PD-1. Our data suggest that the CXCR3 chemokine system is a biomarker for sensitivity to PD-1 blockade and that augmenting the intratumoral function of this chemokine system could improve clinical outcomes.