Abstract
Breast cancer (BC) is one of the most causes of death in women worldwide. It affects Iranian female population approximately a decade earlier than those in other parts of the world. Previous studies have shown that vascular endothelial growth factor (VEGF) gene variants were associated with BC risk. The current study aimed to evaluate the impact of VEGF rs3025039 (+936C>T), rs2010963 (+405C>G), rs833061 (–460T>C), rs699947 (–2578C>A), and rs35569394 (18-bp I/D) polymorphisms on BC risk in an Iranian population in southeast of Iran. This case-control study was done on 250 BC patients and 215 healthy women. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) or PCR was used to genotype the polymorphisms. Our findings showed that VEGF rs699947 variant increased the risk of BC (OR = 1.71, 95% CI = 1.15–2.54, P = 0.009, CA vs CC; OR = 2.12, 95% CI = 1.14–3.93, P = 0.021, AA vs CC; OR = 1.78, 95% CI = 1.22–2.60, P = 0.004, CA+AA vs CC; OR = 1.47, 95% CI = 1.12–1.92, P = 0.005, A vs C). The VEGF rs3025039, rs2010963, rs833061, and rs35569394 variants were not associated with risk/protection of BC. In conclusion, our results proposed that VEGF rs699947 polymorphism may increase the risk of BC development. Furthers studies with larger sample sizes and different ethnicities are necessary to confirm our findings.
Introduction
Breast cancer (BC) is one of the most prevalent malignancies in women worldwide, which affects more than 1 million women annually.1,2 It comprises 21.4% of all cancers among Iranian females. Although the etiology of BC is entirely unknown, there is abundant evidence that genetic factors play key roles in the pathogenesis and progression of BC.3–6
The human VEGF gene is located on the short arm of chromosome 6 (6p12–p21) and consists of eight exons separated by seven introns that exhibit alternative splicing to form a family of proteins. 7 VEGF, also known as VEGFA, is an important regulator of physiological angiogenesis during embryogenesis, skeletal growth, and reproduction functions. It is also involved in pathological angiogenesis associated with tumors.8,9 Angiogenesis is an important step in the development of cancer and is necessary for primary tumor growth, invasiveness, and metastases. 8 Overexpression of VEGF was found in several tumor tissues.10–12 Breast cancer is among the well-known malignancies involving lymphangiogenesis, which is the recruitment of blood and lymphatic vessels, to a growing tumor. 13 Evidence from in vitro and in vivo studies has shown that increased expression of VEGF is associated with tumor growth and metastasis. Moreover, the inhibition of VEGF signaling results in suppression of tumor-induced angiogenesis as well as tumor growth. 14
Single nucleotide polymorphisms (SNPs) in the promoter, intron, exon, or untranslated regions (3′- and 5′-UTR) may affect the production or function of the corresponding protein. A number of SNPs have been designated in the VEGF gene, which have been described to be associated with differential VEGF expression.15–18 Two of these SNPs are located in the VEGF promoter region (–2578 and -1154),15,18 one SNP (+405) in exon 1 of the VEGF gene 16 and a further SNP (+936) located in exon 8, corresponding to the 3′UTR region of the gene. 17 In addition, other SNPs have also been defined, although an association with VEGF expression has not been revealed. 16
Several studies investigated the VEGF genetic polymorphisms in BC in different ethnic groups and led to different conclusions.19–27 To the best of our knowledge, there are no reports concerning the impact of VEGF polymorphisms on BC risk in Iranian women. Hence, this study aimed to evaluate the possible association between VEGF rs3025039 (+936C>T), rs2010963 (+405C>G), rs833061 (–460T>C), rs699947 (–2578C>A), and rs35569394 (18-bp I/D) polymorphisms with risk/protection of BC in an Iranian population.
Patients and Methods
Patients
This population-based case-control study was done on 250 histologically confirmed BC patients admitted to Ali Ebneh Abitaleb Hospital (Iran) and 215 age-matched healthy women with no history of any types of cancer and not related to the patients.
The local ethics committee of Zahedan University of Medical Sciences approved the study, and written informed consent was obtained from all the subjects before enrollment. The research complied with the principles of the Declaration of Helsinki. Blood samples were collected from participants in tubes containing EDTA, and genomic DNA was extracted by the use of salting-out method. 28
Genotyping
Genotyping of VEGF variants was performed using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) or PCR methods. The primers used for detection of polymorphisms are listed in Table 1. PCR was done using commercially available Prime Taq premix (GeNet Bio, South Korea) according to the manufacturer's recommended protocol. In each 0.20-mL PCR reaction tube, 1 μL of genomic DNA (~100 ng/mL), 1 μL of each primer (10 μM), 10 μL of 2X Prime Taq Premix, and appropriate amount of ddH2O were added. The PCR-cycling conditions were 5 minutes at 95°C, followed by 30 cycles of 30 seconds at 95°C, 30 seconds at 62°C for rs2010936, 63°C for rs3025039 and rs833061, 61°C for rs699947, 60°C for rs35569394, and 30 seconds at 72°C with a final extension step of 72°C for 5 min. The 10-μL of PCR product was digested by appropriate restriction enzymes (Table 1). The PCR products were electrophoresed on agarose gel containing 0.5 μg/mL ethidium bromide and visualized on a UV transilluminator.
Primers and restriction enzymes used in PCR-RFLP analysis for detection of VEGF polymorphisms.
To verify genotyping quality for each polymorphism, ~20% of random samples were regenotyped and the results confirmed the preceding genotyping outcomes.
Statistical analysis
The statistical analysis was done with statistical package SPSS version 20. Cases and controls were compared by chi-square test for categorical data or independent Student's t-test for continuous data. The association between each variant and BC was evaluated by computing the odds ratio (OR) and 95% confidence interval (95% CI) from unconditional logistic regression analyses. For the cases, the chi-square test was used to estimate the association between each variant and the clinicopathological characteristics. Haplotype analysis was done using SNPStats software. 29 The statistical level of significance was defined as P < 0.05.
Results
The study group consisted of 250 BC patients with an average age of 48.7 ± 11.0 years and 215 healthy women with a mean age of 50.4 ± 13.1 years. There was no significant difference between the groups concerning age (P = 0.144).
The genotype and allelic frequency of VEGF polymorphisms are listed in Table 2. The findings indicated that VEGF rs699947 variant increased the risk of BC in codominant (OR = 1.71, 95% CI = 1.15–2.54, P = 0.009, CA vs CC; OR = 2.12, 95% CI = 1.14–3.93, P = 0.021, AA vs CC) and dominant (OR = 1.78, 95% CI = 122–2.60, P = 0004, CA+AA vs CC) inheritance models tested. The A allele increased the risk of BC (OR = 1.47, 95% CI = 1.12–1.92, P = 0.005) when compared with C allele. We found that the genotype and allele frequencies of VEGF rs3025039, rs2010963, rs833061, and rs35569394 variants were not associated with risk/protection of BC in our population.
VEGF gene polymorphisms in BC patients and control subjects.
The results of haplotype analysis are given in Table 3. The findings indicated no statistically significant associations of VEGF haplotypes with BC risk.
Haplotype frequencies of VEGF polymorphisms in breast cancer (BC) and control subjects.
The association between clinicopathological characteristics of BC patients with VEGF polymorphisms is shown in Table 4. Patients with rs3025039 CT genotype had a higher risk of developing ER negative (ER–) BC (OR = 2.39, 95% CI = 1.30–4.41, P = 0.007) than those with CC genotype.
Correlation between VEGF polymorphisms and clinicopathological characteristics of breast cancer (BC) patients.
Regarding the rs699947 variant, patients with CA genotype were more prone to develop BC at age >50 years (OR = 2.16, 95% CI = 1.17–4.00, P = 0.0176) than those with CC genotype. Patients with rs699947 AA genotype had a higher risk of BC developing with tumor size >2 cm (OR = 3.88, 95% CI = 1.40–10.74, P = 0.0074) in comparison with CC genotype. Patients with rs699947 CA genotype had a lower risk of developing progesterone receptor negative (PgR-) BC (OR = 0.47, 95% CI = 0.26–0.87, P = 0.020) than those with CC genotype. In addition, patients with rs699947 CA and AA genotype significantly decreased the risk of developing HER2 negative BC (OR = 0.44, 95% CI = 0.24–0.79, P = 0.008 and OR = 0.37, 95% CI = 0.16–0.88, P = 0.031) than those with CC genotype. No significant association between VEGF rs2010963 (+405C>G), rs833061 (–460T>C), and rs35569394 (18-bp I/D) polymorphisms and clinicopathological characteristics of BC patients was found (P > 0.05).
Discussion
The present study investigated the potential impact of VEGF gene polymorphisms on the susceptibility and clinicopathological features of BC in an Iranian population in southeast of Iran. Our findings proposed that VEGF rs699947 variant significantly increased the risk of BC. The results did not support an association between rs3025039, rs2010963, rs833061, and rs35569394 variants of VEGF gene and BC risk.
We performed a stratified analysis by clinicopathological characteristics of BC patients, and our findings proposed that patients with rs3025039 CT genotype had a higher risk of developing ER– BC. Regarding the rs699947 variant, CA genotype increased the risk of developing BC in old individuals (age > 50 years), as well as in patients with tumor size >2 cm. However, patients with CA genotype had a lower risk of developing PgR- and HER negative BC. No other variants were associated with clinicopathological characteristics of the patients.
The –2578 C/C genotype of the VEGF promoter has been associated with higher levels of VEGF in blood mononuclear cells. 18 It has been shown that –2578 C>A variant correlates with disease stages in which angiogenesis plays a critical role.30–32
The results of meta-analysis indicated that rs2010963 (+405C>G) and rs699947 (–2578C>A) variants of VEGF were not associated with BC susceptibility.33,34 The results of a meta-analysis performed by Li and Ju 35 indicated that the functionally important +936 C/T polymorphism of VEGF is not associated with BC risk. The results of meta-analysis indicated that VEGF +936CC genotype and +936C allele increased the risk of BC and tumor growth, especially in Asians. VEGF–634G allele has no effect on BC susceptibility in all subjects, whereas decreases BC susceptibility in Asians and inhibits tumor growth. 36
Kapahi et al 25 investigated the impact of VEGF –2578C/A, –2549I/D, –460T/C, and –7C/T polymorphisms on BC in North Indian population. They found that II genotype of –2549I/D polymorphism significantly increased the risk of BC when compared with DD genotype. VEGF –2578AA genotype and A allele were found to be associated with increased risk of BC. The CC genotype and C allele of VEGF –460T/C variant significantly increased the risk of BC. The +405G>C and –1154G>A variants of VEGF were not associated with BC risk. 20
No significant association between the +936 C/T polymorphism and BC risk was found in Moroccan women. 19 However, the +936 C/T variant was associated with HER2/neu expression. In a meta-analysis performed by Chen et al 21 revealed no significant association between VEGF –2578C/A and the risk of cancer. Subgroup analyses showed that the VEGF –2578C/A polymorphism is associated with colorectal and lung cancers but not with bladder and breast cancers. Furthermore, the polymorphism may decrease the risk of cancer in the Asian population. The results of a meta-analysis showed a significant association between VEGF +936C/T polymorphism and the risk of BC. 26 It has been reported that VEGF –634G/C polymorphism is associated with increased BC risk and aggressiveness. 27
Luo et al 22 investigated the impact of VEGF –634 G/C, +936 C/T, and +1612 G/A polymorphisms on breast cancer in Chinese Han patients. They proposed that individual who carries +936 TT genotypes or 936 T-allele had a protective effect concerning the disease. They found that the -634CC genotype was significantly associated with high tumor aggressiveness. The +1612G/A polymorphism was not associated with BC risk. Rodrigues et al 23 observed that women carriers of +936 CT+TT VEGF genotypes have a protective effect concerning BC. However, rs1109324, rs1547651, and rs833052 variants of VEGF were not associated with BC in Spanish population. The association between VEGF +936C/T, –1154A/G, –2578C/A, –634G/C, and –460T/C polymorphism and risk of BC was evaluated in a meta-analysis performed by Wang et al. 24 The overall results of combined analyses revealed that five polymorphisms of VEGF were not associated with the risk of BC. VEGF –2578 A allele and A carrier genotypes have been shown to be associated with an increased risk of renal cell carcinoma (RCC). 37
As the nature of breast carcinogenesis pathways is complex, there is no clear reason for the discrepancies in different studies. Ethnic, genetic, and/or environmental factors may interact in various ways to affect the risk of BC in different areas.
In conclusion, our findings revealed that VEGF rs699947 (–2578C>A) polymorphism was associated with the risk and clinicopathological characteristics of BC patients.
Author Contributions
Conceived and designed the experiments: MR, MH. Analyzed the data: MH, SS. Wrote the first draft of the manuscript: MR, MAM, MT. Contributed to the writing of the manuscript: MR, MH, SS, MAM, MT. Agree with manuscript results and conclusions: MR, MH, SS, MAM, MT. Jointly developed the structure and arguments for the paper: MR, SS, MAM, MT. Made critical revisions and approved final version: MH. All authors reviewed and approved of the final manuscript.
