Abstract
Background:
Several studies have reported the impact of single nucleotide polymorphisms (SNPs) in vascular endothelial growth factor (VEGF) pathway genes on the efficacy of bevacizumab in metastatic colorectal cancer (mCRC), but results are still inconsistent. The PRODIGE 9 phase III study compared bevacizumab maintenance versus observation alone after induction chemotherapy with FOLFIRI plus bevacizumab.
Objective:
We evaluated the impact of SNPs of VEGF-A, VEGF receptors (VEGFR-1, VEGFR-2), and hypoxia inducible factor-1α (HIF-1α) on tumor control duration (TCD), overall survival (OS), progression-free survival (PFS), and duration of first chemotherapy free-intervals (CFI).
Patients and methods:
We included 314/491 patients from PRODIGE 9 with a DNA blood sample available. Nine SNPs were genotyped on germline DNA using real-time Polymerase Chain Reaction TaqMan TM (Thermo Fisher Scientific, Waltham, MA , USA 02451).
Results:
In the bevacizumab arm, patients with the VEGFR-1 rs9582036 CC genotype (n = 14) had significantly longer TCD [22.4 months (95% confidence interval (CI): 14.75-not reached)] than patients with the AA or CA genotype [14.4 months (95% CI: 11.7–17.1)] (p = 0.036), whereas there was no significant difference in the observation arm. In the bevacizumab arm, no significant difference was found between the CC, and AA or CA genotype for OS [28.2 (95% CI: 18.1–42.8) versus 22.5 (95% CI: 18.6–24.6) months, p = 0.5], PFS [9.4 (95% CI: 7.2–11.3) versus 9.2 (95% CI: 8.71–10.1)], and duration of the first CFI [4.6 (95% CI: 1.6–13.3) versus 4.14 (95% CI: 0.5–29.0) months, p = 0.3].
Conclusion:
Among mCRC patients treated with bevacizumab maintenance, those with the VEGFR-1 rs9582036 CC genotype experienced longer TCD. The presence of this genotype may thus predict a benefit of bevacizumab maintenance in mCRC.
Keywords
Introduction
The prognosis of patients with metastatic colorectal cancer (mCRC) has been significantly improved by the combination of chemotherapy with targeted therapies, antiangiogenic agents (bevacizumab),1,2 or anti-epidermal growth factor receptor agents.3,4
Angiogenesis is a cardinal process that leads to the invasiveness and metastasis of solid tumors. 5 Vascular endothelial growth factor-A (VEGF-A) is a major agent of the angiogenesis pathway. Fixation on VEGF receptor-1 and 2 (VEGFR-1 and VEGFR-2) leads to intracellular pathway activation and promotes vascular permeability, cell proliferation, survival and, migration, tubular morphogenesis and sprouting.5,6 Angiogenesis is a tumor response related to a hypoxic micro-environment. Hypoxia inducible factor-1α (HIF-1α) is a mediator of hypoxia signaling that plays a role in the VEGF-pathway. 7
While the RAS mutation can predict anti-EGFR resistance, 8 there are currently no predictive factors for bevacizumab’s efficacy. A pharmacogenetic study in patients with single nucleotide polymorphisms (SNPs) could be interesting to identify predictive genetic biomarkers on VEGFR genes. Germline SNPs are not dependent on tumor evolution. They are easy to identify and reliable, and can be exploited in a clinical setting with confidence.
Some studies suggest that SNPs on VEGF-A, VEGFR-1, VEGFR-2, or HIF-1α modulate protein expression, and could predict the efficacy of bevacizumab.9–15
The PRODIGE 9 study is a prospective multicenter phase III study comparing maintenance treatment with bevacizumab monotherapy to observation alone after 6 months of induction chemotherapy in mCRC. The results of the study revealed no difference between these two arms for all endpoints. 16
We conducted an ancillary and exploratory study of PRODIGE 9 to evaluate the impact of SNPs on VEGF-A, VEGFR-1, VEGFR-2, and HIF-1α on tumor control duration (TCD), overall survival (OS), progression free survival (PFS), and duration of the first chemotherapy-free interval (CFI).
Patients and methods
Patients’ eligibility criteria
PRODIGE 9 was an open-label, randomized, multicenter, phase III study promoted by the Fédération Francophone de Cancérologie Digestive (FFCD) and the Partenariat de Recherche en Oncologie DIGEstive (PRODIGE) intergroup in 66 French centers. PRODIGE 9 included patients with mCRC without previous chemotherapy or antiangiogenic therapy for metastatic disease. 16 PRODIGE 9 aimed to compare the TCD achieved with first-line chemotherapy followed by either bevacizumab maintenance or observation without treatment during the first CFI. Both arms received induction FOLFIRI plus bevacizumab for 12 cycles. At the end of induction, patients with disease control began a CFI. During the CFI, patients were either treated with bevacizumab monotherapy (maintenance arm) or received no antitumor treatment (observation arm). Chemotherapy was reintroduced at progression for a further eight cycles, followed by one or several new CFI. After progression, the choice of second-line and additional treatments was left to the discretion of the investigator. A blood sample was collected at enrolment of the patient in the trial. All patients from PRODIGE 9 with a DNA blood sample available were eligible for this study.
Follow-up
The tumor response was evaluated by clinicians in accordance with mRECIST 1.1 criteria every 8 weeks.
Polymorphism selection, genotyping
Genes were selected based on their involvement in VEGF-A-dependent angiogenesis. VEGF-A is directly targeted by bevacizumab. Two receptors VEGFR-1 and VEGFR-2 mediate VEGF-A signaling. HIF-1α mediates hypoxic induction of the VEGF-A gene.
All SNPs were shown in previous studies9,13–15,17–19 to play a prognostic and/or predictive role, as well as to have a functional impact and a relative frequency (minor allele frequency) over 5% (Table 1).
SNPs studied, genetic location, ancestral allele and MAF.
SNPs, single nucleotide polymorphism; MAF, minor allele frequency; VEGF, vascular endothelial growth factor; VEGFR, VEGF receptor.
DNA was extracted from blood samples using an automated QIA symphony DSP DNA Mini Kit TM (Qiagen, Hilden, Germany).
DNA was quantified by fluorometry (QubitTM Thermo Fisher Scientific, Waltham, MA, USA 02451) and genotyping was performed with TaqMan® endpoint PCR in a 10 µL reaction volume (15 ng of DNA in 3 µL, 0.5 µL of probe and 5 µL of Master Mix). Negative controls were systematically added, and around 10% of samples were repeated to validate the genotyping procedure.
Outcomes
TCD was defined as the time between randomization and strategy failure, defined as tumor progression during a chemotherapy sequence. 20 Patients whose tumors did not progress during a chemotherapy sequence were censored at their last follow-up within the protocol or censored at the initiation of second-line chemotherapy or another therapeutic strategy. Secondary end points were OS, defined as the time between randomization and death from any cause. Alive patients were censored at the date of the last news. PFS was defined as the time between randomization and first progression or death from any cause. Duration of the first CFI was defined as the time between the end of induction chemotherapy or the date of the last injection for patients in the maintenance arm and the first reintroduction of chemotherapy. From this cohort, all patients for whom a blood sample was available were included. The PRODIGE 9 study was approved by the Committee for the Protection of Persons Ile de France VIII on 12 July 2011 and by the French national agency (AFFSAPS) on the 25 July 2011. The trial was registered on clinicaltrials.gov with the number: NCT00952029. Written informed consent was obtained from all patients before treatment and for the recovery of tumor blocks as requested by the Helsinki declaration (1964) and its amendment (2000).
Statistical analysis
The objectives of the study were to evaluate the impact of SNPs on TCD, OS, PFS, and duration of the first CFI. All SNPs were checked for Hardy–Weinberg equilibrium (HWE) by comparing observed and expected allelic distributions using a chi-square test (χ2 < 3.84 with a significance threshold at 5% and 1 degree of freedom) or Fisher’s exact test for groups smaller than five.
Due to the exploratory design of the study, there was no pre-specified hypothesis tested regarding the effects of a specific genotype and a particular treatment outcome. Multiple statistical tests were performed and no adjustments were made.
Due to small numbers of patients (<4) in some genotype groups: HIF-1 α rs11549465 c.1744C > T and VEGF-A rs3025039 c.*237C > T, patients homozygous for the rare allele and heterozygous were combined for analyses.
Tumor response was analyzed using an univariate correlation test (Fisher’s exact or chi square with a threshold at 0.05).
The survival analysis (TCD, OS, and PFS) were performed using the Kaplan–Meier model with the log rank test (significance threshold 0.05).
We included significant genotypes (p ⩽ 0.05 from univariate analyses) in a single multivariate analysis using the Cox model taking into account independent prognostic factors recognized in the literature, that is, WHO performance status, number of metastatic sites, and age at occurrence of hypertension. As nine SNPs were tested, p values were corrected for multiple testing using the false discovery rate controlling procedure of Benjamini and Liu. 21 Adjusted p values were computed using the discrete module type package for R 3.1 (R Foundation for Statistical Computing, Vienna, Austria). All other statistical tests were performed with IBM SPSS v20.0 (IBM Inc., New York, USA).
Results
From March 2010 to July 2013, 494 patients were randomized in the PRODIGE 9 study. 16 DNA was available for 314 patients: 160 in the bevacizumab maintenance arm and 154 in the observation arm during the first CFI (Figure 1). These 314 patients were comparable to the overall population of PRODIGE 9.

