ABSTRACT

Introduction:

The matrix metalloproteinase-7 (MMP-7) gene -181A>G polymorphism has been reported to be associated with colorectal cancer (CRC) and gastric cancer (GC) susceptibility, yet the results of these previous results have been inconsistent or controversial.

Aim:

To elaborate a meta-analysis to assess the association of -181A>G polymorphism of MMP-7 with CRC and GC risk.

Methods:

Published literature evaluating the association from PubMed, Web of Science, Google Scholar and other databases were retrieved up to April 25, 2018. Pooled odds ratio (OR) and 95% confidence interval (CI) were calculated using random- or fixed-effects model.

Results:

A total of 19 case-control studies, which included eleven studies on CRC (2,169 CRC cases and 2,346 controls) and eight studies on GC (1,545 GC cases and 2,366 controls) were identified. There was a significant association between MMP-7 -181A>G polymorphism and GC risk under the homozygote model (GG vs. AA: OR=1.672, 95% CI 1.161-2.409, p=0.006) and the recessive model (GG vs. GA+AA: OR=1.672, 95% CI 1.319-2.554, p=0.001), but not with CRC. By subgroup analysis based on ethnicity, an increased risk of CRC and GC was found only among Asians.

Conclusions:

This meta-analysis suggests that MMP-7 -181A>G polymorphisms is associated with GC risk, but not with CRC. However, our results clearly showed that the MMP-7 -181A>G polymorphism significantly increased the risk of CRC only in Asians.

HEADINGS:
Matrix metalloproteinase-7; Colorectal neoplasms; Stomach neoplasms; Polymorphism; Single nucleotide; Meta-analysis

RESUMO

Introdução:

O polimorfismo da matriz metaloproteinase-7 (MMP-7) -181A>G tem sido relatado como associado à suscetibilidade dos cânceres colorretal (CRC) e gástrico (GC), mas os resultados desses estudos anteriores foram inconsistentes ou controversos.

Objetivo:

Elaborar metanálise para avaliar a associação do polimorfismo -181A> G da MMP-7 com o risco de CRC e GC.

Métodos:

Revisão da literatura publicada avaliando essa associação no PubMed, Web of Science, Google Acadêmico e outras bases de dados até 25 de abril de 2018. Odds ratio (OR) e o intervalo de confiança de 95% (IC) foram calculados usando dados aleatórios ou modelo de efeitos fixos.

Resultados:

Um total de 19 estudos caso-controle, que incluíram 11 trabalhos sobre CRC (2.169 casos de CCR e 2.346 controles) e oito sobre GC (1.545 casos de GC e 2.366 controles) foram identificados. Houve associação significativa entre o polimorfismo MMP-7 -181A>G e o risco de GC sob o modelo homozigoto (GG vs. AA: OR=1,672, IC 95% 1,161-2,409, p=0,006) e o modelo recessivo (GG vs. GA + AA: OR=1,672, IC 95% 1,319-2,554, p=0,001), mas não com CRC. Por análise de subgrupos com base na etnia, um risco aumentado de CRC e GC foi encontrado apenas entre os asiáticos.

Conclusões:

Esta metanálise sugere que os polimorfismos MMP-7 -181A>G estão associados ao risco de GC, mas não ao CRC. No entanto, estes resultados mostraram claramente que o polimorfismo MMP-7 -181A>G aumentou significativamente o risco de CRC apenas em asiáticos.

DESCRITORES:
Metaloproteinase 7 da matriz; Neoplasias colorretais; Neoplasias gástricas; Polimorfismo de nucleotídeo único; Metanálise

INTRODUCTION

