Rad51 paralogs and the risk of unselected breast cancer: A case-control study

Authors: Peter Grešner aff001;  Ewa Jabłońska aff002;  Jolanta Gromadzińska aff003
Authors place of work: Department of Toxicology and Carcinogenesis, Nofer Institute of Occupational Medicine, Lodz, Poland aff001;  Department of Molecular Genetics and Epigenetics, Nofer Institute of Occupational Medicine, Lodz, Poland aff002;  Department of Biological and Environmental Monitoring, Nofer Institute of Occupational Medicine, Lodz, Poland aff003
Published in the journal: PLoS ONE 15(1)
Category: Research Article
doi: 10.1371/journal.pone.0226976


A case-control study was conducted in which we evaluated the association between genetic variability of DNA repair proteins belonging to the Rad51 family and breast cancer (BrC) risk. In the study, 132 female BrC cases and 189 healthy control females were genotyped for a total of 14 common single nucleotide polymorphisms (SNPs) within Rad51 and Xrcc3. Moreover, our previously reported Rad51C genetic data were involved to explore the nonlinear interactions among SNPs within the three genes and effect of such interactions on BrC risk. The rare rs5030789 genotype (-4601AA) in Rad51 was found to significantly decrease the BrC risk (OR = 0.5, 95% CI: 0.3–1.0, p<0.05). An interaction between this SNP, rs2619679 and rs2928140 (both in Rad51), was found to result in a two three-locus genotypes -4719AA/-4601AA/2972CG and -4719AT/-4601GA/2972CC, both of which were found to increase the risk of BrC (OR = 8.4, 95% CI: 1.8–38.6, p<0.0001), instead. Furthermore, rare Rad51 rs1801320 (135CC) and heterozygous Xrcc3 rs3212057 (10343GA) genotypes were found to respectively increase (OR = 10.6, 95% CI: 1.9–198, p<0.02) and decrease (OR = 0.0, 95% CI: 0.0-NA, p<0.05) the risk of BrC. Associations between these SNPs and BrC risk were further supported by outcomes of employed machine learning analyses. In Xrcc3, the 4541A/9685A haplotype was found to be significantly associated with reduced BrC risk (OR = 0.5; 95% CI: 0.3–0.9; p<0.05). Concluding, our study indicates a complex role of SNPs within Rad51 (especially rs5030789) and Xrcc3 in BrC, although their significance with respect to the disease needs to be further clarified.


Haplotypes – Molecular genetics – machine learning – Variant genotypes – Genetic polymorphism – Introns – DNA repair – breast cancer


Breast cancer (BrC) is known to be the most common malignancy among women, with nearly 1.7 million new cases and more than 520,000 deaths per year worldwide [1]. Its incidence is higher in North America and Western European countries contrary to Asian or African populations. Although not fully elucidated, mechanisms leading to BrC include a number of genetic and environmental factors, family history of the disease, multiparity, early menarche and late menopause [2].

Genetic association and GWAS studies provided valuable insights into genetic factors contributing to BrC risk. In addition to major BrC susceptibility genes including BRCA1 and BRCA2, other high- (such as TP53 and PTEN) and moderate- (CHEK2, ATM, BRIP1, PALB2, and RAD51C) penetrance susceptibility genes were found to play role in the onset of BrC [35]. Both BRCA genes together with the above-mentioned moderate-penetrance BrC susceptibility genes play their roles in homologous recombination (HR) DNA repair pathway involved in repair of DNA double strand breaks (DSB) [58]. It has been proposed that compromised capacity of the HR DNA repair system leads to increased accumulation of DNA damage, mutations and, hence, increased risk of malignancies [911].

The key component of the HR DNA DSB repair pathway is comprised by the Rad51 family proteins including Rad51, a crucial player in the whole HR DNA DSB repair machinery, and its five paralogs—Rad51B, Rad51C, Rad51D, Xrcc2, and Xrcc3. Paralogs interact with each other to form a hetero-tetrameric (Rad51B/Rad51C/Rad51D/Xrcc2; BCDX2) and hetero-dimeric (Rad51C/Xrcc3; CX3) complexes crucial for various processes involved in the HR DNA DSB repair machinery [12].

Rad51 (RecA homolog, Escherichia coli; 15q15.1) is a homolog of bacterial RecA protein forming a nucleoprotein filament on single-stranded DNA which in turn mediates the strand invasion and exchange between the damaged DNA sequence and its undamaged homologue thus facilitating the re-synthesis of damaged DNA region [13]. Xrcc3 (X-ray repair cross-complementing group 3; 14q32.3), on the other hand, has been shown to be crucial with respect to accumulation of Rad51 at sites of DNA DSB in the cell nucleus as well as to enzymatic resolution of the resultant cross-stranded structure (the Holiday junction) [13]. Rad51C (Rad51 homolog C, S. cervisiae; 17q25.1) seems to be required for RAD51/DNA nucleoprotein filament formation as it localizes to DNA DSB sites in early stages of HR [14], but it is also involved in DNA damage response and checkpoint activation [14], migration and resolution of Holiday junction [15], repair of interstrand cross-links [16] and stalled/collapsed replication forks as well as in antioxidant protection of mitochondrial genome [17,18]. Finally, its function as tumor suppressor and cancer susceptibility gene [5,1921] has been proposed.

Most cancer association studies involving Rad51 were focused on two single nucleotide polymorphisms (SNPs) localized in the 5′ untranslated region (5’UTR) of exon 1 of the gene: rs1801320 (c.−98G>C; 135G/C) and rs1801321 (c.−61G>T; 172G/T). Both these SNPs were reported to be associated with altered gene transcription [22,23]. Large meta-analyses have shown that 135C allele increases the general risk as well as the risk of BrC, with a distinct dose-dependent effect [10,24]. The 172T allele-containing genotypes, on the other-hand, were found to be associated with some 25% reduction of general odds of cancer compared to the 172GG wild-type one [10]. Nevertheless, association between 172T allele and the risk of BrC seems to be much more complex, as only limited number of studies are available so far, suggesting both increase [25,26] as well as decrease [27,28] of the disease risk being associated with the allele, not allowing us thus to draw any final conclusion.

In the case of Xrcc3, rs861539 (c.722C>T; 241Thr/Met) in exon 8 is the most frequently tested SNP with respect to cancer risk. The risk of BrC among carriers of the 241Met-containg genotypes was found to be increased compared to wild-type carriers by some 6–10% under various genetic models [2931]. Nevertheless, recent smaller studies conducted in Polish BrC population failed to provide evidence on any unambiguous effect with respect to BrC risk [32,33]. Other studies proposed the 17893G allele (rs1799796; intron 7; c.562-14A>G; 17893A/G) as providing protective effect against BrC (risk reduction of some 10%) [31].

Missense mutations in Rad51C were found to associate with hereditary breast and ovarian cancer (HBOC), which has been further confirmed in several subsequent studies on unselected ovarian cancer (OC). Nevertheless, there is still a considerable amount of studies which failed to find any association between Rad51C mutations and HBOC, what has usually been explained by very rare occurrence of these mutations [21,3438]. Interestingly, none of the above cited studies identified Rad51C mutations associated with the BrC-only families.

The above cited reports prompted us to further contribute to our previous report on associations between genetic variability of Rad51C and the risk of BrC [39] by evaluating the associations between genetic variability of two more enzymes belonging to the Rad51 family—Rad51 and Xrcc3—and the risk of BrC. To this end, a total of 14 common SNPs in Rad51 and Xrcc3 (seven per each gene) were genotyped and tested for significant differences in their distributions between female BrC cases and controls. In addition to conventional analyses (single-site, SNP combinations and haplotype analyses), machine learning (ML) techniques (random forest and multifactor dimensionality reduction) providing increased statistical power and novel approach to cancer association studies [40] were used to explore main and nonlinear (epistatic) interactional effects of SNPs with respect to their association with BrC. In ML analyses, hereby described Rad51 and Xrcc3 genotypic data, supplemented with our previously reported genotypic data obtained for Rad51C [39], were used to broaden the picture of involvement of the Rad51 family in BrC.

Material and methods


In this study, previously described groups of breast cancer patients and control subjects were used [39]. Briefly, the study group consisted of 132 female breast cancer patients of European descent aged 36–86 years (median age at the time of diagnosis of 57 years; interquartile range (IQR): 15 years) hospitalized at the Department of Oncology, Memorial Copernicus Hospital in Lodz, Poland in years 2007 and 2008 with histopathologically confirmed diagnosis of BrC. Only female patients with primary breast cancer tumor without metastases and without any history of previous anti-cancer treatment, undergoing curative resection therapy or chemotherapy, were enrolled. The control group consisted of 189 healthy cancer-free volunteer females of European descent, aged 35–54 years (median age at the time of examination of 43 years; IQR: 6 years), willing to undergo examinations. All subjects enrolled in the study were residents of Lodz district in central Poland.

Additional information on tobacco-smoking habits was collected for both controls and BrC cases and the individuals were classified as either never- or ever-smoker according to the criterion suggested by Pomerleau et al. [41]. According to this criterion, only those subjects who have smoked less than 19 cigarettes (a pack) during their lifetime, were classified as never-smokers, while the others were considered as ever-smokers. No data concerning the alcohol consumption was available for either controls or cases.

Written and informed consent for participation in this study was obtained from each subject enrolled prior to any experiments. The study was performed under the guidelines of the Helsinki Declaration for human research and was approved by the Bioethics Committee in the Nofer Institute of Occupational Medicine (resolution no. 5/2007). Characteristics of the breast cancer group and the control group are summarized in Table 1.