Study flow chart.
The results are presented for the 314 patients enrolled in this study. In the overall population after the first induction chemotherapy, 156 (52.5%) patients had a complete or partial response, 133 (44.8%) patients had stable disease, and eight (2.7%) had disease progression. Patients characteristics are summarized in Table 2.
Population characteristics, TCD, OS, and PFS in each arm (arm A: bevacizumab maintenance; arm B: observation arm without treatment during CFIs).
CFI, chemotherapy-free interval; CI, confidence interval; ECOG PS, Eastern Cooperative Oncology Group performance status; OS, overall survival; PFS, progression-free survival; SD, standard deviation; TCD, tumor control duration.
Median TCD during CFI in the bevacizumab maintenance arm was similar to that in the observation arm. Median TCD, median OS, median PFS, and median duration of the first CFI in the 314 patients were similar to those of the main study.
The genotype distribution of the nine SNPs agreed with those predicted by HWE. Observed genotype frequencies are summarized in Table 3.
Frequencies of genotypes by SNPs.
CFI, chemotherapy-free interval; IQR, interquartile range; OS, overall survival; PFS, progression-free survival; TCD, tumor control duration; VEGFR, vascular endothelial growth factor receptor.
Associations between SNPs and efficacy endpoints
A statistically significant association was found between TCD and VEGFR-1 rs9582036. TCD was significantly longer [median 22.4 months (95% confidence interval (CI) : 14.7–NR)] in patients with the CC genotype than in patients with the AA and CA genotypes [14.4 months (95% CI : 11.7–17.1) (hazard ratio = 0.39; 95% CI : 0.16–0.97) p = 0.04] in the bevacizumab maintenance arm, as shown in Figure 2(a). In the observation arm, no significant association was found between this SNP and TCD, as shown in Table 4.

