Skip to main content

Significant association between ERCC2 and MTHR polymorphisms and breast cancer susceptibility in Moroccan population: genotype and haplotype analysis in a case-control study



Genetic determinants of breast cancer (BC) remained largely unknown in the majority of Moroccan patients. The purpose of this study was to explore the association of ERCC2 and MTHFR polymorphisms with genetic susceptibility to breast cancer in Moroccan population.


We genotyped ERCC2 polymorphisms (rs1799793 (G934A) and rs13181 (A2251C)) and MTHFR polymorphisms (rs1801133 (C677T) and rs1801131 (A1298C)) using TaqMan SNP Genotyping Assays. Genotypes were compared in 151 BC cases and 156 population-matched controls. Allelic, genotypic and haplotype associations with the risk and clinicopathological features of BC were assessed using logistic regression analyses.


ERCC2-rs1799793-AA genotype was associated with high risk of BC compared to wild type genotype (recessive model: OR: 2.90, 95% CI: 1.34–6.26, p = 0.0069) even after Bonferroni correction (p < 0,0125). MTHFR rs1801133-TT genotype was associated with increased risk of BC (recessive model, OR: 2.49, 95% CI: 1.17–5.29, p = 0.017) but the association turned insignificant after Bonferroni correction. For the rest of SNPs, no statistical associations to BC risk were detected.

Significant association with clinical features was detected for MTHFR-rs1801133-TC genotype with early age at diagnosis and familial BC. Following Bonferroni correction, only association with familial BC remained significant. MTHFR-rs1801131-CC genotype was associated with sporadic BC. ERCC2-rs1799793-AA genotype correlated with ER+ and PR+ breast cancer. ERCC2-rs13181-CA genotype was significantly associated large tumors (T ≥ 3) in BC patients. None of these associations passed Bonferroni correction.

Haplotype analysis showed that ERCC2 A-C haplotype was significantly associated with increased BC risk (OR: 3.71, 95% CI: 1.7–8.12, p = 0.0002 and p = 0.0008 before and after Bonferroni correction, respectively) and positive expression of ER and PR in BC patients. ERCC2 G-C haplotype was correlated with PR negative and larger tumor (T4). We did not find any MTHFR haplotypes associated with BC susceptibility. However, the less common haplotype MTHFR T-C was more frequent in young patients and in familial breast cancer, while MTHFR C-C haplotype was associated with sporadic BC form.


Our findings are a first observation of association between ERCC2 SNPs and breast cancer in Moroccan population. The results suggested that ERCC2 and MTHFR polymorphisms may be reliable for assessing risk and prognosis of BC in Moroccan population.

Peer Review reports


Breast cancer (BC) is by far the most frequently diagnosed malignancy among women worldwide. It is a major public health problem in both developed and developing countries. The incidence of BC is steadily increasing over the years. For 2012, there were 1.7 million estimated new cases of BC [1]. More new cases occurred in less developed (883,000 cases) than more developed countries (794,000 cases) [1]. The growing trend has been reported by the Global Burden of Disease (GBD) study in 2015 covering 32 cancer groups in 195 countries. This study ranked the BC as the most common incident cancer for women (2.4 million cases), and as the leading cause of women cancer death (523,000 deaths) [2]. The BC incidence rates are geographically variable with higher rates occurring in Europe and North America than Africa and Asia. The incidence is constantly growing in the Arab countries while remaining below that recorded in Europe or in America [3].

Breast cancer has become the leading cause of malignancy in Moroccan females. The most recent data in Morocco (country of North-western Africa) have described an increasing incidence rate of breast cancer from 39.0 to 49.5 per 100,000 women between 2008 and 2012 [4]. This rate was relatively higher than in other regional countries, but it remained well below the incidence found in Western countries [5].

Breast cancer is a complex and heterogeneous multifactorial disease which is strongly influenced by environmental, lifestyle and genetics risk factors. The effect of rare highly and moderately penetrant alleles located in predisposition genes such as BRCA1, BRCA2, TP53 and DNA repair genes explains only a small percentage of genetic risk of BC. To date, and through multiple previous genome-wide association studies (GWAS), large scale replication studies and meta-analysis studies, more than 90 breast cancer risk SNPs (single nucleotide polymorphisms) have been identified [6,7,8,9,10,11]. Although, individually, these common variants present relatively small increments in BC risk and a modest effect, taken together they may account for about 15–20% of familial clustering and a substantial proportion of sporadic BC susceptibility [12].

The study reported here concerned the population of Morocco. This country is located in the northwestern corner of the African continent (33°, 35’N latitude and 7°, 39’W longitude), bordered by the Mediterranean Sea to the north, the Atlantic Ocean to the west, Algeria to the east and Mauritania to the south. Morocco is host to a number of human populations that are different in their language, culture and ethnic identity. Indeed, this country very coveted since antiquity, and continues to attract peoples coming from the Mediterranean, the near and the Middle East, as well as from sub-Saharan Africa. The overwhelming majority of Moroccan population is composed of Berbers and Arabs. The Berbers, a people of Euro-Asiatic origin are indigenous residents of Morocco since at least 5000 years ago. They were invaded by many civilizations such as Phoenicians, Carthaginians, Romans, Vandals, Byzantines and Arabs. The Arabs came from the Middle East, namely from the Arabian Peninsula, in the 7th Century and conquered the country during the Islamic expansion in North Africa. Other human groups in Morocco are the Africans, Sub-Sahara Africans, Europeans (commonly descended from Spanish or French ancestry), and Sephardic Jews. All of these populations probably have contributed to the genetic diversity of the current population of Morocco.

The genetic basis of BC remains unknown in the majority of Moroccan patients. Identifying genetic factors associated with this prevalent disease is nevertheless of considerable clinical importance. To date, the few genetic studies that have been conducted in Moroccan population have demonstrated an important but complex contribution of genetic factors in BC pathogenesis as reflected by an increased frequency (27.5–31.6%) of BRCA1/2 mutations detected in familial BC cases [13, 14]. Beside BRCA1/2 pathogenic mutations, relatively few single nucleotide polymorphisms (SNPs) have been studied in Moroccan population. Investigations have yielded a small number of suggestive SNPs associated with varying risks of developing BC [15, 16].

Considered collectively, these observations strongly suggested that other loci may be involved in genetic predisposition to BC in Morocco and clear remaining hereditary BC risk. In the present study, we aimed to investigate the role of four common genetic variants (SNPs) in mediating the disease in Moroccan patients. All of them have been associated with BC in different populations. Their positive correlation with BC have been supported by a meta-analysis of 150 published meta-analysis studies grouping 4474 studies for various types of cancers, 2,452,510 cases and 3,091,626 controls [10].

The selected variants rs1799793 and rs13181 are linked to ERCC2 gene (Excision repair cross-complementation group 2) and rs1801131, rs1801133 polymorphisms to MTHFR gene (methylenetetrahydrofolate reductase). These genes are involved in the etiology of cancer through crucial cellular pathways including, DNA repair path (ERCC2) [17], methylation and DNA synthesis (MTHFR) [18].

Herein, we examined the association between these SNPs and BC as well their contribution in modulating major breast cancer clinicopathological traits in Moroccan patients. Moreover, the effect of haplotypes formed by SNPs localized in the same gene was also examined. It should be mentioned that three selected SNPs were considered for the first time in our study in Moroccan BC cases. The MTHFR-rs1801133 polymorphism has been studied before on a group of 96 Moroccan patients [19]. The notable strength of our study, regarding this SNP, was to analyze a large number of patients from a different geographic area of Morocco and to test haplotype associations of MTHFR SNPs.


Study subjects