Tab. 1. Characteristics of the groups of subjects involved in the study.
Characteristics of the groups of subjects involved in the study.

DNA isolation

Peripheral blood leukocytes were used to isolate the genomic DNA using the QIAamp DNA Blood Mini Kit (Qiagen, Germany). Manufacturer’s instructions were followed. RNA contamination was removed by digestion with 1 mg/ml RNase A (Qiagen, Hilden, Germany). DNA quantification, its purity and protein content were assessed using an Eppendorf BioPhotometer (Eppendorf, Hamburg, Germany) instrument. All DNA samples were stored at -80°C until further processing.


Rad51 paralogs show a relatively high degree of conservativeness (for gene and protein conservativeness schemes see Fig 1 and S1 Fig) with missense changes being very rare within these genes, therefore, focus of this study was put on SNPs occurring predominantly in non-coding regions which may plausibly be involved in regulation of gene expression. A total of 14 SNPs with seven of them being localized in Rad51 and additional seven in Xrcc3 were analyzed. In the case of Rad51, SNPs occurring in the promoter, 5′UTR and intron 3 with the minor allele frequency (MAF) in the Caucasian population exceeding 10% (according to the dbSNP database [42]) were selected. In the case of Xrcc3, in addition to frequently analyzed rs1799794 (4541A/G), rs1799796 (17893A/G) and rs861539 (p.Thr241Met), SNPs localized in 5’UTR, intron 5, together with two missense SNPs in exons 6 and 10, for which the minor allele frequency in the Caucasian population was higher than 10%, were selected. Detailed information on SNPs analyzed in this study are provided in Table 2, while additional graphic representation of localization of all analyzed SNPs is provided in Fig 1. Plausible effects of selected SNPs in non-coding regions, as predicted by the PERFECTOS-APE in-silico method for prediction of regulatory functional effect of noncoding SNPs [43], are provided in Table 3. Prediction of classification of non-synonymous SNPs analyzed in the study, obtained by Polyphen-2 [44], can also be found in Table 3.

<i>Rad51</i> and <i>Xrcc3</i> gene and protein structures.
Fig. 1. Rad51 and Xrcc3 gene and protein structures.
Rad51 and Xrcc3 gene and protein structures are shown with graphic representation of the localization of all 14 SNPs investigated in this study. The positions of individual SNPs in respective genes are indicated by arrows annotated with respective working codes used in this study. Localizations of rs3212057 (X4), rs861539 (X6) and rs28903081 (X7) are indicated also within the Xrcc3 protein structure, as these were the only three SNPs localized within coding regions of exons. Rad51 and Xrcc3 are placed on chromosomes 15 and 14 and span 37 and 18 kb, respectively. Their protein products consist of 339 and 346 amino acids, respectively, showing high degree of conservativeness. They both share highly conserved Rad51-like C-terminal domain with approx. 250 amino acids (responsible for DNA binding) within which two subdomains can be distinguished: the RecA-like ATP-binding domain (responsible for ATP binding and hydrolysis) and the RecA-like monomer-monomer interface. The RecA-like ATP-binding domain embodies two highly conserved consensus motifs (Walker A and B) conferring ATP-binding and -hydrolysis activities. Distributed within this domain is also the multimere BRC interface (pink lines in the protein structure) facilitating interaction with DNA. X4, X6 and X7 are all localized within Rad51-like C-terminal domain, while the rs861539 (X6) hits exactly the threonine-241 taking part in interaction with DNA. Rad51 consists additionally of the helicase-hairpin-helicase (HhH) domain responsible for nonspecific DNA binding [42,45].
Tab. 2. Resume of SNPs genotyped in the study.
Resume of SNPs genotyped in the study.
Tab. 3. Predicted effects of SNPs analyzed in the study.
Predicted effects of SNPs analyzed in the study.

All Xrcc3 and all but one (rs1801321) Rad51 SNPs were genotyped using the PCR-restriction fragment length polymorphism (PCR-RFLP) technique on a BioRad’s PTC-200 DNA Engine thermal cycler (BioRad, Hercules, CA, USA) instrument utilizing the Qiagen’s HotStarTaq PCR kit (Qiagen, Hilden, Germany). Used primer sequences, together with their basic characteristics and used PCR conditions, have already been provided in full details elsewhere [46], but are listed again for convenience in Table 4.

Tab. 4. Summary on primer sequences, amplicons, PCR cycling conditions and enzymatic digestion for individual SNPs genotyped in the study [46].
Summary on primer sequences, amplicons, PCR cycling conditions and enzymatic digestion for individual SNPs genotyped in the study [<em class="ref">46</em>].

Rs1801321 in Rad51 was genotyped by the real-time PCR technique using the predesigned commercially available TaqMan SNP Genotyping Assay kit (Life Technologies, Carlsbad, CA, USA) and detailed conditions of the reaction have also been previously provided [46].

In the case of Rad51C, no SNPs within this gene were genotyped in this study. Instead, genotypic data of a set of eight SNPs previously described in our recent study on breast cancer [39] were used. These SNPs included rs302874 (hereby designated RcA), rs12946522 (RcB), rs302873 (RcC), rs16943176 (RcD), rs12946397 (RcE), rs28910276 (RcF), rs17222691 (RcG) and rs28363302 (RcH). More details on genotyping and the results of analysis of associations between these SNPs and BrC can be find elsewhere [39].

Statistical analyses

Single-locus analyses

For all SNPs, both absolute and relative genotypic frequencies are provided. The Hardy-Weinberg equilibrium (HWE) in controls was tested by a goodness-of-fit chi-square test. For each investigated SNP, possible associations with BrC risk at both allelic and genotypic level were sought for by Fisher’s exact test (allelic level) or unconditional logistic regression (genotypic level) in a series of separate univariate (i.e. single-site) analyses. Associations with BrC are expressed as either raw (allelic level) or age-and-smoking-status-adjusted (genotypic level) odds ratios (ORs) with corresponding 95% confidence intervals (95% CI). Dominant, recessive, and additive genetic model, together with direct comparison of variant versus wild-type homozygotes were assumed. Significance was inferred for p<0.05.

Analysis of Rad51/Xrcc3 SNP combinations

Rad51 and Xrcc3 SNP with statistically significant outcomes from single-site analyses were further involved in analysis of association between their mutual combinations and BrC. BrC risk associated with combinations of the so-called high-risk genotypes in these polymorphic sites was estimated by means of age-and-smoking-status-adjusted unconditional logistic regression and expressed as ORs with corresponding 95% CIs.

Analysis of Rad51 and Xrcc3 haplotypes

Linkage disequilibrium (LD) and haplotype reconstruction were performed by means of an expectation-maximization algorithm implemented in the Haploview package [47]. Briefly, “strong LD” blocks were recognized based on normalized measure of allelic association |D’| [48] according to the confidence interval method proposed by Gabriel et al. [49]. As earlier described [39], haplotypes reconstructed within each “strong LD” block were tested for differences in their frequencies between the control and cancer group and the significance of such differences was assessed using a two-sided exact mid-P test. Possible linkage of haplotypes with BrC was expressed by means of OR with corresponding 95% CI and significance was inferred for p<0.05.

Analyses based on machine learning techniques

In these analyses, SNPs for which no variability was found in single-locus analyses (rs45603942 (X2) and rs28903081 (X7) in Xrcc3; rs28910276 (RcF) in Rad51C [39]) were omitted and subjects with any lacks in data were excluded. Therefore, a total of 19 SNPs (7 SNPs in Rad51, 5 SNPs in Xrcc3 and 7 SNPs in Rad51C) were involved.

RF-based analysis of associations between predictors and BrC

In addition to above described single-locus analyses, simple associations between BrC and analyzed SNPs were assessed by means of the random forest (RF) machine learning strategy. The Breiman-Cutler permutation variable importance (VIMP) [50] was used to measure and rank the strength of such associations. In our study, we used the randomForestSRC package for R obtained from the CRAN repository [51], using which we employed a robust strategy allowing us to reliably validate and statistically infer on the ranking of all analyzed SNPs with respect to their ability to accurately predict the BrC/control status. Detailed description of this strategy is provided in Part A of the S1 File.

The whole procedure was performed twice. In the first run, only SNPs were considered as possible predictors of the BrC case/control status, which provided us with “raw” results, while the second one considered SNPs together with subjects’ age (dichotomized with respect to median age) and smoking status (never/ever smoker), providing the VIMP-based ranking of SNPs allowed for interactions with these two common confounders. Levels of significance were obtained by permutation testing (see the Part A in S1 File for further details) and statistical significance was inferred for p<0.05.

Random forest analysis of epistatic interactions

To distinguish between main and interactive (i.e. epistatic) effects of SNPs on BrC, direct analysis of pure epistatic interactions and their association with BrC was performed using a permutation-based machine learning strategy relying on RF methodology termed the permuted random forest (pRF). pRF detects and quantifies pure interaction between selected SNPs and estimates how much it contributes to the model prediction power. We have implemented this method in R using the randomForestSRC [51] and permutations [52] packages obtained from the CRAN repository according to a thorough description of algorithm provided by Li et al. [53] with minor modifications allowing us to perform analysis of interactions between pairs as well as among triplets of SNPs (i.e. 2-way and 3-way interactions). Subjects’ age (dichotomized with respect to age median) and smoking status (never/ever smoker) were also involved in the analysis as possible confounders. The same RF model as the one used to obtain the VIMP-based ranking of simple associations between SNPs and BrC was used in this analysis. All possible 2-way and 3-way combinations of predictors (SNPs and confounders) were analyzed and the strength of associations between individual combinations of predictors and BrC was measured by the so-called differential error (ΔE; see the Part B in the S1 File). The predictor combination with the highest value of ΔE was considered as the best one being in the strongest association with BrC. Detailed description of the algorithm used can be found in Part B of the S1 File.

