Combination anti-CTLA-4 plus anti-PD-1 checkpoint blockade utilizes cellular mechanisms partially distinct from monotherapies.

Immune checkpoint blockade therapy targets T cell-negative costimulatory molecules such as cytotoxic T lymphocyte antigen-4 (CTLA-4) and programmed cell death-1 (PD-1). Combination anti-CTLA-4 and anti-PD-1 blockade therapy has enhanced efficacy, but it remains unclear through what mechanisms such effects are mediated. A critical question is whether combination therapy targets and modulates the same T cell populations as monotherapies. Using a mass cytometry-based systems approach, we comprehensively profiled the response of T cell populations to monotherapy and combination anti-CTLA-4 plus anti-PD-1 therapy in syngeneic murine tumors and clinical samples. Most effects of monotherapies were additive in the context of combination therapy; however, multiple combination therapy-specific effects were observed. Highly phenotypically exhausted cluster of differentiation 8 (CD8) T cells expand in frequency following anti-PD-1 monotherapy but not combination therapy, while activated terminally differentiated effector CD8 T cells expand only following combination therapy. Combination therapy also led to further increased frequency of T helper type 1 (Th1)-like CD4 effector T cells even though anti-PD-1 monotherapy is not sufficient to do so. Mass cytometry analyses of peripheral blood from melanoma patients treated with immune checkpoint blockade therapies similarly revealed mostly additive effects on the frequencies of T cell subsets along with unique modulation of terminally differentiated effector CD8 T cells by combination ipilimumab plus nivolumab therapy. Together, these findings indicate that dual blockade of CTLA-4 and PD-1 therapy is sufficient to induce unique cellular responses compared with either monotherapy.

Immune checkpoint blockade can induce long-term durable responses; however, these remarkable responses are limited to a minority of patients (1–3). One of the promising approaches to improve the scope and efficacy of these therapies is combinatorial therapy. There are over 1,100 clinical trials currently testing the efficacy of programmed cell death-1 (PD-1) blockade agents in combination with other therapies (4). These include over 250 combining therapeutic agents targeting cytotoxic T lymphocyte antigen-4 (CTLA-4) and PD-1, as well as combinations with or of novel agents targeting other T cell costimulatory molecules (e.g., TIM-3, LAG3, OX-40). Additionally, combinations of immune checkpoint blockade with other therapeutic modalities such as radiation, chemotherapy, and targeted therapies are under active investigation. These efforts are hampered, however, by our limited understanding of the underlying biological mechanisms. In particular, it remains unclear in many cases what precise molecular events and cellular populations are modulated by checkpoint blockade. It is clear that mechanistic insights are necessary to guide rational design of combination therapies. We sought to define what information is needed to guide such rational design of immunotherapies. It is possible that combination therapies act through modulation of the same cellular populations as constituent monotherapies. Alternatively, combination therapies, particularly those that have convergent effects at the cellular or molecular level, may utilize mechanisms distinct from either monotherapy. Distinguishing between these possibilities is important to define whether mechanistic understanding of monotherapies is sufficient to identify therapeutic combinations of interest or, alternatively, whether each combination must be investigated at a mechanistic level.

We chose to study these possibilities specifically in the context of anti–CTLA-4 and anti–PD-1 checkpoint blockade. This combination is particularly interesting for several reasons and also presents unique opportunities for investigation. First, our current understanding of CTLA-4 and PD-1 biology indicates that therapies targeting these respective molecules can induce distinct biological responses. This is based on differences between the cellular and molecular mechanisms of CTLA-4– and PD-1–mediated regulation of T cell activation (1, 5). In addition, recent evidence indicates that PD-1 inhibits cluster of differentiation 28 (CD28) signaling in addition to proximal signaling elements of the T cell receptor (TCR) complex (6). Given that CTLA-4 primarily attenuates T cell activation by limiting CD28+ costimulation, this raises the possibility that dual blockade may further enhance CD28 signaling. In addition, modulation of the divergent regulation by PD-1 and CTLA-4 of the phosphatidylinositol 3-kinase/AKT pathway and of proximal TCR signaling, respectively, may act synergistically. Finally, because anti–CTLA-4 (ipilimumab), multiple anti–PD-1 therapies, and combination ipilimumab plus nivolumab (anti–CTLA-4 plus anti–PD-1) are all approved by the US Food and Drug Administration for the treatment of advanced melanoma, this offers an opportunity to interrogate and compare the effects of monotherapies and combination therapy in clinical samples.

In preclinical tumor models, anti–CTLA-4 plus anti–PD-1 blockade significantly enhances responses with greatly increased effector-to-suppressor cell ratios and production of proinflammatory cytokines such as interferon-γ (IFN-γ) and tumor necrosis factor-α (7). Combination ipilimumab plus nivolumab has enhanced efficacy compared with monotherapies in advanced melanoma (8), superior efficacy compared with standard-of-care sunitinib in renal cell carcinoma (9), and higher objective response rates in mismatch repair-deficient/microsatellite instability-high colorectal cancer compared with historical response rates of nivolumab monotherapy (10). Follow-up analyses of the clinical trials are indicative of an increase in overall survival compared with either monotherapy. A remarkable 57% 3-y overall survival of metastatic melanoma patients treated with this combination therapy has been reported (11). Several lines of evidence support the hypothesis that combination therapy utilizes, at least in part, distinct mechanisms. The mechanisms of negative costimulatory molecules, and thus checkpoint blockade therapies, can differ at multiple biological scales such as the molecular, cellular, and anatomic levels (5). Combination anti–CTLA-4 plus anti–PD-1 therapy leads to distinct transcriptional effects with ∼400 genes differentially regulated compared with monotherapies in both mouse and human T cells (12, 13). Clinical studies of peripheral blood of melanoma patients treated with checkpoint blockade reveal that combination blockade of PD-1 and CTLA-4 is associated with specific cellular immune responses such as dramatic enhancement of cytokine production, increases in T cell frequency, and modulation of peripheral B cell populations (13, 14). Previous studies have also revealed that monotherapies targeting PD-1 and CTLA-4 regulate distinct T cell subsets (15), with the main difference being that anti–CTLA-4, but not anti–PD-1, treatment expands T helper type 1 (Th1)-like CD4 effector T cells. These observations raise the possibility that concurrent loss of both CTLA-4 and PD-1 regulation due to antibody blockade may permit noncanonical signaling events. We hypothesized that such dysregulation due to dual blockade would manifest as distinct effects on antitumor T cell populations compared with monotherapies.