A total of 151 pathologically confirmed female breast cancer patients admitted to the Hassan II Regional Oncology Center of Oujda city during 2009–2013 were included in this study. This center covers the entire eastern region of Morocco in term of cancer diagnosis and patients management. The control group consisted of 156 age-matched healthy female with no prior history of any type of cancer, and who were recruited as volunteer blood donors at the Blood Transfusion Center of the same region. All cases and controls were genetically unrelated Moroccans from the same geographical area and were recruited during the same period.

Relevant clinicopathological characteristics recorded for each case were collected by review of patients’ medical files. The recorded information included age at diagnosis, family history of breast cancer, laterality, histology type, Scarff-Bloom-Richardson (SBR) grade, tumor size, lymph node involvement, metastases as well as hormone receptor status including: estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (Her-2).

The following inclusion criteria were used to identify family history (FH) breast cancer cases: single breast cancer diagnosed before the age of 39 years and/or bilaterality; three or more first or second degree relatives with breast cancer in the same side of the family tree; two first degree relatives with breast cancer, with at least one early onset breast cancer case (≤40 years) or male breast cancer case or ovarian cancer case; triple-negative breast cancer diagnosed before 50 years regardless to family history or > 50 years with positive family history of breast cancer; multiple primary cancers in the same individual or in the family. A case was considered sporadic in the absence of the above criteria.

The study protocol was reviewed and approved by the ethics committee of Mohammed VI University Hospital in Marrakech.

Before enrollment in the study and after explaining the procedures, written informed consent for research participation was signed by each participant.

Peripheral venous blood samples were collected from case and control groups into sterile EDTA coated tubes. Genomic DNA extraction was done using a standard salting-out method [20]. The isolated DNA samples were quantitated using the NanoVue Plus™ spectrophotometer (biochrom, Harvard Bioscience Inc. Massachusetts, USA), and stored at − 20 °C until analysis.

Selection of SNPs and TaqMan genotyping

After a review of published literature, we selected four candidate SNPs suggested as significant risk factors for the breast cancer [10], specifically, rs1799793, rs13181 on ERCC2 gene and rs1801133, rs1801131 on MTHFR gene.

Genotyping of the SNPs was performed by allelic discrimination using the TaqMan SNP Genotyping Assays according to the manufacturer’s instructions (Applied Biosystems). Specific primers and FAM/VIC – labeled TaqMan probes were designed and supplied by Applied Biosystems. Briefly, the reaction was performed in a 10 μl final volume containing 5 ng of genomic DNA, 1X TaqMan SNP Genotyping Assay and 1X TaqMan Genotyping master mix (that contains AmpliTaq Gold® DNA Polymerase UP (Ultra Pure), dNTPs without dUTP, and passive internal reference based on proprietary ROX™ dye). All assays were carried out in 24-well plates including positive and negative controls. The PCR conditions were as follows: Initiation at 95 °C for 7 min, followed by 50 cycles of denaturation at 95 °C for 50 s and annealing/extension at 60 °C for 30s. Plates were read on a Thermo Scientific™ PikoReal™ Real-Time PCR System (Thermo Fisher Scientific Oy, Finland), and the alleles were assigned using the PikoReal Software v 2.2.

Statistical analysis

Agreement of genotype frequencies to Hardy-Weinberg expectations was assessed independently among control and case groups for each SNP using the χ2 test analysis with one degree of freedom. Student’s t-test was used to evaluate the differences in mean age at diagnosis between the cases and controls.

The analysis of association between a single variant and breast cancer risk in multiple inheritance models (genotype, dominant, recessive, and additive) was presented in odds ratios (OR) with corresponding 95% confidence intervals (95% CI) using logistic regression. The allelic frequencies of each SNP were also compared between cases and controls. The wide-type genotype was regarded as the reference group. Logistic regression analyses restricted to case group were also performed to compute the odds ratio associating different genotypes with patients’ clinicopathological features.

Data management and statistical analyses were performed using the statistical package SPSS (version 21.0; IBM Corp. Armonk, NY, USA) and SNPStats software ().

Haplotype analysis was restricted to polymorphisms located on the same chromosome: the haplotype rs1801133-rs1801131 (MTHFR), and the haplotype rs1799793-rs13181 (ERCC2). Haplotype frequency distributions were deduced from genotype data and compared between cases and controls using UNPHASED software version 3.1.7 [21] as well as SNPStats program [22]. The most common haplotype was selected as the reference. Odds ratios and 95% CI were calculated to estimate the degree of the association between haplotypes and the risk of breast cancer.

Measure of linkage disequilibrium (LD) between each pair of SNPs, including Lewontin’s standardized disequilibrium coefficient (D’) and the squared correlation coefficient (r2), was computed with Haploview software package version 4.2 [23]. Results were confirmed using UNPHASED and SNPStats programs.

The effect of representative haplotypes on subphenotypes in breast cancer cases was assessed by SNPstats program, and results were presented as odds ratios with 95% CI.

A two-tailed p value less than 0.05 was taken as statistically significant. P values obtained were corrected for multiple testing using Bonferroni correction for the number of tests.


Patient’s clinicopathological characteristics

The clinical characteristics of the breast cancer patients enrolled in this study were summarized in Table 1. The mean age of case and control groups at the time of diagnosis was 43.81 ± 9.74 years, (median = 43 years, range, 26–72 years) and 42.72 ± 9.29 years, (median = 41 years, range, 26–60 years), respectively. No statistically significant difference was observed in term of age between the two groups (p = 0.316).

Table 1 Clinicopathological features of the breast cancer patients included in the study

Out of 151 cases, 39.7% of them had a family history of breast cancer. Almost all patients (97.35%) had unilateral tumor. A large proportion have invasive ductal carcinoma (86.1%), and 71.5% had early stage tumor (Stage T2), while 13.2% had advanced-stage tumor (T3 and T4). Histopathologically, 71.5% of patients had intermediate grade II, and 17.9% presented high grade (grade III). Additionally, 60.9% exhibited lymph node involvement, while metastasis was confirmed in only 12.6% of cases. Regarding hormone receptor expression, most of cases harbored a Her-2 negative (74.8%), ER-positive (70.9%) and PR-positive (58.9%) tumor. Triple negative status was observed in 14.6% of patients.

Associations between SNPs and breast cancer risk

A total of 156 controls and 151 BC subjects were successfully genotyped for the following selected SNPs: ERCC2-rs1799793 (G934A), ERCC2-rs13181 (A2251C), MTHFR-rs1801133 (C677T) and MTHFR-rs1801131 (A1298C).

Genotypes and alleles distributions of the 4 polymorphisms in BC case and control groups are depicted in Table 2. The genotype frequencies of all SNPs were in compliance with Hardy-Weinberg equilibrium in control group (p > 0.05). In patients group, the distribution of ERCC2 rs1799793, ERCC2 rs13181 and MTHFR rs1801133 genotypes did not conform to the HWE (p = 0.002, p = 0.007 and p = 0.013, respectively).

Table 2 Genotype and allele distribution of ERCC2 and MTHFR polymorphisms in breast cancer cases and controls

We investigated the genotypic association between the 4 SNPs and BC risk in five genetic models including co-dominant, dominant, recessive, over-dominant and additive models. All results were age-unadjusted; the age-adjusted model (data not shown) did not diminish the significance of associations.

The AA genotype frequency of ERCC2-rs1799793 polymorphism revealed an association with high risk of breast cancer in both homozygote and recessive models (AA vs. GG, odds ratio (OR): 2.66, 95% confidence interval (CI): 1.20–5.91, p = 0.016; AA vs. GG + GA, OR: 2.90, 95% CI: 1.34–6.26, p = 0.0069), respectively. The association was still significant for homozygote model after Bonferrroni correction (p < 0,0125).

