Genetic variants of DNA repair genes predict the survival of patients with esophageal squamous cell cancer receiving platinum-based adjuvant chemotherapy

Background Adjuvant chemotherapy in patients with resected esophageal squamous cell cancer (ESCC) remains controversial for its uncertain role in improving overall survival (OS). Nucleotide excision repair (NER) removes DNA-adducts in tumor cells induced by the platinum-based chemotherapy and thus may modulate efficacy of the treatment. The present study evaluated if single nucleotide polymorphisms (SNPs) of NER genes were prognostic biomarkers in ESCC patients treated with platinum-based adjuvant chemotherapy (PAC). Methods The analysis included 572 patients, for whom six SNPs of NER genes [i.e., XPC (rs1870134 and rs2228001), ERCC2/XPD rs238406 and ERCC5/XPG (rs2094258, rs2296147 and rs873601)] were detected with the TaqMan assay. Kaplan–Meier analyses and Cox proportional hazards models were used to evaluate their associations with disease free survival (DFS) and OS of these ESCC patients receiving PAC. Receiving operating characteristic curve analysis was used to evaluate the role of the risk genotypes in the DFS and OS. Results We found that ERCC5/XPG rs2094258 and rs873601 and ERCC2/XPD rs238406 SNPs were independently associated with poorer DFS and OS of ESCC patients [ERCC5/XPG rs2094258: CT+TT vs. CC: adjusted hazards ratio (adjHR) = 1.68 and P = 0.012 for DFS; adjHR = 1.99 and P = 0.0001 for OS; ERCC5/XPG rs873601: GA+GG vs. AA: adjHR = 1.59 and P = 0.024 for DFS; adjHR = 1.91 and P = 0.0005 for OS; ERCC2/XPD rs238406: TT vs. GG+GT: adjHR = 1.43 and P = 0.020 for DFS; adjHR = 1.52 and P = 0.008 for OS]. These HRs increased as the number of risk genotypes increased in the combined analysis. The model combining the risk genotypes with clinical characteristics or the TNM stage system was better in predicting outcomes in ESCC patients with PAC. Conclusion SNPs of ERCC2/XPD and ERCC5/XPG may independently and jointly predict survival of ESCC patients treated with PAC in this study population. Further validation in other study populations is warranted. Electronic supplementary material The online version of this article (doi:10.1186/s12967-016-0903-z) contains supplementary material, which is available to authorized users.


