Skip to main content
  • Research article
  • Open access
  • Published:

Oxidative stress in susceptibility to breast cancer: study in Spanish population



Alterations in the redox balance are involved in the origin, promotion and progression of cancer. Inter-individual differences in the oxidative stress regulation can explain a part of the variability in cancer susceptibility.

The aim of this study was to evaluate if polymorphisms in genes codifying for the different systems involved in oxidative stress levels can have a role in susceptibility to breast cancer.


We have analyzed 76 single base polymorphisms located in 27 genes involved in oxidative stress regulation by SNPlex technology. First, we have tested all the selected SNPs in 493 breast cancer patients and 683 controls and we have replicated the significant results in a second independent set of samples (430 patients and 803 controls). Gene-gene interactions were performed by the multifactor dimensionality reduction approach.


Six polymorphisms rs1052133 (OGG1), rs406113 and rs974334 (GPX6), rs2284659 (SOD3), rs4135225 (TXN) and rs207454 (XDH) were significant in the global analysis. The gene-gene interactions demonstrated a significant four-variant interaction among rs406113 (GPX6), rs974334 (GPX6), rs105213 (OGG1) and rs2284659 (SOD3) (p-value = 0.0008) with high-risk genotype combination showing increased risk for breast cancer (OR = 1.75 [95% CI; 1.26-2.44]).


The results of this study indicate that different genotypes in genes of the oxidant/antioxidant pathway could affect the susceptibility to breast cancer. Furthermore, our study highlighted the importance of the analysis of the epistatic interactions to define with more accuracy the influence of genetic variants in susceptibility to breast cancer.

Peer Review reports


Despite breast cancer (BC) being the most frequent cancer in women in western countries and the second cause of cancer death after lung cancer [1], the risk factors that lead to the disease are not completely understood, although is widely accepted that they include a combination of environmental and genetic factors. For genetic approximation, a polygenic model has been proposed in which a combination of common variants, having individually a modest effect, together contribute to BC predisposition [2].

Numerous evidence links carcinogenesis and oxidative stress regulation, including prooxidant and antioxidant defense systems [37]. Oxidative stress is defined as an imbalance in the production of reactive oxygen species (ROS) and reactive nitrogen species (RNS) and their removal by antioxidants. When this imbalance occurs, biomolecules are damaged by ROS and RNS and normal cellular metabolism is impaired, leading to changes of intra- and extracellular environmental conditions. ROS can cause lesions in DNA, such as mutations, deletions, gene amplification and rearrangements, that may lead to malignant transformations and cancer initiation and progression [810]. The effect of ROS and RNS, however, is balanced by the anti-oxidant action of non-enzymatic and anti-oxidant enzymes maintaining cellular redox levels under physiological conditions [4, 11].

Previous studies with knockout animals that lack antioxidant enzymes support the view that ROS contribute to the age-related development of cancer. For instance, mice deficient in the antioxidant enzyme CuZnSOD showed increased cell proliferation in the presence of persistent oxidative damage contributing to hepatocarcinogenesis later in life [12]. Another study showed that mice lacking the antioxidant enzyme Prdx1 had a shortened lifespan owing to the development, beginning at about 9 months, of severe hemolytic anemia and several malignant cancers [13].

In this context, single nucleotide polymorphisms (SNPs) in components of the cellular redox systems can modify the redox balance and take part in both the BC initiation and/or progression, as well as determine possible therapeutic treatments [1417].

Despite the importance of oxidative stress in the development and progression of cancer, few studies have evaluated the relationship between genetic modification in genes coding for enzymes relatives to the redox system and the susceptibility to develop BC. The previous studies had focused mainly on the analysis of genes related to antioxidant defense enzymes [18, 19], but the information about modifications in genes involved in the oxidation process is relatively sparse.

The aim of this study was to evaluate the association between common variants in genes coding for proteins related to the redox system (antioxidant and oxidant systems or proteins) and the susceptibility to develop BC. We hypothesized that common SNPs related to the redox pathway are associated with an altered risk for BC. We chose 76 SNPs on which to perform a two-step study: one first exploratory set and a second, independent, validation set. We also decided to investigate the impact of complex interactions between SNPs at different genes of the stress oxidative pathway. To address this issue, we analyzed the effects of gene-gene interactions by the multifactor dimensionality reduction (MDR) approach. This analysis was carried out in four SNPs that were statistically significant in the combinatorial set.


Study population

The underlying analyses were carried out in a Caucasian Spanish population. The study was carried out in two steps with two population groups. A first group of 1176 samples was composed of 493 female patients diagnosed for BC between the years 1998–2008 at La Paz Hospital and Foundation Jimenez Díaz (Madrid), and 683 healthy women controls recruited at the Hospital of Valladolid (Spain).

Thereupon, we chose the polymorphisms that showed marginally significant association (p-value < = 0.15), and we replicated the procedure in a second independent group (n = 1233) where we included 430 female patients diagnosed for BC between the years 1988–1998 at the Clinic Hospital of Valencia (Spain) and 803 samples from cancer-free women recruited at the blood donor bank at the same Hospital. Blood was collected between 2010 and 2011 during periodical patient visits. The blood from controls was extracted between the years 2009 and 2012. In both groups, the controls were women without pathology or history of cancer. Controls were not matched to cases, but were similar in age. In group 1, cases’ mean age was 57.5 (range 23.5-89.5), and that for donors was 52.7 (21.5-96.5). In group 2, cases’ mean age was 54.1 (20.5-86.5) while in donors, it was 54 (22.5-92.5).

We selected this staged approach because it allowed us to analyze only those polymorphisms with indicative results and reduced the number of genotyping reactions without significantly affecting statistical power [18, 20].

The research protocols were approved by the ethics committee of the INCLIVA Biomedical Research Institute. All the participants in the study were informed and gave their written consent to participate in the study.

Single nucleotide polymorphisms selection and genotyping

