Influence of chemotherapeutic drug-related gene polymorphisms on toxicity and survival of early breast cancer patients receiving adjuvant chemotherapy

Background We investigated whether GSTT1 (“null” allele), GSTM1 (“null”allele), GSTP1 (A313G), RFC1 (G80A), MTHFR (C677T), TS (2R/3R) polymorphisms were associated with toxicity and survival in patients with early breast cancer (EBC) treated with adjuvant chemotherapy (CT). Methods This prospective trial included patients with stage I–III BC subjected to CT with CMF or FEC regimens. PCR-RFLP was performed for MTHFR, RFC1 and GSTP1, while PCR for TS, GSTT1 and GSTM1 genes. Results Among the 244 patients consecutively enrolled, 48.7% were treated with FEC and 51.3% with CMF. Patients with TS2R/3R genotype showed less frequently severe neutropenia (G3/G4) than those with TS2R/2R and 3R/3R genotype (p = 0.038). Patients with MTHFRCT genotype had a higher probability of developing severe neutropenia than those with MTHFR CC genotype (p = 0.043). Patients with RFC1GG or GSTT1-null genotype or their combination (GSTT1-null/RFC1GG) were significantly associated with a shorter disease free survival (DFS) (p = 0.009, p = 0.053, p = 0.003, respectively) and overall survival (OS) (p = 0.036, p = 0.015, p = 0.005, respectively). Multivariate analysis confirmed the association of RFC1GG genotype with a shorter DFS (p = 0.018) and of GSTT1-null genotype of a worse OS (p = 0.003), as well as for the combined genotypes GSTT1-null/RFC1GG, (DFS: p = 0.004 and OS: p = 0.003). Conclusions Our data suggest that TS2R/2R and 3R/3R or MTHFR CT genotypes have a potential role in identifying patients with greater risk of toxicity to CMF/FEC and that RFC1 GG and GSTT1-null genotypes alone or in combination could be important markers in predicting clinical outcome in EBC patients. Electronic supplementary material The online version of this article (doi:10.1186/s12885-017-3483-2) contains supplementary material, which is available to authorized users.