Analysis of epistatic interactions using dimensionality reduction approach

As an alternative approach for elucidating the epistatic interactions among predictors (SNPs and subject’s age and smoking status as confounders) and their associations with BrC, the model-based multifactor dimensionality reduction (MB-MDR) was used.

The algorithm was implemented using the mbmdr package for R obtained from the CRAN repository [54] and is described in more details in Part C of the S1 File. MB-MDR uses a constructive induction technique to merge multi-locus genotypes into a one-dimensional construct, assigning each analyzed combination of genotypes to either “high-risk”, “low-risk”, or a “no-evidence” (or “non-informative”) category. Such new predictive variable with three states (H, L, 0) was then tested for association with the risk of BrC and such association then expressed by means of Wald statistic, OR and respective p value (separately for the “high-risk” and “low-risk” category, where appropriate). Only 2-way and 3-way interactions were analyzed, and permutation testing was used to correct the obtained p-levels for multiple hypotheses testing.


Single-locus analyses of associations between Rad51 and Xrcc3 SNPs and breast cancer

Single site analyses of 14 SNPs within Rad51 and Xrcc3 in terms of their possible associations with BrC were performed. All relevant results including the counts (frequencies) in the BrC and control groups together with logarithmic regression-derived ORs adjusted to age and smoking status are presented in Tables 5A and 6A. For all investigated SNPs, observed genotype frequencies in control groups were in agreement with those predicted by the Hardy-Weinberg law.

Tab. 5. Observed genotype and haplotype frequencies together with respective odds ratios (ORs) and corresponding 95% confidence intervals for Rad51 polymorphisms in breast cancer cases and control subjects.
Observed genotype and haplotype frequencies together with respective odds ratios (ORs) and corresponding 95% confidence intervals for <i>Rad51</i> polymorphisms in breast cancer cases and control subjects.
Tab. 6. Observed genotype and haplotype frequencies together with respective odds ratios (ORs) and corresponding 95% confidence intervals for Xrcc3 polymorphisms in breast cancer cases and control subjects.
Observed genotype and haplotype frequencies together with respective odds ratios (ORs) and corresponding 95% confidence intervals for <i>Xrcc3</i> polymorphisms in breast cancer cases and control subjects.

Out of the 7 Rad51 SNPs analyzed in this study, only rs5030789 (-4601A/G; RB) and rs1801320 (135G/C; RC) showed differences in genotype frequencies between the BrC and control groups.

As long as the rs5030789 is concerned, the rare -4601AA genotype was found to be almost twice less common among BrC cases compared to controls (12.5% vs. 23.4%), hence it was associated with significantly decreased BrC risk under recessive genetic model (OR = 0.5, 95%CI: 0.3–1.0; p<0.05) as well as in direct comparison between wild-type and variant homozygotes (OR = 0.5, 95% CI: 0.2–1.0; p<0.05). According to outcomes of the analysis under additive genetic model, each copy of the -4601A allele was associated with approximately 30% reduction in odds of BrC, although this outcome remained just beyond the edge of statistical significance (OR = 0.7, 95% CI: 0.5–1.0; p = 0.06). It is, however, considerably similar to an outcome obtained from direct comparison of allelic frequencies, in which the -4601A allele was also found to be less frequent among BrC cases compared to controls (40.2% vs. 47.8%) and thus rendering some 30% reduction in BrC risk (OR = 0.7, 95% CI: 0.5–1.0; p = 0.07), yet just beyond the edge of statistical significance, too. No statistically significant outcomes were obtained for this SNP under dominant genetic model (Table 5A).

Concerning the rs1801320, the rare 135CC genotype was found to be significantly more abundant among BrC cases (7.1% vs. 0.8%) hence can be assumed as associated with increased risk of BrC under recessive genetic model (OR = 10.6, 95% CI: 1.9–198; p<0.05) as well as in direct comparison between wild-type and variant homozygotes (OR = 9.8, 95% CI: 1.8–184; p<0.05). However, no statistically significant outcomes were found under dominant or additive genetic models as well as in direct comparison of allelic frequencies for this SNP (Table 5A).

In the case of Xrcc3, two out of seven analyzed SNPs (rs45603942 (4576C/T; X2) and rs28903081 (c.905G>A; X7)) did not show any genetic variability as only the wild-type homozygotes were observed in both the BrC and control groups. Moreover, for rs3212057 (10343G/A; X4) no variant homozygotes were observed in this study. Although the distribution of genotypes for this SNP was found to differ between BrC cases and controls, with heterozygotes being abundant only among controls compared to BrC cases (3.9% vs. 0.0%), logistic regression adjusted to age and smoking-status did not, however, render this difference as statistically significant under either dominant or additive genetic model with respective ORs being unavailable due to zero frequency of heterozygotes among BrC subjects. Recessive genetic model and direct comparison between wild-type and variant homozygotes were not possible for this SNP. The only statistically significant outcome for this SNP was thus observed at allelic level, where the 10343A allele was found to be associated with reduced risk of BrC (OR = 0.0, 95% CI unavailable; p<0.05; Table 6A). No statistically significant outcomes for any other Xrcc3 SNPs were found under any of the models examined.

Linkage disequilibrium and haplotype analysis

Analysis of non-random associations between the investigated Rad51 SNPs revealed a 1-kb long block of “strong LD” spanning from rs1801321 (172G/T; RD) in 5’UTR of exon 1, through rs2619680 (1037A/C; RE) in intron 1, to rs2619681 (1640C/T; RF) in intron 1 (Fig 2A). Within this LD block, four common (172T/1037A/1640C; 172G/1037C/1640C; 172G/1037A/1640C; 172G/1037C/1640T) and two rare (172T/1037C/1640C, 172G/1037A/1640T) haplotypes were reconstructed. Common haplotypes encompassed together 99.4% of all subjects (Table 5B). Nevertheless, no significant associations with BrC risk were found for any of the haplotypes reconstructed.

The map of LD among analyzed SNPs in (A) <i>Rad51</i> and (B) <i>Xrcc3</i>.
Fig. 2. The map of LD among analyzed SNPs in (A) Rad51 and (B) Xrcc3.
In the case of Xrcc3, rs45603942 (X2) and rs28903081 (X7) are not shown as these SNPs were not considered in haplotype analysis due to lack of any genetic variability observed in these loci. For a given pair of SNP, the values in the map are the values of normalized measure of allelic association (|D’|, provided as percentages) and the color scheme represents the corresponding confidence bounds: dark gray—strong evidence of LD; light-gray—inconclusive; white—strong evidence of recombination [48]. Identified LD blocks are indicated by solid lines encompassing respective SNPs.

In the case of Xrcc3, rs45603942 (4576C/T; X2) in 5’UTR of exon 1 and rs28903081 (c.905G>A; X7) in exon 7 were not included in the LD analysis due to lack of any observed variability in our study. The LD analysis revealed two blocks of “strong LD”, one spanning across 5 kb from rs1799794 (4541A/G; X1) in 5’UTR of exon 2 to rs861530 (9685A/G; X3) in intron 5, while the other one spanning over 174 bp from rs1799796 (17893A/G; X5) in intron 7 to rs861539 (18067C/T; X6) in exon 8 (Fig 2B).

Within the first LD block (X1-X3), three common (4541A/9685G, 4541G/9685A, 4541A/9685A) and one rare (4541G/9685G) haplotype were reconstructed, with the common haplotypes encompassing 99.7% of all subjects (Table 6B). Out of these haplotypes, only the 4541A/9685A common haplotype was found to be significantly associated with reduced risk of BrC, as it was significantly less abundant among BrC cases compared to controls, resulting in around 2-fold reduction of the odds of BrC comparing to carriers of all other haplotypes together (5.6% vs. 11.0% among BrC cases and controls, respectively; OR = 0.5; 95% CI: 0.3–0.9; p<0.05). Within the second (X5-X6) LD block, three common (17893G/18067C,17893A/18067T, 17893A/18067C) and one rare (17893G/18067T) haplotypes were reconstructed, with common haplotypes encompassing 99.5% of all subjects. Nevertheless, none of the haplotype reconstructed within this LD block was associated with the risk of BrC (Table 6B).

Associations between Rad51/Xrcc3 SNP combinations and BrC

Only those SNPs for which statistically significant outcomes in single-site analyses were revealed were involved in this analysis (i.e. rs5030789 (-4601A/G; RB), rs1801320 (135G/C; RC) and rs3212057 (10343G/A; X4)). Combinations of respective high-risk genotypes (i.e. rs5030789 (RB) -4601G/G or -4601G/A, rs1801320 (RC) 135C/C and rs3212057 (X4) 10343G/G) were tested for association with BrC against respective genotype combinations encompassing the lowest possible number of high-risk genotypes (which was 0 in the case of RB/RC, RC/X4 combinations and 1 in the case of RB/X4, RB/RC/X4 combinations, as for these combinations no subjects with low-risk-only genotype combinations were found). Obtained results are summarized in Table 7.

