Cytokine network analysis of immune responses before and after autologous dendritic cell and tumor cell vaccine immunotherapies in a randomized trial
Journal of Translational Medicine volume 18, Article number: 176 (2020)
In a randomized phase II trial conducted in patients with metastatic melanoma, patient-specific autologous dendritic cell vaccines (DCV) were associated with longer survival than autologous tumor cell vaccines (TCV). Both vaccines presented antigens from cell-renewing autologous tumor cells. The current analysis was performed to better understand the immune responses induced by these vaccines, and their association with survival.
110 proteomic markers were measured at a week-0 baseline, 1 week before the first of 3 weekly vaccine injections, and at week-4, 1 week after the third injection. Data was presented as a deviation from normal controls. A two-component principal component (PC) statistical analysis and discriminant analysis were performed on this data set for all patients and for each treatment cohort.
At baseline PC-1 contained 64.4% of the variance and included the majority of cytokines associated with Th1 and Th2 responses, which positively correlated with beta-2-microglobulin (B2M), programmed death protein-1 (PD-1) and transforming growth factor beta (TGFβ1). Results were similar at baseline for both treatment cohorts. After three injections, DCV-treated patients showed correlative grouping among Th1/Th17 cytokines on PC-1, with an inverse correlation with B2M, FAS, and IL-18, and correlations among immunoglobulins in PC-2. TCV-treated patients showed a positive correlation on PC-1 among most of the cytokines and tumor markers B2M and FAS receptor. There were also correlative changes of IL12p40 with both Th1 and Th2 cytokines and TGFβ1. Discriminant analysis provided additional evidence that DCV was associated with innate, Th1/Th17, and Th2 responses while TCV was only associated with innate and Th2 responses.
These analyses confirm that DCV induced a different immune response than that induced by TCV, and these immune responses were associated with improved survival.
Trial registration Clinical trials.gov NCT004936930 retrospectively registered 28 July 2009
The development of effective therapeutic cancer vaccines has been an elusive goal for several decades. The Nobel Prize winning research of Ralph Steinman [1, 2] has led to a resurgence of interest in therapeutic dendritic cell vaccines (DCV) [3, 4], especially in patients with metastatic melanoma . Clinical studies utilizing autologous dendritic cells loaded with antigens from autologous tumor cells have been especially promising [6,7,8,9]. A randomized phase II trial tested two vaccines featuring autologous tumor antigens (ATA): injections of autologous dendritic cells loaded ex vivo with antigens from autologous tumor cell lines (DCV), and tumor cell vaccines (TCV) consisting of irradiated autologous proliferating tumor cells [8, 9]. An early analysis showed that DCV was associated with better survival , and this was confirmed when 5-year follow up showed a more than doubling of median survival and 3-year survival rate, and a 70% reduction in the risk of death . This DCV approach is currently being tested in phase II trials in glioblastoma and ovarian cancer, and in a phase IB trial in combination with monoclonal antibodies to programmed death-1 protein (PD-1) in melanoma patients.
The human immune response to pathogens and cancer has many interacting components encompassed by concepts of innate and adaptive immunity [10, 11]. Inducing new immune responses, or enhancing existing weak immune responses to cancer associated antigens, is the goal of anti-cancer vaccines [12, 13]. Innate immunity features natural killer cells and associated cytokines [14,15,16,17], and macrophages which also produce cytokines and their function is affected by cytokines [18,19,20,21]. Adaptive immunity includes Th1 cell-mediate immunity and Th2 humoral immunity . Th1 encompasses cellular responses, notably cytotoxic T lymphocytes and associated cytokines, especially interferon-gamma [23,24,25]. Th2 encompasses humoral immunoglobulin immune responses with class switching of immunoglobulin isotypes orchestrated by B lymphocytes with the assistance of helper T lymphocytes, and these responses are associated with various cytokines [26,27,28,29,30]. There is also a Th17 response that features cells helper T cells that secrete IL17 [31,32,33]. These cells can augment immunosuppressive or immune enhancing effects depending on the microenvironment and local cytokines [34,35,36,37]. At the nexus of innate and adaptive immunity, are antigen-presenting cells, especially dendritic cells , that are characterized by secretion of IL12 when activated [39,40,41,42,43,44]. As part of maintaining homeostasis of the immune system, there are also a number of cytokines including such as PD-1 [45, 46], IL10 [47,48,49], transforming growth factor beta (TGFβ) [50,51,52], beta-2-microglobulin (B2M) [53, 54], and Fas (CD95) , that are associated with suppressed immune responses. In terms of cytokines, IL17 [34, 36, 37], and IL18 [56, 57], can be immune-augmenting or immunosuppressive depending on other signals in the microenvironment.
One of the objectives of vaccine clinical trials is to increase understanding of the immune responses that are induced or enhanced by such vaccines and correlating these responses with survival. In the randomized phase 2 analysis of 110 cytokines showed that variation in specific cytokine groupings were similar at baseline but differed markedly following three vaccinations . A combination of baseline soluble programmed death protein-1 (sPD-1) and changes in sPD-1 after three injections was strongly predictive of 3-year survival in DCV-treated patients, but not TCV-treated . The current analysis was conducted in an effort to better understand the nature of these immune responses in terms of classical concepts of innate and adaptive immune responses as reflected by changes in cytokines in response to vaccine therapy .
Blood samples from melanoma patients enrolled in a randomized phase 2 trial (clinicaltrials.gov NCT00436930) were obtained at week-0, 1 week before the first of 3 weekly injections of DCV or TCV vaccines, and week-4, 1 week after the third weekly injection [8, 9]. The trial was approved by the Western Institutional Review Board (Seattle, WA., WIRB® Protocol #20090753). Patients gave written informed consent for randomization to DCV or TCV, and blood collection and analysis. The protocol and manufacturing procedures were reviewed by the US Food and Drug Administration (BB-IND 5838 and BB-IND 8554). TCV consisted of irradiated autologous tumor cells from a short-term cell line; DCV consisted of autologous dendritic cells incubated with irradiated autologous tumor cells [8, 9]. Both DCV and TCV were admixed in granulocyte-macrophage colony stimulating factor (GM-CSF) just prior to subcutaneous injections scheduled for weeks 1, 2, 3, 8, 12, 16, 20, and 24.
Analysis of serum markers
Cryopreserved 200-microliter serum samples from week-0 and week-4 were analyzed for 110 cytokines, growth factors, proteases, soluble receptors and other proteins as shown in a supplementary table of a previous publication , using a quantitative, multiplex enzyme-linked immunosorbent assay (Quantibody® Cytokine Array, Raybiotech, Inc., Norcross, GA.). Values were expressed as absolute concentration (pg/mL) and as percentage differences above or below the mean value from three normal controls.
Principal component analysis
Principal component analysis (PCA) was performed using IBM SPSS Statistics V26. PCA transforms data into a coordinate system by creating new uncorrelated variables that successively maximize variance [59, 60]. The most important use of PCA is to represent a multivariate data table as smaller sets of variables (summary indices) in order to observe trends, jumps, clusters and outliers. PCA is valuable for exploratory analysis of extensive data and especially useful for integrating genomics and proteomic datasets in effort to understand biological processes [61,62,63].
In a two-component PCA model variables are distributed in a two-dimensional plane so that the largest amount of variance is grouped in principal component-1 (PC-1) along one axis and the next highest is grouped in principal component-2 (PC-2) along an orthogonal axis. The distribution is based on correlation coefficients that vary from 0 to 1 or 0 to − 1 from the origin, which quantitates the strength of positive or negative correlation. Positively correlated variables cluster in the same side of at least one of the components (influential), or in the same quadrant (correlated). Negatively correlated variables are positioned on opposite sides of the origin in diagonally opposed quadrants. After plotting, the component axes can be rotated for optimal delineation of the groups, following established procedures (e.g. Varimax, Quartimax, Promax, etc.) that give a best fit among variables. The rotation changes only the positioning vector of the components to each-other for easier interpretation.
The Kaiser–Meyer–Olkin Measure of Sampling Adequacy (KMO) was used to indicate the proportion of variance that might be attributed to specific variables. High values (close to 1.0) indicate that a factor analysis of the data will be useful as opposed to values less than 0.50. Bartlett’s Test of Sphericity (BTS) was used to test the hypothesis that the correlation matrix is an identity matrix, which would indicate that variables are unrelated and therefore unsuitable for structure detection. Lower significance levels (less than 0.05) indicate that a factor analysis of the data may be useful.
Discriminant analysis (DA) was used to interpret multiple discriminant functions arising from analyses involving more than two groups and more than one variable . DA is very useful for detecting variables that discriminate between different groups and can classify cases into different groups with a better than chance accuracy. It is especially useful for interpreting large data bases such as result from proteomics [65,66,67]. The objective is to develop discriminant functions that are the linear combination of independent variables (i.e. change of cytokines) that will discriminate between the categories of the dependent variable (i.e. survival groups). It enables one to examine whether significant differences exist among the groups based on the predictor variables and evaluates the accuracy of the classification. Several variables are included in order to see which one(s) contribute to discriminating between groups, and then a matrix of total variances and covariances is created as well as a matrix of pooled within-group variances and covariances. Matrices are compared via multivariate F tests to determine whether there are significant differences (with regard to all variables) between groups. This procedure is identical to multivariate analysis of variance, or MANOVA.
In stepwise discriminant function analysis, a model of discrimination is built step-by-step. At each step all variables are reviewed and evaluated to determine which one contributes most to discriminating between groups. That variable is included in the model, and the process starts again. The stepwise procedure is “guided” by F-to-enter and F-to-remove values for statistical-based discrimination between groups; i.e., it is a measure of the extent to which a variable makes a unique contribution to the prediction of group membership. When performing a multiple group discriminant analysis, an optimal combination of variables is automatically determined so that the first function provides the most discrimination between groups, the second provides second most, and so on. The larger the standardized coefficient, the greater the variable’s contribution to discrimination. Another way to determine which variables define a particular discriminant function is to examine the factor structure coefficients, which are the correlations between variables in the model and the discriminant functions. A Bonferroni correction of p values was used in analyses that included multiple comparisons.
Principal cytokine pathway analysis
The pathway analysis considered the findings from both PCA and DA. The methods to investigate the effect of these pathways on survival included linear regression analysis and tests of equality of group means (univariate ANOVA). In order to interpret multiple discriminant functions arising from analyses with more than two groups and more than one variable, we first tested the different functions for statistical significance, and only considered significant functions for further examination. Next, we looked at the standardized b coefficients for each variable for each significant function. The larger the standardized b coefficient, the larger is the respective variable’s unique contribution to the discrimination specified by the respective discriminant function. The classification matrix was used to determine how well the current classification functions predict group membership as defined by survival. The classification matrix shows the number of cases that were correctly classified (on the diagonal of the matrix) and those that were misclassified. To perform the discriminant analysis each treatment arm was sub-divided into three categories: survivors over 60 months, and the remaining patients split at the median survival for each treatment arm.
Paired week-0 and week-4 blood samples were available for 22 of 24 TCV-treated patients and 17 of 18 DCV-treated patients . The two missing TCV-treated patients did not have a week-4 sample because of rapidly progressing metastatic disease and both died within 2 months of enrollment. The missing DCV-treated patient was alive at 5 years but had rescinded permission to study his blood samples.
Principal component analysis at baseline
The PCA of cytokines and immunoglobulins (Ig) at baseline are shown in Fig. 1. Before vaccine treatment was initiated, the distribution of PCA values for TCV-treated and DCV-treated patients were similar. Deviation from normal values for all 39 patients defined two distinctive groups based on the loading plot. The correlative group of cytokines on component 1 (PC1) accounted for 64.4% of variance and a correlative group of immunoglobulins (PC2) accounted for 7.5% of variance (Additional file 1). As shown in Table 1, the majority of the cytokines with correlations greater than 0.5 are those classically associated with Th1, and Th2 immune responses.
This collection of cytokines positively correlated with B2M, PD1 and TGFβ1 suggesting coexistence of pre-existing tumor-associated inflammation and immunosuppression.
Principal component analysis of cytokine changes after 3 weekly injections
Figure 2 shows the PCA of cytokine changes after 3 weekly vaccine injections. The change from baseline for each cytokine was calculated by subtracting the baseline value from the week-4 value and dividing the difference by the smaller of the baseline or the week-4 value to avoid negative numbers when week-4 levels were lower than baseline values. Using this method, a positive or negative correlation vector was calculated for each variable.
Cytokine changes after 3 TCV injections
In TCV-treated patients, PCA showed the contribution of multiple factors with the combination of PC1 and PC2 responsible for only 54% of variance (Fig. 2a, Additional file 2). There was a positive correlation on PC1 for most of the cytokines suggesting the evolution of tumor-associated inflammation, and on PC2 for tumor markers (B2M, FAS receptor) along with Th2-associated cytokines (Fig. 2a, Table 2). Changes in IL12 correlated with both Th1- and Th2-cytokines suggesting an adaptive response mediated by antigen presenting cells (APC) that had been suppressed by TGFβ1. Thus, after three injections, the cytokine changes in the TCV group were most consistent with a Th2 adaptive response and an innate response.
PCA limited to factors with component values above 0.5 (Table 3), revealed a significant association of the variables KMO = 0.701, BTS < 0.001 (Additional file 3). The first two components contributed equally (40.1% and 39.5% of the variance respectively (Additional file 4), thus organizing the variables into two groups. One group contains variables associated with IL6, IL15, and IFNγ, but not associated with the APC-associated cytokine IL12p70, while the other group contains markers associated with APCs. This suggests an adaptive Th2 response mediated through immunoglobulins, and possibly a suppressed Th17/Th1 response based on the association with TGFβ1. Survival correlated negatively with the suppressed Th1/Innate response, but slightly positively with the increased immunoglobulins, suggesting dominance of a Th2 response in TCV-treated patients. TCV-treated patients exhibited changes in Th2 cytokines, suggesting that antigens on the injected irradiated tumor cells primarily elicited a Th2 response mediated by endogenous APCs.
Cytokine changes after 3 DCV injections
In DCV-treated patients there was a correlation between Th1/Th17 cytokines on PC1 and immunoglobulins on PC2 (Fig. 2b) with PC1 and PC2 accounting for 39.1% and 12.6% of variance respectively (Additional file 5). The first two components were used for further analysis and component plotting (Table 4).
Inflammatory cytokines correlated strongly with an APC-driven effect (IL12p70) through Th1 (IFNγ, TNFα) combined with a Th17 (IL17) response, as well as through Th2 driven by IL12p70 in association with IL5 (hypersensitivity) and immunoglobulins on component 2. Survival correlated with the Th1/Th17 factors on PC1. On both PC1 and PC2, survival correlated positively with IgG1 and IgG3 but negatively with IGM, IgG2 and IgG4. The inverse correlation of IgG4/IgG2 with IgG1/IgG3 is consistent with immunoglobulin class-switching suggesting that antigen-specific immunoglobulin responses may have contributed to survival. It is also noteworthy that B2M and FAS receptor are negatively correlated with the immune response and survival. Collectively the data suggests an adaptive response in the DCV-treated group that included immunoglobulin class-switching and a Th17/Th1 response. The cytotoxic Th17 response is not correlated with IL23 [68, 69], a cytokine that is expected to drive naïve CD4+ cells towards Th17 lineages associated with immunosuppression. Instead, IL17 is strongly correlated with TNFα and IFNγ, consistent with conversion of an existing cognate population of Th17 cells from a tolerizing to a cytotoxicity-facilitating phenotype [34, 36, 37].
Closer examination of factors positively correlated on PC1 confirms a significant association of this group, KMO = 0.718 and BTS < 0.001 (Additional file 6). The analysis of the reduced set of factors identifies two components that account for 49.0% and 33.3% of the variance. Within component 1 IL17 has the highest coefficient, associated with a Th1 response driven by APC (IL12p70) (Table 5). The second component contains the innate response associated with Th2 (IL4, TGFβ1) and other pleiotropic factors (IL7, IL8) (Table 5). Both components were associated with IFNγ and TNFα. Thus, after three DCV-injections, the analysis suggests a multifaceted response driven by the Th1/Th17 cytotoxic (Th1-like) pathway and a Th2 immunoglobulin response.
Discriminant analysis by treatment arm and survival
Each treatment arm was sub-divided into three categories based on overall survival (Table 6).
DA was then applied to cytokine data for each patient in each treatment arm with results as displayed in Fig. 3. Regardless of treatment DA correctly classified patients into their appropriate survival subgroup (Additional file 7).
Discriminant analysis of TCV-treated patients
DA based on changes of inflammatory markers among TCV-treated patients yielded a discriminant function for survival (p = 0.015, Additional file 8 that accounted for 100% of variance with 90.2% of variance in the first discriminant function (Additional file 9). DA also accurately classified the three survivor subgroups (Fig. 3a). However, none of the variables reached statistical significance (Additional file 10). Stepwise DA failed to identify any significant variables qualified for further analysis based on F function entry criteria (p < 0.05). Thus, DA of TCV patients was unable to identify a cytokine pattern to explain the survival distribution.
Discriminant analysis of DCV-treated patients
DA based on changes of inflammatory markers among DCV-treated patients yielded a discriminant function for survival (p = 0.020, Additional file 11) that accounted for 100% of variance with 92.9% of variance in the first discriminant function (Additional file 12). DA accurately classified the three survivor subgroups (Fig. 3b). Stepwise DA identified IL17 and TGFβ1 as the most important variables that discriminated among survival. These two variables identified two significant discriminant functions (p = 0.003 and p = 0.023, Additional file 13) that explained 100% of variance with 73.1% in the first function (Additional file 14). The discriminant function of TGFβ1 and IL17 correctly classified 70.6% of the DCV-treated patients (Fig. 3c), All coefficients were less than 0.05 for each variable and each survival group (Table 7).
Thus, persistence of TGFβ1 and a decrease in IL17 was associated with relatively short survival. Persistence or increase of IL17 and an increase of TGFβ1 (consistent with an ongoing response that could be either anti-tumor or tolerizing depending on the TGFβ1 dominance) was associated with intermediate survival. A decrease in TGFβ1 combined with an increase of IL17 was associated with long-term survival.
Immunoglobulins (Th2 pathway)
PCA suggested there was an existing immunoglobulin response at baseline, but also a new immunoglobulin response in each treatment arm after three injections; therefore, this pathway was examined relative to survival. Baseline immunoglobulin values other than IgM were moderately elevated compared to healthy controls. We assumed that the vaccines were inducing a new response; so, we analyzed changes in IgM that might be expected after exposure to new antigens. The correlation between baseline IgM levels and survival is shown for TCV-treated patients in Fig. 4a, and for DCV-treated patients in Fig. 4b. The average IgM of healthy controls was 120.2 mg/dL, which is in the 40–230 mg/dL standard range of normal IgM for ages 45 years and older . The changes in mean and median Ig levels after three injections of either vaccine are also shown in Table 8.
The slightly increased baseline IgM in TCV-treated patients correlated with longer survival (Fig. 4a). The baseline IgM level was normal in DCV-treated patients and not correlated with survival (Fig. 4b). Linear regression analysis was confirmed by Cox’s proportional hazards model for survival using baseline IgM as a cofactor in TCV-treated patients (p = 0.014, Additional file 14). After three injections average IgM increased by 31.5% in the DCV-treated patients (Table 8) with 11 of 17 patients having above-normal values. This suggests that a new IgM response was induced in DCV-treated patients (p = 0.102, Wilcoxon signed rank test, p = 0.090 paired t-test).
The DCV regression graph shows two distinct distributions on either side of 36 months (Fig. 4b). This suggests that long-term survivors may have different mechanisms of response; therefore, we analyzed patients who survived less than or more than 60 months separately in both treatment arms (Table 9).
TCV-treated patients who survived more than 60 months had elevated baseline IgM levels (p = 0.013, ANOVA) that were higher than for TCV-treated patients who survived less than 60 months (p = 0.033 independent sample t-tests, 0.034 Mann–Whitney), higher than for DCV-treated patients who survived less than 60 months (p = 0.008 independent sample t-tests, 0.027 Mann–Whitney), but not higher than in DCV-treated patients who survived 60+ months (p = 0.183 independent sample t-tests, p = 0.175 Mann–Whitney). Perhaps patients who survived 60+ months had an existing anti-tumor Th2 response that contributed to better survival even in the absence of vaccine treatment. After three injections, only DCV-treated patients who survived less than 60 months had an increase in IgM (p = 0.045 paired t-test, p = 0.06 Wilcoxon signed rank), but this did not correlate with survival (Additional file 15).
T-helper cytotoxic pathways (Th1, Th17 pathways)
The Th17 phenotype is induced by a combination of IL6, IL23 and TGFβ1 [31,32,33], and functionally by the amplitude of TGFβ1 . Th17 cells have a dual-state plasticity; a high TGFβ1 induces regulatory or immune-tolerizing effects while low TGFβ1 allows a Th1-like cytotoxic response that includes IFNγ and TNFα [34, 36, 37]. A strong IL12/IL23 response leading to an IL4-mediated response suggests a new T-helper (CD4) response with a new Th17 component . The lack of association with IL23 suggests a helper response that does not require Th17 . Positive correlation with IFNγ suggests a Th1 response; positive correlation with IL10 suggests a Th2 response; and positive correlation with TGFβ1 suggests a regulatory T lymphocyte (Treg) response. In the tumor microenvironment Th17 cells are likely antigen-cognate and can change from a suppressive to cytotoxic state in response to local factors [36, 37]. The sudden increase of IL17 in association with TNFα suggests a pro-inflammatory switch of Th17 cells to a cytotoxic helper function rather than suppression. An IL23 increase could reflect lineage stimulation from naïve CD4+ cells that can be directed to both cytotoxic and regulatory pathways, depending on the local environment. Downstream effectors of Th17 can result in cytotoxic lymphocytes (Th1-like pathway) [36, 37], a Th2 pathway through IgM induction , or suppressive (as in a Treg pathway) .
Because Th17 is a versatile cell type, estimations of Th17 function should include additional factors such as IFNγ, TNF, TGFβ1, and Th17-associated cytokines (IL17, IL21, IL22, IL27, IL31, IL33). All patients were included in linear regression analysis between cytokines and survival, but long-term follow up ended at 60 months; therefore, patients surviving longer were excluded from linear regression analysis of the cytokines selected as impacting survival by PCA and DA. Cytokines identified by DA accurately classified patients into survival groups. The associations between survival, IL17, IL12p70, IFNγ, and TNFα are shown in Fig. 5. The change in IL17 correlated with survival for DCV-treated patients who survived less than 60 months (Fig. 5a). Increases in IL17 were correlated with increased IL12p70, a cytokine produced by DCs (Fig. 5b). IL17 also correlated with IFNγ (Fig. 5c) and TNFα (Fig. 5d), both of which are components of Th1 responses. This suggests that the increased IL17 is associated with a Th1-like response that was triggered by antigen-loaded DC that secreted IL12. To explore the source of IL17, we investigated its association with IL23, another cytokine that can generate IL17-secreting Th17 lymphocytes [73, 74]. IL17 did not correlate with IL23 (p = 0.248, ANOVA); so, the source of IL17 was not a new population of Th17 cells, but possibly resulted from conversion of an existing antigen-cognate Th17 population from a tolerizing-state to a cytotoxicity-inducing state [34, 36, 37].
In the TCV arm there were no post-treatment changes in cytokines that correlated with survival except for IL17, which correlated negatively (Fig. 6a). In TCV-treated patients IL17 positively correlated with TGFβ1 (Fig. 6b), and IFNγ (p = 0.038) and cytokines associated with innate immune responses: IL15 (p = 0.003), IL8 (p = 0.002) or Th2-associated cytokines: IL13 (p = 0.038), IL10 (p = 0.024), IL7 (p = 0.038), IL6 (p = 0.01) IL5 (p = 0.07) and IL2 (p = 0.018). The cytokine changes were not consistent with a cytotoxic Th1 response, but rather with conversion of an existing antigen-cognate Th17 population from a tolerizing-state to a cytotoxicity-inducing state, as would be expected after inoculation with an antigen that is recognized and presented by endogenous APCs in vivo. The association with TGFβ1 suggests that during antigen processing and presentation in vivo, TGFβ1 affected the immune response by promoting Th2 subpopulations and perhaps immunosuppressive Tregs.
In DCV-treated patients DCs were antigen-loaded ex vivo in the absence of suppressive cytokines (such as TGFβ1). After subcutaneous injection, it appears that the large dose of activated DCs produced sufficient amounts of IL12 to trigger a Th1 response. The first response likely came from cross-presentation of antigens directly to cytotoxic T lymphocytes (CTL) and/or from conversion of inherently plastic antigen-cognate Th17 cells into pro-cytotoxic states. It appears this triggering signal was only provided by ex vivo loaded DC leading to increases in IL17, IFNγ, and TNFα and was not observed in TCV-treated patients.
After this first wave of Th17 cytotoxic conversion, a second wave of helper population likely developed based on the positive association between IL12p70 and IL4 in both TCV-treated (ANOVA, p = 0.001) and DCV-treated patients (ANOVA, p = 0.005). This is consistent with de-novo antigen induction of a Th2 response, but in DCV-treated patients this was accompanied by increases in IFNγ (p = 0.001) and TNFα (p < 0.0001), which were not seen following TCV. The association of cytotoxic cytokines in DCV-treated patients suggests new helper cells sustained a persistent Th1-response, while in TCV-treated patients only a Th2-response persisted. Furthermore, in the DCV group there was a new Th2 response evidenced by PCA and DA of the immunoglobulins that showed a correlative increase and a class switching of IgG1 and IgG3 as opposed to IgG2, IgG4, and IgM. This was only observed in the DCV group, suggesting that while both treatment arms respond with a Th2 mechanism, the ex vivo-generated DCs may present antigens at better signal to noise ratio and thereby enhance or induce de-novo immunoglobulin responses.
IL12 is a cytokine produced by APC- when they are activated by exposure to antigens [40, 44]. In terms of antigen processing, the differences between the treatment arms is the source of the APCs: endogenous in situ APC in the case of TCV, and ex vivo antigen loaded DC in the case of DCV. In DCV the antigens are processed by DCs derived ex vivo from peripheral blood monocyte . DCs mature and migrate to lymph nodes after exposure to antigen, where they contact effector T-cells, B-cells and natural killer cells (NKs) [40, 68]. Efficient presentation of the antigen is regulated by the interaction of the MHCs with T-cell receptors (TCR), and by regulatory and costimulatory connectors as well as by cytokine and chemokine signals [69, 75]. In addition to T-cell interaction, the antigen processing and presentation is modified by the state of the APCs. First there is a difference in the maturation stage of the APCs. Studies have shown that DC maturation is accompanied by a marked reorganization of endocytic compartments [76,77,78], and a concomitant inhibition of antigen uptake [79,80,83]. Antigen uptake and processing is limited to immature DCs, which contribute to the functional distinction between immature and mature DCs . Endogenous APCs represent a minor population of cells at the site of TCV injections. In contrast the ex vivo loaded DCV injection contains a massive dose of 1 to 30 million immature DCs that originated from peripheral monocytes, that have not been exposed to additional maturation factors such as lipopolysaccharides. The phenotype and functionality of the ex vivo derived DCs is expected to be substantially different from the in situ APCs. As previously described, the immature DCs generated in vitro are more efficient for cross-presentation initially . In addition to the massive number of simultaneously activated APCs (DCs) in DCV, the maturation differences could explain the initial Th1 cytotoxic response, cytotoxic conversion of Th17, followed by the Th2 response which could be mediated by endogenous APCs responding to apoptotic antigen-loaded DC, as opposed to the predominate Th2 response associated with the smaller number of more mature and less numerous APCs induced by TCV.
DCs are at the nexus of innate and adaptive immunity and have evolved to orchestrate a multi-pronged immune response . Regardless of antigen source, DC are able to present antigen by both MHC I and MHC II pathways to induce both Th1 and Th2 immune responses [83,84,85,86,89], both of which are necessary for an optimal immune response to tumor antigen . DC also can induce Th17 responses . The role of Th17 cells in cancer is of increasing interest, especially because Th17 cells exhibit a plasticity that can result in their differentiation into Treg or Th17/Th1 cells [37, 92, 93]. The latter increasingly lose the ability to secrete IL17, but are able to secrete larger quantities of TNF, IL2, GM-CSF, and IFNγ than classical Th1 cells. It appears that TGFβ drives Th17 cells to become Tregs, and the absence of TGFβ and the presence of cytokines such as IL12 and IL23 is required for conversion to the Th17/Th1 phenotype. Although Th17 cells and Th17/Th1 cells are not believed to be cytotoxic, this Th1 helper phenotype is associated with increased cytotoxic T lymphocytes (CTL) in the tumor microenvironment. The ability of antigen-loaded DC to induce differentiation of Th17 cells into Th1 helper cells associated with anti-tumor effects was demonstrated in animal models . Our analysis suggests that the patient-specific DCVs induced similar immune responses.
Pattern-recognition receptors (PRR) are proteins expressed by immune cells that recognize pathogen-associated molecular patterns as danger signals [95, 96]. IL6, IL8, and TNFα are produced in response to PRR that trigger a response from NK cells that is associated with production of IL15, IL18, and IFNγ [14,15,16,17]. Cytokines that regulate innate immunity are produced primarily by mononuclear phagocytes such as macrophages and DCs, although they can also be produced by T-lymphocytes, natural killer (NK) cells, endothelial cells, and mucosal epithelial cells. PRRs In injured tissue the innate immune system down-regulates effector mechanisms and restores homoeostasis via cytokines such as IL10 and TGFβ1 that are released by macrophages, preferentially the M2 subset, which can induce Tregs, inhibit pro-inflammatory cytokine production, and induce tissue healing by regulating extracellular matrix protein deposition and angiogenesis.
After treatment In the TCV arm, none of the changes in innate cytokines (pro-inflammatory or regulatory) correlated with overall survival or any of the survival subgroups. Th1/Th17 activity may be balancing a regulatory response through TGFβ1 and IL10. In the TCV arm, the lack of significant cytotoxic activity may tip the balance in favor of a regulatory mechanism. Linear regression of post-treatment changes in the NK-secreted cytokine IL15 revealed correlations with IFNγ (p < 0.001), IL2 (p < 0.001), IL5 (p < 0.001), IL6 (p < 0.001), IL7 (p < 0.001), IL8 (p < 0.001), IL10 (p < 0.001), IL13 (p < 0.001), IL17 (p < 0.001), IL21 (p = 0.009), IL23 (p = 0.002), IL12p70 (p < 0.001), TGFβ1 (p < 0.001), consistent with activation of the innate system.
Analysis of IL15 changes in the DCV arm showed a similar activation of the innate system including correlations with IFNγ (p = 0.002), IL2 (p = 0.005), IL5 (p = 0.001), IL6 (p < 0.001), IL7 (p < 0.001), IL8 (p < 0.001), IL10 (p = 0.01), IL13 (p < 0.022), IL17 (p < 0.001), IL12p70 (p < 0.002), TGFβ1 (p < 0.001), and TNFα (p = 0.016). The only difference between TCV and DCV arms is the correlation between changes in levels of IL15 and TNFα. In the DCV arm TNFα (p = 0.014), TGFβ1 (p = 0.027), and IL10 (p = 0.025) correlated with survival in patients who survived less than 60 months. Changes in TNFα, which is thought to be secreted by activated macrophages and T cells, correlated with increased IL12 in DCV-treated patients (Fig. 6c) but did not correlate with IL12 changes in TCV-treated patients (Fig. 6d). This also supports the hypothesis that conversion of the local inflammatory process from regulatory to cytotoxic occurred only in DCV-treated patients in a process that included Th17 and the innate immune system.
This study provides additional insight into differences in immune responses elicited by DCV and TCV. Both vaccines presented autologous tumor antigens but were associated with different immune responses and different survival benefit. While the results in the study are purely correlative, they are suggestive of underlying immunologic mechanisms of action. The major finding of this analysis is that DCV was associated with a multipronged immune response that included innate, Th2, and Th1/Th17 responses while the TCV immune responses were limited to innate and Th2. Direct correlation between Th1/Th17 changes and survival was also demonstrated. The results provide additional evidence that for therapeutic cancer vaccines it may be advantageous to present antigens by DC that were loaded with antigen ex vivo.
The analysis presented herein is complex and involved correlations of multiple variables derived from relatively small sample sizes. In addition to the analysis described in the manuscript, IBM SPSS (build 18.104.22.1688) was used to create a standard model using Automatic Linear Modeling function, targeting survival as an independent variable. A forward stepwise model was selected and AICC criteria, or F statistics (include effects with p < 0.05, remove effects at p > 0.1) was used for entry/removal of variables. Although the results were significant for many predictors, we felt that the reliability of the test might be criticized because of (small sample and large variability. Therefore, we performed a variable reduction with PCA and DA (variables normality verified) based on Bayesian estimation. Later, the synthetical variables obtained by PCA reduction of various groups of variables was examined in a regression model, but that approach did not provide meaningful results. The PCA graphics of the main components were presented herein without a statistical conclusion because representation of the differences between groups was provided with more confidence in the DA results. Nonetheless, these observations are hypothesis generating, or suggestive, clearly more cases would be needed for a robust statistical analysis.
With regard to immunoglobulins, there was evidence of a Th2 response that included increases in IgM in both treatment arms, especially in the DCV arm. The Th2 response was presumably mediated by endogenous, in situ DCs in the TCV-arm, while the ex vivo antigen-loaded DCV likely caused an early cytokine-mediated or cross-presentation response, followed by a new presentation of antigens through a typical helper pathway. In the TCV arm the changes are presumably in response to the additional antigenic stimulation provided by the irradiated tumors cells while in the DCV arm we believe that the initial cross presentation by DCV induced a new Th2 response as evidenced by immunoglobulin class switching. An association with a Th1 T-cell response is evidenced by the increase in IFNγ and TNFα with IL-12 only in the DCV arm, but this was not evident in the TCV arm.
Just prior to each subcutaneous injection, GM-CSF was admixed with tumor cells for TCV and dendritic cells for DCV for its adjuvant effects [97, 98], and specific effects on dendritic cells . GM-CSF was given in the same dose and schedule in both arms; therefore, in the absence of a GM-CSF alone control arm, we cannot identify the specific effects that GM-CSF induced in each arm. GM-CSF was not one of the cytokines measured in the analysis; so, we have not data re its association with other cytokines in the principal component analyses. We know that levels of granulocyte colony stimulating factor (G-CSF) and macrophage colony stimulating factor (M-CSF) did not change after three injections of either TCV or DCV. Thymus and activation regulated chemokine CCL17 (TARC), which is induced by GM-CSF , was elevated significantly and similarly in both arms, and therefore excluded from the analysis.
The clinical trial from which this data was derived was the first randomized study testing therapeutic cancer vaccines in which there was a difference in survival in the treatment arms and for which associated proteomic data has been analyzed extensively for changes in, and correlations, with circulating markers . Trying to decipher immune responses and their relation to clinical outcome in vaccine clinical trials is challenging. In one study in which colorectal cancer patients were treated with autologous dendritic cells loaded with allogeneic tumor cell lysate, plasma and serum samples were collected prior to vaccination and continuously during treatment . Patients classified as having stable disease had increasing levels of IL2, IL5, TNFα, IFNγ, and GM-CSF while increases in carcinoembryonic antigen (CEA) and (TIMP-1) levels were associated with progressive disease. No correlative changes were noted for IL1b, IL4, IL6, IL8, IL10, IL12, macrophage inflammatory protein 1beta (MIP-1β), Interferon-inducible protein 10 (IP-10), or Eotaxin. That study was limited by the lack of a control arm and the limited number of cytokines examined. In another study, immune monitoring was conducted in association with an 815-patient six-arm trial that randomized patients with surgically resected stage 3 and 4 melanoma to peptide vaccines or placebo with GM-CSF or placebo in patients of appropriate HLA-type, and GM-CSF or placebo in patients who were HLA-A2 negative [102, 103]. One challenge for correlative analyses was that none of the treatment variables impacted overall survival compared to placebo . The focus of the immune analysis was primarily on changes in immune cell phenotypes and their recognition of injected antigens rather than on changes in cytokines. There were no vaccine-specific correlations identified, and the cellular and humoral responses did not correlate with survival in the manner predicted .
The major strengths of our study include: (1) the use of data and samples from a randomized clinical trial that tested ATA presentation by two different cell sources (dendritic cells and cancer cells), (2) the availability of paired blood samples obtained at baseline and after 3 weekly injections, (3) the treatments tested were associated with different survivals, enabling a direct correlation between treatment, immune response, and survival, (4) the large number of immune markers tested, (5) the long-term follow-up for correlations with survival, and (6) the power of the statistical tools used to group positively and negatively correlated variables. The limitations of the study include: (1) the relatively small sample size, (2) paired samples were not available for 7.0% of the patients, (3) samples were not available for testing at earlier and later time points other than week-0 and week-4.
DCV induced a more effective immune response than that induced by TCV, and these immune responses were associated with improved survival. DCV was associated with innate, Th1/Th17, and Th2 responses while TCV was only associated with innate and Th2 responses. Figure 7 is an infographic summary of these changes.
Availability of data and materials
The clinical and proteomic datasets analyzed during the current study are available from the corresponding author on reasonable request.
Antigen presenting cell
Autologous tumor antigens
Bartlett’s Test of Sphericity
Cytotoxic T lymphocytes
Dendritic cell vaccine
FS-7-associated surface antigen (CD95)
Granulocyte-macrophage colony stimulating factor
Interferon-inducible protein 10 (CXCL10)
Kaiser–Meyer–Olkin Measure of Sampling Adequacy
Major histocompatibility complex
- MIP-1 β:
macrophage inflammatory protein 1beta (CCL4)
Principal component analysis
Programmed death protein-1
Soluble programmed death protein-1
Tissue metalloproteinase inhibitor 1
Tumor necrosis factor alpha
Transforming growth factor beta
T lymphocyte helper
Regulatory T cell
Tumor cell vaccine
Steinman RM, Cohn ZA. Identification of a novel cell type in peripheral lymphoid organs of mice. J Exp Med. 1973;137:1142–62. https://doi.org/10.1084/jem.137.5.1142.
Steinman RM. Decisions about dendritic cells: past, present, and future. Ann Rev Immunol. 2012;30:1–22. https://doi.org/10.1146/annurev-immunol-100311-102839.
Palucka K, Banchereau J. Cancer immunotherapy via dendritic cells. Nat Rev Cancer. 2012;12:265–77. https://doi.org/10.1038/nrc3258.
Sabado RL, Balan S, Bhardwaj N. Dendritic cell-based immunotherapy. Cell Res. 2017;27:74–95. https://doi.org/10.1038/cr.2016.157.
Dillman RO, Cornforth AN, Nistor GI. Dendritic cell vaccines for melanoma: past, present, and future. Melanoma Manage. 2016;3:273–89. https://doi.org/10.2217/mmt-2016-0014.
Nestle FO, Alijagic S, Gilliet M, Sun Y, Grabbe S, Dummer R, et al. Vaccination of melanoma patients with peptide- or tumor lysate-pulsed dendritic cells. Nat Med. 1998;4:328–32. https://doi.org/10.1038/nm0398-328.
de Rosa F, Ridolfi L, Fiammenghi L, Petrini M, Granato AM, Ancarani V, et al. Dendritic cell vaccination for metastatic melanoma: a 14-year monoinstitutional experience. Melanoma Res. 2017;27:351–7. https://doi.org/10.1097/CMR.0000000000000356.
Dillman RO, Cornforth AN, Depriest C, McClay EF, Amatruda TT, de Leon C, et al. Tumor stem cell antigens as consolidative active specific immunotherapy: a randomized phase II trial of dendritic cells versus tumor cells in patients with metastatic melanoma. J Immunother. 2012;35:641–9. https://doi.org/10.1097/CJI.0b013e31826f79c8.
Dillman RO, Cornforth AN, Nistor GI, McClay EF, Amatruda TT, Depriest C. Randomized phase II trial of autologous dendritic cell vaccines versus autologous tumor cell vaccines in metastatic melanoma: 5-year follow up and additional analyses. J Immunother Cancer. 2018;6:19. https://doi.org/10.1186/s40425-018-0330-1.
Chaplin DD. Overview of the immune response. J Allergy Clin Immunol. 2010;125(2 Suppl 2):S3–23. https://doi.org/10.1016/j.jaci.2009.12.980.
Vesely MD, Kershaw MH, Schreiber RD, Smyth MJ. Natural innate and adaptive immunity to cancer. Annu Rev Immunol. 2011;29:235–71. https://doi.org/10.1146/annurev-immunol-031210-101324.
Dillman RO. Is there a role for therapeutic cancer vaccines in the age of checkpoint inhibitors? Hum Vaccine Immunother. 2017;13:528–32. https://doi.org/10.1080/21645515.2016.1244149.
Sayour EJ, Mitchell DA. Manipulation of innate and adaptive immunity through cancer vaccines. J Immunol Res. 2017;2017:3145742. https://doi.org/10.1155/2017/3145742.
Lin WW, Karin M. A cytokine-mediated link between innate immunity, inflammation, and cancer. J Clin Invest. 2007;117:1175–83. https://doi.org/10.1172/JCI31537.
Marcus A, Gowen BG, Thompson TW, Iannello A, Ardolino M, Deng W, et al. Recognition of tumors by the innate immune system and natural killer cells. Adv Immunol. 2014;122:91–128. https://doi.org/10.1016/B978-0-12-800267-4.00003-1.
Geiger TL, Sun JC. Development and maturation of natural killer cells. Curr Opin Immunol. 2016;39:82–9. https://doi.org/10.1016/j.coi.2016.01.007.
Santana Carrero RM, Beceren-Braun F, Rivas SC, Hegde SM, Gangadharan A, Plote D, et al. IL-15 is a component of the inflammatory milieu in the tumor microenvironment promoting antitumor responses. Proc Natl Acad Sci USA. 2019;116(2):599–608. https://doi.org/10.1073/pnas.1814642116.
Mosser DM, Edwards JP. Exploring the full spectrum of macrophage activation. Nat Rev Immunol. 2008. https://doi.org/10.1038/nri2448.
Arango Duque GA, Descoteaux A. Macrophage cytokines: involvement in immunity and infectious diseases. Front Immunol. 2014;5:491. https://doi.org/10.3389/fimmu.2014.00491.
Gosselin D, Link VM, Romanoski CE, Fonseca GJ, Eichenfield DZ, Spann NJ, et al. Environment drives selection and function of enhancers controlling tissue-specific macrophage identities. Cell. 2014;159:1327–40. https://doi.org/10.1016/j.cell.2014.11.023.
DeNardo DG, Ruffell B. Macrophages as regulators of tumour immunity and immunotherapy. Nat Rev Immunol. 2019;19:369–82. https://doi.org/10.1038/s41577-019-0127-6.
Kaiko GE, Horvat JC, Beagley KW, Hansbro PM. Immunological decision-making: how does the immune system decide to mount a helper T-cell response? Immunology. 2008;123:326–38. https://doi.org/10.1111/j.1365-2567.2007.02719.x.
Ikeda H, Old LJ, Schreiber RD. The roles of IFN gamma in protection against tumor development and cancer immunoediting. Cytokine Growth Factor Rev. 2002;13:95–109. https://doi.org/10.1016/s1359-6101(01)00038-7.
Schroder K, Hertzog PJ, Ravasi T, Hume DA. Interferon-gamma: an overview of signals, mechanisms and functions. J Leukoc Biol. 2004;75:163–89. https://doi.org/10.1189/jlb.0603252.
Alspach E, Lussier DM, Schreiber RD. Interferon γ and its important roles in promoting and inhibiting spontaneous and therapeutic cancer immunity. Cold Spring Harb Perspect Biol. 2019. https://doi.org/10.1101/cshperspect.a028480.
DeKruyff RH, Rizzo LV, Umetsu DT. Induction of immunoglobulin synthesis by CD4+ T cell clones. Semin Immunol. 1993;5:421–30. https://doi.org/10.1006/smim.1993.1048.
Market E, Papavasiliou FN. V(D)J recombination and the evolution of the adaptive immune system. PLoS Biol. 2003;1:E16. https://doi.org/10.1371/journal.pbio.0000016.
Zhu J. T helper 2 (Th2) cell differentiation, type 2 innate lymphoid cell (ILC2) development and regulation of interleukin-4 (IL-4) and IL-13 production. Cytokine. 2015;75:14–24. https://doi.org/10.1016/j.cyto.2015.05.010.
Lloyd CM, Snelgrove RJ. Type 2 immunity: expanding our view. Sci Immunol. 2018. https://doi.org/10.1126/sciimmunol.aat1604.
Gieseck RL 3rd, Wilson MS, Wynn TA. Type 2 immunity in tissue repair and fibrosis. Nat Rev Immunol. 2018;18:62–76. https://doi.org/10.1038/nri.2017.90.
Gaffen SL. An overview of IL-17 function and signaling. Cytokine. 2008;43:402–7. https://doi.org/10.1016/j.cyto.2008.07.017.
Korn T, Bettelli E, Oukka M, Kuchroo VK. IL-17 and Th17 cells. Annu Rev Immunol. 2009;27:485–517. https://doi.org/10.1146/annurev.immunol.021908.132710.
Zuniga LA, Jain R, Haines C, Cua DJ. Th17 cell development from cradle to grave. Immunol Rev. 2013;252:78–88. https://doi.org/10.1111/imr.12036.
Muranski P, Restifo NP. Essentials of Th17 cell commitment and plasticity. Blood. 2013;121:2402–14. https://doi.org/10.1182/blood-2012-09-378653.
Llosa NJ, Geis AL, Thiele Orberg E, Housseau F. Interleukin-17 and type 17 helper T cells in cancer management and research. Immunotarg Ther. 2014;10(3):39–54. https://doi.org/10.2147/ITT.S56529(eCollection 2014).
Asadzadeh Z, Mohammadi H, Safarzadeh E, Hemmatzadeh M, Mahdian-Shakib A, Jadidi-Niaragh F, et al. The paradox of Th17 functions in tumor immunity. Cell Immunol. 2017;322:15–25. https://doi.org/10.1016/j.cellimm.2017.10.015.
Agalioti T, Villablanca EJ, Huber S, Gagliani N. TH17 cell plasticity: the role of dendritic cells and molecular mechanisms. J Autoimmun. 2018;87:50–60. https://doi.org/10.1016/j.jaut.2017.12.003.
Banchereau J, Steinman RM. Dendritic cells and the control of immunity. Nature. 1998;392:245–52. https://doi.org/10.1038/32588.
Macatonia SE, Hosken NA, Litton M, Vieira P, Hsieh CS, Culpepper JA, et al. Dendritic cells produce IL-12 and direct the development of Th1 cells from naive CD4+ T cells. J Immunol. 1995;154:5071–9.
Trinchieri G. Interleukin-12 and the regulation of innate resistance and adaptive immunity. Nat Rev Immunol. 2003;3:1330146. https://doi.org/10.1038/nri1001.
Cooper AM, Khader SA. Il-12p40: an inherently agonistic cytokine. Trends Immunol. 2007;28:33–8. https://doi.org/10.1016/j.it.2006.11.002.
Xu M, Mizoguchi I, Morishima N, Chiba Y, Mizuguchi J, Yoshimoto T. Regulation of antitumor immune responses by the IL-12 family cytokines, IL-12, Il-23, and IL-27. Clin Dev Immunol. 2010. https://doi.org/10.1155/2010/832454.
Johansson U, Walther-Jallow L, Hofmann A, Spetz AL. Dendritic cells are able to produce IL-12p70 after uptake of apoptotic cells. Immunobiology. 2011;216:251–5. https://doi.org/10.1016/j.imbio.2010.04.003.
Zheng H, Ban Y, Wei F, Ma X. Regulation of interleukin-12 production in antigen-presenting cells. Adv Exp Med Biol. 2016;941:117–38. https://doi.org/10.1007/978-94024-0921-5_6.
Okazaki T, Honjo T. PD-1 and PD-1 ligands: from discovery to clinical application. Int Immunol. 2007;19:813–24. https://doi.org/10.1093/intimm/dxm057.
Chamoto K, Al-Habsi M, Honjo T. Role of PD-1 in immunity and diseases. Curr Top Microbiol Immunol. 2017;410:75–97. https://doi.org/10.1007/82_2017_67.
Fiorentino DF, Zlotnik A, Vieira P, Mosmann TR, Howard M, Moore KW, O’Garra A. IL-10 acts on the antigen-presenting cell to inhibit cytokine production by Th1 cells. J Immunol. 1991;146:3444–51.
Moore KW, de Waal Malefyt R, Coffman RL, O’Garra A. Interleukin-10 and the interleukin-10 receptor. Annu Rev Immunol. 2001;19:683–765. https://doi.org/10.1146/annurev.immunol.19.1.683.
Oft M. IL-10: master switch from tumor-promoting inflammation to antitumor immunity. Cancer Immunol Res. 2014;2:194–9. https://doi.org/10.1158/2326-6066.CIR-13-0214.
Jakowlew SB. Transforming growth factor-beta in cancer metastasis. Cancer Metastasis Rev. 2006;25:435–57. https://doi.org/10.1007/s10555-006-9006-2.
Ikushima H, Miyazono K. TGFbeta signaling: a complex web in cancer progression. Nat Rev Cancer. 2010;10:415–24. https://doi.org/10.1038/nrc2853.
Battle E, Massagué J. Transforming Growth Factor-β signaling in immunity and cancer. Immunity. 2019;50:924–40. https://doi.org/10.1016/j.immuni.2019.03.024.
Xie J, Wang Y, Freeman ME 3rd, Barlogie B, Yi Q. Beta 2-microglobulin as a negative regulator of the immune system: high concentrations of the protein inhibit in vitro generation of functional dendritic cells. Blood. 2003;101:4005–12. https://doi.org/10.1182/blood-2002-11-3368.
Li L, Dong M, Wang XG. The implication and significance of beta 2 microglobulin: a conservative multifunctional regulator. Chin Med J. 2016;129:448–55. https://doi.org/10.4103/0366-6999.176084.
Peter ME, Hadji A, Murmann AE, Brockway S, Putzbach W, Pattanayak A, Ceppi P. The role of CD95 and CD95 ligand in cancer. Cell Death Differ. 2015;22:549–59. https://doi.org/10.1038/cdd.2015.3.
Park S, Cheon S, Cho D. The dual effects of Il-18 in tumor progression. Cell Mol Immunol. 2007;4:329–35.
Esmailbeig M, Ghaderi A. Interleukin-18: a regulator of cancer and autoimmune diseases. Eur Cytokine Netw. 2017;28:l127–40. https://doi.org/10.1684/ecn.2018.0401.
Dillman RO, Nistor GI, McLelland BT, Hsieh C, Poole AJ, Cornforth AN, Keirstead HS. Preliminary observations on soluble programmed death protein-1 (sPD1) as a prognostic and predictive biomarker in patients with metastatic melanoma treated with patient-specific autologous vaccines. Oncotarget. 2019;10:5359–71. https://doi.org/10.18632/oncotarget.27164.
Joiliffe IT, Cadima J. Principal component analysis: a review and recent developments. Philos Trans A Math Phys Eng Sci. 2016;374:20150202. https://doi.org/10.1098/rsta.2015.0202.
Giuliani A. The application of principal component analysis to drug discovery and biomedical data. Drug Discov Today. 2017;22:1069–76. https://doi.org/10.1016/j.drudis.2017.01.005.
David CC, Jacobs DJ. Principal component analysis: a method for determining the essential dynamics of proteins. Methods Mol Biol. 2014;1084:193–220. https://doi.org/10.1007/978-1-62703-658-0_11.
Duforet-Frebourg N, Luu K, Laval G, Bazin E, Blum MG. Detecting genomic signatures of natural selection with principal component analysis: application to the 1000 genomes data. Mol Biol Evol. 2016;33:1082–93. https://doi.org/10.1093/molbev/msv334.
Kim S, Kang D, Huo Z, Park Y, Tseng GC. Meta-analytic principal component analysis in integrative omics application. Bioinformatics. 2018;34:1321–8. https://doi.org/10.1093/bioinformatics/btx765.
Stewart S, Ivy MA, Anslyn EV. The use of principal component analysis and discriminant analysis in differential sensing routines. Chem Soc Rev. 2014;43:70–84. https://doi.org/10.1039/c3cs60183h.
Libbrecht MW, Noble WS. Machine learning applications in genetics and genomics. Nat Rev Genet. 2015;16:321–32. https://doi.org/10.1038/nrg3920.
Gromski PS, Muhamadali H, Ellis DI, Correa E, Turner ML, Goodacre R. A tutorial review: metabolomics and partial least squares-discriminant analysis—a marriage of convenience or a shotgun wedding. Anal Chim Acta. 2015;879:10–23. https://doi.org/10.1016/j.aca.2015.02.012.
Kuligowski J, Pérez-Guaita D, Quintás G. Application of discriminant analysis and cross-validation on proteomics data. Methods Mol Biol. 2016;1362:175–84. https://doi.org/10.1007/978-1-4939-3106-4_11.
Burke S, Lakshmikanth T, Colucci F, Carbone E. New views on natural killer cell-based immunotherapy for melanoma treatment. Trends Immunol. 2010;31:339–45. https://doi.org/10.1016/j.it.2010.06.003.
Curtsinger JM, Schmidt CS, Mondino A, Lins DC, Kedl RM, Jenkins MK, et al. Inflammatory cytokines provide a third signal for activation of naive CD4+ and CD8+ T cells. J Immunol. 1999;162:3256–62.
Dati F, Schumann G, Thomas L, Aguzzi F, Baudner S, Bienvenu J, et al. Consensus of a group of professional societies and diagnostic companies on guidelines for interim reference ranges for 14 proteins in serum based on the standardization against the IFCC/BCR/CAP Reference Material (CRM 470). International Federation of Clinical Chemistry, Community Bureau of Reference of the Commission of the European Communities, College of American Pathologists. Eur J Clin Chem Clin Biochem. 1996;34:517–20. https://doi.org/10.1007/978-1-4939-3106-4_11.
Zhang S. The role of transforming growth factor β in T helper 17 differentiation. Immunology. 2018;155:24–35. https://doi.org/10.1111/imm.12938.
Konforte D, Simard N, Paige CJ. IL-21: an executor of B cell fate. J Immunol. 2009;182:1781–7. https://doi.org/10.4049/jimmunol.0803009.
Croxford AL, Mair F, Becher B. IL23: one cytokine in control of autoimmunity. Eur J Immunol. 2012;42:2263–73. https://doi.org/10.1002/eji.201242598.
Li Y, Yu X, Ma Y, Hua S. IL-23 and dendritic cells: what are the roles of their mutual attachment in immune response and immunotherapy? Cytokine. 2019;120:78–84. https://doi.org/10.1016/j.cyto.2019.02.018.
Sansom DM. CD28, CTLA-4 and their ligands: who does what and to whom? Immunology. 2000;101:169–77. https://doi.org/10.1046/j.1365-2567.2000.00121.x.
Boes M, Cuvillier A, Ploegh H. Membrane specializations and endosome maturation in dendritic cells and B cells. Trends Cell Biol. 2004;14:175–83. https://doi.org/10.1016/j.tcb.2004.02.004.
Trombetta ES, Mellman I. Cell biology of antigen processing in vitro and in vivo. Annu Rev Immunol. 2005;2005(23):975–1028. https://doi.org/10.1146/annurev.immunol.22.012703.104538.
Chow AY, Mellman I. Old lysosomes, new tricks: MHC II dynamics in DCs. Trends Immunol. 2005;26:72–8. https://doi.org/10.1016/j.it.2004.11.008.
Wilson NS, Behrens GM, Lundie RJ, Smith CM, Waithman J, Young L, et al. Systemic activation of dendritic cells by Toll-like receptor ligands or malaria infection impairs cross-presentation and antiviral immunity. Nat Immunol. 2006;7:165–72. https://doi.org/10.1038/ni1300.
Winzler C, Rovere P, Rescigno M, Granucci F, Penna G, Adorini L, et al. Maturation stages of mouse dendritic cells in growth factor-dependent long-term cultures. J Exp Med. 1997;185:317–28. https://doi.org/10.1084/jem.185.2.317.
Romani N, Koide S, Crowley M, Witmer-Pack M, Livingstone AM, Fathman CG, et al. Presentation of exogenous protein antigens by dendritic cells to T cell clones. Intact protein is presented best by immature, epidermal Langerhans cells. J Exp Med. 1989;169:1169–78. https://doi.org/10.1084/jem.169.3.1169.
Hickman-Miller HD, Yewdell JW. Youth has its privileges: maturation inhibits DC cross-priming. Nat Immunol. 2006;7:125–6. https://doi.org/10.1038/ni0206-125.
Gil-Torregrosa BC, Lennon-Duménil AM, Kessler B, Guermonprez P, Ploegh HL, Fruci D, et al. Control of cross-presentation during dendritic cell maturation. Eur J Immunol. 2004;34:398–407. https://doi.org/10.1002/eji.200324508.
Constantino J, Gomes C, Falcäo A, Neves BM, Cruz MT. Dendritic cell-based immunotherapy: a basic review and recent advances. J Immunol Res. 2017;65:798–810. https://doi.org/10.1007/s12026-017-8931-1.
Albert ML, Pearce SF, Francisco LM, Sauter B, Roy P, Silverstein RL, Bhardwaj N. Immature dendritic cells phagocytose apoptotic cells via alphavbeta5 and CD36, and cross-present antigens to cytotoxic T lymphocytes. J Exp Med. 1998;188:1359–68. https://doi.org/10.1084/jem.188.7.1359.
Bonaccorsi I, Campana S, Morandi B, Ferlazzo G. Acquisition and presentation of tumor antigens by dendritic cells. Crit Rev Immunol. 2015;35:349–64. https://doi.org/10.1615/critrevimmunol.v35.i5.10.
Campana S, De Pasquale C, Carrega P, Ferlazzo G, Bonaccorsi I. Cross-dressing: an alternative mechanism for antigen presentation. Immunol Lett. 2015;168:349–54. https://doi.org/10.1016/j.imlet.2015.11.002.
Joffre OP, Segura E, Savina A, Amigorena S. Cross-presentation by dendritic cells. Nat Rev Immunol. 2012;12:557–69. https://doi.org/10.1038/nri3254.
Van Endert P. Intracellular recycling and cross-presentation by MHC class I molecules. Immunol Rev. 2016;272:80–96. https://doi.org/10.1111/imr.12424.
Aarntzen EH, De Vries IJ, Lesterhuis WJ, Schuurhuis D, Jacobs JF, Bol K, et al. Targeting CD4(+) T-helper cells improves the induction of antitumor responses in dendritic cell based vaccination. Cancer Res. 2013;73:19–29. https://doi.org/10.1158/0008-5472.CAN-12-1127.
Dhodapkar KM, Barbuto S, Matthews P, Kukreja A, Mazumder A, Vesole D, et al. Dendritic cells mediate the induction of polyfunctional IL17-producing cells (Th17-1 cells) enriched in the bone marrow of patients with myeloma. Blood. 2008;112:2678–85. https://doi.org/10.1182/blood-2008-03-143222.
Guéry L, Hugues S. Th17 cell plasticity and functions in cancer immunity. Biomed Res Int. 2015;2015:314620. https://doi.org/10.1155/2015/314620.
Chen C, Gao F-H. Th17 cells paradoxical roles in melanoma and potential application in immunotherapy. Front Immunol. 2019;10:187. https://doi.org/10.3389/fimmu.2019.00187.
Guéry L, Dubrot J, Lippens C, Brighouse D, Malinge P, Iria M, et al. Ag-presenting CpG activated pDCs prime Th17 cells that induce tumor regression. Cancer Res. 2014;74:6430–40. https://doi.org/10.1158/0008-5472.CAN-14-1149.
Takeuchi O, Akira S. Pattern recognition receptors and inflammation. Cell. 2010;140:805–20. https://doi.org/10.1016/j.cell.2010.01.022.
Brubaker SW, Bonham KS, Zanoni I, Kagan JC. Innate immune pattern recognition: a cell biological perspective. Annu Rev Immunol. 2015;33:257–90. https://doi.org/10.1146/annurev-immunol-032414-112240.
Dranoff G, Jaffee E, Lazenby A, Golumbek P, Levitsky H, Brose K, et al. Vaccination with irradiated tumor cells engineered to secrete murine granulocyte-macrophage colony-stimulating factor stimulates potent, specific, and long-lasting anti-tumor immunity. PNAS USA. 1993;90:3539–43. https://doi.org/10.1073/pnas.90.8.3539.
Nasi ML, Lieberman P, Busam KJ, Prieto V, Panageas KS, Lewis JJ, et al. Intradermal injection of granulocyte-macrophage colony-stimulating factor (GM-CSF) in patients with metastatic melanoma recruits dendritic cells. Cytokines Cell Mol Ther. 1999;5:139–44.
van de Laar L, Coffer PJ, Woltman AM. Regulation of dendritic cell development by GM-CSF: molecular control and implications for immune homeostasis and therapy. Blood. 2012;119:3383–93. https://doi.org/10.1182/blood-2011-11-370130.
Achuthan A, Cook AD, Lee MC, Saleh R, Khiew HW, Chang MW, et al. Granulocyte macrophage colony-stimulating factor induces CCL17 production via IRF4 to mediate inflammation. J Clin Invest. 2016;126:353–66. https://doi.org/10.1172/jc187828.
Burgdorf SK, Claesson MH, Nielsen HJ, Rosenberg J. Changes in cytokine and biomarker levels in patients with colorectal cancer during dendritic cell-based vaccination. Acta Oncol. 2009;48:1157–64. https://doi.org/10.3109/02841860903099964.
Lawson DH, Lee S, Zhao F, Tarhini AA, Margolin KA, Ernstoff MS, et al. Randomized, placebo-controlled, phase III trial of yeast-derived granulocyte-macrophage colony-stimulating factor (GM-CSF) versus peptide vaccination versus GM-CSF Plus peptide vaccination versus placebo in patients with no evidence of disease after complete surgical resection of locally advanced and/or stage IV melanoma: a trial of the Eastern Cooperative Oncology Group-American College of Radiology Imaging Network Cancer Research Group (E4697). J Clin Oncol. 2015;33:4066–76. https://doi.org/10.1200/JCO.2015.62.0500.
Butterfield LH, Zhao F, Lee S, Tarhini AA, Margolin KA, White RL, et al. Immune correlates of GM-CSF and melanoma peptide vaccination in a randomized trial for the adjuvant therapy of resected high-risk melanoma (E4697). Clin Cancer Res. 2017;23:5034–43. https://doi.org/10.1158/1078-0432.CCR-16-3016.
Andrew N. Cornforth, Ph.D., was responsible for the cryopreservation of blood samples that were analyzed in this study.
The Hoag Hospital Foundation provided support for the production of vaccines and conduct of the clinical trial. Funding support for the serum cytokine assays was provided by Caladrius Biosciences, Inc. Funding for statistical analysis was provided by AIVITA Biomedical, Inc. The funding bodies were not involved in the design of the study, the collection, analysis, and interpretation of data, nor in writing the manuscript.
Ethics approval and consent to participate
The randomized clinical trial was conducted per the Declaration of Helsinki after approval by the Western Institutional Review Board (Seattle, WA., ID 2009 0753), and in accord with assurances filed with and approved by the Department of Health and Human Services. All patients gave written informed consent prior to randomization. In addition to Institutional Review Board oversight, conduct of the trial was monitored by an external consultant, and Hoag Hospital’s institutional clinical research oversight committee. The protocol and manufacturing procedures were reviewed by the US Food and Drug Administration in association with BB-IND 5838 for the tumor cell vaccine, 8554 for the dendritic cell vaccine.
Consent for publication
GIN and ROD are both employees of AVITIA Biomedical, Inc.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
All patients: variance of main components for all patients at baseline.
TCV-treated patients: variance of main components.
TCV-treated patients: KMO and Bartlett’s Test.
TCV-treated patients: Variance. Initial Eigenvalues are rotated with Varimax and Kaiser Normalization method.
KMO and Bartlett’s Test.
Variance. Initial Eigenvalues are rotated with Varimax with Kaiser Normalization method.
Wilks’ Lambda test.
First canonical discriminant function explained 90.2% of variance.
Tests of equality of group means (univariate ANOVA).
Wilks’ Lambda is significant for the first function.
First 2 canonical discriminant functions explain 100% of variance.
First 2 canonical discriminant functions used in the analysis explain 100% of variance.
Wilks’ Lambda test of functions.
Cox regression using IgM baseline values.
About this article
Cite this article
Nistor, G.I., Dillman, R.O. Cytokine network analysis of immune responses before and after autologous dendritic cell and tumor cell vaccine immunotherapies in a randomized trial. J Transl Med 18, 176 (2020). https://doi.org/10.1186/s12967-020-02328-6
- Principal component analysis
- Discriminant analysis
- Metastatic melanoma
- Dendritic cells
- Cancer vaccines