Similarly, the TT genotype of MTHFR-rs1801133 polymorphism was found to be associated with increased breast cancer risk in homozygote (TT vs. CC, OR: 2.39, 95% CI: 1.09–5.22, p = 0.028); and recessive (TT vs. CC + CT, OR: 2.49, 95% CI: 1.17–5.29, p = 0.017) models, but the p values could not withstand the Bonferroni correction. For the two remaining SNPs, no significant association was found between the ERCC2-rs13181 and MTHFR-rs1801131 variants and BC in any hereditary model. In addition, the allelic frequencies of all polymorphisms were similar between BC case and control groups.

Subgroup analysis of BC cases according to age at diagnosis and family history

When BC patients were grouped into two categories with regard to age at the diagnosis (age ≤ 40 and age > 40 years) (Table 3), MTHFR rs1801133 revealed a positive correlation with early age at diagnosis (under 40 years); this polymorphism was found to be a BC risk factor among young patients in 4 genetic models (heterozygote: OR: 2.84, 95% CI: 1.35–6.0, p = 0.006; dominant: OR: 2.56, 95% CI: 1.29–5.09, p = 0.0061; over-dominant: OR: 2.34, 95% CI: 1.18–4.67, p = 0.015 and log-additive: OR: 1.62, 95% CI: 1.03–2.55, p = 0.036). However, the association did not meet statistical thresholds of significance after Bonferroni correction.

Table 3 Association of MTHFR-rs1801133, rs1801131, ERCC2-rs1799793 and rs13181 SNPs with age at diagnosis, family history, ER status, PR status and tumor size of breast cancer patients

Interestingly, in the subgroup analysis according to family history, the MTHFR rs1801133 showed a strong association with familial BC in almost all genetic models (heterozygote, OR: 4.64, 95% CI: 2.18–9.86, p = 6.8.10− 5; homozygote: OR: 4.64, 95% CI: 1.73–12.42, p = 0.002; dominant: OR: 4.64, 95% CI: 2.33–9.21, p < 10− 4; over-dominant: OR: 3.11, 95% CI: 1.55–6.24, p = 0.001 and log-additive: OR: 2.61, 95% CI: 1.60–4.25, p < 10− 4). These associations remained significant after Bonferroni correction (p < .0,005) In contrast, there was a significant association between the CC genotype of MTHFR rs1801131 and sporadic form of BC (recessive model: OR: 0.12, 95% CI: 0.01–0.97, p = 0.012). However, it turned insignificant after Bonferroni adjustment.

Otherwise, no significant association was found between the remaining polymorphisms and age at diagnosis or family history.

Association between SNPs and BC clinicopathologic characteristics

We performed further analysis to investigate a possible relationship between the clinicopathological parameters and the distributions of SNPs genotypes in BC group. Positive results of associations are represented in Table 3.

The AA genotype of ERCC2-rs1799793 polymorphism was significantly frequent in ER-positive patients (homozygote model: OR: 0.21, 95% CI: 0.04–0.96, p = 0.045; recessive model: OR: 0.24, 95% CI: 0.05–1.08, p = 0.03; additive model: OR: 0.54, 95% CI: 0.30–0.96, p = 0.028). In addition, patients with the AA genotype of ERCC2-rs1799793 were prone to be PR-positive (recessive model: OR: 0.33, 95% CI: 0.11–1.05, p = 0.042).

Data indicated a trend to higher frequency of ERCC2-rs13181 CC genotype in ER-positive patients (homozygote model: OR: 0.42, 95% CI: 0.11–1.58; recessive model: OR: 0.45, 95% CI: 0.12–1.63). This increased frequency did not reach a significant level. Furthermore, patients with CA genotype showed a 1.83 fold increased risk of developing PR-negative BC (over-dominant model: OR: 1.83, 95% CI: 0.88–3.82).

Stratified analysis by tumor size feature demonstrated that ERCC2-rs13181 CA genotype carriers displayed an elevated risk of BC with large tumor size (T3-T4) in heterozygote model (OR: 4.09, 95% CI: 1.14–14.7, p = 0.03) and overdominant model (OR: 4.38, 95% CI: 1.28–15.1, p = 0.014).

However, these associations did not retain statistical significance after Bonferroni adjustment.

Finally, there was no evidence of significant correlation between all studied loci and other clinicopathological features including histology type, histology grade, lymph node involvement and metastasis in case subjects.

ERCC2 and MTHFR haplotype associations with BC

In this section, we performed association analysis between the risk of BC and SNP haplotypes of ERCC2 gene (rs1799793 - rs13181) on one hand and MTHFR gene (rs1801133 - rs1801131) on the other hand. Haplotypes were reconstructed from the genotypic data and results of their distribution among BC cases and controls were summarized in Table 4.

Table 4 Association of ERCC2 and MTHFR haplotypes and risk of breast cancer

The pairwise linkage disequilibrium is given for each pair of SNPs. The observed low D’ values (0.17 in cases and 0.45 in controls) and low r2 (0.03 in cases and 0.02 in controls) indicated that the studied ERCC2 SNPs were not at high linkage disequilibrium. Likewise, MTHFR were not found in linkage disequilibrium in both cases and controls (cases: D’ = 0.12, r2 = 0.03; controls: D’ = 0.16, r2 = 0.003).

For ERCC2 SNPs, we found all the four expected haplotypes in both cases and controls; the most popular haplotype was G-A, followed by haplotypes A-A, G-C and A-C (cases: 50, 19.5, 16.9 and 13.6%; controls: 51.5, 23.5, 21.2 and 3.8%, respectively). The haplotype containing the two minor alleles A-C was distributed differently between patients and controls. It was significantly associated with about 3.71 fold increase risk of BC when compared to the wild-type haplotype G-A (OR: 3.71, 95% CI: 1.7–8.12, p = 0.0002). These association was maintained after Bonferroni correction (p < 0,0125). The three other haplotypes were distributed similarly between case and control groups.

Considering MTHFR gene, all the four expected haplotypes appeared in our analysis. The most frequent for both BC cases and controls was C-A (rs1801133 C - rs1801131 A) haplotype (48.4 and 52.5%, respectively). The estimated frequencies for the other haplotype were: T-A (26.2 and 21.9%), C-C (18 and 19.6%) and T-C (7.5 and 6%) in cases and controls, respectively. No difference was observed between case and control groups regarding the distribution of all the haplotypes. These findings indicated no statistically significant associations of MTHFR haplotypes with BC risk.

We also considered the association between haplotypes and major clinicopathological features of BC patients. As presented in Table 5, our data suggested the haplotype composed of the 2 alternative alleles of MTHFR SNPs (rs1801133/T - rs1801131/C) was strongly associated with early age at diagnosis of BC (15.4% of early onset BC vs. 5.1% diagnosed after 40 years; OR: 3.52, 95% CI: 1.23–10.06, p = 0.02). It was also found to be significantly more frequent in familial forms of BC than in sporadic (OR: 4.96, 95% CI: 1.5–16.42, p = 0.0097). Contrariwise, the haplotype with one minor allele rs1801131/C (C-C) was more frequent in sporadic BC than familial (24.5% vs. 7.4%; OR: 0.36, 95% CI: 0.14–0.95, p = 0.042). However, none of these associations reached the Bonferroni-adjusted significance threshold (p > 0,005).

Table 5 Association of ERCC2 and MTHFR haplotypes with clinical characteristics of breast cancer patients

With respect to ERCC2 gene, the haplotype composed of the 2 minor alleles (rs1799793/A - rs13181/C) was significantly associated with positive ER and PR expressions in BC tumors (ER: OR: 0.28, 95% CI: 0.09–0.88, p = 0.03; PR: OR: 0.33, 95% CI: 0.12–0.92, p = 0.017). However, patients who carried G-C haplotype tended to be more likely to develop BC with negative-PR tumors (OR: 2.34, 95% CI: 1.1–4.96, p = 0.02) and with large tumor size, T4 (OR: 23.3, 95% CI: 2.03–267.69, p = 0.017) than did positive PR patients and T1 tumor size ones, respectively. No statistical associations were detected after Bonferroni correction.