Results

Comprehensive Profiling of Murine Tumor Immune Infiltrates in the Context of Combination Immune Checkpoint Blockade Therapy.

To identify the specific tumor-infiltrating T cell populations modulated by anti–CTLA-4 plus anti–PD-1 combination therapy, we utilized a mass cytometry-based systems approach to analyze MC38 tumors from mice treated with control isotypes, anti–CTLA-4, anti–PD-1, or combination anti–CTLA-4 plus anti–PD-1 antibodies. MC38 is highly immunogenic and sensitive to both monotherapy and combination immune checkpoint blockade, consistent with induction of robust antitumor immune responses (Fig. 1). As such, this model system can be used to distinguish effects of monotherapy and combination therapy, independent of effects that differentiate effective and ineffective therapies. Although it is likely that combination therapy will provide the most added benefit in less immunogenic tumor types, the effects on T cell populations due to therapy observed in this model are likely to be generalizable, given prior observations suggesting that the type (but not magnitude) of T cell responses is tumor type-independent (15). MC38 tumor infiltrates were then computationally analyzed using PhenoGraph, a well-validated unsupervised clustering approach for the robust classification of immune cell populations (15–20). Changes in the relative frequency of both lymphoid and myeloid populations were observed in response to therapy (); however, we primarily focused on the T cell compartment, given current understanding of the mechanisms of action of anti–CTLA-4 and anti–PD-1 therapies. The frequency of total tumor-infiltrating T cells increased following both monotherapy and combination therapy (Fig. 1); however, we sought to distinguish effects on the relative frequencies of specific T cell subsets. Although relative frequencies do not distinguish in situ expansion and infiltration or necessarily directly correlate with total cell number, such analyses can provide considerable biological insight into the nature of antitumor T cell responses.

Fig. 1.

Combination anti–CTLA-4 plus anti–PD-1 therapy differentially affects MC38 tumor-infiltrating T cell populations. (A) MC38 tumor volumes 13 d after tumor inoculation are displayed from mice treated with control, monotherapy, or combination anti–CTLA-4 plus anti–PD-1 therapy. (B) T cell frequency as a percentage of total MC38 tumor-infiltrating CD45+ cells. (C) Heat map of MC38 tumor-infiltrating T cell metaclusters displaying expression values of individual parameters normalized to the maximum mean value across metaclusters. (D) Relative frequency of CD8 T cell metaclusters. (E) Relative frequency of CD4 effector T cell metaclusters. (F) Relative frequency of Treg metaclusters. (G) Mean intensity of CTLA-4 of Treg metaclusters. Each data point reflects the mean intensity of the individual clusters associated with the metaclusters. Data are pooled from 3 biological replicate cohorts and are displayed on a per mouse basis. Mean and SD are displayed. *P < 0.05, 2-tailed t test with Welch’s correction. (See also .)

To robustly identify tumor-infiltrating T cell populations across multiple replicate cohorts of mice, we utilized a PhenoGraph-driven metaclustering approach (). Briefly, metaclustering group subsets were identified in individual mice to enable robust classification across multiple mice and cohorts, while minimizing technical batch effects. Using this approach, 14 MC38 tumor-infiltrating T cell subsets were identified, including 3 regulatory T cell (Treg), 3 CD4 effector, and 4 CD8 T cell populations (Fig. 1). The effects of anti–CTLA-4 and anti–PD-1 monotherapy on the relative frequencies of T cell subsets were consistent with prior findings using similar approaches (15) (Fig. 1 ). Interestingly, the frequency of multiple T cell subsets was significantly different in mice treated with combination anti–CTLA-4 plus anti–PD-1 therapy compared with those receiving monotherapy. In the CD8 compartment, phenotypically exhausted (PD-1high Lag3+ Tim3+) CD8 T cells expanded in response to anti–PD-1 monotherapy but decreased in frequency following combination therapy (Fig. 1, metacluster 3). In contrast, activated terminally differentiated effector CD8 T cells (PD-1+ Lag3int Tim3int; hereafter referred to as activated effector CD8 T cells) did not increase after anti–PD-1 monotherapy but increased significantly following combination therapy (Fig. 1, metacluster 1).

In the CD4 effector T cell compartment, a Th1-like PD-1+ ICOSint TBET+ CD4 effector T cell population expanded following anti–CTLA-4 but not anti–PD-1 monotherapy (Fig. 1, metacluster 2). Interestingly, combination therapy led to a further increase in frequency compared with anti–CTLA-4 monotherapy. This suggests that in the context of anti–CTLA-4 monotherapy, Th1-like CD4 effector T cells increase in frequency but also up-regulate PD-1, which limits further expansion. Thus, the effect of PD-1 blockade on this population appears to be conditional on initial engagement following anti–CTLA-4 therapy.

