Abstract

Background. There is no definitive conclusion so far on the predictive values of ERCC1 polymorphisms for clinical outcomes of platinum-based chemotherapy in non-small cell lung cancer (NSCLC). We updated this meta-analysis with an expectation to obtain some statistical advancement on this issue. Methods. Relevant studies were identified by searching MEDLINE, EMBASE databases from inception to April 2015. Primary outcomes included objective response rate (ORR), progression-free survival (PFS), and overall survival (OS). All analyses were performed using the Review Manager version 5.3 and the Stata version 12.0. Results. A total of 33 studies including 5373 patients were identified. ERCC1 C118T and C8092A could predict both ORR and OS for platinum-based chemotherapy in Asian NSCLC patients (CT + TT versus CC, ORR: OR = 0.80, 95% CI = 0.67–0.94; OS: HR = 1.24, 95% CI = 1.01–1.53) (CA + AA versus CC, ORR: OR = 0.76, 95% CI = 0.60–0.96; OS: HR = 1.37, 95% CI = 1.06–1.75). Conclusions. Current evidence strongly indicated the prospect of ERCC1 C118T and C8092A as predictive biomarkers for platinum-based chemotherapy in Asian NSCLC patients. However, the results should be interpreted with caution and large prospective studies are still required to further investigate these findings.

1. Introduction

Lung cancer is currently the most common human malignancies and the leading cause of cancer-related mortality in the world. Non-small cell lung cancer (NSCLC) accounts for ~85% of all lung cancer cases, with 30% being squamous cell carcinoma (SCC) and the remaining (~70%) being collectively classified as non-SCC [1]. Despite efforts for early detection, approximately two-thirds of NSCLC are diagnosed at advanced stages with limited surgery options. Recent advances in target therapy in NSCLC have added more choices for non-SCC, but chemotherapy still remains the therapeutic mainstay for SCC and some non-SCC. Platinum combined with the third-generation cytotoxic drugs including gemcitabine and pemetrexed is the standard first-line treatment for advanced NSCLC [2]. However, the chemotherapy response and clinical prognosis for NSCLC patients vary remarkably among individuals. Some patients show significant tumor regression while others may develop intrinsic or acquired resistance to chemical drugs, which highlights the need to tailor the treatment for individuals [3]. Pharmacogenomics has been perceived as a useful tool to predict treatment response as well as clinical prognosis in cancer patients, and numerous studies have suggested the functions of single nucleotide polymorphisms (SNPs) in affecting drug sensitivity by modifying relevant genes [4]. Consequently, it is necessary to identify functional gene polymorphisms as new biomarkers for accurately predicting clinical outcomes of chemotherapy in NSCLC, which could also reduce the side effects and provide the most cost-effective approach for individuals [3, 4].

The cytotoxic mechanism of platinum is that the agents bind to DNA and form platinum-DNA adducts, which consequently block DNA replication and transcription, resulting in cellular apoptosis and growth inhibition. The damaged DNA system could be repaired by many biological processes, leading to a resistance to the platinum. The nucleotide excision repair (NER) pathway plays dominant roles in the DNA repair process, which is responsible for the recognition of DNA damage and removal of the damaged nucleotides [5]. Excision repair cross-complementation group 1 (ERCC1) is a key rate-limiting enzyme in the NER pathway. It has been documented that high expression of ERCC1 could block the platinum efficacy and cause drug resistance [6]. Therefore, ERCC1 has become one of the most promising biomarkers for efficacy of platinum-based chemotherapy. Two most common SNPs of ERCC1 gene are C118T (rs11615) with a C → T substitution at exon 4 and C8092A (rs3212986) with a C→A change in the 3′-untranslated region, both of which have been speculated to regulate the ERCC1 mRNA expression and manage the NER pathway, consequently affecting the efficacy of platinum-based chemotherapy [7].

To date, lots of studies regarding gene variants that may alter response and tolerability to chemotherapy drugs have been carried out, including both retrospective and perspective studies. The ultimate purpose of a pharmacogenomics analysis is to evaluate the potential value of individual’s gene variants to predict the drug efficacy, thus allowing for optimization and personalization of the clinical decision making. However, such an approach seems still a distant goal. The previous results on the relationships between ERCC1 polymorphisms and clinical outcomes of platinum-based chemotherapy in NSCLC are inconclusive and even conflicting to each other, probably due to the complex gene interaction, environmental effects, different detection methods, sample sizes, or study designs. Therefore, we performed a systematic review and meta-analysis of published studies to comprehensively evaluate the predictive values of two promising ERCC1 SNPs for clinical outcomes of platinum-based chemotherapy in NSCLC, with an expectation to provide useful evidence and suggestions for clinical practice and future investigation.