Finally, we did not discover any association with MTHFR, ERCC2 haplotypes and other clinicopathological parameters of BC.


The etiology of breast cancer is complex and multifactorial, as sustained by contribution of various environmental and genetic factors. Beside mutations in predisposition genes, the identification of genetic polymorphisms including SNPs in the genes conferring relatively small increment in BC risk could be beneficial for the understanding of the disease mechanisms. Such information could also be of great interest in identifying high risk individuals and in improving cancer prevention strategies.

Accordingly, in our present case-control study, we investigated whether 4 SNPs of ERCC2 and MTHFR genes affect the pathogenesis of BC in Moroccan population.

The results from this study revealed, for the first time, an association of the four polymorphisms with increased risk of breast cancer and/or with disease sub-phenotypes including age at diagnosis, family history, hormone receptor statuses and tumor size in Moroccan BC patients.

The first polymorphism, i.e., ERCC2-rs1799793 was identified as potential risk factor for BC in this work. Indeed, homozygote carriers of the minor allele (A/A genotype, Asn312Asn) were over-represented in the BC cases compared to controls, which would make them at high risk of developing the disease among Moroccan cases.

The variant ERCC2-rs1799793 is a G > A coding polymorphism causing a codon 312 Asp to Asn amino acid exchange in ERCC2 gene. ERCC2 is an essential gene involved in DNA damage repair pathway and whose product, a DNA helicase, is important in the transcription-coupled nucleotide excision repair process that contribute to preserving integrity and stability of the genome. It is well known that the DNA repair ability is an important determinant of the predisposition toward various malignancies [9, 24, 25]. There is increasing data supporting the hypothesis that genetic polymorphisms (SNPs) in DNA repair genes could lead to disorder in DNA repair machinery resulting in accumulation of mutations, and in turn could contribute to increased susceptibility to various types of cancers including BC [17, 25,26,27]. In particular, the functional SNPs ERCC2-rs1799793 and ERCC2-rs13181 enrolled in our study have been previously associated with specific DNA defects, namely defective repair capacity of ultraviolet light-induced DNA damage [28]. Interestingly, data from a study conducted by wolf et al. [29] indicated that both polymorphisms significantly decreased constitutive ERCC2 mRNA levels in lymphocytes of healthy subjects which consequently reduce ERCC2 protein amounts.

Our findings revealed an association of ERCC2-rs1799793/AA genotype with increased risk of BC. These results are corroborated by previous studies in various populations, such as Russians, Mexicans, Chinese, Egyptian, and Taiwanese [17, 30,31,32,33]. However, the results were not unanimous as other studies of this SNP failed to find positive correlation with BC, especially among Caucasians as well as North American, European sub-populations, Chinese, Portuguese, Poland and Australian [9, 34,35,36,37,38,39]. At the opposite, the recessive genotype was reported to be protective in Asian and Chinese populations [9, 33, 36, 40,41,42]. All of these studies agree to take into account ethnic origin, sample characteristics and environmental factors that interact with that variant in the reading of these results.

The frequency of ERCC2-rs1799793 minor allele (A) reported in our Moroccan control group was 27%. According to 1000 Genomes Project Phase 3 data [43], this value was higher than in East-Asian and African American (5 and 10%, respectively) and lower than that of South Asian and European (34 and 36%, respectively).

The analysis of association between ERCC2-rs1799793 and clinicopathological features showed that women cases with homozygote AA genotype were more likely to have ER-positive and PR-positive breast cancer compared with women carrying the GG genotype. It is believed that an over-expression of ER in BC could be involved in the tumorogenesis by stimulating mammary cells proliferation which leads to uncontrolled cell division and accumulation of DNA mutations. Therefore, one of the therapeutic means relies on the use of ER modulators [44]. Our results suggested that ERCC2-rs1799793 could be a potential risk marker for hormone receptor-positive BC in Moroccan population. Otherwise, inconsistent findings were reported in a prior study showing that Chinese women with heterozygous genotype were more prone to develop PR-negative BC compared to wild type genotype [9].

In the current study, we included the coding SNP, ERCC2-rs13181 due to its functional relevance [28]. This polymorphism changes the charge of the amino acid (nucleotide A to C substitution causing a 751 Lys to Gln amino acid change) and is located in a crucial domain of interaction between ERCC2 protein and p44, its helicase activator within the transcription factor TFIIH complex [45]. Despite the overrepresentation of the alternative homozygote genotype CC (Gln751Gln) in subgroup of Moroccan BC patients (OR: 1.84, 95% CI: 0.86–3.91) and which did not reach significant levels, the results may indicate a possible association with BC risk in view of HWE results. Indeed, there was a clear deviation from the HWE in BC subjects for both SNPs (p = 0.007 for ERCC2-rs13181, and p = 0.002 for ERCC2-rs1799793), while there was an accordance with HWE in controls (p = 0.16 and 0.52, respectively). This indicated that no evolutionary change has occurred affecting the distribution of the normal and alternative alleles in general population. Tupikowski et al. [46] reported similar situation on renal cell carcinoma. Some authors claimed that screening with HWE datasets of affected individuals is relatively efficient to detect genes associated to a disease [46, 47].

In our study, it is possible that the influence of ERCC2-rs13181 CC genotype might be masked by the small size of tested groups. Thus, larger studies are warranted to reveal associations that are not immediately apparent.

As reported by other study populations, there is mixed evidence regarding the contribution of ERCC2-rs13181 polymorphism to the risk of BC. Significant associations with increased risk of BC were found in some populations such as Caucasians, African Americans and Indians [35, 36, 41, 48, 49]. At the opposite, other reports stated that there was no evidence of association for populations of China, North America and Europe [17, 30, 33, 34, 36, 38, 40]. These findings suggested, again, a possible role of the environment, ethnic differences and variable genetics backgrounds in cancer development.

The ERCC2-rs13181 minor allele frequency reported in our study in Moroccan controls (25%) was slightly higher than American (21.5%) and lower than Europeean and Asian frequencies (36.4 and 34.7%, respectively) [43].

In regard to clinicopathological variables, our study showed that ERCC2-rs13181 heterozygote genotype was more prevalent in the BC patients with higher tumor size T3-T4 which is a poor prognostic indicator. These results suggested that ERCC2-rs13181 is more associated with the severity of the disease than its risk and may serve as a biomarker for BC progression in Moroccan population.

In the analysis of association between ERCC2-(rs1799793-rs13181) haplotypes and the risk of BC, we inferred that the haplotype defined by the minor alleles A-C may play a substantial role in increasing the risk of BC. Interestingly, the level of significance was higher (p = 0.0002) than it was when the two SNPs were taken individually (p = 0.0069 for ERCC2-rs1799793 and p = 0.11 for ERCC2-rs13181). This result supports a potential correlation of ERCC2-rs13181 with increased BC risk in Moroccan patients in addition to that of ERCC2-rs1799793. This haplotype may be regarded as susceptibility marker to BC in Moroccan patients. However, this result differed from previous studies. Indeed, the same haplotype was associated with marginal risk of BC in North-Eastern Poland population [38], but failed to exert any effect in African Americans [50]. In the latter population, the haplotype defined by major alleles (G-A) was found more frequently among controls than cases [50], while the G-C combination was considered as the most potent risk-conferring haplotype in German population [51].

Otherwise, it appears that these two polymorphisms have low linkage disequilibrium in both Moroccan cases and controls suggesting that they are located in a haplotype block with high rate of recombination between the two loci. These two SNPs could therefore be regarded as two distinct hereditary units. Previous studies have reported similar results in populations of European and African ancestry based on the HapMap data [50], in contrast to the US and Poland populations where they are in linkage disequilibrium [38, 52]. In our study, the A-C haplotype was found to correlate with ER+ and PR+ expression, whereas the G-C haplotype was connected with higher risk of developing a PR negative and high tumor size BC in Moroccan cases. Accordingly, these two haplotypes could be considered as markers of breast cancer prognosis.