Tab. 7. Associations of Rad51/Xrcc3 SNP combinations with breast cancer.
Associations of <i>Rad51</i>/<i>Xrcc3</i> SNP combinations with breast cancer.

Significant associations with BrC risk were found for all three possible two-way combinations (i.e. rs5030789 (RB) -4601G/G or -4601G/A & rs1801320 (RC) 135C/C, p<0.002; rs5030789 (RB) -4601G/G or -4601G/A & rs3212057 (X4) 10343G/G, p<0.02; rs1801320 (RC) 135C/C & rs3212057 (X4) 10343G/G, p<0.02) as well as for the three-way combination (rs5030789 (RB) -4601G/G or -4601G/A & rs1801320 (RC) 135C/C & rs3212057 (X4) 10343G/G, p<0.002) of high-risk genotypes listed above, with generally higher levels of significance compared to p values obtained when individual SNPs were tested in single-site analyses. Carriers of these combinations of high-risk genotypes were at increased risk of BrC comparing to those carrying the genotype combinations with maximum possible number of low-risk genotypes (RB&RC: OR = 7.3; 95% CI: 2.1–25.8; RB&X4: OR = 2.2; 95% CI: 1.2–4.3; RC&X4: OR = 184.9; 95% CI: 3.1–11172.8; RB&RC&X4: OR = 8.1; 95% CI: 2.3–28.9). Detailed case/control rates are provided in Table 7.

RF-based analysis of associations between analyzed SNPs and BrC

As for rs45603942 (X2) and rs28903081 (X7) in Xrcc3 as well as for rs28910276 (RcF) in Rad51C no variability was found in single-locus analyses, 7 SNPs in Rad51, 5 SNPs in Xrcc3 and 7 SNPs in Rad51C were further included in RF-based analysis of the strength of their associations with BrC. Intermediate outcomes based on which the best RF-based models were selected are provided in the S2 Table. Performance characteristics of these RFs are provided in S3 Table.

When only SNPs were tested for associations with BrC (i.e. without covariates), obtained VIMP values (Table 8) suggest rs5030789 (RB), rs1801321 (RD), rs1801320 (RC), rs861530 (X3), and rs2928140 (RG) as the five best BrC/control predictors with the strongest association with BrC/control status. Moreover, the VIMPs of rs5030789 (RB), rs1801321 (RD), rs1801320 (RC) and rs3212057 (X4) were all found to be statistically significant (p<0.005 for RB; p<0.05 for RD, RC, and X4), while the VIMP value of rs861530 (X3) remained close to the edge of statistical significance (p = 0.074). The results of bootstrapping the RF-based ranking procedure are shown in Fig 3A. It is of note that the rs5030789 (RB) SNP was ranked 1st in the vast majority of all bootstrapped RF models (99.6% of all 10,000 RFs; Part A in S4 Table) and based on derived weighted average ranks it seems to be confirmed as the best predictor with the strongest association with the BrC/control status. Besides that, six predictors with the lowest weighted average ranks in Fig 3A (i.e. the six best ones) were the same as in the ranking obtained based on observed VIMP values (rs5030789 (RB), rs1801321 (RD), rs1801320 (RC), rs861530 (X3), rs2928140 (RG), and rs1799794 (X1)).

Tab. 8. Ranking of all analyzed predictors (SNPs, covariates) according to their VIMP-based association with BrC.
Ranking of all analyzed predictors (SNPs, covariates) according to their VIMP-based association with BrC.
Results of bootstrap analysis of the VIMP-based ranking of predictors.
Fig. 3. Results of bootstrap analysis of the VIMP-based ranking of predictors.
Bootstrap estimate of distribution of the VIMP-based ranks of analyzed predictors obtained either (A) without or (B) with subjects’ age (AgeB) and smoking status (SmoB) as covariates. Individual SNPs and/or covariates are ordered according to their weighted average ranks as revealed by the resampling with replacement procedure inherent to the RF methodology following 10,000 runs. Color of each rectangle indicate the total number of bootstrap samples in which a given rank in the VIMP-based ranking was recorded for respective predictor. The more yellowish the color, the higher the total number of samples with respective rank. (Numerical data are provided in S3 Table).

When subjects’ age and smoking status were added to analysis, obtained VIMP-based ranking (Table 8) indicated rs5030789 (RB) and rs1801321 (RD) again as the two best predictors with the strongest association with the BrC/control status, with rs1801320 (RC), rs861530 (X3), and rs2928140 (RG) running up on consecutive three places, with some minor shuffling compared to their ranks obtained in analysis not involving covariates, though. Again, VIMPs of rs5030789 (RB), rs1801321 (RD), rs1801320 (RC) were statistically significant (p<0.05 for all three SNPs), while those of rs2928140 (RG) and rs3212057 (X4) remained close to the edge of statistical significance (p = 0.060 for RG and p = 0.072 for X4). Fig 3B shows the results of bootstrapping the RF-based ranking procedure. Again, the rs5030789 (RB) SNP was ranked 1st in the vast majority of all bootstrapped RF models (98.5% of all 10,000 RFs; Part B in S4 Table) and based on derived weighted average ranks it was confirmed as the best predictor with the strongest association with the BrC/control status. Predictors which were found among the top based on their observed VIMP values (such as rs1801321 (RD), rs2928140 (RG), rs1801320 (RC), and rs861530 (X3)) were again among the top predictors also based on their weighted average ranks.

Analysis of epistatic interactions

The possible effects of epistatic interactions among SNPs on BrC/control status were analyzed by the pRF and MB-MDR strategies.

In the case of 2-way interactions, the MB-MDR strategy revealed the rs2619679/rs2928140 (-4719A/T / 2972C/G; RA/RG) Rad51 SNP interaction as the one with the strongest association with BrC (Table 9A). The -4719AA/2972CG and -4719AT/2972CC two-locus genotypes were identified as the “high-risk” genotypes with both of them being more frequent among BrC cases compared to controls (p = 0.071 and p = 0.015, respectively). Analyzing this non-linear interaction effect on BrC as a whole, carriers of these “high-risk” two-locus genotypes were found to be significantly more frequent among cases (cases vs. controls: 10.9% vs. 1.2%), thus indicating an increased risk of BrC among these subjects (OR = 11.3; 95% CI: 2.5–49.5; p<0.005). Moreover, this interaction effect on BrC risk remained statistically significant following the correction for multiple testing using the 10,000 random permutations test (p = 0.0001). Interestingly, the pRF strategy confirmed the results obtained by MB-MDR also identifying the rs2619679/rs2928140 (-4719A/T / 2972C/G; RA/RG) Rad51 SNP interaction as being in the strongest association with the BrC/control status (Table 9A). In this strategy, omitting such SNP interaction in RF-based classification model increased the BrC/control classification differential error (ΔE) by 1.81%.

Tab. 9. Epistatic SNP interactions with the strongest association with BrC.: Results of MB-MDR and pRF analysis.
Epistatic SNP interactions with the strongest association with BrC.: Results of MB-MDR and pRF analysis.

In the case of 3-way interactions, the MB-MDR strategy revealed the rs2619679/rs5030789/rs2928140 (-4719A/T / -4601A/G / 2972C/G; RA/RB/RG) Rad51 SNP triplet as the one with the strongest association with BrC (Table 9B). The -4719AA/-4601AA/2972CG and -4719AT/-4601GA/2972CC three-locus genotypes were found to be more frequent among BrC cases compared to controls (p = 0.087 and p = 0.013, respectively) and thus identified as the “high-risk” combination of genotypes. Carriers of these “high-risk” triplets were found to be significantly more frequent among BrC cases (cases vs. controls 8.6% vs. 1.1%) suggesting an increased risk of BrC among these subjects (OR = 8.4; 95% CI: 1.8–38.6; p<0.005) which remained significant following the correction for multiple hypotheses testing (p<<0.0001). The significant association between the rs2619679/rs5030789/rs2928140 (-4719A/T / -4601A/G / 2972C/G; RA/RB/RG) Rad51 SNP interaction and the BrC/control status was confirmed also in pRF strategy, according to which involvement of this interaction effect in RF-based classification models improves their differential error (ΔE) by 3.57% (Table 9B). Detailed distribution of subjects with respect to genotypes carried at SNP loci involved in analyzed 2-way or 3-way multi-locus genotypes is presented in S5 Table.


In the hereby presented study we examined the role of genetic variability of two proteins belonging to the HR DSB DNA repair pathway—Rad51 and Xrcc3—as a risk factor for breast cancer in Polish population. In total we investigated 14 common single nucleotide polymorphisms in the genes encoding the above mentioned enzymes, seven SNPs per each protein (rs2619679, rs5030789, rs1801320, rs1801321, rs2619680, rs2619681, and rs2928140 in Rad51; rs1799794, rs45603942, rs861530, rs3212057, rs1799796, rs861539, and rs28903081 in Xrcc3).

