The Relation between eNOS −786 C/T, 4 a/b, MMP-13 rs640198 G/T, Eotaxin 426 C/T, −384 A/G, and 67 G/A Polymorphisms and Long-Term Outcome in Patients with Coronary Artery Disease
Aim. The purpose of this study is to determine the association between eotaxin 426 C/T, −384 A/G, 67 G/A, eNOS −786 T/C, 4 a/b, and MMP-13 rs640198 G/T and prognosis of patients with known CAD. Methods. From total of 1161 patients referred to coronary angiography, 532 patients with angiographically confirmed CAD were selected. Their long-term outcome was followed up using hospital database. Subsequent events were assessed in this study: death or combined endpoint-myocardial infarction, unstable angina pectoris, revascularization, heart failure hospitalization, and cardioverter-defibrillator implantation. Results. The multivariate Cox regression model identified age, smoking, and 3-vessel disease as significant predictors of all-cause death. Further analysis showed that eotaxin 67 G/A (GA + AA versus GG) and eotaxin −384 A/G (GG versus GA + AA) were significant independent prognostic factors when added into the model: HR (95% CI) 2.81 (1.35–5.85), ; HR (95% CI) 2.63 (1.19–5.83), ; eotaxin −384 A/G was significantly associated with the event-free survival, but it did not provide the prognostic information above the effect of two- or three-vessel disease. Conclusion. The A allele in eotaxin 67 G/A polymorphism is associated with worse survival in CAD patients.
Cardiovascular diseases, particularly coronary artery disease (CAD), are the leading cause of mortality in developed countries. The aetiology of CAD is multifactorial and various well-known risk factors have been described in previous decades. Association of CAD prevalence and deoxyribonucleic acid (DNA) polymorphisms in human genes has been extensively studied [1, 2], but relation of these polymorphisms to long-term prognosis of patients was reported less frequently [3–6]. Some of the polymorphisms did not show any significant association with the prognosis of CAD or the results were controversial. On the other hand, several papers reported promising results; for example, genome-wide association study (GWAS) identified polymorphisms in loci 1p13.3 (rs599839), 1q41 (rs17465637), and 3q22.3 (rs9818870) , or the variants of paraoxonase 1 gene (rs854560)  and phospholipase A2 gene (rs1805017) . These data were repetitively confirmed by following studies [8–14]. However, there are many variants in human genome influencing the presence of CAD, of which their relationship to long-term outcome in CAD patients is still not elucidated as the information based on prospective studies is lacking.
Matrix metalloproteinases are zinc-dependent endopeptidases cleaving proteins of extracellular matrix . Matrix-metalloproteinase 13 (MMP-13) preferentially cleaves type II collagen and is produced by macrophages and smooth muscle cells in atherosclerotic lesions . This contributes to restructuring of atherosclerotic plaque in animal models . Polymorphisms of MMP-13 gene influencing its production have been associated with an extent of atherosclerosis , including rs640198 G/T intronic single nucleotide polymorphism (SNP) .
Eotaxin (CCL11) is a chemokine with pleiotropic proinflammatory effects produced by several cellular types present in atherosclerotic lesion . Studies in human subjects have found eotaxin overexpression in atherosclerotic plaques  or its elevated plasmatic concentration in the presence of coronary atherosclerosis [21, 22]. Genetic variation of eotaxin gene, SCYA11, has been found to contribute to CAD susceptibility, most notably the first exon (67 G/A) polymorphism leading into amino acid change Ala/Thr in site 23, where secreted eotaxin is cleaved off from the signalling peptide . The association of eotaxin (CCL 11) Ala/Ala variant and development of coronary artery disease was reported by Wang et al. . The authors performed a prospective study of 1297 Chinese patients with diabetes mellitus; median follow-up was 8 years. The variant Ala/Ala of eotaxin Ala23Thr polymorphism showed a hazard ratio of 1.7 (95% CI 1.10–2.61, ). Other functionally relevant polymorphisms were identified in the promoter of the gene, as −384 A/G and −426 C/T , which has been associated with eotaxin levels and asthma phenotype [25, 26].
Endothelial nitric oxide synthase (eNOS) is an enzyme that constitutively produces nitric oxide from L-arginine, mediating vasodilation and having other beneficial effects on vessel wall . The rate of NO synthesis by eNOS is modulated by variability in eNOS gene, which is the case of 27-bp tandem repetition in intron 4, where an “A allele” (consisting of four repetitions) is associated with higher NO production than a “B allele” (consisting of five repetitions) . Other polymorphisms, −786 T/C promoter SNP, have been linked to changes in eNOS promoter activity and coronary spasms .
The aim of our study is to assess the relationship between polymorphisms in eotaxin gene (67 G/A, −384 A/G, and −426 C/T), MMP- 13 gene (rs640198), and endothelial NO-synthase (eNOS) gene (−786 C/T and 4 a/b) and a long-term outcome in Caucasian patients initially referred to coronary angiography.
2. Patients and Methods
2.1. Patient Cohort and Genetic Analyses
The study comprised of 532 subjects suffering of CAD that were recruited from the total number of 1161 consecutive patients referred to coronary angiography after applying the following exclusion criteria: known malignancy, advanced renal dysfunction (blood creatinine level above 200 μmol/L), smooth coronary arteries or insignificant atherosclerosis, or life expectancy below 1 year. All patients gave their informed consent and study was approved by the institutional ethic committee. The personal and family history were obtained and all patients underwent physical examination, electrocardiography, haematological and biochemical blood tests, and coronary angiography. Coronary artery disease was defined as ≥50% stenosis of at least one coronary artery. Blood samples for DNA analysis from peripheral blood leucocytes were also taken. The selection of DNA polymorphism, analyses, results, and impact on incidence and severity of CAD were described in our previous publications [19, 30, 31].
The clinical outcome of patients was assessed in 2014 using hospital database. The following events were recorded: death, acute myocardial infarction (MI), unstable angina pectoris, necessity of myocardial revascularization (either percutaneous or coronary artery bypass graft), hospitalization for heart failure, and cardioverter/defibrillator (ICD) implantation. Patients who were not under clinical observation in our hospital, patients without angiographically confirmed significant CAD, or those with unsuccessful DNA analysis were excluded from the study.
2.3. Statistical Analysis
Multivariate Cox regression model was used for estimating the contribution of genetic polymorphisms and other risk factors to survival. Genetic variants had been preselected out of the total number of 8 polymorphisms in candidate genes using Kaplan-Meier method and log-rank tests in dominant, recessive, and codominant mode of allelic expression. value 0.05 was used as a cut-off for an inclusion of variable in the Cox regression analysis. Hardy-Weinberg equilibrium was calculated for each polymorphism using test. Nongenetic variables included in the model were the following: age on admission, sex, presence of acute coronary syndrome, MI in personal history, ejection fraction, number of diseased vessels, obesity, hypertension, hyperlipidaemia, diabetes mellitus, concomitant presence of cardiomyopathy, significant valvular disease, and smoking status. Two endpoints were observed: all-cause death and cardiac endpoint consisting of MI, unstable angina, coronary revascularization, and rehospitalization for heart failure. Supposing unknown dominance of alleles, three modes of gene expressions were used: dominant, where minor allele carriers were compared with major allele homozygotes; recessive, where major allele carriers were compared with minor allele homozygotes; and additive, gene dose-based model, where the values 0, 1, and 2 were attributed to major allele homozygote, heterozygote, and minor allele homozygote. Stepwise construction of a Cox regression model with value to include 0.05 and value to remove 0.051 was finally employed to determine all significant independent factors (clinical and genetic) contributing to patients’ survival or the incidence of cardiac event. STATISTICA software (StatSoft, version 12) was used for statistical analysis. MIDAS software version 1.0  was used for Hardy-Weinberg and linkage disequilibrium determination. Generally, values ≤ 0.05 were considered to be statistically significant.
Five hundred thirty-two patients with significant CAD were finally comprised into analysis. The basic characteristics are listed in Table 1. Twenty-nine patients died and 271 cardiac events were recorded (26 patients with unstable angina, 44 acute MI, 154 myocardial revascularization, 26 hospitalization for heart failure, and 21 ICD implantation) during follow-up, median follow-up period was 77 months. Kaplan-Meier curves presenting proportion of survival and event-free survival can be seen in Figures 1 and 2. The list of studied polymorphisms and their genotype distribution are shown in Table 2. The preselection of candidate polymorphisms was based on log-rank values. The eNOS 4 a/b, eotaxin −426 C/T, eotaxin −384 A/G, and eotaxin 67 G/A were found to be significantly associated with death, as well as eotaxin −384 A/G with cardiac events. There was also nonsignificant trend for eotaxin 67 G/A to influence the event-free survival.
Complete data are shown in Tables 3 and 4. Subsequently, the multivariate Cox regression model was used to reveal potential effects of selected variants on survival and cardiac endpoint-free survival. The eotaxin 67 G/A and eotaxin −384 A/G were considered significant for survival and event-free survival, respectively (Tables 5 and 6). The multivariate stepwise Cox regression model identified age, smoking, and three-vessel disease (compared to one vessel disease) as significant predictors of all-cause death. Together with those factors, eotaxin 67 G/A was a significant independent prognostic factor when added into the model (Table 7, Figure 3). When eotaxin −384 A/G was added into the stepwise Cox regression model instead of 67 G/A, the model gave similar results with the same prognostic factors; however, hazard ratio for the given polymorphism was slightly lower (GG versus GA + AA; HR = 2.63; 95% CI = 1.19–5.83; ). As shown below, both polymorphisms were in significant linkage disequilibrium, and the “risky” alleles of both SNPs were frequently inherited together. In the multivariate stepwise Cox regression analysis of cardiac endpoints, a number of diseased coronary arteries (two versus one: HR = 1.86; 95% CI = 1.23–2.79; , three versus one: HR = 4.80; 95% CI = 3.32–6.94; ) was the only significant independent risk factor. Adding eotaxin −384 A/G into the model did not change the results; that is, although the −384 A/G polymorphism was associated with the hazard of cardiac endpoint both with and without adjustments, it does not provide the prognostic information above the effect of two- or three-vessel disease. Finally, it was determined, whether the variants in eotaxin and eNOS gene were or were not in linkage equilibrium. The analysis revealed modest, but significant linkage disequilibrium (LD) between polymorphisms −426 C/T and −384 A/G in eotaxin gene promoter (; ; ), as well as between −384 A/G and 67 G/A in the same gene (; ; ). In this polymorphism, the risky A allele of 67 G/A was more often inherited together with the risky G allele of 384 A/G, suggesting the possible role of LD in resulting effects. On the other hand, −426 C/T and 67 G/A were in linkage equilibrium with each other. In eNOS gene, the eNOS −786 C/T and eNOS 4 a/b polymorphisms were also in significant LD with each other (; ; ). The B allele of eNOS 4 a/b polymorphism tended to be inherited with the T allele of −786 C/T, while the A allele of the former polymorphism was associated with the C allele of the latter in most cases.
Many studies have been published about the influence of plasma eotaxin levels and variants of eotaxin gene on presence and severity of atherosclerosis but the results are controversial. Several authors reported an association of higher plasma levels with atherosclerotic burden [21, 22, 33]; on the other hand, other studies did not prove this [34, 35]. Higher eotaxin concentration was found in local blood samples from patients with acute coronary syndrome . Among other effects, eotaxin was proved to stimulate migration of smooth muscle cells from media to intima of injured arterial wall . This mechanism is considered to have an impact on atherosclerotic plaque progress and restenosis formation [38, 39]. Influence of eotaxin gene polymorphisms on cardiac events in Chinese diabetic patients was reported in the abovementioned paper ; the main result was increased morbidity in patients homozygote for the G allele of 67 G/A SNP in the first exon of the gene corresponding to alanine in resulting protein; this allele is associated with higher production and concentration of eotaxin [40, 41]. Likewise, in our previous study , there was greater susceptibility to acute coronary syndromes in subjects with homozygous 67 G/A GG genotype. This seemingly contradicts our current results, where GG genotype showed protective effect. However, given all subjects in our study suffered from CAD, our study group consisted of highly selected patients and their treatment may have been responsible for the results. Actually, the influence of A allele on prognosis can be explained by an effect of statin therapy. There is some evidence from both human patients and animal models of substantial chemokine levels reduction after statin treatment, including eotaxin [42–44]. Because A allele carriers have lower eotaxin levels , the protective effect of statin therapy may be reduced in comparison with GG homozygotes. This mechanism seems to be likely in our population of CAD patients, where the prevalence of statin therapy was very high (514/532, 96.6%).
In a large study with 14,916 American men, Zee et al.  studied the same polymorphism with different results. Authors compared data from 523 patients with myocardial infarction and 2092 persons who were free from cardiovascular disease during a follow-up (mean period 13.2 years). They reported that T23 (67 A) homozygotes had significantly higher risk of myocardial infarction: odds ratio (OR) 1.86; 95% confidence interval (CI), 1.15–3.01; ; the effect remains significant even after adjusting for body mass index, history of hypertension, the presence of diabetes, and randomized treatment assignment. There was no significant difference between GG homozygotes and heterozygotes. Since there was no information about a treatment in that study, its impact cannot be excluded. The recessive model (AA versus GA + GG) in our study had relatively low power in all cases, there were only 16 AA homozygotes, to detect eventual risky effect of AA homozygote genotype. The linkage disequilibrium between 67 G/A and −384 A/G must be taken into account as well. Unlike the former one, the latter SNP was associated also with cardiac events and the “protective” and “risky” alleles were frequently inherited together.
The influence of genetic variants of different matrix metalloproteinase genes, such as MMP-9 [47, 48] or MMP-3 , on clinical outcome of CAD was described in previous studies. In our current study, no effect of rs640198 polymorphism in MMP-13 gene on either cardiac events or all-cause mortality was found both with and without adjustments. As no other study associating MMP-13 genetic variants and CAD exists so far, we can thus conclude that MMP-13 does not predict the prognosis of CAD, even when in our previous study based on the same population we found significant association of rs640198 with triple vessel disease. However, triple vessel disease itself was a significant predictor of both cardiac events and all-cause mortality. The gene-prognosis relation may be modified by commonly used treatment, as in the case of eotaxin.
The association of eNOS polymorphisms with CAD was reported in many studies. The large meta-analysis involving a total sample size of 69,235 subjects  confirmed the association of the three NOS3 gene polymorphisms (Glu298Asp, −786 T/C, and 4 a/b) with the presence of CAD, but no significance was found for 4 a/b polymorphism in European population. No impact of −786 C/T or 4 a/b on CAD prognosis was found in our study.
In conclusion, a significant association between A allele in eotaxin 67 G/A single nucleotide polymorphism (SNP) in the first exon of the gene and lower survival rate in patients with coronary artery disease was found with and without adjustments for other factors. Similar results were observed for −384 A/G polymorphism in the same gene, which was also associated with cardiac events. Further studies will be needed to elucidate the mechanism of this effect.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This study was supported by European Regional Development Fund, Project FNUSA-ICRC (no. CZ.1.05/1.1.00/02.0123).
S. Blankenberg, H. J. Rupprecht, C. Bickel et al., “Common genetic variation of the cholesteryl ester transfer protein gene strongly predicts future cardiovascular death in patients with coronary artery disease,” Journal of the American College of Cardiology, vol. 41, no. 11, pp. 1983–1989, 2003.View at: Publisher Site | Google Scholar
K. L. Ellis, C. M. Frampton, A. P. Pilbrow et al., “Genomic risk variants at 1p13.3, 1q41, and 3q22.3 are associated with subsequent cardiovascular outcomes in healthy controls and in established coronary artery disease,” Circulation: Cardiovascular Genetics, vol. 4, no. 6, pp. 636–646, 2011.View at: Publisher Site | Google Scholar
A. Brautbar, C. M. Ballantyne, K. Lawson et al., “Impact of adding a single allele in the 9p21 locus to traditional risk factors on reclassification of coronary heart disease risk and implications for lipid-modifying therapy in the atherosclerosis risk in communities study,” Circulation: Cardiovascular Genetics, vol. 2, no. 3, pp. 279–285, 2009.View at: Publisher Site | Google Scholar
K. J. Haley, C. M. Lilly, J.-H. Yang et al., “Overexpression of eotaxin and the CCR3 receptor in human atherosclerosis: using genomic technology to identify a potential novel pathway of vascular inflammation,” Circulation, vol. 102, no. 18, pp. 2185–2189, 2000.View at: Publisher Site | Google Scholar
Y. Wang, A. O. Y. Luk, R. C. W. Ma et al., “Independent predictive roles of eotaxin Ala23Thr, paraoxonase 2 Ser311Cys and β3-adrenergic receptor Trp64Arg polymorphisms on cardiac disease in type 2 diabetes—an 8-year prospective cohort analysis of 1297 patients,” Diabetic Medicine, vol. 27, no. 4, pp. 376–383, 2010.View at: Publisher Site | Google Scholar
L. Rigoli, L. Caminiti, C. Di Bella et al., “Investigation of the eotaxin gene −426C → T, −384A → G and 67G → A single-nucleotide polymorphisms and atopic dermatitis in Italian children using family-based association methods,” Clinical and Experimental Dermatology, vol. 33, no. 3, pp. 316–321, 2008.View at: Publisher Site | Google Scholar
J. Máchal, A. Vašků, V. Kincl et al., “Association between three single nucleotide polymorphisms in eotaxin (CCL 11) gene, hexanucleotide repetition upstream, severity and course of coronary atherosclerosis,” Journal of Applied Genetics, vol. 53, no. 3, pp. 271–278, 2012.View at: Publisher Site | Google Scholar
V. Kincl, A. Vasků, J. Meluzín, R. Panovský, J. Seménka, and L. Groch, “Association of the eNOS 4a/b and -786T/C polymormphisms with coronary artery disease, obesity and diabetes mellitus,” Folia Biologica, vol. 55, pp. 187–191, 2009.View at: Google Scholar
B. V. Loughrey, A. McGinty, I. S. Young, D. R. McCance, and L. A. Powell, “Increased circulating CC chemokine levels in the metabolic syndrome are reduced by low-dose atorvastatin treatment: evidence from a randomized controlled trial,” Clinical Endocrinology, vol. 79, no. 6, pp. 800–806, 2013.View at: Publisher Site | Google Scholar
Y. Sheikine, B. Olsen, B. Gharizadeh, K. Jatta, P. Tornvall, and M. Ghaderi, “Influence of eotaxin 67G>A polymorphism on plasma eotaxin concentrations in myocardial infarction survivors and healthy controls,” Atherosclerosis, vol. 189, no. 2, pp. 458–463, 2006.View at: Publisher Site | Google Scholar
T. B. Opstad, H. Arnesen, A. Å. Pettersen, and I. Seljeflot, “The MMP-9 -1562 C/T polymorphism in the presence of metabolic syndrome increases the risk of clinical events in patients with coronary artery disease,” PLoS ONE, vol. 9, no. 9, Article ID e106816, 2014.View at: Publisher Site | Google Scholar