In the Treg compartment, both anti–CTLA-4 and anti–PD-1 monotherapies led to decreased relative frequencies of specific subsets and combination therapy led to additive reductions in Treg frequencies (Fig. 1). These decreases likely constitute both Fc-mediated cell depletion and relative dilution, as anti–CTLA-4 treatment leads to Treg depletion in murine tumor models (21–23). Here, the anti–CTLA-4 antibody clone 9H10, which can mediate Treg depletion in murine tumor models, was utilized in both monotherapy and combination therapy. Similar experiments in murine tumor models using nondepleting anti–CTLA-4 antibodies may help to further define the contributions of direct effects on T effector cells and Treg depletion. Notably, however, while specific anti–CTLA-4 antibody clones can mediate Treg depletion in murine tumor models, recent evidence indicates that ipilimumab therapy does not lead to Treg depletion in humans (24). Interestingly, however, the frequency of only specific Treg populations decreased following both monotherapy and combination therapy (Fig. 1). This may be explained by differential sensitivity of these populations to antibody-dependent cellular cytotoxicity (ADCC) efficiency or, alternatively, differential effects on the proliferation of specific Treg populations due to checkpoint blockade. Notably, these are not mutually exclusive possibilities but do reflect fundamentally different mechanistic explanations.

We investigated whether different Treg populations expressed varying levels of CTLA-4 as a potential explanation for varying sensitivities to blockade. Treg metacluster 8 expressed the highest levels of CTLA-4 and displayed the most dramatic decrease in frequency following CTLA-4 blockade (Fig. 1 ). Interestingly, however, Treg metacluster 9 expressed significant but slightly lower levels of CTLA-4 than metacluster 8, yet did not decrease in frequency following either anti–CTLA-4 or combination therapy. Furthermore, Treg metacluster 12 did not express detectable levels of CTLA-4 but did decrease in frequency following anti–CTLA-4 and combination therapies (Fig. 1 ). This suggests that modulation of metacluster 12 is mediated by indirect effects or by low amounts of CTLA-4 below the limit of detection.

Together, these data reveal that relative to monotherapies, the effects of combination therapy fall into discrete groups, including opposing effects on phenotypically exhausted CD8 T cells, additive expansion of Th1-like CD4 effector T cells, and decreases in specific types of Tregs. Similar effects on the relative frequency of T cell subsets as a frequency of total CD45+ tumor-infiltrating cells support the notion that these changes reflect increases in population, rather than simply relative changes contained only within the T cell compartment (). In addition, these data raise the possibility that combination therapy may have unique biological effects on T cell function. Notably, the analysis of total T cells was not sufficient to detect differences between combination and monotherapies (Fig. 1), while unsupervised phenotypic classification enabled the detection of relevant biological differences (Fig. 1 ). From a mechanistic perspective, such subsets can arise due to modulated proliferation, but also changes in differentiation. Consistent with this possibility, negative costimulation restricts the boundaries of T cell phenotypes achieved through peripheral differentiation (25). In addition, checkpoint blockade can lead to proliferative expansion of specific T cell subsets, but these cells can differentiate and display altered phenotypes when detected at later times. For example, PD-1 blockade leads to expansion of TCF1+ CD8 T cells with memory and stem-like properties, which contributes to the increased frequency of differentiated PD-1+ CD8 T cells in the periphery (26–28).

The Frequencies of Specific Tumor-Infiltrating T Cell Populations Correlate with Response in the Contexts of Monotherapy and Combination Checkpoint Blockade Therapy.

We next sought to gain insight into whether differential effects on T cell subset frequency due to monotherapy and combination therapy reflect functional changes with respect to tumor rejection. We assessed the correlation between individual populations and tumor volume. In the case of most T cell subsets, findings were similar to prior observations using similar methods (15) (Fig. 2 and ). For example, the frequency of Th1-like CD4 effector T cells negatively correlated with tumor volume, consistent with prior findings and a functional role in tumor rejection (15) (Fig. 2, metacluster 2). It is unclear whether combination therapy differentially modulates the activity of this subset compared with anti–CTLA-4 monotherapy (which also leads to an increase in frequency). The functional contribution of this population is not well defined, but may contribute to tumor rejection through several mechanisms such as enhancing CD8 T cell infiltration, increasing antigen presentation, or direct tumor lysis.

Fig. 2.

Correlative relationships between MC38 tumor growth and T cell subset frequency. (A) Relative frequency of denoted MC38 tumor-infiltrating T cell metaclusters plotted as a function of final tumor volume on day 13 postinoculation. Data from all treatment groups are displayed together. Data are pooled from the 3 biological replicate cohorts analyzed in Fig. 1 and are displayed on a per mouse basis. Data points are color-coded by treatment: control (blue), anti–CTLA-4 (red), anti–PD-1 (green), and combination (purple). n.s., not significant. (B) Relative frequency of Th1-like CD4 effector cells (metacluster 2) plotted as a function of final tumor volume. Data from all treatment groups are displayed together. (C) Relative frequency of regulatory T cells (metacluster 8) plotted as a function of final tumor volume. Data from all treatment groups are displayed together. (D) Frequency of IdU+ T cells is displayed for each MC38 tumor-infiltrating T cell metacluster on a per mouse basis. The mean and SD are displayed. Linear regression lines are displayed along with Spearman r and P values. (See also .)

The frequency of Treg metacluster 8 positively correlated with tumor volume, consistent with a protumor functional role (Fig. 2). The frequency of activated effector CD8 T cells (metacluster 1) negatively correlated with tumor volume (Fig. 2). Thus, this population has a functional role independent of therapy, but only combination therapy leads to an increase in frequency. All T cell populations had a similar correlative relationship with tumor volume in comparison to prior findings in the context of monotherapies, with the notable exception of metacluster 3. The frequency of phenotypically exhausted PD-1high Lag3+ Tim3+ CD8 T cells (metacluster 3) did not significantly correlate with tumor volume. At first, this appears to contrast with the notion that exhausted CD8 T cells retain functional capacity as well as prior findings that this population negatively correlates with tumor volume in preclinical tumor models (15). However, upon further inspection, the correlations between the frequency of metacluster 3 and tumor volume diverge in the contexts of monotherapies and combination therapy (). While a strong negative correlation was observed in the context of monotherapy, consistent with prior observations, a positive correlation was observed in the context of combination therapy. This suggests that metacluster 3 has distinct effects on antitumor immune responses in the context of monotherapy and combination therapy.

The striking differences in the frequencies and functional associations of phenotypically exhausted and activated effector CD8 T cells (metaclusters 3 and 1, respectively) in the context of monotherapy versus combination therapy raise several interesting nonmutually exclusive possibilities. First, combination therapy may be able to prevent CD8 T cells from acquiring a fully exhausted phenotype. It is also possible that combination therapy may be sufficient to reverse T cell exhaustion, although this is unlikely, given the kinetics of these experiments and recent evidence that epigenetic programs enforce an exhausted cell state (29, 30). Second, the enhanced CD28 costimulation enabled by CTLA-4 blockade may enhance priming and repriming of newly activated as well as already activated T cells, delaying or perhaps preventing induction of T cell exhaustion. Finally, combination therapy may have distinct effects on the proliferative capacity of specific T cell populations that manifest as differences in relative frequency.

We next investigated whether differences in the efficacy of therapy (and thus the “T cell inflammation status”) could explain the apparent switch from expansion of phenotypically exhausted cells to expansion of newly activated cells. Taking advantage of the biological variance in this model system, we compared approximately size-matched tumors of monotherapy and combination therapy groups, observing similar effects on CD8 T cell populations in response to monotherapies and combination therapy (). This suggests that differences in efficacy do not underlie the differential effects on CD8 T cell metacluster 1 and 3 subsets, although it remains possible that the kinetics of antitumor responses may be regulated in a treatment-dependent manner. To investigate whether monotherapies and combination therapy differentially modulate proliferation of tumor-infiltrating T cell populations, we assessed short-term iodo-deoxyuridine (IdU) incorporation to identify differences in cellular proliferation between tumor-infiltrating T cell subsets (). Consistent with previous findings (15), tumor-infiltrating T cell subsets displayed a wide range of proliferative activity (Fig. 2). No differences in the proliferation of phenotypically exhausted CD8 T cells were observed between treatment groups (), which indicates that the decrease in frequency in response to combination therapy is not due to modulation of proliferation. Interestingly, despite displaying a highly exhausted phenotype, nearly 20% of PD-1high Lag3+ Tim3+ CD8 T cells (metacluster 3) were actively proliferating regardless of therapy. These data support a model in which phenotypically exhausted CD8 T cells still functionally contribute to antitumor responses (Fig. 2 and ).

The shift from phenotypically exhausted to activated effector CD8 T cell effects following combination therapy may be exerted through multiple mechanisms, including modulation of T cell trafficking to the tumor, altered T cell survival, or preventing the adoption of a fully exhausted state. Inhibition of lymphocyte egress from secondary lymphoid organs leads to a loss of therapeutic efficacy in the context of anti–PD-1 with intratumoral adjuvants (31), consistent with the necessity of engagement of peripheral T cells in secondary lymphoid organs for therapeutic efficacy. Notably, however, inhibition of lymphocyte egress attenuates responses to immune checkpoint blockade combination therapies when performed during early phases of responses to therapy but not when performed later (32). This suggests that combination immune checkpoint blockade enhances T cell infiltration and also modulates the functional responses of T cells already within the tumor microenvironment (TME).

Comparison of T Cell Populations in Peripheral Blood from Melanoma Patients Treated with Monotherapies versus Ipilimumab Plus Nivolumab Combination Therapy.

We next sought to understand whether such mechanistic differences underlie responses to monotherapies and combination immune checkpoint blockade therapy in humans. Prior studies indicate that significant changes in the frequencies of CD8 T cell and other immune cell populations can be observed in peripheral blood of patients receiving anti–PD-1 therapy (33–35). Identifying peripheral blood immune cell populations that are modulated by monotherapy and combination therapy might provide candidate biomarkers of response or pharmacokinetics, regardless of whether such phenotypes are present within the tumor. Based on our preclinical studies, we sought to determine whether analyses of peripheral blood could provide insights into the cellular effects of combination therapy. We analyzed peripheral blood samples from a retrospective cohort of metastatic melanoma patients treated with ipilimumab, anti–PD-1 (nivolumab or pembrolizumab), or combination ipilimumab plus nivolumab (). A small number of normal donor peripheral blood samples were also analyzed in each batch as reference samples to enable the assessment and correction of technical batch effects ( and ).

First, we investigated whether immune checkpoint blockade therapies preferentially engage already activated T cells. Manual gating revealed a strong enrichment of Ki-67+ proliferative cells in PD-1+ versus PD-1− CD8 T cells following anti–PD-1 monotherapy (Fig. 3), consistent with prior observations (33, 34). This aligns with the notion that only already activated T cells expressing PD-1 would be sensitive to anti–PD-1 therapy. Interestingly, however, a similar effect was observed in patients who had received ipilimumab monotherapy (Fig. 3). It is unclear whether this effect of anti–CTLA-4 therapy is due to modulation of already activated PD-1–expressing T cells or activation of resting T cells that subsequently up-regulate PD-1. Most notably, the frequencies of Ki-67+ cells were comparably elevated in both PD-1− and PD-1+ CD8 T cells following combination ipilimumab plus nivolumab therapy, consistent with our preclinical findings (Fig. 3). This suggests that monotherapies are sufficient to induce the proliferation of already activated PD-1+ CD8 T cells and combination therapy is able to induce the proliferation of PD-1+ and PD-1low CD8 T cells. In addition, manual gating revealed trends in the frequency of total CD4 and CD8 T cells consistent with engagement of CD8 and CD4 by anti–PD-1 and anti–CTLA-4, respectively (). Combination therapy led to an increase in the overall frequency of CD8 T cells compared with anti–PD-1 monotherapies. In contrast, the frequency of CD4 T cells was increased in the blood of patients treated with ipilimumab compared with those treated with anti–PD-1 therapy. Interestingly, the frequency of Ki-67+ CD4 T cells in peripheral blood was not different between treatments, in contrast to the observed effects on CD8 T cells ().

Fig. 3.

Mass cytometry analysis of peripheral blood from patients treated with checkpoint blockade reveals treatment-specific effects. (A) Frequency of Ki-67+ cells within PD-1− and PD-1+ fractions of CD8 T cells in peripheral blood from patients treated with checkpoint blockade therapies (anti–CTLA-4, ipilimumab [Ipi], anti–PD-1 monotherapy, and ipilimumab plus nivolumab [Ipi + Nivo]) and normal donors is displayed. *P < 0.05, 1-way ANOVA with Sidak’s multiple testing correction. n.s., not significant. (B) tSNE plot of T cells from peripheral blood from patients treated with checkpoint blockade therapies. The plot is color-coded by treatment. Equal numbers of events per treatment group are displayed. (C) tSNE plots of T cells from peripheral blood overlaid with the expression levels of CD4, CD8, and TBET as heat maps. (D) Frequency of CD8 T cell metaclusters identified in peripheral blood from normal donors and patients treated with monotherapy and combination checkpoint blockade therapy. *P < 0.05 Tukey’s 2-way ANOVA with multiple testing correction. The mean and SD are displayed for each frequency plot. ND, normal donor. (E) Frequency of CD4 effector T cell metaclusters identified in peripheral blood from normal donors and patients treated with monotherapy and combination checkpoint blockade therapy. *P < 0.05, Tukey’s 2-way ANOVA with multiple testing correction. (See also .)

To precisely delineate T cell populations that respond to therapy, a similar metaclustering approach was utilized. PD-1 and CTLA-4 were not used as parameters for clustering to exclude the possibility that staining can be inhibited by receptor occupancy following therapy, as has been observed previously in the context of anti–PD-1 therapy (33). Twenty-four T cell metaclusters were identified in peripheral blood samples, including 8 CD8 and 15 CD4 T cell subsets (, Fig. 3 , and ). Although anti–CTLA-4 treatment can induce ADCC in murine models, treatment with ipilimumab did not appear to lead to decreased Treg relative frequencies (Fig. 3, metacluster 7). This may be attributed to differences in the effects on peripheral blood, but is consistent with a recent report that ipilimumab does not deplete Tregs (24). How these findings can be reconciled with prior findings that high-affinity Fc receptor allelic variants are associated with improved responses to anti–CTLA-4 therapy remains to be clarified (36). Recent findings that engagement of anti–CTLA-4 antibodies by Fcγ receptors can modulate T cell responses to therapy independent of Tregs may explain these findings in part (37).

Treatment-dependent associations with subset frequencies were observed in the CD8 and CD4 effector compartments. The frequency of a terminally differentiated TBET+ EOMES+ CD8 T cell subset was significantly increased in blood samples from patients receiving combination therapy compared with those receiving monotherapies (Fig. 3, metacluster 1). Furthermore, metacluster 1 is more frequent in posttherapy samples versus in unmatched pretherapy samples (). Surprisingly, no significant changes in the frequency of PD-1+ CD8 T subsets (e.g., metacluster 16) were observed. These observations suggest that while monotherapies are sufficient to increase the proliferation of activated CD8 T cells (Fig. 3), combination therapy is sufficient to also increase the frequency of terminally differentiated effector CD8 T cells. Supportive of this, monotherapies only led to increased Ki-67+ frequencies in PD-1+ CD8 T cells, whereas combination therapy led to similar Ki-67+ frequencies in PD-1+ and PD-1− CD8 T cells (Fig. 3).

Three CD4 effector populations were elevated following therapies containing anti–CTLA-4. The frequencies of activated CD45RO+ and CD45RO+ TBET+ EOMES+ CD4 effector T cell subsets were higher following ipilimumab compared with anti–PD-1 monotherapy (Fig. 3, metaclusters 3 and 8, respectively). In addition, the frequency of a CD45RO+ ICOS+ PD-1+ CD4 effector T cell subset was higher in patients receiving ipilimumab than in patients receiving anti–PD-1 monotherapy or combination therapy (Fig. 3, metacluster 19). This observation is consistent with prior observations that anti–CTLA-4 therapy leads to expansion of ICOS+ CD4 effector T cells (38–41); however, it remains unclear why a similar elevated frequency is not observed in the peripheral blood of combination therapy-treated patients.

Systematic Characterization of the Relationship Between T Cell Populations in the TME versus Those in Peripheral Blood.

A key question is how closely changes in cellular frequencies in peripheral blood reflect changes in T cell frequencies within the TME. It is generally recognized that the phenotypes of tumor-infiltrating and peripheral blood lymphocyte populations are quite different; however, the extent of phenotypic and functional overlap remains ill-defined. We sought to apply a mass cytometry-based approach to systematically characterize relationships between the phenotypes and frequencies of such T cell populations in each tissue. To investigate the degree to which phenotypically defined T cell populations in peripheral blood reflect intratumoral T cell populations, we analyzed a publicly available mass cytometry dataset of matched lung tumors, normal adjacent lungs, and peripheral blood (20). This matched analysis approach enables assessment of the relative anatomical distribution of phenotypically defined T cell populations. A total of 21 T cell metaclusters were identified with 8 CD8 and 9 CD4 T cell populations (Fig. 4).

Fig. 4.

Specific subsets of phenotypically defined T cell populations within the TME are detectable in peripheral blood of treatment-naive patients. (A) Heat map of T cell metaclusters identified across matched human lung tumors, normal adjacent lung, and peripheral blood. Expression values of individual parameters are normalized to the maximum mean value across metaclusters. (B) Mean frequency of T cell metaclusters in blood plotted as a function of their respective frequencies in tumor tissue. Interpatient variability is displayed within each tissue as the SE. Metacluster clusters are annotated. (C) Mean frequency of T cell metaclusters in normal adjacent lung plotted as a function of their respective frequencies in tumor tissue. (See also .)

