A pharmacogenetic interaction analysis of bevacizumab with paclitaxel in advanced breast cancer patients
- Select a language for the TTS:
- UK English Female
- UK English Male
- US English Female
- US English Male
- Australian Female
- Australian Male
- Language selected: (auto detect) - EN

Play all audios:

To investigate pharmacogenetic interactions among VEGF-A, VEGFR-2, IL-8, HIF-1α, EPAS-1, and TSP-1 SNPs and their role on progression-free survival (PFS) in metastatic breast cancer (MBC)
patients treated with bevacizumab plus first-line paclitaxel or with paclitaxel alone. Analyses were performed on germline DNA, and SNPs were investigated by real-time PCR technique. The
multifactor dimensionality reduction (MDR) methodology was applied to investigate the interaction between SNPs. The present study was an explorative, ambidirectional cohort study: 307
patients from 11 Oncology Units were evaluated retrospectively from 2009 to 2016, then followed prospectively (NCT01935102). Two hundred and fifteen patients were treated with paclitaxel and
bevacizumab, whereas 92 patients with paclitaxel alone. In the bevacizumab plus paclitaxel group, the MDR software provided two pharmacogenetic interaction profiles consisting of the
combination between specific VEGF-A rs833061 and VEGFR-2 rs1870377 genotypes. Median PFS for favorable genetic profile was 16.8 vs. the 10.6 months of unfavorable genetic profile (p =
0.0011). Cox proportional hazards model showed an adjusted hazard ratio of 0.64 (95% CI, 0.5–0.9; p = 0.004). Median OS for the favorable genetic profile was 39.6 vs. 28 months of
unfavorable genetic profile (p = 0.0103). Cox proportional hazards model revealed an adjusted hazard ratio of 0.71 (95% CI, 0.5–1.01; p = 0.058). In the 92 patients treated with paclitaxel
alone, the results showed no effect of the favorable genetic profile, as compared to the unfavorable genetic profile, either on the PFS (p = 0.509) and on the OS (p = 0.732). The
pharmacogenetic statistical interaction between VEGF-A rs833061 and VEGFR-2 rs1870377 genotypes may identify a population of bevacizumab-treated patients with a better PFS.
The treatment of metastatic breast cancer (MBC) patients with hormone-receptors positive (HR+) and human epidermal receptor 2 negative (HER2−) tumors is dramatically changed over the years.
In this setting, cyclin-dependent kinase 4/6 inhibitors (CDK4/6i), such as palbociclib, ribociclib and abemaciclib, in combination with aromatase inhibitors or fulvestrant represent today
the first and later lines of therapy1.
However, chemotherapy-based treatment is still a therapeutic choice when hormone resistance occurs, in triple-negative tumor or in case of visceral crisis2,3. In this scenario, the humanized
monoclonal antibody bevacizumab, in combination with paclitaxel, is a treatment option compared to chemotherapy alone4. Although a significant improvement in progression-free survival (PFS)
was observed from three comparative studies, the US Food and Drug Administration (FDA), but not the European Medicines Agency (EMA), revoked the initial approval of bevacizumab for the
first-line treatment of MBC patients because of the lack of benefit in terms of overall survival (OS). However, it has been theorized that when a long survival post first-line progression is
expected after a first-line chemotherapy, such as in breast cancer, the lack of an apparent benefit in OS could not mean a lack of improvement in OS for the first line of
treatment4,5,6,7,8.
Different strategies have been investigated to find possible predictive biomarkers and select those patients with the best chance of response to bevacizumab. Indeed, the PFS improvement due
to bevacizumab was identical for magnitude in all subgroups of patients with different clinical and pathological characteristics9, and therefore new selective biomarkers should be needed to
identify those patients who can have a major advantage in terms of outcome. Despite many attempts have been done, no validated biomarkers are currently available in the clinical practice and
the prospective MERiDiAN trial failed to demonstrate a possible role of VEGF-A baseline in predicting the response to bevacizumab in breast cancer patients10,11,12,13,14,15.
Germline and somatic polymorphisms of genes involved in the angiogenic pathways have also been widely investigated in this research area to predict bevacizumab outcome, with contrasting
results12,16,17,18. Due to the retrospective nature of these studies and to their inconclusive results, the role of single nucleotide polymorphisms (SNPs) as predictive markers remains to
define19. Therefore, the current approach of correlating the bevacizumab response to a single SNP may be replaced by a genetic analysis of the interaction between SNPs, defined as non-linear
interaction or epistasis. Moore and colleagues have established and validated a methodology, called multifactor dimensionality reduction (MDR) analysis, to identify a genetic profile with
the ability to predict the drug response20. To test this hypothesis, we conducted a retrospective study on 113 MBC patients to assess the ability of MDR methodology to identify a favorable
pharmacogenetic profile associated to PFS in patients treated with bevacizumab, combined with first-line paclitaxel. The MDR analysis provided two pharmacogenetic interaction profiles
consisting of the combination between specific VEGFR-2 rs11133360 and IL-8 rs4073 genotypes. The median PFS was 14.1 months (95% CI, 11.4–16.8) and 10.2 months (95% CI, 8.8–11.5) for the
favorable and the unfavorable genetic profile, respectively (HR = 0.44; 95% CI, 0.29–0.66; p C and TYMS-TSER polymorphisms, instead of an individual polymorphism, seemed to predict the
CAPOX-B (capecitabine, oxaliplatin, and bevacizumab combination) response in terms of PFS, suggesting a paradigm shift from SNPs to a more complex interaction gene analysis able to predict
response to antitumor agents.
The current study was planned to evaluate the effects of the combination of paclitaxel with bevacizumab versus paclitaxel alone on MBC patients harboring different genetic profiles,
exploring the possibility to predict the best favorable profile in terms of PFS. The second step was to test if the eventual seen PFS advantage could be maintained also in terms of OS in
these patients. The previous study on VEGFR-2 and IL-8 genetic interaction analysis21, the favorable profile in terms of PFS was not predictive of OS benefit. In the present study, the seen
advantage in PFS was indeed confirmed in OS (an adjusted hazard ratio of 0.71) but with a p = 0.058, a value very close to a statistically significance, but not significant. However, the
reported data, although statistically negative, seem to suggest that the favorable profile in terms of PFS may probably be maintained also in terms of OS and undoubtedly merits further
investigations in a validation prospective study. Indeed, evaluation and confirmation of these findings in an independent cohort is critical because of the exploratory nature of our
ambidirectional trial.
The analyses conducted with the MDR methodology in this unselected population of MBC patients revealed more than a genetic interaction profile, consisting of the combination between specific
genotypes, but, due to nature of the MDR methodology, we investigated the genetic profile with a benefit in terms of both PFS and OS. The analysis conducted revealed a genetic interaction
profile, consisting of the combination between specific genotypes of VEGF rs833061 and VEGFR-2 rs1870377. Particularly, two genetic profiles were identified in patients, as reported in Table
5. The first one was associated with a greater both PFS and OS benefit compared to the second one. However, this model considered all the candidates and allows for any and all combination
of SNPs to correlate with outcome. Thus, there are other significant or borderline permutations. Indeed, we have also included, as an example in the supplementary data (Supplementary Table 2
and Supplementary Fig. 2), another interesting genetic profile with a significant advantage in term of PFS but without any advantage in OS (not even a tendency).
Therefore, in our study we demonstrated, through the MDR methodology, a statistical interaction between VEGF-A and VEGFR-2 gene SNPs that potentially relates to bevacizumab efficacy on both
PFS and OS. The two genes, and, consequently, the two proteins belong to the same signaling pathway, and it has been clearly demonstrated that VEGF-A stimulates VEGFR-2 phosphorylation and
tumor angiogenesis24. Based on these premises, it is conceivable to hypothesize that, in patients carrying the favorable genetic profile, the tumor angiogenesis is successfully inhibited in
the presence of bevacizumab. The pharmacological inhibition of the angiogenic process by bevacizumab could be effective because of the physiological (not increased) production of VEGF-A due
to the presence of VEGF-A rs833061 CC genotype or C allele. Indeed, for this SNP VEGF-A rs833061 C>T it has been described an increased promoter activity due to the T allele25 that may
explain an eventual resistance to the treatment. Moreover, the VEGFR-2 rs1870377 is a nonsynonymous SNP substituting glycine with histidine (Q472H) located in the extracellular ligand
binding region of the receptor, potentially impacting VEGFR-2 degradation26. The VEGFR-2 rs1870377 TT genotype or T allele present in the favorable profile synthetize a receptor not modified
in its structure, suggesting that it is not abnormally activated or degraded. Therefore, it might be plausible that the genetic background characterized by a physiological activation of the
VEGF-A pathway may be responsible, in part, for the positive effect of bevacizumab maintenance therapy in these MBC patients. In contrast, in patients with an unfavorable genetic profile,
the microenvironment conditions due to the different genotype combinations may result in an increase of the VEGF-A production and/or the presence of an altered VEGFR-2 on tumor endothelial
cells which may be capable to proliferate, migrate or survive because the VEGF action is not completely blocked by bevacizumab.
The absence of any advantage in terms of efficacy in the patients treated with chemotherapy without bevacizumab could suggest a possible predictive role of the favorable genetic profile for
bevacizumab response, but the exploratory nature of this ambidirectional study may limit this hypothesis. However, the main findings of our analyses support the conclusion that a genetic
profile may identify a group of patients with longer PFS and OS, predicting the response to bevacizumab in combination with paclitaxel.
The MDR approach is a major reason for differences between our trial and other studies on bevacizumab biomarkers such as E210016. There are additional aspects between the E2100 US patients
and the Italian population of our study that may account for different results. First of all, Italian patients were of Caucasian origin and no patients of African origin were represented.
Secondly, although the frequencies of our studied VEGF-A and VEGFR-2 SNPs were superimposable with the ones of the Caucasian patients published in the article by Schneider and colleagues16,
there was an exception regarding the VEGFR-2 889A/G (rs2071559). In that case, the frequency of the minor allele A in our population was 0.49 whereas in the E2100 study was 0.09.
New pharmacogenetic favorable biomarkers of bevacizumab-combined therapies could be retrieved from a genetic analysis of the interaction among SNPs rather than from the examination of a
single SNP of a single gene. Surely, a multigene-risk biomarkers may be more beneficial from a comprehensive agnostic approach using genome-wide association studies (GWAS) rather than a
candidate gene approach as the one that we have used in our study. However, some challenges have been faced when scientists tried the scaling of MDR to big data, as the one from GWAS, such
as the necessity to filter the data prior to MDR analysis27, also using biological knowledge through tools such as BioFilter28. Moreover, our work can definitively be strengthened by the
biological characterization of the VEGF expression in the pre-treatment tissue. Indeed, since rs833061 is located in the promoter region of VEGF-A, the difference in expression levels of
VEGF-A in tumors of patients harboring the favorable vs. unfavorable profile could be an important strategy to confirm our statistical findings.
In conclusion, the MDR methodology could be successfully used as witnessed by the experience in this unselected MBC patients where the investigation of an interaction between VEGF-A rs833061
and VEGFR-2 rs1870377 gene polymorphisms resulted in the identification of a genetic profile associated with a longer PFS.
This is an explorative, ambidirectional cohort study, meaning that eligible patients were enrolled and evaluated retrospectively from January 2009 until September 2016 and then followed
prospectively. The oncology units, all located in the north or center of Italy, were selected based on their clinical experience in the use of the combination of paclitaxel and bevacizumab
as first-line therapy in histologically confirmed HER-2-negative MBC patients. Two-hundred fifteen patients from 11 Italian divisions of Medical Oncology, with histologically confirmed
HER2-negative MBC, were treated with a first-line therapy including bevacizumab 10 mg/m2 i.v. on days 1 and 15 combined with first-line paclitaxel 90 mg/m2 i.v. on days 1, 8, and 15, every 4
weeks, and they were enrolled for the present pharmacogenetic study. Ninety-two MBC patients treated with a first-line chemotherapy including paclitaxel without bevacizumab, during the same
period, were also included into the study as a control group. The patients enrolled in the previously published study21 have been also included in the present analysis. Basal and
pathological characteristics recorded from both groups were the following: age (≤ or >65 years); Eastern Cooperative Oncology Group (ECOG) performance status (0 or 1–2); hormonal-receptor
status (positive or negative); previous adjuvant chemotherapy (none, anthracycline or anthracycline plus taxanes); previous hormonal therapy (adjuvant or metastatic); disease-free interval
from the first diagnosis of breast cancer (≤ or >12 months); extent of disease (≤ or >3 sites); location of disease (viscera or bone); disease evaluation (measurable or non-measurable).
Patients with human epidermal growth factor receptor type 2 (HER2)-positive, were excluded from the present study. The characteristics of the patients are summarized in Table 1.
The treatment with chemotherapy was continued until either disease progression occurred or unacceptable toxicities registered, or it was stopped for medical choice. The bevacizumab
maintenance was continued, and hormone therapy added for both groups when indicated, until disease progression or unacceptable toxicities occurred.
Sites of metastatic disease were radiologically re-evaluated according to the RECIST criteria 1.1, in patients with measurable disease, every 2 months. In patients without measurable
lesions, progression of disease was defined when new lesions appeared or when existing lesions evolved. Likewise, in the case of non-measurable lesions, deterioration of clinical condition
not due to treatment toxicity, was defined as progression of disease.
PFS was defined as the period from the beginning of the treatment to the first observation of disease progression as above described, or death from any cause. OS was defined as the period
from the beginning of the treatment to death from any cause. All patients were assessed for response, PFS and OS. Each patient entering the study signed the informed consent. The disease
assessment was conducted by the investigators based on the approved protocol and all the oncology units followed the same assessment schedule and criteria for the prospective follow-up. The
protocol was approved by ethic committee of Azienda Ospedaliera-Universitaria Pisana (CESM-AOUP 3077/2010; clinicaltrials.gov identifier NCT01935102) for Pisa, Livorno, Lucca, Massa Carrara,
Versilia, and Pontedera Hospitals, and by the ethic committees of all participating centers.
Blood samples (3 ml) were collected in EDTA tubes and stored at −80 °C. Genes and polymorphisms, involved in the angiogenesis pathways, were selected for the present analyses based on our
previous study21. In the Table 9, the selected polymorphisms are reported. Germline DNA extraction was performed using QIAamp DNA Blood Mini Kit (Qiagen, Valencia, CA, USA). Allelic
discrimination of genes was performed using an ABI PRISM 7900 SDS (Applied Biosystems, Carlsbad, CA, USA) and with validated TaqMan® SNP genotyping assays (Table 9; Applied Biosystems). PCR
reactions were carried out according to the manufacturer’s protocol. Genotyping was not performed until an adequate number of events (>80% on study population) was reported in terms of PFS.
All the samples were analyzed twice to replicate the obtained genotype.
The investigators responsible for data analysis were blinded to which samples were from patients treated with paclitaxel alone and paclitaxel plus bevacizumab.
The aim of the present study was to identify a favorable genetic profile in terms of PFS in MBC patients treated with bevacizumab in association with paclitaxel. The corresponding OS in
these patients remained a secondary endpoint as well as response rate. All polymorphisms were analysed for deviation from the Hardy–Weinberg Equilibrium (HWE) by means of comparison between
observed allelic distributions with those expected from the HWE by on χ2 test (see Supplementary Tables 3 and 4).
Any association between gene polymorphisms and response rate was analysed by the two-sided Fisher’s exact test. The association between each individual polymorphism and the most relevant
clinical-pathological characteristics with PFS and OS was tested using a Cox proportional hazards model. In these analyses we used Bonferroni’s correction and the p value