Nowadays, gastrointestinal-related cancers especially gastric cancer (GC) and colorectal cancer (CRC) are major causes of cancer-related death worldwide ⁷7 Forat-Yazdi M, Gholi-Nataj M, Neamatzadeh H, Nourbakhsh P, Shaker-Ardakani H. Association of XRCC1 Arg399Gln polymorphism with colorectal cancer risk: A huge meta analysis of 35 studies. Asian Pacific J Cancer Prev [Internet]. 2015;16(8).^,2222 Namazi A, Forat-Yazdi M, Jafari MA, Foroughi E, Farahnak S, Nasiri R, et al. Association between polymorphisms of ERCC5 gene and susceptibility to gastric cancer: A systematic review and meta-analysis. Asian Pacific J Cancer Prev. 2017;18(10).^,2323 Namazi A, Forat-Yazdi M, Jafari M, Farahnak S, Nasiri R, Foroughi E, et al. Association of interleukin-10 -1082 A/G (rs1800896) polymorphism with susceptibility to gastric cancer: meta-analysis of 6,101 cases and 8,557 controls. Arq Gastroenterol. 2018 Mar;55(1):33-40.. Globally, CRC and GC are the third and fourth leading cause most common cancers, making up 10% and 7% of cases all new cancer cases, respectively ²³23 Namazi A, Forat-Yazdi M, Jafari M, Farahnak S, Nasiri R, Foroughi E, et al. Association of interleukin-10 -1082 A/G (rs1800896) polymorphism with susceptibility to gastric cancer: meta-analysis of 6,101 cases and 8,557 controls. Arq Gastroenterol. 2018 Mar;55(1):33-40.^,2727 Tak DH, Moon HS, Kang SH, Sung JK, Jeong HY. Prevalence and Risk Factors of Gastric Adenoma and Gastric Cancer in Colorectal Cancer Patients. Gastroenterol Res Pract. 2016 Dec 26;2016:1-7.. It is well-known that the development of CRC and GC can be induced by the interactions of multiple genetic and environmental factors in complex ways. However, the pathogenesis of CRC differs from that of GC in that it is affected by ethnic background, regional, life style, environmental factors and molecular pathogenesis²⁷27 Tak DH, Moon HS, Kang SH, Sung JK, Jeong HY. Prevalence and Risk Factors of Gastric Adenoma and Gastric Cancer in Colorectal Cancer Patients. Gastroenterol Res Pract. 2016 Dec 26;2016:1-7.. Currently, several genes have been reported to be associated with CRC and GC, and the matrix metalloproteinases (MMPs) genes has received increasing attention^{10, 29,31}.

MMPs are classified as a large family of zinc-containing proteases, which involved in normal physiological and pathological processes such as degradation and remolding of extracellular matrix, embryonic development, reproduction and cancer⁶6 Fanjul-Fernández M, Folgueras AR, Cabrera S, López-Otín C. Matrix metalloproteinases: Evolution, gene regulation and functional analysis in mouse models. Biochim Biophys Acta - Mol Cell Res. 2010 Jan 1;1803(1):3-19.^,1313 Klein T, Bischoff R. Physiology and pathophysiology of matrix metalloproteases. Amino Acids. 2011 Jul;41(2):271-90.. MMP-7, the smallest member of MMP family, is an endopeptidases with broad substrate specificity, which break down extracellular matrix (ECM) by degrading macromolecules including casein, type I, II, IV, and V gelatins, fibronectin, and proteoglycan¹⁸18 Mittal R, Patel AP, Debs LH, Nguyen D, Patel K, Grati M, et al. Intricate Functions of Matrix Metalloproteinases in Physiological and Pathological Conditions. J Cell Physiol. 2016 Dec 1;231(12):2599-621.. In addition, MMP-7 is one of the main regulatory enzymes involved in apoptosis by releases the Fas ligand (FasL) from the membrane then induces apoptosis of neighboring cells, or decreases cancer-cell apoptosis²¹21 Musial K, Zwolinska D. Matrix metalloproteinases and soluble Fas/FasL system as novel regulators of apoptosis in children and young adults on chronic dialysis. Apoptosis. 2011 Jul;16(7):653-9.. Thus, MMP-7 promotes cell survival by resisting apoptosis through cleaving FasL. MMP7 is potentially involved in tumor metastasis and inflammatory processes¹⁰10 Ke P, Wu Z-D, Wen H-S, Ying M-X, Long H-C, Qing L-G. Current evidence on associations between the MMP-7 (-181A>G) polymorphism and digestive system cancer risk. Asian Pac J Cancer Prev. 2013;14(4):2269-72.^,3131 Yang X, Liu Y, Yang Y, Li B. Update meta-analysis on MMP-7 -181A>G polymorphism and cancer risk: Evidence from 25 studies. Gene. 2013 Jun 1;521(2):252-8..

