Chronic obstructive pulmonary disease (COPD) is characterized by decreased air flow and is associated with abnormal chronic inflammation in the airways and extensive tissue remodeling. Matrix metalloproteinase-7 (MMP7) is produced primarily by the epithelium of many organs, including the lungs. A functional MMP7 –181A>G (rs11568818) promoter polymorphism influences the binding of nuclear regulatory proteins modulating the transcription of the gene. In this study, we genotyped 191 patients with COPD for MMP7 –181A>G single nucleotide polymorphism (SNP) and 215 control subjects using the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method and explored the role of that polymorphism as a risk factor for COPD. There were no differences in the genotype and allele distribution of the MMP7 –181A>G SNP between the COPD patients and control groups (p = 0.341 and p = 0.214). However, the carries of the G allele (AG and GG genotypes), appeared to develop COPD significantly earlier than those with the AA genotype (61.01 ± 10.11 vs. 64.87 ± 9.00 years, p = 0.032). When the genotype distribution was studied only in the groups of patients (n = 76) and controls (n = 106) younger than 60 years, we found significantly higher frequency of the carriers of the G allele in COPD patients than in the controls, determining about a 3-fold higher risk for COPD [odds ratio (OR) –3.33, 1.36-8.14, p = 0.008 for GG, and OR = 2.91, 1.38-6.13, p = 0.005 for AG+GG]. Based on our results, the MMP7 –181A>G promoter variant may influence early development of COPD. This effect could be attributed to the increased production of the enzyme resulting in enhanced airway wall protein degradation and injury.
Chronic obstructive pulmonary disease (COPD) is a preventable and treatable disease with some significant extrapulmonary effects that may contribute to its severity in individual patients. Its pulmonary component is char-acterized by airflow limitation that is not fully reversible. The airflow limitation is associated with an abnormal inflammatory response of the lung to noxious particles or gases .
Smoking is one of the main risk factors for COPD, but since not all smokers develop COPD, as well as the fact that the disease often develops in middle age, it is suggested that other factors may play a role in the pathogenesis such as genetic factors . Inhalation of cigarette smoke, organic and/or inorganic dust, chemical agents and particle matters increase the risk of developing COPD. The presence of these irritants, may lead to chronic inflammation and structural changes in the lung due to repeated injury and repair . Pathological changes characteristic for COPD are found in the proximal airways, peripheral airways, lung parenchyma and pulmonary vasculature .
One of the main roles of the epithelial cells is to provide a barrier against pathogens and to release antimicrobial products. By producing chemoattractants and adhesion molecules, epithelial cells contribute to the migration of the inflammatory cells to injury sites [5,6]. Epithelial damage is an important characteristic of several pulmonary diseases including COPD. Numerous enzymes, proteins and peptides are involved in the process of tissue repair and remodeling .
The matrix metalloproteinases (MMPs) family is composed of more than 25 zinc-dependent proteases that cleave the extracellular matrix and cell-surface proteins to regulate wound healing, physiological angiogenesis and immune response . Matrix metalloproteinases can activate and increase the bioavailability of a variety of non matrix proteins, including cytokines, chemokines, recep-tors and antimicrobial peptides [5,9]. Matrilysin 1 (MMP-7), unlike many MMPs, is expressed by non injured, non inflamed mucosal epithelia in most adult human tissues (7). Besides extracellular matrix (ECM) components, MMP-7 processes cell surface molecules such as pro-defensin, Fas-ligand, pro-tumor necrosis factor (TNF), and E-cadherin [10,11].
The production of the enzyme is highly up-regulated by injury or exposure to bacteria, stimulating cell migration and coordinating the inflammatory response [12, 13, 14]. Thus, MMP-7 participates in the processes of defense, repair and inflammation.
In the promoter of the gene of MMP-7 an A>G transition has been identified as a functional polymorphism (MMP7 –181A>G; rs11568818). The polymorphism has been shown to modulate transcriptional activity by influencing the binding of nuclear regulatory proteins. The G allele has greater basal transcriptional activity than the A allele in vitro experiment in the human monocyte/ macrophage cell line U937 . The nuclear proteins bind with higher affinity to the G allele than to the A allele, resulting in significant increase in promoter activity [16,17].
So far, there are only very limited studies concerning the role of MMP-7 and its genetic variants in lung diseases such as lung cancer , pneumoconiosis , bron-chiolitis obliterans syndrome (BOS) in patients after lung transplantation , and idiopathic pulmonary fibrosis (IPF) . Concerning COPD, there was only one report showing increased serum level of MMP-7 in patients exposed to biomass and tobacco smoke compared with non smoking healthy controls, and the levels were negatively correlated with spirometric index of lung function (FEV1 %pr.) . However, there has only been one study exploring the functional variants of MMP-7 in the development of COPD .
In this respect we aimed to evaluate the possible role of the MMP 7 –181A>G promoter polymorphism as predisposing factor for COPD in a population from the region of Stara Zagora, Bulgaria.
Material and methods
Patients and Controls
We have genotyped 191 patients with COPD and 215 healthy volunteers or individuals unaffected by lung or cancer diseases. The inclusion criteria for COPD were as follows: age higher than 40 years; forced expiratory volume in one second (FEV1) of <80.0%; forced expiratory volume in one second (FEV1)/ forced vital capacity (FVC) ratio of ≤70.0%; FEV1 reversibility after inhalation of 400 μg Salbutamol of <12.0%.
In both groups, the age of inclusion in the study and smoking status were noted; in the patients’ group: age of diagnosis, the spirometric indexes, duration and the stages of the disease (GOLD stages) were also reported. The available demographic and clinical data are presented in Table 1. Informed consent was obtained from patients and controls before the beginning of the study.
Demographic and clinical data of patients with chronic obstructive pulmonary disease and controls.
|Characteristics||Patients with COPD (%)||Controls (%)||p Value|
|males||146 (76.4)||109 (50.7)||<0.001a|
|females||45 (23.6)||106 (49.3)|
|Age (years) at inclusion in the study median (range)||67 (36-88)||59 (18-80)||<0.001b|
|Age (years) at diagnosis of the disease median (range)||64 (30-86)||–|
|Smoking status||n = 186||n = 164|
|non-smokers||54 (29.0)||97 (59.0)||<0 001a|
|ex-smokers||89 (48.0)||25 (15.0)|
|current smokers||43 (23.0)||42 (26.0)|
|FEV1 %pr (mean ± SD) (range)||50.60 ± 14.10 (15-79)||93.40 ± 11.90 (82-114)||<0.001b|
|FEV1/FVC % (mean ± SD (range||61.1 ± 8.6 (27.1-70.0)||80.4 ± 7.2 (76.3-93.0)||<0.001b|
|COPD staging (GOLD 2009)|
|stage II (moderate)||102 (55.0%)||–|
|stage III (severe)||73 (39.0)||–|
|stage IV (very severe)||11 (67.0)||–|
DNA Isolation and Genotyping
Genomic DNA was isolated from 0.2 mL of whole blood using a commercial kit for isolation of genomic DNA from blood (GenElute™ Mammalian Genomic DNA Miniprep Kit, Sigma-Aldrich, St. Louis, MO, USA).
The genotyping for the MMP7–181A>G (rs11568818) was performed by the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP)-based method. The final volume of each reaction was 15 μL, containing 0.5 U Dream Taq Polymerase (Fermentas, Waltham, MA, USA), 1.5 μL 10 × PCR buffer (with 1.5 mM MgCl2), 0.6 μL dNTPs (Sigma-Aldrich) in a final concentration of 200 μM for each of the four dNTPs, 0.3 μL of each primer in concentration of 20 pmol/μL (MMP7F: 5’-TGG TAC CAT AAT GTC CTG AAT G-3’; MMP7R: 5’-TCG TTA TTG GCA GGA AGC ACA CAA TGA ATT- 3’) and distilled water to the end volume.
The temperature profile of the PCR reactions included primary denaturing of the template DNA for 3 min. at 94°C, followed by 30 cycles of denaturation for 30 seconds at 94°C, annealing for 30 seconds at 53.6°C and poly-merization for 30 seconds at 72°C. The PCR reaction was terminated by a final extension for 5 min. at 72°C.
The restriction reaction was performed with 5U EcoRI in final volume of 5 μL for 16 hours at 37°C. The obtained restriction products were analyzed by 4.0% agarose gel stained with ethidium bromide. The results were documented by the Gel documentation system (Syngene; Synoptics Ltd., Cambridge, Cambridgeshire, UK).
Statistical analyses were performed using the Statistical Package for the Social Sciences (SPSS), version 16.0 for Windows (SPSS Inc., Chicago, IL, USA). Continuous variables were analyzed for normality of the distribution using the Kolmogorov-Smirnov test (One-Sample Kolmogorov-Smirnov D-Test in SPSS, version 16; SPSS Inc.). When the level of signif-icance in this test was lower than 0.05 (p <0.05), the hypothesis for normal distribution was rejected. The continuous variables with normal distribution were com-pared between two or more independent groups by the Student t-test or one-way analysis of variance (ANOVA) test with least significant difference (LSD) post hoc ana-lysis, while those with an abnormal distribution were analyzed with the Mann-Whitney U or Kruskal-Wallis tests. The frequencies of distribution in the contingency tables were analyzed using χ2 test, or Fisher’s exact test, when needed. The odds ratio (OR) and 95% confidence interval (95% CI) were calculated by binary logistic regression with age and sex as covariates. The Hardy-Weinberg equilibrium (HWE) was calculated by an inter-active calculation tool for χ2 tests of goodness of fit and independence . Factors with a p value of <0.05 were considered to be statistically significant.
The PCR product amplified with the primers for MMP7 –181A>G was 150 bp in length. The EcoRI digestion resulted in 150 bp for the AA genotype, 150, 120 and 30 bp for the AG genotype, and 120 and 30 bp for the GG genotype (Figure 1).
The genotype distribution in both controls and patients did not deviate from HWE (p = 0.998 and p = 0.999). After statistical processing of the data, we found no differences in the genotype and allele distribution of the MMP7 –181A>G polymorphism between COPD patients and controls (p = 0.341 and p = 0.214) (Table 2). In the group of patients, there was no association between the genotype and gender (p = 0.358), smoking habit (p = 0.587), stage of the disease (p = 0.924) or spirometric indexes for lung function (for FEV1: p = 0.598 and for FEV1/FVC: p = 0.3, Kruskal-Wallis test). Association of the MMP-7 G allele with severity of COPD was also not found (p = 0.876).
Genotype and allele frequencies of the MMP7 –181A>G gene polymorphism in patients with chronic obstructive pulmonary disease and controls.
|MMP7 –181A>G||Patients||Controls||OR (95% CI)||p Value|
|n = 191||Frequency||n = 215||Frequency|
|–181 A||185||0.484||227||0.528||1.0 (referent)|
|–181 G||197||0.516||203||0.472||1.191 (0.904-1.569)||0.232|
However, the carriers of the G allele (AG and GG genotypes) appeared to develop COPD significantly early compared to those with the AA genotype (61.01 ± 10.11 vs. 64.87 ± 9.00 years, p = 0.032, Student t-test). When the genotype distribution was studied only in the patients (n = 76) and controls (n = 106) younger than 60 years, significant difference was found. The frequency of the carriers of the G allele (AG and GG genotypes), was higher in COPD patients than in controls, showing a 3-fold higher risk for the development of COPD before the age of 60 years (Table 3).
Genotype and allele and genotype frequencies of the MMP7 –181A>G gene polymorphism in the patients with chronic obstructive pulmonary disease and controls younger than 60 years.
|MMP7 –181A>G||Patients||Controls||OR (95% CI)||p Value||OR (95% CI)||p Value|
|n = 76||Frequency||n = 106||Frequency|
|AA||11||0.145||35||0.330||1.0 (referent)||1.0 (referent)|
|AG||43||0.566||50||0.471||2.736 (1.252-5.962)||0.016a||4.742 (1.426-15.768)||0.011a|
|GG||22||0.289||21||0.198||3.333 (1.363-8.148)||0.008a||3.194 (0.862-11.830)||0.082|
|AG+GG||65||0.855||71||0.670||2.913 (1.380-6.133)||0.005a||4.122 (1.318-12.890)||0.015a|
|–181 A||65||0.428||120||0.566||1.0 (referent)||–||–||–|
|–181 G||87||0.572||92||0.434||1.746 (1.147-2.657)||0.009a||–||–|
A common feature of COPD is the chronic in inflammation in the airways and the development of extensive tissue remodeling during the course of the disease process . The matrix metalloproteinases are a family of zinc-containing enzymes with proteolytic activity against a wide range of extracellular proteins . Under normal physiological conditions, the activities of MMPs are precisely regulated at the level of transcription, activation of the precursor zymogens and inhibition by endogenous inhibitors . Due to their activity, MMPs participate in many physiological and pathological processes in the body such as development, involution, inflammation, tumor growth, and repair [27, 28].
Matrix metalloproteinases play an important role in the turnover of almost all extracellular matrix molecules and thus, participate in the pathogenesis of COPD [22, 29]. An A>G substitution in the promoter of the MMP7 gene has been shown to affect the promoter activity, as the G allele determines higher basal transcriptional activity in vitro in human monocyte/macrophage cell line U937 . In this study, we did not find any significant difference in geno-type and allele distribution when we studied the entire COPD patient and control groups. However, the carries of the G allele (AG and GG genotypes) appeared to develop COPD significantly early compared to those with the AA genotype. Moreover, the G allele determines about a 3-fold higher risk for COPD before the age of 60 years.
So far, only one study has been found in the literature exploring the role of MMP7 –181A>G (rs11568818) in COPD. In the study of Mogulkoc et al. , the MMP-7 AA genotype has been found to be associated with an increased risk of COPD. On the contrary, in our study we found no association of the A allele with COPD, which might be explained by the difference in Bulgarian and Turkish populations (e.g., age, smoking, diagnosis, etc.).
It has been shown that the same variant is not a risk factor for IPF, but influences the plasma level of MMP-7 in patients, as carriers of the AA genotype had higher concentrations than carriers of other genotypes . MMP7 –181A>G has shown no association with the risk for lung cancer . However, an increased risk for devel-opment of bronchiolitis obliterans syndrome has been reported for carriers of the AA genotype in patients after lung transplantation .
In COPD patients, the serum levels of MMP-7, as well as of some other MMPs, MMP-1 and MMP-9, has been found to be significantly higher in exposed to biomass and tobacco smoke when compared with unexposed healthy controls. Moreover, the levels of those enzymes have been found in negative correlation with the lung function indexes (FEV1 %pr) .
Matrilysin 1 has also been reported to associate with moderate panlobular emphysema as well as with severe and moderate centrilobular emphysema . Matrilysin 1 efficiently cleaves the basement membrane protein entac-tin, which bridges laminin and collagen type IV, and sug-gested a potentially important role for MMP-7 in the dis-ruption of basement membranes by inflammatory cells . Matrilysin 1 is produced by the epithelium of several uninjured, uninflamed tissues, such as lung, liver and breast. Except in intact tissues, matrilysin is expressed in migrating epithelium in injured airways .
It has been shown that MMP-7 mediates shedding of E-cadherin from alveolar epithelium during progression of bleomycin-induced pulmonary fibrosis . In specimens of emphysema, strong immunoreactive signal for matrilysin protein has been detected in epithelial cells lining damaged alveoli, especially in cells bordering denuded epithelium . The higher promoter activity of the G allele may contribute to ongoing epithelial activation by mediating persistent shedding of the E-cadherin ectodomain, altering cell-cell interactions. By cleaving the Fas ligand, matrilysin can regulate apoptosis and it may promote local coagulation by cleaving tissue factor pathway inhibitor . Altogether, these data suggest that the carriers of the G allele may lead to disruption of ECM and cell-cell interactions, impairment of the process of tissue repair or fibrosis, and thus, to early development of COPD.
In conclusion, the results of our study suggest that MMP7 –181A>G (rs11568818) promoter polymorphism might affect the risk for COPD, as the carriers of the G allele (AG and GG genotypes) could be considered as predisposing factors for early onset of COPD. This effect could be attributed to the increased production of the enzyme resulting in enhanced airway wall protein degra-dation and injury due to its direct ability for ECM degra-dation. In this respect, the latter suggestion is warranted to be proven by analyzing the possible association of geno-types with serum MMP-7 levels in patients with COPD.
This study was performed with the support of projects 23/2015 MF and 4/2016 MF, Trakia University, Bulgaria and COST Action CA15129 “Diagnosis, Monitoring and Prevention of Exposure-Related Noncommunicable Diseases (DiMoPEx).”
Declaration of Interest. The authors report no conflicts of interest. The authors alone are responsible for the content and writing of this article.
GOLD 2017. Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease 2017 report. [Available from http://goldcopd.org/gold-2017-global-strategy-diagnosis-management-prevention-copd/.]
Echave-Sustaeta JM, Comeche Casanova L, Cosio BG, Soler-Cataluna JJ, Garcia-Lujan R, Ribera X. Comor-bidity in chronic obstructive pulmonary disease. Related to disease severity? Int J Chron Obstruct Pulmon Dis. 2014; 9: 1307-1314.
Parks WC. Matrix metalloproteinases in lung repair. Eur Respir J Suppl. 2003; 44: 36s-8s.
Han JC, Li XD, Du J, Xu F, Wei YJ, Li HB, et al. Elevated matrix metalloproteinase-7 expression promotes metastasis in human lung carcinoma. World J Surg Oncol. 2015; 13(5): 1477-7819.
Jormsjo S, Whatling C, Walter DH, Zeiher AM, Hamsten A, Eriksson P. Allele-specific regulation of matrix metalloproteinase-7 promoter activity is associated with coronary artery luminal dimensions among hyper-cholesterolemic patients. Arterioscler Thromb Vasc Biol. 2001; 21(11): 1834-1839.
Liu D, Guo H, Li Y, Xu X, Yang K, Bai Y. Association between polymorphisms in the promoter regions of matrix metalloproteinases (MMPs) and risk of cancer metastasis: A meta-analysis. PLoS One. 2012;7(2):14.
Kastelijn EA, van Moorsel CH, Ruven HJ, Karthaus V, Kwakkel-van Erp JM, van de Graaf EA, et al. Genetic polymorphisms in MMP7 and reduced serum levels associate with the development of bronchiolitis obliterans syndrome after lung transplantation. J Heart Lung Transplant. 2010; 29(6):680-686.
Richards TJ, Park C, Chen Y, Gibson KF, Peter Di Y, Pardo A, et al. Allele-specific transactivation of matrix metalloproteinase 7 by FOXA2 and correlation with plasma levels in idiopathic pulmonary fibrosis. Am J Physiol Lung Cell Mol Physiol. 2012; 302(8): L746-L754.
Lu L, Gunja-Smith Z, Woessner JF, Ursell PC, Nissen T, Galardy RE, et al. Matrix metalloproteinases and collagen ultrastructure in moderate myocardial ischemia and reperfusion in vivo. Am J Physiol Heart Circ Physiol. 2000; 279(2): H601-H609.
Brajer B, Batura-Gabryel H, Nowicka A, Kuznar-Kaminska B, Szczepanik A. Concentration of matrix metalloproteinase-9 in serum of patients with chronic obstructive pulmonary disease and a degree of airway obstruction and disease progression. J Physiol Pharmacol. 2008; 59(Suppl 6): 145-152.
Manjari KS, Jyothy A, Kumar PS, Prabhakar B, Nallari P, Venkateshwari A. Association of matrix metalloproteinase-7 (-181A/G) promoter polymorphism in chronic pancreatitis. Indian J Med Res. 2014; 140(5): 609-615.
Belaaouaj AA, Li A, Wun TC, Welgus HG, Shapiro SD. Matrix metalloproteinases cleave tissue factor pathway inhibitor. Effects on coagulation. J Biol Chem. 2000; 275(35): 27123-27128.
Determined by the χ2 test.
Determined by the χ2 test.