The other most relevant result of the current work was the association of MTHFR- rs1801133 (C677T) polymorphism with increased susceptibility for BC in Moroccan patients. We have detected more homozygote carriers (TT) among patients than controls. The genotype distributions of MTHFR-rs1801133 had a somewhat deviation from the HWE in BC cases, but not in healthy controls, giving further evidence of its role in increased BC risk.

MTHFR is the gene encoding methylenetetrahydrofolate reductase enzyme which is involved in folate metabolism. The enzyme assists the irreversible conversion of 5,10-methylenetetrahydrofolate (5,10-MTHF) to 5-methyltetrahydrofolate (5-MTHF). The 5-MTHF, the predominant circulatory form of folate, plays an integral role in DNA synthesis, DNA methylation and DNA repair and maintenance. Deficiency of folate has been shown to result in DNA damage, DNA hypomethylation and reduced DNA repair leading to an increased risk of chromosomal breaks. It has been appreciated that depletion of folate might be linked to the carcinogenesis in multiple cancer conditions through the process cited above. Consequently, the potential influence of MTHFR activity on folate availability makes the MTHFR gene an attractive candidate for cancer predisposition [53].

Accordingly, the most common studied functional MTHFR variants, rs1801133/C677T (Ala222Val) and rs1801131/A1298C (Glu429Ala) have been shown to generate therm-olabile and less active enzyme, with C677T having a higher effect than A1298C variant [54,55,56].

In the case of breast cancer, rs1801133 and rs1801131 have been widely assessed for their implication in increasing breast cancer risk in numerous epidemiological studies. However, conflicting results have been reported depending on the study group. For several studies, there was a significant association between rs1801133 SNP and high risk of BC [53, 57,58,59,60,61], whereas a number of others failed to detect any association [62,63,64,65,66,67,68].

When taking the age at diagnosis and family history into account, we found a strong association of MTHFR-rs1801133 with young age (< 40 years) and with familial form of BC in Moroccan patients. Campbell et al. [69] and Semenza et al. [70] reported similar results of significant association between MTHFR-rs1801133 variant and early onset of breast cancer (before 40 years of age) in English population. Another study showed that the risk estimates were maintained in group of women diagnosed at or before 50 years [71]. Early age of diagnosis is well accepted as a prognostic indicator associated to more aggressive form of BC [72]. Moreover, family history is a particularly major factor associated with increased risk of BC. Our finding are in accordance to an earlier study conducted in Jewish population showing high frequency of hereditary BC in individuals with TT genotype [73] and in Italian BRCA1 mutations carriers [74]. Therefore, our results suggested MTHFR-rs1801133 polymorphism as a real risk modifier in overall Moroccan BC cases, especially in young and familial subgroups.

Otherwise, the MTHFR-rs1801133 had no statistically significant association with clinicopathologic features. This result was supported by a recent report in Indian subjects [59] and by previous findings in Brazilian and Austrian cases. Meanwhile Huang et al. [75] reported only a weak correlation of MTHFR C/T or TT genotype distribution with RE positive status in Taiwanese population.

The frequency of the MTHFR-rs1801133 variant allele amongst healthy Moroccan population (28%) was similar to East-Asian one (29.6%) but lower than those reported in American (47.4%) and European (36.5%) [43].

Regarding the second SNP, MTHFR-rs1801131, while there was no significant association with overall BC risk in Moroccan patients, this polymorphism could be a potential marker for sporadic BC subphenotype and marginally for patients aged over 40 years. Similarly, high risk of sporadic breast cancer associated with MTHFR-rs1801131 was previously reported in Turkish women carrying homozygote variant genotype [76]. Likewise, the presence of the alternative C allele confered an increased risk of breast cancer in sporadic cases of Italian population [74].

In the light of these findings, and although both MTHFR polymorphisms affect the total enzymatic activity, it is not excluded that their dissimilar contributions in modulating breast cancer risk direction and extend may depend on interactions with other still unknown endogenous or exogenous factors.

We further evaluated the contribution of MTHFR haplotypes generated by MTHFR-rs1801133 and MTHFR-rs1801131 SNPs to BC risk. The current findings displayed a frequency of 6 and 7.5% for the haplotype T-C in Moroccan healthy population and in BC cases. Different data were reported in previous studies depending on ethnic origin. The estimated frequency was reported to be zero in German, Spanish and Japanese populations [67, 77, 78]. Nevertheless, a recent Arabic study showed slightly higher frequency in Jordanian population (8.3%) and lower frequency (3.6%) in matched BC cases [79].

Haplotype analysis inferred that none of the MTHFR haplotypes was significantly associated with overall BC risk in our population, although we have noticed that carriers of the haplotype defined by the minor alleles (T-C) were 1.36 times more likely to have BC. However, this haplotype exhibited a positive correlation with familial form and with early onset of the disease. At the opposite, the C-C haplotype showed higher representation in sporadic BC subgroup. Other clinical conditions of BC were independent of MTHFR haplotypes in Moroccan patients.

Previous studies that investigated the contributory role of MTHFR haplotypes in BC development have produced inconclusive results. A borderline line significant protection was observed for the C-C haplotype in German and East asian populations [67, 80], while the C-A haplotype was protective in South-Eastern European population [54]. Carriers of the T-C haplotype in Caucasians [80] and of the T-A haplotype in Jordanian [79] were more prone to develop BC.

Interestingly, the Lewontin’s estimate was consistent with no linkage disequilibrium in both cases and controls, suggesting that both SNP were independent of each other. It seems likely that there is a specific linkage disequilibrium pattern in Moroccan population for these MTHFR SNPs, suggesting that they act independently to affect BC susceptibility. LD patterns at these loci appear to be population dependant. LD was strong in populations of Europe, Brazil, Pakistan and china [78, 81, 82] and much smaller in Mexican and African populations [78]. In contrast, these two variants are genetically independent in Russian and Puerto Rican populations [82, 83].


To the best of our knowledge, this is the first study assigning increased BC risk to ERCC2-rs1799793 (Asn312Asn) polymorphism and the corresponding haplotype determined by Asn312-Gln751 codons in Moroccan population.

The other finding of special interest is the association between MTHFR-rs1801133 (Val222Val) with increased risk of BC. Our results suggested that ERCC2-rs1799793 and MTHFR-rs1801133 represent suitable tool for assessing susceptibility to breast cancer in Moroccan population and prognosis.

For the two other SNPs investigated in this study, it is likely that either they do not contribute to BC risk or, more likely, their influence is small and can be detected only in larger samples. Thus, it is strongly recommended to reproduce these results on a larger number of participants.



Breast cancer


Confidence interval


Estrogen receptor


Family history


Human epidermal growth factor receptor 2


Hardy-Weinberg equilibrium


Linkage disequilibrium




Progesterone receptor