Kaplan–Meier curves according to study treatment for VEGFR-1 rs9582036. (a) TCD. (b) OS.
Efficacy endpoint for the VEGFR-1 rs9582036.
CFI, chemotherapy-free interval; IQR, interquartile range; OS, overall survival; PFS, progression-free survival; TCD, tumor control duration; VEGFR, vascular endothelial growth factor receptor.
In the bevacizumab maintenance arm, there was a trend toward better OS in patients with the CC genotype for VEGFR-1 rs9582036 versus those with the AC + AA genotype, although the difference was not significant [28.2 (95% CI : 18.1-42.8) versus 22.5 (95% CI : 18.6–24.6) months, p = 0.5), as shown in Figure 2(b). In the observation arm, the median OS was similar in both the CC and AC + AA genotype subgroups.
Moreover, no statistically significant difference was observed for either PFS or first CFI duration for CC patients versus AC and AA patients (Table 4).
There was no significant association between any genotype and objective response after the first induction chemotherapy. However, there was a trend toward a better tumor response in patients with the VEGFR-1 rs9582036 CC genotype than that in AC + AA patients, as shown in Table 5.
Best response on the first induction treatment for VEGFR-1 rs9582036.
VEGFR, vascular endothelial growth factor receptor.
There was no other statistically significant association between other SNPs evaluated and any endpoints, as described in Supplemental Tables 1 and 2.
Discussion
Bevacizumab is one of the main biotherapies for mCRC. In the era of ‘personalized medicine’, it has become essential to determine predictive factors of response to bevacizumab so as to identify the subgroups of patients that may benefit from this treatment without added toxicity.
In our study, we studied nine SNPs involved in the VEGF-pathway in a large series of 314 patients included in the PRODIGE 9 phase III study comparing bevacizumab maintenance versus observation alone after induction first-line chemotherapy with FOLFIRI plus bevacizumab.
Despite the randomized design of this trial set up to determine whether patients benefited or not from bevacizumab maintenance therapy, no clinical factors consistently predicted bevacizumab efficacy in this exploratory study. However, we did find that the VEGFR-1 rs9582036 genotype had an impact on survival in patients treated with maintenance bevacizumab. TCD was significantly longer in patients with the CC genotype than in those with the CA and AA genotype. We also found a trend toward better OS in patients with the CC genotype. Nonetheless, no prolongation of the CFI was observed, suggesting that the maintenance of antiangiogenic pressure could be pointless for this genotype. To date, this is the largest genotyping study of SNPs involved in the VEGF pathway in mCRC patients treated with maintenance bevacizumab as a single agent, and the first to report long-term follow-up data regarding differences in outcomes between the CC and the CA + AA genotype for VEGFR-1 rs9582036. However, the VEGFR-1 rs9582036 CC genotype was rare, found in 14 patients in the bevacizumab arm and 16 patients in the observation arm. The potential clinical value is unquestionable; but the data should still be validated in an independent series.
Hansen et al. reported a similar low frequency of the VEGFR-1 rs9582036 CC genotype, present in 22 of the 218 mCRC patients analyzed. However, they showed an inverse association, with a better response in patients with the VEGFR-1 rs9582036 AA genotype versus CC + AC. 22 Chionh et al. reported no predictive association between VEGFR-1 rs9582036 and survival in patients with mCRC receiving bevacizumab plus chemotherapy. In their study, however, tumor tissue and not blood was used for the SNP analysis. 23 Loupakis et al. also reported no predictive impact of the VEGFR-1 rs9582036 genotype on survival outcomes.13,24
Several studies have suggested some associations between remaining SNPs and bevacizumab efficacy in mCRC, but no correlation was strong enough to be used in clinical practice.13–15,17,22,24–27 A prognostic impact was found for VEGF-A rs1570360, 15 VEGF-A rs2010963, 15 VEGFR-2 rs125057758, 13 and VEGFR-2 rs2305948. 14 Previous data suggested the possible role of the VEGF-A rs833061 genotype in predicting the efficacy of bevacizumab,24,27 but failed to show a predictive impact. 13 VEGF-A rs699947 and VEGF-A rs3025039 have been found to be associated with a lower prevalence of bevacizumab-induced hypertension. 25
In other kinds of cancers, several studies have linked VEGFR-1 rs9582036 with clinical outcomes in chemotherapy plus bevacizumab-treated patients. In a retrospective study of 77 patients with pancreatic cancer, the VEGFR-1 rs9582036 CC genotype was a predictor of shorter survival in patients treated with bevacizumab, and correlated with increased VEGFR-1 expression. 28 The same was shown for patients with non-small-cell lung cancer treated with bevacizumab in two small studies: shorter survival correlated with higher baseline plasma levels of VEGF-A and the presence of VEGFR-1 rs9582036 CC.29,30 Conversely, some studies showed that clinical outcomes in patients with VEGFR-1 rs9582036 AA were better than outcomes in those with CC + AC, especially in patients with breast cancer treated with bevacizumab-based therapy. 31 In patients with renal cancer treated with sunitinib, a multitarget antiangiogenic agent, VEGFR-1 rs959582036 CC was associated with lower objective tumor response, PFS and OS. 32
More studies are needed to understand the molecular mechanism of VEGFR-1 polymorphism in response to VEGF blockage.
VEGFR-1 rs9582036 has been shown to be in high linkage disequilibrium with rs7993418. Fine-mapping experiments of the VEGFR-1 locus identified rs7993418, a synonymous SNP affecting tyrosine 1213 in the VEGFR-1 tyrosine-kinase domain, as the functional variant underlying the association. This SNP causes a shift in codon usage, leading to increased VEGFR-1 expression and downstream VEGFR-1 signaling. We could infer that CC patients in our study could be more sensitive to bevacizumab because of that variation of the VEGR-1 level. An approach with simultaneous VEGF-A dosing could be interesting to corroborate this hypothesis.
Some limitations must be considered. First, similar to other investigations in this field and due to the exploratory design, no pre-specified hypothesis was tested regarding the effects of a specific genotype and a particular treatment outcome, we conducted multiple statistical tests, and no adjustments were made.
Then, bevacizumab’s mechanisms of action are not fully known. Bevacizumab not only blocks VEGF-A, but also plays a role in the tumor microenvironment, including the inhibition of bone marrow-derived progenitors, the normalization of vessel structure, vascular ‘constriction’, disruption of the cancer stem cells niche, a direct effect on tumor cells, and interactions with the host immune system. This may explain the difficulty to identify a single ‘candidate SNP’. A broader genome-wide approach to identify more relevant SNPs could be interesting. This method was used by Kim, 19 who conducted genome-wide SNP screening on high-speed DNA chips (SNP array) in blood samples from colon cancer patients. Surprisingly, genes unrelated to angiogenesis were found to be significantly linked to the bevacizumab response. However, the population studied was heterogeneous: patients were treated with bevacizumab or cetuximab, an anti-EGFR which acts via a different mechanism.
Other prognostic factors have been identified in the PRODIGE 9 study, such as the BRAF V600E mutation, 33 baseline splenic volume, 34 and radiomic signature, 35 but these were not evaluated in correlation with SNP in this study. In recent studies, the side of the colorectal cancer, left or right, has also been reported as a prognostic factor and a predictor of response to treatment,36,37 with a worst prognosis in right-side colorectal cancer.36,38 This prognostic effect was not confirmed in other prospective trials4,39 including PRODIGE 9 16 and should be interpreted with caution.
It is important to note that our analyses were limited by the small number of VEGFR-1 rs9582036 CC patients, as in other studies. A larger validation cohort with more CC patients is needed to confirm the predictive value of VEGFR-1 rs9582036 for bevacizumab treatment. Moreover, the pharmacokinetics of bevacizumab was not taken into account, even though it has been suggested that clinical outcomes vary according to the concentration of bevacizumab. 40 A global approach with pharmacokinetic dosing could be interesting to take into account this possible cofounding factor. In addition, a VEGFR-1 expression analysis seems necessary to understand how VEGFR-1 rs9582036 CC affects protein function.
In the PRODIGE 9 trial, bevacizumab monotherapy did not impair tumor progression during CFI. The application of this findings regarding the VEGFR-1 s9582036 CC genotype in patients treated with bevacizumab maintenance in clinical practice seems to be extremely limited. Nevertheless, a pooled analysis of several trials including PRODIGE 9. That evaluated bevacizumab monotherapy during CFI reported an improvement in PFS in patients treated with bevacizumab monotherapy. 41 It would be interesting to evaluate VEGFR-1 rs9582036 in these patients to select patients that could potentially benefit from bevacizumab monotherapy during CFI.
In conclusion, our study found that VEGR-1 rs9582036 CC was associated with prolonged TCD for patients with mCRC on bevacizumab-based maintenance therapy. This result needs to be confirmed in larger studies, which should include pharmacokinetics and expression analyses to understand and strengthen the potential predictive relevance of VEGR-1 rs9582036 in patients treated with bevacizumab.
Supplemental Material
sj-docx-1-tam-10.1177_17588359221141307 – Supplemental material for Predictive value of vascular endothelial growth factor polymorphisms for maintenance bevacizumab efficacy in metastatic colorectal cancer: an ancillary study of the PRODIGE 9 phase III trial
Supplemental material, sj-docx-1-tam-10.1177_17588359221141307 for Predictive value of vascular endothelial growth factor polymorphisms for maintenance bevacizumab efficacy in metastatic colorectal cancer: an ancillary study of the PRODIGE 9 phase III trial by Bernadette de Rauglaudre, Camille Sibertin-Blanc, Aurélie Fabre, Karine Le Malicot, Jaafar Bennouna, François Ghiringhelli, Julien Taïeb, Valérie Boige, Olivier Bouché, Thierry Chatellier, Roger Faroux, Eric François, Stéphane Jacquot, Dominique Genet, Claire Mulot, Sylviane Olschwang, Jean-François Seitz, Thomas Aparicio and Laetitia Dahan in Therapeutic Advances in Medical Oncology
Supplemental Material
sj-docx-2-tam-10.1177_17588359221141307 – Supplemental material for Predictive value of vascular endothelial growth factor polymorphisms for maintenance bevacizumab efficacy in metastatic colorectal cancer: an ancillary study of the PRODIGE 9 phase III trial
Supplemental material, sj-docx-2-tam-10.1177_17588359221141307 for Predictive value of vascular endothelial growth factor polymorphisms for maintenance bevacizumab efficacy in metastatic colorectal cancer: an ancillary study of the PRODIGE 9 phase III trial by Bernadette de Rauglaudre, Camille Sibertin-Blanc, Aurélie Fabre, Karine Le Malicot, Jaafar Bennouna, François Ghiringhelli, Julien Taïeb, Valérie Boige, Olivier Bouché, Thierry Chatellier, Roger Faroux, Eric François, Stéphane Jacquot, Dominique Genet, Claire Mulot, Sylviane Olschwang, Jean-François Seitz, Thomas Aparicio and Laetitia Dahan in Therapeutic Advances in Medical Oncology
Footnotes
Acknowledgements
The authors thank all participating patients and their families. We thank all the study groups, investigators from the participating countries. The authors also thank the team from the FFCD data center.
Declarations
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
