Germline Genetic Variants in TEK, ANGPT1, ANGPT2, MMP9, FGF2 and VEGFA Are Associated with Pathologic Complete Response to Bevacizumab in Breast Cancer Patients

Abstract

Background

We previously reported improved pathologic complete response (pCR) in a prospective phase II study using neoadjuvant bevacizumab in combination with chemotherapy compared to chemotherapy alone in breast cancer patients (41% vs. 25%, p = 0.0291). In this study, we queried germline single-nucleotide polymorphisms (SNPs) in angiogenesis-related genes for their impact on pCR and overall survival (OS).

Methods

DNA for genotyping was available from 34 subjects who received bevacizumab in addition to chemotherapy and 29 subjects who did not. Using Illumina® technology, we queried 504 SNPs with a minor allele frequency (MAF) of at least 5%, located in 10 angiogenesis-related genes, for their effect on pCR via logistic regression with an additive-inheritance model while adjusting for race and bevacizumab treatment. SNPs that showed significant associations with pCR were selected for additional characterization.

Results

After adjusting for race and tumor type, patients who had bevacizumab added to their neoadjuvant therapy were found to experience a significantly improved rate of pCR compared to patients who did not (adjusted OR 8.40, 95% CI 1.90–37.1). When patients were analyzed for SNP effects via logistic regression with race and bevacizumab treatment included as covariates, two SNPs in angiopoietin 1 (ANGPT1), six in ANGPT2, three in fibroblast growth factor 2 (FGF2), four in matrix metalloproteinase 9 (MMP9), three in tyrosine kinase, endothelial (TEK) and two in vascular endothelial growth factor A (VEGFA) were associated with pCR (P<0.05). However, when overall survival was considered, there was no difference between treatment groups or between genotypes.

Conclusion

Genetic variability in TEK, ANGPT1, ANGPT2, FGF2, MMP9 and VEGFA is associated with pCR in bevacizumab-treated patients. Consistent with other studies, adding bevacizumab to standard chemotherapy did not impact OS, likely due to other factors and thus, while SNPs in TEK, ANGPT1, ANGPT2, FGF2, MMP9 and VEGFA were associated with pCR, they were not predictive of OS in this patient population.

Trial Registration

ClinicalTrials.gov NCT00203502

About the Authors:

Issam Makhoul

* E-mail: [email protected]

Affiliation: Division of Hematology/Oncology Division, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

ORCID http://orcid.org/0000-0001-8646-5658

Valentina K. Todorova

Affiliation: Division of Medical Genetics, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Eric R. Siegel

Affiliation: Department of Biostatistics, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Stephen W. Erickson

Affiliation: Department of Biostatistics, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Ishwori Dhakal

Affiliation: Department of Biostatistics, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Vinay R. Raj

Affiliation: Division of Medical Genetics, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Jeannette Y. Lee

Affiliation: Department of Biostatistics, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Mohammed S. Orloff

Affiliation: Department of Epidemiology, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Robert J. Griffin

Affiliation: Department of Radiation Oncology, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Ronda S. Henry-Tillman

Affiliation: Division of Breast Surgical Oncology, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Suzanne Klimberg

Affiliations Division of Breast Surgical Oncology, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America, Department of Pathology, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Laura F. Hutchins

Affiliation: Division of Hematology/Oncology Division, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Susan A. Kadlubar

Affiliation: Department of Radiation Oncology, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Introduction

Bevacizumab, a humanized monoclonal antibody that was designed to target vascular endothelial growth factor A (VEGF-A), was first approved for breast cancer in the metastatic setting in 2008 based on the results of the E2100 Intergroup phase III trial comparing paclitaxel with or without the addition of bevacizumab as first-line therapy. [1] The approval was contingent on the results of further studies, which ultimately did not demonstrate significant improvements in overall survival. Thus, FDA approval of bevacizumab for breast cancer was officially withdrawn in 2010 [2], although additional studies of bevacizumab in the adjuvant setting were ongoing. [3] In the neoadjuvant setting, the addition of bevacizumab to chemotherapy results in increased pathologic complete response [(pCR); reviewed in [4]]. While these results are promising, it is apparent that clinical response to bevacizumab is highly variable. For this reason, biomarkers are urgently needed to identify patients most likely to benefit.

