Immunotherapy (IT) and targeted therapy (TT) are both effective against melanoma, but their combination is frequently toxic. Here, we investigated whether the sequence of IT (anti–PD-1)→ TT (ceritinib–trametinib or dabrafenib–trametinib) was associated with improved antitumor responses in mouse models of BRAF- and NRAS-mutant melanoma. Mice with NRAS-mutant (SW1) or BRAF-mutant (SM1) mouse melanomas were treated with either IT, TT, or the sequence of IT→TT. Tumor volumes were measured, and samples from the NRAS-mutant melanomas were collected for immune-cell analysis, single-cell RNA sequencing (scRNA-seq), and reverse phase protein analysis (RPPA). scRNA-seq demonstrated that the IT→TT sequence modulated the immune environment, leading to increased infiltration of T cells, monocytes, dendritic cells and natural killer cells, and decreased numbers of tumor-associated macrophages, myeloid-derived suppressor cells, and regulatory T cells. Durable responses to the IT→TT sequence were dependent on T-cell activity, with depletion of CD8+, but not CD4+, T cells abrogating the therapeutic response. An analysis of transcriptional heterogeneity in the melanoma compartment showed the sequence of IT→TT enriched for a population of melanoma cells with increased expression of MHC class I and melanoma antigens. RPPA analysis demonstrated that the sustained immune response induced by IT→TT suppressed tumor-intrinsic signaling pathways required for therapeutic escape. These studies establish that upfront IT improves the responses to TT in BRAF- and NRAS-mutant melanoma models.

Melanoma is the deadliest form of skin cancer. For many years, it remained refractory to all available therapies. The discovery that approximately 50% of all cutaneous melanomas harbor activating mutations in the serine/threonine kinase, BRAF, led to the development of BRAF inhibitors and later the combination of BRAF and MEK inhibitors (1). These targeted therapies revolutionized the landscape of melanoma therapy, delivering impressive responses in patients whose melanomas harbor activating BRAF mutations. One of the most exciting developments of recent years was the development of effective immunotherapies for melanoma. These strategies, which use therapeutic antibodies to block inhibitory immune checkpoints, enable tumor-reactive T cells to overcome negative regulation and mount effective antitumor responses (2, 3). One of the most successful immunotherapy approaches thus far has been the targeting of programmed cell death (PD)-1, a receptor that maintains peripheral immune tolerance by fine-tuning T-cell responses (2). In the clinic, anti–PD-1 therapy has proven effective in >30% of patients with advanced melanoma, irrespective of tumor genotype, with 70% to 80% of these patients maintaining a response at three years (4).

The success of both targeted therapy and immunotherapy led to attempts to combine these two therapeutic modalities. Constitutive MAPK signaling in melanoma cells allows for immune escape through mechanisms including the recruitment of regulatory T cells (Treg), decreased antigen presentation (via downregulation of MHC class I), and inhibition of IFNγ, IL2, and TNFα release (5–7). Inhibition of BRAF reverses these processes and can potentially restore tumor-immune recognition (8). In preclinical studies, BRAF inhibition leads to increased CD40L expression and interferon-γ release from CD4+ T cells, reduces accumulation of myeloid-derived suppressor cells (MDSC) and Tregs, and decreases IL1, IL6, IL10, and VEGF levels (9–12). In coculture studies of melanoma cells and dendritic cells (DC), BRAF inhibition restores IL12 and TNFα expression and increases levels of T-cell stimulatory molecules, including CD80, CD83, and CD86 (13). In transgenic mouse melanoma models, BRAF inhibition improves the ratio of CD8+ T cells to MDSCs in the tumor (14). In preclinical studies, BRAF–MEK inhibitor combinations synergize with concurrent immune-checkpoint blockade, an effect that associates with decreased macrophage and Treg accumulation, as well as improved IFNγ release and antigen presentation (10).

Initial attempts to develop targeted therapy/immunotherapy combinations clinically (particularly with ipilimumab) have not been successful due to severe toxicity (15). More success has been seen when BRAF and MEK inhibitors are combined with anti–PD-1 therapy (16, 17). Mechanistic studies show this combination to be associated with enhanced CD8+ T-cell accumulation and increased expression of MHC I and II (16, 17). Although promising, these targeted therapy/immunotherapy combinations frequently lead to serious off-target effects (58% grade 3–5 toxicity). It therefore appears that although antitumor immune responses can be improved by combining targeted therapy and immunotherapy, toxicity is a problem that may limit the widespread use of these regimens. One potential strategy is the development of sequential schedules of immunotherapy and targeted therapy, which could deliver more durable antitumor responses with reduced levels of toxicity (under investigation in the clinical trial NCT 03149029). Equivalent approaches for BRAF-wild-type melanoma have not yet been explored. In the current study, we build upon previous work from our group that identified the combination of trametinib and ceritinib as being effective against cell culture models of BRAF-mutant and BRAF-wild-type melanoma cell lines and xenograft models of BRAF-wild-type melanoma (18). As immunotherapy is the mainstay of therapy for BRAF-wild-type melanoma, we asked whether this same drug combination could deliver more durable responses when used in sequence with an anti–PD-1 immune-checkpoint inhibitor.

Cell lines and mice