2. Materials and Methods

2.1. Literature Search

We searched for original studies investigating the associations between ERCC1 SNPs and platinum-based chemotherapy in NSCLC by using MEDLINE and EMBASE databases. The following terms were combined variously for search: “ERCC1,” “platinum or cisplatin or carboplatin,” “polymorphism or variant or mutant,” and “lung cancer.” The literature search was last updated on April 2015. References from retrieved articles and previous meta-analysis were further screened manually for additional qualified studies.

2.2. Inclusion Criteria

The inclusion criteria were as follows: (1) patients should be pathologically confirmed to have NSCLC; (2) studies should assess the relationships between ERCC1 polymorphisms and clinical outcomes of platinum-based chemotherapy; (3) polymorphisms should be genotyped: ERCC1 C118T (rs11615), C8092A (rs3212986); and (4) clinical outcomes including the objective response rate (ORR), overall survival (OS), and progression-free survival (PFS) should be reported. The most recent and comprehensive data were included if there were duplications.

2.3. Quality Assessment

The quality of each study was assessed independently by two reviewers (Han and Liu) using the Newcastle-Ottawa Scale (NOS), and disputations were settled by a third reviewer (Sun). The NOS consists of 3 parameters for the quality of case-control/cohort study: selection, comparability, and exposure/outcome. Studies with NOS scores >6 were considered with high quality [8].

2.4. Data Extraction

Data extraction was performed independently by two investigators (Han and Liu) and disagreements were adjudicated by a third reviewer (Sun). For each study, general characteristics such as the first author, publication year, country and ethnicity of patients, sample size, tumor stages, follow-up time, chemotherapy drugs and treatment line, SNP allele frequency, and genotyping method were collected. Results data included odds ratio (OR) for the ORR, hazard ratios (HR), and 95% confidence interval (CI) for the OS and PFS. If HR and 95% CI were not available directly, the estimated value was derived indirectly using Tierney’s methods [9].

2.5. Statistic Analysis

The ORR was used to measure the efficacy of chemotherapy regimens. Patients were divided into the responder group including the complete responders (CR) or partial responders (PR) and the nonresponders group including the stable disease (SD) or progressive disease (PD) according to the Response Evaluation Criteria in Solid Tumors criteria [10]. The pooled OR and 95% CI were calculated by (CR + PR)/(SD + PD). The OS and PFS on behalf of prognosis of NSCLC patients were evaluated by calculating pooled Cox proportional HR and 95% CI as relevant effect measures. The associations between SNPs and ORR/OS/PFS were examined by both the dominant (heterozygote or homozygote variant versus wild type) and codominant models (heterozygote variant versus wild type, homozygote variant versus wild type). Heterogeneity among studies was tested by the Chi-square-based test and statistics. Values of for the test or % were considered as lack of heterogeneity and a fixed-effect model was used; otherwise, a random-effect model was used. Subgroup analysis by ethnicity was used to detect and reduce potential source of heterogeneity among studies. The publication bias was investigated by the inverted funnel plot, Begg’s test, and Egger’s test. Values of were indicative of statistically significant publication bias. All values were two sided. All analyses were performed using the Review Manager version 5.3 (Oxford, England) and the Stata version 12.0 (Stata Corporation, College Station, TX).

3. Results

3.1. Study Identification

As shown in the flow chart for the study selection process (Figure 1), 126 publications were initially retrieved from the PUBMED and EMBASE databases. After screening by checking title and abstract, 57 articles being reviews, meta-analyses, or molecular experiments and other irrelevant studies were excluded. Then, we checked full text for the remaining 69 studies, and subsequently 36 studies were excluded for the following reasons: evaluating the association between gene expression levels and clinical outcome, including patients with SCLC or other tumors, including patients treated with chemotherapy regimens without platinum, conducting the genotype conjoint analysis, and reporting no detailed or extractable survival data. Finally, a total of 33 studies including 5373 patients were identified to be eligible for this meta-analysis. All of our included studies had a high quality with NOS scores >6. The baseline characteristics of the included studies were shown in Table 1 and all the results were summarized in Table 2.

3.2. ERCC1 C118T
3.2.1. Objective Response Rate