Pharmacogenomics studies have reported associations between single nucleotide polymorphisms (SNPs) in VEGF and VEGFR-2 and response in metastatic breast cancer [5–7], but studies of the impact of SNPs on pCR in the neoadjuvant setting have just started being performed. A recent analysis of the randomized phase III GeparQuinto study examined genetic variants within the VEGF pathway genes in relation to response to neoadjuvant bevacizumab. [8] In that study, a total of 15 genes were examined, but the only two common genes investigated in both that study and the current report were KDR and VEGFA. However, a comprehensive examination of germline SNPs in other genes involved in angiogenesis is lacking.

We have previously conducted and reported on a Phase II neoadjuvant trial examining the combination of chemotherapy with bevacizumab. [9] In that trial, we found significant differences in pCR according to tumor types (52% for ductal carcinoma compared to only 10% in non-ductal disease, p = 0.021), receptor status (59% in estrogen/progesterone receptor negative compared to only 27% in patients positive for at least one receptor, p = 0.047) and race (pCR in African Americans (AA) 75% compared to 28% in European Americans (EA), p = 0.007). Subsequent investigations in this patient population revealed that elevated serum levels of tyrosine kinase with immunoglobulin and epidermal growth factor-homology domains 2 (Tie-2) and basic fibroblast growth factor (bFGF) at baseline were seen in subjects who attained pCR. [10] In the current study, we combined two cohorts in order to assess the effects of the SNPs in the angiogenic pathway while adjusting for the effect of bevacizumab.

Materials and Methods

Study populations

Cancer patients–Data from two independent clinical studies were combined for the analyses in this report. One was a Phase II neoadjuvant clinical trial of bevacizumab added to chemotherapy, for which the study design, sample collection, and protein analysis have been described elsewhere. [9,10] The other was a prospective observational study investigating doxorubicin toxicity in breast-cancer patients who followed the same neoadjuvant chemotherapy regimen as the Phase II trial, but without the addition of bevacizumab; these patients were enrolled between 12/16/2010 and 9/1/2013. [11] Both studies were conducted at the University of Arkansas for Medical Sciences, and both had IRB approval. Serial blood samples were collected at baseline and at each visit from patients in both studies. From the Phase II trial, we excluded one Asian and one Hispanic to leave 37 patients who were either AA or EA. All 38 subjects from the observational study were AA or EA.

Healthy, cancer-free subjects–Healthy subjects were those recruited as controls in case-control studies of breast, colorectal and prostate cancer conducted at UAMS from 1997–2007. These studies were approved by the IRB at UAMS. All participants have signed a written consent form that was approved by the IRB of UAMS. Germline DNA was isolated from buffy coat samples provided by the healthy individuals.

Genotyping assays

High-density genotyping–To test for the association between SNPs in candidate genes and pCR, DNA was extracted from buffy coats of 63 study participants (34 receiving bevacizumab/29 not). Genotyping was performed using the Illumina® Whole Genome Genotyping Infinium chemistry on both cohorts, with the HumanOmni2.5-8v1-1 BeadChip for the cohort receiving bevacizumab, and the HumanOmni5-4 BeadChip for the cohort not receiving bevacizumab. Individual samples with genotype call rates less than 99% and SNPs with call rates less than 95% were removed. The modest sample size precluded the ability to detect genome-wide associations, so even though many other genotypes were generated, we used a candidate gene/SNP approach for analysis.

Candidate SNP genotyping–Specific SNPs were selected for genotyping due to their degree of linkage with other SNPs in the gene under investigation. Genotype was determined using TaqMan® SNP Genotyping Assays (Applied Biosystems) for rs2515462 and rs1960669 using 7900 HT Fast Real-Time PCR system (Applied Biosystems, Life Technologies, Carlsbad, California). Reactions were heated to 95°C for 10 minutes, and then subjected to 40 cycles of amplification at 95°C for 10 seconds and 60°C for 1 minute. PCR amplification was followed by allelic discrimination plate read and analysis.

Statistical analysis

SNP Selection: SNPs were targeted for selection if they resided in one of 10 genes (ANGPT1, ANGPT2, FGF2, IL1A, KDR, MMP9, PDGFB, PECAM1, TEK, and VEGFA) that code for the proteins of our previous study. [10] This covered 833 SNPs on the HumanOmni2.5-8v1-1 BeadChip, and 1460 SNPs on the HumanOmni5-4 BeadChip system. SNPs were eligible for inclusion if they had matching SNP ID names plus matching chromosome + base positions on both chip systems. Differences in reported strand usage between cohorts were reconciled to match the so-called “dbsnp format–forward strand coding scheme” employed for the cohort that received bevacizumab. Eligible SNPs were included if they had a 100% call rate and an observed minor-allele frequency (MAF) >5% across all 63 genotyped samples; 504 SNPs qualified for inclusion into the initial genotype screen.