Background
Breast cancer (BC) currently accounts for 20% of all female cancers worldwide and is the most frequent malignancy occurring in women [1]. There is convincing evidence that adjuvant systemic chemotherapy (AC) increases survival of patients with BC [2]. AC imparted a statistically significant reduction in the risk of BC relapse and death at 5 years of follow-up (with a hazard reduction of approximately 25%), and combination chemotherapy was found to be significantly more effective than singleagent therapy [3]. Trials included more than 15 years of follow-up and led to the conclusion that AC conferred benefit to both premenopausal and postmenopausal patients and also to node-positive and node-negative patients [4]. In general, approximately one of every four recurrences and one of seven deaths is avoided annually by adjuvant chemotherapy [5].
Among the treatments used in this adjuvant setting, the combination of cyclophosphamide (CP), methotrexate (MTX) and 5-fluorouracil (5-FU) (CMF treatment) or the combination of 5-FU, anthracycline-based chemotherapy (adriamycin or its analogue epirubicin) and CP (FAC/FEC treatment) are the most commonly used. Although the benefit of BC chemotherapy has been demonstrated, these drugs have shown the ability to induce DNA damage in eukaryotic cells [6,7] and, consequently, chemotherapy treatment involves a risk of provoking DNAdamage even in proliferative non-cancer cells [8] therefore leading to a marked toxicity state. Adverse events represent an important physical, psychological and financial burden for the patient and society since up to 15% of the patients receiving FEC will experience at least one serious adverse event [9,10]. Besides toxicity, another major clinical problem encountered during adjuvant CMF or FEC treatments is BC recurrence of therapeutically resistant disease and thus affecting the long-term outcome of the patient. Significant variability in drug response may occur among cancer patients treated with the same medications [11].
Germline genetic variation in drug metabolizing enzymes and transporters is thought to contribute to the observed inter-individual variation in treatment toxicity and/ or efficacy [12]. Recently, pharmacogenomic studies have elucidated the inherited nature of these differences in drug disposition and effects, thereby providing a stronger scientific basis for optimizing drug therapy according to each patient's genetic constitution. Candidate genes are thymidylate synthase (TS), 5, 10-methylenetetrahydrofolate reductase (MTHFR), the reducer folate carrier (RFC1) and glutathione-S-transferases (GSTs), involved in CMF or FEC adjuvant chemotherapies transport and/or metabolism, or being targets of such drugs, as it is shown in Fig. 1. TS is an enzyme implicated in the conversion of deoxyuridine monophosphate (dUMP) into deoxythymidine monophosphate (dTMP), which is essential in DNA synthesis. The human TS gene (hTS) is polymorphic with either double (2R) or triple (3R) tandem repeats of a 28 base-pair sequence downstream of the cap site in the 5′ terminal regulatory region [13]. In vitro studies, the activity of a reporter gene linked to the 5′ terminal fragment of the hTS gene with triple (3R) tandem repeats was 2.6 times higher than that with double (2R) tandem repeats [14]. Thus, this polymorphic region TS 2R/3R appears to be functional and may modulate TS gene expression. MTHFR is an enzyme responsible for the Fig. 1 Metabolism of chemotherapeutic drugs-related gene polymorphisms. In cancer cells 5-FU is converted to 5-fluorodeoxyuridine monophosphate (5-FdUMP). 5-FdUMP inhibits the DNA synthesis by competing with deoxyuridine monophosphate (dUMP) for binding to thymidylate synthase (TS) in a complex that is stabilized by the reduced folate 5,10-methylene tetrahydrofolate. 5-FU can also inhibit RNA synthesis in a pathway that involves its metabolism to 5-fluorouridinemonophosphate (5-FUMP) and subsequent conversion to 5-fluorouridine triphosphate (5-FUTP) via 5-fluorouridine diphosphate (5-FUDP). The main effect of cyclophosphamide is due to its metabolite phosphoramide mustard that forms DNA crosslinks both between and within DNA strands at guanine N-7 positions (known as interstrand and intrastrand crosslinkages, respectively). This is irreversible and leads to cell apoptosis. Anthracyclines inhibit DNA and RNA synthesis by intercalating between base pairs of the DNA/RNA strand, thus preventing the replication of rapidly growing cancer cells. In addition, they can generate reactive oxygen species (ROS) damaging DNA, proteins and cell membranes. Glutathione S-transferases (GSTs) catalyse the detoxification of alkylating agents used in chemotherapy and/or ROS metabolization of vitamin B9 (folate), which is required for DNA synthesis. A known MTHFR gene polymorphism consists of a 677C > T transition, in exon 4, which results in an alanine to valine substitution in the predicted catalytic domain of MTHFR. This substitution renders the enzyme thermolabile, and homozygotes and heterozygotes have about 70 and 35% reduced enzyme activity, respectively [15]. RFC1 is a major MTX transporter whose impaired function has been recognized as a frequent mechanism of antifolate resistence [16]. Different gene alterations affecting RFC1 transport properties were found in cell lines selected for antifolate resistance [17]. A polymorphism G > A at position 80 in exon 2 of RFC1 gene which replaces His by Arg at position 27 of the RFC1 protein was identified. A recent study implied an effect of G > A80 in combination with C > T677 in MTHFR on plasma folate levels and homocysteine pools [18]. It is known that the mechanism of cytotoxicity with chemotherapy is through the generation of reactive oxygen species (ROS) and their by-products. The reactive molecules responsible for cytotoxicity of these therapies are subject to enzymatic removal, and variability of cells in sensitivity to therapy could depend, at least in part, on the availability and activity of specific metabolizing enzymes. GSTs enzymes are an important cellular defence system that protects cells from chemical injury by catalyzing conjugation of reactive electrophilic molecules with glutathione (GSH). GSTs catalyze the detoxification of some alkylating agents used in chemotherapy and detoxification of products of reactive oxidation [19]. GSTs M1 and T1 have been shown to have activity toward lipid hydroperoxides [20], and individuals lacking each of these enzymes (null allele) may have reduced removal of secondary organic oxidation products produced by cancer therapy and thus may have better prognoses. The pi-class human GST (GSTP1) besides playing a role in protection from oxidative damage was shown to catalyze GSH conjugation of reactive cyclophosphamide metabolites in vitro assays [21]. The present study aimed at investigating the association between TS 2R/3R, MTHFR C677T, RFC1 G80A and GSTT1 null, GSTM1 null or GSTP1 A313G polymorphisms with toxicity, disease free survival (DFS) and overall survival (OS) in Caucasian patients with early BC treated with CMF or FEC regimens.