MMP-7 gene (also known as matrilysin) is localized on chromosome 11q21-q22 and contains 13 exons³⁵35 Zhang Y, Chen Q. Relationship between matrix metalloproteinases and the occurrence and development of ovarian cancer. Brazilian J Med Biol Res = Rev Bras Pesqui medicas e Biol. 2017 May 18;50(6):e6104.. The single nucleotide polymorphism (SNP) -181A>G in promoter region of MMP-7 gene has been considered to be a candidate SNP for various conditions including gastrointestinal related malignancies¹1 Adachi Y, Yamamoto H, Itoh F, Hinoda Y, Okada Y, Imai K. Contribution of matrilysin (MMP-7) to the metastatic pathway of human colorectal cancers. Gut. 1999 Aug 1;45(2):252-8.^,1010 Ke P, Wu Z-D, Wen H-S, Ying M-X, Long H-C, Qing L-G. Current evidence on associations between the MMP-7 (-181A>G) polymorphism and digestive system cancer risk. Asian Pac J Cancer Prev. 2013;14(4):2269-72.^,3131 Yang X, Liu Y, Yang Y, Li B. Update meta-analysis on MMP-7 -181A>G polymorphism and cancer risk: Evidence from 25 studies. Gene. 2013 Jun 1;521(2):252-8.. There is clear evidence that MMP-7 gene up-regulation is significantly related to the promoter activity variation of the -181A>G16. Molecular epidemiological studies have reported the association of MMP-7 -181A>G polymorphism with CRC and GC risk, but the results remain conflicting rather than conclusive. Several studies previously have performed on the association of MMP-7 -181A>G polymorphism with CRC and GC risk¹⁰10 Ke P, Wu Z-D, Wen H-S, Ying M-X, Long H-C, Qing L-G. Current evidence on associations between the MMP-7 (-181A>G) polymorphism and digestive system cancer risk. Asian Pac J Cancer Prev. 2013;14(4):2269-72.^,3131 Yang X, Liu Y, Yang Y, Li B. Update meta-analysis on MMP-7 -181A>G polymorphism and cancer risk: Evidence from 25 studies. Gene. 2013 Jun 1;521(2):252-8.. However, these studies had opposite observations and additional case-control studies with larger sample sizes have been published since then. Hence, the association of MMP-7 -181A>G polymorphism remains unknown.

Therefore, we have performed a meta-analysis of all eligible studies to derive more precise estimation of the association of MMP-7 -181A>G polymorphism with CRC and GC risk.

METHODS

Literature search

Following PRISMA guidance, we searched the electronic literature databases including PubMed, EMBASE, Elsevier, Science Direct, Wan Fang, Chinese National Knowledge Infrastructure (CNKI) and Chinese Biomedical Literature for all relevant articles published up to April 25, 2018. The search strategies were based on combinations of the following key words: (‘’matrix metalloproteinase-7’’ OR ‘’MMP-7’’ OR ‘’matrilysin’’ OR uterine metalloproteinase OR ‘’pump-1 protease’’ OR ‘’PUMP-1’’) AND (‘’-181A>G OR ‘’rs11568818’’) AND (‘’colorectal cancer’’ OR ‘’CRC’’) AND (‘’gastric adenocarcinoma’’ OR ‘’stomach cancer’’ OR ‘’gastric cancer’’ OR ‘’GC’’) AND (‘’gene” or ‘’allele” or ‘’genotype” or ‘’mutation” or ‘’variant” or ‘’variation” or ‘’polymorphism”), without any restriction on language. Review articles were hand-searched to find additional eligible studies and only published studies with full-text articles were included. We excluded studies that were not full-length publications articles or letters in peer-reviewed English journals. When the same patient population was included in different articles, the one with the largest population of participants or the most recent one was selected.

Inclusion and exclusion criteria

Inclusion criteria was defined as follows: 1) published studies and contained original data; 2) case-control studies; 3) evaluating the association of MMP-7 -181A>G polymorphism with CRC and GC risk; and 4) sufficient published data available to estimate an odds ratio (OR) with 95% confidence interval (CI). Major reasons for exclusion of studies were as follows: 1) only case population, family based or linkage studies; 2) studies that could not offer the number of cases and controls or other essential data; 3) reviews, abstracts, letters to editor, case reports or animal studies; 4) duplicate of previous publication or studies with overlapping patient populations; and 5) studies without histologically confirmed of CRC and GC. For more than two studies with overlapping data, the study with the most subjects or newest published data was selected.

Data extraction

Data were carefully extracted from all eligible publications by two of the authors independently. If the study provided stratum information, the data coming from similar stratum were added up to make full use of the data. Disagreements between the two authors were resolved by discussing the results with a third one. For each study, the following variables were collected: first author’s name, year of publication, country, ethnicity of participants, number of cases and controls, genotyping methods, and allele numbers and genotype distributions in cases and controls, minor allele frequencies (MAFs) in control subjects, and the results of Hardy-Weinberg equilibrium (HWE) test. Different ethnicities were categorized as Asian, Caucasian and Latinos (mixed). Study designs were stratified to population-based studies and hospital-based studies.

Statistical analysis

