- Original investigation
- Open Access
Plasma matrix metalloproteinases are associated with incident cardiovascular disease and all-cause mortality in patients with type 1 diabetes: a 12-year follow-up study
Cardiovascular Diabetology volume 16, Article number: 55 (2017)
Altered regulation of extracellular matrix remodeling by matrix metalloproteinases (MMPs) and tissue inhibitor of metalloproteinase (TIMP) may contribute to vascular complications in type 1 diabetes. We investigated associations between plasma MMP-1, -2, -3, -9, -10 and TIMP-1, and cardiovascular events and all-cause mortality in type 1 diabetic patients.
We prospectively followed 337 type 1 diabetic patients [mean age 41.4 years (9.6), 39% female], 170 with and 167 without diabetic nephropathy, with median follow-up of 12.3 years. Survival analyses were applied to investigate differences in plasma MMP-1, -2, -3, -9, -10, and TIMP-1-levels in patients with and without a cardiovascular event and in those who died vs survivors. All analyses were adjusted for age, sex, duration of diabetes, HbA1c, nephropathy and for other conventional cardiovascular risk factors.
After adjustment for potential confounders, higher MMP-2 plasma levels were significantly associated with higher incidence of cardiovascular events [HR 1.49 (95% CI 1.11; 1.99)], and higher plasma levels of MMP-1 [1.38 (1.07; 1.78)], MMP-2 [1.60 (1.19; 2.15)] and MMP-3 [1.39 (1.05; 1.85)] were associated with all-cause mortality. All associations were independent of low-grade inflammation and endothelial dysfunction as estimated by plasma markers. Associations between MMP-2 and cardiovascular events and between MMP-3 and mortality were attenuated after further adjustment for eGFR and changes in eGFR.
Higher levels of MMP-2 are associated with CVD and higher MMP-1, -2 and -3 with all-cause mortality. In addition, associations between MMP-2 and CVD, and MMP-3 and mortality were attenuated after adjustment for eGFR while both MMPs were associated with eGFR decline, indicating a possible mediating role of eGFR.
The pathophysiological mechanisms leading to cardiovascular disease (CVD) in type 1 diabetes  have only been partly elucidated. Recent data suggest that altered regulation of extracellular matrix (ECM) remodeling by matrix metalloproteinases (MMPs) could play a role . MMPs are proteases, which degrade ECM components in normal and pathological conditions. Increased plasma MMP-2, for example, has been associated with macrovascular disease by increasing vascular remodeling as well as thrombus formation [3, 4]. MMP production may be enhanced by hyperglycemia, pro-inflammatory mediators, reactive oxygen species and aldosterone ; MMP activity is dependent on the number of MMPs and their inhibitors, such as tissue inhibitor of metalloproteinase (TIMP) and α2-macroglobulin . In turn, increased activity of MMP-2 and MMP-9 is associated with low-grade inflammation (LGI)  and endothelial dysfunction (ED) , and increased MMP-9 levels have been associated with renal dysfunction . Therefore, LGI, ED and renal dysfunction constitute potential pathophysiological mechanisms through which MMPs could lead to CVD.
Interestingly, higher plasma levels of MMP-2 and -10 have been observed in patients with type 1 diabetes compared to non-diabetic controls [9, 10]. In cross-sectional studies, individuals with type 1 diabetes and retinopathy or nephropathy showed significantly higher circulating levels of MMP-9 and -10 compared to individuals without these complications [10, 11]. These data support the concept that MMPs may play an important role in the development of CVD. In addition, several prospective studies in (mainly) non-diabetic individuals have shown independent associations between plasma MMP-2 , MMP-9  and TIMP-1  and CVD, as well as associations between MMP-1 , MMP-2 [12, 16, 17] and -3 [18, 19] and all-cause mortality. However, to our knowledge, no prospective studies have been reported on the associations between plasma MMPs or TIMP levels and CVD or all-cause mortality in type 1 diabetes.
In view of these considerations, we investigated associations, in type 1 diabetes, between plasma MMP-1, -2, -3, -9 and -10, and TIMP-1 and incident non-fatal and fatal cardiovascular events, as well as all-cause mortality, in a prospective study with more than 12 years of follow-up. In addition, we investigated the extent to which LGI, ED and renal dysfunction contributed to these associations.
Study population and design
In 1993, 199 type 1 diabetic patients with diabetic nephropathy and older than 18 years were enrolled in a prospective observational study at the outpatient clinic at Steno Diabetes Center. Diabetic nephropathy was defined as persistent macroalbuminuric [urinary albumin excretion (UAE) >300 mg/24 h] in at least two out of three previous consecutive 24-h urine collections, presence of retinopathy, and absence of other kidney or urinary tract disease. In addition, 192 normoalbuminuric diabetic patients (UAE <30 mg/24 h), matched for age, sex, and duration of diabetes, were enrolled for comparison .
Patients with prior CVD (n = 24), who were lost to follow-up (n = 17), with end-stage renal disease at baseline (n = 7) and those with missing data on MMPs/TIMP-1 (n = 6) were excluded from the present analysis resulting in a total study population of 337 individuals. Patient selection and inclusion criteria have been described in detail elsewhere . The study was approved by the local ethics committee, in accordance with the Helsinki Declaration, and all patients gave their written informed consent.
Commercially available ELISA kits were used to measure levels of MMP-1, MMP-2, MMP-3, MMP-9, MMP-10 and TIMP-1 in EDTA plasma samples [Human MMP 3-Plex Kit (for MMP-1, -3 and -9), Human MMP-2-Plex Kit (for MMP-2 and -10) and Human TIMP-1 Kit, MSD, Rockville, USA], according to the manufacturer’s protocol. The MMPs were detected in both pro- and active form. TIMP-1 was detected in the active form. The intra- and inter-assay coefficients of variation were 6.6 and 9.7% for MMP-1, 3.7 and 7.1% for MMP-2, 7.2 and 14.1% for MMP-3, 4.6 and 13.7% for MMP-9, 2.9 and 7.6% for MMP-10, and 5.1 and 7.4% for TIMP-1, respectively.
Other baseline measurements
Measurements of other biomarkers and risk factors have been described in detail elsewhere . All clinical investigations were performed in the morning after an overnight fast. In 88% of the normoalbuminuric and 24% of macroalbuminuric patients antihypertensive medication was never prescribed. Patients using antihypertensive medication were asked to stop their antihypertensive and diuretic treatment eight days before the examination, but 15% of the 337 patients (27% in the nephropathy group) considered in our analyses did not comply with this recommendation.
Arterial blood pressure was measured twice, with an appropriate cuff size, following at least a 10-min rest in supine position. Mean arterial pressure (MAP) was calculated as [systolic blood pressure + (2*diastolic blood pressure)]/3. Body mass index (BMI) was defined as weight (kg) divided by height (m) squared. An ELISA was used to measure UAE from 24-h urine collections. Serum creatinine was determined by a kinetic Jaffé method and we estimated the glomerular filtration rate (eGFR) using the chronic kidney disease epidemiology Collaboration (CKD-EPI) equation , after creatinine was IDMS standardized . Patients were divided into three groups according to their smoking status: never, former or current smokers.
An in-house ELISA was used to determine high-sensitivity C-reactive protein (CRP)  and commercially available ELISA kits were used to measure plasma levels of secreted phospholipase A2 (sPLA2), interleukin-6 (IL-6), soluble vascular cell adhesion molecule-1 (sVCAM-1) and soluble intercellular adhesion molecule-1 (sICAM-1).
Follow-up and study end-points
All patients were followed up to the last visit at Steno Diabetes Center, until September 1st 2006, death (n = 83) or emigration (n = 3). All patients were traced through the national register during autumn 2006. When a patient died before 1 September 2006, the date of death was recorded and the primary cause of death was obtained from the death certificate, which was reviewed by two independent observers. Additionally, information from necropsy reports was included, when available. All deaths were classified as cardiovascular unless an unequivocal non-cardiovascular cause was established. At the end of follow up, non-fatal events were retrieved from patients’ records at Steno Diabetes Center or other hospital records. The primary endpoint of this study was fatal or non-fatal CVD (i.e. myocardial infarction, percutaneous coronary intervention, coronary bypass grafting, amputation due to ischemia, vascular surgery for peripheral atherosclerotic disease and stroke). Ischemic heart disease was examined through a WHO questionnaire combined with an ECG . In several patients, who were alive at follow-up and had more than one cardiovascular event during follow-up, the first event was used in the analyses. All-cause mortality was the secondary end-point.
All analyses were performed with SPSS version 20 for Windows (SPSS, Chicago, IL, USA). Triglycerides, UAE, CRP, IL-6, sPLA2, MMP-1, MMP-2, MMP-3, MMP-9, MMP-10 and TIMP-1 showed a skewed distribution and were log℮ transformed prior to further analyses. Student’s t or Chi square tests were performed for comparisons of baseline characteristics between groups, when appropriate. Cox proportional hazards regression models were applied to investigate associations between lnMMP-1, lnMMP-2, lnMMP-3, lnMMP-9, lnMMP-10 and lnTIMP-1 and study endpoints with adjustment for sex, age, duration of diabetes, HbA1c and nephropathy-no nephropathy status (model 1). Further adjustments were performed for BMI, MAP, total cholesterol, smoking status, antihypertensive treatment and whether or not patients withheld their medication prior to baseline examinations (model 2). We used z-scores (in SD) of lnMMPs and lnTIMP-1 to enable direct comparisons between the strengths of the associations (if any) between the various markers and outcome measures.
An LGI score was calculated as the z-score of the average of the z-scores of lnIL-6, lnCRP, sICAM-1, and lnsPLA2 [2, 25,26,27]. Similarly, an ED score was calculated as the z-score of the average of the z-scores of sICAM-1 and sVCAM-1 [2, 26, 27]. We examined the cross-sectional associations between lnMMP-1, lnMMP-2, lnMMP-3, lnMMP-9, lnMMP-10 and lnTIMP-1, and renal function (i.e. eGFR and ln-UAE), the LGI score and the ED score using linear regression analyses. Further analyses with adjustment of the associations between MMPs and TIMP-1 and study endpoints for eGFR, and ln-UAE (model 3), LGI (model 4) and ED (model 5) were performed to evaluate the extent to which these pathophysiological processes could explain the associations between MMPs/TIMP-1 and study outcomes. Additionally, in case renal function does explain part of the association between circulating MMPs and TIMP-1 and outcome measures, the possible additional influence of decrease in eGFR over time will also be investigated. Finally, we investigated whether or not any such associations differed between patients with and without nephropathy or between men and women by adding interaction terms between nephropathy or sex and MMPs or TIMP-1 to the models. Whenever a p value <0.1 was observed, results were presented for both groups separately.
Additional file 1: Figure S1 depicts the flowchart concerning patient selection. During a median follow-up of 12.3 years (interquartile range 7.7–12.5), 86 patients had a non-fatal (n = 54) or fatal (n = 49) cardiovascular event. Eighty-three patients died from various causes: cardiovascular (n = 49), diabetes (n = 4), infection (n = 13), cancer (n = 6), suicide (n = 4) and other (n = 7). Patients with a cardiovascular event or who died during follow-up were significantly older, had a longer duration of diabetes, higher levels of HbA1c and a higher prevalence of other cardiovascular risk factors (Table 1). Additionally, MMP-1, MMP-2, MMP-3, MMP-9, MMP-10 and TIMP-1 were higher in this group compared to patients without cardiovascular events or surviving the period of follow-up. Additional file 1: Table S1 shows the baseline characteristics according to nephropathy-no nephropathy status. Additional file 1: Table S2 shows the correlations between five MMPs and TIMP-1.
Associations of MMPs and TIMP-1 with incident CVD
Additional file 1: Figure S2 shows the cumulative hazard plots of cardiovascular events according to tertiles of lnMMPs and lnTIMP-1. After adjustment for age, sex, duration of diabetes, HbA1c and nephropathy, higher levels of MMPs 1, 2, 3 and 9, and TIMP-1 were significantly associated with incidence of cardiovascular events, with a hazard ratio (HR) of 1.51 (95% CI 1.20; 1.91) per 1 SD higher lnMMP-1; 1.65 (1.28; 2.13) for lnMMP-2; 1.47 (1.13; 1.92) for lnMMP-3; 1.44 (1.17; 1.78) for lnMMP-9; and 1.87 (1.40; 2.49) for lnTIMP-1, respectively (Table 2, model 1). LnMMP-10 was not significantly associated with incident CVD. After further adjustment for other cardiovascular risk factors, antihypertensive treatment and continuation of antihypertensive treatment at baseline, higher lnMMP-2 [1.49 (1.11; 1.99)] remained significantly associated with the incidence of cardiovascular events (Table 2, model 2). This HR remained significant after further adjustment for markers of LGI and ED (Table 2, models 4 and 5). However, further adjustment for eGFR and lnUAE (Table 2, model 3) decreased the HR from 1.49 (1.11; 1.99) to 1.34 (0.96; 1.85). This decline was explained mainly by the adjustment for eGFR, which decreased the HR to 1.36 (0.98; 1.88), whereas adjustment for lnUAE decreased the HR to 1.46 (1.08; 1.98).
Associations of MMPs and TIMP-1 with all-cause mortality
Additional file 1: Figure S3 shows the cumulative hazard plots of all-cause mortality according to tertiles of lnMMPs and lnTIMP-1. After adjustment for age, sex, duration of diabetes, HbA1c and nephropathy, higher levels of MMPs 1, 2, 3, 9 and 10, and TIMP-1 were significantly associated with all-cause mortality, with HRs of 1.62 (1.28; 2.06) per 1 SD higher lnMMP-1; 1.93 (1.48; 2.51) for lnMMP-2; 1.82 (1.38; 2.41) for lnMMP-3; 1.39 (1.14; 1.69) for lnMMP-9; 1.33 (1.06; 1.67) for lnMMP-10; and 2.10 (1.55; 2.85) for lnTIMP-1, respectively (Table 3, model 1). After adjustment for other cardiovascular risk factors, antihypertensive treatment and continuation of antihypertensive treatment at baseline, higher lnMMP-1 [1.38 (1.07; 1.78)], lnMMP-2 [1.60 (1.19; 2.15)] and lnMMP-3 [1.39 (1.05; 1.85)] remained significantly associated with all-cause mortality (Table 3, model 2). These associations were independent of eGFR, UAE, LGI and ED (Table 3, models 3–5), except for the association between lnMMP-3 and all-cause mortality, where the HR decreased from 1.39 (1.05; 1.85) to 1.32 (0.95; 1.83) after adjustment for eGFR and lnUAE (Table 3, model 3). This decrease was explained mainly by the adjustment for eGFR, which decreased the HR to 1.32 (0.96; 1.84), whereas adjustment for lnUAE decreased the HR to 1.36 (1.02; 1.82). The associations between MMP-1, MMP-2 and MMP-3 and all-cause mortality did not materially change after mutual adjustments for each of the other MMPs, except for the association between MMP-3 and all-cause mortality, which decreased after adjustment for MMP-1 [to 1.29 (0.96; 1.73)] and MMP-2 [to 1.20 (0.88; 1.62)].
Associations of MMPs and TIMP-1 with eGFR, lnUAE, LGI and ED
Cross-sectionally and after full adjustment for cardiovascular risk factors, lnMMP-1, -2, -3, -10 and lnTIMP-1 were significantly and inversely associated with eGFR at baseline (Table 4, model 2). LnMMP-2, -3, -10 and lnTIMP-1 were significantly and positively associated with lnUAE (Table 4, model 2). LGI was significantly associated with lnMMP-1, -9 and lnTIMP-1, and ED was significantly associated with lnMMP-1 and lnTIMP-1 (Table 4, model 2).
MMP-2 and -3 and eGFR decline during follow-up
The significant associations between lnMMP-2 and incident CVD and between lnMMP-3 and all-cause mortality were attenuated after adjustment for eGFR. To further investigate the role of eGFR in the association between plasma levels of MMPs and TIMP-1 on the one hand and incident CVD and all-cause mortality on the other, we investigated to what extent baseline MMPs and TIMP-1 levels were associated with changes in eGFR during follow-up. In addition, we investigated to what extent changes in eGFR (in addition to baseline eGFR) explained the associations between plasma markers and outcome measures when added to model 2. Follow-up data on eGFR, during follow-up measurements in 2003 , were available in 317/337 patients; 20 patients were excluded because they underwent kidney transplantation during follow-up and no data were available on eGFR pretransplantation (n = 14) or because no follow-up eGFR measurement was available (n = 6). Median follow-up duration was 9.9 years [9.1; 10.5]. The median decrease in eGFR was −1.67 ml/min/1.73 m2/year (interquartile range −3.18 to −0.64). After adjustment for age, sex, duration of diabetes, HbA1c, smoking, MAP, BMI, total cholesterol, antihypertensive medication and eGFR at baseline, lnMMP-2 and -3 (per 1 SD increase) were associated with decline in eGFR of −0.24 (95% CI −0.60; 0.12) and of −0.55 (−0.94; −0.16) ml/min/1.73 m2/year, respectively. Additional file 1: Figure S4 shows the decrease in eGFR per year according to tertiles of all MMPs and TIMP-1 after adjustment for the potential confounders used in model 2 plus baseline eGFR. The decline in eGFR was significantly greater in patients with nephropathy compared to those without (p < 0.001).
In the 317 patients with available data on eGFR during follow-up, the HR for developing CVD per 1 SD higher lnMMP-2 decreased from 1.60 (1.17; 2.18) to 1.40 (1.00; 1.96) after additional adjustment for baseline eGFR. Further adjustment for the decline in eGFR further decreased the HR to 1.36 (0.96; 1.92). This indicates that MMP-2-associated decrease of eGFR may partly explain the association between plasma MMP-2 and incident CVD. Next, in these 317 patients, the HR for all-cause mortality per 1 SD increase in lnMMP-3 decreased from 2.04 (1.45; 2.88) to 1.80 (1.22; 2.65) after additional adjustment for baseline eGFR. Further adjustment for the decrease in eGFR did not materially change the HR [1.79 (1.22; 2.62)].
Stratified analyses showed that, after full adjustment for cardiovascular risk factors (model 2), one SD higher lnTIMP-1 was associated with an increase in all-cause mortality in patients without [HR 3.65 (1.49; 8.92)], but not in those with nephropathy at baseline [HR 1.15 (0.75; 1.76)] (p for interaction was 0.009). In addition, after full adjustment for cardiovascular risk factors (model 2), one SD higher lnMMP-2 was associated with an increased HR for all-cause mortality in men [HR 1.99 (1.34; 2.95)], but not in women [0.85 (0.49; 1.46)] (p for interaction was 0.025).
Associations between plasma lnMMPs and lnTIMP-1 on the one hand and eGFR, lnUAE, LGI and ED on the other differed significantly between patients with and without nephropathy (Additional file 1: Tables S3, S4). The associations between plasma lnMMPs and lnTIMP-1 and eGFR and UAE were generally stronger in patients with diabetic nephropathy. Associations between markers of LGI and ED on the one hand and MMPs and TIMP-1 on the other did not differ consistently between patients with or without nephropathy.
The present study in type 1 diabetic patients, with a median follow-up of 12.3 years, shows significant associations between higher levels of plasma MMP-2 and incident CVD. In addition, higher plasma levels of MMP-1, -2 and -3 were significantly associated with all-cause mortality. These associations were independent of cardiovascular risk factors, LGI and ED (as expressed by z-scores), but partly dependent on eGFR and changes in eGFR. Despite the widespread correlations between the various MMPs and TIMP-1 (Additional file 1: Table S2), we showed specific associations as mentioned previously, which will be further discussed both with regards to current literature, but also on possible pathophysiological mechanisms.
MMPs, TIMP-1 and cardiovascular disease
Higher MMP-2-levels were independently associated with incident cardiovascular events, which is supported by the few studies available in mainly non-diabetic subjects [12, 28, 29]. In patients with atrial fibrillation, higher baseline plasma MMP-2 levels were independently associated with cardiovascular events or death during a mean follow-up of 28 months . In a Mexican population, a significant association between the specific MMP-2-1575 A/G polymorphism and occurrence of myocardial infarction has been established . Additionally, in patients with ST-elevation myocardial infarction undergoing percutaneous coronary intervention, higher plasma MMP-2-levels at baseline were associated with increased myocardial infarct size and decreased left ventricular function after four months of follow-up .
Various mechanisms may contribute to MMP-2-mediated macrovascular disease. We hypothesize that the occurrence of macrovascular disease can be enhanced by increased extracellular MMP-2 activity leading to plaque rupture, thrombus formation and vasoconstriction caused by MMP-2-mediated cleavage of vasoactive peptides. This hypothesis is supported by current literature [4, 30]. Additionally, altered MMP-2 activity may trigger cardiovascular pathology by vascular remodeling, i.e. increased migration of smooth muscle vascular cells to the intima, increased fibrosis and decrease of elastin content [3, 31]. Furthermore, MMP-2 leads to sustained thrombus formation by increased platelet activation, adhesion and aggregation .
In our study plasma MMP-1, -3, -9, -10 and TIMP-1-levels were not associated with incident CVD, which is in line with several prospective studies in the general population [32,33,34] and in patients with rheumatoid arthritis . In contrast, one study in 1127 patients with stable and unstable angina showed significant associations between higher MMP-9 levels and cardiovascular mortality with a HR of 2.41 (1.26; 4.60) in the upper quartile . The discrepancies with our study can possibly be explained by the differences in study population and by the fact that we used a combined endpoint for fatal and non-fatal CVD. For the combined endpoint of cardiovascular death and nonfatal myocardial infarction in the study by Blankenberg et al. , the unadjusted HR was 1.30 (1.11; 1.52, p < 0.002) with increasing MMP-9 quartiles, but it is not clear whether or not this remained significant after adjustment for potential confounders. In our study we did not show an association between TIMP-1-levels and occurrence of CVD. This is in contrast to a recent study which showed a significant association between plasma TIMP-1-levels and recurrent cardiovascular events in patients suffering from previous myocardial infarction or unstable angina 3–36 months before enrollment . This discrepancy might be explained by the difference in study populations, as in our study patients with prior CVD were excluded. In addition, in our study more adjustments for potential confounders were made compared with the latter two studies.
MMPs, TIMP-1 and all-cause mortality
Baseline MMP-1, -2 and -3-levels were significantly associated with all-cause mortality during follow-up. Stratification according to sex showed that the association between MMP-2 and all-cause mortality was restricted to men. However, the present study can neither explain this sex interaction nor can it exclude the possibility of a chance finding. Thus far, no prospective studies have been published regarding the association between plasma MMPs and all-cause mortality in type 1 diabetic patients. In different study populations higher circulating levels of MMP-1 , MMP-2 [12, 16, 17] and MMP-3 [18, 19] were strongly associated with all-cause mortality in patients with colorectal cancer , myocardial infarction [16, 18] or heart failure [17, 19]. MMP-1 is thought to be involved in cancer progression and invasion , whereas the mechanism by which MMP-3 contributed was not clarified [18, 19].
Plasma TIMP-1 levels did not show a significant association with all-cause mortality in the total study group and in patients with nephropathy at baseline. In contrast, plasma TIMP-1 was positively associated with all-cause mortality in patients without nephropathy at baseline. However, the observed interaction between TIMP-1 and nephropathy-status should be interpreted with caution as we cannot exclude the play of chance. Previous studies on the association between plasma TIMP-1 and mortality have been conflicting, both no  and positive associations [34, 37, 38] have been described. Positive associations between TIMP-1 and mortality were encountered in the general population , patients who underwent coronary angiography  and in patients with heart failure . These diverse results do not allow us to draw a solid conclusion on the association between TIMP-1 and all-cause mortality for all study populations. However, in type 1 diabetes, the present study is the first to report on this association.
Certain genetic polymorphisms coding for MMP-1 and MMP-3 have been associated with all-cause mortality . The MMP-1 2G/2G polymorphism has been associated with increased mortality in a 3-year follow-up study in 99 hemodialysis patients and 133 matched controls, HR 2.96 (1.29–6.80). In this study, the MMP-3 6A/6A polymorphism was non-significantly associated with all-cause mortality, HR 3.01 (0.88–10.3) and combined homozygosity of MMP-1 and MMP-3 polymorphisms lead to a HR of 4.69 (1.72–12.8). It is unknown, however, how these polymorphisms influence MMP expression and/or activity and the underlying mechanism of this association thus remains to be further investigated.
Potential mediation through eGFR
The association between plasma MMP-2 and CVD was attenuated after adjustment for baseline eGFR, which explained ~10% of the association. The HR decreased further after adjustment for decline in eGFR during follow-up. This may indicate that eGFR and decline in eGFR play a mediating role in the association between MMP-2 and CVD. In support, MMP-2 has been shown to induce complex alterations characterizing renal tubular epithelial-mesenchymal transformation leading to increased renal fibrosis , which could lead to decreased renal function. Higher plasma MMP-2 levels may also be the result of decreased renal function. However, MMP-2 levels are only minimally dependent on renal clearance due to the molecular weight of MMP-2, which is higher than that of albumin.
The association between MMP-3 and all-cause mortality was also attenuated by eGFR and this explained ~5% of the association. Although MMP-3 was significantly associated with eGFR decline during follow-up, the association between MMP-3 and mortality did not materially change after adjustment for eGFR decline. This makes a mediating role for eGFR in the association between MMP-3 and all-cause mortality less likely. Taken together, these results show that increased plasma levels of MMP-2 and -3, possible due to, but to a greater extent independent of, decreased eGFR are positively associated with CVD and all-cause mortality, respectively.
There are several limitations to our study. Firstly, plasma levels of MMP-1, -2, -3, -9, -10 and of TIMP-1 were only determined at baseline. Changes in plasma levels over time, revealing more details about outcome measures, are not known. Secondly, we do not know to what extent plasma levels reflect the local pathological situation at tissue level. Thirdly, we only determined five different MMPs and one of the four known TIMPs, considered to be the most important MMPs and TIMP to contribute to our outcome measures, but we cannot exclude that other MMPs and TIMPs are involved. Fourthly, the non-significant associations between the MMPs/TIMP-1 and outcome measures could be due to lack of statistical power, because of the relatively small number of events. Patients who died during follow-up were classified as cardiovascular unless another cause was established and this could have biased our results towards the null, resulting in underestimation of associations. Finally, the associations between MMPs and TIMP-1 with eGFR changes over time might be weakened by excluding those with renal transplantation during follow-up. Besides these limitations, the current study, with extensive adjustments for potential confounders and mediators, provides additional and consistent information on ECM remodeling by MMPs and TIMP-1 preceding incident CVD and all-cause mortality in a large cohort of patients with type 1 diabetes. These results are generalizable to normo- and macroalbuminuric type 1 diabetic patients. Although microalbuminuric patients were not included in the current study, we have no reason to expect different associations in this group.
In patients with type 1 diabetes followed for a median of 12 years, higher plasma MMP-2 levels are associated with incident CVD and higher plasma MMP-1, MMP-2, MMP-3 are associated with all-cause mortality. These associations are independent of cardiovascular risk factors, low-grade inflammation and endothelial dysfunction. However, baseline eGFR and decline in eGFR during follow-up attenuated the association between MMP-2 and CVD as well as the association between MMP-3 and all-cause mortality; eGFR may thus, in part, mediate these associations.
estimated glomerular filtration rate
tissue inhibitor of metalloproteinase
urinary albumin excretion
de Ferranti SD, de Boer IH, Fonseca V, Fox CS, Golden SH, Lavie CJ, Magge SN, Marx N, McGuire DK, Orchard TJ, et al. Type 1 diabetes mellitus and cardiovascular disease: a scientific statement from the American Heart Association and American Diabetes Association. Diabetes Care. 2014;37(10):2843–63.
Peeters SA, Engelen L, Buijs J, Chaturvedi N, Fuller JH, Schalkwijk CG, Stehouwer CD, Group EPCS. Plasma levels of matrix metalloproteinase-2, -3, -10, and tissue inhibitor of metalloproteinase-1 are associated with vascular complications in patients with type 1 diabetes: the EURODIAB Prospective Complications Study. Cardiovasc Diabetol. 2015;14:31.
Chen Q, Jin M, Yang F, Zhu J, Xiao Q, Zhang L. Matrix metalloproteinases: inflammatory regulators of cell behaviors in vascular formation and remodeling. Mediat Inflamm. 2013;2013:928315.
Gresele P, Falcinelli E, Loffredo F, Cimmino G, Corazzi T, Forte L, Guglielmini G, Momi S, Golino P. Platelets release matrix metalloproteinase-2 in the coronary circulation of patients with acute coronary syndromes: possible role in sustained platelet activation. Eur Heart J. 2011;32(3):316–25.
Visse R, Nagase H. Matrix metalloproteinases and tissue inhibitors of metalloproteinases: structure, function, and biochemistry. Circ Res. 2003;92(8):827–39.
Chase AJ, Newby AC. Regulation of matrix metalloproteinase (matrixin) genes in blood vessels: a multi-step recruitment model for pathological remodelling. J Vasc Res. 2003;40(4):329–43.
Giebel SJ, Menicucci G, McGuire PG, Das A. Matrix metalloproteinases in early diabetic retinopathy and their role in alteration of the blood-retinal barrier. Lab Invest. 2005;85(5):597–607.
Ebihara I, Nakamura T, Shimada N, Koide H. Increased plasma metalloproteinase-9 concentrations precede development of microalbuminuria in non-insulin-dependent diabetes mellitus. Am J Kidney Dis. 1998;32(4):544–50.
Thrailkill KM, Bunn RC, Moreau CS, Cockrell GE, Simpson PM, Coleman HN, Frindik JP, Kemp SF, Fowlkes JL. Matrix metalloproteinase-2 dysregulation in type 1 diabetes. Diabetes Care. 2007;30(9):2321–6.
Toni M, Hermida J, Goni MJ, Fernandez P, Parks WC, Toledo E, Montes R, Diez N. Matrix metalloproteinase-10 plays an active role in microvascular complications in type 1 diabetic patients. Diabetologia. 2013;56(12):2743–52.
Jacqueminet S, Ben Abdesselam O, Chapman MJ, Nicolay N, Foglietti MJ, Grimaldi A, Beaudeux JL. Elevated circulating levels of matrix metalloproteinase-9 in type 1 diabetic patients with and without retinopathy. Clin Chim Acta. 2006;367(1–2):103–7.
Ehrlich JR, Kaluzny M, Baumann S, Lehmann R, Hohnloser SH. Biomarkers of structural remodelling and endothelial dysfunction for prediction of cardiovascular events or death in patients with atrial fibrillation. Clin Res Cardiol. 2011;100(11):1029–36.
Blankenberg S, Rupprecht HJ, Poirier O, Bickel C, Smieja M, Hafner G, Meyer J, Cambien F, Tiret L. Plasma concentrations and genetic variation of matrix metalloproteinase 9 and prognosis of patients with cardiovascular disease. Circulation. 2003;107(12):1579–85.
West MJ, Nestel PJ, Kirby AC, Schnabel R, Sullivan D, Simes RJ, Pollicino C, Lubos E, Munzel TF, White HD, et al. The value of N-terminal fragment of brain natriuretic peptide and tissue inhibitor of metalloproteinase-1 levels as predictors of cardiovascular outcome in the LIPID study. Eur Heart J. 2008;29(7):923–31.
Tahara K, Mimori K, Iinuma H, Iwatsuki M, Yokobori T, Ishii H, Anai H, Kitano S, Mori M. Serum matrix-metalloproteinase-1 is a bona fide prognostic marker for colorectal cancer. Ann Surg Oncol. 2010;17(12):3362–9.
Dhillon OS, Khan SQ, Narayan HK, Ng KH, Mohammed N, Quinn PA, Squire IB, Davies JE, Ng LL. Matrix metalloproteinase-2 predicts mortality in patients with acute coronary syndrome. Clin Sci (Lond). 2010;118(4):249–57.
George J, Patal S, Wexler D, Roth A, Sheps D, Keren G. Circulating matrix metalloproteinase-2 but not matrix metalloproteinase-3, matrix metalloproteinase-9, or tissue inhibitor of metalloproteinase-1 predicts outcome in patients with congestive heart failure. Am Heart J. 2005;150(3):484–7.
Kelly D, Khan S, Cockerill G, Ng LL, Thompson M, Samani NJ, Squire IB. Circulating stromelysin-1 (MMP-3): a novel predictor of LV dysfunction, remodelling and all-cause mortality after acute myocardial infarction. Eur J Heart Fail. 2008;10(2):133–9.
Narayan H, Ng L, Squire I, Mohammed N, Quinn P. 51 Plasma MMP-3: a novel indicator of left ventricle remodelling and adverse outcomes in patients with acute heart failure. Heart. 2014;100(Suppl 3):A29–30.
Astrup AS, Tarnow L, Rossing P, Pietraszek L, Riis Hansen P, Parving HH. Improved prognosis in type 1 diabetic patients with nephropathy: a prospective follow-up study. Kidney Int. 2005;68(3):1250–7.
Tarnow L, Cambien F, Rossing P, Nielsen FS, Hansen BV, Lecerf L, Poirier O, Danilov S, Parving HH. Lack of relationship between an insertion/deletion polymorphism in the angiotensin I-converting enzyme gene and diabetic nephropathy and proliferative retinopathy in IDDM patients. Diabetes. 1995;44(5):489–94.
Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF 3rd, Feldman HI, Kusek JW, Eggers P, Van Lente F, Greene T, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–12.
Carrillo S. Updated information for IDMS traceable VITROS® chemistry products CREA slides. Rochester: Ortho clinical diagnostics (OCD); 2008.
De Jager J, Kooy A, Lehert P, Bets D, Wulffele MG, Teerlink T, Scheffer PG, Schalkwijk CG, Donker AJ, Stehouwer CD. Effects of short-term treatment with metformin on markers of endothelial function and inflammatory activity in type 2 diabetes mellitus: a randomized, placebo-controlled trial. J Intern Med. 2005;257(1):100–9.
Schram MT, Chaturvedi N, Schalkwijk CG, Fuller JH, Stehouwer CD, Group EPCS. Markers of inflammation are cross-sectionally associated with microvascular complications and cardiovascular disease in type 1 diabetes–the EURODIAB Prospective Complications Study. Diabetologia. 2005;48(2):370–8.
Nin JW, Jorsal A, Ferreira I, Schalkwijk CG, Prins MH, Parving HH, Tarnow L, Rossing P, Stehouwer CD. Higher plasma soluble receptor for advanced glycation end products (sRAGE) levels are associated with incident cardiovascular disease and all-cause mortality in type 1 diabetes: a 12-year follow-up study. Diabetes. 2010;59(8):2027–32.
van Bussel BC, Henry RM, Ferreira I, van Greevenbroek MM, van der Kallen CJ, Twisk JW, Feskens EJ, Schalkwijk CG, Stehouwer CD. A healthy diet is associated with less endothelial dysfunction and less low-grade inflammation over a 7-year period in adults at risk of cardiovascular disease. J Nutr. 2015;145(3):532–40.
Perez-Hernandez N, Vargas-Alarcon G, Martinez-Rodriguez N, Martinez-Rios MA, Pena-Duque MA, Pena-Diaz Ade L, Valente-Acosta B, Posadas-Romero C, Medina A, Rodriguez-Perez JM. The matrix metalloproteinase 2-1575 gene polymorphism is associated with the risk of developing myocardial infarction in Mexican patients. J Atheroscler Thromb. 2012;19(8):718–27.
Nilsson L, Hallen J, Atar D, Jonasson L, Swahn E. Early measurements of plasma matrix metalloproteinase-2 predict infarct size and ventricular dysfunction in ST-elevation myocardial infarction. Heart. 2012;98(1):31–6.
Berry E, Bosonea AM, Wang X, Fernandez-Patron C. Insights into the activity, differential expression, mutual regulation, and functions of matrix metalloproteinases and a disintegrin and metalloproteinases in hypertension and cardiac disease. J Vasc Res. 2013;50(1):52–68.
Lan TH, Huang XQ, Tan HM. Vascular fibrosis in atherosclerosis. Cardiovasc Pathol. 2013;22(5):401–7.
Jefferis BJ, Whincup P, Welsh P, Wannamethee G, Rumley A, Lennon L, Thomson A, Lawlor D, Carson C, Ebrahim S, et al. Prospective study of matrix metalloproteinase-9 and risk of myocardial infarction and stroke in older men and women. Atherosclerosis. 2010;208(2):557–63.
Welsh P, Whincup PH, Papacosta O, Wannamethee SG, Lennon L, Thomson A, Rumley A, Lowe GD. Serum matrix metalloproteinase-9 and coronary heart disease: a prospective study in middle-aged men. QJM. 2008;101(10):785–91.
Velagaleti RS, Gona P, Sundstrom J, Larson MG, Siwik D, Colucci WS, Benjamin EJ, Vasan RS. Relations of biomarkers of extracellular matrix remodeling to incident cardiovascular events and mortality. Arterioscler Thromb Vasc Biol. 2010;30(11):2283–8.
Vuilleumier N, Bas S, Pagano S, Montecucco F, Guerne PA, Finckh A, Lovis C, Mach F, Hochstrasser D, Roux-Lombard P, et al. Anti-apolipoprotein A-1 IgG predicts major cardiovascular events in patients with rheumatoid arthritis. Arthritis Rheum. 2010;62(9):2640–50.
Hansson J, Vasan RS, Arnlov J, Ingelsson E, Lind L, Larsson A, Michaelsson K, Sundstrom J. Biomarkers of extracellular matrix metabolism (MMP-9 and TIMP-1) and risk of stroke, myocardial infarction, and cause-specific mortality: cohort study. PLoS ONE. 2011;6(1):e16185.
Cavusoglu E, Ruwende C, Chopra V, Yanamadala S, Eng C, Clark LT, Pinsky DJ, Marmur JD. Tissue inhibitor of metalloproteinase-1 (TIMP-1) is an independent predictor of all-cause mortality, cardiac mortality, and myocardial infarction. Am Heart J. 2006;151(5):e1101–8.
Frantz S, Stork S, Michels K, Eigenthaler M, Ertl G, Bauersachs J, Angermann CE. Tissue inhibitor of metalloproteinases levels in patients with chronic heart failure: an independent predictor of mortality. Eur J Heart Fail. 2008;10(4):388–95.
Cozzolino M, Biondi ML, Galassi A, Turri O, Brancaccio D, Gallieni M. Matrix metalloproteinase-1 and matrix metalloproteinase-3 gene promoter polymorphisms are associated with mortality in haemodialysis patients. Nephrol Dial Transpl. 2009;24(7):2207–12.
Cheng S, Lovett DH. Gelatinase A (MMP-2) is necessary and sufficient for renal tubular cell epithelial-mesenchymal transformation. Am J Pathol. 2003;162(6):1937–49.
AJ collected the follow-up data. SP and LE researched data; SP, LE, CGS and CDS wrote the manuscript. JB, AJ, H-HP, LT and PR reviewed and edited the manuscript. CGS and CDS had equal contributions, contributed to the discussion and reviewed and edited the manuscript. All authors read and approved the final manuscript.
We would like to thank all the participants who took part in the study.
The authors declare that they have no competing interests.
Availability of data and materials
The dataset used during the current study is available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The study was approved by the local ethics committee, in accordance with the Helsinki Declaration, and all patients gave their written informed consent.
About this article
Cite this article
Peeters, S.A., Engelen, L., Buijs, J. et al. Plasma matrix metalloproteinases are associated with incident cardiovascular disease and all-cause mortality in patients with type 1 diabetes: a 12-year follow-up study. Cardiovasc Diabetol 16, 55 (2017). https://doi.org/10.1186/s12933-017-0539-1
- Type 1 diabetes
- Cardiovascular disease
- All-cause mortality
- Matrix metalloproteinase
- Tissue inhibitor of metalloproteinase
- Estimated glomerular filtration rate