Twenty-three studies including 3272 patients were eligible for this analysis [1133]. Overall, the T allele showed a significant association with a poorer ORR (CT + TT versus CC: OR = 0.82, 95% CI = 0.71–0.96, %, and ) (Figure 2(a)). For the ethnicity-specific subgroup analysis, the T allele was found to be correlated with a decreased ORR in Asian patients, but not in Caucasian patients (CT + TT versus CC: for Asians, OR = 0.80, 95% CI = 0.67–0.94, %, and ; for Caucasians: OR = 0.93, 95% CI = 0.66–1.30, %, and ) (Figure 2(a)). Similar results were found by the codominant model (CT versus CC: for total, OR = 0.73, 95% CI = 0.61–0.89, %, and ; for Asians, OR = 0.69, 95% CI = 0.55–0.85, %, and ; for Caucasians, OR = 0.93, 95% CI = 0.62–1.38, %, and ) (TT versus CC: for total, OR = 0.63, 95% CI = 0.47–0.85, %, and ; for Asians, OR = 0.64, 95% CI = 0.44–0.93, %, and ; for Caucasians, OR = 0.63, 95% CI = 0.39–1.01, %, and ) (Table 2).

3.2.2. Overall Survival

Seventeen studies including 2926 patients were available for this evaluation [11, 14, 15, 21, 22, 25, 27, 3241]. Generally, patients with CT or TT genotype had a poorer OS when compared to those with CC genotype (CT + TT versus CC: HR = 1.25, 95% CI = 1.05–1.50). However, significant heterogeneity (%, ) was found in this analysis and the random effect model was applied (Figure 2(b)). In the ethnic subgroup analysis, CT or TT genotype was found to be associated with an unfavorable OS in Asians but not in Caucasians (CT + TT versus CC: for Asians, HR = 1.24, 95% CI = 1.01–1.53, %, and ; for Caucasians, HR = 1.27, 95% CI = 0.88–1.85, %, and ) (Figure 2(b)). Obvious heterogeneity still existed and sensitivity analysis indicated no study to be deleted. In the codominant model analysis with only Asians, the homozygote showed a more significant trend to be associated with survival than the heterozygote (CT versus CC: HR = 1.09, 95% CI = 0.86–1.38, %, and ) (TT versus CC: HR = 1.49, 95% CI = 1.18–1.88, %, and ) (Table 2).

3.2.3. Progression-Free Survival

Nine studies with 943 patients were included in this analysis [11, 14, 19, 25, 27, 29, 36, 40, 41]. There was no evidence to support an association between the ERCC1 C118T polymorphism and PFS in the dominant model (CT + TT versus CC: HR = 1.11, 95% CI = 0.84–1.47, %, and ) (Figure 2(c)). The ethnic subgroup analysis also demonstrated no relationship between the ERCC1 C118T polymorphism and PFS in Asians or Caucasians (CT + TT versus CC: for Asians, HR = 1.15, 95% CI = 0.86–1.54, %, and ; for Caucasians, HR = 1.06; 95% CI = 0.64–1.75, %, and ) (Figure 2(c)). Only one study conducted this evaluation in the codominant model (CT versus CC: HR = 1.47, 95% CI = 0.91–2.36; TT versus CC: HR = 1.86, 95% CI = 1.13–3.07), and the homozygote showed a more obvious trend to be associated with PFS than the heterozygote (Table 2).

3.3. ERCC1 C8092A
3.3.1. Objective Response Rate

The data from eleven studies with 1607 patients were included in this evaluation [15, 16, 19, 24, 28, 29, 3133, 42, 43]. In total, no statistically significant association was found between ERCC1 C8092A and ORR (CA + AA versus CC: OR = 0.87, 95% CI = 0.71–1.07, %, and ) (Figure 3(a)). Stratified analysis by ethnicity showed an obvious association in Asians, but not in Caucasians (CA + AA versus CC: for Asians, OR = 0.76, 95% CI = 0.60–0.96, %, and ; for Caucasians, OR = 1.51, 95% CI = 0.95–2.39, %, and ) (Figure 3(a)). Positive results were also found by the codominant model in Asian patients (CA versus CC: OR = 0.71, 95% CI = 0.53–0.96, %, and ) (AA versus CC: OR = 0.42, 95% CI = 0.26–0.70, %, and ) (Table 2).

3.3.2. Overall Survival and Progression-Free Survival