The strength of the association of MMP-7 -181A>G polymorphism with CRC and GC risk was measured using odds ratios (ORs) with 95% confidence intervals (CIs). The statistical significance of the pooled OR was assessed with the Z-test and p<0.05 was considered significant. The pooled ORs were performed under five genetic models, i.e., allele (G vs. A), heterozygote (GA vs. AA), homozygote (GG vs. AA), dominant (GG+GA vs. AA) and recessive (GG vs. GA+AA). The between-study heterogeneity was evaluated by a chi-square-based Q test, which p value for the Q-test less than 0.10 indicates existing heterogeneity among studies. In addition, the I² statistics was used to quantify the proportion of the total variation across studies due to heterogeneity. A high value of I2 indicated a higher probability of the existence of heterogeneity (I²=0% to 25%, no heterogeneity; I²=25% to 50%, moderate heterogeneity; I²=50% to 75%, large heterogeneity; and I²=75% to 100%, extreme heterogeneity). A random-effects (DerSimonian- Laird method) or fixed-effects (Mantel-Haenszel method) model was used to calculate pooled effect estimates in the presence or absence of heterogeneity. HWE of genotype distribution in the controls of included studies was conducted using by Pearson’s x ²2 Banday MZ, Sameer AS, Mir AH, Mokhdomi TA, Chowdri NA, Haq E. Matrix metalloproteinase (MMP) -2, -7 and -9 promoter polymorphisms in colorectal cancer in ethnic Kashmiri population - A case-control study and a mini review. Gene. 2016 Sep 1;589(1):81-9. test, in which p-value less than 0.05 was considered significantly deviating from HWE. Subgroup analyses were performed by ethnicity, genotyping method, HWE status, source of controls and cancer type (CRC and GC). In addition, to consider the possible sources of heterogeneity, the studies we stratified. To validate the reliability of the results, sensitivity analysis was performed though omitting one case-control study each time, as well as limiting this meta-analysis to studies which were conformed to HWE. Funnel plots and Egger’s linear regression test were used to diagnose potential publication bias (p<0.05). All analyses were performed with the comprehensive meta-analysis (CMA) 2.0 software (Biostat, USA). Two-sided p-values<0.05 were considered statistically significant.

RESULTS

Extraction process and study characteristics

The flow diagram of literature search was given in Figure 1. The initial search of databases yielded 103 relevant publications based on our literature search strategy, and an additional one study was identified through hand searching. However, 36 of them were ruled out because of duplicate results obtained from multiple databases, 68 articles remained. In addition, after the titles and abstracts of the 68 articles were reviewed, 49 full-text irrelevant studies were excluded. Finally, 19 eligible case-control studies with 3,714 cases and 4,712 controls were included in this meta-analysis. The characteristics of studies included in the current meta-analysis are shown in Table 1. Among these studies, eleven studies with 2,169 CRC cases and 2,346 controls were on CRC²2 Banday MZ, Sameer AS, Mir AH, Mokhdomi TA, Chowdri NA, Haq E. Matrix metalloproteinase (MMP) -2, -7 and -9 promoter polymorphisms in colorectal cancer in ethnic Kashmiri population - A case-control study and a mini review. Gene. 2016 Sep 1;589(1):81-9.