Although many phenotypically defined T cell populations were identified in all tissues, several metaclusters were tissue-specific (Fig. 4 ). Tumor-infiltrating T cell frequencies were more closely associated with T cell frequencies in normal adjacent lung than in peripheral blood. Although this observation is expected, it highlights a notable caveat to the identification and interpretation of tumor-infiltrating T cell populations. Although particular T cell populations may be present in tumor and not blood, they may not necessarily be directly relevant to the antitumor response but rather represent “bystander” tissue-resident populations. For example, the frequency of metacluster 7 in the tumor correlates with that in the normal adjacent lung, while it does not correlate with that in blood (and is present at low frequencies) (). While resident T cell populations may play a significant role in antitumor immunity, this observation also raises the possibility that lung tissue-resident T cell populations may be present or even expand within the TME as “passengers” due to indirect effects.

Of interest were 1 CD8 and 2 CD4 effector populations that were present at significant frequencies in tumor but undetectable in blood (Fig. 4 , metaclusters 8, 14, and 18). The phenotypes of these subsets suggested active involvement in the antitumor response. The CD8 T cell subset displayed a PD-1+ LAG3int TIM3int phenotype (Fig. 4, metacluster 18). This metacluster also displayed expression of CD39, which has been reported to be a marker of exhausted CD8 T cells in humans (42). Of the 2 CD4 T cell populations, metacluster 8 displayed an activated PD-1hi ICOSint CD39int phenotype, while metacluster 14 displayed a PD-1int CD127+ phenotype (Fig. 4 ). These observations are consistent with the notion that T cells are exposed to chronic antigen stimulation within the TME and suggest that specific aspects of endogenous antitumor T cell responses can only be detected in tumor tissue. The presence of metacluster 8 but not metacluster 18 in normal adjacent lung tissue raises the possibility that metacluster 8 T cells are tissue resident-derived, while metacluster 18 arises from a systemic immune response. In addition, because these analyses are of treatment-naive patients, the extent to which therapy-responsive populations can be detected remains unclear. These analyses indicate that while the bulk of T cell phenotypes are detectable in both peripheral blood and tumor samples, only a minority of subsets have a correlation between their frequencies in tumor and blood. These analyses also identify several populations for which frequencies in blood reflect intratumoral frequencies. These include an activated antigen D-related human leukocyte antigen-positive (HLA-DR+) CXCR3+ CD4 effector population and a CD45RAint CD8 subset (Fig. 4, metaclusters 1 and 3, respectively). Together, these findings support the utility of peripheral blood analyses, albeit limited in scope in treatment-naive patients, but also indicate that analyses of peripheral blood alone cannot provide a complete characterization of antitumor immune responses.

Discussion

Here, we sought to determine whether combination anti–CTLA-4 plus anti–PD-1 therapy utilizes similar or distinct cellular mechanisms compared with monotherapies. These studies indicate that combination therapy effects on the frequencies of tumor-infiltrating and peripheral T cell populations are, in large part, similar to those of the constituent monotherapies, but also reveal distinct effects on phenotypically exhausted CD8 T cell populations. This indicates that dual blockade of CTLA-4 and PD-1 allows for engagement of divergent biological pathways in CD8 T cells. These findings are consistent with recent evidence that combination anti–CTLA-4 plus anti–PD-1 therapy leads to an expansion of both CD4 and CD8 CD45RO+ EOMES+ CD69+ T cells within the tumor (43). Importantly, our findings further define differential effects of combination therapy on phenotypically exhausted PD-1high Lag3+ Tim3+ and activated effector PD-1+ Lag3int Tim3int CD8 T cells. These findings are in addition to additive effects of combination therapy on multiple CD4 and CD8 T cell subsets, including Th1-like CD4 effector T cells.

Overall, these findings provide insight into what mechanistic insights are required to guide rational design of combination immunotherapies. These findings suggest it is necessary to investigate the mechanisms of action of combination therapies rather than indirectly inferring such biology from the effects of monotherapies. Likewise, the known pharmacodynamics and biomarkers of response for monotherapies may not accurately reflect those of combination therapy. Larger and prospective clinical studies, as well as broader preclinical investigation, are warranted to further investigate these possibilities. Such mechanistic differences may make combination therapy more efficacious in the less immunogenic tumor types refractory to monotherapies. The enhanced efficacy of combination therapy in the context of renal cell carcinoma versus sunitinib supports this notion (9); however, the relative contribution of ipilimumab to the response rate in this context remains ill-defined, given emerging evidence that anti–PD-1 monotherapy has a relatively high response rate in the setting of renal cell carcinoma (44, 45).

A key outstanding question is what specific signaling pathways are differentially modulated by combination anti–CTLA-4 plus anti–PD-1 therapy in order to lead to the increase in frequency of activated effector PD-1+ Lag3int Tim3int CD8 T cells. Prior transcriptional analyses of murine tumor-specific CD8 T cells revealed that combination anti–CTLA-4 plus anti–PD-1 treatment leads to regulation of genes associated with effector function (12). Although the brief timing of our experimental model likely does not allow for the induction of a completely exhausted T cell state, as defined by adoption of exhaustion markers and decrease of functional output, the differential modulation of activated and phenotypically exhausted CD8 T cells by combination therapy raises several interesting possibilities. For example, combination therapy may be sufficient to attenuate the induction of a T cell exhaustion phenotype and/or associated reductions in functional output. Prior studies, however, reveal that the epigenetic state of exhausted CD8 T cells prevents a stable reversion by anti–PD-1 monotherapy (29, 30). If combination therapy can prevent, attenuate, or potentially reverse the induction of T cell exhaustion, this would have significant therapeutic implications.