SW1 (passage #10) and SM1 (passage #5) melanoma cells were obtained from Dr. Eric Lau (Moffitt Cancer Center) in 2017. The NRAS #5 melanoma cell line (passage #3) was derived in 2018 by Dr. Kenneth Tsai from tumors initiated in the laboratory of Dr. Christin Burd (Ohio State) as previously described (19, 20). Cells were kept for a maximum of 10 passages (total). Cell lines were maintained in RPMI-1640 (cat. #MT-10-040-CM; Fisher Scientific) + 10% FBS (cat. #F0926; Sigma). Cell lines were routinely tested for Mycoplasma (every three months) and were authenticated (6-month intervals) by STR authentication (last date of stock testing: 1December 18, 2019). Female, 4-week-old immunocompetent C3H/HeNCrl (Charles River Laboratories) and C57BL/6J (The Jackson Laboratory) mice were observed daily, and all the protocols were reviewed and approved by Institutional Animal Care and Use Committee at University of South Florida (approval #4882R).

In vivo procedures

Mice were subcutaneously injected with 1.5 × 106 cells in Matrigel (cat. #CB40234; Fisher Scientific). The tumors were allowed to grow approximately to ∼50 mm3 before initiation of drug dosing. For immunotherapy, mice received two intraperitoneal doses of anti–PD-1 antibody (200 μg/100 μL; clone RMP1-14; cat. #BE0146; Bio X Cell) or IgG2a isotype control (200 μg/100 μL; clone 2A3; cat. #BE0089; Bio X Cell) every 5 days. On day 10, the immunotherapy was stopped, and the mice received ceritinib (25 mg/kg; cat. #CT-LDK378; Chemietek) or trametinib (1 mg/kg; cat. #CT-GSK212; Chemietek) alone or in combination via oral gavage. A solution of 0.5% carboxy methyl cellulose (cat. #C5678; Sigma-Aldrich) and 0.1% Tween 80 (cat. #AC278632500; Acros Organics) was used as vehicle control. In one study, mice with BRAF-mutant SM1 mouse melanomas were treated with IgG control or anti–PD-1 before being switched onto either control chow (cat. #D10001i; Research Diets) or chow containing dabrafenib (150 mg/kg; cat. #S2807; Selleck Chemicals) and trametinib (1.5 mg/kg; cat. #S2673; Selleck Chemicals; cat. #D20052202i; Research Diets). All the chow diets were manufactured and irradiated by Research Diets, Inc. Tumor size was measured twice weekly. In the in vivo sequence of targeted therapy followed by immunotherapy, mice received a combination of ceritinib and trametinib. Then, on day 10, the targeted therapy was stopped and the mice were switched to intraperitoneal doses of anti–PD-1 antibody or IgG isotype control every 5 days until the experimental endpoint. The tumors were collected at the endpoint, weighed (Supplementary Table S1), and processed for single-cell RNA sequencing (scRNA-seq) or flow cytometry analysis (see below).

CD4+ and CD8+ T-cell depletion

CD4-specific antibody (clone YTS191; cat. #BE0003-1; Bio X Cell) and CD8a-specific antibody (clone YTS169.4; cat. #BE0117; Bio X Cell) were from Bio X Cell were used to deplete CD4+ T cells and CD8+ T cells, respectively. The C3H/HeNCrl mice were administered anti-CD4 and anti-CD8a (100 μg/100 μL) via intraperitoneal injections three days before injection with SW1 cells and then every 4 days thereafter. When the tumors reached 50 to 70 mm3 in size, the mice were treated with two doses of immunotherapy followed by a combination of ceritinib and trametinib every day. Tumor size was measured twice weekly, with CD4+ T-cell or CD8+ T-cell depletion measured by flow cytometry at termination of the experiment.

Flow cytometry

Tumors were harvested at the endpoint under sterile conditions and weighed. Single-cell suspensions were prepared by enzymatic digestion, using a MACS tumor dissociation kit (cat. #130-095-929; Miltenyi Biotec). Numbers of viable cells were counted. To analyze immune-cell populations, 1 × 106 cells were blocked with purified mouse CD16/32 antibody (1:100 dilution; cat. #101301; BioLegend) for 5 minutes on ice. The cells were then incubated with antibody cocktail of Live/Dead Near IR antibody (cat. #L10119; Thermo Fisher Scientific), anti-CD45-BUV395 (clone 30-F11; cat. #565967; BD Biosciences), anti-CD3-BUV737 (clone 17A2, cat. #564380; BD Biosciences), anti-CD4-BUV 496 (clone GK1.5; cat. #564667; BD Biosciences), anti-CD8-BUV805 (clone 53-6.7; cat. #564920; BD Biosciences), CD127-BV711 (clone SB/199; cat. #565490; BD Biosciences), CD69-AF488 (clone H1.2F3; cat. # 104516; BioLegend), CD44-APCR700 (clone IM7; cat. #565480; BD Biosciences), CD62L-BV650 (clone MEL-14; cat. #564108; BD Biosciences), PD-1-BV785 (clone 29F.1A12; cat. #135225; BioLegend), CTLA4-BV421 (clone UC10-4B9; cat. #106311; BioLegend), TIM3-PECF594 (clone B8.2C12; cat. #134013; BioLegend), LAG3-PE (clone C9B7W; cat. #552380; BD Biosciences) for T-cell analysis and Live/Dead Near IR antibody (cat. #L10119; Thermo Fisher Scientific), anti-CD45-BUV395 (clone 30-F11; cat. #565967; BD Biosciences), anti-CD3-BUV737 (clone 17A2, cat. #564380; BD Biosciences), CD11b-BB700 (M1/70, cat. #566417; BD Biosciences), Gr-1-PE-Cy7 (RB6-8C5; cat. #108415; BioLegend), anti-Ly6C-BV421 (clone HK1.4; cat. #128031; BioLegend), anti-Ly6G-APC (clone 1A8; cat. #127613; BioLegend), anti-CD11c-BV605 (clone N418; cat. #117333; BioLegend), anti-MHC II-BB515 (clone 2G9; cat. #565254; BD Biosciences), anti-F4/80-BV785 (clone BM8; cat. #123141; BioLegend), and anti-CD103-PE (clone M290; cat. #561043; BD Biosciences) for myeloid cells analysis. Most of the antibodies were used at 1:50 to 1:100 dilutions. The cells were incubated with the antibody cocktail for 20 minutes at 4°C in dark. For FOXP3 staining, the cells were fixed/permeabilized overnight at 4°C in dark using eBioscience FOXP3 transcription factor staining buffer set (cat. #00-5523-00; Thermo Fisher) and FOXP3 monoclonal antibody (1:25 dilution; clone FJK-16s; cat. #17-5773-80; Thermo Fisher). All the washings were done with PBS (cat. #SH30256FS; Fisher Scientific) + 2% FBS (cat. #F0926; Sigma). Flow cytometry acquisition was performed on the BD FACS Symphony or LSR II. The data analysis was carried out using FlowJo software. To detect melanoma antigen expression on tumor cells or cell lines treated with kinase inhibitors, the cells were incubated with anti-Tyrp1-APC (clone TA99; cat. # NBP2-34720APC; Novus Biologicals) at 1:50 dilution for 20 minutes at 4°C in dark. Flow cytometry acquisition was performed on the BD FACS LSR II. The data analysis was carried out using FlowJo software.

scRNA-seq

Tumors were harvested at the endpoint under sterile conditions and weighed. Single-cell suspensions were prepared by enzymatic digestion, using a MACS tumor dissociation kit (cat. #130-095-929; Miltenyi Biotec). Cells were strained through MACS strainer (cat. # 130-098-458; Miltenyi Biotec). The cell count and viability was analyzed by staining the cells with AO/PI stain on the Nexcelom Cellometer K2. The cells were then resuspended at a concentration of 500 cells/μL in PBS (cat. #SH30256FS; Fisher Scientific) + 0.4% nonacetylated BSA (cat. #BP1605100; Fisher Scientific). The samples were then loaded onto 10X Genomics Chromium Single-Cell Controller (10X Genomics) to prepare scRNA-seq libraries. Around 50,000 to 1,000,000 mean sequencing reads per cell were generated on Illumina NextSeq 500 instrument using v2.5 flow cells. 10X Genomics CellRanger software was used for demultiplexing, barcode processing, alignment, and gene counting. Finally, the analysis of single-cell data sets was performed using Interactive Single-Cell Visual Analytics (ISCVA; see below).

ISCVA

To facilitate the rapid analysis of single-cell data sets, we developed a new computational tool consisting of two major components. The first component is comprised of a collection of Bash and R scripts (utilizing many of the widely used algorithms in the single-cell community, including Seurat for general processing (21), SingleR for cell-type recognition (22), and single-cell signature explorer for gene set signature scoring; ref. 23) that processed the scRNA-seq data offline, and a second Web-based component (based upon state-of-the-art technologies, including react.js from Facebook, tensorflow.js from Google and Plotly.js) that allows convenient real-time interactive exploration and ad hoc analysis. The heterogeneity analyses implemented in SinCHet (24) were also performed as part of the analytical modules. A node.js backend was also created to serve the on-demand queries of the Web application, allowing for real-time interactive investigation of genes expressed in selected samples or subsets of cells. Cells with high mitochondria content were not filtered as these may reflect cell populations going through apoptosis. Data are available through Gene Expression Omnibus (GEO; GSE165582).

Western blot analysis

Western blotting was performed as previously described (25). Primary antibodies for phospho-ERK (clone D13.14.4E; cat. #4370S), total ERK (clone 137F5; cat. #4695S), phospho-AKT (clone D9E; cat. #4060S), total AKT (clone 40D4; cat. #2920S), phospho-Met (clone D26; cat. #3077S), and total Met (clone D1C2; cat. #8198S) from Cell Signaling Technology.

Reverse phase protein analysis

Frozen cell pellets obtained from enzymatic digestion of mouse melanoma tumors, using a MACS tumor dissociation kit (cat. #130-095-929; Miltenyi Biotec) were used to extract proteins for reverse phase protein analysis (RPPA). Proteins were extracted using RIPA lysis buffer containing beta-mercaptoethanol (cat. #M7522; Sigma-Aldrich) and cOmplete, Mini Protease Inhibitor Cocktail (cat. #11836153001, Millipore Sigma). Protein concentration was determined using a Pierce BCA protein assay kit (cat. #23227; Thermo Scientific). The lysates were submitted at a concentration of 80 μg/80 μL to Functional Proteomics RPPA core facility, MD Anderson. All the further sample processing and data analysis for RPPA was carried out as previously described (26).

MTT assay

The cells were seeded in 96-well plates at a concentration of 2,000 cells/100 μL/well overnight. Next day, the cells were treated with increasing doses of ceritinib or trametinib. The cells were allowed to incubate with the drugs for 72 hours. MTT stock solution was prepared from MTT reagent (cat. #M5655; Sigma-Aldrich) at a concentration of 5 mg/mL in water. MTT solution (100 μL) was added to the cells per well after aspirating the medium. The cells were incubated with MTT solution for 3 to 4 hours. Followed by complete removal of MTT solution after 3 to 4 hours, 100 μL of DMSO was added to dissolve the crystals. The plates were read at 490 nm using Synergy H1 microplate reader from BioTek.

Statistical analysis

One-way ANOVA in Microsoft Excel version 15.40 was used to compare the results between different groups with a single independent variable. The mean of three independent experiments ± SEM is shown for each data set. Wilcoxon rank sum test results within each cluster were used to compare results between groups for scRNA-seq. Results with P values ≤ 0.05 were considered statistically significant.

The IT→TT sequence is more effective than either IT or TT alone in mouse models of NRAS- and BRAF-mutant melanoma

We began by investigating the effects of continuously dosed anti–PD-1 (this immunotherapy is referred to in the Results section as IT) therapy in a syngeneic mouse model of NRAS-mutant melanoma (SW1 cells) and noted a delay in tumor growth but no tumor regression (Supplementary Fig. S1). It therefore seemed that this model was at least partly resistant to IT. In previous studies, we noted that the combination of ceritinib–trametinib (targeted therapy referred to in the results section as TT) led to impressive levels of tumor suppression in BRAF-mutant and BRAF-wild-type melanoma cell lines (18). Thus, we next asked whether use of IT and TT in sequence would lead to improved therapeutic responses compared with either therapy alone and designed two schedules based on 10 days of either IT or TT followed by a switch to TT or IT, respectively (Fig. 1A). Evaluation of trametinib and ceritinib both alone and in combination in the SW1 mouse melanoma model demonstrated that trametinib led to initial tumor shrinkage, but then tumor growth resumed (Fig. 1B). Use of the ceritinib–trametinib combination (TT) led to more durable responses than the monotherapy, with some recovery of tumor growth (Fig. 1B). We next evaluated the effects of ceritinib alone, trametinib alone, and the TT combination following two doses of IT (Fig. 1C). It was noted that upfront use of IT improved responses to both mono- and combination TT therapies (Fig. 1C). Longer-term treatment of up to 41 days demonstrated that the IT→TT sequence gave more durable responses than just TT alone, with no evidence of treatment failure in the SW1 model (Fig. 1D). Repeat of this treatment schedule in a syngeneic BRAF-mutant mouse melanoma model (SM1 cells) showed similar responses and demonstrated that initial IT dramatically increased the efficacy of TT compared with the IgG→TT sequence (Fig. 1E). In this instance, IT alone was more effective than it was in the SW1 model. Reponses to TT alone were observed, but resistance rapidly occurred (Fig. 1E). Similarly, more durable antitumor responses were seen when the BRAF-mutant SM1 melanoma model was treated with IT→dabrafenib–trametinib compared with vehicle (veh), IT alone, or dabrafenib–trametinib alone (Supplementary Fig. S2). Investigation of the IT→TT sequence in a second mouse model of NRAS-mutant melanoma (NRAS #5) revealed a similar potential to improve TT responses (Supplementary Fig. S3). As the IT→TT sequence delayed the onset of TT resistance, we next asked whether the reverse sequence of TT→IT would be equally effective (see Fig. 1A for scheme). It was noted that treatment with the TT combination for 10 days until tumor regression followed by IT improved responses compared with TT alone, but did not result in the rapid tumor regressions seen to the IT→TT sequence in the SW1 model (Fig. 1F). Together, these data suggested that initial use of IT enhanced subsequent responses to TT in multiple mouse melanoma models.

Figure 1.

The sequence of IT→TT suppresses the growth of NRAS-mutant (SW1) and BRAF-mutant (SM1) tumors and limits relapse. A, Dosing schema for the IT→TT and TT→IT sequences. B, TT alone causes tumor suppression up to day 30 of treatment. Mice received two doses of IgG control (200 μg/100 μL), followed by vehicle, ceritinib (25 mg/kg), trametinib (1 mg/kg), or a combination of ceritinib and trametinib for 30 days. C, Treatment with IT→TT yields durable antitumor responses. Mice received two doses of anti–PD-1 (200 μg/100 μL), followed by vehicle, ceritinib, trametinib, or ceritinib and trametinib combination therapy for 30 days. D, IT→TT sequence leads to long-term antitumor responses. Mice were treated with IgG or IT, followed by vehicle or combination TT therapy for a period of 41 days. E, Treatment with IT→TT in BRAF-mutant melanoma mouse model system (SM1) shows antitumor responses without relapse. Mice were treated as per the SW1 melanoma in C for 41 days. F, The reverse sequence TT→IT was associated with less tumor suppression in SW1 NRAS-mutant melanoma. Mice received 10 doses of vehicle or a combination of ceritinib and trametinib every day, followed by six doses of anti–PD-1 or IgG control (200 μg/100 μL) every 5 days for a period of 35 days. The results are represented as average ± SEM. For all panels, eight mice per group were used except for E and F, where five mice per group were used. Statistical significance was assessed with one-way ANOVA test (**, 0.0001 ≤ P ≤ 0.001; ***, P ≤ 0.0001).

Figure 1.

The sequence of IT→TT suppresses the growth of NRAS-mutant (SW1) and BRAF-mutant (SM1) tumors and limits relapse. A, Dosing schema for the IT→TT and TT→IT sequences. B, TT alone causes tumor suppression up to day 30 of treatment. Mice received two doses of IgG control (200 μg/100 μL), followed by vehicle, ceritinib (25 mg/kg), trametinib (1 mg/kg), or a combination of ceritinib and trametinib for 30 days. C, Treatment with IT→TT yields durable antitumor responses. Mice received two doses of anti–PD-1 (200 μg/100 μL), followed by vehicle, ceritinib, trametinib, or ceritinib and trametinib combination therapy for 30 days. D, IT→TT sequence leads to long-term antitumor responses. Mice were treated with IgG or IT, followed by vehicle or combination TT therapy for a period of 41 days. E, Treatment with IT→TT in BRAF-mutant melanoma mouse model system (SM1) shows antitumor responses without relapse. Mice were treated as per the SW1 melanoma in C for 41 days. F, The reverse sequence TT→IT was associated with less tumor suppression in SW1 NRAS-mutant melanoma. Mice received 10 doses of vehicle or a combination of ceritinib and trametinib every day, followed by six doses of anti–PD-1 or IgG control (200 μg/100 μL) every 5 days for a period of 35 days. The results are represented as average ± SEM. For all panels, eight mice per group were used except for E and F, where five mice per group were used. Statistical significance was assessed with one-way ANOVA test (**, 0.0001 ≤ P ≤ 0.001; ***, P ≤ 0.0001).

Close modal

scRNA-seq defines the effect of each therapy sequence on the immune microenvironment

We next used scRNA-seq to determine how different therapy sequences affected the immune–tumor landscape. A high-level overview of cellular composition identified multiple cell types in each tumor, including melanoma cells, T cells, natural killer cells, myeloid cells, granulocytes, monocytes, endothelial cells, and fibroblasts (Fig. 2A and B). Each therapy sequence caused marked changes to the melanoma compartment (Fig. 2B). Treatment-dependent alterations in the extent of immune infiltration and the presence of myeloid cells were also evident. We saw that TT increased the proportion of fibroblasts in the tumor, IT→Veh increased T-cell accumulation, and the IT→TT sequence led to the largest immune-cell influx (including monocytes, granulocytes, and T cells) along with higher numbers of fibroblasts and endothelial cells (Fig. 2C). Detailed cell curation and parallel flow cytometry were used to better understand the constituent cell types in the tumor microenvironment (Fig. 2D–F; Supplementary Fig. S4). The greatest number of CD4+ and CD8+ T cells were seen in flow cytometry data from tumors treated with the IT→TT sequence (Fig. 2E and F). Furthermore, T cells from IT→TT-treated tumors had higher CD69 and IFNG levels, as measured by scRNA-seq (Supplementary Fig. S5). The IT→TT sequence was also associated with the decreased infiltration of regulatory CD4+ T cells (Treg) and immune suppressive (CD11b+GR1+) MDSCs (Fig. 2G and H). It was further noted that the IgG→TT sequence also reduced Treg and MDSC numbers, whereas the IT→Veh did not (Fig. 2G and H). We next examined the tumors from the reverse TT→IT sequence (Supplementary Fig. S6) and observed fewer tumor-infiltrating CD4+ or CD8+ T cells (∼100-fold less than IT→TT; Supplementary Fig. S6A). Although the TT→IT sequence had some impact on Treg numbers, this was less than that seen with the IT→TT sequence (Supplementary Fig. S6B). MDSC numbers either did not change or slightly increased with TT→IT compared with the dramatic decreases seen after IT→TT administration (Supplementary Fig. S6C).

Figure 2.

IT→TT increases T-cell infiltration and reduces the burden of Tregs and MDSCs in NRAS-mutant tumors. A, t-SNE plots showing the major cell types and subtypes identified in SW1 tumors from the indicated treatment groups. NK, natural killer. B, t-SNE plots showing the distribution of all cell subtypes in SW1 tumors from the indicated treatment groups. C, Number (top) and proportion (bottom) of each cell type in tumors from the indicated treatment groups. D, Heat map showing the expression of defining markers across the 19 identified cell subclusters. E, Total number of tumor-infiltrating CD4+ and CD8a+ T cells as measured by flow cytometry. F, Number of tumor-infiltrating CD4+ and CD8a+ T cells after normalization to tumor weight (see E). The results are represented as average ± SEM of three mice per group. G, Percentage of tumor-infiltrating CD4+ T cells that stained positive for the Treg marker FoxP3. H, Percentage of live cells staining positive for the MDSC markers CD11b and Gr-1. The results are represented as average ± SEM of three mice per group. One representative experiment of three is shown. Statistical significance was assessed with one-way ANOVA test (*, 0.001 ≤ P ≤ 0.05; **, 0.0001 ≤ P ≤ 0.001).

Figure 2.

IT→TT increases T-cell infiltration and reduces the burden of Tregs and MDSCs in NRAS-mutant tumors. A, t-SNE plots showing the major cell types and subtypes identified in SW1 tumors from the indicated treatment groups. NK, natural killer. B, t-SNE plots showing the distribution of all cell subtypes in SW1 tumors from the indicated treatment groups. C, Number (top) and proportion (bottom) of each cell type in tumors from the indicated treatment groups. D, Heat map showing the expression of defining markers across the 19 identified cell subclusters. E, Total number of tumor-infiltrating CD4+ and CD8a+ T cells as measured by flow cytometry. F, Number of tumor-infiltrating CD4+ and CD8a+ T cells after normalization to tumor weight (see E). The results are represented as average ± SEM of three mice per group. G, Percentage of tumor-infiltrating CD4+ T cells that stained positive for the Treg marker FoxP3. H, Percentage of live cells staining positive for the MDSC markers CD11b and Gr-1. The results are represented as average ± SEM of three mice per group. One representative experiment of three is shown. Statistical significance was assessed with one-way ANOVA test (*, 0.001 ≤ P ≤ 0.05; **, 0.0001 ≤ P ≤ 0.001).

Close modal

Responses to the IT→TT sequence are dependent upon an active CD8+ T-cell response

We next determined whether the durable responses seen with the IT→TT sequence were dependent on T-cell activity. CD4+ or CD8+ T cells were depleted using blocking antibodies prior to the injection of tumor cells. Once tumors had formed, the mice were treated with either IgG→Veh, IgG→TT, or IT→TT. It was found that depletion of CD8+ T cells, but not CD4+ T cells, was required for the activity of TT, even without prior treatment with IT (Fig. 3A and B). To better understand the T-cell phenotypes required for durable responses to the IT→TT sequence, we performed a detailed analysis of the scRNA-seq data and identified four populations of T cells, including one subcluster of CD4+ T cells and three clusters of CD8+ T cells. The CD4+ T cells identified expressed multiple activation markers including Cd40lg, Icos, Cd5, and Maf (Fig. 3C and D; Supplementary Fig. S7). In the CD8+ T-cell clusters, cluster #1 was activated CD8+ T effector cells, expressing Ifngγ, multiple granzymes (Gzm), Ccl4, and Ccl5, cluster #2 was a subset of exhausted CD8+ T cells expressing multiple inhibitory markers including ApoD, Msp, Sparc, and Mt1, as well as some activation markers including Gzmk and Cd38 (Fig. 3C and D; Supplementary Fig. S7). Cluster #3 was characteristic of a rapidly proliferating population of CD8+ T cells that expressed genes associated with memory T cells (increased Birc5 and Hist1h3c; lower levels of Sell and Tcf7; Fig. 3C; Supplementary Fig. S7). Pathway analysis of these three CD8+ T-cell subsets confirmed these predictions (Supplementary Fig. S8A). Quantification of the T-cell composition following treatment with each therapy sequence demonstrated that the IT→TT sequence was associated with the largest accumulation of T cells and B cells (Fig. 3E). The T-cell proportions seen following IT→TT were similar to the IgG→Veh group and the IgG→TT sequence (albeit with far fewer total T cells seen for the latter two sequences). By contrast, the IT→Veh sequence was associated with the most dramatic change in T-cell composition, with an increased accumulation of cluster #2 CD8+ T cells (Fig. 3E). Orthogonal validation by flow cytometry analysis demonstrated that the IT→TT sequence was associated with the highest levels of infiltrating CD44+CD62L+, CD44+CD62L, and CD69+ effector T cells, and reduced expression of exhaustion markers TIM3 on the tumor-infiltrating CD8+ T cells (Fig. 3F and G; Supplementary Fig. S8B and S8C). We also examined the tumors from the reverse TT→IT sequence for CD44+CD62L+, CD44+CD62L, and CD69+ effector T cells (Supplementary Fig. S8D and S8E). The TT→IT sequence was associated with virtually no CD44+CD62L CD8+ T cells and drastically reduced numbers of CD44+CD62L+ and CD69+ effector T cells compared with the IT→TT sequence (Supplementary Fig. S8D and S8E).

Figure 3.

The IT→TT sequence increases the accumulation of activated CD8+ T cells. A and B, Responses to IT→TT are dependent on CD8+ T cells. Mice were treated with CD4- or CD8-specific antibodies (i.p., 100 μg/100 μL) every 4 days before beginning the IT→TT sequence for 23 days. The results are represented as average ± SEM of five mice per group. C, Unsupervised hierarchical clustering of scRNA-seq data from T and natural killer (NK) cells derived from vehicle- and drug-treated SW1 tumors. D, Heat map showing expression of genes associated with T-cell activation and exhaustion. E, Different types of lymphoid cells (T cells and NK cells) identified in tumors with different treatment condition by cell-type count (top) and cell-type proportion (bottom). F, T-cell activation as shown by CD44 and CD62 L staining. G, Quantification of TIM3 expression on tumor-infiltrating T cells after each therapy sequence. Histograms and bar graphs show Tim3 mean fluorescence intensity (MFI). The results are represented as average ± SEM of three mice per group for F and G. One representative experiment of three is shown. Statistical significance was assessed with one-way ANOVA test (*, 0.001 ≤ P ≤ 0.05; ***, P ≤ 0.0001).

Figure 3.

The IT→TT sequence increases the accumulation of activated CD8+ T cells. A and B, Responses to IT→TT are dependent on CD8+ T cells. Mice were treated with CD4- or CD8-specific antibodies (i.p., 100 μg/100 μL) every 4 days before beginning the IT→TT sequence for 23 days. The results are represented as average ± SEM of five mice per group. C, Unsupervised hierarchical clustering of scRNA-seq data from T and natural killer (NK) cells derived from vehicle- and drug-treated SW1 tumors. D, Heat map showing expression of genes associated with T-cell activation and exhaustion. E, Different types of lymphoid cells (T cells and NK cells) identified in tumors with different treatment condition by cell-type count (top) and cell-type proportion (bottom). F, T-cell activation as shown by CD44 and CD62 L staining. G, Quantification of TIM3 expression on tumor-infiltrating T cells after each therapy sequence. Histograms and bar graphs show Tim3 mean fluorescence intensity (MFI). The results are represented as average ± SEM of three mice per group for F and G. One representative experiment of three is shown. Statistical significance was assessed with one-way ANOVA test (*, 0.001 ≤ P ≤ 0.05; ***, P ≤ 0.0001).

Close modal

IT→TT alters the myeloid cell landscape

We next investigated how each therapy sequence modulated the myeloid cell compartment. Our analyses identified one population of granulocytes, three populations of macrophages, and four populations of cells with characteristics of both monocytic cells and DCs (Fig. 4A; Supplementary Fig. S9). Detailed cell curation suggested that macrophage clusters #1 and #2 had features of tumor-associated macrophage (TAM) and M2 macrophage phenotypes, and cluster #3 represented a hybrid phenotype with characteristics of M1 and M2 macrophages. In the monocytic cell–DC subsets, cluster #1 was characteristic of monocytes differentiating to DCs, cluster #2 contained circulating monocytes, and cluster #3 represented cDC1s as they expressed DC development markers such as Irf8 and Cd80. DCs in cluster #3 also expressed Integrin αE and Ccr7, which are necessary to direct DCs to tumor-draining lymph nodes for antigen cross presentation (Fig. 4A). Cluster #4 contained cells characteristic of activated immune stimulatory cDC2s, which express markers necessary for DC–T-cell cross-talk and the activation of cytotoxic T cells (including Ifngr2, Cd40, Tlr2, and Cd86; Fig. 4A).

Figure 4.

IT→TT increases DC infiltration and upregulates MHC class I expression across multiple cell types. A, Unsupervised hierarchical clustering of granulocytes, monocytes, and macrophages from vehicle- and drug-treated SW1 tumors (top). Heat map showing the expression of DC activation markers across monocyte clusters (bottom). B, Curation of the myeloid cell types identified in tumors following each therapy regimen by cell-type count (top) and cell-type proportion (bottom). C, IT→TT increases the tumor infiltration of antigen-presenting DCs (DC1s and DC2s). D, Average number of DC1s and DC2s normalized to tumor weight. The cells that stained positive for CD11c and CD103, and negative for CD11b and F4/80 were designated as DC1. The results are represented as average ± SEM of three mice per group. E, t-SNE plots showing major MHC class I markers, B2m and H2-D1, identified in immune and tumor cell clusters with different treatment conditions (top). Expression of MHC class I markers in each immune cell and melanoma cell cluster identified by violin plots (bottom). F, IT→TT increases melanoma cell-surface MHC class I expression. SW1 tumors were digested, and cells were stained with MHC I antibody before being analyzed by flow cytometry. The results are represented as average ± SEM of three mice per group. One representative experiment of three is shown. Statistical significance was assessed with one-way ANOVA test. Wilcoxon rank sum test results within each cluster were used to compare results between groups for violin plots (*, 0.001 ≤ P ≤ 0.05; **, 0.0001 ≤ P ≤ 0.001; ***, P ≤ 0.0001).

Figure 4.

IT→TT increases DC infiltration and upregulates MHC class I expression across multiple cell types. A, Unsupervised hierarchical clustering of granulocytes, monocytes, and macrophages from vehicle- and drug-treated SW1 tumors (top). Heat map showing the expression of DC activation markers across monocyte clusters (bottom). B, Curation of the myeloid cell types identified in tumors following each therapy regimen by cell-type count (top) and cell-type proportion (bottom). C, IT→TT increases the tumor infiltration of antigen-presenting DCs (DC1s and DC2s). D, Average number of DC1s and DC2s normalized to tumor weight. The cells that stained positive for CD11c and CD103, and negative for CD11b and F4/80 were designated as DC1. The results are represented as average ± SEM of three mice per group. E, t-SNE plots showing major MHC class I markers, B2m and H2-D1, identified in immune and tumor cell clusters with different treatment conditions (top). Expression of MHC class I markers in each immune cell and melanoma cell cluster identified by violin plots (bottom). F, IT→TT increases melanoma cell-surface MHC class I expression. SW1 tumors were digested, and cells were stained with MHC I antibody before being analyzed by flow cytometry. The results are represented as average ± SEM of three mice per group. One representative experiment of three is shown. Statistical significance was assessed with one-way ANOVA test. Wilcoxon rank sum test results within each cluster were used to compare results between groups for violin plots (*, 0.001 ≤ P ≤ 0.05; **, 0.0001 ≤ P ≤ 0.001; ***, P ≤ 0.0001).

Close modal

An analysis of the myeloid cell numbers under each therapy sequence demonstrated a marked increase in cell infiltration following IT→Veh and IT→TT. Important differences were seen between the two regimens, with IT→Veh being associated with greater infiltration of macrophage clusters #1 and #2 (Fig. 4B). Cluster #1 macrophages expressed Apoe, Ms4a7, Trem2, Cxcl14, and Cd72, which are the markers for M2 macrophages and TAMs. Cluster #2 macrophages expressed Sparc, which is secreted by TAMs, as well as Serpin1, Cryab, and Cpe, which are also markers for TAMs. In contrast, the IT→TT sequence was associated with increases in granulocytes and cells from monocytic/DC clusters #1, #2, #3, #4, and in particular the highest total accumulation of cDC1 and cDC2 (Fig. 4B). Significantly, the IT→TT sequence was the only therapy regimen to be associated with large numbers of cDC2s. Tumors from mice treated with IT→TT also had fewer macrophages, which correlated with elevated numbers of cDC2s and, to a lower extent, expansion of cDC1s (Fig. 4B). Tumors treated with IgG→TT had fewer macrophages than IT→Veh, and slightly more DCs. Together, these data suggest that priming with IT reduces the accumulation of immune suppressive macrophages and increases the number of antigen-presenting DCs in TT-treated animals. Flow cytometry confirmed that the IT→TT sequence increased the number of tumor-associated cDC1s and cDC2s (Fig. 4C and D). As the increased expression of MHC I on DCs plays a major role in the activation of protective antitumor CD8+ T-cell responses, we interrogated our single-cell data for the expression of MHC class I mRNAs. The IT→Veh sequence increased the expression of mRNAs such as B2m, H2-D1, and H2-K1, and these transcripts were further upregulated following the IT→TT sequence (Fig. 4E; Supplementary Fig. S10A and S10B). Interestingly, the increased expression of B2m was observed in melanoma cells as well as multiple immune-cell populations (T cells, monocytes, and granulocytes) following IT→TT treatment. A similar, but less pronounced, trend was seen for H2-D1 (Fig. 4E). Flow cytometry validated these findings, showing that the IT→TT sequence dramatically increased cell-surface expression of MHC class I in the tumor (Fig. 4F). These data support the idea that IT alone upregulates MHC class I in the tumor microenvironment and that these effects can be augmented by subsequent TT. We examined the tumor-associated cDC1 and cDC2 by flow cytometry in tumors from the reverse TT→IT sequence (Supplementary Fig. S10C). Our results confirmed reduced number of cDC1 and cDC2 infiltrations in these tumors compared with the IT→TT sequence.

IT→TT enriches for melanoma cells with immune-responsive signatures

Antitumor immune responses are dependent on the recognition of tumor antigens by antigen-presenting DCs. We reasoned that this could be mediated in part by increased antigen expression in the melanoma cells following IT→TT. To explore this at a single-cell resolution, we determined how each therapy sequence modulated the transcriptional heterogeneity of the tumor cells (Fig. 5A). These analyses utilized the single-cell heterogeneity (SinCHet) software platform developed by our group (24). Using this approach, we used the minimum change point at the lowest cluster level to define 13 distinct transcriptional states in the SW1 melanoma tumors (Fig. 5A and B). The transcriptional clusters we identified had gene signatures that predicted discrete cellular processes, with each therapy sequence altering the transcriptional heterogeneity in a unique manner (Fig. 5B). Treatment with TT led to a massive shift in the melanoma cell transcriptional profiles, which then differed if the tumors were pretreated with IT. In particular, we observed that IgG→TT was primarily associated with an increase in cluster #2, whereas the IT→TT sequence enriched for cluster #5. Both sequences showed an increase in cluster #11. A detailed analysis revealed expression of genes implicated in immune regulation and antigen expression such as Dct, Tyrp1, and Sox10 in clusters #2, #5, and #11 (Fig. 5B–D). Notably, cluster #5 was also associated with increased expression of B2m and H2-D1, suggesting this may represent a particularly immunogenic cluster of tumor cells (Fig. 4E). In contrast, the major clusters observed in the Veh→IgG-treated tumors (such as #1) were mostly associated with glycolysis, metastasis, and drug resistance (Fig. 5B and C). Treatment with IT→Veh enriched for cluster #4, which was associated with glycolysis, metabolism, innate immunity, and transcription/translation. An analysis of SW1 mouse tumor samples by flow cytometry confirmed that the IT→TT sequence dramatically increased the cell-surface expression of TYRP1 and the melanoma lineage marker SOX10 (Fig. 5E). These effects were recapitulated in SW1 melanoma cells treated with trametinib (and to a lesser extent ceritinib) for 48 to 72 hours in vitro (Supplementary Fig. S11).

Figure 5.

IT→TT enriches for melanoma cells with increased antigen expression. A, UMAP plots showing how each therapy sequence alters heterogeneity in the melanoma compartment. SinCHET analysis identified 13 melanoma subclusters, allowing the transcriptional composition of the melanoma compartment to be defined following each treatment regimen. B, Unsupervised hierarchical clustering showing the genes that differentiate each melanoma cluster. C, STRING analysis of genes differentially expressed in clusters 1 and 2 (vehicle cluster) and clusters 4 and 5 (combination cluster). Analysis was performed using STRING: functional protein association networks database (http://string-db.org/). D, Violin plot showing expression of the melanoma antigens Sox10, Tyrp1, and Dct in each melanoma cluster identified. E, Quantification of TYRP1 and SOX10 expression in SW1 tumors following each therapy regimen. Single-cell suspensions of tumor cells were stained with anti-TYRP1 or anti-SOX10 and quantified by flow cytometry. The results are represented as average ± SEM of three mice per group. Statistical significance was assessed with one-way ANOVA test. Wilcoxon rank sum test results within each cluster were used to compare results between groups for violin plots (*, 0.001 ≤ P ≤ 0.05).

Figure 5.

IT→TT enriches for melanoma cells with increased antigen expression. A, UMAP plots showing how each therapy sequence alters heterogeneity in the melanoma compartment. SinCHET analysis identified 13 melanoma subclusters, allowing the transcriptional composition of the melanoma compartment to be defined following each treatment regimen. B, Unsupervised hierarchical clustering showing the genes that differentiate each melanoma cluster. C, STRING analysis of genes differentially expressed in clusters 1 and 2 (vehicle cluster) and clusters 4 and 5 (combination cluster). Analysis was performed using STRING: functional protein association networks database (http://string-db.org/). D, Violin plot showing expression of the melanoma antigens Sox10, Tyrp1, and Dct in each melanoma cluster identified. E, Quantification of TYRP1 and SOX10 expression in SW1 tumors following each therapy regimen. Single-cell suspensions of tumor cells were stained with anti-TYRP1 or anti-SOX10 and quantified by flow cytometry. The results are represented as average ± SEM of three mice per group. Statistical significance was assessed with one-way ANOVA test. Wilcoxon rank sum test results within each cluster were used to compare results between groups for violin plots (*, 0.001 ≤ P ≤ 0.05).

Close modal

IT→TT suppresses the emergence of a resistance signature in NRAS-mutant melanoma

Our in vivo experiments demonstrated that the IT→TT sequence was associated with improved antitumor responses in the SW1 model. We performed RPPA on tumor samples collected at days 5 and 15 to determine if upfront IT limited subsequent TT failure by altering signaling in the tumor. Key differences were noted in the signaling patterns of tumors from mice treated with IgG→TT versus IT→TT (Fig. 6A). Pathways that showed significant differences included those involved in cell death and signaling adaptation (Fig. 6B and C). Of note, tumors treated with the IgG→TT sequence had a signaling signature associated with therapeutic escape including increased levels of phospho-RAF, AKT, MAPK, and multiple receptor tyrosine kinases (including c-MET, IGF1R, HER1, and HER2; ref. Fig. 6B). A potential link between the recovery of mitogenic signaling in the IgG→TT-treated tumors and immune escape was suggested by studies of cells generated from SW1 tumors that had failed IgG→TT therapy. These cells, which showed resistance to trametinib, but not ceritinib (Supplementary Fig. S12A and S12B), exhibited little change in AKT or MET signaling following TT treatment (Fig. 6D) and had lower baseline expression of melanoma antigens compared with their drug-naïve counterparts (Fig. 6E). Upon further drug treatment, only minor increases in antigen expression were seen in the TT-resistant SW1 cells compared with the drug-naïve SW1 cells (Fig. 6E). A role for increased MAPK signaling in the immune escape was suggested by studies demonstrating that treatment of the resistant cells with the ERK inhibitor SCH772984 (which is known to limit MEK inhibitor resistance in NRAS-mutant melanoma cells; refs. 27, 28) increased antigen expression, albeit to a lesser degree than in the drug-naïve SW1 cells (Fig. 6F). Together, these results suggest that escape from targeted therapy is accompanied by a decrease in immune recognition, further supporting a role for the immune system in mediating continued targeted therapy responses.

Figure 6.

IT→TT suppresses therapeutic escape in the tumor compartment. A, RPPA analysis of SW1 tumors following each therapy sequence identified changes in resistance-associated signaling following the IT→TT sequence compared with TT alone. B, Detailed view of RPPA highlighting phospho-proteins involved in therapy escape and tumor progression. C, Pathway analysis of the genes from the RPPA heat map (B) using GeneGo MetaCore software, highlighting the upregulation of pathways involved in cell proliferation and therapy resistance in tumors treated with IgG→TT. D, Western blot analysis showing increased AKT and MET signaling in SW1 cells with acquired resistance to TT. Cells were treated with the ceritinib–trametinib combination (500 and 1 nmol/L, respectively, 72 hours) and blotted for the expression of phospho-AKT and phospho-MET. E, Quantification of TYRP1 antigen expression on drug-naïve and TT-resistant SW1 cells. The cells were treated with ceritinib (500 nmol/L), trametinib (1 nmol/L), or combination therapy (72 hours) before being stained with TYRP1 surface antibody. F, Quantification of TYRP1 antigen expression on drug-naïve and TT-resistant SW1 cells. The cells were treated with the ERK inhibitor SCH772984 (500 nmol/L) for 72 hours before being stained with TYRP1 surface antibody. The results are represented as average ± SEM of three independent experiments. Statistical significance was assessed with one-way ANOVA test (*, 0.001 ≤ P ≤ 0.05; **, 0.0001 ≤ P ≤ 0.001; ***, P ≤ 0.0001).

Figure 6.

IT→TT suppresses therapeutic escape in the tumor compartment. A, RPPA analysis of SW1 tumors following each therapy sequence identified changes in resistance-associated signaling following the IT→TT sequence compared with TT alone. B, Detailed view of RPPA highlighting phospho-proteins involved in therapy escape and tumor progression. C, Pathway analysis of the genes from the RPPA heat map (B) using GeneGo MetaCore software, highlighting the upregulation of pathways involved in cell proliferation and therapy resistance in tumors treated with IgG→TT. D, Western blot analysis showing increased AKT and MET signaling in SW1 cells with acquired resistance to TT. Cells were treated with the ceritinib–trametinib combination (500 and 1 nmol/L, respectively, 72 hours) and blotted for the expression of phospho-AKT and phospho-MET. E, Quantification of TYRP1 antigen expression on drug-naïve and TT-resistant SW1 cells. The cells were treated with ceritinib (500 nmol/L), trametinib (1 nmol/L), or combination therapy (72 hours) before being stained with TYRP1 surface antibody. F, Quantification of TYRP1 antigen expression on drug-naïve and TT-resistant SW1 cells. The cells were treated with the ERK inhibitor SCH772984 (500 nmol/L) for 72 hours before being stained with TYRP1 surface antibody. The results are represented as average ± SEM of three independent experiments. Statistical significance was assessed with one-way ANOVA test (*, 0.001 ≤ P ≤ 0.05; **, 0.0001 ≤ P ≤ 0.001; ***, P ≤ 0.0001).

Close modal

Recent years have seen the development of targeted therapies and immunotherapies that have revolutionized the treatment of advanced melanoma. Although there has been much interest in using immunotherapy and targeted therapy in combination, a good mechanistic basis for this strategy has been lacking. In the present study, we focused on studying different immunotherapy and targeted therapy sequences for both NRAS- and BRAF-mutant melanoma. As NRAS-mutant melanoma currently lacks any FDA-approved targeted therapy, we focused on the ceritinib–trametinib combination. This combination, which we previously showed to be effective against multiple melanoma genotypes, is currently being explored in a phase I clinical trial (NCT03501368; ref. 18). In the present study, we demonstrated that the sequence of anti–PD-1 followed by ceritinib–trametinib (IT→TT) was more effective than either treatment alone in the SW1 NRAS-mutant and SM1 BRAF-mutant mouse melanoma models. Similar results were also seen in BRAF-mutant mouse melanoma models following treatment with an IT→dabrafenib–trametinib sequence.

The underlying premise of targeted therapy is based on the inhibition of the growth/survival signals that drive uncontrolled cancer growth. To date, the majority of targeted therapy preclinical studies have been performed in human cancer cell lines grown as xenografts in immunocompromised mice (29–31). Although useful for understanding the tumor-intrinsic effects of these drugs, these approaches fail to adequately account for the role of the tumor-immune microenvironment in the therapeutic responses. There is emerging evidence that TT responses can differ when immune-competent mouse models are used. As one example, BRAF–MEK inhibitor responses are more durable for mouse melanoma cell lines grown in immunocompetent C57/BL6 mice versus immunocompromised NOD/SCID mice. In the immune-competent mice, BRAF–MEK inhibition is associated with vigorous infiltration of CD4+ and CD8+ T cells, and reduced levels of protumorigenic macrophages (32). Analysis of clinical specimens from patients with BRAF-mutant melanoma has also demonstrated that BRAF inhibitor therapy is associated with immune infiltration (9, 33). Other studies using syngeneic models of BRAF-mutant melanoma have shown that concurrent anti–PD-1 therapy improves the durability of responses to BRAF–MEK inhibition (10). This has led to multiple clinical trials combining BRAF–MEK inhibition with anti–PD-1 therapy, which have yielded longer durations of progression-free survival and treatment response compared with BRAF–MEK inhibition alone, at the same time demonstrating significant levels of toxicity (17, 34). There are also trials exploring targeted therapy and immunotherapy in sequence, with some studies waiting until tumor progression for treatment switch (NCT02224781), whereas others focus on using BRAF–MEK inhibitor therapy first (for 2–3 weeks) followed by a switch to anti–PD-1 (NCT03149029).

We here demonstrate that the IT→TT sequence significantly improved responses to targeted therapy by increasing levels of T-cell infiltration, DC accumulation, and upregulation of the tumor antigen presentation machinery. This then sensitized the immune environment to the increase in tumor antigen expression that occurred following the administration of TT. In addition to increasing T-cell numbers, it was noted that the IT→TT sequence increased the numbers of T cells that expressed lower levels of immune checkpoints and exhaustion markers, compared with the numbers observed following IT→Veh. The effects of TT were mediated by both the tumor microenvironment and the tumor and were associated with increased accumulation of multiple myeloid cell types and T cells, and decreased levels of inhibitory immune populations such as Tregs and MDSCs.

An analysis of tumor heterogeneity following IT→TT treatment demonstrated that the initial immunotherapy altered the transcriptional composition of the tumors in vivo, enriching for different immune transcriptional profiles than seen with TT alone. One of the critical changes in gene expression following TT was the increase in melanoma antigen expression. Of note, the IT→TT sequence was also found to enrich for a subcluster of melanoma cells with increased expression of both melanoma antigens and MHC class I (including B2m, H2-D1, and H2-K1). Increased expression of MHC class I was also seen in multiple immune subsets too, suggesting a role for anti–PD-1 therapy in “priming” the immune system for the increase in antigen that followed TT administration (Supplementary Fig. S13). DCs are the major cell type involved in the presentation of tumor antigen to T cells. An analysis of the immune microenvironment demonstrated the IT→TT sequence to be associated with a marked accumulation of a diverse array of DCs. At least three subtypes of DCs were identified in this response, including one subset of activated DC2s that expressed CD80, TLR2, IFNγR2, and CD86 and were not detected in significant quantities following any other therapy sequence. It is likely that this DC2 population played a key role in the sustained CD8+ T-cell responses that were associated with durable therapeutic responses.

One surprise finding was the observation that IT→TT sequence was sufficient to prevent the emergence of a resistance-associated signaling signature in the melanoma cells. One possible explanation for this is that the immune system preferentially recognized the melanoma cells with reduced MEK signaling, due to their increased antigen expression. It is also likely that once the drug-resistant cells emerged, their increased resistance-associated signaling led to downregulation of antigen expression, setting up a vicious cycle of increased tumor growth and immune evasion. Recent work has also suggested that increased tumor burden suppresses systemic immune responses, indicating that once tumors reach a critical size immune function is impaired, setting up a scenario of uncontrolled tumor growth in the face of immune evasion (35). Analysis of TT-resistant SW1 cells supported this idea and showed a reduced inhibition of AKT and MET signaling following drug treatment, along with decreased antigen expression. Analysis of melanoma patient specimens shows that the CD8+ T-cell accumulation that follows BRAF–MEK inhibition declines as a patient fails therapy, and this is accompanied by increased numbers of suppressive immune cells including MDSCs, Tregs, and TAMs (33, 36). In mouse models, depletion of MDSCs sensitizes BRAF inhibitor–resistant tumors to subsequent immune-checkpoint blockade (14). Escape from immunotherapy is also associated with impaired immune responses, characterized by T-cell exhaustion, increased numbers of Tregs, MDSCs, TAMs, and metabolic changes within the tumor microenvironment (37). A potential link between resistance to immunotherapy and targeted therapy has been suggested by transcriptional analyses of tumors from patients who did not respond to anti–PD-1 therapy (38).

Our data suggest that in addition to direct effects upon the tumor, targeted therapy drugs may also improve the durability of the immune response through the suppression of Treg and MDSC accumulation (Supplementary Fig. S13). A number of therapeutic strategies to limit MDSC accumulation have already been proposed, including the use of all-trans retinoic acid to increase MDSC differentiation (39) and CBP/EP300 bromodomain inhibitors to downregulate Arg1 and iNOS in MDSCs (40). There is also evidence that inhibition of TAM-family kinases (Axl, MERTK, TYRO3) can reverse the prooncogenic activities of MDSCs (41). Other studies have demonstrated that ceritinib can suppress MDSC recruitment through the inhibition of CCR2 expression (42). It is therefore likely that the sequential use of immunotherapy and targeted therapy could improve the durability of patient responses. The mechanisms underlying the improved responses to the IT→TT sequence are likely to be complex and go beyond a simple “immune priming” resulting from the initial PD1 dosing. Recent work has shown that PD1, in addition to being expressed on activated/exhausted T cells, is also found on myeloid cells (43, 44). Inhibition or silencing of PD1 in the myeloid compartment leads to marked change in tumor-driven myelopoiesis, characterized by a shift in cell fate away from immature myeloid cells and toward differentiated monocytes and macrophages (43). Our results support these observations and show that treatment with anti–PD-1 induced partial antitumor responses and a large infiltration of macrophages. One of the key changes associated with improved therapeutic responses to the IT→TT combination was the dramatic accumulation of DC1s and DC2s. It is likely that this could have resulted from an amplified immune response, which was induced by the “brake release” provided by the combined effects of anti–PD-1 upon myelopoiesis plus the immunogenic tumor cell death and release of immunostimulatory damage-associated molecular patterns triggered by the targeted therapy (45, 46). Achieving durable responses to targeted therapy seems to be predicated by maintaining a sustained immune response, with the onset of resistance marked by both a recovery of tumor cell growth and immune evasion. Together, our data demonstrate that sustained antitumor responses to targeted therapy are dependent upon a vigorous, sustained immune response and that upfront use of immunotherapy can augment this.

M.A. Davies reports grants and personal fees from Roche/Genentech, GlaxoSmithKline, and Sanofi Aventis and personal fees from Array, Pfizer, Novartis, Bristol-Myers Squibb, Vaccinex, Apexigen, and ABM outside the submitted work. I. Smalley reports grants from NIH during the conduct of the study. D.R. Duckett is a cofounder of Cadw Therapeutics. P.A. Forsyth reports grants from U.S. Department of Defense, CDMRP, NIH/NCI, Pfizer, and State of Florida and other support from AbbVie, Boehringer Ingelheim, NCI, NCRI, NIH, Novellus, Ziopharm, Bayer, BTG, Inovio, and Novocure outside the submitted work. Z. Eroglu reports personal fees from Array, Genentech, OncoSec, Regeneron, SunPharma, and Natera and grants and personal fees from Novartis outside the submitted work. K.S.M. Smalley reports grants from NIH, Florida Department of Health, Lesa Kennedy and Bill Christie Foundation, and U.S. Department of Defense during the conduct of the study, as well as personal fees from Elsevier outside the submitted work. No disclosures were reported by the other authors.

M.S. Phadke: Formal analysis, validation, investigation, writing–original draft. Z. Chen: Data curation, software, formal analysis. J. Li: Data curation, software, formal analysis. E. Mohamed: Investigation, methodology. M.A. Davies: Resources, methodology. I. Smalley: Data curation, formal analysis. D.R. Duckett: Writing–review and editing. V. Palve: Formal analysis, investigation. B.J. Czerniecki: Investigation, writing–review and editing. P.A. Forsyth: Writing–review and editing. D. Noyes: Investigation, methodology. D.O. Adeegbe: Investigation, methodology, writing–review and editing. Z. Eroglu: Writing–review and editing. K.T. Nguyen: Investigation. K.Y. Tsai: Resources, writing–review and editing. U. Rix: Formal analysis, investigation. C.E. Burd: Resources, writing–review and editing. Y.A. Chen: Data curation, software, formal analysis, supervision. P.C. Rodriguez: Formal analysis, methodology, writing–review and editing. K.S.M. Smalley: Conceptualization, resources, supervision, funding acquisition, writing–original draft, project administration.

This work was supported by NIH grants P50 CA168536, R21 CA198550 (to K.S.M. Smalley), and K99 CA226679 (to I. Smalley); U.S. Department of Defense grant W81XWH1810268 (to K.S.M. Smalley); Live Like Bella Grant (9LA03) from the State of Florida (to K.S.M. Smalley); and an Innovative Core Pilot Project funded by a gracious gift from DPR Construction (to Y.A. Chen). The authors thank Lesa Kennedy and Bill Christy for their generous support of the Moffitt Melanoma and Skin Cancer Center of Excellence that significantly contributed to these studies. The Molecular Genomics and Flow Cytometry Core at Moffitt is supported in part by the NCI through a Cancer Center Support Grant (P30-CA076292). M.A. Davies is supported by NIH/NCI (2T32CA009666-21) and the Dr. Miriam and Sheldon G. Adelson Medical Research Foundation. The RPPA analysis was made possible through NCI grant CA16672.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Larkin
J
,
Ascierto
PA
,
Dreno
B
,
Atkinson
V
,
Liszkay
G
,
Maio
M
, et al
Combined vemurafenib and cobimetinib in BRAF-mutated melanoma
.
N Engl J Med
2014
;
371
:
1867
76
.
2.
Wei
SC
,
Duffy
CR
,
Allison
JP
. 
Fundamental mechanisms of immune checkpoint blockade therapy
.
Cancer Discov
2018
;
8
:
1069
86
.
3.
Ott
PA
,
Hodi
FS
,
Robert
C
. 
CTLA-4 and PD-1/PD-L1 blockade: new immunotherapeutic modalities with durable clinical benefit in melanoma patients
.
Clin Cancer Res
2013
;
19
:
5300
9
.
4.
Larkin
J
,
Chiarion-Sileni
V
,
Gonzalez
R
,
Grob
JJ
,
Cowey
CL
,
Lao
CD
, et al
Combined nivolumab and ipilimumab or monotherapy in untreated melanoma
.
N Engl J Med
2015
;
373
:
23
34
.
5.
Sumimoto
H
,
Imabayashi
F
,
Iwata
T
,
Kawakami
Y
. 
The BRAF-MAPK signaling pathway is essential for cancer-immune evasion in human melanoma cells
.
J Exp Med
2006
;
203
:
1651
6
.
6.
Shabaneh
TB
,
Molodtsov
AK
,
Steinberg
SM
,
Zhang
P
,
Torres
GM
,
Mohamed
GA
, et al
Oncogenic BRAF(V600E) governs regulatory T-cell recruitment during melanoma tumorigenesis
.
Cancer Res
2018
;
78
:
5038
49
.
7.
Leslie
C
,
Bowyer
SE
,
White
A
,
Grieu-Iacopetta
F
,
Trevenen
M
,
Iacopetta
B
, et al
FOXP3+ T regulatory lymphocytes in primary melanoma are associated with BRAF mutation but not with response to BRAF inhibitor
.
Pathology
2015
;
47
:
557
63
.
8.
Wilmott
JS
,
Long
GV
,
Howle
JR
,
Haydu
LE
,
Sharma
RN
,
Thompson
JF
, et al
Selective BRAF inhibitors induce marked T-cell infiltration into human metastatic melanoma
.
Clin Cancer Res
2012
;
18
:
1386
94
.
9.
Boni
A
,
Cogdill
AP
,
Dang
P
,
Udayakumar
D
,
Njauw
CN
,
Sloss
CM
, et al
Selective BRAFV600E inhibition enhances T-cell recognition of melanoma without affecting lymphocyte function
.
Cancer Res
2010
;
70
:
5213
9
.
10.
Hu-Lieskovan
S
,
Mok
S
,
Homet Moreno
B
,
Tsoi
J
,
Robert
L
,
Goedert
L
, et al
Improved antitumor activity of immunotherapy with BRAF and MEK inhibitors in BRAFV600E melanoma
.
Sci Transl Med
2015
;
7
:
279ra41
.
11.
Ho
PC
,
Kaech
SM
. 
BRAF-targeted therapy alters the functions of intratumoral CD4(+) T cells to inhibit melanoma progression
.
Oncoimmunology
2014
;
3
:
e29126
.
12.
Ho
PC
,
Meeth
KM
,
Tsui
YC
,
Srivastava
B
,
Bosenberg
MW
,
Kaech
SM
. 
Immune-based antitumor effects of BRAF inhibitors rely on signaling by CD40L and IFNgamma
.
Cancer Res
2014
;
74
:
3205
17
.
13.
Hajek
E
,
Krebs
F
,
Bent
R
,
Haas
K
,
Bast
A
,
Steinmetz
I
, et al
BRAF inhibitors stimulate inflammasome activation and interleukin 1 beta production in dendritic cells
.
Oncotarget
2018
;
9
:
28294
308
.
14.
Steinberg
SM
,
Shabaneh
TB
,
Zhang
P
,
Martyanov
V
,
Li
Z
,
Malik
BT
, et al
Myeloid cells that impair immunotherapy are restored in melanomas with acquired resistance to BRAF inhibitors
.
Cancer Res
2017
;
77
:
1599
610
.
15.
Ribas
A
,
Hodi
FS
,
Callahan
M
,
Konto
C
,
Wolchok
J
. 
Hepatotoxicity with combination of vemurafenib and ipilimumab
.
N Engl J Med
2013
;
368
:
1365
6
.
16.
Ribas
A
,
Lawrence
D
,
Atkinson
V
,
Agarwal
S
,
Miller
WH
 Jr
,
Carlino
MS
, et al
Combined BRAF and MEK inhibition with PD-1 blockade immunotherapy in BRAF-mutant melanoma
.
Nat Med
2019
;
25
:
936
40
.
17.
Ascierto
PA
,
Ferrucci
PF
,
Fisher
R
,
Del Vecchio
M
,
Atkinson
V
,
Schmidt
H
, et al
Dabrafenib, trametinib and pembrolizumab or placebo in BRAF-mutant melanoma
.
Nat Med
2019
;
25
:
941
6
.
18.
Verduzco
D
,
Kuenzi
BM
,
Kinose
F
,
Sondak
VK
,
Eroglu
Z
,
Rix
U
, et al
Ceritinib enhances the efficacy of trametinib in BRAF/NRAS-wild-type melanoma cell lines
.
Mol Cancer Ther
2018
;
17
:
73
83
.
19.
Hennessey
RC
,
Holderbaum
AM
,
Bonilla
A
,
Delaney
C
,
Gillahan
JE
,
Tober
KL
, et al
Ultraviolet radiation accelerates NRas-mutant melanomagenesis: a cooperative effect blocked by sunscreen
.
Pigment Cell Melanoma Res
2017
;
30
:
477
87
.
20.
Burd
CE
,
Liu
W
,
Huynh
MV
,
Waqas
MA
,
Gillahan
JE
,
Clark
KS
, et al
Mutation-specific RAS oncogenicity explains NRAS codon 61 selection in melanoma
.
Cancer Discov
2014
;
4
:
1418
29
.
21.
Butler
A
,
Hoffman
P
,
Smibert
P
,
Papalexi
E
,
Satija
R
. 
Integrating single-cell transcriptomic data across different conditions, technologies, and species
.
Nat Biotechnol
2018
;
36
:
411
.
22.
Aran
D
,
Looney
AP
,
Liu
L
,
Wu
E
,
Fong
V
,
Hsu
A
, et al
Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage
.
Nat Immunol
2019
;
20
:
163
72
.
23.
Pont
F
,
Tosolini
M
,
Fournie
JJ
. 
Single-Cell Signature Explorer for comprehensive visualization of single cell signatures across scRNA-seq datasets
.
Nucleic Acids Res
2019
;
47
:
e133
.
24.
Li
J
,
Smalley
I
,
Schell
MJ
,
Smalley
KSM
,
Chen
YA
. 
SinCHet: a MATLAB toolbox for single cell heterogeneity analysis in cancer
.
Bioinformatics
2017
;
33
:
2951
3
.
25.
Smalley
KS
,
Brafford
P
,
Haass
NK
,
Brandner
JM
,
Brown
E
,
Herlyn
M
. 
Up-regulated expression of zonula occludens protein-1 in human melanoma associates with N-cadherin and contributes to invasion and adhesion
.
Am J Pathol
2005
;
166
:
1541
54
.
26.
Vashisht Gopal
YN
,
Gammon
S
,
Prasad
R
,
Knighton
B
,
Pisaneschi
F
,
Roszik
J
, et al
A novel mitochondrial inhibitor blocks MAPK pathway and overcomes MAPK inhibitor resistance in melanoma
.
Clin Cancer Res
2019
;
25
:
6429
42
.
27.
Rebecca
VW
,
Alicea
GM
,
Paraiso
KH
,
Lawrence
H
,
Gibney
GT
,
Smalley
KS
. 
Vertical inhibition of the MAPK pathway enhances therapeutic responses in NRAS-mutant melanoma
.
Pigment Cell Melanoma Res
2014
;
27
:
1154
8
.
28.
Sullivan
RJ
,
Infante
JR
,
Janku
F
,
Wong
DJL
,
Sosman
JA
,
Keedy
V
, et al
First-in-class ERK1/2 inhibitor ulixertinib (BVD-523) in patients with MAPK mutant advanced solid tumors: results of a phase I dose escalation and expansion study
.
Cancer Discov
2018
;
8
:
184
95
.
29.
Haarberg
HE
,
Paraiso
KH
,
Wood
E
,
Rebecca
VW
,
Sondak
VK
,
Koomen
JM
, et al
Inhibition of Wee1, AKT, and CDK4 underlies the efficacy of the HSP90 inhibitor XL888 in an in vivo model of NRAS-mutant melanoma
.
Mol Cancer Ther
2013
;
12
:
901
12
.
30.
Phadke
MS
,
Sini
P
,
Smalley
KS
. 
The novel ATP-competitive MEK/Aurora kinase inhibitor BI-847325 overcomes acquired BRAF inhibitor resistance through suppression of Mcl-1 and MEK expression
.
Mol Cancer Ther
2015
;
14
:
1354
64
.
31.
Abel
EV
,
Basile
KJ
,
Kugel
CH
 3rd
,
Witkiewicz
AK
,
Le
K
,
Amaravadi
RK
, et al
Melanoma adapts to RAF/MEK inhibitors through FOXD3-mediated upregulation of ERBB3
.
J Clin Invest
2013
;
123
:
2155
68
.
32.
Erkes
DA
,
Cai
W
,
Sanchez
IM
,
Purwin
TJ
,
Rogers
C
,
Field
CO
, et al
Mutant BRAF and MEK inhibitors regulate the tumor immune microenvironment via pyroptosis
.
Cancer Discov
2020
;
10
:
254
69
.
33.
Cooper
ZA
,
Reuben
A
,
Spencer
CN
,
Prieto
PA
,
Austin-Breneman
JL
,
Jiang
H
, et al
Distinct clinical patterns and immune infiltrates are observed at time of progression on targeted therapy versus immune checkpoint blockade for melanoma
.
Oncoimmunology
2016
;
5
:
e1136044
.
34.
Gutzmer
R
,
Stroyakovskiy
D
,
Gogas
H
,
Robert
C
,
Lewis
K
,
Protsenko
S
, et al
Atezolizumab, vemurafenib, and cobimetinib as first-line treatment for unresectable advanced BRAFV600 mutation-positive melanoma (IMspire150): primary analysis of the randomised, double-blind, placebo-controlled, phase 3 trial
.
Lancet North Am Ed
2020
;
395
:
1835
44
.
35.
Allen
BM
,
Hiam
KJ
,
Burnett
CE
,
Venida
A
,
DeBarge
R
,
Tenvooren
I
, et al
Systemic dysfunction and plasticity of the immune macroenvironment in cancer models
.
Nat Med
2020
;
26
:
1125
34
.
36.
Cooper
ZA
,
Juneja
VR
,
Sage
PT
,
Frederick
DT
,
Piris
A
,
Mitra
D
, et al
Response to BRAF inhibition in melanoma is enhanced when combined with immune checkpoint blockade
.
Cancer Immunol Res
2014
;
2
:
643
54
.
37.
Jenkins
RW
,
Barbie
DA
,
Flaherty
KT
. 
Mechanisms of resistance to immune checkpoint inhibitors
.
Br J Cancer
2018
;
118
:
9
16
.
38.
Hugo
W
,
Zaretsky
JM
,
Sun
L
,
Song
C
,
Moreno
BH
,
Hu-Lieskovan
S
, et al
Genomic and transcriptomic features of response to anti-PD-1 therapy in metastatic melanoma
.
Cell
2016
;
165
:
35
44
.
39.
Tobin
RP
,
Jordan
KR
,
Robinson
WA
,
Davis
D
,
Borges
VF
,
Gonzalez
R
, et al
Targeting myeloid-derived suppressor cells using all-trans retinoic acid in melanoma patients treated with Ipilimumab
.
Int Immunopharmacol
2018
;
63
:
282
91
.
40.
de Almeida Nagata
DE
,
Chiang
EY
,
Jhunjhunwala
S
,
Caplazi
P
,
Arumugam
V
,
Modrusan
Z
, et al
Regulation of tumor-associated myeloid cell activity by CBP/EP300 bromodomain modulation of H3K27 acetylation
.
Cell Rep
2019
;
27
:
269
81
.
41.
Holtzhausen
A
,
Harris
W
,
Ubil
E
,
Hunter
DM
,
Zhao
J
,
Zhang
Y
, et al
TAM family receptor kinase inhibition reverses MDSC-mediated suppression and augments anti-PD-1 therapy in melanoma
.
Cancer Immunol Res
2019
;
7
:
1672
86
.
42.
Hu
J
,
Zhang
W
,
Liu
Y
,
Yang
Y
,
Tan
C
,
Wei
X
, et al
LDK378 inhibits the recruitment of myeloid derived suppressor cells to spleen via the p38/GRK2/CCR2 pathway in mice with sepsis
.
Immunol Cell Biol
2019
;
97
:
902
15
.
43.
Strauss
L
,
Mahmoud
MAA
,
Weaver
JD
,
Tijaro-Ovalle
NM
,
Christofides
A
,
Wang
Q
, et al
Targeted deletion of PD-1 in myeloid cells induces antitumor immunity
.
Sci Immunol
2020
;
5
:
eaay1863
.
44.
Gordon
SR
,
Maute
RL
,
Dulken
BW
,
Hutter
G
,
George
BM
,
McCracken
MN
, et al
PD-1 expression by tumour-associated macrophages inhibits phagocytosis and tumour immunity
.
Nature
2017
;
545
:
495
9
.
45.
Galluzzi
L
,
Buque
A
,
Kepp
O
,
Zitvogel
L
,
Kroemer
G
. 
Immunogenic cell death in cancer and infectious disease
.
Nat Rev Immunol
2017
;
17
:
97
111
.
46.
Wculek
SK
,
Cueto
FJ
,
Mujal
AM
,
Melero
I
,
Krummel
MF
,
Sancho
D
. 
Dendritic cells in cancer immunology and immunotherapy
.
Nat Rev Immunol
2020
;
20
:
7
24
.