3 Dziki L, Przybylowska K, Majsterek I, Trzcinski R, Mik M, Sygut A. A/G Polymorphism of the MMP-7 Gene Promoter Region in Colorectal Cancer. Polish J Surg. 2011 Jan 1;83(11):622-6.^-44 Fang W-L, Liang W-B, He H, Zhu Y, Li S-L, Gao L-B, et al. Association of Matrix Metalloproteinases 1, 7, and 9 Gene Polymorphisms with Genetic Susceptibility to Colorectal Carcinoma in a Han Chinese Population. DNA Cell Biol. 2010 Nov;29(11):657-61.^,88 Ghilardi G, Biondi ML, Erario M, Guagnellini E, Scorza R. Colorectal carcinoma susceptibility and metastases are associated with matrix metalloproteinase-7 promoter polymorphisms. Clin Chem. 2003 Nov;49(11):1940-2.^,1616 Lièvre A, Milet J, Carayol J, Le Corre D, Milan C, Pariente A, et al. Genetic polymorphisms of MMP1, MMP3 and MMP7gene promoter and risk of colorectal adenoma. BMC Cancer. 2006 Dec 24;6(1):270.^,1919 Moreno-Ortiz JM, Gutierrez-Angulo M, Partida-Pérez M, Peregrina-Sandoval J, Ramírez-Ramírez R, Muñiz-Mendoza R, et al. Association of MMP7-181A/G and MMP13-77A/G polymorphisms with colorectal cancer in a Mexican population. Genet Mol Res. 2014 Feb 14;13(2):3537-44.^,2020 Motovali-Bashi M, Sadeghi Z, Hemati S. Relationship between -181 (A/G) region single nucleotide polymorphisms of matrilysin gene promoter and the onset and prevalence of colorectal cancer using tetra-primer ARMS PCR and RFLP-PCR techniques [Internet]. Vol. 7, Basic & Clinical Cancer Research. Tehran University of Medical Sciences, Cancer Institute of Iran, Cancer Research Center; 2015. 20-31.^,2424 Ohtani H, Maeda N MY. Functional polymorphisms in the promoter regions of matrix metalloproteinase-2,-3,-7,-9 and TNF- alpha genes, and the risk of colorectal neoplasm in Japanese. Yonago Acta Med. 2009;52:47-56.^,2828 Woo M, Park K, Nam J, Kim JC. Clinical implications of matrix metalloproteinase-1, -3, -7, -9, -12, and plasminogen activator inhibitor-1 gene polymorphisms in colorectal cancer. J Gastroenterol Hepatol. 2007 Jul;22(7):1064-70.^,3333 Yueh T-C, Wu C-N, Hung Y-W, Chang W-S, Fu C-K, Pei J-S, et al. The Contribution of MMP-7 Genotypes to Colorectal Cancer Susceptibility in Taiwan. Cancer Genomics Proteomics. 2018 Apr 24;15(3):207-12., and eight studies with 1,545 GC cases and 2,366 controls were on GC⁵5 Fang W-L, Liang W-B, Gao L-B, Zhou B, Xiao F-L, Zhang L. Genetic polymorphisms in Matrix Metalloproteinases -1 and -7 and susceptibility to gastric cancer: an association study and meta-analysis. Iran J Allergy Asthma Immunol. 2013 Jul 9;12(3):203-10.^,1111 Kesh K, Subramanian L, Ghosh N, Gupta V, Gupta A, Bhattacharya S, et al. Association of MMP7 -181A?G Promoter Polymorphism with Gastric Cancer Risk. J Biol Chem. 2015 Jun 5;290(23):14391-406.^,1212 Kim JH, Pyun JA, Lee KJ, Cho SW, Kwack KB. [Study on association between single nucleotide polymorphisms of MMP7, MMP8, MMP9 genes and development of gastric cancer and lymph node metastasis]. Korean J Gastroenterol. 2011 Nov 25;58(5):245-51.^,1414 Kubben FJGM, Sier CFM, Meijer MJW, van den Berg M, van der Reijden JJ, Griffioen G, et al. Clinical impact of MMP and TIMP gene polymorphisms in gastric cancer. Br J Cancer. 2006 Sep 29;95(6):744-51.^,1515 Li J, Tian M, Zhao A. T2028 Polymorphism in the Promoter Region of the Metalloproteinase-7 Increases Susceptibility and Risk of Metastasis of Gastric Adenocarcinoma. Gastroenterology. 2008 Apr 1;134(4):A-603.^,1717 Malik MA, Zargar SA, Mittal B. Role of the metalloproteinase-7 (181A&gt;G) polymorphism in gastric cancer susceptibility: a case control study in Kashmir valley. Asian Pac J Cancer Prev. 2011;12(1):73-6.^,2626 Sugimoto M, Furuta T, Kodaira C, Nishino M, Yamade M, Ikuma M, et al. Polymorphisms of matrix metalloproteinase-7 and chymase are associated with susceptibility to and progression of gastric cancer in Japan. J Gastroenterol. 2008 Oct 29;43(10):751-61.^,3434 Zhang J, Jin X, Fang S, Wang R, Li Y, Wang N, et al. The functional polymorphism in the matrix metalloproteinase-7 promoter increases susceptibility to esophageal squamous cell carcinoma, gastric cardiac adenocarcinoma and non-small cell lung carcinoma. Carcinogenesis. 2005 Oct 1;26(10):1748-53.. By ethnics, there were 13 studies of Asians (countries: Korea, Japan, China, Iran, Kashmir, Taiwan, and India), four studies of Caucasians (countries: Italy, France, Poland and Netherland), and two studies of Latinos (countries: Brazil and Mexico). According to the control source, ten studies were hospital-based, eight studies were population-based and one study was not clear. The studies used four different genotyping methods including direct sequencing, TaqMan, polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis and tetra-primer amplification refractory mutation system-polymerase chain (ARMS-PCR). All of the studies indicated that the distribution of genotypes in the controls was consistent with HWE except for two studies (Table 1).