Background
Esophageal cancer (EC), more than 90 % of which are esophageal squamous cell carcinoma (ESCC), is the fourth leading cause of cancer-related deaths in China [1]. To date, surgery remains the standard treatment for resectable ESCC in China. But, for those patients who received an esophagectomy alone, their 5-year survival is still in a disappointing range of 15-40 % [2]. As a consequence, the surgery combined with adjuvant treatment has been employed to improve patients' survival.
In recent years, neoadjuvant chemotherapy or neoadjuvant chemoradiotherapy has been widely introduced as a standard regimen by various guides. Unfortunately, the results of some major multicenter prospective randomized controlled trials (MPRCTs) were controversial. Besides, most of these trials were conducted in western countries, and more than 50 % of the patients included in these trials were diagnosed with an adenocarcinoma. Shapiro et al. [3] reported some overall survival (OS) benefits when neoadjuvant chemoradiotherapy was added to surgery, while others all reported opposite conclusions [4][5][6]. Furthermore, the only one MPRCT conducted in Hong Kong Chinese patients with ESCC reported a negative conclusion [7]. Of the four MPRCTs employing neoadjuvant chemotherapy, two reported benefits in OS [8,9] but the other two did not [10,11].
Adjuvant chemotherapy is not taken seriously in clinical practice, because it is not thought as effective as it should be, and it may have impaired patients' functions as a result of esophagectomy and prolonged convalescence that hamper patients' timely administration of adjuvant therapy [12]. To date, only three MPRCTs about the efficacy of platinum-based adjuvant chemotherapy (PAC) for EC have been published, among which two were conducted in Japanese ESCC patients. Although one did not yield the expected results [13], the other larger trial found an enhanced 5-year disease free survival (DFS) for patients with lymph node metastasis [14]. Traditionally, the anatomic and pathologic staging has been the most commonly used prognostic factors in EC patients, but it did not provide sufficient information for evaluating the efficacy of PAC, because it did not account for host factors, because genetic variants may interact with PAC and thus play a role in determining clinical response and prognosis of the patients.
In China, PAC is still the preference of many oncologists to treat EC patients [15]. Platinum compounds produce DNA adducts by reacting with DNA to form both intrastrand and interstrand cross-links in tumor cells, mainly with the N 7 atom of guanine. These adducts result in a bulky distortion of the DNA helix, inhibit DNA replication, and eventually lead to cell death, if not repaired [16]. The amount of DNA adducts accumulated in tumor cells were correlated with the efficacy of platinum therapy and had an impact on clinical outcome of the patients [17].
Nucleotide excision repair (NER), which participates in two path ways of DNA repair, global genomic repair and transcription-coupled repair, plays an important role in detection and repair of DNA damage caused by UV, tobacco-related carcinogens and other carcinogenic chemicals [18,19], including DNA adducts formed by platinum [20]. Several xeroderma pigmentosum (XP) group genes are involved in NER, including group C (XPC), group D (ERCC2/XPD) and group G (ERCC5/XPG) [21][22][23][24]. Several studies suggested that single nucleotide polymorphisms (SNPs) in these three genes might be responsible for the variation in DNA repair capacity, leading to individual variation in cancer susceptibility and treatment response [25][26][27][28]. However, few reported studies have studied the SNPs' roles in the prognosis of EC patients treated with PAC.
Therefore, we hypothesize that potentially functional SNPs of the three NER genes may modulate prognosis of EC patients treated with PAC. In this study, we selected six well-studied potentially functional SNPs of the NER genes, including three from ERCC5/XPG (rs2094258 at 5′ near gene, rs22961475 at 5′ untranslated region (UTR) and rs873601 at 3′UTR), one from ERCC2/XPD (rs238406 at codon Arg156Arg) and two from XPC (rs1870134 at codon Val16Leu and rs2228001 at codon Lys939Gln) and evaluated their roles in survival of ethnic Han Chinese patients in eastern China.

Study population
The present study was done in a retrospective patient cohort in Fudan University Shanghai Cancer Center (Shanghai, China), and the research protocol was approved by the Institutional Ethics Review Board. Written informed consents were obtained from all patients before blood samples were obtained for genotype testing. Patients with perioperative mortality, defined as a death within 30 days of the operation or during the same hospitalization period, were excluded from the analysis. As a result, a cohort of 572 patients of ethnic Han Chinese in eastern China, who received an esophagectomy and had pathologically confirmed ESCC in the Department of Thoracic Surgery between March 2009 and December 2010, were included in the present study. Of these patients, additional 228 patients were excluded for the following reasons: 159 patients without undergoing postoperative chemotherapy for stage I disease, 35 patients without complete follow-up information, 7 patients for neoadjuvant chemotherapy and 27 patients for postoperative chemoradiotherapy. Therefore, the final analysis included 344 patients who completed four cycles of adjuvant chemotherapy (Oxaliplatin 135 mg/ m 2 d1 or cisplatin 40 mg/m 2 d1-3 plus 5-Fu 750 mg/m 2 d1-5).
Demographic and clinical information of the patients was abstracted from the medical records. All patients were staged according to the 7th edition of the American Joint Committee on Cancer staging system. Survival data were obtained through the follow-up in outpatient clinics or by telephone calls quarterly upto Oct. 31, 2014. The DFS was defined as the time interval between the date of surgical resection and the first confirmed detection of local recurrence or the appearance of new metastases. The OS duration of a patient was defined as the interval between surgical resection and the date of the latest follow-up or the death of the patients from any cause.