Statistical Analysis: Patient demographics and tumor characteristics were compared for differences between cohorts using the Wilcoxon rank-sum and chi-square tests. Prior to SNP assessment, the multivariate effect of cohort, race, and tumor type on pCR were evaluated via logistic regression, since race and tumor type showed significant associations with pCR in the Phase II trial of bevacizumab. [9] Because this evaluation showed that tumor type was not independently significant (P = 0.16) in the presence of cohort and race, it was ignored in the subsequent SNP assessments. For the initial genotype screen, the 504 SNPs were individually assessed for their effect on the odds of achieving pCR using a logistic-regression model with cohort and race included as covariates. The model assumed an additive mode of inheritance for the number of minor alleles carried by each subject. A SNP was selected for further analysis if the logistic regression returned an unadjusted P<0.05 for its effect on pCR, although its FDR-adjusted P-value was calculated for reporting purposes. Selected SNPs, in turn, had their minor-allele counts similarly assessed at P<0.05 via logistic regression for association with suitably dichotomized tumor characteristics (stage, grade, ER/PR status, HER2-Neu expression, triple-negative status, and histology). These logistic regressions included Race as an adjustment factor for every characteristic except histology because of the 100% prevalence of ductal carcinoma among African Americans in our two study cohorts. To analyze the selected SNPs for MAF imbalance between cohorts, subjects were subgrouped by race, and the MAF of each SNP was compared between cohorts within race via Fisher’s exact tests. Finally, for two of the 20 selected SNPs, the genotypes of both cohorts were combined and compared via chi-square tests for MAF differences with a race-matched cohort of healthy, cancer-free subjects.

Results

Study population demographics and characterization of pCR

The characteristics of the study population are detailed in Table 1. There were 37 patients who received additional bevacizumab with neoadjuvant chemotherapy, 33 of whom had genotype information available. Similarly, there were 38 patients who had standard neoadjuvant therapy, 29 of whom also had genotypes. The two groups differed significantly at P<0.05 with respect to race, age, tumor size, tumor type, HER2/Neu status, and lymph-node positivity (Table 1). Additionally, the number (%) of patients staged to a Stage III subgroup in Table 1 was 18 (49%) in the chemo + bevacizumab cohort compared to only 8 (21%) in the standard-chemotherapy cohort (P = 0.012). Table 1 demonstrates that patients in the bevacizumab cohort attained pCR much more frequently than patients in the standard chemo cohort. As shown in Fig 1, tumor-related characteristics associated with pCR among genotyped study subjects included hormone receptor negativity, but not ductal tumor type. Table 2 shows the multivariate effect of cohort, race, and tumor type on the odds of attaining pCR. The effect of cohort was independently significant (P = 0.0050), and the effect of race was marginally so (P = 0.051), but the effect of tumor type was insignificant (P = 0.16) in the presence of the other two covariates. As a consequence, tumor type was removed from the multivariate model; Table 3 shows the result. Both cohort and race had independently significant associations with pCR, and both were accordingly included as covariates to adjust for in the initial genotype screen.

[Figure omitted. See PDF.]

Fig 1. Rates of pCR by Treatment cohort, Race, Hormone-receptor status, and Tumor type among genotyped subjects.

Stacked barchart showing how the number and proportion of pCRs varies by treatment cohort, by race, by ER/PR Status, and by tumor type. Group labels appear below their bars. P-values are from chi-square tests comparing groups for pCR differences.

https://doi.org/10.1371/journal.pone.0168550.g001

[Figure omitted. See PDF.]

Table 1. Patient demographics,tumor characteristics, and response to treatment.

https://doi.org/10.1371/journal.pone.0168550.t001

[Figure omitted. See PDF.]

Table 2. Multivariate evaluation of cohort, race, and tumor type on pCR.

https://doi.org/10.1371/journal.pone.0168550.t002

[Figure omitted. See PDF.]

Table 3. Effect of cohort and race on pCR after removing tumor type.

https://doi.org/10.1371/journal.pone.0168550.t003

Initial genotyping screen

