Abstract
Matrix metalloproteinase-7 (MMP-7) is a small secreted proteolytic enzyme with broad substrate specificity. Its expression is associated with tumor invasion, metastasis, and survival for a number of cancers. However, data from published studies with individually low statistical power are conflicting. Here, we performed a meta-analysis of 14 publications (16 case–control studies) to better assess the purported relationship. Eligible studies were identified by searching the Pubmed database. Odds ratios (ORs) with 95% confidence intervals (CIs) were estimated to assess the association. Overall, we found that the −181 G allele increased cancer risk in East Asians (G-allele vs. A-allele, OR=1.35, 95% CI: 1.25–1.46, P heterogeneity=0.01; GA vs. AA, OR=1.40, 95% CI: 1.16–1.69, P heterogeneity=0.04; GG+GA vs. AA, OR=1.52, 95% CI: 1.30–1.78, P heterogeneity=0.00). Similarly, in the stratified analysis by cancer type and source of control, significantly increased cancer risk was indicated. Our study showed evidence that MMP7 −181A/G polymorphism may increase cancer risk in the East Asian population. Future studies with larger sample size are warranted to further evaluate this association in greater detail.
Introduction
MMP7 (also known as matrilysin or punctuated metalloproteinase-1, PUMP-1), which is encoded on chromosome 11q21-22, degrades elastin, proteoglycans, fibronectin, and type IV collagen (Wilson and Matrisian, 1999). In a transient transfection study, the promoter activity of the MMP7 −181 (rs11568818) G allele is 2- to 3-fold higher than that of the A allele, which may induce elevation of MMP7 mRNA transcription and subsequently increase its protein levels (Jormsjö et al., 2001). The substitution of A by the G allele implies higher activity of the promoter region, thereby increasing the levels of MMP7, and this is a possible mechanism by which MMP7 −181A/G polymorphism is related to an increase in the development of cancer.
Numerous studies on the association of this polymorphism with cancer susceptibility have been carried out; however, the results remain inconclusive. A quantitative synthesis to accumulate data from different studies is needed to provide better evidence on the relationship. In this article, we performed a meta-analysis of 14 publications (Zhang et al., 2005; Kubben et al., 2006; Li et al., 2006; Lu et al., 2007; Vairaktaris et al., 2007; Woo et al., 2007; Qiu et al., 2008; Singh et al., 2008; Sugimoto et al., 2008; Fang et al., 2010; Malik et al., 2010; Srivastava et al., 2010; Yi et al., 2010; Wu et al., 2011) to estimate the effect of MMP7 −181A/G polymorphism on cancer risk.
Materials and Methods
Literature search and data extraction
We conducted searches on the PubMed database (
Inclusion and exclusion criteria
Eligible studies had to meet the following criteria: (a) the study assessed the correlation between global cancer and MMP7 −181A/G polymorphism; (b) case–control studies; (c) control subjects matched with case patients for age and gender; (d) only full-text manuscripts were included; and (e) all the individuals were East Asians. Major exclusion criteria were (a) no control population; (b) no available genotype frequency; and (c) duplication of the previous publications.
Data extraction
Data were collected on the first author's last name, year of publication, country of origin, cancer type, number of cases and controls, study design, subjects of genotype in cases and controls, source of control, and Hardy–Weinberg equilibrium (HWE) of controls and G% in the control groups.
Statistic analysis
Crude odds ratios (ORs) with 95% confidence intervals (CIs) were used to measure the strength of the association between MMP7 −181A/G polymorphism and cancer risk based on the genotype frequencies in cases and controls. Subgroup analysis stratified by race was performed first. If one cancer type contained only one individual study, then it was combined into the “other cancers” subgroup. Source of control subgroup analysis was performed on two classifications: population-based and hospital-based.
The fixed-effects model and the random-effects model were used to calculate the pooled OR. The statistical significance of the summary OR was determined with the Z-test. Heterogeneity assumption was evaluated with a chi-square-based Q test among the studies. A p-value of>0.05 for the Q-test indicated a lack of heterogeneity among the studies. If the p-value for heterogeneity was<0.05 using Q-test, then the random-effects model (DerSimonian and Laird method) was adopted (DerSimonian and Laird, 1986); in contrast, the fixed-effects model (Mantel–Haenszel method) was chosen (Mantel and Haenszel, 1959). The OR and 95% CI value were calculated under different models. To better evaluate the extent of heterogeneity between studies, the I 2 test was also used. As a guide, I 2 values of<25% may be considered low, values of 25%–50% may be considered moderate, and values of>75% may be considered high (Higgins et al., 2003). For MMP7 −181A/G polymorphism, we investigated the association between genetic variants and cancer risk in the dominant genetic model (GG+GA vs. AA), heterozygote comparison (GA vs. AA), and allelic contrast (G-allele vs. A-allele). The funnel plot asymmetry was assessed with Egger's test. Publication bias was assessed with Egger's test, and p<0.05 was considered statistically significant (Egger et al., 1997). The departure of frequencies of MMP7 −181A/G polymorphism from expectation under HWE was assessed by χ 2 test in controls by using the Pearson chi-square test, and p<0.05 was considered significant. All statistical tests for our meta-analysis were performed with Stata software (version 10.0; StataCorp LP, College Station, TX).
Results
Eligible studies
Overall, we identified 14 articles (16 case–control studies; 14 different first authors) to evaluate the association of MMP7 −181A/G polymorphism with cancer risk. One study (Lu et al., 2007) was excluded, because the same authors published several reports on the same patients and controls, and only the most complete study was selected. The analysis includes 16 studies with 3099 cases and 4280 controls concerning the MMP7 −181A/G polymorphism and cancer risk. Characteristics of studies of MMP7 −181A/G polymorphism are summarized in Table 1. Cases in most of the studies were histologically diagnosed, and controls were free from cancer.
HB, hospital-based; PB, population-based; P HWE, Hardy–Weinberg equilibrium of control; G%, the G allele frequency in control group.
For the MMP7 −181A/G polymorphism in East Asians, there were three studies of gastric cancer, two of cervical cancer, two of colorectal cancer, and seven of other cancers. Hospital-based cases were used in 10 studies. Zhang et al. (2005) investigated three types of cancer in one publication. All the included studies used the polymerase chain reaction restriction fragment length polymorphism method to analysis genotypes. Except for three case–control studies (Kubben et al., 2006; Vairaktaris et al., 2007; Wu et al., 2011), all the data were consistent with the HWE.
Quantitative synthesis and test of heterogeneity
In the overall analysis, significantly increased association could be observed between cancer risk and the variant genotypes of MMP7 −181A/G polymorphism for East Asians, but not for Caucasians: in the allelic contrast (OR=1.35, 95% CI=1.25–1.46, P heterogeneity=0.01, p<0.01, I 2=48.7), the heterozygote comparison (OR=1.40, 95% CI=1.16–1.69, P heterogeneity=0.04, p<0.01, I 2=43.4), and the dominant models (OR=1.52, 95% CI=1.30–1.78, P heterogeneity<0.01, p<0.01, I 2=65.7) (Table 2).
P h, the value for heterogeneity in each subgroup; P, the statistical value; OR, odds ratio.
In the stratified analysis by cancer type in East Asians, a significant association was found between MMP7 –181A/G polymorphism and gastric, cervical, and other cancers (Gastric cancer: OR=1.36, 95% CI=1.16–1.59, P heterogeneity=0.32, p<0.01, I 2=12.0 for G-allele vs. A-allele; OR=1.64, 95% CI=1.19–2.27, P heterogeneity=0.02, p<0.01, I 2=74.5 for GG+GA vs. AA. Cervical cancer: OR=1.17, 95% CI=1.02–1.34, P heterogeneity=0.89, p=0.02, I 2=0.0 for G allele vs. A allele. Other cancer: OR=1.76, 95% CI=1.46–2.11, P heterogeneity=0.01, p<0.01, I 2=62.9 for G allele vs. A allele; OR=1.60, 95% CI=1.26–2.04, P heterogeneity=0.01, p<0.01, I 2=67.8 for GA vs. AA; OR=1.73, 95% CI=1.41–2.12, P heterogeneity<0.01, p<0.01, I 2=74.6 for GG+GA vs. AA). Similarly, in the subgroup of source of control, significantly increased association was found in hospital-based or population-based groups in several genetic models (Table 2).
Sensitivity analysis and publication bias
Sensitivity analysis was performed to assess the influence of each individual study on the pooled OR by sequential removal of individual studies. The results suggested that no individual study significantly affected the overall OR in a dominant manner. The Begg's funnel plot and Egger's test were performed to access the publication bias of literatures. The shape of the funnel plots seemed asymmetrical in allele comparison for MMP7 –181A/G polymorphism, thus suggesting the presence of publication bias. Then, Egger's test was used to provide statistical evidence of funnel plot symmetry. As expected, the results have not shown an obvious evidence of publication bias in overall and East Asian analysis (Table 3).
Discussion
This analysis, based on 16 case–control studies with 3099 cases and 4280 controls, explored the association between MMP7 −181A/G polymorphism and cancer risk. The results showed that there was a significant association between MMP7 single nucleotide polymorphism (SNP) −181 and cancer risk in East Asians. Moreover, the variant genotypes were associated with significantly increased gastric, cervical, and other cancers. A similar trend was detected among both population-based and hospital-based studies. SNPs are the most abundant form of genetic variation in the human genome. By convention, a point mutation is referred to as a SNP when the frequency of the minor (rarer) allele exceeds 1% in at least one population (Matsumura et al., 2005). It is thought that SNPs contribute to inter-individual variability in susceptibility to common diseases such as cancer (Sugimoto et al., 2006). MMP7 is a protease with a broad substrate specificity. It is not only able to degrade elastin, proteoglycans, fibronectin, and type IV collagen, but also can cleave nonmatrix substrates from the cell surface, including VE-cadherin, pro-tumor necrosis factor α, insulin-like growth factor binding proteins, and Fas ligand. Overexpression of MMP7 has been shown to occur in wide types of cancers, such as tumors of esophagus, colorectum, and kidney. Further, MMP7 expression has been shown to be associated with metastasis, disease progression, and decreased survival among patients with esophageal, gastric, colorectal, and renal cancer (Tanioka et al., 2003; Kitoh et al., 2004; Miyata et al., 2006; Wang et al., 2006).
The main conclusion of our meta-analysis is that significantly increased associations were found between whole cancer and MMP7 −181A/G polymorphism in East Asians. Genes may play various roles in different cancers. For this reason, we divided the cancers into four different cancer groups, including gastric cancer, cervical cancer, colorectal cancer, and other cancers. In our study, MMP7 −181A/G allele represents a risk factor in gastric cancer and cervical cancer. Malik et al. (2011) reported that people with the MMP7 −181 GG genotype had a 2.13-fold increased risk of gastric cancer. Singh et al. (2008) reported that individuals with the −181 GG genotype or −181 G allele had a 1.94-fold or 1.37-fold higher risk of cervical cancer, respectively. Our data analysis supports these findings. Our study has several limitations. First of all, the number of published studies included in our meta-analysis was not sufficient for a comprehensive analysis. Since only studies that were indexed by the selected database were included for data analysis, some relevant published studies in other databases may be missed, which may have biased our results. Third, publication bias might have occurred, and our Egger's test results may have a substantial risk of being affected by such bias, although the funnel plots as well as Egger's linear regression tests indicated no remarkable publication biases in this meta-analysis. Fourth, lack of individual data from each study prevents more detailed analysis, such as joint effects of SNP-SNP, gene-environment factors, age, sex, and metastasis/differentiation status. Fifth, one group of authors studied three kinds of cancer using the same healthy people as the control group; however, we calculated the control number three times (Zhang et al., 2005).
In summary, this meta-analysis suggests that MMP7 −181A/G polymorphisms may increase risk for overall cancer or some types of cancers in East Asian populations, but not for Caucasians. We encourage future evaluations on the contribution of this polymorphism to cancer risk in large sample sizes, more comprehensive and well-designed association studies.
Footnotes
Disclosure Statement
No competing financial interests exist. The authors are fully responsible for all the content and editorial decisions and have not received any financial support or other form of remuneration related to the development of this article.