In the case of Rad51, no associations with BrC were found for rs2619679, rs1801321, rs2619680, rs2619681, and rs2928140 under any of genetic models assumed in this study. Contrary to this, our study provides some very interesting outcomes concerning the rare rs5030789 (-4601A/G; RB) -4601AA genotype which can be associated with reduced BrC risk under recessive genetic model (OR = 0.5; p<0.05). Outcomes of the analysis under additive genetic model and direct analysis of allele frequencies seem to provide further support in favor of this conclusion suggesting a 30% BrC risk reduction being associated with the variant -4601A allele, yet just beyond the edge of statistical significance (p = 0.06 and p = 0.07, respectively). In our previous study [46] we have found the rare rs5030789 -4601AA genotype to confer some protective effect against head and neck cancer (HNC) among men, thus the present study is yet another report suggesting protective effect of this SNP against cancer. It has to be, however, stated that rs5030789 has not yet been studied in relation to cancer risk by any other group and thus no other reports on its involvement in cancer risk modulation are available. The protective anti-cancer effect of this SNP needs to be thus treated with caution and verified in a larger case-control study. Either way, the rs5030789 in Rad51 seems to be a plausible cancer risk-reducing SNP.

Recent huge meta-analyses involving several tens of thousands of subjects provided solid evidence that the variant rs1801320 (135G/C; RC) 135C Rad51 allele localized in the 5’UTR the gene increases the overall risk of cancer [10,24]. The same effect was suggested also in the case of BrC, with the odds ratio of BrC under the recessive genetic model being estimated at 1.7 [24] and 3.3 [10], respectively. The hereby presented outcomes of our study are in line with those cited above, as we have also found the 135CC genotype to be more frequent among BrC cases as compared to healthy controls (7.1% vs. 0.8%, respectively), resulting in around 10-times higher odds of BrC among rare homozygotes. Of note might be the fact, that the frequencies and odds reported by our study differ somewhat from the ones reported by earlier study conducted in Polish population, in which this variant genotype was found to be present in as much as almost 70% of BrC cases and 20% of healthy controls [55]. Discrepancies among BrC cases may at least partially be explained by the fact that the study [55] was focused on triple-negative BrC cases only, while we applied no filtering of the BrC cases based on their estrogen, progesterone and HER-2 receptor status. Large difference in frequencies among controls is, however, strange and difficult to be explained. Nevertheless, the study [55] has also provided quite high values of odds ratio for triple-negative BrC being associated with the rare 135CC Rad51 genotype (OR = 6.0), analogically to our study.

Of note is the fact that our study seems to provide no evidence in favor of any risk-modifying effect of rs1801321 (172G/T; RD) Rad51 SNP, yet another SNP in Rad51 frequently investigated in relation to cancer risk. This SNP was shown to be located in the P300/CBP transcription factor binding site leading to increased activity of Rad51 promoter and increased capacity of DSB DNA repair pathway [22,5658]. Even though the 172T allele-containing genotypes were found to be associated with reduced risk of cancer in general [24], their effect on BrC risk could not be unambiguously specified so far. In Polish population, there are only two relatively small studies available so far, delivering results not allowing to draw any conclusion, as well [25,27]. Our study does not provide any evidence useful in clarifying this conundrum, thus the role of 172T Rad51 SNP in BrC development remains unresolved.

Concerning the SNPs in Xrcc3 investigated in this study, only rs3212057 (10343G/A, X4) localized in exon 6 of the gene resulting in 94Arg/His mismatch mutation can plausibly be associated with BrC. Variant 10343A allele was found among the control subjects only, which suggests its protective effect against BrC. Unfortunately, observed level of significance of such effect was close to marginal 0.05 (p = 0.0453) and we were unable to find any supportive outcome in terms of distribution of genotypic frequencies compared under any of the genetic models assumed. Therefore, this outcome needs to be interpreted cautiously mainly due to the fact that this SNP turned out to be quite rare and our study surely lacked the statistical power needed to provide additional reliable outcomes.

Concerning the SNPs in Xrcc3 investigated in this study, only rs3212057 (10343G/A, 94Arg/His, X4) in exon 6 of the gene can plausibly be associated with BrC. Variant 10343A allele was found among the control subjects only, which suggests its protective effect against BrC. Unfortunately, observed level of significance of such effect was close to marginal 0.05 (p = 0.0453) and we were unable to find any supportive outcome in terms of distribution of genotypic frequencies compared under any of the genetic models assumed. Therefore, this outcome needs to be interpreted cautiously mainly due to the fact that this SNP turned out to be quite rare and our study surely lacked the statistical power needed to provide additional reliable outcomes. Nevertheless, it should be emphasized that the 94Arg/His SNP finds itself in the Xrcc3 gene segment encoding amino acids 63 to 346 of the protein, a region directly involved in Xrcc3/Rad51C heterodimer formation [59]. Owing to crucial role played by this heterodimer in several steps of HR (binding of DNA, resolution of Holiday junctions) [12], a polymorphism affecting its formation may thus influence the capacity of the whole HR DNA repair machinery and play its role in cancer development. Evidence in favor of this hypothesis is, still, rather scarce, as in addition to our previous study in which the Xrcc3 94His allele was found to be associated with increased risk of HNC [46], only four other studies investigating the association between 94Arg/His SNP and risk of various cancers exist [6063]. These studies were, however, conducted in Taiwanese population completely lacking any variability at this locus, so they failed to provide any valuable information on plausible role of rs3212057 in cancer risk modulation. Taken all together, the hereby suggested cancer type- and population-specific association between 94Arg/His SNP in Xrcc3 and cancer risk needs to be further verified.

We failed to find any statistically significant association between rs861539 (241Thr/Met) in exon 6 of Xrcc3 (the most frequently investigated Xrcc3 SNP) and BrC risk. A huge meta-analysis has shown that this SNP provides slight but statistically significant BrC risk increase [64], mainly due to altered protein function, increased genetic instability and DNA DSB accumulation [23]. Recently performed studies on Polish population, however, failed to confirm such effect on either unselected or triple-negative BrC [32,33]. It has to be, however stated, that these two Polish studies, alike the hereby presented one, were relatively small in their sizes (up to 200 BrC subjects only), thus the effect revealed by a meta-analysis involving almost 10,000 BrC cases may have not simply been provable.

For those SNPs identified in single-site analyses as significantly associated with BrC risk, we further analyzed the effect of all possible combinations of respective high-risk genotypes on BrC. Based on obtained outcomes, it is obvious that the relationship between BrC risk and the number of possessed high-risk genotypes is not linear nor additive. One may, however notice that the combinations containing Rad51 -4601GG or -4601GA (RB) seem to confer a higher BrC risk compared to the risk associated with the RB SNP itself, while, on the other hand, the combinations containing Rad51 135CC (RC) SNP tend not to present such distinct one-directional change of BrC risk when compared to RC-only-associated BrC risk. Interpretation of how the resultant disease risk changes due to such combinations is even more difficult due to our inability to exactly count the BrC risk associated with the Xrcc3 10343G/A (X4) SNP (due to zero BrC cases carrying the low-risk genotype). Even though it seems unable to draw any further conclusions on the influence of combinations of Rad51 and Xrcc3 high-risk genotypes on BrC risk, it is worth mentioning that in spite of the hereby presented outcomes being obtained based on relatively small number of subjects (see Table 7) in an analysis adjusted to confounders, the outcomes related to genotype combinations are characterized with considerably higher levels of significance compared to p-levels obtained in respective single-site analyses. It thus clearly shows that analyzing SNP combinations instead of individual SNPs may lead to increased statistical power and that it may indeed be the way of how to move forward in search for links between genetic variability and complex diseases, especially in the case of SNPs within closely related proteins from a given pathway.

In haplotype analysis, none of the 6 haplotypes reconstructed within the rs1801321-rs2619680- rs2619681 (i.e. RD-RE-RF) LD block in Rad51 were found to be associated with the BrC risk, which seems to be in line with the results of above discussed single-site analyses of these three SNPs. Contrary to this, outcomes of haplotype analysis for rs5030789 (RB) and rs1801320 (RC), which were found to be associated with BrC risk in single-site and SNP combinations analyses, however suggest that these SNPs are not part of any LD block and most probably recombine during meiotic chromatid segregation and modulate the BrC risk independently of each other. In the case of Xrcc3, rs3212057 (X4), significantly associated with BrC risk in single-site analysis, was not the part of any of two LD blocks identified within the Xrcc3 gene. Instead, the 4541A/9685A haplotype within the LD block spanning from rs1799794 (X1) to rs861530 (X3) was found to confer a significantly lower BrC risk. Such an interesting observation of two SNPs not associated with BrC risk in single-site analyses revealing significant risk modulating effect in haplotype analysis may simply indicate that the possession of certain variant in one SNP locus might be not enough to impose the risk-modulating effect, and that certain “configuration” of two or several other SNPs is required for the effect to take place. Nevertheless, the haplotype analyses in relation to cancer risk are still very scarce, rendering thus discussion of such outcomes quite challenging. The very few studies analyzing possible cancer risk-modifying effects of Rad51 and/or Xrcc3 haplotypes are, unfortunately, related to different type of cancers [46,65].

In addition to above described analyses, we investigated the effect of nonlinear SNP-SNP (i.e. epistatic) interactions between selected SNPs on the risk of BrC as well. Pure nonlinear epistatic interactions are believed to affect functionality of ternary and quaternary structures involved in specific biological processes without any directly observable main effects of interacting SNPs [53]. Since Xrcc3 and Rad51C were shown to form a heterodimer which directly interacts with Rad51 facilitating the functionality of HR DNA DSB repair machinery, we hypothesized that epistatic interactions between these three proteins might at least partially influence the system’s DNA repair capacity and affect the resultant BrC risk. However, due to the fact that conventional statistical methods (linear regression, logistic regressions, chi-square tests, etc.) turned out to be ineffective and not able to deal with some very specific challenges faced in such kind of analyses (such as the so-called curse of dimensionality, large computational burden, etc.), machine learning techniques have recently became increasingly implemented in order to uncover associations between complex diseases (including cancer) and such otherwise “hidden” interactions. These techniques are becoming a hallmark of the so-called post-GWAS era, in which it is known that risk profiles generated by common low and/or moderate susceptibility loci put together in a simple additive model provide only limited usefulness with respect to complex diseases risk prediction [66,67].

