Prognostic or predictive value of circulating cytokines and angiogenic factors for initial treatment of multiple myeloma in the GIMEMA MM0305 randomized controlled trial

Abstract

Background

Several new drugs are approved for treatment of patients with multiple myeloma (MM), but no validated biomarkers are available for the prediction of a clinical outcome. We aimed to establish whether pretreatment blood and bone marrow plasma concentrations of major cytokines and angiogenic factors (CAFs) of patients from a phase 3 trial of a MM treatment could have a prognostic and predictive value in terms of response to therapy and progression-free and overall survival and whether these patients could be stratified for their prognosis.

Methods

Blood and bone marrow plasma levels of Ang-2, FGF-2, HGF, VEGF, PDGF-β, IL-8, TNF-α, TIMP-1, and TIMP-2 were determined at diagnosis in MM patients enrolled in the GIMEMA MM0305 randomized controlled trial by an enzyme-linked immunosorbent assay (ELISA). These levels were correlated both reciprocally and with the type of therapy and patients’ characteristics and with a group of non-MM patients as controls.

Results

No significant differences were detected between the blood and bone marrow plasma levels of angiogenic cytokines. A cutoff for each CAF was established. The therapeutic response of patients with blood plasma levels of CAFs lower than the cutoff was better than the response of those with higher levels in terms of percentage of responding patients and quality of response.

Conclusion

FGF-2, HGF, VEGF, and PDGF-β plasma levels at diagnosis have predictive significance for response to treatment. The stratification of patients based on the levels of CAFs at diagnosis and their variations after therapy is useful to characterize different risk groups concerning outcome and response to therapy.

Trial registration

Clinical trial information can be found at the following link: NCT01063179

Background

Multiple myeloma (MM) is the second most common hematologic cancer after non-Hodgkin’s lymphoma, with a higher incidence in the elderly. Patients older than 70 years account for 56% of new cases and for 73% of all deaths from MM [1]. Combined melphalan-prednisone has been the standard of care for more than 40 years and has been found to be associated with a median survival of 29 to 37 months [2,3,4]. Today, the availability of novel agents, such as the first-in-class proteasome inhibitor bortezomib and the immunomodulatory drugs thalidomide and lenalidomide, has significantly improved the clinical outcome of these patients [5,6,7,8,9,10,11,12,13].

Accurate identification of high-risk patients and risk stratification are crucial in improving outcomes for MM patients, but considerable heterogeneity exists in their overall survival. Although a large number of prognostic markers have been described, including disease burden (Durie-Salmon staging system, International Staging System, magnetic resonance imaging, (18F)fluorodeoxyglucose positron emission tomography, presence of extramedullary disease or plasma-cell leukemia), host factors (age, performance status, and renal function), tumor biology (proliferation rate, conventional cytogenetics, interphase fluorescence in situ hybridization, and gene expression profiling), and depth of response to therapy [14,15,16,17,18], none of them completely explains the heterogeneity seen in this tumor.

To further complicate matters, some of the new treatments appear to overcome the high risk defined by one or more of these prognostic factors [19, 20]. With the increased treatment options, the ability of some treatments to overcome certain risk factors, and the availability of markers to define risk categories, risk stratification in the management of MM is becoming an important issue [21]. The achievement of a uniform risk stratification system would also allow a better comparison of patient groups across different trials [22].

Angiogenesis is a constant hallmark of MM progression and has prognostic potential. The pathophysiology of MM-induced angiogenesis involves both direct production of angiogenic cytokines by plasma cells and their induction within the bone marrow microenvironment [23]. Moreover, inhibitors of the vascular endothelial growth factor (VEGF) pathway, including bevacizumab and tyrosine kinase inhibitors such as sorafenib, sunitinib, and pazopanib, have been shown to prolong progression-free survival (PFS) and overall survival (OS) [24,25,26,27,28] and are in fact approved for the treatment of solid cancer.

It has been previously demonstrated that the plasma levels of cytokines and angiogenic factors (CAFs) decrease after therapy in patients with cancer, and this may be relevant for treatment response and PFS [29,30,31,32]. Here, we demonstrate that high levels of CAFs are negative prognostic factors in patients with MM and seem to be predictive of relative benefit from therapy. Moreover, the stratification of patients based on CAF levels at diagnosis is useful to detect different risk groups for outcome and response to therapy.