SNP selection and genotyping
We selected six potentially functional SNPs from the NCBI dbSNP database (http://www.ncbi.nlm.nih. gov/) and the SNPinfo (http://snpinfo.niehs.nih.gov/). Genomic DNA was extracted from the buffy-coat fraction of the blood samples using the Qiagen Blood DNA Mini Kit (Qiagen Inc., Valencia, CA). All the six SNPs were genotyped using the Taqman real-time PCR method with a 7900 HT sequence detector system (Applied Biosystems, Foster City, CA). The primers used in genotyping for these SNPs are listed in Additional file 1: Table S1. To ensure high genotyping accuracy, strict quality control procedures were implemented, and four duplicated positive controls and four negative controls (no DNA) were used in each of 384-well plates. Approximately 5 % of the samples were repeatedly genotyped, and the results were 100 % concordant.

Statistical methods
Cox proportional hazards regression analysis was used to evaluate the effect of genotypes and clinicopathological variables on patients' DFS and OS, calculated as hazards ratios (HRs) with their corresponding 95 % confidence intervals (CIs). Kaplan-Meier analysis was used to present the visual effects of clinicopathological and genetic variables on the cumulative probability of DFS and OS. Receiver operating characteristic (ROC) analysis was used to compare sensitivity and specificity of predicting overall survival by the parameters. Statistical significance of the improvement in area under the receiver operator characteristic curve (AUC) after adding an explanatory factor was calculated by Delong's test [29]. All reported P values were two-sided, and P < 0.05 was considered statistically significant. All analyses were performed using SAS software (version 9.2; SAS Institute, Cary, NC).

Demographics and clinicopathological characteristics of ESCC patients and their associations with DFS and OS
The final analysis included 344 ESCC patients who received esophagectomy and PAC (Table 1) and had complete data on demographics, clinical characteristics, genotyping, DFS and OS. These patients aged between 37 and 77 years at the time of diagnosis with a mean of 58.43 years and a standard deviation of 8.03 years. More patients were men than women (85.8 vs. 14.2 %), with 33.7 % of stage II and 66.3 % of stage III diseases, among whom 81.7 and 18.3 % underwent radical operation through two-field and three-field lymphadenectomy, respectively. The median follow-up time was 36.13 months, during which 196 (57.0 %) patients died at the last follow-up. In multivariate analysis, three variables, i.e., TNM stage [adjusted hazards ratio (adjHR) = 1.55 and 95 % CI 1.13-2.12 for III vs. II], vessel invasion (adjHR = 1.44 and 95 % CI 1.07-1.94 for yes vs. no), and lymphadenectomy (adjHR = 1.42 and 95 % CI 1.03-1.97 for three fields vs. two fields), were significantly associated with DFS (P < 0.05). The two variables, TNM stage (adjHR = 1.49 and 95 % CI 1.06-2.08 for III vs. II) and vessel invasion (adjHR = 1.56 and 95 % CI 1.14-2.12 for yes vs. no) remained to be independent prognostic factors for OS (P < 0.05), but smoking (adjHR = 1.44 and 95 % CI 1.01-2.07 for yes vs. no) instead of lymphadenectomy became the third independent prognostic factor for OS.
To evaluate the collective effect of the significant SNPs on patients' DFS and OS, we combined the risk genotypes of ERCC5/XPG rs2094258CT/TT and rs873601 GA/GG and ERCC2/XPD rs238406 TT for DFS and OS into a genotype score as the number of risk genotypes    Fig. 1a). Similarly, with the increasing NRG, risk of death increased correspondingly (P trend < 0.0001) ( Table 2; Fig. 1c).

Stratified analysis between the risk genotypes and survival of ESCC patients
We performed stratified analysis to assess whether the combined effect of risk genotypes (HG vs. LG) on DFS and OS was modified by some important demographic and clinicopathological factors listed in Table 1. For DFS or OS, we found that ESCC patients tended to exhibit an increased risk for disease progression or death in the subgroups with younger age (<60), of male, with smoking and drinking history, with a relatively earlier stage, without vessel and neural invasion and with two-field lymphadenectomy (P < 0.05) ( Table 3 in multivariate analysis and Additional file 2: Table S2 in univariate  analysis).

Discussion
In this study, we reported that some SNPs of the NER genes, such as ERCC5/XPG rs2094258 and rs873601 and ERCC2/XPD rs238406, may independently or jointly influence the prognosis of ESCC patients treated with PAC in eastern China. These genetic variants or genotypes, combined with some demographic and clinicopathological factors, once validated by others, may provide an improved prognostic tool for ESCC patients treated with PAC (Additional file 3: Table S3).
In the present study, we found that ERCC5/XPG rs2094258 CT/TT genotypes were associated with a decreased DFS and OS in ESCC patients treated with PAC. Although one previous study of 84 patients with squamous cell carcinomas did not find an association between the ERCC5/XPG rs2094258 SNP and response to PAC in non-small-cell lung cancer (NSCLC) [30], another study of 433 patients with advanced NSCLC did find an association between the ERCC5/XPG rs2094258 SNP and outcome of PAC, with the conclusion consistent with the present study [31]. The rs2094258 SNP is located in the 5′ UTR of ERCC5/XPG, which is a putative transcription factor binding site. Although genetic variants in gene promoters may alter gene expression levels and thus likely exert some influence on clinic outcome [25,32,33], there has been no report about biological or functional validation for this polymorphic site, which warrants additional mechanistic studies.
In the present study, we found that the ERCC2/XPD rs238406 TT genotype was associated with a reduced DFS and OS in ESCC patients, but few studies have reported its role in prognosis of cancer patients. One Taiwan study found that ERCC2/XPD rs238406 CC (or GG of its antisense) instead of the AA (or TT of its antisense) genotype of 185 ESCC patients with neoadjuvant chemoirradiation followed by esophagectomy could additively increase risk of death and disease progression in cisplatin-based neoadjuvant concurrent chemoradiation therapy [34]. Contrary to their results, the TT genotype was associated with the worse DFS and OS in the present study. Another study showed that rs238406 AA carriers had less efficiency of DNA adduct formation We also found that patients with the ERCC5/XPG rs873601 GA/GG genotypes had an increased risk of progression and death of ESCC after PAC, which was not reported before. One study reported that the rs873601 G allele was associated with better PFS and OS of patients with advanced NSCLC, which may be disease-specific and need to be validated in future studies [37].
It is likely that the effect of a single SNP on clinical outcome may be much restricted, but the combined effect of several SNPs in the same or different genes could be much greater. Indeed, we found that the collective effect of the risk genotypes identified in the present study better predicted DFS and OS of the patients. Compared with some clinical factors, genetic variants such as SNPs may have a weaker effect on prognosis. As shown by ROC curves in the present study, SNPs had almost the same effect on prognosis as the TNM stage, although it served a relatively inferior role in prognosis, compared with other clinical characteristics that included clinical characteristics such as smoking, vessel invasion, neural invasion, TNM stage and lymphadenectomy.
In the stratified analyses, we found that the genotypesurvival association was more evident for a mild status of clinical characteristics, such as without vessel and neural invasion, II stage, two-field lymphadenectomy, which is consistent with what were reported in Lee's study for EC with neoadjuvant chemoradiation and esophagectomy [34]. It is likely that the severe effects of a later TNM stage and vessel invasion as poor prognosis factors [38][39][40] on the survival may have masked those benefit form genetic factors.
There were some limitations in the present study. First, patients included in the analyses were from one hospital in eastern China, which may not represent the general population. Second, only six putatively functional SNPs of three NER genes were tested in the study, and there were other genetic variants in these genes or other NER pathway genes that may affect the prognosis of ESCC. Finally, the present study was retrospective instead of a prospective or randomized design, thus the bias caused by other possible factors, such as standardization of dose and the judgment of disease progress, could not have been completely excluded.

Conclusions
In summary, we identified that ERCC5/XPG rs2094258 CT/TT and rs873601 GA/GG and ERCC2/XPD rs238406 TT genotypes may independently or jointly affect survival of ESCC patients treated with PAC. These findings, once validated in future prospective studies with large sample sizes and better study designs, will provide some promising guidance for personalized treatment for ESCC patients in the adjuvant setting in China.