Twenty of the 504 SNPs assessed in the initial genotype screen yielded unadjusted P-values <0.05. They included two SNPs in ANGPT1, six in ANGPT2, three in FGF2, four in MMP9, three in TEK and two in VEGFA (Table 4). Although none of the 20 SNPs retained significance after adjusting for multiple comparisons (FDR-adjusted P = 0.88), they were selected for further analysis. Tables 5 and 6 show the minor-allele frequencies of the 20 selected SNPs in the two breast-cancer cohorts receiving neoadjuvant chemotherapy. There were no significant differences in the distribution of alleles between the Bevacizumab cohort and the chemotherapy-alone cohort among AA subjects (Table 5), but significant differences in FGF2 and TEK allele frequencies were evident in EA subjects (Table 6). However, there was no significant difference in pCR for EA patients between the two cohorts (6 pCR /22 EA with bevacizumab versus 2 pCR /24 EA with chemotherapy alone; Fisher’s exact P = 0.13), suggesting that differences in allele frequency were due to chance, given the small number of subjects. None of the other 16 SNPs were significantly different between cohorts in EA subjects. We also examined these 20 SNPs in relation to tumor characteristics (Table 7). One SNP in TEK (kgp3874749) was negatively associated with higher stage of disease while another SNP in FGF2 (kgp11559855) was positively associated with both double-negative ER/PR status and triple-negative ER/PR/HER2-Neu status. Two closely linked SNPs in VEGFA were positively associated with ductal carcinoma. No SNP was associated with Grade or with HER2-Neu expression as a single factor.

[Figure omitted. See PDF.]

Table 4. Odds ratios (95% CI) for pCR by Genotype.

https://doi.org/10.1371/journal.pone.0168550.t004

[Figure omitted. See PDF.]

Table 5. SNP allele distribution by cohort in African American (AA) subjects.

https://doi.org/10.1371/journal.pone.0168550.t005

[Figure omitted. See PDF.]

Table 6. SNP Allele Distribution by cohort in European American (EA) subjects.

https://doi.org/10.1371/journal.pone.0168550.t006

[Figure omitted. See PDF.]

Table 7. SNPs associated with patient characteristics.

https://doi.org/10.1371/journal.pone.0168550.t007

Allele frequency distribution between cancer patients and healthy subjects

Response to Bevacizumab in the clinical trial varied significantly by race. Thus, we examined differences in allele frequency by race (data not shown). We also selected two representative SNPs to examine whether they displayed differences in allele frequencies between the breast cancer patients and healthy individuals from the same population. There were no differences in allele frequencies compared to healthy individuals, suggesting that the two selected SNPs are involved with treatment response and are not risk alleles for breast cancer (Tables 8 and 9).

[Figure omitted. See PDF.]

Table 8. Comparison in European Americans (EA) of allele frequencies between cancer patients and healthy population.

https://doi.org/10.1371/journal.pone.0168550.t008

[Figure omitted. See PDF.]

Table 9. Comparison in African Americans (AA) of allele frequencies between cancer patients and healthy population.

https://doi.org/10.1371/journal.pone.0168550.t009

Discussion

Since the introduction of anti-angiogenic therapy targeting VEGF to clinical practice in the late 1990s-early 2000s, noted differences in efficacy have been observed between cancer types and between patients within the same cancer type. The efficacy of bevacizumab as an anti-angiogenic agent for breast cancer has generated mixed findings. A recent meta-analysis, however, reports that in the neoadjuvant setting, the addition of Bev to standard chemotherapy significantly improves pCR, particularly in patients who are HER2-negative and hormone receptor-negative. [12] While these tumor characteristics are good predictors of pCR, the incorporation of other variables, such as germline SNPs in the angiogenic pathway, could potentially improve prediction models.

Our overarching hypothesis was that the interaction of the type of cancer and the host’s angiogenic profile explains differences in clinical response to VEGF inhibition. In order for a cancer to be sensitive to the interruption of angiogenesis, it has to be undergoing active blood vessel remodeling at the time of treatment. The maturation of blood vessels in many cancers is so advanced by the time of the diagnosis that targeting the machinery for angiogenesis may, at best, lead to stabilization of the tumor rather than its eradication. Clearly, tumor characteristics associated with response to the addition of bevacizumab to chemotherapy are hormone-receptor negativity and ductal type. These cancers are usually highly proliferative, high-grade cancers that are responsive to chemotherapy. Defining the host’s angiogenic profile is more problematic. The last four decades of angiogenesis research have unraveled the multiplicity and redundancy of pathways involved in this crucial process for multicellular organisms. [13] Hence, targeting one pathway when the host may offer the tumor alternative pathways to rescue angiogenesis is not likely to succeed. Alternatively, success is more likely if the targeted pathway is the only one–or the main one–for a specific host.