Quantitative synthesis

Overall study

Table 2 listed the main results of the meta-analysis of MMP-7 -181A>G polymorphism with CRC and GC risk. We pooled all the 19 case-control studies together to assess the overall association of MMP-7 -181A>G polymorphism with CRC and GC risk. Overall, no significant main effects on CRC and GC susceptibility were observed in the overall population under all the five genetic models, i.e., allele (G vs. A: OR=1.049, 95% CI 0.889-1.239, p=0.570, Figure 2A), heterozygote (GA vs. AA: OR=1.083, 95% CI 0.813-1.443, p=0.586), homozygote (GG vs. AA: OR=0.982, 95% CI 0.701-1.375, p=0.915), dominant (GG+GA vs. AA: OR=1.061, 95% CI 0.869-1.296, p=559) and recessive (GG vs. GA+AA: OR=1.084, 95% CI 0.786-1.495, p=0.622).

Colorectal cancer

Table 2 also listed the main results of the meta-analysis of MMP-7 -181A>G polymorphism with CRC risk. When all the 11 eligible studies were pooled into the meta-analysis of MMP-7 -181A>G polymorphism, we have not found evidence of a significant MMP-7 -181A>G polymorphism with CRC risk under all the five genetic models. In the subgroup analysis by ethnicity, significantly increased risk of CRC was observed in Asians under three genetic models, i.e., allele (G vs. A: OR=0.798, 95% CI 0.661-0.964, p=0.019, Figure 2A), homozygote (GG vs. AA: OR=0.490, 95% CI 00.286-0.838, p=0.009) and recessive model (GG vs. GA+AA: OR=0.530, 95% CI 0.340-0.826, p=0.005), but not in Caucasians and Latinos populations (Table 2).

We also performed subgroup analyses based on the source of control and genotyping method, when it was available (Table 3). The hospital-based subgroup analysis revealed that the presence of the MMP-7 -181A>G polymorphism was related to a higher risk of CRC under the homozygote model (GG vs. AA: OR=0.671, 95% CI 0.484-0.951, p=0.023). In the PCR-RFLP group, significantly increased association between MMP-7 -181A>G polymorphism and CRC risk was found under the homozygote model (GG vs. AA: OR=0.680, 95% CI 0.486-0.950, p=0.024, Table 3).

Gastric cancer

Table 2 also listed the main results of the meta-analysis of MMP-7 -181A>G polymorphism with GC risk. There was a significant association between MMP-7 -181A>G polymorphism and GC risk under two genetic models, i.e., homozygote (GG vs. AA: OR=1.672, 95% CI 1.161-2.409, p=0.006, Fig 2B) and recessive (GG vs. GA+AA: OR=1.835, 95% CI 1.319-2.554, p=0.001). Similarly, when stratified by ethnicity, a significant association between MMP-7 -181A>G polymorphism and increased risk of GC was detected among Asians under the homozygote model (GG vs. AA: OR=1.975, 95% CI 1.331-2.934, p=0.006) and the recessive model (GG vs. GA+AA: OR=2.022, 95% CI 1.416-2.886, p=0.001).

The studies were further stratified on the basis of source of controls (Table 3). When stratifying by source of control, a significant association between MMP-7 -181A>G polymorphism and increased risk of GC was detected in population-based studies under the recessive model (GG vs. GA+AA: OR=1.819, 95% CI 1.173-2.819, p=0.008), and in hospital-based studies under two genetic models, i.e., homozygote (GG vs. AA: OR=1.884, 95% CI 1.107-3.204, p=0.020) and recessive (GG vs. GA+AA: OR=1.858, 95% CI 1.124-3.069, p=0.016).

Heterogeneity analysis

Heterogeneity was detected among studies under all the five genetic models, i.e., allele (G vs. A: I2=70.90%, PH=0.001), heterozygote (GA vs. AA: I2=82.36%, PH=0.001), homozygote (GG vs. AA: I2=54.65%, PH=0.006), dominant (GG+GA vs. AA: I2=65.43%, PH=0.001) and recessive (GG vs. GA+AA: I2=65.4%, PH=0.001). Thus, to explore the potential sources of heterogeneity across studies, we assessed the pooled ORs via stratification by cancer type, ethnicity, genotyping method, HWE status and source of controls. The results showed that the heterogeneity effectively removed by subgroup analyses based on ethnicity among studies on CRC. Therefore, we found that genotyping method, HWE status and source of controls did not contribute to substantial heterogeneity among the meta-analysis.

