Abstract
Introduction
Materials and Methods
Study population
A study group consisting of 101 women with histologically confirmed breast cancer, diagnosed at Child and Women Health Department, Gynecologic Oncology Division, Hospital Centre of Cova da Beira, Covilhã, Portugal, was enrolled between May 2008 and March 2009. A nonmalignant control group, 121 healthy female blood donors with no previous history of cancer and without family histories of cancers, was also studied. All subjects were Caucasian and had the same geographic area of residence. Both cases and controls completed a questionnaire that assessed parameters such as age, family history of breast cancer, parity, ethnicity, reproductive and menstrual history, and history of breast diseases and procedures. Informed consent was obtained from both patients and controls before entering in the study. The study was approved by the Institutional Review Board of Hospital Centre of Cova da Beira, Covilhã, Portugal.
DNA extraction
Genomic DNA of all cases and controls was isolated from either frozen or fresh blood samples using Wizard® Genomic DNA Purification Kit (Promega) according to the manufacturer's instruction, and stored at 4°C.
Genotyping
CYP19A1 codon 39 Trp/Arg (T/C) polymorphism, GSTM1 and GSTT1 homozygous deletions, and GSTP1 codon 105 Ile/Val (A/G) polymorphism analyses were performed by polymerase chain reaction (PCR)–based methods. CYP19A1 genotyping was performed using PCR with confronting two-pair primers, modified from the procedure described by Hirose et al. (2004). Briefly, each PCR mixture was carried out in a total volume of 25 μL that contained 10 pmol of each primer, 1.5 mM of MgCl2, 100 nM of each deoxynucleotide triphosphate, 1 unit of Taq DNA polymerase (Promega), and 100 ng of genomic DNA, using MyCycler Thermal Cycler (Bio-Rad). Reaction mixtures were preincubated for 10 min at 95°C. PCR conditions were 1 min at 95°C, 1 min at 56°C, and 1 min at 72°C, for 30 cycles. The final extension was at 72°C for 5 min. The amplified DNA was electrophoresed through 2% agarose gels stained with ethidium bromide. Genotypes were distinguished by the presence of a 200-bp band for T allele, a 264-bp band for C allele, and a 427-bp common band. GSTM1, GSTT1, and GSTP1 genotyping was performed as previously published by our group (Ramalhinho et al., 2011). Results were confirmed by re-genotyping 10% of randomly selected samples, and all the results were in agreement with the ones obtained previously.
Statistical analysis
To examine the association between genotypes and the development of breast cancer, we calculated odds ratios (ORs) and 95% confidence intervals (95% CIs) as estimates of relative risk, using logistic regression analysis with computer software SPSS for Windows (version 16.0). Chi-square test was applied to compare the allelic frequencies between normal controls and breast cancer patients. p-Values <0.05 were considered statistically significant.
Results
The mean age was 55.5 years (SD, 20.3 years; range, 17–84 years) for the controls and 62.8 years (SD, 12.9 years; range, 34–85 years) for the cases. Demographic data from cases and controls are shown in Table 1. The majority of study participants were postmenopausal (62.0% controls and 69.3% cases). The risk of breast cancer was higher for women with first degree family history of breast cancer (OR=2.532; 95% CI=2.133–3.005) and for women with history of benign breast disease (OR=1.959; 95% CI=1.146–3.350). Also, women with body mass index (BMI) superior to 25 seemed to have increased risk for breast cancer (OR=1.770; 95% CI=1.036–3.024). CYP19A1 codon 39 genotype results in cases and controls are presented in Table 2. Significant statistical association of the TC/CC genotypes combined with breast cancer risk, with reference to TT genotype, was documented (OR=1.770; 95% CI=1.036–3.024; p=0.036). Additionally, we analyzed the allele frequency in cases and controls (Table 3) and a positive correlation between C allele carriers (homozygous and heterozygous) and breast cancer risk was found (OR=1.540; 95% CI=1.013–2.341; p=0.043). Results regarding the genotype distribution of GSTM1, GSTT1, and GSTP1 in cases and controls in this population were previously published by our group, as well as the results of the association in two- and three-way combinations to evaluate the impact of gene–gene interaction (Ramalhinho et al., 2011). To investigate whether profiles of GSTs and CYP19A1 genotypes might be associated with the risk of breast cancer, we examined combinations of genotypes. The reference group consisted of individuals with all putative low-risk genotypes, that is, the presence of GSTM1 and GSTT1, the homozygous Ile/Ile genotype for GSTP1, and the homozygous TT genotype for CYP19A1. Heterozygous and homozygous individuals for the Val allele and also heterozygous and homozygous individuals for C allele were grouped for this analysis, due to the low number of individuals with homozygous Val allele and homozygous C allele, and so to increase the statistical power. We analyzed the two-way combination of GSTM1 and CYP19A1 genotypes (Table 4) and found that when GSTM1 is present, there is no significant association with breast cancer risk (OR=1.784; 95% CI=0.786–4.050; p=0.164) but we found a significant increase of breast cancer risk for carriers of GSTM1 null genotype along with CYP19A1 TT wild-type genotype (OR=3.115; 95% CI=1.362–7.124; p=0.006), and this risk is higher for carriers of both variants, GSTM1 null and CYP19A1 C allele, simultaneously (OR=6.158; 95% CI=2.676–14.171; p<0.001). We could find the same relation in the two-way analysis of GSTT1 and CYP19A1 genotypes (Table 5). There is no significant increase in breast cancer risk for carriers of GSTT1 present along with CYP19A1 C allele (OR=1.841; 95% CI=0.939–3.609; p=0.074), but for carriers of GSTT1 null genotype the risk is significantly higher in individuals with simultaneous CYP19A1 TT genotype (OR=4.599; 95% CI=1.792–11.804; p=0.001) or CYP19A1 TC/CC genotype (OR=4.870; 95% CI=2.216–10.700; p<0.001). Regarding the two-way analysis of CYP19A1 and GSTP1 genotypes (Table 6), we found no association of any genotype combination and risk of breast cancer. By analyzing the three-way combination of CYP19A1, GSTM1, and GSTT1 polymorphisms (Table 7), we found a positive association with breast cancer risk for women who carry C allele in CYP19A1 and null genotypes in GSTM1 (OR=5.000; 95% CI=1.491–16.771; p=0.007) or GSTT1 (OR=6.333; 95% CI=2.195–18.271; p=0.001). We also found that when both GSTs are deleted, the risk of breast cancer is higher, independently of CYP19A1 genotype, and the combination of all putative high-risk genotypes, CYP19A1 TC/CC, GSTM1 null, and GSTT1 null, carries a strong association with breast cancer risk (OR=11.429; 95% CI=3.590–36.385; p<0.001).
OR, odds ratio; CI, confidence interval; BMI, body mass index.
Discussion
High estradiol levels have been consistently associated with an increased breast cancer risk (Dorgan et al., 1996; Key et al., 2002; Missmer et al., 2004; Zeleniuch-Jacquotte et al., 2004), as estrogen promotes cellular growth and contributes to tumor growth by promoting the proliferation of cells with existing mutations or perhaps by increasing the opportunity for mutations (Zhu and Conney, 1998). Also, ROS have been related to the etiology of cancer, as they are known to be mitogenic to a variety of cells, and therefore capable of tumor promotion (Mitrunen and Hirvonen, 2003). Interindividual differences in estrogen biosynthetic and metabolic pathways may define subpopulations of women with higher lifetime exposure to hormone-dependent growth promotion or to cellular damage from estrogens or their metabolites. Such variation could explain a portion of the cancer susceptibility associated with reproductive events and hormone exposure. In most studies, polymorphisms in estrogen biosynthesis and metabolism genes are considered separately (Helzlsouer et al., 1998; Kristensen et al., 1998; Haiman et al., 2000; Healey et al., 2000; Millikan et al., 2000; Miyoshi et al., 2000; Mitrunen et al., 2001; Baxter et al., 2001; Gudmundsdottir et al., 2001; Hirose et al., 2004; Raskin et al., 2009; Ramalhinho et al., 2011) but here we hypothesize that SNPs in these low penetrance genes can exhibit synergistic effects on modulating individual susceptibility to breast cancer. Furthermore, specific associations of polymorphisms in estrogen biosynthesis and metabolism genes could result in a high-risk profile, by influencing lifetime levels of estrogen that could influence breast cancer risk. The four genes analyzed in this study are considered as low penetrating genes. These genes, alone or in association, may identify subjects as poor or fast synthesizers, like subjects as poor or fast estrogen metabolizers. While fast synthesizers are prolonging their exposure to estrogen, with potential increased carcinogenic activity, poor metabolizers are more exposed to the formation of carcinogen-DNA adducts and/or mutations, which confers them higher susceptibility to complex genetic disorders such as cancer. Thus, women having the “fast synthesizer/poor metabolizer” profile may possess higher risk for developing breast cancer than women with other profiles. In this study we investigated, in a population-based, case–control study of Caucasian women from a Portuguese region, breast cancer risk in association with polymorphism in a gene involved in estrogen biosynthesis (CYP19A1) along with polymorphisms in other three genes implicated in estrogen metabolism (GSTM1, GSTT1, and GSTP1). To avoid potential biases, we selected cases and controls from a single institution that only accepts patients from the same geographic area, Central Eastern Portugal. Patients and controls were matched by ages and belonged to the same ethnic group (Caucasian). Regarding CYP19A1 codon 39 genotype distribution, our results were in agreement with the few previously published reports. Miyoshi et al. (2000) identified this polymorphism and performed a case–control study in order to evaluate its association with breast cancer risk in a Japanese population. They concluded that homozygous and heterozygous carriers of the variant Arg (C) allele at codon 39 had a decreased risk of breast cancer compared with the noncarriers. In our study, breast cancer patients had higher relative frequency of C allele and lower prevalence of Trp (T) allele than controls. Thus, for developing breast cancer, C allele carriers apparently have higher risk, while T allele carriers have lower risk, consequently T allele seems to be protective for carriers. A similar study, also in a Japanese population, was conducted by Hirose et al. (2004). They found, in homozygous and heterozygous carriers of the C allele, a significantly increased risk of breast cancer among premenopausal women with a late age at first full term pregnancy or a high BMI. In our study, we found positive statistical association of the TC/CC genotypes (C allele carriers) with breast cancer risk with no association with any particular demographic data. More recently, Tuzuneri et al. (2010) also found an association between C allele-bearing and breast cancer risk, and they propose TT haplotype as a protective haplotype. No published reports, regarding functional effects of CYP19 codon 39 Trp/Arg amino acid change, are of our knowledge. To explain our data, we can speculate about deleterious effect of Trp to Arg substitution in codon 39 of CYP19A1 gene; it seems that the encoded protein P450 aromatase variant tends to have higher activity, thereby raising estrogens levels, what leads to the association of C allele with breast cancer risk. Considering GSTs polymorphisms, our group has already demonstrated the significant increase of breast cancer risk associated with GSTM1 and GSTT1 null polymorphisms, both alone or in combination, in the same population, while GSTP1 genotypes seem to have no influence in breast cancer susceptibility (Ramalhinho et al., 2011). In that article we reported that GSTM1 null genotype was significantly more common among breast cancer cases compared with controls, so as GSTT1, and no significant increase of breast cancer risk associated with GSTP1 genotypes was found. We also observed an eight-fold increased breast cancer risk in women who carry null genotype both in GSTM1 and GSTT1, and we found that Ile/Ile genotype, so as presence of Val allele, seemed to be associated with risk of breast cancer when combined with both GSTM1 and GSTT1 null genotypes; thus, it appears that the increase of breast cancer risk was mainly linked to both GSTM1 and GSTT1 deletion, rather than to any genotype of GSTP1. Grouping CYP19A1 and GSTM1 genotypes, breast cancer susceptibility was not altered when GSTM1 was present. However, breast cancer risk for carriers of GSTM1 null genotype along with CYP19A1 TT wild-type genotype is increased, so as for simultaneous carriers of both variants, GSTM1 null and CYP19A1 C allele. Similar association was observed in the analysis of CYP19A1 and GSTT1 genotypes; independently of CYP19A1 polymorphism, breast cancer risk is lower for carriers of GSTT1. Again, no association between both CYP19A1 and GSTP1 genotypes taken together and risk of breast cancer was found. To our knowledge, this is the first study to analyze the potential role of CYP19A1 codon 39 polymorphism in Portuguese women, and its possible combined effect with GSTM1, GSTT1, and GSTP1 polymorphisms in the development of breast cancer. We determined that C allele, in homozygosity or heterozygosity, is significantly associated with an increased risk of breast cancer in the population studied. Furthermore, it seems that the effects of CYP19A1 T/C polymorphism in estradiol biosynthesis appear to be modulated by the presence of GSTM1 and GSTT1 in estrogen metabolic pathway, because breast cancer susceptibility is lower in carriers of GSTM1 and GSTT1, independently of CYP19A1 genotype. The results support the accepted hypothesis that estrogen formation and metabolism, which can be dependent of genetic and environmental variation, are determinant in defining breast cancer risk. Further studies are needed to evaluate the influence of the expression of these genes, especially in association with the expression of other known low penetrance genes in defining breast cancer risk.
Footnotes
Disclosure Statement
The authors declare that we have no conflict of interest.
Acknowledgments
The authors would like to thank all the technical staff from Centro Hospitalar Cova da Beira, Covilhã, Portugal, for their kindly cooperation in the collection of the blood samples and all the volunteers who agreed to participate in this study.