Analysis of angiogenesis-related proteins in our patients’ serum at baseline showed that high serum Tie2 and high serum bFGF were associated with pCR. [10] Baseline levels of serum ANG1 and ANG2 and their ratios did not show a significant association with pCR. Their levels remained stable during treatment regardless of the response, which suggested that these proteins were related to the host and not the tumor. [10]

In this study, we found that specific SNPs in some of the host angiogenic genes were associated with pCR in the neoadjuvant setting in breast cancer patients receiving bevacizumab and chemotherapy. In our previous study, some of the proteins encoded by these genes were elevated in patients who achieved pCR after receiving bevacizumab; this was the case for bFGF and Tie2, proteins encoded by bFGF and TEK genes, respectively. However, it is interesting that some of these proteins (Tie2) and genes (TEK, ANGP1 and 2) are components of the same pathway. Other pathways may also be involved, including MMP9 and VEGFA. Results from the randomized phase III GeparQuinto study identified five SNPs in the host VEGF-pathway genes that were associated with improved pCR, particularly in patients with triple negative disease. [8] However, there was no overlap in the 15 genes tested in that study with variants investigated in our study except for the KDR and VEGFA genes. The strongest association with pCR in the GeparQuinto study was reported for rs833058, which is located upstream of the VEGFA promoter, whereas in our study, rs833068 and rs833069 were significantly associated with pCR. These SNPs are located within approximately 10,000bp of each other and are likely to be linked. Taken together, these reports indicate that genetic variation in VEGFA can influence response to neoadjuvant bevacizumab. This difference between our study and the GeparQuinto study could be due to differences in patient populations and to the size of the studies. Additionally, the GeparQuinto study reported interactions between VEGFA SNPs and hormone receptor status, while our study indicated associations between VEGFA SNPs and ductal carcinoma, and FGF2 SNPs and receptor status. Again, these differences are likely due to sample sizes and population differences.

Results from the two studies combined for analysis in this report indicate that pCR is improved with the addition of bevacizumab. When overall survival was examined, there was no difference between the cohorts of our study. While this may appear as a shortcoming for this work, in reality, it is likely due to the fact that survival was extremely high (>95%) in both groups (data not shown). What is intriguing in our study is the observation that AA women had much improved pCR compared to women of European ancestry. Whether this finding could be due to a race-related difference in the tumor’s responsiveness to bevacizumab is deserving of further study to fully evaluate if ANGPT1, ANGPT2, TEK, FGF2, MMP9 and VEGFA are promising targets for future research. Genotyping analysis results suggest that the interaction of the host’s angiogenic profile and the type of cancer may contribute to differences in clinical response to VEGF inhibition. Therefore, it is possible that variants of the ANGPT1, ANGPT2, TEK, FGF2, MMP9 and VEGFA genes or their respective proteins might render certain tumors more susceptible to targeting of VEGF. These findings are hypothesis-generating and should be validated in a larger ethnically diverse population. The authors have initiated the process to access DNA from SWOG S0800, a randomized open-label phase II clinical trial that compared the combination of weekly nab-paclitaxel and bevacizumab followed by dose-dense doxorubicin and cyclophosphamide (AC) with nab-paclitaxel followed or preceded by AC as neoadjuvant treatment for HER2-negative locally advanced breast cancer. The primary endpoint was pCR. A total of 215 patients were accrued including 32% with TNBC. The addition of bevacizumab significantly increased the pCR rate overall (36 vs. 21%; p = 0.019) and in TNBC (59 vs. 29%; p = 0.014), but not in hormone receptor-positive disease (24 vs. 18%; p = 0.41). [14] Baseline and on-treatment blood samples were collected. The results of this trial are consistent with our findings and make it a perfect setting to validate our hypothesis.

Supporting Information

[Figure omitted. See PDF.]

S1 File. Dictionary.

AngioSNPs data dictionary. This Excel file contains the Data Dictionary for the ten other Excel files included as Supporting Information; they are described below.

https://doi.org/10.1371/journal.pone.0168550.s001