Single nucleotide polymorphism


  1. 1.

    Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in (GLOBOCAN) 2012. Int J Cancer. 2015;136:E359–86.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Global Burden of Disease Cancer Collaboration, Fitzmaurice C, Allen C, Barber RM, Barregard L, Bhutta ZA, Brenner H, Dicker DJ, Chimed-Orchir O, Dandona R, et al. Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 32 cancer groups, 1990 to 2015: a systematic analysis for the global burden of disease study. JAMA Oncol. 2017;3:524–48.

    Article  Google Scholar 

  3. 3.

    Chouchane L, Boussen H, Sastry KSR. Breast cancer in Arab populations: molecular characteristics and disease management implications. Lancet Oncol. 2013;14:417–24.

    Article  Google Scholar 

  4. 4.

    Fondation Lalla Salma Prévention et traitement des Cancers. Registre des cancers de la région du Grand Casablanca pour la période 2008–2012. 2016.

    Google Scholar 

  5. 5.

    Khalis M, Rhazi K, Charaka H, Chajès V, Rinaldi S, Nejjari C, Romieu I, Charbotel B. Female breast cancer incidence and mortality in Morocco: comparison with other countries. Asian Pac J Cancer Prev. 2016;17:5211–6.

    PubMed  PubMed Central  Google Scholar 

  6. 6.

    Hamdi Y, Soucy P, Adoue V, Michailidou K, Canisius S, Lemaçon A, Droit A, Andrulis IL, Anton-Culver H, Arndt V, et al. Association of breast cancer risk with genetic variants showing differential allelic expression: identification of a novel breast cancer susceptibility locus at 4q21. Oncotarget. 2016;7:80140–63.

    Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Shi M, O'Brien KM, Sandler DP, Taylor JA, Zaykin DV, Weinberg CR. Previous GWAS hits in relation to young-onset breast cancer. Breast Cancer Res Treat. 2017;161:333–44.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Véron A, Blein S, Cox DG. Genome-wide association studies and the clinic: a focus on breast cancer. Biomark Med. 2014;8:287–96.

    Article  PubMed  Google Scholar 

  9. 9.

    Wang T, Wang H, Guo H, Yang S, Zhu G, Guo H, et al. Polymorphisms in the DNA repair gene ERCC2/XPD and breast cancer risk: a HapMap-based case–control study among Han women in a chinese less-developed area. Genet Test Mol Biomarkers. 2014;18:703–10.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Lanara Z, Giannopoulou E, Fullen M, Kostantinopoulos E, Nebel J-C, Kalofonos HP, et al. Comparative study and meta-analysis of meta-analysis studies for the correlation of genomic markers with early cancer detection. Hum Genomics. 2013;7:14.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Michailidou K, Hall P, Anna GN, Ghoussaini M, Dennis J, Milne RL, et al. Large-scale genotyping identifies 41 new loci associated with breast cancer risk. Nat Genet. 2013;45:353–61.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Michailidou K, Beesley J, Lindstrom S, Canisius S, Dennis J, Lush MJ, et al. Genome-wide association analysis of more than 120,000 individuals identifies 15 new susceptibility loci for breast cancer. Nat Genet. 2015;47:373–80.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Jouhadi H, Tazzite A, Azeddoug H, Naim A, Nadifi S, Benider A. Clinical and pathological features of BRCA1/2 tumors in a sample of high-risk Moroccan breast cancer patients. BMC Res Notes. 2016;9:248.

    Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Laraqui A, Uhrhammer N, Idriss L-A, Rhaffouli HEL, Baghdadi J, Dehayni M, et al. Mutation screening of the BRCA1 gene in early onset and familial breast/ovarian cancer in Moroccan population. Int J Med Sci. 2013;10:60–7.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Rahoui J, Sbitti Y, Touil N, Laraqui A, Ibrahimi A, Rhrab B, et al. The single nucleotide polymorphism +936 C/T VEGF is associated with human epidermal growth factor receptor 2 expression in Moroccan breast cancer women. Med Oncol. 2014;31:336.

    Article  PubMed  Google Scholar 

  16. 16.

    Marouf C, Göhler S, Filho M, Hajji O, Hemminki K, Nadifi S, et al. Analysis of functional germline variants in APOBEC3 and driver genes on breast cancer risk in Moroccan study population. BMC Cancer. 2016;16:165.

    Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Gómez-Díaz B, DE LA Luz Ayala-Madrigal M, Gutiérrez-Angulo M, Valle-Solis AE, Linares-González LM, González-Guzmán R, Cruz-Guillén D, Cedeño-Garcidueñas AL, Canto P, López-Hernández LB. Analysis of ERCC1 and ERCC2 gene variants in osteosarcoma, colorectal and breast cancer. Oncol Lett. 2015;9:1657–61.

    Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Xie SZ, Liu ZZ, Yu J, Liu L, Wang W, Xie DL, et al. Association between the MTHFR C677T polymorphism and risk of cancer: evidence from 446 case–control studies. Tumor Biol. 2015;36:8953–72.

    CAS  Article  Google Scholar 

  19. 19.

    Diakite B, Tazzite A, Hamzi K, Jouhadi H, Nadifi S. Methylenetetrahydrofolate reductase C677T polymorphism and breast cancer risk in Moroccan women. Afr Health Sci. 2012;12:204–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988;16:1215.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Dudbridge F. Likelihood-based association analysis for nuclear families and unrelated subjects with missing genotype data. Hum Hered. 2008;66:87–98.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Solé X, Guinó E, Valls J, Iniesta R, Moreno V. SNPStats: a web tool for the analysis of association studies. Bioinform Oxf Engl. 2006;22:1928–9.

    Article  Google Scholar 

  23. 23.

    Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–5.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Katarzyna P-S, Scott RJ, Górski B, Cybulski C, Kurzawski G, Dymerska D, et al. Polymorphisms in nucleotide excision repair genes and susceptibility to colorectal cancer in the polish population. Mol Biol Rep. 2015;42:755–64.

    Article  Google Scholar 

  25. 25.

    Srivastava K, Srivastava A, Mittal B. Polymorphisms in ERCC2, MSH2, and OGG1 DNA repair genes and gallbladder cancer risk in a population of Northern India. Cancer. 2010;116:3160–9.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Chu H, Gu D, Xu M, Xu Z, Gong Y, Gong W, et al. A genetic variant in ERCC2 is associated with gastric cancer prognosis in a Chinese population. Mutagenesis. 2013;28:441–6.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Goode EL, Ulrich CM, Potter JD. Polymorphisms in DNA repair genes and associations with cancer risk. Cancer Epidemiol Biomark Prev. 2002;11:1513–30.

    CAS  Google Scholar 

  28. 28.

    Au WW, Salama SA, Sierra-Torres CH. Functional characterization of polymorphisms in DNA repair genes using cytogenetic challenge assays. Environ Health Perspect. 2003;111:1843–50.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Wolfe KJ, Wickliffe JK, Hill CE, Paolini M, Ammenheuser MM, Abdel-Rahman SZ. Single nucleotide polymorphisms of the DNA repair gene XPD/ERCC2 alter mRNA expression. Pharmacogenet Genomics. 2007;17:897–905.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Shadrina AS, Ermolenko NA, Boyarskikh UA, Sinkina TV, Lazarev AF, Petrova VD, et al. Polymorphisms in DNA repair genes and breast cancer risk in Russian population: a case-control study. Clin Exp Med. 2016;16:21–8.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Hussien YM, Gharib AF, Awad HA, Karam RA, Elsawy WH. Impact of DNA repair genes polymorphism (XPD and XRCC1) on the risk of breast cancer in Egyptian female patients. Mol Biol Rep. 2012;39:1895–901.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Wang HC, Liu CS, Wang CH, Tsai RY, Tsai CW, Wang RF, et al. Significant association of XPD Asp312Asn polymorphism with breast cancer in Taiwanese patients. Chin J Physiol. 2010;53:130–5.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Zhang L, Zhang Z, Yan W. Single nucleotide polymorphisms for DNA repair genes in breast cancer patients. Clin Chim Acta. 2005;359:150–5.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    He BS, Xu T, Pan YQ, Wang HJ, Cho WC, Lin K, Sun HL, Gao TY, Wang SK. Nucleotide excision repair pathway gene polymorphisms are linked to breast cancer risk in a Chinese population. Oncotarget. 2016;7:84872–82.

    PubMed  PubMed Central  Google Scholar 

  35. 35.

    Yan Y, Liang H, Light M, Li T, Deng Y, Li M, et al. XPD Asp312Asn and Lys751Gln polymorphisms and breast cancer susceptibility: a meta-analysis. Tumor Biol. 2014;35:1907–15.

    CAS  Article  Google Scholar 

  36. 36.

    Pabalan N, Ofelia FP, Sung L, Jarjanazi H, Ozcelik H. Meta-analysis of two ERCC2 (XPD) polymorphisms, Asp312Asn and Lys751Gln, in breast cancer. Breast Cancer Res Treat. 2010;124:531–41.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Costa S, Pinto D, Pereira D, Rodrigues H, Jorge C-T, Medeiros R, et al. DNA repair polymorphisms might contribute differentially on familial and sporadic breast cancer susceptibility: a study on a Portuguese population. Breast Cancer Res Treat. 2007;103:209–17.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Dębniak T, Scott RJ, Huzarski T, Byrski T, Masojć B, van de Wetering T, et al. XPD common variants and their association with melanoma and breast cancer risk. Breast Cancer Res Treat. 2006;98:209–15.

    Article  PubMed  Google Scholar 

  39. 39.

    Kuschel B, Georgia CT, Spurdle AB, Chen X, Hopper JL, Giles GG, et al. Common polymorphisms in ERCC2 (Xeroderma pigmentosum D) are not associated with breast cancer risk. Cancer Epidemiol Biomark Prev. 2005;14:1828–31.

    CAS  Article  Google Scholar 

  40. 40.

    Zhao R, Ying MF. Association between ERCC1 and ERCC2 polymorphisms and breast cancer risk in a Chinese population. Genet Mol Res. 2016;

  41. 41.

    Jiang Z, Li C, Xu Y, Cai S, Wang X. Associations between XPD polymorphisms and risk of breast cancer: a meta-analysis. Breast Cancer Res Treat. 2010;123:203–12.

    Article  PubMed  Google Scholar 

  42. 42.

    Yao L, Qiu LX, Yu L, Yang Z, Yu X-J, Zhong Y, et al. The association between ERCC2 Asp312Asn polymorphism and breast cancer risk: a meta-analysis involving 22,766 subjects. Breast Cancer Res Treat. 2010;123:227–31.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    IGSR: The international genome sample resource. Providing ongoing support for the 1000 genome project data.

  44. 44.

    McCafferty MP, McNeil RE, Miller N, Kerin MJ. Interactions between the estrogen receptor, its cofactors and microRNAs in breast cancer. Breast Cancer Res Treat. 2009;116:425–32.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Coin F, Marinoni J-C, Rodolfo C, Fribourg S, Pedrini A, Egly J-M. Mutations in the XPD helicase gene result in XP and TTD phenotypes, preventing interaction between XPD and the p44 subunit of TFIIH. Nat Genet. 1998;20:184–8.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Tupikowski K, Partyka A, Kolodziej A, Dembowski J, Debinski P, Halon A, et al. CTLA-4 and CD28 genes’ polymorphisms and renal cell carcinoma susceptibility in the polish population – a prospective study. Tissue Antigens. 2015;86:353–61.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Lee WC. Searching for disease-susceptibility loci by testing for Hardy-Weinberg Disequilibrium in a Gene Bank of affected individuals. Am J Epidemiol. 2003;158:397–400.

    Article  PubMed  Google Scholar 

  48. 48.

    Qiu LX, Yao L, Zhang J, Zhu XD, Zhao XM, Xue K, et al. XPD Lys751Gln polymorphism and breast cancer susceptibility: a meta-analysis involving 28,709 subjects. Breast Cancer Res Treat. 2010;124:229–35.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Mitra AK, Singh N, Garg VK, Chaturvedi R, Sharma M, Rath SR. Statistically significant association of the single nucleotide polymorphism (SNP) rs13181 (ERCC2) with predisposition to squamous cell carcinomas of the head and neck (SCCHN) and breast cancer in the north Indian population. J Exp Clin Cancer Res. 2009;28:104.

    Article  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Smith TR, Levine EA, Freimanis RI, Akman SA, Allen GO, Hoang KN, et al. Polygenic model of DNA repair genetic polymorphisms in human breast cancer risk. Carcinogenesis. 2008;29:2132–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Justenhoven C, Hamann U, Pesch B, Harth V, Rabstein S, Baisch C, Vollmert C, Illig T, Ko YD, Brüning T, et al. ERCC2 genotypes and a corresponding haplotype are linked with breast cancer risk in a German population. Cancer Epidemiol Biomark Prev. 2004;13:2059–64.

    CAS  Google Scholar 

  52. 52.

    Jiao L, Hassan MM, Bondy ML, Abbruzzese JL, Evans DB, Li D. The XPD Asp312Asn and Lys751Gln polymorphisms, corresponding haplotype, and pancreatic cancer risk. Cancer Lett. 2007;245:61–8.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Gonzales MC, Yu P, Shiao PSK. MTHFR gene polymorphism mutations and air pollution as risk factors for breast cancer: a metaprediction study. Nurs Res. 2017;66:152–63.

    Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Papandreou CN, Doxani C, Zdoukopoulos N, Vlachostergios PJ, Hatzidaki E, Bakalos G, et al. Evidence of association between methylenetetrahydrofolate reductase gene and susceptibility to breast cancer: a candidate gene association study in a South-Eastern European population. DNA Cell Biol. 2012;31:193–8.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Weisberg IS, Jacques PF, Selhub J, Bostom AG, Chen Z, Ellison CR, et al. The (1298A→C) polymorphism in methylenetetrahydrofolate reductase (MTHFR): in vitro expression and association with homocysteine. Atherosclerosis. 2001;156:409–15.

    CAS  Article  PubMed  Google Scholar 

  56. 56.

    Yamada K, Chen Z, Rozen R, Matthews RG. Effects of common polymorphisms on the properties of recombinant human methylenetetrahydrofolate reductase. Proc Natl Acad Sci. 2001;98:14853–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    He L, Shen Y. MTHFR C677T polymorphism and breast, ovarian cancer risk: a meta-analysis of 19,260 patients and 26,364 controls. Onco Targets Ther. 2017;10:227–38.

    Article  PubMed  PubMed Central  Google Scholar 

  58. 58.

    Naushad S, Divya C, Ramaiah JM, Hussain T, Alrokayan SA, Kutala V. Population-level diversity in the association of genetic polymorphisms of one-carbon metabolism with breast cancer risk. J Community Genet. 2016;7:279–90.

    Article  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Waseem M, Hussain SR, Kumar S, Serajuddin M, Mahdi F, Sonkar SK, Bansal C, Ahmad MK. Association of MTHFR (C677T) gene polymorphism with breast cancer in North India. Biomark Cancer. 2016;8:111–7.

    Article  PubMed  PubMed Central  Google Scholar 

  60. 60.

    Lu Q, Jiang K, Li Q, Ji YJ, Chen WL, Xue XH. Polymorphisms in the MTHFR gene are associated with breast cancer risk and prognosis in a Chinese population. Tumor Biol. 2015;36:3757–62.

    CAS  Article  Google Scholar 

  61. 61.

    Qi X, Ma X, Yang X, Fan L, Zhang Y, Zhang F, et al. Methylenetetrahydrofolate reductase polymorphisms and breast cancer risk: a meta-analysis from 41 studies with 16,480 cases and 22,388 controls. Breast Cancer Res Treat. 2010;123:499–6.

    CAS  Article  PubMed  Google Scholar 

  62. 62.

    Kaya EF, Karakus N, Ulusoy AN, Özaslan C, Kara N. Association of the MTHFR gene (C677T) polymorphism with breast cancer in a Turkish population. Oncol Res Treat. 2016;39:534–8.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Song A, Zhao L, Li Y, Wu L, Li Y, Liu X, et al. Haplotypes of the MTHFR gene are associated with an increased risk of breast cancer in a Han Chinese population in Gansu province. IUBMB Life. 2016;68:526–34.

    CAS  Article  PubMed  Google Scholar 

  64. 64.

    Kakkoura MG, Demetriou CA, Loizidou MA, Loucaides G, Neophytou I, Marcou Y, et al. Single-nucleotide polymorphisms in one-carbon metabolism genes, Mediterranean diet and breast cancer risk: a case–control study in the Greek-Cypriot female population. Genes Nutr. 2015;10:5.

    Article  PubMed Central  Google Scholar 

  65. 65.

    Qiu LX, Zhang J, Li WH, Zhang QL, Yu H, Wang BY, et al. Lack of association between methylenetetrahydrofolate reductase gene A1298C polymorphism and breast cancer susceptibility. Mol Biol Rep. 2011;38:2295–9.

    CAS  Article  PubMed  Google Scholar 

  66. 66.

    Ziva Cerne J, Stegel V, Gersak K, Novakovic S. Lack of association between methylenetetrahydrofolate reductase genetic polymorphisms and postmenopausal breast cancer risk. Mol Med Rep. 2011;4:175–9.

    PubMed  Google Scholar 

  67. 67.

    Justenhoven C, Hamann U, Pierl CB, Rabstein S, Pesch B, Harth V, et al. One-carbon metabolism and breast cancer risk: no association of MTHFR, MTR, and TYMS polymorphisms in the GENICA study from Germany. Cancer Epidemiol Prev Biomarkers. 2005;14:3015–8.

    CAS  Article  Google Scholar 

  68. 68.

    Langsenlehner U, Krippl P, Renner W, Babak Y-B, Wolf G, Wascher TC, et al. The common 677C>T gene polymorphism of methylenetetrahydrofolate reductase gene is not associated with breast cancer risk. Breast Cancer Res Treat. 2003;81:169–72.

    CAS  Article  PubMed  Google Scholar 

  69. 69.

    Campbell IG, Baxter SW, Eccles DM, Choong DYH. Methylenetetrahydrofolate reductase polymorphism and susceptibility to breast cancer. Breast Cancer Res. 2002;4:R14.

    Article  PubMed  PubMed Central  Google Scholar 

  70. 70.

    Semenza JC, Delfino RJ, Ziogas A, Hoda AC. Breast cancer risk and methylenetetrahydrofolate reductase polymorphism. Breast Cancer Res Treat. 2003;77:217–23.

    CAS  Article  PubMed  Google Scholar 

  71. 71.

    Kotsopoulos J, Zhang WW, Zhang S, David M, Trudeau M, Zhang P, et al. Polymorphisms in folate metabolizing enzymes and transport proteins and the risk of breast cancer. Breast Cancer Res Treat. 2008;112:585–93.

    CAS  Article  PubMed  Google Scholar 

  72. 72.

    Kaabi B, Belaaloui G, Benbrahim W, Hamizi K, Sadelaoud M, Toumi W, et al. ADRA2A germline gene polymorphism is associated to the severity, but not to the risk of breast cancer. Pathol Oncol Res. 2016;22:357–65.

    CAS  Article  PubMed  Google Scholar 

  73. 73.

    Gershoni-Baruch R, Dagan E, Israeli D, Kasinetz L, Kadouri E, Friedman E. Association of the C677T polymorphism in the MTHFR gene with breast and/or ovarian cancer risk in Jewish women. Eur J Cancer. 2000;36:2313–6.

    CAS  Article  PubMed  Google Scholar 

  74. 74.

    Pepe C, Guidugli L, Sensi E, Aretini P, D'Andrea E, Montagna M, et al. Methyl group metabolism gene polymorphisms as modifier of breast cancer risk in Italian BRCA1/2 carriers. Breast Cancer Res Treat. 2007;103:29–36.

    CAS  Article  PubMed  Google Scholar 

  75. 75.

    Huang MY, Wang YH, Chen FM, Lee SC, Fang WY, Cheng TL, et al. Multiple genetic polymorphisms of GSTP1313AG, MDR1 3435CC, and MTHFR 677CC highly correlated with early relapse of breast cancer patients in Taiwan. Ann Surg Oncol. 2007;15:872–80.

    Article  PubMed  Google Scholar 

  76. 76.

    Ergul E, Sazci A, Utkan Z, Canturk NZ. Polymorphisms in the MTHFR gene are associated with breast cancer. Tumor Biol. 2004;24:286–90.

    Article  Google Scholar 

  77. 77.

    Guillem VM, Collado M, Terol MJ, Calasanz MJ, Esteve J, Gonzalez M, et al. Role of MTHFR (677, 1298) haplotype in the risk of developing secondary leukemia after treatment of breast cancer and hematological malignancies. Leukemia. 2007;21:1413–22.

    CAS  Article  PubMed  Google Scholar 

  78. 78.

    Shi M, Caprau D, Romitti P, Christensen K, Murray JC. Genotype frequencies and linkage disequilibrium in the CEPH human diversity panel for variants in folate pathway genes MTHFR, MTHFD, MTRR, RFC1, and GCP2. Birth Defects Res Part Clin Mol Teratol. 2003;67:545–9.

    CAS  Article  Google Scholar 

  79. 79.

    Awwad N, Yousef AM, Abuhaliema A, Abdalla I, Yousef M. Relationship between genetic polymorphisms in MTHFR (C677T, A1298C and their haplotypes) and the incidence of breast cancer among Jordanian females - case-control study. Asian Pac J Cancer P. 2015;16:5007–11.

    Article  Google Scholar 

  80. 80.

    Zhong S, Chen Z, Yu X, Li W, Tang J, Zhao J. A meta-analysis of genotypes and haplotypes of methylenetetrahydrofolate reductase gene polymorphisms in breast cancer. Mol Biol Rep. 2014;41:5775–85.

    CAS  Article  PubMed  Google Scholar 

  81. 81.

    Fan S, Yang B, Zhi X, Wang Y, Zheng Q, Sun G. Combined genotype and haplotype distributions of MTHFR C677T and A1298C polymorphisms: a cross-sectional descriptive study of 13,473 Chinese adult women. Medicine (Baltimore). 2016;95:1–6.

    Article  Google Scholar 

  82. 82.

    Huang T, Tucker KL, Lee YC, Crott JW, Parnell LD, Shen J, et al. Methylenetetrahydrofolate reductase variants associated with hypertension and cardiovascular disease interact with dietary polyunsaturated fatty acids to modulate plasma homocysteine in Puerto Rican adults. J Nutr. 2011;141:654–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  83. 83.

    Trifonova EA, Spiridonova MG, Gabidulina TV, Urnov FD, Puzyrev VP, Stepanov VA. Analysis of the MTHFR gene linkage disequilibrium structure and association of polymorphic gene variants with coronary atherosclerosis. Russ J Genet. 2012;48:1035–47.

    CAS  Article  Google Scholar 

Download references


We would like to thank all the patients and control subjects for their participation in this study. We gratefully acknowledge all the Blood Transfusion Center staff in Oujda city for their valuable help.


Internal support for this study was provided by the Mohammed First University-Oujda, Morocco. The funding body did not play any role in the study design, collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author information




NB conceived and coordinated the study and drafted manuscript. RM performed the statistical analyses, interpreted the results and wrote the manuscript. HH carried out SNP molecular genotyping. ZB and TE were involved in the recruitment of the patients and provided the clinical data. SA contributed to clinical data management. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Rahma Melki.

Ethics declarations

Ethics approval and consent to participate

The study protocol was reviewed and approved by the Ethics Committee of Mohammed VI University Hospital Center of Marrakech in Morocco. Informed and written consent was obtained from all individual participants included in the study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Hardi, H., Melki, R., Boughaleb, Z. et al. Significant association between ERCC2 and MTHR polymorphisms and breast cancer susceptibility in Moroccan population: genotype and haplotype analysis in a case-control study. BMC Cancer 18, 292 (2018).

Download citation


  • Breast cancer
  • ERCC2
  • Genetics
  • Single nucleotide polymorphism
  • Haplotype
  • Moroccan population