Seven studies consisting of 1078 patients were selected for this analysis [15, 24, 3234, 37, 42]. Overall, no correlation was observed between ERCC1 C8092A and survival (CA + AA versus CC: HR = 1.18, 95% CI = 0.77–1.79), and high level of heterogeneity was also observed in this analysis (%, ) (Figure 3(b)). After excluding the only one study for Caucasians, A allele was observed to be associated with a worse survival in the remaining four Asian studies (HR = 1.37, 95% CI = 1.06–1.75, %, and ). Significant association was also found by the codominant model in Asian patients (CA versus CC: HR = 1.38, 95% CI = 1.03–1.86, %, and ) (AA versus CC: HR = 2.03, 95% CI = 1.19–3.47, %, and ) (Table 2). Only one Chinese study with 90 patients examined the relationship between ERCC1 C8092A and PFS, and no association was observed (HR = 1.08, 95% CI = 0.72–1.63) (Table 2).

3.4. Publication Bias

The publication bias was performed in the evaluation of the relationship between ERCC1 C118T and ORR by the dominant model, since this comparison included the most studies (). The funnel plot revealed no obvious publication bias, with a symmetrical distribution of study results around the pooled measurement of effect. No evidence of publication bias was detected according to Begg’s test () and Egger’s test ().

4. Discussion

The optimal treatment for cancer patients has been evolving for decades. Lots of studies have demonstrated the significant association between gene characteristics and chemotherapy efficacy in individuals. With respect to the gene characteristics, gene polymorphisms are considered more accessible for detection and valuable for evaluation than gene mRNA levels. Single nucleotide polymorphisms (SNPs) have been widely suggested to be capable of affecting the drug sensitivity by modifying functions of relevant genes, thus being predictive biomarkers for treatment efficacy in cancer patients. ERCC1 has been well known to be a promising screening tool in the application of platinum-based chemotherapy in NSCLC. Numerous studies have been carried out to address this question and quite a few meta-analyses evaluating this question have been published. However, to date, no definitive conclusion has been reached yet. The latest meta-analysis was conducted by Yang and Xian who did the literature search to July 31, 2013 [44]. We updated this meta-analysis to evaluate the issue more comprehensively with an expectation to get some statistical advancement. We hope our efforts to select gene biomarkers could make a positive contribution to the progress of optimization and personalization of chemotherapy in NSCLC.

Yang and Xian demonstrated merely a trend of ERCC1 C118T in predicting ORR with no statistical significance (CT + TT versus CC: OR = 0.94, 95% CI = 0.72–1.23), and stratified analysis by ethnicity also failed to show any valuable results. Besides, they got an opposite trend when analyzing the relevance of C8092A with ORR by the dominant model, probably due to the heterogeneity resulting from ethnicity (CA + AA versus CC: for total: OR = 1.05, 95% CI = 0.83–1.32; for Asians: OR = 0.97, 95% CI = 0.75–1.25; for Caucasians: OR = 1.55, 95% CI = 0.87–2.74). The assessment by the codominant model did not reveal any significant results either. Notably, our updated analysis with more data exactly showed a significant relevance of T allele and A allele with a worse ORR compared with C allele in the Asian subgroup, by both the dominant and codominant models. Our statistically significant results highlight a good prospect of these two ERCC1 SNPs as predictive biomarkers for the response of platinum-based chemotherapy in NSCLC.

Furthermore, when assessing the predictive values of ERCC1 C118T/C8092A for prognosis of NSCLC patients receiving platinum-based chemotherapy, we also obtained some statistical advancements. Yang and Xian observed an association between C118T and OS in totals and in Asians (CT + TT versus CC, for total: HR = 1.26, 95% CI = 1.02–1.55, %, and ; for Asians: HR = 1.35, 95% CI = 1.04–1.75, %, and ; for Caucasians: HR = 1.12, 95% CI = 0.74–1.68, %, and ), while no apparent correlation between C8092A and OS in both ethnic groups was observed (CA + AA versus CC, for total: HR = 1.26, 95% CI = 0.81–1.95, %, and ; for Asians: HR = 1.50, 95% CI = 0.89–2.55, %, and ; for Caucasians: HR = 0.87, 95% CI = 0.30–2.56, %, and ). The great heterogeneity in Yang’s analysis to some extent impacted the credibility of their results. Importantly, in our analysis for C8092A, A allele exhibited a significant association with a poorer OS in the Asian subgroup analysis with smaller heterogeneity (CA + AA versus CC, for Asians: HR = 1.37, 95% CI = 1.06–1.75, %, and ), and consistent results were observed by the codominant model, which indicated a possibility that C8092A may serve as a predictor for clinical prognosis of platinum-based chemotherapy in NSCLC. Nevertheless, in the C118T analysis, we revealed similar results with Yang’s analysis and significant heterogeneity still existed by the dominant model. Notably, in our evaluation by the codominant model which only included Asian patients, the homozygote showed obvious association with an unfavorable OS while the heterozygote only showed a slight trend without statistical significance. This phenomenon presents a possibility that the number of T alleles may be inversely proportional to the clinical outcome of platinum-based chemotherapy.

Fewer studies reported data of PFS, and most of them performed analysis by the dominant model. Combining the extracted data, we could not identify a predictive role of ERCC1 C118T/C8092A for PFS of platinum-based chemotherapy in NSCLC. However, we also observed a slight trend that the T allele may correlate with a poorer PFS, as the HRs in various comparisons of mutant genotype with wild genotype all showed greater than 1 (CT + TT versus CC, HR = 1.11; CT versus CC, HR = 1.47; and TT versus CC, HR = 1.86). Only one study reported the relevance of C8092A with PFS, which also revealed no obvious association. In consideration of the positive results obtained in the evaluation for ORR and OS, there still exists the possibility that minor variants in ERCC1 SNPs can as well predict PFS. More studies in the future are expected to confirm this possibility.

In particular consideration of the influence of race on chemosensitivity, we conducted the stratified analysis by ethnicity. Notably, the results revealed a significant distinction between the Asian and the Caucasian subgroups, with ERCC1 SNPs predicting both ORR and OS prominently in Asian patients but not in Caucasian patients. The significant ethnic discrepancy in the predictive roles of ERCC1 SNPs may be due to the gene-gene and gene-environment interactions, different genetic background, and lifestyle, which is not yet clearly defined.

Our work has several strengths. First, we comprehensively evaluated gene effects on the standard first-line treatment in NSCLC, platinum-based chemotherapy, and the indicators covered both treatment efficacy and clinical prognosis (ORR, OS, and PFS); second, we made a statistical progress when evaluating the predictive values of ERCC1 polymorphisms for platinum-based chemotherapy in NSCLC; third, all of our included studies were of high quality and no evident publication bias was observed. All these points significantly increased the statistical power of our analysis.

Despite the strengths of our work, some limitations should be taken into consideration. Above all, an obvious heterogeneity was observed in the analysis of the ERCC1 C118T for OS and PFS and still existed after stratified analysis by ethnicity. It is speculated that the heterogeneity mainly derived from heterogeneous samples and various clinical trial designs: firstly, the baseline characteristics of populations in different trials were impossible to coincide with each other, such as age, sex, histology type, and tumor stage; secondly, clinical trial designs varied in many aspects, such as the coadministration of other drugs, drug administration mode, treatment cycle, treatment line with or without surgery or radiotherapy, and genotyping method. Most of the included studies employed the platinum-based chemotherapy as a first-line treatment; only a few studies did not specialize the treatment line. And the concurrent drugs combined in the platinum-based regimen also to some extent affected the clinical outcome, as the efficacy of gemcitabin and pemetrexed also varied among individuals and could be influenced by certain gene polymorphisms as well, such as ribonucleotide reductase M1(RRM1) and thymidylate synthase (TS). All of these potential discrepancies interfere with the standardization of data acquisition; however, we were unable to conduct a further stratified analysis since most studies did not provide detailed information about these factors. In addition, there inevitably existed some inaccuracies in our data extraction and analysis: firstly, some researches in the literature did not report concrete time-to-event data so that we had to obtain calculable data indirectly by reading Kaplan–Meier curves or conducting data transformation; secondly, most studies reported unadjusted estimates and only a few studies are supplied with adjusted estimates; moreover, those estimates were not always adjusted for the same potential confounders.

We still have a long way to go to ultimately turn to ERCC1 polymorphisms for our clinical decision making, since various factors contribute to the individual variation in drug response. Polymorphism of genes encoding proteins involved in the transport, metabolism, and action of drugs influences the clinical outcome of chemotherapies prominently, and other factors such as age, sex, physical condition, hepatic and renal function, tolerance to toxicity, and alcohol and tobacco use also affect the drug efficacy. All the factors should be taken into account comprehensively when making clinical decisions. Although it is a distant goal, our work has made an important step along this path, indicating a necessity to conduct future prospective studies with large sample sizes and better study designs to validate these conclusions and prove the feasibility of the customization approach to our clinical decision making.

5. Conclusion

In conclusion, the prospect of optimal chemotherapy in NSCLC based on validated biomarkers was further confirmed with the latest data. We provided statistical evidence that both ERCC1 C118T and C8092A could be useful predictive biomarkers for platinum-based chemotherapy in Asian NSCLC patients. However, considering the limitations, heterogeneity and biases existed within our analysis, our conclusions should be interpreted with caution, and large prospective studies are still required to further validate these findings.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work was supported by grants from the National Natural Science Foundation of China (81372334).