Methods

Patients

Patient characteristics are reported in Table 1. This study has been carried out on MM patients enrolled in the multicenter clinical trial GIMEMA-MM0305, with the participation of 61 centers in Italy from May 2006 to January 2009. The study compared the combination bortezomib-melphalan-prednisone-thalidomide followed by maintenance with bortezomib-thalidomide (VMPT-VT) with bortezomib-melphalan-prednisone (VMP) administered for nine cycles without maintenance. The details and results of the trial have been published previously [33,34,35]. Clinical protocol and informed consent documents were approved by the participating local institution’s review boards, and the trial was undertaken in accordance with the International Conference on Harmonization Guidelines for Good Clinical Practice and the amended Declaration of Helsinki. Patients without MM or other tumors (patients with stage I arterial hypertension without organ damage and without other diseases) who gave their consent were used as controls.

Table 1 Patients’ characteristics

Methods

Before starting the treatment, peripheral blood and bone marrow plasma (the initial 1 ml of bone marrow aspirate) samples were collected into EDTA-containing tubes. Both blood and bone marrow plasma samples from 124 of 511 patients enrolled in the study (24%) were available for analysis. Plasma was separated by centrifugation (2,000 rpm for 20 min at 4 °C) within 1 h from blood drawing and aliquoted into multiple cryovials. Plasma samples were stored at − 80 °C until use. Before analysis, plasma samples were thawed slowly in an ice bath and all analyses were done from a one-off thaw sample. CAFs were measured by using Q-Plex™ Array Human Angiogenesis Antigen (Quansys Biosciences, Logan, Utah) allowing the simultaneous quantification of the following factors: angiopoietin-2 (ANG-2), fibroblast growth factor-2 (FGF-2), hepatocyte growth factor (HGF), interleukin-8 (IL-8), platelet-derived growth factor-BB (PDGF-BB), tissue inhibitor of matrix metalloproteinase-1 and 2 (TIMP-1, TIMP-2), tumor necrosis factor-alpha (TNF-α), and vascular endothelial growth factor (VEGF), according to the manufacturer’s instructions. Secreted levels of CAFs were quantified through Q-View Software (Quansys Biosciences, Logan, Utah) in triplicate samples, and the mean results were used in biomarker analysis.

Statistical analysis

In the first step, CAF levels were measured in the blood and bone marrow plasma samples of MM patients to assess whether significant differences could be detected in the two compartments by Student’s t test (p values less than 0.05 was considered significant). Since no significant differences between the two compartments were detected (see the “Results” section), subsequent analyses were carried out using plasma samples from the peripheral blood. With the Student t test, the CAF levels in the plasma samples of the MM patients were compared with those of controls.

In the second phase, plasma levels of CAFs were measured as an independent variable to predict binary response status (≥ VGPR vs < VGPR) by logistic regression analysis. CAF plasma levels were also correlated as a continuous variable with tumor response by linear regression and logarithmic transformation to normalize CAF values. The correlation between log-CAFs and tumor response was approximately linear. Selection of individual CAF markers from the screening phase was done on the basis of results of median cutoff, ROC curve estimation of cutoff, and logistic regression analysis between dichotomized tumor response and CAF plasma levels. We assessed the association between CAF plasma levels and progression-free survival (PFS) with the Cox proportional hazard model.

To establish whether plasma levels of any CAF might have prognostic or predictive significance, the Kaplan-Meier method was used to analyze PFS and OS. We used the Cox regression model to verify significant differences noted in Kaplan-Meier curves for both treatment groups between a high- and a low-CAF subgroup, defined by the respective median CAF plasma levels. Sensitivity analyses confirmed that median cutoff achieved the most significant segregation of clinical benefit. To assess the potential differential effects of baseline CAF concentrations between two treatment groups, a treatment versus CAF status interaction term was included in the Cox model analysis for PFS, with treatment group and CAF status as two additional independent variables. CAF plasma levels with a significant interaction value with treatment were regarded as predictive. A post hoc analysis was done to adjust for multiple testing of CAF markers with the Bonferroni test. Exploratory analyses included correlation between CAF plasma levels and Eastern Cooperative Oncology Group (ECOG) performance status, Durie & Salmon stage, International Staging System stage, Cytogenetic risk, and age and sex hierarchical clustering (unweighted pair group method with arithmetic mean, unweighted average, and Euclidean distance for similarity measure) to assess a multiple-CAF signature association with PFS or OS (Kaplan-Meier method for PFS and OS, and Cox regression models to assess differences). All statistical analyses were done with SPSS software.

