Abstract
Alveolar macrophages play an important role in chronic obstructive pulmonary disease via production of matrix metalloproteinases (MMPs) and cathepsins as well as their inhibitors, tissue inhibitors of metalloproteinases and cystatin C. We hypothesised that expression levels of these molecules by alveolar macrophages at baseline and after stimulation would be influenced by genotype and associated with chronic obstructive pulmonary disease phenotypes.
Quantitative PCR and ELISAs/gelatine zymography were used to investigate expression levels of mRNA and protein, respectively. The relationships of expression with genotype, pulmonary function and emphysema were analysed.
The results showed that basal expression level of MMP12 mRNA was inversely related to the diffusing capacity of the lung for carbon monoxide/alveolar volume and to forced expiratory volume in 1 s/forced vital capacity after correction for multiple comparisons. The expression level of MMP12 protein stimulated with lipopolysaccharide was also inversely related to the diffusing capacity of the lung for carbon monoxide/alveolar volume and was positively related to the extent of emphysema. The basal expression of MMP1 mRNA was positively correlated with the extent of emphysema. Cathepsin L protein level was positively associated with forced expiratory volume in 1 s % predicted.
We conclude that increased MMP12 and MMP1 expression may play a role in the pathogenesis of emphysema. Cathepsin L and MMP9 may be involved in the development of airflow limitation.
Abstract
MMP12 expression may play a role in the pathogenesis of emphysema and airflow limitation http://ow.ly/pICmy
Introduction
Chronic obstructive pulmonary disease (COPD) is characterised by airflow limitation that is not fully reversible and is a major cause of chronic morbidity and mortality throughout the world [1]. The hallmark of this disease is an inflammatory process within the lungs in response to noxious particles originating mainly from cigarette smoke, with concomitant destruction and repair of the lung architecture, leading to tissue remodelling. Although the major trigger for the development of COPD is cigarette smoke, only a minority of chronic heavy smokers develop symptomatic COPD [2]. The clustering of COPD in families suggests that there is a genetic component [3] and may explain why only a subset of smokers develops this disease.
A proteinase:antiproteinase imbalance is thought to play a central role in the pathogenesis of COPD [4]. An important source of both proteinases and antiproteinases within the lung is the alveolar macrophage, which produces matrix metalloproteinases (MMPs) and their inhibitors, tissue inhibitors of metalloproteinases (TIMPs), as well as cathepsins and a potent inhibitor of cathepsins, cystatin C (CST3).
MMPs are a family of proteolytic enzymes that degrade collagen or elastin and play an essential role in tissue remodelling and repair associated with inflammation [5]. MMP1, MMP9 and MMP12 are produced by alveolar macrophages, and several studies in animals and humans have provided evidence that these MMPs are important in airway inflammation and the development of emphysema [6, 7]. Cathepsin L (CTSL) and cathepsin S (CTSS) are lysosomal cysteine proteinases produced by alveolar macrophages [6] that may also be involved in COPD through degradation of elastin [8].
As chronic inflammatory processes are thought to be important in the pathology of COPD, cytokines such as interleukin (IL)-1β and tumour necrosis factor (TNF)-α may play a critical role in the disease process [9]. Microbial infections induce exacerbations through bronchial inflammation in COPD patients, which causes tissue damage and could contribute to the development of emphysema [10]. Thus, it is important to also examine the expression change of proteinases and antiproteinases induced by inflammatory cytokines or lipopolysaccharide (LPS), a potent secretagogue for a variety of inflammatory cytokines [10].
We hypothesised that higher expression levels (basal and stimulated by IL-1β, TNF-α or LPS) of proteinases (MMP1, MMP9, MMP12, CTSL and CTSS) and/or lower levels of antiproteinases (TIMP1, TIMP2, TIMP3 and CST3) in human alveolar macrophages would be associated with lower lung function (airflow limitation and diffusing capacity of the lung for carbon monoxide (DLCO)) and greater extent of emphysema. We also hypothesised that polymorphisms of these genes would alter the level of expression of their mRNA and protein.
Materials and methods
Study subjects
109 subjects who underwent lung resection for small (<3 cm diameter), TNM (tumour, node, metastasis) stage I or II, peripheral tumours at Vancouver General Hospital (Vancouver, BC, Canada) were recruited for this study. Partial results from 54 of these subjects were included in a previous study [11]. Preoperative pulmonary function tests were performed including forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and single breath DLCO/alveolar volume (VA), in accordance with American Thoracic Society guidelines [12]. Computed tomography (CT) scans were available for assessment of emphysema in 75 subjects. The characteristics of the subjects are shown in table 1.
This study was approved by the University of British Columbia/Providence Health Care and Vancouver Hospital and Health Sciences Research Ethics Boards (Vancouver, BC, Canada). All subjects provided written informed consent for the study.
Computed tomographic analysis
Details of a semi-quantitative emphysema visual score, and a quantitative computerised analysis of the percentage of lung with emphysema, are given in the online supplementary materials and methods.
Bronchoalveolar lavage and cell culture
The resected lung or lobe was obtained immediately post-operatively and bronchoalveolar lavage was performed. Details of the lavage, macrophage separation and subsequent culture are given in the online supplementary materials and methods.
RNA extraction, cDNA synthesis and quantitative real-time PCR
Details of cDNA preparation and quantification by real-time PCR (online supplementary table S1) are given in the online supplementary materials and methods.
Protein quantification
Details of the protein quantification by ELISAs or zymography are given in the online supplementary materials and methods.
Genotyping
We selected polymorphisms (online supplementary table S2) in MMPs, cathepsins, TIMPs and CST3 that have been previously associated with COPD and/or with regulation of gene expression (rs1799750 in MMP1, D20S838 and rs3918242 in MMP9, rs2276109 in MMP12, and rs2277698 in TIMP2 [13–17]. In addition, we selected single nucleotide polymorphisms (SNPs) from the HapMap Project (www.hapmap.org) to identify a haplotype block covering the promoter region of each gene (details are given in the online supplementary materials and methods).
Statistical analysis
Statistical analyses were performed using JMP software version 5.1 (SAS Institute Inc., Cary, NC, USA). Linear regression analysis was used to evaluate mRNA or protein expression versus pulmonary function tests and quantitative CT scan data. We also used linear regression to evaluate genotype versus mRNA or protein expression with an additive model, and p-values were calculated in the linear model based on a Wald test statistic. Ordinal regression was used to examine the association between semiquantitative CT scores and gene expression levels. Adjustment was performed for age, sex, race and smoking history (pack-years). If the values were not normally distributed, they were log-transformed before analyses.
p-values <0.05 were considered nominally significant. Due to the nature of the analysis and the multiple-testing of independent and dependent variables, modified Bonferroni correction for multiple comparisons was done as described by Nyholt [18] and details of this correction procedure are provided in the online supplementary tables S3 and S4. The power of this study design was estimated using the many proportions function and linear regression function in PASS 2005 [19].
Results
Proteinase and antiproteinase mRNA expression
The mRNA expression levels of the proteinases (MMP1, MMP9, MMP12, CTSL and CTSS) and antiproteinases (TIMP1, TIMP2, TIMP3 and CST3) in alveolar macrophages under the different experimental conditions (basal, i.e. uncultured macrophages, cultured without stimulation, cultured with stimulation by IL-1β, TNF-α or LPS) are shown in online supplementary figures S1–S9. The alveolar macrophages cultured without stimulation showed a significant increase in expression of all of the genes (except TIMP3 and CST3) relative to the basal gene expression level. With LPS stimulation, the expression levels of all the proteases (except CTSS), as well as TIMP1, were significantly increased. Conversely, the expression levels of TIMP2 and CST3 were significantly decreased after LPS stimulation. Following stimulation with either IL-1β or TNF-α the expression levels of MMP1, MMP9 and TIMP1 were significantly increased. Stimulation with either IL-1β or TNF-α resulted in significantly decreased TIMP2 expression.
There were no significant effects of age, sex or race on the basal expression of the genes. The effect of smoking status on basal gene expression is shown in online supplementary table S5. MMP12 basal expression was higher and TIMP2 basal expression was lower in current smokers compared with former smokers. TIMP3 basal expression was positively correlated with pack-years while CTSS expression was inversely correlated with pack-years.
Proteinase and antiproteinase genotypes versus basal and stimulated gene expression
As genetic variation could affect COPD phenotypes through modulation of the proteinase and antiproteinase mRNA and protein expression profiles, we investigated whether polymorphisms altered expression levels of their respective gene products. All of the polymorphisms we investigated were in Hardy–Weinberg equilibrium (p≥0.18) in the Caucasian subgroup (n=87), as shown in online supplementary table S2.
A summary of the results of association between genotypes and basal gene expression and expression change is shown in online supplementary table S6. MMP1 mRNA basal expression and protein expression induced by LPS were both inversely associated with increasing number of the 2G allele of MMP1 rs1799750 (online supplementary figs S10 and S11, respectively). CTSL protein expression without stimulation was inversely associated with increasing number of the T allele of CTSL rs2274611 (online supplementary fig. S12). The most significant association of genotype with gene expression was observed for the basal expression of CST3 mRNA, which was inversely associated with increasing number of the A allele of CST3 rs6036478 (online supplementary fig. S13). All of these associations remained significant after correction for multiple comparisons.
Proteinase and antiproteinase mRNA expression and COPD phenotypes
We analysed the relationship of COPD phenotypes and basal mRNA expression levels using linear regression analyses (table 2 and online supplementary table S7). The basal expression level of MMP12 in alveolar macrophages showed a significant inverse association with DLCO/VA % predicted that remained significant after correction for multiple comparisons. A scatter plot of the association between MMP12 basal mRNA expression and DLCO/VA % pred is shown in figure 1. There was also an inverse association between basal MMP12 and FEV1/FVC % (fig. 2). The data also suggest that there was a positive association (p=0.0052, uncorrected for multiple comparisons) between basal expression of MMP1 and the semi-quantitative visual assessment of extent of emphysema of all CT scans (total n=75); this association was no longer significant after correction for multiple comparisons, as the level of significance required was p<0.0012. There was no significant association (uncorrected p=0.120) with the computerised quantitative assessment of emphysema by CT scanning (n=45), which may be due to the lower number of subjects with these data available.
We also tested for associations between changes in mRNA expression induced by LPS, TNF-α and IL-1β and COPD phenotypes (table 3 and online supplementary table S8). The most significant result was that expression change of MMP9 mRNA stimulated by TNF-α was inversely associated with FEV1 % pred (p=0.0005), but this did not survive correction for multiple comparisons, since the level of significance required was p<0.00042.
The relationship between mRNA and protein expression
As we postulated that gene expression would be associated with COPD phenotypes through the level of protein expression, we investigated the association between mRNA and protein expression for the genes that showed at least nominally significant associations with the clinical outcomes. We compared the mRNA levels at 24 h with the protein concentrations of MMP1, MMP9, MMP12, CTSL, TIMP2 and CST3 in media after 48 h of culture. The levels of protein are shown in online supplementary table S9. We observed significant correlations between the mRNA and protein levels for MMP1 stimulated with LPS, MMP9 with and without LPS, and MMP12 with and without LPS (online supplementary table S10). The correlations between MMP12 mRNA and protein expression of MMP12 in alveolar macrophages are shown in figure 3. The basal mRNA expression of MMP12 was also positively correlated with its protein expression level stimulated with LPS (p=1×10−5).
We also investigated the correlation between the fold change in mRNA and protein levels after culture with the various stimuli (online supplementary table S11). Significant relationships were observed for changes in MMP1 level after stimulation with LPS and IL-β. The changes in expression of MMP9 mRNA and protein after culture with all the stimuli were also significantly correlated.
Protein expression of proteinases and antiproteinases and COPD phenotypes
The protein expression of the proteinases and antiproteinases was investigated in relation to COPD phenotypes (table 4 and online supplementary table S12). The protein expression of MMP12 stimulated with LPS was inversely associated with DLCO/VA % predicted and positively associated with quantitative score for emphysema. The protein expression of CTSL with or without LPS was also positively associated with FEV1 % predicted. However, none of these associations with protein expression were significant after correction for multiple comparisons.
Discussion
In this study, we found that the genotypes of MMP1, CTSL and CST3 showed significant associations with their expression at the mRNA and/or protein level. The MMP1 1G/2G polymorphism (rs1799750) has previously been associated with gene expression [20], but this is the first description of an effect on transcription in alveolar macrophages. There are no previous studies of the rs2274611 CTSL polymorphism, but a CTSL promoter SNP (rs3118869) has been associated with reporter gene expression in vitro [21]. The C allele of rs3118869 showed greater expression than the A allele (p=3.36×10−6) in luciferase assays and was associated (p=0.0143) with a higher level of CTSL protein in our study (online supplementary table S6). The C allele of rs3118869 is in high linkage disequilibrium (r2=0.8648) with the C allele of rs2274611, which was associated (p=0.0007) with higher CTSL protein in our study (online supplementary table S6). Thus, this is the first study to confirm the functional effect of rs3118869 on transcription and we also demonstrated a novel association with CTSL protein level. There are no previous studies of the CST3 rs6036478 polymorphism in the literature, but this SNP is part of a haplotype of SNPs in high linkage disequilibrium (all r2>0.92), which has been shown to influence serum/plasma CST3 levels [22–25]. Our observation that the A (minor) allele of rs6036478 was associated with lower basal CST3 mRNA (online supplementary table S6) is consistent with these previous data and is novel in the context of COPD.
The basal level of MMP12 mRNA expression was inversely related to DLCO/VA % pred and FEV1/FVC % (table 2). MMP12 protein level after stimulation with LPS was also inversely related to DLCO/VA % pred and was positively associated with the quantitative score for emphysema (table 4 and online supplementary table S12). Although these latter associations were not significant after correction for multiple comparisons, the concordance of direction of association suggests a true pathophysiological role for this proteinase.
Our novel findings of associations between MMP12 expression and lung diffusing capacity, as well as the extent of emphysema, suggest that MMP12 has a role in the development of emphysema in humans. It has been reported that genetic variation in MMP12 is associated with COPD and asthma [26]. However, there are no reports on the direct association between emphysema and the expression of MMP12 from alveolar macrophages in humans. MMP12 knockout mice do not develop emphysema following exposure to cigarette smoke [27], suggesting that the presence of MMP12 is critical in smoke-induced emphysema in mice. Our study further suggests a role for MMP12 in the pathogenesis of emphysema in humans. The underlying mechanism for this role may relate to MMP12 stimulation of neutrophils through activation of TNF-α, resulting in the production of MMP9, neutrophil elastase and cathepsin G, all of which degrade elastin and may lead to emphysema [28]. In addition, MMP9 and MMP12 degrade specific collagens [29, 30]. MMP12 is increased in the lungs [31] and sputum [32] of COPD patients.
We also showed that the basal level of MMP1 mRNA was positively correlated with a semiquantitative assessment of the extent of emphysema (p=0.0052), as in our previous report, based on the first 44 cases with CT scans [11]. Thus, we have extended our previous observations by demonstrating that not only MMP12 mRNA but also protein level is associated with DLCO/VA and emphysema.
Cigarette smoke and other irritants stimulate alveolar macrophages to produce MMPs [33, 34] and cathepsins [35]. MMP1 degrades collagen whereas MMP9, MMP12, CTSL and CTSS degrade elastin [36]. Although human MMP1 is a collagenase, its expression in the lungs of transgenic mice caused emphysema [37] and expression of MMP1 was increased in the lungs of patients with emphysema [38]. Smokers who have airway obstruction show increased expression of MMP1 and MMP9 compared with smokers without COPD and nonsmokers [39].
We observed a novel association of TNF-α induced MMP9 mRNA expression with FEV1 % pred (table 3). This observation is consistent with previous data that demonstrated MMP9 expression increased in the lung parenchyma surrounding small airways as FEV1 % pred decreased [40]. There is also supportive evidence that MMP9 contributes to the development of emphysema from mouse models [41, 42]. Our data suggest that it is the level of MMP9 in response to specific pro-inflammatory stimuli that plays a critical role in the development of COPD, in agreement with studies showing that increases in TNF-α are associated with reductions in lung function [40] and that TNF-α plays an important role in cigarette smoke-induced emphysema [43].
CTSL protein level was associated with FEV1 % pred, which suggests that this proteinase has a role in development of COPD. However, individuals whose macrophages secreted more CTSL had better lung function. We did not observe this relationship in our previous study [11] but the same trends were present (i.e. higher CTSL with or without LPS associated with higher lung function) with similar estimates from the regression analysis. However, the p-values were not significant due to the smaller sample size. It is unclear what mechanism underlies this “paradoxical” relationship. However, CTSL expression in thymocytes is essential for natural killer T-cell development [44], and natural killer cells, which have an essential role in immunity against viral infection [45], are decreased in the peripheral blood in COPD [46]. Similarly, CTSL is involved in antigen processing [47, 48] and T-helper cell development [49, 50]. Thus, CTSL protein expression could reinforce immunological capacity, reduce exacerbations and prevent progression of COPD. CSTL also plays a prominent role in the production of neuropeptides [51], such as neuropeptide Y, which is decreased in the airway epithelium in COPD patients compared with unaffected smokers [52], suggesting a protective role for this molecule. In addition, CSTL is involved in the production of adrenocorticotropic hormone, and upregulation of this hormone is attenuated in COPD patients following exercise [53].
In this study, there were significant differences in gene expression in macrophages cultured without stimulation compared with primary macrophages, as we have previously shown [11]. This phenomenon has also been observed in murine macrophages [54] and other cell types following culture [55, 56]. While it was not surprising that nonphysiological conditions (i.e. cell culture) altered gene expression in our study, whether this was due to adherence of the gcells to the plastic substrate in the cell culture dishes or to components of the culture medium is not known.
A recent study investigated changes in gene expression following instillation of LPS into the lung in healthy individuals (n=8) compared with instillation of saline into the contralateral lung [57]. In that in vivo study [57], TIMP1 was upregulated 2.4-fold (p=3×10−4) following challenge with LPS, mirroring the increase in TIMP1 mRNA that we observed after culture with LPS (online supplementary fig. S6). Similarly, TIMP2 mRNA was downregulated following instillation of LPS (p=3×10−3), and we also observed a significant decrease in TIMP2 mRNA following culture with LPS (online supplementary fig. S7).
In comparison with our previous study [11], we have doubled the sample size to increase the power of the analyses and have assessed additional genes (i.e. the antiproteinases TIMP1, TIMP2, TIMP3 and CST3). We also measured the levels of the proteins corresponding to six of the genes whereas our previous manuscript only included mRNA data [11].
There are some limitations of our study with respect to gene expression and COPD or emphysema phenotypes. With a total sample of 109 subjects we have 80% power to detect associations between gene expression and COPD phenotypes only where the r2 value is ≥0.3. As we have a smaller number of subjects for DLCO/VA % pred and the CT scan data, the power associated with these phenotypes is less. In addition, there was a small proportion of severe COPD patients in our study (59% of our subjects had FEV1 values within the normal range). In addition, the use of alveolar macrophages from patients with lung cancer adds a potential bias to our study. Specifically, several MMP polymorphisms have been associated with lung cancer [58] and, therefore, the risk alleles may have been selected for in our study group. However, all SNPs in our study showed good agreement with expectations under Hardy–Weinberg equilibrium (online supplemetary table S2) suggesting that there was no strong selection bias. In addition, the allele frequencies are generally similar to the HapMap data (online supplemetary table S2). Expression levels of proteinase and antiproteinase genes may also have been affected by the presence of tumours in the lungs. However, as noted above, our results are similar to those of Reynier et al. [57] in subjects without cancer, suggesting that the observations in this study are applicable to noncancer populations.
In conclusion, among proteinases and antiproteinases, our data further support previous studies demonstrating that MMP12 and MMP1 may have an important role in the pathogenesis of emphysema. In addition, CTSL may positively influence the development of airflow limitation through protein expression and is affected by genetic variation.
Acknowledgments
The authors gratefully acknowledge the technical assistance of Hong Li, Rosa Garcia, Beth Whalen, Sebastian Cogswell, Barbara Moore, Mark Elliot and Steve Kalloger (University of British Columbia, Vancouver, BC, Canada). We also gratefully acknowledge Ken Evans and John Yee (University of British Columbia) for providing us with lung specimens from their patients undergoing lung resection surgery, as well as Paola Nasute Fauerbach and Nestor Müller (University of British Columbia) for assigning the semiquantitative emphysema scores to the CT scans.
Peter D. Paré is a Michael Smith Distinguished Scholar and the Jacob Churg Distinguished Researcher.
Footnotes
This article has supplementary material available from www.erj.ersjournals.com
Support statement: This study was supported by a grant from the Canadian Institutes of Health Research. A.J. Sandford is the recipient of a Michael Smith Foundation for Health Research Senior Scholarship award (grant number CI-SSH-01874(07-1)). T. Ishii was the recipient of a Pfizer Fellowship from the Japanese Respiratory Society, and fellowship support from the BC Lung Association, and the CIHR/HSFC IMPACT Strategic Training Program.
Conflict of interest: None declared.
- Received October 30, 2012.
- Accepted April 30, 2013.
- ©ERS 2014