To elucidate plausible effects of nonlinear epistatic interactions between Xrcc3, Rad51 and Rad51C SNPs on BrC risk, two machine learning techniques, presenting two different approaches were used: MB-MDR which is based on a simple model (logistic regression), and pRF built upon a model-free CART analysis. Quite surprisingly, despite different approaches, both these methods pointed out the same interaction as the one with the strongest effect on BrC risk, irrespectively of whether 2-way or 3-way interactions were considered. In the case of 2-way interactions, a simple yet interesting model was proposed by MB-MDR, according to which the highest BrC risk is dependent on inheriting just one heterozygous genotype from two different loci—either the rs2619679 (RA) -4719AT or the rs2928140 (RG) 2972CG, but not both. We cannot currently provide any biologically plausible interpretation of this outcome, as details of molecular interaction among Rad51-family proteins and/or DNA are still unknown. We can only speculate that these SNPs may be localized in gene regions crucial with respect to biological/functional properties of Rad51. Of note is, however, the fact that this interactional effect on BrC risk has been identified as the pure epistatic effect, with no main effects observable for neither rs2619679 (RA) or rs2928140 (RG) in single-site analyses, nor they have been suggested as being involved in haplotypes affecting the BrC risk. According to pRF outcomes, involving this interaction in a classification model results in increased classification “correctness” by some 2% (Table 9), which doesn’t seem to be a considerably high value, but it is quite astonishing when we realize that it is a result of allowing the model for a combination of just 2 out of some 10 million SNPs found in human genome. It is just as striking that subjects belonging to high-risk group identified within this model are conferred with quite considerable BrC risk increase (OR = 11.3), which furthermore turned out to be highly statistically significant (p10000 = 0.0001, corrected for multiple hypotheses testing). Of course, this model cannot be yet considered as clinically useful, as more studies aimed to verify and validate its clinical relevancy in larger experimental setup involving also higher-order interaction models are needed. The importance of the order of such interaction is argued for by the fact that involving 3-way epistatic interactions among SNPs in a classification model almost doubled the gain in classification correctness in this study. Here, both methods indicated the rs2619679/rs5030789/rs2928140 (i.e. RA/RB/RG) interaction as the one with the strongest association to BrC. Adding the rs5030789 (RB) -4601G/A SNP into previously identified RA/RG best 2-way interaction provides basically more detailed specification of the high-risk group, showing that high risk of BrC in subjects carrying either rs2619679 (RA) -4719AT or the rs2928140 (RG) 2972CG genotypes is dependent on co-occurrence of either rs5030789 (RB) -4601GA or the rs5030789 (RB) -4601AA genotypes, respectively (Table 9). It all somehow suggests that it takes two more rare RB alleles (i.e. -4601AA) for the -4719AA/2972CG carriers to exert the BrC risk-increasing effect, while only one such RB allele (i.e. -4601GA) is sufficient for BrC risk increase among -4719AT/2972CC carriers. Even though the high-risk group was narrowed based on additional third locus, it still retained its relatively high increase of BrC risk (OR = 8.5) as well as its high level of statistical significance (p10000 = 0.0001 corrected for multiple hypotheses testing). The correctness of classification model was also improved by almost 3.6%. Nevertheless, the biological/functional background of such effects are unknown.

To sum all outcomes up, Table 10 briefly highlights all outcomes obtained in the study. Of note here is especially the rs5030789 (RB) SNP, for which evidence suggesting its involvement in BrC risk modulation was provided by four out of five different analytical methods (single-site analysis, VIMP-based ranking, MB-MDR and pRF). Moreover, it was the only SNP found to be associated with BrC in both single-site and epistatic interaction analysis, although the predictions of this SNP’s effect on BrC risk seems to be contradictory. It might be a bit puzzling to understand how both protective and detrimental effects on disease risk may be exerted by a single SNP, nevertheless, it has to be kept in mind, that, by definition, epistatic effects based on interaction with other SNPs may lead to completely different effects compared to those observed for individual interacting SNPs. Considering the fact that this SNP is placed within the crucial player of the whole HR DSB DNA repair machinery, it seems as a strong rationale for rs5030789 (RB) being indeed of importance with respect to BrC development, although complete understanding of its role requires further studies.

Tab. 10. Summary of inferred effects of individual SNPs on BrC risk.
Summary of inferred effects of individual SNPs on BrC risk.

The rare rs1801320 (RC) genotype, on the other hand, seems to exert detrimental effect on BrC risk and the conclusion in based on outcomes of single-site analysis and high place in the VIMP-based ranking (3rd or 4th place). Such combination of outcome suggests rather considerable main effect of this genotype on BrC risk, without any additional epistatic effects, as simple VIMP-based rankings cannot distinguish between main and interactional effects [68]. Importance of this SNP with respect to BrC risk seems to be, however, in line with outcomes reported previously by meta-analyses, in which rs1801320 135C allele was found to be associated with significant cancer (including BrC) risk increase [10,24].

Outcomes for rs3212057 (X4) also suggest that heterozygous 10343GA genotype in this locus may protect against BrC, as some indications of such effect were found in single-site analysis and further supported by statistically significant yet rather low (11th and 12th place) place in VIMP-based ranking. This conclusion should be, however, interpreted with caution as the amount of currently available evidence on involvement of rs3212057 in carcinogenesis is still limited.

Outcomes for rs1799794 (X1) and rs861530 (X3) (statistically significant protective effect of X1/X3 haplotype), rs1801321 (RD) (2nd place in VIMP-based rankings with statistically significant VIMP values) as well as rs2619679 (RA) and rs2928140 (RG) (both involved in risk-increasing three-locus genotype revealed by both MB-MDR and pRF) these are relatively novel outcomes definitely requiring further verification (Table 10).

It is worth mentioning that the analysis based on machine learning techniques pointed out an interaction between SNPs localized within the only one gene. This conclusion is further supported by the VIMP-based ranking, where Rad51 SNPs took four out of five best ranks. Although this underlines the importance of Rad51 SNPs with respect to BrC, we cannot say whether it would be so also in the case of higher–order interactions.

Among minor limitations our study admittedly suffers from, one has to mention the lack of proper case-control matching and limited study size. While the lack of proper matching was partially solved by considering certain subjects’ characteristics as confounders in each analysis enabling to do so, the limited study size was balanced out by employing novel machine learning approach specifically designed to overcome such limitation. Nevertheless, the low study size still prevented us from being able to verify our outcomes in a validation group, an approach often used as a golden standard with the hereby used techniques. It is also possible that some other machine learning techniques could have led to slightly different outcomes, nevertheless, the goal of this study was to investigate the role of genetic variability of Rad51 family members in BrC development and compare outcomes obtained by conventional analytical methods to those obtained by the most popular novel machine learning techniques used in the field, rather than provide comprehensive comparison of all possible approaches. Last but not least, data on some of well-known BrC risk factors interacting with HR DSB DNA repair pathway (such as BRCA1 and BRCA2 mutations, alcohol consumption) would possibly increase the significance of outcomes and help better understand the above discussed relationships, but such data were unavailable. Either way, we still argue that our study provides some valuable novel outcomes which may successfully provide clues for future fruitful research.

To sum all up, our study provides evidence that the genetic variability of Xrcc3 and Rad51 may be of relevance with respect to BrC risk modulation. Especially the rs5030789 -4601G/A Rad51 SNP seems to be of importance in this regard, as it was found to independently predict the disease risk as well as to co-participate (together with specific rs2619679 -4719A/T and rs2928140 2972C/G genotypes) in a BrC risk-modulating epistatic interactions, suggesting its possible complex role in BrC development. Important roles in BrC risk modulation were suggested also for rs1801320 135CC Rad51 genotype and rs3212057 10343A (94His) Xrcc3 allele.

Supporting information

S1 Fig [d1]
Multiple sequence alignment of Rad51, Xrcc3 and Rad51C amino acid sequences across 16 different species.

S1 File [docx]
Thorough description of algorithms used in analyses based on machine learning techniques.

S1 Table [xls]

S2 Table [xls]
Results of the systematic examination of the impact of crucial parameters affecting the resultant RF ability to accurately predict the BrC/control status.

S3 Table [xls]
Basic characteristics of RFs used in analyses aimed to obtain the VIMP-based ranking of predictors and in RF-based analysis of epistatic interactions.

S4 Table [xls]
Bootstrap estimates of distribution of the VIMP-based ranks of analyzed predictors.

S5 Table [xls]
Distribution of BrC and control subjects with respect to genotypes carried at SNP loci involved in analyzed 2-way or 3-way multi-locus genotype.


1. Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA. Cancer J. Clin. 2015;65(2):87–108 doi: 10.3322/caac.21262 25651787

2. Sekhar D, Pooja S, Kumar S, Rajender S. RAD51 135G>C substitution increases breast cancer risk in an ethnic-specific manner: a meta-analysis on 21,236 cases and 19,407 controls. Sci. Rep. 2015;5(1):11588

3. Stratton MR, Rahman N. The emerging landscape of breast cancer susceptibility. Nat. Genet. 2008;40(1):17–22 doi: 10.1038/ng.2007.53 18163131