(XLSX)

S2 File. Clinical.

Clinical Data for Both Cohorts. This Excel file contains the clinical information on the 63 genotyped cohorts in the study. It includes each subject’s cohort, race, tumor type, receptor status, and whether or not she attained pCR.

https://doi.org/10.1371/journal.pone.0168550.s002

(XLSX)

S3 File. Bev Samples.

BEV Ten-gene Sample Genotypes. This Excel file contains the genotype information for the 34 subjects in the Bevacuzumab cohort. The genotype information consists of the alleles determined for every SNP inside the 10 target genes that were included on Illumina’s HumanOmni2.5-8v1-1 BeadChip.

https://doi.org/10.1371/journal.pone.0168550.s003

(XLSX)

S4 File. Bev Anno.

BEV Ten-gene SNP Annotations. This Excel file contains the SNP annotations for all SNPS inside the 10 target genes that were included on Illumina’s HumanOmni2.5-8v1-1 BeadChip.

https://doi.org/10.1371/journal.pone.0168550.s004

(XLSX)

S5 File. Dox Samples.

DOX Ten-gene Sample Genotypes. This Excel file contains the genotype information for the 29 subjects in the Doxorubicin cohort. The genotype information consists of the alleles determined for every SNP inside the 10 target genes that were included on Illumina’s HumanOmni5-4 BeadChip.

https://doi.org/10.1371/journal.pone.0168550.s005

(XLSX)

S6 File. Dox Anno.

DOX Ten-gene SNP Annotations. This Excel file contains the SNP annotations for all SNPS inside the 10 target genes that were included on Illumina’s HumanOmni5-4 BeadChip for all SNPS inside the 10 target genes.

https://doi.org/10.1371/journal.pone.0168550.s006

(XLSX)

S7 File. Dox Strand.

DOX Ten-gene SNP Strand Info. This Excel file contains the SNP strand-usage information for the genotype data in the Doxorubicin cohort (S5 File). It is needed to reconcile the strand usages in the Doxorubicin cohort with those in the Bevacizumab cohort after combining the genotype data for both cohorts.

https://doi.org/10.1371/journal.pone.0168550.s007

(XLSX)

S8 File. Reconciled.

Both Cohorts Reconciled Strands. This Excel file contains the results after combining the genotype data from the two cohorts and making strand usage in the Doxorubicin cohort match that in the Bevacizumab cohort at any SNP where strand usage originally differed between cohorts. These results were used as inputs to generate the data shown in the S9, S10, and S11 Files.

https://doi.org/10.1371/journal.pone.0168550.s008

(XLSX)

S9 File. AlleleFreqs.

Allele Frequencies of All SNPs. This Excel file contains the observed allele frequencies of all SNPs of interest in the two cohorts combined. The results in this file were used to filter the SNPs, i.e., to exclude from further analysis any SNP whose observed minor-allele frequency was less than 5% among the 63 subjects in the analysis.

https://doi.org/10.1371/journal.pone.0168550.s009

(XLSX)

S10 File. FinalSNPs.

Both Cohorts Final SNPs (bases). This Excel file contains the genotype data (i.e., the 126 base calls per SNP) for the SNPs whose observed minor-allele frequencies were greater than 5% among the 63 subjects in the analysis.

https://doi.org/10.1371/journal.pone.0168550.s010

(XLSX)

S11 File. MinorCounts.

Minor-Allele Counts per subject. This Excel file contains, for each SNP in the analysis, the number of minor alleles carried by each subject in the analysis. This file and S2 File Clinical together contain the data used in the logistic regressions.

https://doi.org/10.1371/journal.pone.0168550.s011

(XLSX)

S12 File. Logistic Clinical.

Logistic Regressions vs Cohort, Race, and Tumor Chars. This Word document contains the results of logistic regression of pCR on Cohort, Race, and each of several tumor characteristics taken one at a time, using the 63 subjects who had been genotyped. No SNPs were included.

https://doi.org/10.1371/journal.pone.0168550.s012

(DOCX)

Author Contributions

1. Conceptualization: IM SAK.

2. Data curation: ERS SWE ID VRR.

3. Formal analysis: ERS SWE ID VRR JYL MSO RJG.

4. Funding acquisition: IM VKT.

5. Investigation: VKT SAK.

6. Methodology: IM SAK.

7. Project administration: IM.

8. Resources: RSHT SK LFH.