Another related open question is which functional aspects of the antitumor response are modulated by combination therapy. Specifically, it is unclear whether combined checkpoint blockade enhances the magnitude of responses (e.g., increased killing capacity of CD8 effector T cells), alters the nature of antitumor immune responses (e.g., types of cytokines), or both, relative to checkpoint blockade monotherapies. Our studies here indicate that anti–CTLA-4 and combination therapy, but not anti–PD-1 monotherapy, leads to an expansion of CD4 effector T cells. This does not indicate that CD4 effector T cells are unnecessary for effective antitumor T cell responses in the context of PD-1 blockade alone, but that PD-1 blockade does not enhance the endogenous CD4 effector response. The direct functional consequence of Th1-like CD4 effector T cell expansion following anti–CTLA-4 monotherapy and combined anti–CTLA-4 and anti–PD-1 therapy remains unclear. Recent evidence indicates that enhancement of CD4 effector function by combination blockade leads to increased activation of CD103+ dendritic cells (46). In addition, based on insights from viral systems, Th1 CD4 helper T cells may enhance CD8 T cell infiltration (47), enhance antibody penetration (48), enhance T cell memory formation, or have direct cytolytic capability as has been previously observed with tumor-specific CD4 T cells in the context of adoptive cell transfer (49). Another intriguing possibility is that CD4 T cells may employ cytokine-mediated killing of tumor cells rather than engaging in direct cytolysis. In addition, recent findings reveal that immune checkpoint blockade induces profound changes in specific myeloid subsets. These changes are most striking in the context of combination anti–CTLA-4 plus anti–PD-1 therapy and are partially dependent on IFN-γ (50). This raises the possibility that indirect modulation of tumor-infiltrating myeloid subsets by IFN-γ–producing CD8 and Th1-like CD4 effector T cells may represent part of the mechanism of combination therapy. Consistent with this notion, loss of negative costimulation releases limits on T cell differentiation, with genetic loss of CTLA-4 leading to emergence of noncanonical CD4 T cell phenotypes with elevated polyfunctional cytokine expression, including IFN-γ (25).

Understanding whether combination therapies utilize similar mechanisms as monotherapies is further important for determining whether combination therapies exhibit drug independence. Prior reports have suggested that most drug combinations, including ipilimumab plus nivolumab, derive their benefit through drug independence rather than additivity (51). These analyses were based on short-term follow-up of phase III clinical trial data, and further investigation in the context of long-term follow-up is clearly warranted. The contributions of drug independence and drug additivity are not completely mutually exclusive possibilities, and may both contribute to the enhanced efficacy of combination therapy. These possibilities highlight the need for further mechanistic investigation of monotherapies and combination therapies in both preclinical models and clinical samples. In summary, our findings here indicate that combination immune checkpoint blockade therapy induces distinct biological effects that manifest in different cellular changes in T cell populations.

Methods

Mice and In Vivo Murine Tumor Experiments.

Seven-week-old female C57BL/6J mice were purchased from The Jackson Laboratory and allowed to acclimate in the housing facility for at least 1 wk prior to the initiation of tumor experiments. A previously empirically defined experimental system was utilized for MC38 tumor experiments (15). Briefly, 4 × 105 MC38 tumor cells suspended in 100 μL of phosphate-buffered saline (PBS) were injected subcutaneously on the flank of mice. Immune checkpoint blockade antibodies (polyclonal hamster and rat IgG2a isotype [BioXCell], anti–CTLA-4 [clone 9H10; BioXCell], anti–PD-1 [RMP1-14; BioXCell], or combination anti–CTLA-4 plus anti–PD-1) were injected intraperitoneally 5, 8, and 11 d following tumor injection. Hamster polyclonal and anti–CTLA-4 antibodies were administered using doses of 200 μg for the initial injection and 100 μg for subsequent injections. Anti–PD-1 and rat isotype antibodies were administered at a dose of 250 μg per injection. Mice were randomized to treatment groups prior to the first treatment dose. For short-term IdU experiments, 100 μL of 5 mg/mL IdU (Sigma–Aldrich) in PBS was administered by intraperitoneal injection 30 min prior to euthanasia. Tumors were measured by digital caliper 2 to 3 times per week, with tumor volume calculated as volume = length × width × height. Mice were euthanized using CO2 and subsequent cervical dislocation at the time of euthanasia for analysis. All animal experiments were performed in accordance with The University of Texas MD Anderson Cancer Center (MDACC) Institutional Animal Care and Use Committee (IACUC) guidelines in an MDACC American Association for Accreditation of Laboratory Animal Care (AAALAC)-accredited barrier facility vivarium.

Cell Lines.

MC38 was grown as previously described (15, 52). Briefly, MC38 murine colon carcinoma cells were originally obtained from N. Restifo, National Cancer Institute, Bethesda, MD; cultured in Dulbecco’s modified Eagle’s medium supplemented with 10% fetal bovine serum and penicillin and streptomycin; and periodically tested for Mycobacterium contamination. MC38 was derived from a female C57BL6 mouse. Cell lines were previously analyzed using whole-exome sequencing to interrogate mutational load (15), but have not been further authenticated by other approaches.

Human Subjects.

Peripheral blood samples were from patients treated at The University of Texas MDACC between December 2011 and May 2017. All samples were obtained with patient informed consent, deidentified, and then analyzed under The University of Texas MDACC Institutional Review Board-approved protocols and in accordance with the Declaration of Helsinki. Clinical annotation data are displayed in . Normal donor blood samples were collected from the Gulf Coast Blood Bank. Protected health information (PHI) is protected under The University of Texas MDACC PHI policy ADM0396. Patient age and gender information are not disclosed accordingly. Because normal donor blood samples were primarily used as reference samples (biological and technical) to ensure consistency across multiple acquisition batches, the donor samples are the same as in prior cohorts analyzed (15).

Mass Cytometry Antibodies.

Metal-conjugated antibodies were either purchased from Fluidigm or conjugated in-house using X8 polymer kits per the manufacturer’s protocol as previously described (15). Antibody, metal conjugate, and clone information is detailed in .

Mass Cytometry Sample Preparation and Data Acquisition.