4. Turnbull C, Rahman N. Genetic predisposition to breast cancer: past, present, and future. Annu. Rev. Genomics Hum. Genet. 2008;9:321–45 doi: 10.1146/annurev.genom.9.081307.164339 18544032

5. Levy-Lahad E. Fanconi anemia and breast cancer susceptibility meet again. Nat.Genet. 2010;42(5):368–9 doi: 10.1038/ng0510-368 20428093

6. Rahman N, Seal S, Thompson D, Kelly P, Renwick A, Elliott A, et al. PALB2, which encodes a BRCA2-interacting protein, is a breast cancer susceptibility gene. Nat. Genet. 2007;39(2):165–7 doi: 10.1038/ng1959 17200668

7. Seal S, Thompson D, Renwick A, Elliott A, Kelly P, Barfoot R, et al. Truncating mutations in the Fanconi anemia J gene BRIP1 are low-penetrance breast cancer susceptibility alleles. Nat. Genet. 2006;38(11):1239–41 doi: 10.1038/ng1902 17033622

8. Ding S-L, Yu J-C, Chen S-T, Hsu G-C, Kuo S-J, Lin YH, et al. Genetic variants of BLM interact with RAD51 to increase breast cancer susceptibility. Carcinogenesis 2009;30(1):43–9 doi: 10.1093/carcin/bgn233 18974064

9. Sun H, Bai J, Chen F, Jin Y, Yu Y, Jin L, et al. RAD51 G135C polymorphism is associated with breast cancer susceptibility: a meta-analysis involving 22,399 subjects. Breast Cancer Res. Treat. 2011;125(1):157–61 doi: 10.1007/s10549-010-0922-z 20454923

10. Zhang B-B, Wang D-G, Xuan C, Sun G-L, Deng K-F. Genetic 135G/C polymorphism of RAD51 gene and risk of cancer: a meta-analysis of 28,956 cases and 28,372 controls. Fam. Cancer 2014;13(4):515–26 doi: 10.1007/s10689-014-9729-0 24859942

11. Ripperger T, Gadzicki D, Meindl A, Schlegelberger B. Breast cancer susceptibility: current knowledge and implications for genetic counselling. Eur. J. Hum. Genet. 2009;17(6):722–31 doi: 10.1038/ejhg.2008.212 19092773

12. Masson J, Tarsounas MC, Stasiak AZ, Stasiak A, Shah R, Michael J, et al. Identification and purification of two distinct complexes containing the five RAD51 paralogs. Genes Dev. 2001;15(24):3296–307 doi: 10.1101/gad.947001 11751635

13. Thacker J. The RAD51 gene family, genetic instability and cancer. Cancer Lett. 2005;219(2):125–35 doi: 10.1016/j.canlet.2004.08.018 15723711

14. Badie S, Liao C, Thanasoula M, Barber P, Hill MA, Tarsounas M. RAD51C facilitates checkpoint signaling by promoting CHK2 phosphorylation. J.Cell Biol. 2009;185(4):587–600 doi: 10.1083/jcb.200811079 19451272

15. Liu Y, Masson JY, Shah R, O’Regan P, West SC. RAD51C is required for Holliday junction processing in mammalian cells. Science. 2004;303(5655):243–6 doi: 10.1126/science.1093037 14716019

16. Somyajit K, Subramanya S, Nagaraju G. Distinct roles of FANCO/RAD51C protein in DNA damage signaling and repair: implications for Fanconi anemia and breast cancer susceptibility. J. Biol. Chem. 2012;287(5):3366–80 doi: 10.1074/jbc.M111.311241 22167183

17. Saleh-Gohari N, Bryant HE, Schultz N, Parker KM, Cassel TN, Helleday T. Spontaneous homologous recombination is induced by collapsed replication forks that are caused by endogenous DNA single-strand breaks. Mol. Cell. Biol. 2005;25(16):7158–69 doi: 10.1128/MCB.25.16.7158-7169.2005 16055725

18. Sage JM, Gildemeister OS, Knight KL. Discovery of a novel function for human Rad51: maintenance of the mitochondrial genome. J. Biol. Chem. 2010;285(25):18984–90 doi: 10.1074/jbc.M109.099846 20413593

19. Meindl A, Hellebrand H, Wiek C, Erven V, Wappenschmidt B, Niederacher D, et al. Germline mutations in breast and ovarian cancer pedigrees establish RAD51C as a human cancer susceptibility gene. Nat. Genet. 2010;42(5):410–4 doi: 10.1038/ng.569 20400964

20. Vuorela M, Pylkas K, Hartikainen JM, Sundfeldt K, Lindblom A, von Wachenfeldt WA, et al. Further evidence for the contribution of the RAD51C gene in hereditary breast and ovarian cancer susceptibility. Breast Cancer Res.Treat. 2011;130(3):1003–10. doi: 10.1007/s10549-011-1677-x 21750962

21. De Leeneer K, Van Bockstal M, De Brouwer S, Swietek N, Schietecatte P, Sabbaghian N, et al. Evaluation of RAD51C as cancer susceptibility gene in a large breast-ovarian cancer patient population referred for genetic testing. Breast Cancer Res. Treat. 2012;133(1):393–8 doi: 10.1007/s10549-012-1998-4 22370629

22. Hasselbach L, Haase S, Fischer D, Kolberg HC, Sturzbecher HW. Characterisation of the promoter region of the human DNA-repair gene Rad51. Eur. J Gynaecol. Oncol. 2005;26(6):589–98. 16398215

23. Nowacka-Zawisza M, Wiśnik E, Wasilewski A, Skowrońska M, Forma E, Bryś M, et al. Polymorphisms of Homologous Recombination RAD51, RAD51B, XRCC2, and XRCC3 Genes and the Risk of Prostate Cancer. Anal. Cell. Pathol. 2015;2015:1–9

24. Zhao M, Chen P, Dong Y, Zhu X, Zhang X. Relationship between Rad51 G135C and G172T Variants and the Susceptibility to Cancer: A Meta-Analysis Involving 54 Case-Control Studies. He B, editor. PLoS One 2014;9(1):e87259 doi: 10.1371/journal.pone.0087259 24475258

25. Michalska MM, Samulak D, Romanowicz H, Smolarz B. Single Nucleotide Polymorphisms (SNPs) of RAD51-G172T and XRCC2-41657C/T Homologous Recombination Repair Genes and the Risk of Triple- Negative Breast Cancer in Polish Women. Pathol. Oncol. Res. 2015;21(4):935–40 doi: 10.1007/s12253-015-9922-y 25743260

26. Al-Zoubi MS, Mazzanti CM, Zavaglia K, Al Hamad M, Armogida I, Lisanti MP, et al. Homozygous T172T and Heterozygous G135C Variants of Homologous Recombination Repairing Protein RAD51 are Related to Sporadic Breast Cancer Susceptibility. Biochem. Genet. 2016;54(1):83–94 doi: 10.1007/s10528-015-9703-z 26650628

27. Sassi A, Popielarski M, Synowiec E, Morawiec Z, Wozniak K. BLM and RAD51 genes polymorphism and susceptibility to breast cancer. Pathol. Oncol. Res. 2013;19(3):451–9 doi: 10.1007/s12253-013-9602-8 23404160

28. Tulbah S, Alabdulkarim H, Alanazi M, Parine NR, Shaik J, Pathan AAK, et al. Polymorphisms in RAD51 and their relation with breast cancer in Saudi females. Onco. Targets. Ther. 2016;9:269–77 doi: 10.2147/OTT.S93343 26834486

29. Chai F, Liang Y, Chen L, Zhang F, Jiang J. Association between XRCC3 Thr241Met Polymorphism and Risk of Breast Cancer: Meta-Analysis of 23 Case-Control Studies. Med. Sci. Monit. 2015;21:3231–40 doi: 10.12659/MSM.894637 26498491

30. Mao C-F, Qian W-Y, Wu J-Z, Sun D-W, Tang J-H. Association between the XRCC3 Thr241Met polymorphism and breast cancer risk: an updated meta-analysis of 36 case-control studies. Asian Pac. J. Cancer Prev. 2014;15(16):6613–8 doi: 10.7314/apjcp.2014.15.16.6613 25169497

31. He X-F, Wei W, Su J, Yang Z-X, Liu Y, Zhang Y, et al. Association between the XRCC3 polymorphisms and breast cancer risk: meta-analysis based on case–control studies. Mol. Biol. Rep. 2012;39(5):5125–34 doi: 10.1007/s11033-011-1308-y 22161248

32. Smolarz B, Makowska M, Samulak D, Michalska MM, Mojs E, Wilczak M, et al. Association between single nucleotide polymorphisms (SNPs) of XRCC2 and XRCC3 homologous recombination repair genes and triple-negative breast cancer in Polish women. Clin. Exp. Med. 2015;15(2):151–7 doi: 10.1007/s10238-014-0284-7 24728564

33. Romanowicz H, Pyziak Ł, Jabłoński F, Bryś M, Forma E, Smolarz B. Analysis of DNA Repair Genes Polymorphisms in Breast Cancer. Pathol. Oncol. Res. 2017;23(1):117–23 doi: 10.1007/s12253-016-0110-5 27571987

34. Zheng Y, Zhang J, Hope K, Niu Q, Huo D, Olopade OI. Screening RAD51C nucleotide alterations in patients with a family history of breast and ovarian cancer. Breast Cancer Res.Treat. 2010;124(3):857–61 doi: 10.1007/s10549-010-1095-5 20697805