Two public databases were used to collect information about SNPs in oxidative pathway genes: NCBI ( and HapMap ( The selection of polymorphisms was performed by SYSNP [20] and by a literature search in PubMed, Scopus and EBSCO databases using the terms “breast cancer and polymorphisms and oxidative”, along with additional terms such as “SNPs and oxidative pathway and susceptibility”, and their possible combinations. The following criteria were used to select the SNPs: functional known or potentially functional effect, location in promoter regions, minor allele frequency (MAF) over 0.1 in Caucasian populations analyzed previously, localization and distribution along the gene (including upstream and downstream regions) and low described linkage disequilibrium between candidate polymorphisms. We included variants with potential influence in the gene and protein function, as well as the most important variants described in the literature.

Finally, we select a total of 76 polymorphisms located in 27 genes related to the redox system: 17 were classified as antioxidant genes (CAT, GCLC, GCLM, GNAS, GPX6, GSR, GSS, M6PR, MSRB2, OGG1, SOD1, SOD2, SOD3, TXN, TXN2, TXNRD1, TXNRD2) and 10 as reactive species generators (mainly NADPH oxidase-related genes CYBB, NCF2, NCF4, NOS1, NOS2A, NOX1, NOX3, NOX4, NOX5 and XDH). Reference names and characteristics of the selected SNPs are provided in Table 1.

Table 1 Summary of the 76 selected SNPs in 27 genes

Experimental procedures

The blood samples remained frozen until the DNA extraction was performed. Genomic DNA was extracted from blood samples using DNA Isolation Kit (Qiagen, Izasa, Madrid, Spain) following the manufacturer’s protocol, but a final elution volume of 100 μl used. DNA concentration and quality were measured in a NanoDrop spectrophotometer. Each DNA sample was stored at -20°C until analysis, which in all cases was performed within a year of the DNA extraction.

Genotyping analysis in both sets was performed by SNPlex technology (Applied Biosystems, Foster City, California, USA) according to the manufacturer’s protocol [21]. This genotyping system, based on oligation assay/polymerase chain reaction and capillary electrophoresis, was developed for accurate genotyping, high sample throughput, design flexibility and cost efficiency. It has validated its precision and concordance with genotypes analyzed using TaqMan probes-based assays. The sets of SNPlex probes were reanalyzed in about 10% of the samples with a reproducibility of over 99%. Those polymorphisms and samples with genotyping lower than 85% in the first set were excluded from further analysis.

Statistical and MDR analyses

Statistical analysis was performed using SNPstats software [22], a free web-based tool, which allows the analysis of association between genetic polymorphisms and diseases. The proper analysis of these studies can be performed with general purpose statistical packages, but this software facilitates the integration of data. The association with disease is modeled as binary; the application assumes an unmatched case–control design and unconditional logistic regression models are used. The statistical analyses are performed in a batch call to the R package ( SNPStats returns a complete set of results for the analysis. SNPstats provides genotype frequencies, proportions, odds ratios (OR) and 95% confidence intervals (CI), and p-values for multiple inheritance models. The lowest Akaike’s Information Criterion and Bayesian Information Criterion values indicate the best inheritance genetic model for each specific polymorphism. All the analyses were adjusted by age. Only SNPs with no significant deviation from Hardy-Weinberg equilibrium (HWE) in controls and a MAF exceeding 5% were retained for the association analysis (Table 1).

To identify gene-gene interactions, MDR was used. It is a non-parametric and a genetic model-free approach that uses a data reduction strategy [2325]. This method considers a single variable that incorporates information from several loci that can be divided into high risk and low risk combinations. This new variable can be evaluated for its ability to classify and predict outcome risk status using cross validation and permutation testing. Both were used to prevent over-fitting and false-positives from the multiple testing. With n-fold cross-validation, the data are divided into n equal size pieces. An MDR model is fit using (n-1)/n of the data (the training set) and then evaluated for its generalizability on the remaining 1/n of the data (the testing set). The fitness of a MDR model is assessed by estimating accuracy in the training set and the testing set. Moreover, it estimates the degree to which the same best model is discovered across n divisions of the data, referred to as the cross-validation consistency (CVC). The best MDR model is the one with the maximum testing accuracy. Statistical significance is determined using permutation testing. We used 10-fold cross-validation and 1000-fold permutation testing. MDR results were considered statistically significant at the 0.05 level. The advantages of this method are that there are no underlying assumptions about the independence or biological relevance of SNPs or any other factor. This is important for diseases as sporadic BC where the etiology is not completely known. We used the MDR software (version 2.0 beta 8.4) which is freely available (


Single nucleotide polymorphisms and susceptibility to breast cancer

Set 1: To determine the possible association of polymorphisms related to oxidative stress genes and BC we analyzed 76 polymorphisms in 27 genes of the redox system in 493 cases and 683 controls (Table 1). Seven SNPs (rs3749930, rs2036343, rs34990910, rs17881274, rs17011353, rs17011368, rs17323225) with MAF <0.05 in controls, as along with two SNPs (rs725521 and rs231954) not showing Hardy-Weinberg equilibrium, were excluded from the association analysis (Table 1). A total of 67 SNPs were successfully genotyped and analyzed.

Our association analysis in set 1 pointed out four nominally statistically significant results (p < 0.05). Table 2 shows the results found in the selected polymorphisms. Polymorphisms rs974334, rs1805754 rs4135225 and rs207454 showed an association with modifications in the risk for BC. All the results were adjusted by age.

Table 2 Comparison of genotype frequencies between breast cancer patients and controls (Set 1)

Set 2: Subsequently, and in order to better identify those polymorphisms that could be associated with BC, we replicated the 10 SNPs with a p-value equal to or lower than 0.15 in group 1 [rs3736729, OR: 0.74 (0.54-1.01); rs406113, OR:1.26 (0.98-1.62); rs974334, OR:2.01 (1.07-3.80); rs1805754, OR: 1.31 (1.02-1.68); rs1052133, OR: 1.76 (1.00-3.10); rs2284659, OR: 1.30 (0.92-1.84); rs2301241, OR: 0.80 (0.60-1.07); rs4135179, OR: 1.27 (0.97-1.66); rs4135225, OR: 0.66 (0.45-0.96); rs207454, OR: 4.98 (1.28-19.34)] in a second independent set.

Set 1 + Set 2: Finally, we analyzed the 10 polymorphisms in the global population set 1 + set 2 (n = 2409; cases = 923, controls = 1486). The results are listed in Table 3. From the 10 polymorphisms analyzed in both samples, 6 presented a statistically significant association with increased risk when the combined data were analyzed: rs406113 [OR: 1.23 (1.04-1.46)], rs974334 [OR: 1.73 (1.09-2.73)], rs1052133 [OR:1.82 (1.31-2.52)], rs2284659 [OR:1.33 (1.05-1.67), rs4135225 [OR: 0.77 (0.60-0.99)], rs207454 [OR: 2.12 (1.11-4.04)]. Of these polymorphisms, the rs105213 on the OGG1 gene maintained the statistical significance (p-value = 0.0004) after the Bonferroni correction.

Table 3 Genotype frequencies of relevant polymorphisms in different Sets

Gene-gene interactions in breast cancer patients

There is growing evidence that epistasis interactions between genes may play a role in cancer risk, and different variable selection approaches have been developed to analyze the potential gene-gene and gene-environment interactions [25]. The four most significantly associated polymorphisms in set 1 + set 2 with susceptibility to BC were selected for this analysis: rs406113 [OR: 1.23 (1.04-1.46)], rs974334 [OR: 1.73 (1.09-2.73)], rs1052133 [OR:1.82 (1.31-2.52)] and rs2284659 [OR:1.33 (1.05-1.67)]. Data from 1182 samples (controls and patients) from both groups were used. The combination was performed grouping the genotypes according to the model predicted for the four polymorphisms: recessive model for rs1052133 (CC and CG were grouped into a single block), dominant model for rs406113 (CC and AC genotypes were grouped into a single block), recessive model for rs974334 (CC and CG genotypes were grouped into a single block) and recessive model for rs2284659 (GG and GT genotypes were grouped into a single block). For a two-loci interaction, the combination of polymorphisms rs406113 (GPX6) and rs1052133 (OGG1) was the most significant (p = 0.041). The best three-loci model included rs406113 on the GPX6 gene, rs1052133 on the OGG1 gene and rs2284659 on the SOD3 gene, and it showed statistical significance (p < 0.0007) with an OR = 1.82 and 95% CI = 1.28-2.58. A four-way interaction found that between rs406113 on the GPX6 gene, rs974334 on the GPX6 gene, rs1052133 on the OGG1 gene and rs2284659 on the SOD3 gene predicts breast cancer with a testing balance accuracy of 0.5267. This four-loci model had a chi-square value of 11.284 (p = 0.0008) and an OR of 1.75 [95% CI = 1.26-2.44]. The four polymorphism combinatory model showed a higher predisposition to BC than the polymorphisms rs406113, rs974334 and rs2284659 did individually (ORX2 = 1.23, ORX3 = 1.73, ORX6 = 1.33) and had values similar to the ones of polymorphism rs1052133 (ORX5 = 1.82). The summary of the multi-factor dimensionality results are listed in Table 4.

Table 4 Summary of Multi-factor Dimensionality (MDR) results

The combined genotype AA for rs406113, CC/CG for rs974334, CC/CG for rs1052133 and GG/ GT for rs2284659 showed a higher risk for BC, which is consistent with the models described for the polymorphisms individually. Figure 1 summarizes the four-loci genotype combinations associated with high and low risk and with the distribution of cases and controls.

Figure 1
figure 1

The polymorphisms rs406113, rs974334, rs1052133 and rs2284659, showing the highest statistical significance in the combinatorial set 1+set 2, were chosen for the gene-gene interaction analysis. The MDR analysis was done with the genotypes collapsed according to the genetic models selected: rs1052133 (OGG1) recessive model; CC/CG vs.GG, rs406113 (GPX6) dominant model; CC/AC vs. AA, rs974334 (GPX6) recessive model; CC/CG vs. GG, rs2284659 (SOD3) recessive model; GG/GT vs. TT. The figure shows the summary of four-loci genotype combinations associated with high and low risk. Cases: left bars, controls: right bars. The epistatic gene-gene interaction corresponds to the high risk combinations (darkest color).


Genetic association studies involving SNPs and their possible interactions have become increasingly important for the study of human diseases. The present study has focused on genes encoding for proteins of the redox system. It is long proven that they are clearly involved in extensive damage to DNA, which in turn leads to gene mutations and, finally, carcinogenesis. The functionality of polymorphisms in relation to oxidative stress has been proven in several cases. For instance, the polymorphism in exon 2 of the superoxide dismutase 2 (SOD2) gene A16V (C/T) (rs4880) led to structural alterations in the domain responsible to target the mitochondria, giving a reduction in the antioxidant potential [26]. Furthermore, a functional polymorphism in exon 9 of the CAT gene and other polymorphisms in endothelia NO synthase (eNOs) that seem relevant for their activity have been documented [2628]. Therefore, it is clear that a single oligonucleotide modification can lead to structural changes, modifications in the affinity to bind proteins or in activity that may be relevant to the redox system. Our hypothesis was that variations in genes from the stress oxidative pathway, that have shown to have a possible linkage to BC, can be associated with predisposition to this disease. Indeed, genetic variations in these pathways have shown to modify the risk for BC [29].

In the present epidemiological study, we have assessed the effect of 76 SNPs in 27 genes in a case–control study in a Spanish population. Genotype distributions in the controls did not differ significantly from those expected under HWE. The study was performed in two independent sets of patients and controls, first, to select the relevant polymorphisms and second, to check the reproducibility and significance of these preliminary results.

Six SNPs (rs406113 and rs974334 on the glutathione peroxidase 6 (GPX6) gene, rs1052133 on the 8-oxoguanine DNA glycosylase (OGG1) gene, rs2284659 on the superoxide dismutase 3 (SOD3) gene, rs4135225 on the thioredoxin (TXN) gene and rs207454 on the xanthine dehydrogenase (XDH) gene) are associated with variations in the predisposition to BC.

The rs406113 (c.39 T > G; p.F13L) and rs974334 (c.242-12G > C) polymorphisms on the GPX6 gene had not been studied previously; in fact, there was no information available in the literature about polymorphisms on the GPX6 gene even though they can have a functional effect. Genetic variants in other genes of the GPX family have been associated with BC [3032].

Thioredoxin (TXN) is overexpressed in BC, and it is related to tumor grade [33], being a crucial element in redox homeostasis [34]. Studies of polymorphisms in the TXN gene, encoding thioredoxin, are few in cancer. Seibold et al. [19], evaluated the influence of common variants on TXN, thioredoxin reductase 1 (TXNRD1) and thioredoxin 2 (TXN2) genes and the risk of BC after menopause, including seven of the SNPs analyzed in our study. Rs2301241 and rs2281082 were not significantly related to BC risk in our study, however, Seibold et al. found a limited association of rs2301241 with BC risk when comparing rare homozygote vs. common homozygote. Other studies found a borderline significance [18]. In the case of rs2281082, the borderline association of the Seibold study was not confirm in other publications [19]. In our study population, we found that carriers of one T allele on rs4135225 (c.196-192C > T) were associated with lower risk for BC development (OR = 0.77 [95% CI; 0.60-0.99] p = 0.041). Seidbol et al. found a predisposition to BC for this polymorphisms (OR = 1.22 [95% CI; 1.06-1.41]. One must take into account that the Seidbold study is focused in postmenopausal women, unlike ours. Still, in their analysis, they compared only two of the three possible genotypes (heterozygotes vs. common homozygotes). In any case, their results showed borderline significance.

The xanthine dehydrogenase (XDH) is an important enzyme involved in the first-pass metabolism of 6-mercaptopurine [35]. Polymorphisms in the XDH gene have been related to cancer. The rs1884725 polymorphism has been identify as a genetic variant associated with disease risk and outcomes in multiple myeloma [36]. In our study, one of the thirteen polymorphisms evaluated on this gene showed an association with BC risk. Carriers of one A allele of rs207454 displayed 2.12 times ([95% CI; 1.11-4.04], p = 0.024) more risk to develop the illness than did non-carriers. To our knowledge, there are no studies of these polymorphisms in the literature. The results presented here suggest an association with the development of BC, although further confirmatory studies would be needed to confirm it.

The polymorphism on the SOD3 gene (rs2284659) analyzed in our study showed a trend to the predisposition for BC in the global analysis. There was no information about this polymorphism in the literature. Other polymorphisms in this gene, like rs2536512 and rs699473, have been associated in BC patients with the incidence of tumor and poorer progression-free survival (PFS) [37]. Moreover, some results suggest that rs699473 may influence brain tumor risk [38].

The variant rs1052133 (Ser326Cys) in the OGG1 gene has the same tendency to predisposition for BC in both sets, separately and in the combined data set. Concerning this polymorphism, previous studies had conflicting results [3945]. Three meta-analyses have attempted to summarize the results [39, 41, 46]. In one study, the authors analyzed this polymorphism in relation to several cancers founding only significant association with the risk for lung cancer [46]. The others two meta-analyses are focused on BC, and the results are contradictory. Yuan et al. found an association just in the European population subgroup [41], while Gu et al. did not show any association, even when stratifying the analysis by ethnicity or menopausal status [39]. These differences may have arisen from the different number of studies included in the European group.

We found an increment for the risk to develop BC in the carriers of at least one Ser allele (recessive model) if we consider the sets both separately and together ((OR = 1.82 [95% CI 1.31-2.52]) and p-value = 0.0004). Our results are in concordance with the meta-analysis by Yuan and collaborators that suggests that the hOGG1 326 Cys allele provides a significant protective effect for BC in European women [41]. The importance of this SNP rests in the role of the 8-oxoguanine DNA glycosylase, encoded by OGG1[47]. This enzyme can excise the 8-hydroxy-2´-deoxyguanosine (8-OHdG) modifications occurring in the DNA as a result of hydroxyl radical interaction [41, 48, 49]. An incorrect expression of the protein could interfere with the suitable repair of the genetic material. Other polymorphisms in the OGG1 gene, like rs2304277, and recently described by Osorio et al., have been associated with ovarian cancer risk in BRCA1 mutation carriers [50]. This data certainly support the importance of genetic changes in the OGG1 gene in relation to the predisposition to cancer.

The epistatic analysis of the four most significant polymorphisms in relation to the susceptibility to BC was performed with the MDR method. This is a reliable approach that has been widely used [23, 25, 5154]. The combination was performed grouping the genotypes according to the model predicted for the four polymorphisms in Tables 2 and 3. The result obtained was an OR = 1.75 [95% CI = 1.26-2.44; p-value = 0.0008], a value similar to that obtained for rs1052133 (OR = 1.82 [95% CI = 1.31-2.52; p-value = 0.0004]). The previous study of Cebrian et al. in antioxidant defence enzymes and BC susceptibility has twelve common SNPs with our study. Two showed discrepancies with our data: Rs511895 in the CAT gene was not significant in our analysis, but it presented a borderline tendency in the Cebrian et al. study. Moreover, they found a significant difference in genotype distribution between cases and controls in rs4135179 (TXN). We, however, were unable to confirm this in our global analysis, although we detected a marginal significance in set 1. The reason for this discrepancy can be found in the population’s characteristics, in the superior age of the population included in the Cebrian study [18].

Our study has several limitations to take into consideration. Firstly, there is no data available about the lifestyle of controls and patients that could be related to oxidative stress, such as diet, exercise and the consumption of tobacco and alcohol. Secondly, polymorphisms that were not explored in our study may affect the risk to develop BC and should be taken into account in the analysis of our data and in further studies. Nevertheless, the association between SNPs and risk for BC is reliable since that power exceeded 95% in all the cases. All samples are from the same country and ethnicity, and the adjustment for age reduces variability.

Additionally, MDR has 80% statistical power to detect true interactions in two-, three-, and four-way gene-gene interactions, even with a small number of cases and controls [24]. Furthermore, several associations detected in these data involved SNPs occurring in non-coding regions. However, variations in the intronic structure have been proposed to influence cancer susceptibility via regulation of gene expression, gene splicing or mRNA stability. It is also possible that these polymorphisms are in linkage disequilibrium with other functional polymorphisms that may affect BC susceptibility.

Despite these considerations, our work, as far as we know, is the largest study in the Spanish population that analyzes the influence of polymorphisms in oxidative genes in susceptibility to BC. Overall, our data, together with that published in the bibliography [18, 19, 29, 37, 41, 45, 5562], suggest a role of stress-response gene variants in the susceptibility to BC.


Our results suggest that different genotypes in genes of the oxidant/antioxidant pathway could affect the susceptibility to breast cancer. We have found six polymorphisms in OGG1, GPX6, SOD3, TXN and XDH genes significantly associated with predisposition to breast cancer. These associations have not been described previously, except for rs1052133 (OGG1). Concerning this polymorphism the published results in breast cancer were contradictory, and some authors found only a significant association with the risk of developing lung cancer. We have found an increment in the risk of developing breast cancer in the carriers of at least one Ser allele (recessive model) in concordance with a meta-analysis of breast cancer susceptibility in European women. In this particular case, an incorrect expression of the protein encoded by the OGG1 gene could interfere with the suitable repair of the genetic material. Furthermore, our study highlighted the importance of the analysis of the epistatic interactions in order to define the influence of genetic variants in susceptibility to breast cancer more precisely. Further studies on the relevance of these and other polymorphisms in the development of breast cancer should be performed.



Breast cancer


Breast cancer 1


Confidence interval




Cytochrome b-245, beta polypeptide


Glutamate-cysteine ligase, catalytic subunit


Glutamate-cysteine ligase, modifier subunit


GNAS complex locus


Glutathione peroxidase 6


Glutathione reductase


Glutathione synthetase


Hardy-Weinberg Equilibrium


Minor Allele Frequency


Multifactor dimensionality reduction


mannose-6-phosphate receptor


Methionine sulfoxide reductase B2


Neutrophil cytosolic factor 2


Neutrophil cytosolic factor 4


Nitric oxide synthase 1


Nitric oxide synthase 2


NADPH oxidase 1


NADPH oxidase 3


NADPH oxidase 4


NADPH oxidase 5


8-oxoguanine DNA glycosylase


Odds ratio


Reactive Nitrogen Species


Reactive Oxygen Species


Superoxide dismutase 1


Superoxide dismutase 2


Superoxide dismutase 3


Single Nucleotide Polymorphisms




Thioredoxin 2


Thioredoxin reductase 1


Thioredoxin reductase 2


Xanthine dehydrogenase.


  1. Dumitrescu RG, Cotarla I: Understanding breast cancer risk – where do we stand in 2005?. J Cell Mol Med. 2005, 9 (1): 208-221. 10.1111/j.1582-4934.2005.tb00350.x.

    Article  CAS  PubMed  Google Scholar 

  2. Stratton MR, Rahman N: The emerging landscape of breast cancer susceptibility. Nat Genet. 2008, 40 (1): 17-22. 10.1038/ng.2007.53.

    Article  CAS  PubMed  Google Scholar 

  3. Sun Y: Free radicals, antioxidant enzymes, and carcinogenesis. Free Radic Biol Med. 1990, 8 (6): 583-599. 10.1016/0891-5849(90)90156-D.

    Article  CAS  PubMed  Google Scholar 

  4. Valko M, Leibfritz D, Moncol J, Cronin MT, Mazur M, Telser J: Free radicals and antioxidants in normal physiological functions and human disease. Int J Biochem Cell Biol. 2007, 39 (1): 44-84. 10.1016/j.biocel.2006.07.001.

    Article  CAS  PubMed  Google Scholar 

  5. Valko M, Rhodes CJ, Moncol J, Izakovic M, Mazur M: Free radicals, metals and antioxidants in oxidative stress-induced cancer. Chem Biol Interact. 2006, 160 (1): 1-40. 10.1016/j.cbi.2005.12.009.

    Article  CAS  PubMed  Google Scholar 

  6. Parri M, Chiarugi P: Redox molecular machines involved in tumor progression. Antioxid Redox Signal. 2013, 19 (15): 1828-1845. 10.1089/ars.2012.5040.

    Article  CAS  PubMed  Google Scholar 

  7. Scott TL, Rangaswamy S, Wicker CA, Izumi T: Repair of oxidative DNA damage and cancer: recent progress in DNA base excision repair. Antioxid Redox Signal. 2014, 20 (4): 708-726. 10.1089/ars.2013.5529.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Cerutti PA: Prooxidant states and tumor promotion. Science. 1985, 227 (4685): 375-381. 10.1126/science.2981433.

    Article  CAS  PubMed  Google Scholar 

  9. Trush MA, Kensler TW: An overview of the relationship between oxidative stress and chemical carcinogenesis. Free Radic Biol Med. 1991, 10 (3–4): 201-209.

    Article  CAS  PubMed  Google Scholar 

  10. Costa A, Scholer-Dahirel A, Mechta-Grigoriou F: The role of reactive oxygen species and metabolism on cancer cells and their microenvironment. Semin Cancer Biol. 2014, 25: 23-32.

    Article  CAS  PubMed  Google Scholar 

  11. Trachootham D, Lu W, Ogasawara MA, Nilsa RD, Huang P: Redox regulation of cell survival. Antioxid Redox Signal. 2008, 10 (8): 1343-1374. 10.1089/ars.2007.1957.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Elchuri S, Oberley TD, Qi W, Eisenstein RS, Jackson Roberts L, Van Remmen H, Epstein CJ, Huang TT: CuZnSOD deficiency leads to persistent and widespread oxidative damage and hepatocarcinogenesis later in life. Oncogene. 2005, 24 (3): 367-380. 10.1038/sj.onc.1208207.

    Article  CAS  PubMed  Google Scholar 

  13. Neumann CA, Krause DS, Carman CV, Das S, Dubey DP, Abraham JL, Bronson RT, Fujiwara Y, Orkin SH, Van Etten RA: Essential role for the peroxiredoxin Prdx1 in erythrocyte antioxidant defence and tumour suppression. Nature. 2003, 424 (6948): 561-565. 10.1038/nature01819.

    Article  CAS  PubMed  Google Scholar 

  14. Acharya A, Das I, Chandhok D, Saha T: Redox regulation in cancer: a double-edged sword with therapeutic potential. Oxid Med Cell Longev. 2010, 3 (1): 23-34. 10.4161/oxim.3.1.10095.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Fang J, Seki T, Maeda H: Therapeutic strategies by modulating oxygen stress in cancer and inflammation. Adv Drug Deliv Rev. 2009, 61 (4): 290-302. 10.1016/j.addr.2009.02.005.

    Article  CAS  PubMed  Google Scholar 

  16. Ivanova D, Bakalova R, Lazarova D, Gadjeva V, Zhelev Z: The impact of reactive oxygen species on anticancer therapeutic strategies. Adv Clin Exp Med. 2013, 22 (6): 899-908.

    PubMed  Google Scholar 

  17. Milaeva ER: Metal-based antioxidants–potential therapeutic candidates for prevention the oxidative stress-related carcinogenesis: mini-review. Curr Top Med Chem. 2011, 11 (21): 2703-2713. 10.2174/156802611798040741.

    Article  CAS  PubMed  Google Scholar 

  18. Cebrian A, Pharoah PD, Ahmed S, Smith PL, Luccarini C, Luben R, Redman K, Munday H, Easton DF, Dunning AM, Ponder BA: Tagging single-nucleotide polymorphisms in antioxidant defense enzymes and susceptibility to breast cancer. Cancer Res. 2006, 66 (2): 1225-1233. 10.1158/0008-5472.CAN-05-1857.

    Article  CAS  PubMed  Google Scholar 

  19. Seibold P, Hein R, Schmezer P, Hall P, Liu J, Dahmen N, Flesch-Janys D, Popanda O, Chang-Claude J: Polymorphisms in oxidative stress-related genes and postmenopausal breast cancer risk. Int J Cancer. 2011, 129 (6): 1467-1476. 10.1002/ijc.25761.

    Article  CAS  PubMed  Google Scholar 

  20. Hirschhorn JN, Daly MJ: Genome-wide association studies for common diseases and complex traits. Nat Rev Genet. 2005, 6 (2): 95-108.

    Article  CAS  PubMed  Google Scholar 

  21. Tobler AR, Short S, Andersen MR, Paner TM, Briggs JC, Lambert SM, Wu PP, Wang Y, Spoonde AY, Koehler RT, Peyret N, Chen C, Broomer AJ, Ridzon DA, Zhou H, Hoo BS, Hayashibara KC, Leong LN, Ma CN, Rosenblum BB, Day JP, Ziegle JS, De La Vega FM, Rhodes MD, Hennessy KM, Wenz HM: The SNPlex genotyping system: a flexible and scalable platform for SNP genotyping. J Biomol Tech. 2005, 16 (4): 398-406.

    PubMed  PubMed Central  Google Scholar 

  22. Sole X, Guino E, Valls J, Iniesta R, Moreno V: SNPStats: a web tool for the analysis of association studies. Bioinformatics. 2006, 22 (15): 1928-1929. 10.1093/bioinformatics/btl268.

    Article  CAS  PubMed  Google Scholar 

  23. Moore JH, Gilbert JC, Tsai CT, Chiang FT, Holden T, Barney N, White BC: A flexible computational framework for detecting, characterizing, and interpreting statistical patterns of epistasis in genetic studies of human disease susceptibility. J Theor Biol. 2006, 241 (2): 252-261. 10.1016/j.jtbi.2005.11.036.

    Article  PubMed  Google Scholar 

  24. Ritchie MD, Hahn LW, Moore JH: Power of multifactor dimensionality reduction for detecting gene-gene interactions in the presence of genotyping error, missing data, phenocopy, and genetic heterogeneity. Genet Epidemiol. 2003, 24 (2): 150-157. 10.1002/gepi.10218.

    Article  PubMed  Google Scholar 

  25. Ritchie MD, Hahn LW, Roodi N, Bailey LR, Dupont WD, Parl FF, Moore JH: Multifactor-dimensionality reduction reveals high-order interactions among estrogen-metabolism genes in sporadic breast cancer. Am J Hum Genet. 2001, 69 (1): 138-147. 10.1086/321276.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Banerjee M, Vats P: Reactive metabolites and antioxidant gene polymorphisms in type 2 diabetes mellitus. Indian J Hum Genet. 2014, 20 (1): 10-19. 10.4103/0971-6866.132747.

    Article  PubMed  PubMed Central  Google Scholar 

  27. dos Santos KG, Canani LH, Gross JL, Tschiedel B, Souto KE, Roisenberg I: The catalase -262C/T promoter polymorphism and diabetic complications in Caucasians with type 2 diabetes. Dis Markers. 2006, 22 (5–6): 355-359.

    Article  PubMed  Google Scholar 

  28. Tarnai I, Csordas M, Sukei E, Shemirani AH, Kaplar M, Goth L: Effect of C111T polymorphism in exon 9 of the catalase gene on blood catalase activity in different types of diabetes mellitus. Free Radic Res. 2007, 41 (7): 806-811. 10.1080/10715760701381778.

    Article  CAS  PubMed  Google Scholar 

  29. Ambrosone CB, Freudenheim JL, Thompson PA, Bowman E, Vena JE, Marshall JR, Graham S, Laughlin R, Nemoto T, Shields PG: Manganese superoxide dismutase (MnSOD) genetic polymorphisms, dietary antioxidants, and risk of breast cancer. Cancer Res. 1999, 59 (3): 602-606.

    CAS  PubMed  Google Scholar 

  30. Cox DG, Tamimi RM, Hunter DJ: Gene x Gene interaction between MnSOD and GPX-1 and breast cancer risk: a nested case–control study. BMC Cancer. 2006, 6: 217-10.1186/1471-2407-6-217.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Ravn-Haren G, Olsen A, Tjonneland A, Dragsted LO, Nexo BA, Wallin H, Overvad K, Raaschou-Nielsen O, Vogel U: Associations between GPX1 Pro198Leu polymorphism, erythrocyte GPX activity, alcohol consumption and breast cancer risk in a prospective cohort study. Carcinogenesis. 2006, 27 (4): 820-825.

    Article  CAS  PubMed  Google Scholar 

  32. Tsai SM, Wu SH, Hou MF, Chen YL, Ma H, Tsai LY: Oxidative stress-related enzyme gene polymorphisms and susceptibility to breast cancer in non-smoking, non-alcohol-consuming Taiwanese women: a case–control study. Ann Clin Biochem. 2012, 49 (Pt 2): 152-158.

    Article  CAS  PubMed  Google Scholar 

  33. Cha MK, Suh KH, Kim IH: Overexpression of peroxiredoxin I and thioredoxin1 in human breast carcinoma. J Exp Clin Cancer Res. 2009, 28: 93-10.1186/1756-9966-28-93.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Penney RB, Roy D: Thioredoxin-mediated redox regulation of resistance to endocrine therapy in breast cancer. Biochim Biophys Acta. 2013, 1836 (1): 60-79.

    CAS  PubMed  Google Scholar 

  35. Bostrom B, Erdmann G: Cellular pharmacology of 6-mercaptopurine in acute lymphoblastic leukemia. Am J Pediatr Hematol Oncol. 1993, 15 (1): 80-86. 10.1097/00043426-199302000-00010.

    Article  CAS  PubMed  Google Scholar 

  36. Van Ness B, Ramos C, Haznadar M, Hoering A, Haessler J, Crowley J, Jacobus S, Oken M, Rajkumar V, Greipp P, Barlogie B, Durie B, Katz M, Atluri G, Fang G, Gupta R, Steinbach M, Kumar V, Mushlin R, Johnson D, Morgan G: Genomic variation in myeloma: design, content, and initial application of the Bank On A Cure SNP Panel to detect associations with progression-free survival. BMC Med. 2008, 6: 26-10.1186/1741-7015-6-26.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Hubackova M, Vaclavikova R, Ehrlichova M, Mrhalova M, Kodet R, Kubackova K, Vrana D, Gut I, Soucek P: Association of superoxide dismutases and NAD(P)H quinone oxidoreductases with prognosis of patients with breast carcinomas. Int J Cancer. 2012, 130 (2): 338-348. 10.1002/ijc.26006.

    Article  CAS  PubMed  Google Scholar 

  38. Rajaraman P, Hutchinson A, Rothman N, Black PM, Fine HA, Loeffler JS, Selker RG, Shapiro WR, Linet MS, Inskip PD: Oxidative response gene polymorphisms and risk of adult brain tumors. Neuro Oncol. 2008, 10 (5): 709-715. 10.1215/15228517-2008-037.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Gu D, Wang M, Zhang Z, Chen J: Lack of association between the hOGG1 Ser326Cys polymorphism and breast cancer risk: evidence from 11 case–control studies. Breast Cancer Res Treat. 2010, 122 (2): 527-531. 10.1007/s10549-009-0723-4.

    Article  CAS  PubMed  Google Scholar 

  40. Romanowicz-Makowska H, Smolarz B, Makowski M, Polac I, Pertynski T: Ser326Cys polymorphism in DNA repair genes hOGG1 in breast cancer women. Pol J Pathol. 2008, 59 (4): 201-204.

    CAS  PubMed  Google Scholar 

  41. Yuan W, Xu L, Feng Y, Yang Y, Chen W, Wang J, Pang D, Li D: The hOGG1 Ser326Cys polymorphism and breast cancer risk: a meta-analysis. Breast Cancer Res Treat. 2010, 122 (3): 835-842. 10.1007/s10549-009-0722-5.

    Article  CAS  PubMed  Google Scholar 

  42. Sangrajrang S, Schmezer P, Burkholder I, Waas P, Boffetta P, Brennan P, Bartsch H, Wiangnon S, Popanda O: Polymorphisms in three base excision repair genes and breast cancer risk in Thai women. Breast Cancer Res Treat. 2008, 111 (2): 279-288. 10.1007/s10549-007-9773-7.

    Article  CAS  PubMed  Google Scholar 

  43. Synowiec E, Stefanska J, Morawiec Z, Blasiak J, Wozniak K: Association between DNA damage, DNA repair genes variability and clinical characteristics in breast cancer patients. Mutat Res. 2008, 648 (1–2): 65-72.

    Article  CAS  PubMed  Google Scholar 

  44. Cai Q, Shu XO, Wen W, Courtney R, Dai Q, Gao YT, Zheng W: Functional Ser326Cys polymorphism in the hOGG1 gene is not associated with breast cancer risk. Cancer Epidemiol Biomarkers Prev. 2006, 15 (2): 403-404. 10.1158/1055-9965.EPI-05-0868.

    Article  CAS  PubMed  Google Scholar 

  45. Choi JY, Hamajima N, Tajima K, Yoo KY, Yoon KS, Park SK, Kim SU, Lee KM, Noh DY, Ahn SH, Choe KJ, Han W, Hirvonen A, Kang D: hOGG1 Ser326Cys polymorphism and breast cancer risk among Asian women. Breast Cancer Res Treat. 2003, 79 (1): 59-62. 10.1023/A:1023305826726.

    Article  CAS  PubMed  Google Scholar 

  46. Wei B, Zhou Y, Xu Z, Xi B, Cheng H, Ruan J, Zhu M, Hu Q, Wang Q, Wang Z, Yan Z, Jin K, Zhou D, Xuan F, Huang X, Shao J, Lu P: The effect of hOGG1 Ser326Cys polymorphism on cancer risk: evidence from a meta-analysis. PLoS One. 2011, 6 (11): e27545-10.1371/journal.pone.0027545.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Hirano T: Repair system of 7, 8-dihydro-8-oxoguanine as a defense line against carcinogenesis. J Radiat Res. 2008, 49 (4): 329-340. 10.1269/jrr.08049.

    Article  CAS  PubMed  Google Scholar 

  48. David SS, O’Shea VL, Kundu S: Base-excision repair of oxidative DNA damage. Nature. 2007, 447 (7147): 941-950. 10.1038/nature05978.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Friedman JI, Stivers JT: Detection of damaged DNA bases by DNA glycosylase enzymes. Biochemistry. 2010, 49 (24): 4957-4967. 10.1021/bi100593a.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Osorio A, Milne R, Kuchenbaecker K, Vaclová T, Pita G, Alonso R, Peterlongo P, Blanco I: DNA glycosylases involved in Base Excision Repair may be associated with cancer risk in BRCA1 and BRCA2 mutation carriers. Plos Genet. 2014, 10 (4): e1004256-10.1371/journal.pgen.1004256.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Andrew AS, Karagas MR, Nelson HH, Guarrera S, Polidoro S, Gamberini S, Sacerdote C, Moore JH, Kelsey KT, Demidenko E, Vineis P, Matullo G: DNA repair polymorphisms modify bladder cancer risk: a multi-factor analytic strategy. Hum Hered. 2008, 65 (2): 105-118. 10.1159/000108942.

    Article  CAS  PubMed  Google Scholar 

  52. Gui J, Andrew AS, Andrews P, Nelson HM, Kelsey KT, Karagas MR, Moore JH: A robust multifactor dimensionality reduction method for detecting gene-gene interactions with application to the genetic analysis of bladder cancer susceptibility. Ann Hum Genet. 2011, 75 (1): 20-28. 10.1111/j.1469-1809.2010.00624.x.

    Article  PubMed  Google Scholar 

  53. Ihsan R, Chauhan PS, Mishra AK, Yadav DS, Kaushal M, Sharma JD, Zomawia E, Verma Y, Kapur S, Saxena S: Multiple analytical approaches reveal distinct gene-environment interactions in smokers and non smokers in lung cancer. PLoS One. 2011, 6 (12): e29431-10.1371/journal.pone.0029431.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Su MW, Tung KY, Liang PH, Tsai CH, Kuo NW, Lee YL: Gene-gene and gene-environmental interactions of childhood asthma: a multifactor dimension reduction approach. PLoS One. 2012, 7 (2): e30694-10.1371/journal.pone.0030694.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Chen X, Wang J, Guo W, Liu X, Sun C, Cai Z, Fan Y, Wang Y: Two functional variations in 5’-UTR of hoGG1 gene associated with the risk of breast cancer in Chinese. Breast Cancer Res Treat. 2011, 127 (3): 795-803. 10.1007/s10549-010-1284-2.

    Article  CAS  PubMed  Google Scholar 

  56. He C, Tamimi RM, Hankinson SE, Hunter DJ, Han J: A prospective study of genetic polymorphism in MPO, antioxidant status, and breast cancer risk. Breast Cancer Res Treat. 2009, 113 (3): 585-594. 10.1007/s10549-008-9962-z.

    Article  CAS  PubMed  Google Scholar 

  57. Oestergaard MZ, Tyrer J, Cebrian A, Shah M, Dunning AM, Ponder BA, Easton DF, Pharoah PD: Interactions between genes involved in the antioxidant defence system and breast cancer risk. Br J Cancer. 2006, 95 (4): 525-531. 10.1038/sj.bjc.6603272.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. McCullough LE, Santella RM, Cleveland RJ, Bradshaw PT, Millikan RC, North KE, Olshan AF, Eng SM, Ambrosone CB, Ahn J, Steck SE, Teitelbaum SL, Neugut AI, Gammon MD: Polymorphisms in oxidative stress genes, physical activity, and breast cancer risk. Cancer Causes Control. 2012, 23 (12): 1949-1958. 10.1007/s10552-012-0072-1.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Rossner P, Terry MB, Gammon MD, Zhang FF, Teitelbaum SL, Eng SM, Sagiv SK, Gaudet MM, Neugut AI, Santella RM: OGG1 polymorphisms and breast cancer risk. Cancer Epidemiol Biomarkers Prev. 2006, 15 (4): 811-815. 10.1158/1055-9965.EPI-05-0659.

    Article  CAS  PubMed  Google Scholar 

  60. Li Y, Ambrosone CB, McCullough MJ, Ahn J, Stevens VL, Thun MJ, Hong CC: Oxidative stress-related genotypes, fruit and vegetable consumption and breast cancer risk. Carcinogenesis. 2009, 30 (5): 777-784. 10.1093/carcin/bgp053.

    Article  CAS  PubMed  Google Scholar 

  61. Crawford A, Fassett RG, Geraghty DP, Kunde DA, Ball MJ, Robertson IK, Coombes JS: Relationships between single nucleotide polymorphisms of antioxidant enzymes and disease. Gene. 2012, 501 (2): 89-103. 10.1016/j.gene.2012.04.011.

    Article  CAS  PubMed  Google Scholar 

  62. Weiss JM, Goode EL, Ladiges WC, Ulrich CM: Polymorphic variation in hOGG1 and risk of cancer: a review of the functional and epidemiologic literature. Mol Carcinog. 2005, 42 (3): 127-141. 10.1002/mc.20067.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references


This work has been supported by grants PS09/01700 and PI12/01421 (Ministerio de Ciencia y Tecnología-Fondo de Investigación Sanitaria del Instituto de Salud Carlos III-FEDER) and RD12/0036/0070 (RTICC) to A.LL. and by grants GE-004/09 and ACOMP/2009/201 (Consellería de Sanidad Valenciana). JF was funded from the RTICC RD12/0036/0070, and PR has been hired under the Santiago Grisolia program. CEGEN is funded by the Instituto de Salud Carlos III.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Pilar Eroles.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

PR processed the material, extracted DNA, prepared samples for SNPlex, analyzed all the data, performed statistical analyses, carried out literature searches and contributed in the drafting of the manuscript. GdM processed samples, optimized the technique used, contributed to the analysis of the data and revised the manuscript. JF processed samples and helped in the data analyses, interpretation of the results and in drafting the manuscript. MLM performed the SNPlex protocol, analyzed the data and contributed in the drafting of the manuscript. MPA and AG-N helped in the sample processing and interpretation of the results and revised the final manuscript. JC M-E and JB performed the selection of the patients, provided intellectual content and revised the final manuscript; ALL performed the conceptual design, performed the selection of the patients, supervised the study, revised the manuscript and participated in the acquisition of funding. FJC performed the conceptual design of the study, revised the manuscript and provided intellectual content. PE conceived and supervised the study, contributed to the interpretation of the data, gave intellectual support, contributed to writing the manuscript and helped in the acquisition of funding. All authors have given final approval to the version to be published.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rodrigues, P., de Marco, G., Furriol, J. et al. Oxidative stress in susceptibility to breast cancer: study in Spanish population. BMC Cancer 14, 861 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: