Abstract
The prognosis for lung cancer patients treated with chemotherapy is poor. Single nucleotide polymorphisms (SNPs) in matrix metalloproteinase (MMP) genes could influence treatment outcome by altering apoptotic pathways. Eight SNPs with known or suspected phenotypic effect in six genes (MMP1, MMP2, MMP3, MMP7, MMP9 and MMP12) were investigated.
For 349 Caucasian patients with primary lung cancer, receiving first-line chemotherapy, three different endpoints were analysed: response after the second cycle, progression free survival (PFS) and overall survival (OS).
The prognostic value of the SNPs was analysed using multiple logistic regression for all patients and histology-, stage- and treatment-specific subgroups. Hazard ratio estimates for PFS and OS were calculated using Cox regression methods.
None of the investigated polymorphisms modified response significantly in the whole patient population.
However, tumour stage IIIB variant allele carriers of MMP2 C-735T showed a significantly worse response. PFS was significantly prolonged in MMP1 G-1607GG variant allele carriers and OS in small cell lung cancer patients carrying the MMP12 A-82G variant allele.
In conclusion, this study identified SNPs in MMP1, MMP2, MMP7 and MMP12 for further investigation as possible predictors of chemotherapy outcome in lung cancer patients.
With about 1.35 million new cases diagnosed per year, lung cancer is the most common cancer worldwide 1. It is a disease with major morbidity and continuingly bad prognosis. While early stage nonsmall cell lung cancer (NSCLC) can be treated by surgery, late stage NSCLC and small cell lung cancer (SCLC) cases receive chemotherapy as the treatment of choice 2. Polymorphisms, which reduce or inhibit apoptosis, can cause chemotherapy resistance 3. Cohort studies are used in order to elucidate correlations between biomarkers in the host genome, shown to be relevant for patient outcome, and therapy response.
Matrix metalloproteinase (MMP) expression is associated with the development of an extensive list of diseases especially various malignant tumours. Their involvement in promotion of metastasis, chemotherapy resistance of cancer and a bad treatment outcome has been shown in previous studies 4, 5.
High expression of MMPs is usually associated with all steps of cancer initiation and progression. MMP3, MMP7 and MMP9, have also been reported to influence the Fas/FasL-mediated extrinsic, as well as the p53/PKC mediated intrinsic apoptotic pathway 6, 7. By cleavage of plasminogen and collagen XVIII, MMP12 is one of the most effective producers of the angiogenesis inhibitors angiostatin and endostatin 8. In addition, it has an influence on the plasmin levels and subsequent activation of MMPs like MMP3 and MMP9 9. MMP2, also known as gelatinase A, is involved in migration and invasion processes and also seems to have influence on chemotherapy response. It has been shown in rat models that platinum-based chemotherapy has higher response rates when the animals are co-treated with prinomastat, a specific MMP2 inhibitor 10. MMP1 is the most highly expressed interstitial collagenase, degrading fibrillar collagen. Overexpression of MMP1 shown in tumour tissues has been suggested to be associated with tumour invasion and metastasis 11.
Polymorphisms in MMP genes have been shown to influence the expression pattern of the protein. One polymorphism in the promoter region of MMP2 (C-735T) was reported to lead to a lower expression of the protein due to a disruption of a SP1-binding site 12. As far as MMP3 is concerned, the variant allele of the ins/del polymorphism (6A-1171 5A) is associated with a higher transcription rate 13. For MMP7, two single nucleotide polymorphisms (SNPs) (C-153T and A-181G) have been functionally characterised and, in both cases, the variant alleles lead to an enhanced expression 14. Another SNP in the promoter of the MMP9 gene (C-1562T) increases transcription by the disruption of a repressor binding site 12. Polymorphisms in the MMP12 gene (A-82G; A1082G) also have been shown to influence the outcome of cancer patients. The promoter polymorphism is responsible for lowering the transcription by disruption of an AP1-binding site 15, whereas the polymorphism A1082G has not been shown to have a phenotypic effect. The GG allele of the functional MMP1 G-1607GG polymorphism is associated with higher expression levels of MMP1 and with increased susceptibility to head and neck and lung cancer 16, 17.
This study focuses on the relationship between polymorphisms in MMP genes and three defined endpoints of clinical outcome: 1) overall response rate (ORR) to chemotherapy after the second cycle; 2) overall survival (OS); and 3) progression free survival (PFS). We hypothesised that genetic background, such as polymorphisms in gDNA can influence therapy outcome. We also assumed that this background has a more immediate effect on an early endpoint, such as ORR, than the clinically more relevant late endpoints, OS and PFS, due to tumour-specific accumulation of genetic variations by clonal selection during many cycles of therapy.
Eight SNPs in six MMP genes were analysed in a cohort of 349 lung cancer patients, consisting of 187 NSCLC, 161 SCLC and one patient with a mixed histology, receiving first-line chemotherapy.
PATIENTS AND METHODS
Study cohort
349 patients of Caucasian origin with histologically confirmed primary lung cancer eligible for first-line chemotherapy were recruited between March 1999 and October 2004 at the Thoraxklinik in Heidelberg. The cut-off date for the follow-up was March 2005. All patients had never received antineoplastic chemotherapy nor had they previously been diagnosed with another malignancy. Where possible, tumour stage at the time of diagnosis was determined according to the clinical tumour, node, metastasis (cTNM) of the Union Internationale contre le Cancer 18 using hospital records. For some SCLC patients, the tumour stage was classified as limited or extensive disease based on the Veterans Administration Lung Cancer Study Group criteria 19. For statistical analysis, limited and extensive disease were regarded as stage III and IV, respectively (table 1⇓).
Main clinical characteristics for the nonsmall cell lung cancer(NSCLC) and small cell lung cancer (SCLC) patient groups
Tumour response was assessed after the second cycle of first-line chemotherapy as complete remission (CR), partial remission (PR), stable disease (SD) or progressive disease (PD) according to RECIST (response evaluation criteria in solid tumours) 20. Progression free survival (PFS) is defined as the time interval (in days) between start of chemotherapy and documented progression (uncensored observation). In case no progression was documented, PFS was calculated until the last progression-free examination (censored observation) irrespective of whether that patient was lost to follow up or whether death occurred later. OS was defined as the time interval (in days) between start of chemotherapy and the documented date of death (uncensored observation) or, when date of death was unknown, the last date when the patient was still alive (censored observation).
60% of the 187 NSCLC cases died during the observation period, the median OS time was 410 days (95% CI 353–514). The median PFS time was 216 days (95% CI 189–251), and the median observation time was 812 days (95% CI 499–1120). 63% of the 161 SCLC patients died during observation, the median OS time was 405 days (95% CI 347–557). The median PFS time was 271 days (95% CI 239–314), and the median observation time was 777 days (95% CI 648–838). All participants signed informed consent forms and most study subjects completed a self-administered questionnaire with detailed information on pre-treatment, smoking habits and occupational exposure. The study was approved by the ethical committee of the Medical Faculty of the University of Heidelberg (Heidelberg, Germany) (Ref.-Nr.199/2001). Blood samples were collected prior to chemotherapy.
DNA extraction
Buffy coat from 5 mL venous blood in EDTA was archived at -80°C. Genomic DNA was isolated using either the QIAamp DNA blood midi kit (Qiagen, Hilden, Germany), or by an automated DNA extraction protocol on the MagNA Pure System (Roche, Basel, Switzerland) according to the manufacturer's instructions.
Genotyping
For the detection of the polymorphisms MMP3 6A-1171 5A, MMP12 A-82G and MMP12 A1082G, assays using fluorescence-based melting curve analysis (LightCycler480) were designed (shown in table 2⇓). Each PCR with a final volume of 10 μL was performed in 96-well plates (Roche, Mannheim, Germany) with ∼20 ng of genomic DNA as template. After 2 min of 95°C initial incubation, 45 subsequent cycles of 5 s at 95°C, 10 s for annealing at 56°C and 15 s at 72°C for elongation were carried out. Concentrations for the multiplex reaction for MMP12 A-82G and MMP12 A1082G genotyping were 1× PCR Buffer, 5 mM MgCl2, 200 μM dNTPs, 0.5 U Qiagen Polymerase, 0.7 μM forward primer (A-82G), 0.2 μM reverse primer (A-82G), 0.1 μM forward primer (A1082G), 0.7 μM reverse primer (A1082G) and for all four probes 0.2 μM each. Concentrations for MMP3 6A-1171 5A genotyping were 0.1 μM forward primer and 0.5 μM reverse primer, 2.5 mM MgCl2, 1× PCR buffer, 200 μM dNTPs, 0.5 U Taq Polymerase (Eurolone, England) and both probes 0.2 μM each. All primers and probes (TIB MOLBIOL, Berlin, Germany) are detailed in table 2⇓.
Primers and probes
For validation and quality control of the LightCycler480 method, PCR restriction fragment length polymorphism (RFLP) was performed for MMP12 A-82G 21. Another PCR-RFLP analysis, for MMP12 A1082G was newly developed and carried out in a total volume of 20 μL with 400 nM of each primer (forward primer 5′-GATGACAAATACTGGTTAATTAGGA-3′; reverse primer 5′-CTGGTTATCTACAAAGAAGT-3′), 200 μM dNTPs, 1× PCR Buffer, 0.8 U Taq-Polymerase and 5 mM magnesium. Cycling conditions were as follows: initial denaturation at 95°C for 3 min followed by 35 cycles of denaturation at 94°C for 30 s, annealing at 59°C for 30 s, elongation at 72°C for 30 s with a temperature transition rate of 20°C per second and a final elongation step at 72°C for 2 min. Results were compared for at least 200 samples and showed 100% concordance. All genotyping was carried out blinded before data analysis. 10% of samples were genotyped twice for quality control and results showed no discrepancies. PCR-RFLP methods were employed as previously published for MMP1 G-1607GG, MMP3 6A-1171 5A, MMP2 C-735T, MMP9 C-1562T, MMP7 C-153T and MMP7 A-181G 12, 14, 17, 22.
Linkage analysis of the polymorphisms situated on the long arm of chromosome 11 was calculated using the HAPReg software. Calculation of haplotypes and the logistic regression for evaluation of haplotypic effect was conducted with the Thesias software 23.
Statistical analysis
The allele frequencies and their standard deviation for the group of all patients were calculated and the genotype distribution was tested for deviations from Hardy–Weinberg equilibrium (HWE) using the Chi-squared test. The prognostic value of the SNPs was analysed using multivariable logistic regression. Odds ratios were calculated by comparing genotype frequencies in responders (CR and PR) and non-responders (SD and PD) after the second cycle of first-line chemotherapy.
Hazard ratio (HR) estimates with 95% confidence intervals were calculated using Cox proportional hazard models for PFS (adjusted for tumour stage) and for OS (adjusted for tumour stage and sex). Kaplan–Meier survival curves were plotted and the log-rank test was used to test for differences in both PFS and OS. In univariate analysis, tumour stage had a statistically significant influence on the response after the second cycle of NSCLC patients (p = 0.0012), and on PFS and OS of all lung cancer patients (p<0.001) in this study. Therefore, we included stage as adjustment factor in all analyses. Sex had a significant influence on OS (p = 0.044) but not on PFS (p = 0.982) of all lung cancer patients and also not on chemotherapy response. Therefore, it was included as an adjustment factor only for OS. There was no significant influence of age in all three clinical outcomes (for response p = 0.59, for OS p = 0.88, and for PFS p = 0.82). The same was observed for performance status measured by ECOG (for response p = 0.61, for OS p = 0.81, and for PFS p = 0.13). Thus, neither of these two parameters was used as confounding factor.
This analysis was performed for all 349 lung cancer patients. In addition, it was performed separately for the two histological sub-populations of SCLC and NSCLC. The wild-type genotype was considered as reference. We analysed therapy-based subsets including the patients receiving gemcitabine- (143; all NSCLC), etoposide- (161; 158 SCLC, 2 NSCLC and 1 mixed) and platinum-based therapy (256; 117 SCLC and 139 NSCLC). Additionally we analysed early- (IIB+IIIA) and late-stage (IIIB+IV) patients as well as IIIB and IV separately in the histological subgroups.
Allele frequencies were calculated using the website of the International Helmholtz Gemeinschaft 24. All other calculations were performed using the statistical software package SAS version 9.1.3 (SAS Institute, Cary, NC, USA).
RESULTS
The total of 349 lung cancer patients, who received first-line chemotherapy were included into this study. Three different endpoints were analysed: response after the second cycle (R), PFS and OS, each for the whole patient group, as well as for the two subgroups, SCLC and NSCLC.
The distribution of the genotypes of the SNPs analysed is given in table 3⇓. For all SNPs, the distribution of genotypes was within Hardy–Weinberg equilibrium (MMP1 G-1607GG p = 0.65; MMP2 C-735T p = 0.57; MMP3 6A-11715A p = 0.39; MMP7 C-153T p = 0.11; A-181G p = 0.18; MMP9 C-1562T p = 0.71; MMP12 A-82G p = 0.97; and A1082G p = 0.99). The calculated frequencies of the wild-type alleles were MMP2 C-735T 0.88±0.015, MMP3 6A-11715A 0.49±0.024, MMP7 C-153T 0.93±0.012, A-181G 0.46±0.022, MMP9 C-1562T 0.86±0.016, MMP12 A-82G 0.85±0.017, and A1082G 0.93±0.012. These frequencies are comparable to those previously published for Caucasians by other groups 5, 12–14.
Association of matrix metalloproteinase(MMP) single nucleotide polymorphisms and the response after the second cycle of chemotherapy for all lung cancer patients
While for SCLC, chemotherapy is the treatment of choice for all stages, early stage NSCLC patients are preferentially treated surgically. This is reflected in table 1⇑ where response rates of SCLC patients are also much higher compared to NSCLC patients.
Odds ratios and their 95% confidence intervals describing the risk of being a nonresponder for variant allele carriers after multivariate analysis are shown in table 3⇑.
In the group of all lung cancer patients, none of the examined polymorphisms significantly modified the chemotherapy response after the second cycle; only a borderline significant effect was observed for the MMP7 C-153T polymorphism (OR 1.79, 95% CI 0.99–3.20; p = 0.05). Among stage IIIB patients, the risk of being a nonresponder was significantly increased in the aggregated groups of variant allele homozygotes and heterozygotes for MMP2 C-735T (OR 2.79, 95% CI 1.06–7.41; p = 0.02) and MMP3 6A-1171 5A (OR 10.37, 95% CI 1.32–81.46; p = 0.01) (table 3⇑).
An enhanced risk for being a nonresponder was also found for IIB-IIIA patients among the SCLC cohort carrying the MMP2 C-735T polymorphism (OR 12.75, 95% CI 1.39–167.39; p = 0.03) as well as for stage IIIB NSCLC patients carrying the variant allele of the MMP3 6A-1171 5A polymorphism (OR 7.78, 95% CI 1.35–147.84; p = 0.03) (tables 4⇓ and 5⇓).
Association of matrix metalloproteinase(MMP) single nucleotide polymorphisms and the response after the second cycle of chemotherapy for small cell lung cancer patients
Association of matrix metalloproteinase(MMP) single nucleotide polymorphisms and the response after the second cycle of chemotherapy for nonsmall cell lung cancer patients
PFS was significantly affected in the group of NSCLC patients, who received gemcitabine as first-line chemotherapeutic drug and are carriers of the variant allele of the MMP12 A1082G polymorphism (HR 1.88, 95% CI 1.02–3.46; p = 0.04). Individuals with the variant genotype of MMP9 C-1562T in the whole cohort among late stage patients (HR 1.36, 95% CI 1.00–1.84; p = 0.05, data not shown) and also among late-stage SCLC patients (adjusted HR 1.72, 95% CI 1.02–2.89; p = 0.04) (table 6⇓) had a significant shorter PFS than those with the wild-type genotype. The Kaplan–Meier curve describing the progression-free survival for all patients also shows a significant modulation by the MMP9 polymorphism (fig. 1⇓). Individuals in the whole cohort carrying the GG allele of the MMP1 G-1607GG allele had a significantly longer PFS compared with the reference genotype carriers (HR 0.72, 95% CI 0.55–0.95; p = 0.04) and a longer PFS and OS in early stage SCLC (tables 6⇓ and 7⇓).
Selected Kaplan–Meier curves showing different patient groups in relation to their genotypes for matrix metalloproteinase (MMP) 9 and the promoter polymorphisms of MMP12. a) MMP9 C-1562T, progression free survival of all patients (n = 349). ——: CC; - - - -: TT+CT. p = 0.0460. b) MMP12 A-82G, overall survival in small cell lung cancer patients (n = 161). ——: AA; - - - -: GG+AG. p = 0.0321.
Association of matrix metalloproteinase(MMP) single nucleotide polymorphisms and progression free survival (PFS) or overall survival (OS) in early and late-stage small cell lung cancer (SCLC) and nonsmall cell lung cancer (NSCLC)
Association of matrix metalloproteinase(MMP) single nucleotide polymorphisms and progression free survival (PFS) or overall survival (OS)
DISCUSSION
The choice of therapy for the treatment of lung cancer therapy is currently dependent on stage, histology and performance status of the patient. Given that there are a number of chemotherapy treatments to choose from, host genetic polymorphisms which affect therapy success should be taken into account. Genotyping prior to therapy would be feasible because it is a quick, reliable and minimally invasive procedure as only a blood sample is required. Individually tailored treatment could improve patient survival and quality of life and is thus very desirable for the future.
Therefore, we genotyped seven polymorphisms in five different matrix metalloproteinase genes in a hypothesis-driven approach and analysed their impact on three endpoints: chemotherapy response after the second cycle, PFS and OS. Response rates in SCLC and NSCLC are known to differ strongly and are also dependent on therapy received. Therefore, we stratified data for histology, therapy protocol and stage. Stage has been shown to correlate strongly with MMP expression.
A particular strength of this study lies in the analysis of three endpoints, chemotherapy response, PFS and OS. Numbers are adequate for the separate analysis of SCLC and NSCLC. A big limitation for this study is the lack of adequate sample numbers in the stage subgroups. This is illustrated by the fact that certain ORs and HRs are not calculable for the MMP3 polymorphism in the SCLC cohort due to a lack of wild-type carriers in the non-responder group. Small sample sizes in stage subgroups also lead some to inconsistencies in the results between the whole cohort and the histological subgroups; e.g. the nonsignificant point estimate for chemotherapy response for carriers of the MMP2 variant allele points towards a protective effect among NSCLC stage IIIA patients, but is associated with a significant risk of non-response among all Stage IIB and IIIA and Stage IIIB patients. However a number of effects appear consistent, if not significant, across subgroups.
MMP1 G-1607GG variant allele carriers have a longer PFS than individuals homozygous for the G allele. This result is not in concordance with the hypothesis that the GG allele, representing higher MMP1 expression, has harmful effects during therapy and in the later endpoints. We observed a significant protective effect in the whole cohort and especially in early stage SCLC. This result has to be validated in further studies of preferably larger cohorts.
The variant allele of MMP2 correlates significantly with a poorer response in IIIB and SCLC early-stage patients and the variant allele of MMP3 correlates significantly, with a poorer response in IIIB patients, especially NSCLC, but not stage IV patients after the second cycle of chemotherapy.
Regarding MMP2, the result does not reflect the hypothesised benefit of the T allele, which is expected to result in low expression. This is not in accordance with Liu et al. 10, in which rats treated with platinum-based therapy in combination with prinomastat, a specific MMP2 inhibitor, had a higher response. Why this effect is seen in the logistic regression analysis is not clear and should be confirmed in other studies.
The significant results for MMP3 and MMP7 show large confidence intervals, which means the point estimate of the effect size is much less precise. However, they are in accordance with a publication by Blons et al. 6 which suggests modulation of chemotherapy response by altered MMP3 expression. The reason for this could be enhanced Fas ligand cleavage from the cell surface by MMP7 and MMP3 25, leading to inhibition of the extrinsic apoptotic pathway.
A possible explanation for why this effect, in the case of the MMP7C-153T polymorphism is observed only in early-stage patients is that the genotype-dependent modulation of treatment outcome is expected to be stronger before extensive metastasis. Both MMP3 and MMP7 are well known for their involvement in metastasis. In stage IV patients where metastasis has already occurred, other resistance mechanisms which are not specific for MMP3 and MMP7 could be more important for response after the second cycle of chemotherapy. The results presented here, which will require confirmation in further pharmacogenetic studies, suggest a genotype specific modulation of therapy protocol in earlier stages of lung cancer.
The Kaplan–Meier estimate for all lung cancer patients shows a significant harmful effect of the T allele, situated in the promoter region of MMP9 and representing high expression 26, on the progression free survival. This is in accordance with previous publications, which have shown high MMP9 expression as relevant to metastasis and patient outcome. However, this result could not be confirmed by corresponding hazard ratios for the MMP9 C-1562T analysis.
The MMP7 A-181G polymorphism did not show a synergetic effect with C-153T in the logistic regression. SCLC patients carrying the variant allele even have a longer lifetime after starting chemotherapy. This is not in accordance with the assumption, that both polymorphisms have similar effects on the transcription 14. However, haplotypic analysis shows a harmful effect of the GT haplotype, which is due to a combination of the high expression alleles. The reason for this is not clear and further studies are required to validate this result.
Due to the antiangiogenic features of MMP12 27 it was predicted that the hypothetically low expressing G allele has a harmful effect on overall survival and progression free survival. This could not be confirmed by our data. In contrast, it showed a significant protective tendency in the analysis for the overall survival for SCLC patients.
Low expression level also represents a higher plasmin level and a consequently higher activation rate of MMP9 in a MMP3 dependent manner, which might promote tumour growth and shorten PFS and survival 28. Several studies indicate a benefit of MMP12 expression in various cancers 8, 29. However, it was also shown that NSCLC patients who have high MMP12 expression have earlier relapse 30. A reason for this effect could be that there is a link between high MMP12 expression, especially in the macrophages in the lung 31, COPD, inflammation and cancer progression 32–34 and that a lower expression indicates a better outcome for the patients. This would also indicate a tissue dependency especially in men, where cases of COPD appear more frequently.
The function of the exon polymorphism of MMP12 is still unknown, but in this study it associated with bad outcome for NSCLC patients who received gemcitabine as chemotherapy. An explanation for this effect remains speculative but a change in the amino acid sequence of the hemopexin-like domain could alter its affinity or specificity.
Of course, the overall survival and the PFS are parameters which are most clinically relevant. However, in the long run, studies which measure the effect of genotypes on chemotherapy response could, if the results are confirmed by larger independent study cohorts, contribute to better treatment protocols. To our knowledge, this is the first pharmacogenetic study on MMPs and lung cancer that complements analysis of tumour progression and patient survival with an evaluation of the early tumour response. The results prioritise MMP1, MMP2, MMP7, MMP9 and MMP12 as candidate genes to be further investigated as possible predictors of the clinical outcome of chemotherapy in lung cancer patients.
Support statement
This work was supported, in part, by the Cooperation Program in Cancer Research of the Deutsches Krebsforschungszentrum (DKFZ), by the “Deutsche Krebshilfe”, and by the “Verein zur Förderung der Krebsforschung in Deutschland e. V.”
Statement of interest
None declared.
Acknowledgments
The authors are grateful to all patients and staff at the Thoraxklinik GmbH (Heidelberg, Germany) who helped to make this study possible. They also thank those who were involved with sample and/or data collection and archiving, O. Landt (TIB MOLBIOL, Berlin, Germany) for helping with the LightCycler probe design, R. Rausch (Dept of Biostatistics, German Cancer Research Center, Heidelberg) for assisting with statistical analysis as well as M. Timofeeva and L. Knoefel (both Dept of Epigenomics and Cancer Risk Factors, German Cancer Research Center, Heidelberg) for valuable discussions.
- Received August 14, 2008.
- Accepted July 13, 2009.
- © ERS Journals Ltd