The Chromosome 9p21 Variant Not Predicting Long-Term Cardiovascular Mortality in Chinese with Established Coronary Artery Disease: An Eleven-Year Follow-Up Study
Introduction. We examined whether the variant at chromosome 9p21, rs4977574, was associated with long-term cardiovascular mortality in Han Chinese patients with coronary artery disease (CAD). Methodology. Subjects who underwent coronary angiography for chest pain were consecutively enrolled. Fasting blood samples were collected for laboratory and genotype assessments. The information was correlated with data collected from the national death database. Results. There were 925 cases with CAD and 634 without CAD enrolled in the present study. The G allele conferred a significant increase in risk of CAD (odds ratio = 1.47, in the dominant model; odds ratio = 1.36, in the recessive model). During a median of 11 years (inter-quartile range between 5.2 and 12.5 years) of follow-up, neither the total nor the cardiovascular mortality was different among CAD subjects with different genotypes. Using Cox regression analysis, genotypes of rs4977574 still failed to predict cardiovascular mortality (hazard ratio = 1.25, in the dominant model; hazard ratio = 1.05, in the recessive model). Conclusions. The rs4977574 at chromosome 9p21 is associated with presence of CAD in Han Chinese. However, rs4977574 could not predict cardiovascular mortality in these CAD subjects during the eleven-year period of the study.
Coronary artery disease (CAD) is the leading cause of death worldwide . Heredity plays an important role in CAD susceptibility . The chromosome 9p21 region has been shown to be associated with CAD, and several single nucleotide polymorphisms (SNPs) selected from Caucasians have been replicated in Chinese [3–5]. However, SNPs in chromosome 9p21 contribute diverse intensities of risk to cardiovascular disease, and it has been reported that the effect size of identified SNPs is dependent on the chosen population [6–8].
Among the tag SNPs, rs4977574 has been initially demonstrated to be associated with risk of CAD in the Myocardial Infarction Genetics Consortium (MIGen) study . The association between rs4977574 and CAD has also been found in several populations [10–13]. However, a small study comprising 334 subjects did not show any significant association between CAD risk and rs4977574 in a Chinese population .
Although rs4977574 was recently reported to be associated with cardiovascular incidence in the subjects without CAD at baseline , results from longitudinal assessments between 9p21 variants and cardiovascular mortality in subjects with established CAD were inconsistent. In the Global Registry of Acute Coronary Events (GRACE) genetic study, the 9p21 variant showed a significant association with cardiovascular death . However, Virani et al.  recently reported that the 9p21 variant was not associated with cardiovascular outcome. In order to clarify the long-term impact of 9p21 variant on cardiovascular mortality in CAD patients, we therefore examined the polymorphism of rs4977574 in a Taiwanese cohort of Han Chinese descent in an eleven-year follow-up duration.
2. Material and Methods
The participants were enrolled during hospitalization in the Cardiovascular Department in Taichung Veterans General Hospital between March 1999 and March 2000. The study comprised both case-control and follow-up assessments for subjects undergoing coronary angiography due to angina and clinical suspicion of ischemic heart disease. Blood samples were collected after the subjects had fasted overnight. Mortality data up to 2011 were collected for the follow-up assessment. The study protocol was approved by the Institutional Review Board of Taichung Veterans General Hospital, Taichung, Taiwan.
Genomic DNA was extracted from patients’ peripheral blood leukocytes using a QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany). Genotypes of rs4977574 were determined using TaqMan allelic discrimination assay (Custom TaqMan SNP Genotyping Assay). The genotyping reaction was amplified and detected on a StepOnePlus Real-Time PCR (Applied Biosystems, Foster City, USA).
Plasma glucose was determined using the glucose oxidase-peroxidase method (Wako Diagnostics, Tokyo, Japan). Serum cholesterol and triglyceride levels were determined by the enzymatic method using commercial kits (Wako Diagnostics, Tokyo, Japan). High-density lipoprotein (HDL) cholesterol was determined after the precipitation of apo B-containing lipoproteins. Low-density lipoprotein (LDL) cholesterol concentration was calculated using the formula of Friedewald et al.  if fasting serum triglyceride was less than 4.5 mmol/L. Otherwise, LDL cholesterol was determined after separation of very-low-density lipoprotein (VLDL) from serum by ultracentrifugation and precipitation of apo B-containing particles.
2.3. Exposure Information
The subjects were assigned to the CAD group if they fulfilled one or both of the following criteria: (1) a history of myocardial infarction or coronary revascularization or (2) 50% or more luminal narrowing in any coronary artery according to angiography. Subjects who did not fulfill any of these criteria were assigned to the non-CAD group. Diabetes mellitus was defined as (1) a history of diabetes mellitus or use of antidiabetic medications or (2) fasting blood glucose levels equal to or higher than 7.0 mmol/L. Hypertension was defined as (1) a history of hypertension or use of antihypertensive medications or (2) systolic blood pressure higher than 140 mmHg or diastolic blood pressure higher than 90 mmHg.
The coronary artery lesions were assessed at the angiography viewing workstation using software for quantitative analysis (Philips Inturis Suite, R2.2, Philips Medical Systems, Eindhoven, Netherlands). The baseline CAD severity was evaluated by Gensini score, which was calculated by the sum of obstructive severity score times the lesion location score system .
The mortality data were provided by the Collaboration Center of Health Information Application (CCHIA), Department of Health, Executive Yuan, Taiwan. The causes of death were determined according to the International Classification of Disease, 9th Revision, Clinical Modification (ICD-9 CM) diagnostic criteria. Cardiovascular disease was defined as CAD (402, 404, 410, 411, 412, 413, 414, 426, 427, 428, and 429), cerebrovascular disease (430, 431, 432, 433, 434, 435, 436, 437, and 785.9), and peripheral artery disease (440.2, 443.8, 443.9, and 444.2) .
2.4. Statistical Analysis
Data are presented as mean ± standard deviation (SD). The genotype distribution was tested for Hardy-Weinberg equilibrium by the goodness-of-fit test. Chi-square test was used to assess the differences in gender, genotypes, diabetes, hypertension, and smoking status. The differences in continuous variables at baseline were analyzed by independent -test. Multivariate logistic regression analyses were used to assess the genetic effects on CAD at baseline. Multivariable Cox proportional-hazard regression analyses were used to examine the cardiovascular mortality between genotypes. Statistical analysis was performed using SPSS 19.0 (SPSS Inc., Chicago, IL, USA).
3.1. Case-Control Assessments at Baseline
A total of 1559 subjects were enrolled with 925 in the CAD group and 634 in the non-CAD group (Figure 1). The distribution of alleles was in agreement with the Hardy-Weinberg equilibrium. The clinical features of the two study groups are shown in Table 1. The CAD group was significantly older with male dominance compared with the non-CAD group (both values less than 0.001). There were significantly higher proportions of current smoking status (), hypertension (), and diabetes mellitus () in the CAD group than those in the non-CAD group. Fasting serum concentrations of LDL cholesterol and total cholesterol were also significantly higher in the CAD group than those in the non-CAD group (both values less than 0.001). The genotype distributions showed that G allele of rs4977574 was overrepresented in the CAD group compared to the non-CAD group ( for genotype distribution and for allele distribution, resp.). The risk allele remained independently associated with CAD after adjusting for age, gender, smoking, hypertension, diabetes, and LDL cholesterol in the logistic regression model (Table 2).
3.2. Longitudinal Cardiovascular and Total Mortality for CAD Patients
Among the 925 patients with CAD at baseline, there were 248 with GG genotype, 479 with GA genotype, and 198 with AA genotype (Table 3). The CAD severity, based on the Gensini score, was not significantly different among all genotypes (). During the follow-up period (median of 11.0 years, interquartile range between 5.2 and 12.5 years), there was no significant difference in cardiovascular mortality or total mortality among the CAD patients with different genotypes (Figure 2). Using Cox regression analysis, the risk G allele was not associated with cardiovascular mortality and total mortality in either the dominant model or the recessive model after adjustment for age, gender, smoking status, past history of myocardial infarction, statins treatment, presence of hypertension, and presence of diabetes (Table 4).
Our major findings were that the G allele of SNP rs4977574 was associated with the presence of CAD but failed to predict cardiovascular death in the Han Taiwanese patients with established CAD. The case-control findings of our study are in-line with results reported in other populations [9–13]. SNP rs4977574 is located at a locus with a gene encoding a large antisense noncoding RNA (ANRIL) [21, 22]. ANRIL might be involved in atherosclerosis by being expressed in the endothelial cells and monocyte-derived cells [21, 23]. Upstream from the locus, the cyclin-dependent kinase inhibitor 2A (CDKN2A) and 2B (CDKN2B) genes are associated with transcriptional regulation of ANRIL [24–26]. It is also important to note that the baseline Gensini scores were not significantly different among the CAD patients with three different genotypes. In accordance with our observations, Chen et al.  reported that the 9p21 variant is not associated with CAD severity in subjects with established CAD. Similar findings, reported in Chinese subjects with myocardial infarction, showed that severity of coronary lesion was not significantly different among different genotypes of the 9p21 variant . Further investigations are warranted to clarify the underlying pathogenesis contributing to the observed association between chromosome 9p21 variants and coronary lesions.
In the present study, rs4977574 failed to predict cardiovascular or total mortality in subjects with CAD during the eleven-year period of follow-up. Our finding in CAD Chinese was different from results in subjects without CAD . Furthermore, the GRACE genetic study showed the 9p21 variant was significantly associated with recurrent myocardial infarction within only six months , but Peng et al.  reported that the 9p21 variant could not predict total mortality in subjects following myocardial infarction in a 2.5-year study. Virani et al.  also found no significant association between the 9p21 variant and cardiac death in subjects with established CAD in a 3.2-year follow-up. One of the main strengths of the present study was that long-term mortality in a follow-up period of 11 years was analyzed in CAD patients, using death records as the endpoint to test our hypothesis.
It is well acknowledged that the current established risk factors of atherosclerosis are not able to fully account for cardiovascular death . The discrepancy probably results from different pathophysiological mechanisms. Endothelial dysfunction and monocyte activation play roles in inflammation and subsequent development of coronary lesions [30, 31]. However, the atherosclerotic biomarkers and monocyte activity were poor predictors of longitudinal cardiovascular events in patients with established CAD . In fact, cardiovascular death might be associated with the critical mechanisms of plaque stabilization and arrhythmia in CAD patients [33–35]. In large-scale genome-wide association studies (GWAS), susceptibility variants were determined based on case-control studies, and thus it might not be possible to prospectively predict the outcomes of diseases . Previous authors have asserted that the effects of the 9p21 genetic polymorphisms might be only associated with coronary lesion but not capable of predicting cardiovascular mortality [37, 38]. Furthermore, cardiovascular death might be associated with many common diseases, which in turn may be linked to both genetic and environmental factors [39, 40]. The association between chromosome 9p21 variants and sequential cardiovascular mortality might be abated by the medical treatment or procedural intervention during the follow-up period . Along the same lines, it has been reported that the 9p21 variants could not predict the cardiovascular outcomes in subjects with coronary lesions after drug-eluting stent placement .
Cavender et al.  reported a total mortality of 37% in subjects with CAD at five years. In the present study, the accumulative total mortality of about 55% was observed in a median 11-year period of follow-up. Ethnic differences in mortality were reported in CAD subjects, with a higher mortality in Asians compared with that in Caucasians [44, 45]. However, the prognosis of CAD seems to be better in Chinese than that in south Asians, although the mechanism is not well understood . In a recent report from a Vienna cohort study with a similar duration of follow-up , a lower total mortality rate (39%) in the CAD subjects was observed, although the proportion of cardiovascular cause (57%) was similar to that in the present study (56.6%). Because the majority of deaths in the present study were cardiovascular causes, the total mortality findings were close to those of cardiovascular mortality.
There were some limitations in the current study. First, the cardiovascular mortality data was collected from the national mortality registry which covers 99.6% of the total population in Taiwan . However, nonfatal cardiovascular events could not be traced in the follow-up assessment as these data are not provided in the database. Second, we assessed cardiovascular mortality rather than CAD mortality due to indistinct classifications of death causes; for example, some types of CAD might be classified to cardiovascular disease (ex. ICD-9 code of 429.2) . Furthermore, we could not collect information on the effect of environmental factors and treatments during the follow-up period.
In conclusion, the SNP rs4977574 at chromosome 9p21 was an independent risk factor for CAD in Taiwanese individuals of Han Chinese descent. However, this SNP could not predict the total or cardiovascular mortality in these subjects, who had CAD at baseline, during the eleven-year period of follow-up.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was supported by Grants from Taichung Veterans General Hospital, Taichung, Taiwan, [TCVGH-1013501B and TCVGH-103DHA0500405] and National Science Council, Taiwan [NSC 102-2314-B-075A-001].
R. Lozano, M. Naghavi, K. Foreman et al., “Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010,” The Lancet, vol. 380, no. 9859, pp. 2095–2128, 2013.View at: Google Scholar
J. Guo, W. Li, Z. Wu, X. Cheng, Y. Wang, and T. Chen, “Association between 9p21. 3 genomic markers and coronary artery disease in East Asians: a meta-analysis involving 9, 813 cases and 10, 710 controls,” Molecular Biology Reports, vol. 40, no. 1, pp. 337–343, 2013.View at: Publisher Site | Google Scholar
L. Qi, J. M. Li, H. Sun et al., “Association between gene polymorphisms and myocardial infarction in Han Chinese of Yunnan province,” Chinese Journal of Medical Genetics, vol. 29, no. 4, pp. 413–419, 2012.View at: Google Scholar
I. Buysschaert, K. F. Carruthers, D. R. Dunbar et al., “A variant at chromosome 9p21 is associated with recurrent myocardial infarction and cardiac death after acute coronary syndrome: the GRACE genetics study,” European Heart Journal, vol. 31, no. 9, pp. 1132–1141, 2010.View at: Publisher Site | Google Scholar
S. S. Virani, A. Brautbar, V.-V. Lee et al., “Chromosome 9p21 single nucleotide polymorphisms are not associated with recurrent myocardial infarction in patients with established coronary artery disease,” Circulation Journal, vol. 76, no. 4, pp. 950–956, 2012.View at: Publisher Site | Google Scholar
W. T. Friedewald, R. I. Levy, and D. S. Fredrickson, “Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge,” Clinical Chemistry, vol. 18, no. 6, pp. 499–502, 1972.View at: Google Scholar
G. G. Gensini, “A more meaningful scoring system for determining the severity of coronary heart disease,” The American Journal of Cardiology, vol. 51, no. 3, p. 606, 1983.View at: Google Scholar
L. N. Tseng, Y. H. Tseng, Y. D. Jiang et al., “Prevalence of hypertension and dyslipidemia and their associations with micro- and macrovascular diseases in patients with diabetes in Taiwan: an analysis of nationwide data for 2000–2009,” Journal of the Formosan Medical Association, vol. 111, no. 11, pp. 625–636, 2012.View at: Publisher Site | Google Scholar
J. G. Canto, C. I. Kiefe, W. J. Rogers et al., “Atherosclerotic risk factors and their association with hospital mortality among patients with first myocardial infarction (from the National Registry of Myocardial Infarction),” The American Journal of Cardiology, vol. 110, no. 9, pp. 1256–1261, 2012.View at: Publisher Site | Google Scholar
K.-W. Liang, W.-J. Lee, W.-L. Lee, C.-T. Ting, and W. H.-H. Sheu, “Decreased ratio of high-molecular-weight to total adiponectin is associated with angiographic coronary atherosclerosis severity but not restenosis,” Clinica Chimica Acta, vol. 405, no. 1-2, pp. 114–118, 2009.View at: Publisher Site | Google Scholar
D. Ardissino, C. Berzuini, P. A. Merlini et al., “Influence of 9p21.3 genetic variants on clinical and angiographic outcomes in early-onset myocardial infarction,” Journal of the American College of Cardiology, vol. 58, no. 4, pp. 426–434, 2011.View at: Google Scholar
J. B. Muhlestein and J. L. Anderson, “The 9p21.3 genetic region and coronary heart disease where do we go from here?” Journal of the American College of Cardiology, vol. 58, no. 4, pp. 435–437, 2011.View at: Google Scholar
P. Wilkinson, J. Sayer, K. Laji et al., “Comparison of case fatality in south Asian and white patients after acute myocardial infarction: observational study,” British Medical Journal, vol. 312, no. 7042, pp. 1330–1333, 1996.View at: Google Scholar
G. Goliasch, M. E. Kleber, B. Richter et al., “Routinely available biomarkers improve prediction of long-term mortality in stable coronary artery disease: the Vienna and Ludwigshafen Coronary Artery Disease (VILCAD) risk score,” European Heart Journal, vol. 33, no. 18, pp. 2282–2289, 2012.View at: Publisher Site | Google Scholar
C. P. Wen, S. P. Tsai, and W.-S. I. Chung, “A 10-year experience with universal health insurance in Taiwan: measuring changes in health and health disparity,” Annals of Internal Medicine, vol. 148, no. 4, pp. 258–267, 2008.View at: Google Scholar
S. De Henauw, P. De Smet, W. Aelvoet, M. Kornitzer, and G. De Backer, “Misclassification of coronary heart disease in mortality statistics. Evidence from the WHO-MONICA Ghent-Charleroi study in Belgium,” Journal of Epidemiology and Community Health, vol. 52, no. 8, pp. 513–519, 1998.View at: Google Scholar