Sensitivity analysis

To evaluate the effect of individual study on the pooled ORs and stability of the meta-analysis results, we excluded one study at a time. However, the omission of any single study made no significant difference, suggesting that the results of this meta-analysis were stable. Moreover, sensitivity analysis was performed after excluding HWE-violating studies, and the corresponding pooled ORs were not qualitatively altered (data not shown).

Publication bias

Both Begg’s funnel plot and Egger’s test were performed to assess the publication bias of literature. Begg’s funnel plots did not reveal any evidence of obvious asymmetry under all five genetic models in the overall meta-analysis. For example, the shape of the funnel plot did not indicate any evidence of obvious asymmetry under the allele model (Figure 3), and the Egger’s test suggested the absence of publication bias (PBeggs=0.293 and PEggers=0.483). However, the results of Egger’s regression test showed evidence of publication bias among Caucasian’s studies on CRC under the recessive model (GG vs. GA+AA: PBegg’s=0.296, PEggers=0.001).

Minor allele frequencies (MAFs)

The minor allele frequencies (MAFs) of the MMP-7 -181A>G polymorphism are shown in Table 2. The allele and genotype distributions of MMP-7 -181A>G polymorphism exhibited ethnic variations. The MMP-7 -181A>G polymorphism MAF in overall populations, Asians, Caucasians and Latinos were 37.4% (4.0%-70.8%), 37.4% (4.0%-70.8%), 43.65% (40.0%-47.30%), and 37.75% (37.1%-38.4%), respectively.

DISCUSSION

MMP-7 gene is implicated in cancer susceptibility and metastasis in a variety of gastrointestinal-related cancers ^{19, 11}. In the recent decade, several molecular epidemiological studies have been performed to assess the association of MMP-7 -181A>G polymorphism with CRC and GC risk. However, the results were conflicting. Thus, we conducted a comprehensive meta-analysis involving published data, to assess the strength of association of MMP-7 -181A>G polymorphism CRC and GC risk. In current meta-analysis, a total of 19 case-control studies including 3,714 cases and 4,712 controls were recruited. The subgroup analysis by cancer type, showed that there was significant association between the MMP-7 -181A>G polymorphism and increased risk of GC in overall estimations, but not with CRC. In a stratified analysis by ethnicity, our results indicated that MMP-7 -181A>G polymorphism was associated with a significantly increased risk of CRC and GC in Asians. Moreover, stratified analysis according to source of controls and genotyping method revealed a significantly increased risk of CRC and GC in participants with the MP-7 -181A>G polymorphism in those studies involving PCR-RFLP, population-based and hospital-based (Table 3).

Our results were consistent with the previous meta-analyses on MMP-7 -181A>G polymorphism and GC risk. Thus, MMP-7 -181A>G polymorphism might serve as a susceptibility marker to GC risk. However, our results are inconsistent with the previous meta-analyses on MMP-7 -181A>G polymorphism and CRC risk. In 2013, two meta-analyses by Ke et al and Yang et al were conducted to examine the association between MMP-7 -181A>G polymorphism and CRC risk ¹⁰10 Ke P, Wu Z-D, Wen H-S, Ying M-X, Long H-C, Qing L-G. Current evidence on associations between the MMP-7 (-181A&gt;G) polymorphism and digestive system cancer risk. Asian Pac J Cancer Prev. 2013;14(4):2269-72.^,3131 Yang X, Liu Y, Yang Y, Li B. Update meta-analysis on MMP-7 -181A&gt;G polymorphism and cancer risk: Evidence from 25 studies. Gene. 2013 Jun 1;521(2):252-8.. Ke et al in meta-analysis of seven case-control have found a significant association between MMP-7 -181A>G polymorphism and CRC under the homozygote model (GG vs. AA, OR=1.13, 95% CI=1.01-1.26) ¹⁰10 Ke P, Wu Z-D, Wen H-S, Ying M-X, Long H-C, Qing L-G. Current evidence on associations between the MMP-7 (-181A&gt;G) polymorphism and digestive system cancer risk. Asian Pac J Cancer Prev. 2013;14(4):2269-72.. Similarly, Yang et al in a meta-analysis of seven studies with 1,502 CRC cases and 1,602 controls found significant association between MMP-7 -181A>G polymorphism and CRC under the homozygote model (OR=1.31, 95% CI 1.02-1.69) ³¹31 Yang X, Liu Y, Yang Y, Li B. Update meta-analysis on MMP-7 -181A&gt;G polymorphism and cancer risk: Evidence from 25 studies. Gene. 2013 Jun 1;521(2):252-8.. However, their findings about MMP-7 -181A>G polymorphism and CRC risk essentially remains an open field, as the number of studies is considerably smaller than that needed to yield a robust conclusion. In addition, the previous meta-analyses did not perform stratified analysis by ethnicity to identify possible association of MMP-7 -181A>G polymorphism with CRC among different ethnic groups.