Study population
This prospective study was conducted in patients with a histological diagnosis of stage I-III BC treated with conservative surgery or mastectomy, and subjected to adjuvant chemotherapy with CMF or FEC regimens. Tumor staging followed the TNM-AJCC classification [22] and the pTNM was obtained after classical pathological examination. Patients with metastatic disease and with other previous tumors were excluded from this study. Recorded clinical and pathological features for each patient included: age, menopausal status, histology, grade, stage, estrogen receptors (ER) and progesterone receptor (PgR) status, Ki67, p53, HER2 and medical adjuvant therapy. ER, PgR, Ki67, p53 and HER2 status were assessed at the time of surgery on formalin-fixed paraffin-embedded tissue blocks of the primary tumor in the Pathology Department of the University of Perugia. We used the following cut-off for considering Ki 67 positive >14%, [23] p53 positive ≥ 1%, Her2 positive IHC 3+ or IHC 2+ and FISH amplified. Written informed consent was obtained by all patients and the study was reviewed and approved by the institution's Ethics Committee in accordance with the principles established in the Helsinki declaration.

Chemotherapy regimen
Treatment combined regimen was as follows: CMF (cyclophosphamide 600 mg/m 2 , MTX 40 mg/m 2 and 5fluorouracil 600 mg/m 2 ) administered on day 1 and 8 each 4 weeks, for 6 cycles; FEC (5-fluorouracil 600 mg/ m 2 , 4-epirubicin 90 mg/m 2 and cyclophosphamide 600 mg/m 2 ) administered on day 1, every 21 days, for 6 cycles. Physical examination and a full blood counts were performed after each chemotherapy cycle. Hepatic and renal function tests were assessed at baseline and repeated before each cycle of treatment. All patients who had received at least one course of chemotherapy were evaluated for toxicity. Toxicity was scored every 3 weeks according to the Common Toxicity Criteria of the National Cancer Institute (NCI-CTC, version 2.0) [24].
We defined "severe toxicity" as hematological or gastrointestinal toxicity of grade 3-4.

Genotyping analysis
Genomic DNA was extracted from 200 μL of whole blood using the Qiamp blood kit (Qiagen, Milan, Italy) according to the manufacturer's instructions. Polymorphisms were characterized using the PCR-RFLP for genotyping analyses of MTHFR, RFC1 and GSTP1, while PCR was used for TS polymorphism determination. Multiplex PCR was used to simultaneously amplify GSTT1 and GSTM1, with albumin as a control gene. All primers used in this study were designed by using Primer express 2.0 software (Applied Biosystems, Italy). The primer sequences, restriction enzymes and PCR conditions used in the study are shown in Additional file 1: Table S1.