35. Akbari MR, Tonin P, Foulkes WD, Ghadirian P, Tischkowitz M, Narod SA. RAD51C germline mutations in breast and ovarian cancer patients. Breast Cancer Res. 2010;12(4):404 doi: 10.1186/bcr2619 20723205

36. Wong MW, Nordfors C, Mossman D, Pecenpetelovska G, Avery-Kiejda KA, Talseth-Palmer B, et al. BRIP1, PALB2, and RAD51C mutation analysis reveals their relative importance as genetic susceptibility factors for breast cancer. Breast Cancer Res.Treat. 2011;127(3):853–9 doi: 10.1007/s10549-011-1443-0 21409391

37. Pang Z, Yao L, Zhang J, Ouyang T, Li J, Wang T, et al. RAD51C germline mutations in Chinese women with familial breast cancer. Breast Cancer Res.Treat. 2011;129(3):1019–20 doi: 10.1007/s10549-011-1574-3 21597919

38. Clague J, Wilhoite G, Adamson A, Bailis A, Weitzel JN, Neuhausen SL. RAD51C germline mutations in breast and ovarian cancer cases from high-risk families. PLoS.One. 2011;6(9):e25632 doi: 10.1371/journal.pone.0025632 21980511

39. Gresner P, Gromadzinska J, Jablonska E, Stepnik M, Zambrano Quispe O, Twardowska E, et al. Single nucleotide polymorphisms in noncoding regions of Rad51C do not change the risk of unselected breast cancer but they modulate the level of oxidative stress and the DNA damage characteristics: a case-control study. Woloschak GE, editor. PLoS One 2014;9(10):e110696 doi: 10.1371/journal.pone.0110696 25343521

40. Moore JH, Gilbert JC, Tsai CT, Chiang FT, Holden T, Barney N, et al. 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–61 doi: 10.1016/j.jtbi.2005.11.036 16457852

41. Pomerleau CS, Pomerleau OF, Snedecor SM, Mehringer AM. Defining a never-smoker: results from the nonsmokers survey. Addict.Behav. 2004;29(6):1149–54 doi: 10.1016/j.addbeh.2004.03.008 15236816

42. Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, et al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29(1):308–11 doi: 10.1093/nar/29.1.308 11125122

43. Vorontsov IE, Kulakovskiy I V., Khimulya G, Nikolaeva DD, Makeev VJ. PERFECTOS-APE: Predicting regulatory functional effect of SNPs by approximate P-value estimation. Bioinforma. 2015—6th Int. Conf. Bioinforma. Model. Methods Algorithms, Proceedings; Part 8th Int. Jt. Conf. Biomed. Eng. Syst. Technol. BIOSTEC 2015 2015;102–8

44. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat. Methods. 2010. p. 248–9 doi: 10.1038/nmeth0410-248 20354512

45. The UniProt Consortium. UniProt: a worldwide hub of protein knowledge. Nucleic Acids Res. 2019;47(D1):D506–15 doi: 10.1093/nar/gky1049 30395287

46. Gresner P, Gromadzinska J, Polanska K, Twardowska E, Jurewicz J, Wasowicz W. Genetic variability of Xrcc3 and Rad51 modulates the risk of head and neck cancer. Gene. 2012;504(2):166–74 doi: 10.1016/j.gene.2012.05.030 22613844

47. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–5 doi: 10.1093/bioinformatics/bth457 15297300

48. Lewontin RC. The detection of linkage disequilibrium in molecular sequence data. Genetics. 1995;140(1):377–88 7635301

49. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, et al. The structure of haplotype blocks in the human genome. Science. 2002;296(5576):2225–9 doi: 10.1126/science.1069424 12029063

50. Breiman L. Random Forests. Mach. Learn. 2001;45:5–32

51. Ishwaran H, Kogalur U. Fast Unified Random Forests for Survival, Regression, and Classification (RF-SRC). R package version 2.9.2. 2019

52. Hankin RKS. Package “permutations” R package version 1.0.5 2017

53. Li J, Malley JD, Andrew AS, Karagas MR, Moore JH, Hirschhorn J, et al. Detecting gene-gene interactions using a permutation-based random forest method. BioData Mining; 2016;9(1):14

54. Calle ML, Urrea V, Malats N, Van Steen K. mbmdr: an R package for exploring gene–gene interactions associated with binary or quantitative traits. Bioinformatics 2010;26(17):2198–9. doi: 10.1093/bioinformatics/btq352 20595460

55. Smolarz B, Zadrożny M, Duda-Szymańska J, Makowska M, Samulak D, Michalska MM, et al. RAD51 genotype and triple-negative breast cancer (TNBC) risk in Polish women. Pol. J. Pathol. 2013;64(1):39–43 doi: 10.5114/pjp.2013.34602 23625599

56. Lu J, Wang L-E, Xiong P, Sturgis EM, Spitz MR, Wei Q. 172G>T variant in the 5’ untranslated region of DNA repair gene RAD51 reduces risk of squamous cell carcinoma of the head and neck and interacts with a P53 codon 72 variant. Carcinogenesis 2007;28(5):988–94 doi: 10.1093/carcin/bgl225 17118968

57. Flygare J, Falt S, Ottervald J, Castro J, Dackland AL, Hellgren D, et al. Effects of HsRad51 overexpression on cell proliferation, cell cycle progression, and apoptosis. Exp.Cell Res. 2001;268(1):61–9 doi: 10.1006/excr.2001.5265 11461118

58. Yoo S, McKee BD. Overexpression of Drosophila Rad51 protein (DmRad51) disrupts cell cycle progression and leads to apoptosis. Chromosoma. 2004;113(2):92–101 doi: 10.1007/s00412-004-0300-x 15257466

59. Kurumizaka H, Enomoto R, Nakada M, Eda K, Yokoyama S, Shibata T. Region and amino acid residues required for Rad51C binding in the human Xrcc3 protein. Nucleic Acids Res. 2003;31(14):4041–50 doi: 10.1093/nar/gkg442 12853621

60. Liu J-C, Tsai C-W, Hsu C-M, Chang W-S, Li C-Y, Liu S-P, et al. Contribution of double strand break repair gene XRCC3 genotypes to nasopharyngeal carcinoma risk in Taiwan. Chin. J. Physiol. 2015;58(1):64–71 doi: 10.4077/CJP.2015.BAD279 25687493

61. Chen H-J, Chang W-S, Hsia T-C, Miao C-E, Chen W-C, Liang S-J, et al. Contribution of Genotype of DNA Double-strand Break Repair Gene XRCC3, Gender, and Smoking Behavior to Lung Cancer Risk in Taiwan. Anticancer Res. 2015;35(7):3893–9 26124335

62. Chang W-S, Tsai C-W, Wang J-Y, Ying T-H, Hsiao T-S, Chuang C-L, et al. Contribution of X-Ray Repair Complementing Defective Repair in Chinese Hamster Cells 3 (XRCC3) Genotype to Leiomyoma Risk. Anticancer Res. 2015;35(9):4691–6 26254358

63. Su C-H, Chang W-S, Hu P-S, Hsiao C-L, Ji H-X, Liao C-H, et al. Contribution of DNA Double-strand Break Repair Gene XRCC3 Genotypes to Triple-negative Breast Cancer Risk. Cancer Genomics Proteomics 2015;12(6):359–67 26543082

64. He X-F, Wei W, Li J-L, Shen X-L, Ding D, Wang S-L, et al. Association between the XRCC3 T241M polymorphism and risk of cancer: Evidence from 157 case–control studies. Gene 2013;523(1):10–9 doi: 10.1016/j.gene.2013.03.071 23562721

65. Rollinson S, Smith AG, Allan JM, Adamson PJ, Scott K, Skibola CF, et al. RAD51 homologous recombination repair gene haplotypes and risk of acute myeloid leukaemia. Leuk. Res. 2007;31(2):169–74 doi: 10.1016/j.leukres.2006.05.028 16890287

66. Spitz MR, Amos CI, D’Amelio A, Dong Q, Etzel C, Etzel C. Re: Discriminatory accuracy from single-nucleotide polymorphisms in models to predict breast cancer risk. J. Natl. Cancer Inst. 2009;101(24):1731–2 doi: 10.1093/jnci/djp394 19903803

67. Moore JH, Williams SM. Epistasis and Its Implications for Personal Genetics. Am. J. Hum. Genet. 2009;85(3):309–20 doi: 10.1016/j.ajhg.2009.08.006 19733727

68. Jiang R, Tang W, Wu X, Fu W. A random forest approach to the detection of epistatic interactions in case-control studies. BMC Bioinformatics 2009;10 Suppl 1:S65

69. Kulakovskiy I V., Vorontsov IE, Yevshin IS, Sharipov RN, Fedorova AD, Rumynskiy EI, et al. HOCOMOCO: Towards a complete collection of transcription factor binding models for human and mouse via large-scale ChIP-Seq analysis. Nucleic Acids Res. 2018;46(D1):D252–9 doi: 10.1093/nar/gkx1106 29140464

Článok vyšiel v časopise


2020 Číslo 1
Najčítanejšie tento týždeň
Najčítanejšie v tomto čísle
Zabudnuté heslo

Zadajte e-mailovú adresu, s ktorou ste vytvárali účet. Budú Vám na ňu zasielané informácie k nastaveniu nového hesla.


Nemáte účet?  Registrujte sa