Results

Biological samples (blood and bone marrow) from 124 MM patients randomly assigned to receive VMPT-VT (53 pts) or VMP (71 pts) in the GIMEMA MM0305 phase III clinical trial and blood plasma samples from 54 control subjects were available for this study. Baseline demographic and disease characteristics are reported in Table 1.

We evaluated the concentration of ANG-2, FGF-2, HGF, IL-8, PDGF-BB, TIMP-1, TIMP-2, TNF-α, and VEGF that are the major cytokines involved in angiogenesis in MM and other cancers [36] and, as previously demonstrated [37] in MM patients, directly correlate with disease activity and increase with progression. Moreover, plasma levels of CAFs are directly related to disease response to therapy in hematologic and solid tumors [29,30,31,32].

Our results showed that there were no differences in the levels of the studied CAFs between the peripheral blood and bone marrow plasma samples of MM patients (Fig. 1, Additional file 1: Table S1), indicating that the concentrations of circulating cytokines well reflect those of the bone marrow and could be used for all subsequent analyses. As expected, the plasma levels of CAFs in MM patients were significantly higher compared to controls (p < 0.0001 for all CAFs). In addition, the levels of CAFs were found to be significantly related to MM response to therapy (Fig. 2) with the exception of TIMP-1 and TIMP-2. More precisely, low levels of ANG-2 (p < 0.05), FGF-2 (p < 0.005), HGF (p < 0.05), IL-8 (p < 0.05), PDGF-BB (p < 0.005), TNF-α (p < 0.05), and VEGF (p < 0.005) were indicative of more profound response (very good partial response [VGPR] or better in all patients, with no evident differences between the two therapeutic regimens (VMPT-VT vs VMP: p = 0.1).

Fig. 1

Analysis of the CAF levels in blood and bone marrow plasma samples of MM patients. No differences are evident in their concentration between peripheral blood and bone marrow samples. Significantly higher levels of CAFs are detected in blood and bone marrow samples of MM patients as compared with control subjects (p < 0.0001 for all cytokines)

Fig. 2

Response rate of MM patients based on CAF levels. The blood levels of CAFs significantly correlate with MM response to therapy. Lower levels of ANG-2 (p < 0.05), FGF-2 (p < 0.005), HGF (p < 0.05), IL-8 (p < 0.05), PDGF-BB (p < 0.005), TNF-α (p < 0.05), and VEGF (p < 0.005) are indicative of more profound response, VGPR or better, in all patients, with no evident differences between the two therapy regimens (VMPT-VT vs VMP: p = 0.1)

On the basis of ROC curve estimation of cutoff and logistic regression analysis between dichotomized tumor response and CAF plasma levels, a cutoff for each cytokine was established that could be used to discriminate the probability of response to therapy of MM patients with high sensitivity and specificity (Additional file 1: Table S2). Seventy-four patients had the best degree of tumor response (≥ VGPR) and were termed good responders, whereas 50 had the smallest degree of tumor response (< VGPR) and were termed poor responders. Among the seven CAFs that were shown to be significantly related to MM response to therapy (Fig. 2), low (relative to median) concentrations of ANG-2 (p < 0.0001), FGF-2 (p < 0.0001), HGF (p < 0.0001), IL-8 (p < 0.0001), PDGF-BB (p < 0.0001), TNF-α (p < 0.001), and VEGF (p < 0.0001) highly correlated with best response (Additional file 1: Table S2).

In terms of survival, PFS correlated with CAF levels as a dichotomous variable, and low levels of FGF-2 (p < 0.0001), HGF (p < 0.05), IL-8 (p < 0.05), TNF-α (p < 0.05), and VEGF (p < 0.005) were associated with better PFS. Instead, only low levels of FGF-2 (p < 0.001) and VEGF (p < 0.004) were associated with prolonged OS. Similar results were obtained when CAF levels were related to PFS and OS according to the combination therapy administered. A better PFS was in fact related to low levels of FGF-2 (p < 0.0001), HGF (p < 0.0001), TNF-α (p < 0.005), and VEGF (p < 0.005) in the VMP arm and FGF-2 (p < 0.005) and VEGF (p < 0.05) in the VMPT-VT arm. When OS was taken into consideration, only FGF-2 (p < 0.001 in the VMP arm, p < 0.005 in the VMPT-VT arm), and VEGF (p < 0.05 for both arms) were demonstrated (Additional file 1: Table S3).

Hierarchical clustering analysis (Fig. 3) showed three distinct risk groups of patients, based on the concentrations of two CAF levels (FGF-2 and VEGF). Patients with elevation of both CAF levels had a worse prognosis with significantly shorter PFS and OS (high risk) compared with those with high level of only one (p < 0.0001, intermediate risk) and those of low blood levels of both CAFs (p < 0.0001, low risk) (Fig. 3). No significant differences were detected depending on the type of therapy received. Trying to construct a survival model, a three-stage system using FGF-2 and VEGF blood levels provided the highest statistically significant results (Table 2). Median survivals of the risk groups were as follows: low risk: PFS 44 months, OS 70 months; intermediate risk: PFS 23.5 months, OS 62 months; and high risk: PFS 14.5 months, OS 34 months (p < 0.0001 for differences). Patient numbers were well distributed across the three groups (low risk, 40%; intermediate risk, 29%; and high risk, 31%).

Fig. 3

Progression-free and overall survival analysis in MM patients based on the peripheral blood plasma concentrations of FGF-2 and VEGF. The hierarchical clustering analysis of MM patients shows three distinct risk groups of patients based on the concentrations of FGF-2 and VEGF. High risk: patients who present both high FGF-2 and VEGF plasma levels showing a worse prognosis with significantly shorter PFS and OS; intermediate risk: patients who present high plasma levels in only one of the two cytokines; low risk: patients who show low blood levels in both angiogenic cytokines. Again, no evident differences between the two therapy regimens were detected

Table 2 Patient stratification based on CAF circulating levels

By multivariate analysis (Table 3), the other variables significantly associated with better outcome were age < 65 years in the intermediate- and high-risk groups (HR 0.61, CI 0.44–0.85, p = 0.0047 and HR 0.60, CI 0.38–0.89, p = 0.0391, respectively) for OS; age < 65 years and best response to induction therapy (≥ VGPR) only in the high-risk group (HR 0.68, CI 0.48–0.96, p = 0.0342) for PFS. Durie and Salmon (D&S) stage 1 or 2 impacted only on PFS for the high-risk group (HR 0.64, CI 0.49–0.89, p = 0.0169). No significant differences were demonstrated regarding both PFS and OS for sex, isotype of the M-component, renal failure, and type of induction therapy. The administration of maintenance therapy significantly impacted on PFS in all the risk groups (low: HR 0.70, CI 0.44–0.89, p = 0.0210; intermediate: HR 0.54, CI 0.38–0.80, p = 0.0008; high: HR 0.71, CI 0.56–0.92, p = 0.0249) and OS (low: HR 0.70, CI 0.48–0.88, p = 0.0210; intermediate: HR 0.74, CI 0.57–1.01, p = 0.0330; high: HR 0.52, CI 0.36–0.75, p = 0.0003).

Table 3 Multivariate analysis of risk stratification of patients based on the CAFs FGF-2 and VEGF

As regards the cytogenetic risk, because many data relative to the cytogenetic characteristics of the patients were missing in the database of the trial, the statistical power of the relative analysis was low. For this reason, it is not possible to reach a correct conclusion on the value of CAFs stratification risk in correlation to this risk parameter.

Discussion

The emergence of new treatment options for MM has extended the patients’ survival [5] and the need to prospectively identify those patients who are likely to benefit from a specific treatment and understand the mechanisms underlying therapeutic resistance. Several adverse prognostic factors have been identified in MM at diagnosis and before initiation of treatment [38, 39], including an advanced stage in the international staging system (ISS) based on plasma albumin [40]. The dual activity of the new and newest drugs active both on MM plasma cells and bone marrow stromal cells, and in particular on angiogenesis [41,42,43], obviously indicates that any new prognostic markers cannot ignore the angiogenesis aspect. Previous studies have indeed suggested that the plasma levels of CAFs might be used to identify prognostic and predictive markers in solid tumors [44,45,46,47,48,49] and are indicative of the response to antineoplastic therapy [29,30,31,32].

Seven CAFs of the nine evaluated (ANG-2, FGF-2, HGF, IL-8, PDGF-BB, TNF-α, and VEGF) initially emerged as being related to disease activity, but further testing showed that only FGF-2 and VEGF were significantly associated with PFS and OS and were therefore evaluated for prognostic stratification of patients. We also assessed whether determination of the studied CAFs could add prognostic information to D&S [50] and International Staging Systems (ISS) [40] or whether these systems were, in fact, predictive of benefit. D&S and ISS staging were associated with prognosis based on PFS in both groups; however, they were not strong prognostic parameters as FGF-2 and/or VEGF. Thus, the plasma levels of the two angiogenic cytokines provided prognostic and, more importantly, predictive value beyond that of standard clinical staging.

Our analyses of survival were based on assignment at initial randomization. High FGF-2 and VEGF plasma levels were negative prognostic markers and were associated with lower relative OS in both harms. The associations reported here with FGF-2 and VEGF were akin to previous studies [29,30,31,32], in which patients with the highest reduction of circulating CAF concentrations obtained the greatest benefit from anticancer therapy. However, other studies assessing the use of VEGF as a predictive marker of benefit from VEGF-targeted therapies in renal-cell carcinoma have yielded inconsistent results [51].

Previous studies have shown that groups of related angiogenic or inflammatory factors are often correlated [44, 46, 52]. Hierarchical clustering of six circulating CAFs showed a strong correlation among many of them, including osteopontin and VEGF [46, 52]. Patients defined by high concentrations of these CAF and inflammatory or immunomodulatory factors had a significantly worse prognosis, but derived a greater relative OS benefit from therapy. The circulating CAFs identified in these studies might themselves have important biological roles or might be markers for alternative pathways or mechanisms affecting treatment responses (e.g., hypoxia-induced factor-1α or NF-κB pathways). This association between factors suggests that common mechanisms might regulate their production. Studying a cohort of patients with metastatic renal-cell carcinoma, a classification of the disease based on the higher expression of angiogenic versus inflammatory circulating CAFs defined by a six-cytokine signature was established [46]. In that study, the angiogenic group derived greater benefit from sorafenib alone, whereas the other group benefited from the combination of sorafenib and IFN-α. In the present study, a similar correlation of treatment benefit with circulating CAF signature was observed in MM patients.

All the evidences provided by the studies on the role of microenvironment [3, 17], and particularly of the angiogenic process [23, 37] in myeloma progression as well as in cancer cell protection mediated by microenvironment components [37], indicated that the response to therapy is also related to the activity of anticancer drugs on tumor microenvironment [42, 43]. Therefore, the inhibition of cytokine production which mediates the interaction between cancer cells and their microenvironment represents one of the major goals of the modern therapeutic approaches. On these bases, the evaluation of CAF levels is indicative of the response to therapy and may represent a good indicator of refractoriness in cancer patients.

Conclusion

Overall, our findings support the use of circulating CAF profiling to define biologically distinct subgroups of MM patients whose tumors have a greater angiogenic drive. As such, these patients might have a more aggressive disease course but are likely to derive relative benefit from inhibition of angiogenic pathways. Circulating CAF profiling might be particularly well suited for angiogenesis inhibitors and other drugs targeting the tumor microenvironment, in which both circulating host-derived and tumor-derived factors could affect response. Such an approach may have important advantages, including straightforward and relatively non-invasive sample collection, availability of robust analytical platforms, and the ability to monitor changes during treatment or disease progression, which can help identify markers of resistance.

The limitation of this study is that the results have been obtained in older patients not-eligible for ASCT, whereas the major strength of the study is the homogeneous stratification and longer follow-up of patients. Further studies on a greater cohort of subjects, including young patients eligible for ASCT before and after treatment, will be needed to evaluate the variation of CAFs after therapy, to study their power as an indicator of minimal residual disease and, then, of risk of relapse, and to evaluate the value of the risk stratification based on CAFs in the new drugs era, with the aim of definitively establish the value of this approach in the application of precision, personalized therapy for patients with MM.