Between-studies heterogeneity plays an important role when performing a meta-analysis ³²32 Yazdi MM, Jamalaldini MH, Sobhan MR, Jafari M, Mazaheri M, Zare-Shehneh M, et al. Association of ESRa Gene Pvu II T&gt;C, XbaI A&gt;G and BtgI G&gt;A Polymorphisms with Knee Osteoarthritis Susceptibility: A Systematic Review and Meta-Analysis Based on 22 Case-Control Studies. Arch bone Jt Surg. 2017 Nov;5(6):351-62.. Heterogeneity could result from study design, genotyping error, selection bias, population stratification, sample size, allelic heterogeneity, or chance ⁹9 Gohari M, Neámatzadeh H, Jafari MA, Mazaheri M, Zare-Shehneh M, Abbasi-Shavazi E. Association between the p53 codon 72 polymorphism and primary open-angle glaucoma risk: Meta-analysis based on 11 case-control studies. Indian J Ophthalmol. 2016 Oct;64(10):756-61.^,2525 Sobhan MR, Mahdinezhad-Yazdi M, Aghili K, Zare-Shehneh M, Rastegar S, Sadeghizadeh-Yazdi J, et al. Association of TNF-a-308?G?&gt;?A and -238G?&gt;?A polymorphisms with knee osteoarthritis risk: A case-control study and meta-analysis. J Orthop. 2018 Sep 1;15(3):747-53.. Therefore, finding the source of heterogeneity is very important for the final result of meta-analysis. Through performing sub group analysis, we found that the heterogeneity could not be explained by genotyping method, HWE status and source of controls in this meta-analysis. However, the results showed that the heterogeneity effectively removed by subgroup analysis based on ethnicity among studies on CRC, indicating that studies among Asian populations regarding CRC might be a source of the heterogeneity in our meta-analysis.

This meta-analysis had three main strengths. First, this is the biggest and most recent meta-analysis of the association of MMP-7 -181A>G polymorphism with CRC and GC risk, and it was more powerful than previous single case-control studies. Second, this is the first meta-analysis by subgroup analysis showed that the MMP-7 -181A>G polymorphism was associated with CRC risk in Asians. Third, a comprehensive searching strategy from several electronic databases with manual search made the eligible studies included as much as possible.

Despite the clear strengths of this meta-analysis, limitations of our meta-analysis should be noted. First, although all the eligible studies were included to this meta-analysis, the sample size of the included studies was not large enough, which could increase the likelihood of type I and type II errors. Therefore, there was a lack of statistical power to better evaluate the association of MMP-7 -181A>G polymorphism with CRC and GC risk. Second, most of included studies in the present meta-analysis mainly provided data in Asians. In addition, the sample size was relatively small for stratified analyses by ethnicity and might not have provided sufficient power to estimate the association of MMP-7 -181A>G polymorphism among different ethnic groups. Third, although the funnel plot and Egger’s test did not show evidence of publication bias in overall estimations, the influence of bias in the present analysis could not be completely excluded. For example, the negative findings are usually difficult to get published, or in this meta-analysis we have included only studies published in English, which produced selection bias at the start of our study. Fourth, the present meta-analysis was based primarily on unadjusted effect estimates, because most studies did not provide the adjusted OR and 95%CI controlling for potential confounding factors, thus the effect estimates were relatively imprecise. If individual data were available, adjusted ORs could be obtained to give a more precise analysis. Finally, it is well known that CRC and GC are multifactor conditions; however, the effects of gene-gene and gene-environment interactions were not addressed in the current meta-analysis.

CONCLUSION

This meta-analysis indicated that the MMP-7 -181A>G polymorphism might be a risk factor for susceptibility to GC in overall estimations and in Asians. However, our results clearly showed that the MMP-7 -181A>G polymorphism significantly increased the risk of CRC only in Asians. Considering the limited sample size and ethnicities, well-designed studies taking into consideration of gene-gene and gene-environment interactions should be performed to confirm our results.

REFERENCES

Publication Dates

History