Statistical analysis
Allele and genotype frequencies for each polymorphism were calculated and tested as to whether they were distributed according to the Hardy-Weinberg equilibrium. A chi-square test for deviation from Hardy-Weinberg equilibrium was used to estimate differences in allele frequencies. The association of each polymorphism and clinical-pathological features of the patients was assessed by means of a chi-square test. A univariate logistic regression model was used to assess the effect of the same variables, included as dummy variables on incidence of toxicity (0-1-2 grade vs. [3][4], expressing results as odds ratios (OR) and relative 95% confidence intervals (95% CIs). Disease free survival (DFS) was defined as the time from the treatment start up to the date of first progression or death from any cause, whichever came first. Patients who had not died or had disease progression at the date of analysis were censored at the last available information on status. Overall survival (OS) was defined as the time from the treatment start to the date of death from any cause. Time-to-event data were described by the Kaplan-Meier curves. Cox proportional hazards models were used for univariate and multivariate analyses to estimate and test clinical-pathological features and polymorphisms for their associations with DFS and OS. Variables statistically significant at univariate analysis (at a level of p < 0.10) were included in the multivariate models. Results were expressed as hazard ratio (HRs) and their 95% CIs. Due to the explorative nature of the study, no adjustment of the significance level to make allowance for multiple tests has been made. Statistical significance was set at p < 0.05. All statistical analyses were carried out using SAS version 9.2 (SAS Institute, Cary, NC).

Patient characteristics
From June 2000 to September 2005 a total of 244 consecutive Caucasian patients with conservative surgery or mastectomy for primary BC, referred to the Breast Unit Surgical Department of the University of Perugia, Italy, were recruited. Histological diagnosis was confirmed by a pathologist at the Institute of Pathology, University of Perugia. The main clinical-pathological characteristics of the patients are summarized in Table 1.

Frequencies and associations among the polymorphisms and clinical-pathological features
The associations between genetic polymorphisms and the patient clinical-pathological features are reported in Additional file 2: Table S2.
The frequencies of genotypes GSTT1-null e GSTM1null were 20.5% and 54.1%, respectively and GSTM1-null allele was significantly higher in stage I than the GSTM1present allele (p = 0.042). The frequencies of the genotypes GSTP1 AA, AG, and GG were 59.4%, 39.3%, and 1.2%, respectively. GSTP1 AA genotype was significantly higher in stage III, in positive lymph nodes and in negative p53, than the GSTP1 AG or GG genotype (p = 0.006, p = 0.027 and p = 0.033, respectively). For MTHFR the frequencies of CC, CT, and TT were 27.5%, 47.5%, and 25.0%, respectively and the MTHFR CT or TT genotypes were significantly higher in stage III or in positive lymph nodes than the MTHFR CC genotype (p = 0.025 and p = 0.011, respectively). For the RFC1 polymorphism, the frequencies of GG, GA, and AA were 30.3%, 46.3%, and 23.4%, respectively. The frequencies of TS tandem repeat genotype distribution were 32.8% in 3R3R, 35.2% in 3R2R, and 32.0% in 2R2R. There was no statistically significant association among genotype distributions and tumor size, grading, ER, PgR, Ki67 and HER2 status. The genotype  Table 4). According to genotypes of GSTT1 and RFC1 genes we classified patients in three groups: the first with GSTT1-present and RFC1-AA (group1), the second with GSTT1-present and RFC1-GA/RFC1-GG or GSTT1-null and RFC1-GA/RFC1-AA (group2), and the third with GSTT1-null and RFC1-GG (group3).
Kaplan-Meier curves for DFS and OS are reported in Fig. 2e and f, respectively. At univariate analysis, confirmed at multivariate analysis (

Discussion
In the present study, we demonstrated that among BC patients who received CMF or FEC, those possessing the TS 2R/3R variant showed a significantly lower risk of severe toxicity (grade 3-4) for neutropenia and, when considering toxicity of any grade (G1-4), the same variant conferred a lower probability of developing oral mucositis. Our data are in agreement with previously published  studies [25][26][27] confirming a significant inverse association of TS 2R/3R polymorphism and severity toxicity. However, whereas in the study by Lecomte et al. patients with the 2R/2R genotype were 20 times more likely to have severe toxicity compared with 3R/3R carriers, this effect was much less pronounced in our study and more similar to the results of Schwab's study [28]. However, the role of other 5-FU catabolism-involved polymorphisms, such as dihydropyrimidine dehydrogenase (DPYD), should be explored to improve prediction of 5-FU toxicity [29]. At present, the real predictive value of MTHFR C677T polymorphism on MTX and 5-FU toxicity is not completely established. In our study, we found that the patients with MTHFR CT genotype had a higher probability of developing severe neutropenia than patients with MTHFR CC genotype. Some recent studies have shown increased toxicity in 677 T-carriers treated with methotrexate [30][31][32], although other studies did not confirm such an association [33,34]. Different methotrexate doses and schemes as well as diverse nutritional/folate status might account, at least in part, for these discrepant results. Probably, the heterozygous effects of MTHFR CT and TS 2R/3R genotypes as compared to each homozygous effect might be justified by considering that exogen factors, environmental conditions, dietary habits and lifestyle might play an important role [25-27, 35, 36]. No other significant differences in toxicity were found with respect to the other polymorphisms. There are a few studies on the role of GSTs isoenzymes on mortality in BC survivors drawn from community practice. The majority of these studies have small sample sizes, are based on participants diagnosed prior to 1999 and on women undergoing chemotherapy and/or radiotherapy. In addition, most of them examined only one GST gene (usually GSTP1). In our study, we showed that genotype GSTT1-null was associated with worse DFS and OS in EBC patients. This association was maintained in the multivariate model only for OS independently of age and other traditional predictors of prognosis. Our results are based on the assumption that the individuals with GSTT1-null genotype, that is associated with an absence of enzyme activity, are considered to be at increased risk for malignancies due to reduced efficiency in protection against environmental carcinogens [37,38]. Conversely, Ambrosone et al. [39], showed that Our results on GSTM1genotype are in agreement with those of Lizard-Nacol et al. [40] who, showed no effect of GSTM1-null genotype on DFS or OS among 92 women with advanced BC who had received cyclophosphamide, doxorubicin, and 5-FU. Whereas, Kristensen et al. [41] found that patients with GSTM1-null allele had a significantly shorter OS. Moreover, Yu Ke-Da et al. [42] showed a more complicated role for GSTM1 that should be considered in breast cancer risk prediction. The results of this study indicated a U-shaped association of GSTM1 with breast cancer, which challenges the linear genedosage effect of GSTM1 that was previously proposed. This effect was due to a new SNP, rs412543 (−498C > G) located in the promoter region that decreased gene transcription by 30-40% via reducing the DNA-binding affinity of AP-2. In contrast to these previous studies, our study is the only one to examine adjuvant therapy in a population of patients with a relatively uniform recurrence risk, with a longer follow-up (9.2 years), providing a homogeneous patient population in which to study treatment related genotypes and outcomes. Genetic background differences among races account for differences in the frequencies of allelic variants so that the association of polymorphic variants with a disease risk can significantly vary among populations. As far as we know, scanty information is available on the association of chemotherapeutic drug-related gene polymorphisms on toxicity and survival of breast cancer patients in non Caucasian populations. The results of Yang et al. showed no association between any of the GSTM1 or GSTT1 genotypes in patients with breast carcinoma who were treated with chemotherapy [43]. RFC1 genotypes, as predictors of BC treatment efficacy, have not been previously reported. Recent evidence suggests that G80A polymorphism in RFC1 is associated with altered folate/antifolate levels and may influence the efficacy of therapy with MTX [39]. Data suggest that subjects carrying the homozygous mutant AA genotype tend to have higher plasma folate and MTX levels and higher erythrocyte polyglutamate levels compared with those with the wild type or heterozygous genotype. In our study, for the first time to our knowledge, we showed that patients with RFC1 GG genotype had a shorter DFS and OS than carriers of the AA genotype. These observations are in keeping with previous studies on rheumatoid arthritis (RA). The work of Drozdzik et al. [44] showed that patients with RFC1 AA genotype responded to the therapy more effectively than carriers of AG and GG genotypes. The remission of RA symptoms was significantly higher (3.32-fold) in AA carriers in comparison to GG individuals. In contrast to RA patients, the study on acute lymphoblastic leukemia of Laverdiere et al. [45] showed children with AA genotype had worse prognoses than patients with GG genotype, and AA genotype was associated with higher plasma levels of MTX than other genotypes. Moreover, we showed, in an explorative analysis, that the combined genotypes (GSTT1-null/ RFC1-GG) had a negative prognostic effect on DFS and OS. This subgroup of tumors could have a more aggressive clinical course and the availability of a non-invasive, repeatable and reproducible technique to detect polymorphisms in the blood appears to be a useful tool for identifying high-risk BC patients. Therefore, further large sample size and well designed studies are greatly needed to confirm these preliminary results. Limitations of our study include relatively small sample size and low number of events, thus we were not able to evaluate the association with outcome by subgroups, such as menopausal status. Nevertheless, the association between GST polymorphisms and BC survival, showed by our results seems to be in agreement with those of the literature [39,40]. The cohort was established before some current treatments, such as aromatase inhibitors, and Her2/neu targeted therapies were available. Therefore, we cannot estimate what associations GST isoenzymes might have with survival in women using these treatments. However, our study has a larger sample size than most prior studies examining the association between GST polymorphisms and survival and it is the first study to evaluate RFC1 genotypes as predictors of BC treatment efficacy.