9. Software: ERS.

10. Supervision: IM.

11. Validation: IM ERS SAK.

12. Visualization: IM.

13. Writing – original draft: IM.

14. Writing – review & editing: IM ERS SAK.

Citation: Makhoul I, Todorova VK, Siegel ER, Erickson SW, Dhakal I, Raj VR, et al. (2017) Germline Genetic Variants in TEK, ANGPT1, ANGPT2, MMP9, FGF2 and VEGFA Are Associated with Pathologic Complete Response to Bevacizumab in Breast Cancer Patients. PLoS ONE 12(1): e0168550. https://doi.org/10.1371/journal.pone.0168550

References

1. Miller KD, Chap LI, Holmes FA, Cobleigh MA, Marcom PK, Fehrenbacher L, et al. Randomized phase III trial of capecitabine compared with bevacizumab plus capecitabine in patients with previously treated metastatic breast cancer. J Clin Oncol. 2005;23: 792–799. pmid:15681523

2. Lyman GH, Burstein HJ, Buzdar AU, D'Agostino R, Ellis PA. Making genuine progress against metastatic breast cancer. J Clin Oncol. 2012;30: 3448–3451. pmid:22927527

3. Westbrook K, Stearns V. Pharmacogenomics of breast cancer therapy: an update. Pharmacol Ther. 2013;139: 1–11. pmid:23500718

4. Kumler I, Christiansen OG, Nielsen DL. A systematic review of bevacizumab efficacy in breast cancer. Cancer Treat Rev. 2014;40: 960–973. pmid:24909311

5. Schneider BP, Miller KD. Angiogenesis of breast cancer. J Clin Oncol. 2005;23: 1782–1790. pmid:15755986

6. Schneider BP, Radovich M, Miller KD. The role of vascular endothelial growth factor genetic variability in cancer. Clin Cancer Res. 2009;15: 5297–5302. pmid:19706811

7. Schneider BP, Shen F, Miller KD. Pharmacogenetic biomarkers for the prediction of response to antiangiogenic treatment. Lancet Oncol. 2012;13: e427–36. pmid:23026828

8. Hein A, Lambrechts D, von Minckwitz G, Haberle L, Eidtmann H, Tesch H, et al. Genetic variants in VEGF pathway genes in neoadjuvant breast cancer patients receiving bevacizumab: Results from the randomized phase III GeparQuinto study. Int J Cancer. 2015;137: 2981–2988. pmid:26100253

9. Makhoul I, Klimberg VS, Korourian S, Henry-Tillman RS, Siegel ER, Westbrook KC, et al. Combined Neoadjuvant Chemotherapy With Bevacizumab Improves Pathologic Complete Response in Patients With Hormone Receptor Negative Operable or Locally Advanced Breast Cancer. Am J Clin Oncol. 2013.

10. Makhoul I, Griffin RJ, Siegel E, Lee J, Dhakal I, Raj V, et al. High-circulating Tie2 is Associated With Pathologic Complete Response to Chemotherapy and Antiangiogenic Therapy in Breast Cancer. Am J Clin Oncol. 2014.

11. Todorova VK, Makhoul I, Siegel ER, Wei J, Stone A, Carter W, et al. Biomarkers for Presymptomatic Doxorubicin-Induced Cardiotoxicity in Breast Cancer Patients. PLoS One. 2016;11: e0160224. pmid:27490685

12. Cao L, Yao GY, Liu MF, Chen LJ, Hu XL, Ye CS. Neoadjuvant Bevacizumab plus Chemotherapy versus Chemotherapy Alone to Treat Non-Metastatic Breast Cancer: A Meta-Analysis of Randomised Controlled Trials. PLoS One. 2015;10: e0145442. pmid:26717149

13. Yadav L, Puri N, Rastogi V, Satpute P, Sharma V. Tumour Angiogenesis and Angiogenic Inhibitors: A Review. J Clin Diagn Res. 2015;9: XE01–XE05.

14. Nahleh ZA, Barlow WE, Hayes DF, Schott AF, Gralow JR, Sikov WM, et al. SWOG S0800 (NCI CDR0000636131): addition of bevacizumab to neoadjuvant nab-paclitaxel with dose-dense doxorubicin and cyclophosphamide improves pathologic complete response (pCR) rates in inflammatory or locally advanced breast cancer. Breast Cancer Res Treat. 2016;158: 485–495. pmid:27393622