MC38 murine tumors were isolated from mice treated with control or checkpoint blockade antibodies and killed 13 d after tumor inoculation. Tumor tissue was isolated, processed, and stained with mass cytometry antibodies as previously described (15). Briefly, following manual dissociation, tumor samples were enzymatically digested with Liberase TL (Roche) and DNase I (Roche) for 30 min at 37 °C. Single-cell suspensions were generated by mashing digested tumors through 70-μm filters and purified using a Histopaque-1119 (Sigma–Aldrich) discontinuous gradient on RPMI cell culture media centrifuged at 2,000 rpm for 20 min. Up to 2.5 × 106 cells were washed with fluorescence-activated cell sorting (FACS) buffer and then incubated with 2% of each bovine, murine, rat, hamster, and rabbit serum in PBS with 25 μg/mL 2.4G2 antibody at 4 °C for 10 min prior to surface staining with an antibody mixture at 4 °C for 30 min in a 50-μL volume. A 50-μL volume of 2× solution of 194Pt monoisotopic cisplatin (Fluidigm) in PBS was added directly to the samples for a final concentration of 2.5 μM for 1 min. Samples were immediately washed twice with FACS buffer and labeled using palladium barcoding per the manufacturer’s protocol. Samples were prepared for intracellular staining using a FOXP3 fix and permeabilization kit per the manufacturer’s protocol (Thermo Fisher Scientific) and then stained for 30 min at room temperature. Following 2 washes with permeabilization buffer and FACS buffer, samples were fixed again using 1.6% paraformaldehyde in PBS supplemented with 100 nM iridium nucleic acid intercalator (Fluidigm). Following fixation, cells were washed twice with 0.5% bovine serum albumin (BSA) PBS and 0.1% BSA water, filtered, and resuspended in 0.1% BSA water. Barcoded samples were then analyzed using a Helios mass cytometer (Fluidigm).

For analysis of human peripheral blood, cryopreserved samples were thawed and immediately washed prior to additional processing. Samples were then purified on a Histopaque gradient prior to staining using a similar protocol as described above for murine tumor infiltrating lymphocytes (TILs). Surface and intracellular antibody staining panels are detailed in .

Mass Cytometry Data Analysis.

Mass cytometry flow cytometry standard (FCS) data files were concatenated and bead-normalized using Helios software (Fluidigm), debarcoded using Helios software or the cell debarcoder in MATLAB (53), and exported for downstream analysis. Total CD45+ or T cell populations were manually identified and exported using negative and positive gating strategies of lineage markers in FlowJo. T cell-gated data were not subsampled and were subjected to an arcsinh transformation with a cofactor of 4. Clustering and metaclustering were performed using the well-validated algorithm PhenoGraph (16). A metaclustering approach was used to identify MC38 tumor-infiltrating T cell populations using a Euclidean distance metric with k = 30 for initial cluster identification at the per mouse level and a cosine distance metric with k = 15 for metacluster assignment across cohorts. A similar metaclustering approach with these variable values was used for identification of T cell populations in publicly available human lung tumor mass cytometry data and human peripheral blood data. For all clustering approaches, samples with fewer than 1,000 events were excluded from the analysis.

The human peripheral blood mass cytometry data were acquired in 4 batches consisting of analytical samples and technical controls (repeated sampling of cryopreserved normal donors). Comparison of controls across cohorts (runs) revealed a significant batch effect (, Left). When we performed metaclustering agnostic of multiple parameters (including TIM3, HLA-DR, CD28, and CD8a) the observed batch effect was reduced (, Center and Right), suggesting that these specific channels contributed to most of the technical batch effect. Upon further inspection, these 4 channels did not appear to have any significant flaw other than a batch effect (an apparent shift in the background, and thus average signal). This initial approach to mitigate the batch effect present in the data ameliorated much of the batch effect but precluded the use of these parameters in the assignment of clusters across cohorts.

In order to more systematically correct the observed batch effect across cohorts, we chose to measure z-score protein expression within each cohort separately. More precisely, we subtracted the cohort mean expression of all proteins from each cell (within that cohort) and divided the result by the cohort level SD of respective protein expression. This procedure results in the protein expression distribution across cohorts having the same mean (0) and the same variance (1) but does not alter it within cohort biological variances; therefore, this minimizes any technical batch effect between cohorts. For this, we utilized the zscore function in MATLAB.

Subsequent analyses such as tSNE and PhenoGraph were performed using similar default parameters. In the case of PhenoGraph clustering of human samples, k = 30 was used to construct the graph. To determine whether this procedure minimized the technical batch effect between cohorts, replicate normal donor samples between cohorts were compared. tSNE overlays between these samples, indicating that this procedure removed the technical batch effects ().

Heat maps were generated using MATLAB. Individual parameters were normalized to the maximal mean expression of each parameter across clusters such that for each parameter, 1 cluster displayed maximal expression. To help in interpretation of the data, only clusters over 1% frequency were plotted.

Analysis of Publicly Available Data.

For the analysis of previously published mass cytometry data (20), FCS files were downloaded from the Flow Repository (https://flowrepository.org; repository ID: FR-FCM-ZY9M). Data were analyzed using PhenoGraph and FlowJo as described above.

Datasets.

Mass cytometry data have been deposited in the Flow Repository (murine TIL data; repository ID: FR-FCM-ZYQQ) (54). Human peripheral blood data from patients being treated with immune checkpoint blockade therapies were also deposited in the Flow Repository (repository ID: FCM-FR-ZYQR) (55).

Statistics.

Statistics were performed in Prism 7.0 (GraphPad). Tumor growth and cellular frequencies between 2 groups were compared using 2-tailed t tests with Welch’s correction or 1-way ANOVA with Sidak’s multiple testing correction. Cluster frequencies were compared using 2-way ANOVA with Tukey’s multiple testing correction. Correlations were displayed with linear regression lines with Spearman’s rank correlation.