Email updates

Keep up to date with the latest news and content from BMC Medical Genetics and BioMed Central.

Open Access Highly Accessed Research article

Assessment of the 9p21.3 locus in severity of coronary artery disease in the presence and absence of type 2 diabetes

Natalia V Rivera1*, Robert Carreras-Torres2, Roberta Roncarati13, Chiara Viviani-Anselmi1, Francesca De Micco1, Alessandra Mezzelani4, Werner Koch5, Petra Hoppmann5, Adnan Kastrati5, Alexandre FR Stewart6, Li Chen6, Robert Roberts6, Lennart C Karssen7, Najaf Amin7, Valentina Trimarco8, Raffaele Izzo9, Guido Iaccarino110, Gerolama Condorelli11, Annibale A Puca110, Paolo Pagnotta12, Flavio Airoldi1, Bruno Trimarco9, Cornelia M van Duijn7, Gianluigi Condorelli123 and Carlo Briguori13

Author Affiliations

1 IRCCS Multimedica, Via Fantoli 16/15, 20138, Milan, Italy

2 Facultat de Biologia, Universitat de Barcelona, Barcelona, Spain

3 Institute of Genetic and Biomedical Research (IRGB), National Research Council (CNR), Milan, Italy

4 Institute of Biomedical Technologies (ITB), National Research Council (CNR), Milan, Italy

5 Deutsches Herzzentrum München, Munich, Germany

6 Ruddy Canadian Cardiovascular Genetics Centre, University of Ottawa Heart Institute, Ottawa, Canada

7 Department of Genetic Epidemiology, Erasmus MC, Rotterdam, The Netherlands

8 Department of Neuroscience, Federico II University, Federico II University Hospital, Naples, Italy

9 Department of Clinical Medicine and Cardiovascular Sciences, Federico II University Hospital, Naples, Italy

10 School of Medicine, University of Salerno, Baronissi, Salerno, Italy

11 Department of Biology, Cellular and Molecular Pathology, University Federico II, Naples, Italy

12 Humanitas Clinical and Research Center, Milan, Rozzano, Italy

13 Laboratory of Interventional Cardiology, Department of Cardiology, Clinica Meditterranea of Naples, Via Orazio 2, 80121, Naples, Italy

For all author emails, please log on.

BMC Medical Genetics 2013, 14:11  doi:10.1186/1471-2350-14-11


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2350/14/11


Received:30 March 2012
Accepted:22 October 2012
Published:23 January 2013

© 2013 Rivera et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The 9p21.3 locus is strongly associated with the risk of coronary artery disease (CAD) and with type 2 diabetes (T2D). We investigated the association of 9p21.3 variants with severity of CAD (defined by the number of vessel diseased [VD]) in the presence and absence of T2D.

Methods

We tested 11 9p21.3-variants for association in a white Italian study (N = 2,908), and carried out replication in 2 independent white populations, a German study (N = 2,028) and a Canadian Study (N=950). SNP association and permutation analyses were conducted.

Results

We identified two 9p21.3-variants, rs4977574 (P < 4×10-4) and rs2383207 (P < 1.5×10-3) that were associated with severity of CAD in subjects without T2D. Association of rs4977574 with severity of CAD was confirmed in the Canadian Study. Results from subgroup analysis among patients with T2D showed an interaction between rs10738610 and T2D with P = 4.82×10-2. Further investigation showed that rs10738610 (P < 1.99×10-2) was found to be significantly associated with severity of CAD in subjects with T2D.

Conclusions

The 9p21.3 locus is significantly associated with severity of CAD. The number of associations of 9p21.3 variants with severity of CAD is variable to the presence and absence of T2D. In a CAD-susceptible region of 115 kb, there is only one variant associated with the severity of coronary vessel disease in the presence of type 2 diabetes.

Keywords:
Severity of CAD; Coronary artery disease; Diabetes mellitus; T2D; 9p21.3; Genetics; Single nucleotide polymorphism

Background

Coronary artery disease (CAD) concomitant with type 2 diabetes (T2D) is the leading cause of cardiovascular (CV) morbidity and mortality. Severity of CAD is higher in diabetics than in non-diabetics, often presenting itself as a 3-vessel disease with a more diffuse coronary artery involvement. Thus, diabetics have higher risks for complications frequently leading to death [1].

The combination of genetics, intermediate phenotypes, and environmental factors (such life-style) play a major role in the pathogenesis of CAD and T2D, including determining the conditions and stability of these diseases. The genetic architectures of CAD and T2D are dependent on multiple genetic and environmental triggers that may overlap or may be specific to each disease. The estimated heritability of CAD ranges from 30 to 60% [2,3], whereas of T2D is from 40 to 90% [4,5].

To date, the 9p21.3 locus is the most consistently replicated genetic locus for CAD [6-11] and T2D [12-14]. The association with this locus is independent of cardiovascular risk factors, including age, gender, obesity, smoking, hypertension, and hyperlipidemia. Moreover, the 9p21.3 locus is phenotypically heterogeneous, since it is also associated with various phenotypes, including myocardial infarction [15], arterial stiffness [16], aortic aneurysm [17], ischemic stroke [18,19], glioma [20], and malignant melanoma [21]. The role of the susceptible region for CAD and T2D at the 9p21.3 locus, a region located on two adjacent haplotype blocks, a 44-kb (from rs10116277 to rs1333049) and a 4-kb block (from rs10965243 to rs10757283) [22], has been a major question of research for determining its implication in the pathogenesis of these diseases. Clearly, identifying genetic variants associated with severity of CAD in the presence and absence of T2D is crucial for understanding to what extent genetic variations at the 9p21.3 locus interplay in the pathogenesis of these conditions.

This study was undertaken to investigate the association of 9p21.3-variants with severity of CAD in patients without and with T2D. Thus, first, we sought to identify whether 9p21.3 variants have a role in determining the severity of CAD; and second, we assessed the impact of T2D in the associations with severity of CAD.

Methods

Subjects and phenotype characterization

Participants were of white European ancestry. All studies had protocols approved by local institutional boards. Participants provided written consent and gave permission to use their DNA for research purposes. Baseline characteristics for all participating studies are presented in Additional file 1: Table S1. Baseline measures of clinical and demographic characteristics were measured at each medical facility at the time of the participants’ enrollment in the study.

Additional file 1. Table S1. Demographics and comorbidities of participants. A.Italian study (N=2,908). B.German study (N=2,028). C.Canadian study (N=950). Table S2. Patient characteristics by 9p21.3 (rs4977574, rs2383207, and rs10738610) genotypes for severity of CAD. Table S3. Genotyping information of genotyped and imputed SNP in the Italian study (N=2,908). Table S4. A.Observed genotypes and allele counts (proportions) in the Italian study (N=2,908). B.Observed genotypes and allele counts (proportions) in the German study (N=2,028). C. Observed genotypes and allele counts (proportions) in the Canadian study (N=950). Table S5. Association results derived from linear regression between severity of CAD and 9p21.3 SNPs with interaction term (T2D * SNP) performed in the Italian study (N=2,908). Table S6. Association results of rs4977574, rs2383207, and rs10738610 with other quantifying scores for severity of CAD in the Italian (N=2,908) study. Table S7. Linkage disequilibrium (LD) among the 11 9p21.3 SNPs in the Italian study (N=2,908).

Format: DOC Size: 345KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

In all participating studies, severity of CAD was defined according to angiograms showing number of vessels disease ≥ 1 (namely, 1-vessel, 2-vessel, and 3-vessel disease) with > 50% lumen diameter stenosis and in accordance with clinical standards and guidelines for coronary angiography. T2D was defined as treatment with insulin or oral hypoglycemic drugs, or as elevated (> 126 mg/dL or > 7.0 mmol/L) levels of fasting, non-stressed blood glucose on at least two separate occasions in conjunction with adhering to ongoing dietary measures to control glucose level. Early-onset CAD was defined as onset ≤ 55 years of age and late-onset CAD as otherwise. Type 1 diabetics were excluded in all studies.

The Italian study consists of 2,908 subjects (67.2% male) aged 20 to 90 years diagnosed with and without severity of CAD. Specifically, 1,696 subjects were recruited from both the Clinica Mediterranea in Naples and the Multimedica Hospital in Milan and 1,212 were recruited from the Campania Salute (CS) Network at the University Federico II in Naples. Subjects from the Clinica Mediterranea and the Multimedica Hospital were scheduled for first-time elective coronary angiography because of clinical suspicion and/or diagnoses of CAD. Herein, an angiogram was called “normal” (or “unaffected without severity of CAD”) if it showed number of vessels <1 and smooth lumens with stenosis < 50%, whereas if it showed otherwise, the angiogram was called “affected with severity of CAD” and was further characterized by objective coronary angiography and in accordance with the ACC/AHA guidelines for coronary angiography [23]. Within the 1,696 subjects, 1,541 were “affected with severity of CAD” and 155 were “unaffected without severity of CAD”. The remaining 1,212 subjects were drawn from the Campania Salute (CS) Network, an open registry that collects information from a network of general practitioners and hospitals linked with a coordinating center (University Federico II in Naples). These subjects were used as “controls” (or “unaffected without severity of CAD”) and were selected based on their medical history and absent diagnosis for cardiovascular disease predefined by the exclusion criteria of the open registry study design. Exclusion criteria of the CS Network included history of previous myocardial infarction or angina or procedures of coronary revascularization, stroke, or transitory ischemic attack, congestive heart failure, and chronic kidney disease (defined as GFR < 30 ml/min per 1.73 m2). Additionally, participants enrolled in the CS network were subjected to an active follow-up program including periodically cardiac and carotid ultrasounds aimed to the prevention and the reduction of cardiovascular events in a population of hypertensives. Patients who presented signs of preclinical atherosclerosis or CAD were subjected to further exams, including coronary angiography. Further details on the design of the CS Network can be found elsewhere [24].

Combining all subjects from the three centers, we had a total sample size of 2,908 participants. Within this cohort, 19% of participants were diabetics. Additionally, 15% of patients had early-onset CAD and 31% had late-onset CAD according to available data records.

The German study included 2,028 subjects (78% males) aged 27 to 92 years characterized with ischemic symptoms or evidence of myocardial ischemia in the presence of 50% de novo stenosis located in native coronary vessels. Specifically, 1,917 patients were “affected with severity of CAD” and 111 were “unaffected without severity of CAD”. Severity of vessel disease was recorded as qualitative morphological lesion conditions by means of standard criteria. Baseline, post-procedural, and follow-up angiograms were digitally recorded and assessed off-line in the quantitative angiographic core laboratory (ISAR Center) with an automated edge-detection system (CMS version 7.1; MEDIS Imaging Systems, Leiden, the Netherlands) by two independent operators. Herein, 28% of subjects were diabetics. Further details of this study are described elsewhere [25].

The Canadian study was part of the Ottawa Heart Genomics Study [26] and consisted of 950 subjects (71% males) aged < 55 (males) and < 65 (female) with at least one epicardial stenosis > 50% who were considered as early onset cases of CAD. All subjects were “affected with severity of CAD”. The coronary angiograms were assessed by two independent angiographers blind to the results of the genotype. For a vessel to be scored, stenosis > 50% had to be noted in an epicardial coronary vessel of interest or in one of its major branches. In the event of discordance on the number of vessel scored between two reviewers, angiograms were scored by a third independent reviewer. Modified Gensini and modified Duke coronary scores were derived in subjects who were not previously bypassed. Subjects who had previously undergone PCI for a lesion for which the degree of stenosis was not known were not analyzed. Herein, all subjects were non-diabetics. Exclusion of diabetic patients was part of the study design and their exclusion criteria of the Ottawa Heart Genomics Study. Further details of this study and subject characteristics are described in detailed elsewhere [27].

Specimen collection, SNP selection, and genotyping

A blood sample was drawn from each participant during the regular routing examination for coronary artery angiogram. Genomic DNA was extracted from peripheral blood lymphocytes using the DNA Qiamp Midi kit (QIAGEN Inc. Valencia, CA) and genotyped accordingly.

In the Italian study, 11 SNPs (rs7044859, rs10965215, rs564398, rs7865618, rs10116277, rs4977574, rs2383207, rs10738610, rs10757278, rs1333049, rs10811661) were selected from published genome-wide association studies (GWAS) of CAD [9,11,15] and T2D [13,14] - that is, 7 SNPs were chosen at P = 1×10-8 and 4 SNPs at P < 1×10-6. Herein, genotyping was carried out in two stages; in the first stage 1,696 subjects (from Clinica Mediterranea and Multimedica Hospital) were genotyped using TaqMan® Pre-Designed SNP Genotyping Assays (Applied Biosystems), whereas in the second stage 1,212 subjects (from the CS Network) were genotyped with Illumina Human1M-Duo BeadChip array. In the array, 7 of 11 selected SNPs were present and thus extracted after quality control filtering, which included minor allele frequency (MAF) ≤1%, Hardy-Weinberg equilibrium (HWE) P < 1×10-6, and SNP and subject call rates ≤ 97%. The missing SNPs (rs10116277, rs10738610, rs10757278, rs1333049) were imputed with reference to HapMap CEU release 22 build 36 using the maximum likelihood method implemented in MACH (version 1.0.15). Further details on genotyping and imputation are given in Additional file 1: Table S3. Genotypes of Italian participants were pooled together, resulting in 2,908 subjects and 11 SNPs to analyze.

In the replication stage: in the German study, SNP rs1333049 was available for replication and selection for studying the severity of CAD. Genotyping of this SNP in 2,028 subjects was performed using the TaqMan allelic discrimination assay. Further details of genotyping and quality control were available elsewhere [25], whereas in the Canadian study, SNP rs4977574 was available for replication and selection for studying the severity of CAD. This SNP was extracted from the Affymetrix (Santa Clara, California) 500 K array after GWAS quality control in 950 subjects. Details of genotyping and quality control were available elsewhere [27].

Observed genotypes and allele counts for all participating studies are shown in Additional file 1: Table S4.

Statistical analysis

Statistical analysis within studies

Cardiovascular risk factors (i.e., LDL-cholesterol, HDL-cholesterol, triglycerides, type 2 diabetes mellitus, hypertension, etc.) obtained at baseline were analyzed using SPSS (version 17, SPSS, Chicago, Illinois, USA) and expressed as mean ± SD. Continuous variables were presented as means ± SD and categorical variables as absolute numbers and proportions (%) with one-way analysis of variance and chi-squared tests. These tests were also used to determine difference by genotype for some 9p21.3 variants. Clinical variables deviating from normal distribution were transformed by taking the natural logarithm of their value. Additional quantifying scores for severity of CAD, such as the Gensini score [28] (or modified Gensini [29]), the CAD Prognostic Index (also called the Duke score) [30], modified Duke [31], Dahlen score [32], and percent stenosis diameter (defined by the modified ACC/AHA classification) were assessed with linear regression. Hardy-Weinberg equilibrium (HWE) test was carried out for each SNP in the control group. SNPs with HWE P < 5×10-2 were excluded from the analysis. For the analysis of genetic data, we used the GenABEL [33] library of the R package software.

Association analysis

In the Italian study, we conducted association analyses using linear regression method under an additive genetic model and executed the analyses in two stages; first, an association analysis on the 11 SNPs regressed on severity of CAD (kept as a quantitative variable) using age, sex, and T2D as covariates; and second, an association analysis stratified by T2D on the 11 SNPs regressed on severity of CAD using age and sex as covariates after performing interaction tests between 9p21.3 SNPs and T2D. In both analytical stages, we supported nominal P value significance with empirical P values (i.e., Pemp), which were computed by performing 10,000 permutations on the 9p21.3 region that included the 11 SNPs tested. Association and permutation analyses were performed using methods implemented in the GenABEL [33] library of the R software.

In the replication stage: in the German study, association analysis was carried out using linear regression method on severity of CAD under an additive model and adjusted for age, sex, and T2D. Additionally, association in stratified groups by T2D was performed after testing for interaction between rs1333049 and T2D. Association was performed using GenABEL package [33]. In the Canadian study, the association analysis was conducted using a logistic regression method. Severity of CAD in the early-onset CAD group (as defined by the study design) was dichotomized in two vessel disease groups: (a) 1-VD versus 2- and 3-VD, and (b) 1- and 2-VD versus 3-VD. Logistic regression was performed under an additive genetic model, adjusted for age and sex.

Results

Study samples

The characteristics of the patients participating in these studies are shown in Additional file 1: Table S1. Additionally, baseline characteristics by 9p21.3 genotypes for severity of CAD are shown in Additional file 1: Table S2. In the Italian study, 17% of participants had single-vessel (1VD), 20% had double-vessel (2VD), 16% had triple-vessel (3VD); the remaining 47% had normal coronary arteries or a stenosis diameter of < 50%. Furthermore, 64.5% had severity of CAD without T2D and the remaining percent had severity of CAD with T2D. In the German study, 83.4% of participants had 1VD, 11% had 2VD, and 5.5% had normal coronary vessels. Herein, 68.5% of participants had severity of CAD without T2D and the remaining percent had severe CAD with T2D. In the Canadian study, 30.1% of participants had 1VD, 32.9% had 2VD, and 36.9% had 3VD. In the early onset CAD, 19.8% of participants had 3VD and 35.6% had 1VD and 2VD. In the late onset CAD, 17.1% of participants had 3VD and 27.5% had 1VD and 2VD. Within this study, all participants had severity of CAD without T2D. Results (Table 1) of the association analysis performed on the 11 SNPs in the 2,908 Italian subjects showed that two SNPs (rs4977574 and rs2383207) in high LD (r2 > 0.93) were significantly (Pemp < 0.05) associated with severity of CAD after correcting for multiple testing. Results from testing interaction between 9p21.3 SNPs and T2D are shown in Additional file 1: Table S5. Results (Table 2) of the association analysis stratified by T2D revealed that the association between these two SNPs and severity of CAD continued to be present in subjects without and with T2D. Additionally, a third SNP rs10738610 was found to be significantly (Pemp < 1.99×10-2) associated with severity of CAD in the subjects with T2D. This SNP was imputed in 1,212 controls and had an imputation quality of R2 = 0.96. Imputation quality of imputed SNPs is shown in Additional file 1: Table S3. The linkage disequilibrium (LD, denoted by r2) among the 11 SNPs is shown in Additional file 1: Table S7.

Table 1. Association results derived from linear regression between 11 9p21.3-variants and severity of disease (defined as quantitative trait) in the Italian study (N=2,908)

Table 2. Association results derived from linear regression between 11 9p21.3-variants and severity of disease (defined as quantitative trait) in the Italian study (N=2,908) stratified by T2D

Association tests employing other quantifying angiographic scores, including Duke scores, Dahlen scores, and percent diameter stenosis showed significant evidence of association for rs4977574 and rs2383207 with severity of CAD, as shown in Additional file 1: Table S6.

Results from the replication stage are shown in Tables 3 and 4. In the Canadian study, we observed an association between rs4977574 and severity of CAD in the early-onset CAD patients without T2D, which validated our result for this SNP in this particular group. No genotype data for this SNP was available in the late-onset group. In the German study, we observed no association between rs1333049 and severity of CAD. In addition, no association was detected in the stratified groups by T2D. This finding is in agreement with what we observed for this SNP.

Table 3. Replication results; Association results derived from logistic regression between rs4977574 and severity of CAD without T2D in Italian (N=2,343) and Canadian studies (N=950)

Table 4. Replication results; Association results derived from linear regression between rs1333049 and severity of CAD in the German study (N=2,028)

Discussion

In this study, our aim was to investigate the association between 9p21.3-variants and the severity of CAD in the presence and absence of T2D. Results from our analysis showed that two 9p21.3-variants (rs4977574 and rs2383207) from the known CAD-associated haplotype, are also associated with severity of CAD. We also showed that these variants are associated with severity of CAD in subjects with and without T2D. Additionally, a third SNP rs10738610, located inside the same CAD-associated haplotype, was found to be significantly associated with severity of CAD in subjects with T2D. The significance of the association for this SNP in the pooled analysis was slightly above the threshold; however, this heightened significance was strengthened after stratification by T2D. Association of rs4977574 with severity of CAD was confirmed in an independent study in subjects without T2D. The results of the association for SNP rs1333049, which is a well-documented variant in the association with incidence of CAD and is in high LD with rs494457 (r2 = 0.89), rs2383207 (r2 = 0.89), and rs10738610 (r2 = 0.93), showed no association with severity of CAD in the Italian and German studies. However, Dandona et al. [34] in the Canadian study demonstrated that this SNP predicted the severity of coronary atheromatous burden in patients without T2D.

From our results, we hypothesized that although the 9p21.3 locus does not contain any genes, there are two different pathological mechanisms harboring the same haplotype, one that shares a common genetic platform between molecular mechanisms of CAD and T2D and the other, an independent one that is linked to the pathogenesis of CAD alone. Additionally, disagreement on the associations with severity of CAD of well-documented 9p21.3 variants associated with CAD variants may be due to various factors, including phenotype heterogeneity and difference on risk-allele frequencies among participating studies. The mechanisms of severity of CAD and T2D in patients with acute coronary syndrome (ACS) are complex, and they are not so clear at present. While many genetic studies, including GWAS, for CAD and T2D, have identified the 9p21.3 locus as a common risk region (determined by two adjacent haplotypes), there are no functional variants (non-coding variants) within this region. Therefore, the link between the 9p21.3 locus and the genetic backbones of these two diseases still remains poorly understood. The 9p21.3 locus lies in proximity to CDKN2A and CDKN2B. These cyclin-dependent kinase inhibitors (CDKNs) generate several transcript variants that are involved in diverse biological processes, including cell proliferation, cellular senescence, cell cycle arrest, induction of apoptosis, and activation of caspase activity, to name a few [35,36]. Some of these biological processes are presumably implicated in the pathogenesis of atherosclerosis as well as in the other complex phenotypes (i.e., myocardial infarction, aortic aneurysm, ischemic stroke, T2D, glioma, and malignant melanoma). Recent works have speculated that the pathogenic role of the 9p21.3 locus in atherosclerosis is influenced by gene expression [37,38] and a transcriptional and translational control mechanism presumably caused by cis- and trans-acting elements located in the vicinity of this region. Yet, no cis- and/or trans-acting mechanisms have been identified. Other works have elaborated many discussions about the analysis of the 9p21.3 locus and the appropriate phenotypic characterization. Thus, genotypic variability among 9p21.3 variants may be due to wrong phenotyping rather than sequence variations influenced by gene expression. Moreover, it seems that the significance of variants found at this locus is changeable to the phenotype of interest, its severity, and the presence of combined risk and environmental factors that correlated with the phenotype.

In our study, we have a well-characterized phenotype for severity of CAD in subjects with and without T2D. We performed association analyses in three different white populations of European ancestry, totaling 5,886 subjects. This enabled us to study the relationship between the 9p21.3 locus, including CAD- and T2D-associated variants, and the severity of CAD in subjects with and without T2D. Moreover, to reduce phenotype-genotype bias, we used additional scoring methods for assessing extension of CAD. Then, we conducted association analyses following two approaches for assessing the relationship between severity of CAD and 9p21.3-variants. Our finding showed (1) that rs4977574 and rs2383207 at the 9p21.3 locus are associated with severity of CAD regardless of T2D status and (2) that rs10738610 is associated with severity of CAD comorbid with T2D.

This study has limitations: (1) in order to characterize the 9p21.3 locus more effectively in subjects with severity of CAD and T2D, more variants need to be included in the analysis and hence tested for association with this phenotype; (2) observed differences among genetic effects of 9p21.3 variants across populations may be due to the quantitative method in which angiographic scores were calculated in each participating study; (3) lack of availability of variants rs2383207 and rs10738610 in replication stage demands validation of their associations with severity of CAD in other white population studies with similar characteristics; (4) phenotype heterogeneity among participating studies is likely to be present since each study was previously designed to answer its own research question and had predefined inclusion and exclusion criteria; (5) In the Italian study, although we cannot discard the presence of atherosclerosis in the 1,212 controls, preliminary testing between 9p21.3 variants and incidence of CAD obtained from follow-up records showed no association with CAD; (6) lack of availability of 3-vessel disease subjects in the German study may had underestimated the genetic effect of SNP rs1333049 attained in the replication stage; and (7) lack of availability of genotype data of SNP rs497754 in late-onset CAD subjects with severity of CAD in the Canadian study requires a replication in another cohort with similar characteristics as to strengthen the significance of this SNP with severity of CAD regardless on the onset of the disease.

Conclusions

The 9p21.3 locus is significantly associated with severity of CAD. The number of associations of 9p21.3 variants with severity of CAD is variable to the presence and absence of T2D. In a CAD-susceptible region of 115 kb, there is only one variant associated with the severity of coronary vessel disease in the presence of type 2 diabetes.

Abbreviations

CABG: Coronary artery bypass graft surgery; CAD: Coronary artery disease; T2D: Type 2 diabetes mellitus; eGFR: Estimated glomerular filtration rate; GWAS: Genome-wide association study; HDL-C: High density lipoprotein cholesterol; HWE: Hardy-Weinberg equilibrium; LD: Linkage disequilibrium; LDL-C: Low density lipoprotein cholesterol; Pemp: Empirical P-value; SNP: Single nucleotide polymorphism; kb: kilo base pairs.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the experiments: CB, NVR, GC, RR, and CVA. Performed the clinical work: CB and FA. Performed the biological experiments: RR, CVA, and AM. Analyzed the data: NVR, FDM, LCK, and RCT. Performed study in German population: WK, PH, and AK. Performed study in Canadian population: AFRS, LC, RR. Contributed reagents/materials/analysis tools: RR, CVA, AM, AAP, NVR, P.P. and FDM. Wrote the paper: NVR, RCT, CMVD, NA, AFRS , CB, and GC. Other: drafting article NVR, CB, GC, RR, and CVA. All authors read and approved the final manuscript.

Acknowledgements

Italian Population Study This study was supported by grants from the Italian Ministry of Health and Research (MIUR), Italian Ministry of Health and Region of Lombardy to GC.

Imputations were carried out on the Genetic Cluster Computer (http://www.geneticcluster.org) which is financially supported by the Netherlands Scientific Organization (NWO 480-05-003) along with a supplement from the Dutch Brain Foundation and the VU University Amsterdam.

Canadian Population Study The Ottawa Heart genomics Study was supported by the Canada Foundation for Innovation and the Ontario Research Fund.

German Population Study This study was funded by an institutional grant from the Deutsches Herzzentrum München.

We thank Dr. John Bruce Cantrell and Dr. Michael Latronico for the nice and kind editing of this work.

References

  1. Herlitz J, Wognsen GB, Emanuelsson H, Haglid M, Karlson BW, Karlsson T, Albertsson P, Westberg S: Mortality and morbidity in diabetic and nondiabetic patients during a 2-year period after coronary artery bypass grafting.

    Diabetes Care 1996, 19(7):698-703. PubMed Abstract | Publisher Full Text OpenURL

  2. Zdravkovic S, Wienke A, Pedersen NL, Marenberg ME, Yashin AI, De Faire U: Heritability of death from coronary heart disease: a 36-year follow-up of 20 966 Swedish twins.

    J Intern Med 2002, 252:247-254. PubMed Abstract | Publisher Full Text OpenURL

  3. Schunkert H, Konig IR, Kathiresan S, Reilly MP, Assimes TL, Holm H, Preuss M, Stewart AF, Barbalic M, Gieger C, et al.: Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease.

    Nat Genet 2011, 43(4):333-338. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Zeggini E, Scott LJ, Saxena R, Voight BF, Marchini JL, Hu T, de Bakker PI, Abecasis GR, Almgren P, Andersen G, et al.: Meta-analysis of genome-wide association data and large-scale replication identifies additional susceptibility loci for type 2 diabetes.

    Nat Genet 2008, 40(5):638-645. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Voight BF, Scott LJ, Steinthorsdottir V, Morris AP, Dina C, Welch RP, Zeggini E, Huth C, Aulchenko YS, Thorleifsson G, et al.: Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis.

    Nat Genet 2010, 42(7):579-589. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Assimes TL, Knowles JW, Basu A, Iribarren C, Southwick A, Tang H, Absher D, Li J, Fair JM, Rubin GD, et al.: Susceptibility locus for clinical and subclinical coronary artery disease at chromosome 9p21 in the multi-ethnic ADVANCE study.

    Hum Mol Genet 2008, 17(15):2320-2328. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Meng W, Hughes AE, Patterson CC, Belton C, Kee F, McKeown PP: Chromosome 9p21.3 is associated with early-onset coronary heart disease in the Irish population.

    Dis Markers 2008, 25(2):81-85. PubMed Abstract | Publisher Full Text OpenURL

  8. Schunkert H, Götz A, Braund P, McGinnis R, Tregouet DA, Mangino M, Linsel-Nitschke P, Cambien F, Hengstenberg C, Stark K, et al.: Repeated replication and a prospective meta-analysis of the association between chromosome 9p21.3 and coronary artery disease.

    Circulation 2008, 117(13):1675-1684. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. The Wellcome Trust Case Control Consortium: Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls.

    Nature 2007, 447(7145):661-678. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Samani NJ, Deloukas P, Erdmann J, Hengstenberg C, Kuulasmaa K, McGinnis R, Coronary Artery Disease Consortium, et al.: Large scale association analysis of novel genetic loci for coronary artery disease.

    Arterioscler Thromb Vasc Biol 2009, 29(5):774-780. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. McPherson R, Pertsemlidis A, Kavaslar N, Stewart A, Roberts R, Cox DR, Hinds DA, Pennacchio LA, Tybjaerg-Hansen A, Folsom AR, et al.: A common allele on chromosome 9 associated with coronary heart disease.

    Science 2007, 316(5830):1488-1491. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Broadbent HM, Peden JF, Lorkowski S, Goel A, Ongen H, Green F, Clarke R, Collins R, Franzosi MG, Tognoni G, et al.: Susceptibility to coronary artery disease and diabetes is encoded by distinct, tightly linked SNPs in the ANRIL locus on chromosome 9p.

    Hum Mol Genet 2008, 17(6):806-814. PubMed Abstract | Publisher Full Text OpenURL

  13. Scott LJ, Mohlke KL, Bonnycastle LL, Willer CJ, Li Y, Duren WL, Erdos MR, Stringham HM, Chines PS, Jackson AU, et al.: A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants.

    Science 2007, 316(5829):1341-1345. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Zeggini E, Weedon MN, Lindgren CM, Frayling TM, Elliott KS, Lango H, Timpson NJ, Perry JR, Rayner NW, Freathy RM, et al.: Replication of Genome-Wide Association Signals in UK Samples Reveals Risk Loci for Type 2 Diabetes.

    Science 2007, 316(5829):1336-1341. PubMed Abstract | Publisher Full Text OpenURL

  15. Helgadottir A, Thorleifsson G, Manolescu A, Gretarsdottir S, Blondal T, Jonasdottir A, Jonasdottir A, Sigurdsson A, Baker A, Palsson A, et al.: A common variant on chromosome 9p21 affects the risk of myocardial infarction.

    Science 2007, 316(5830):1491-1493. PubMed Abstract | Publisher Full Text OpenURL

  16. Björck HM, Länne T, Alehagen U, Persson K, Rundkvist L, Hamsten A, Dahlström U, Eriksson P: Association of genetic variation on chromosome 9p21.3 and arterial stiffness.

    J Intern Med 2009, 265(3):373-381. PubMed Abstract | Publisher Full Text OpenURL

  17. Helgadottir A, Thorleifsson G, Magnusson KP, Grétarsdottir S, Steinthorsdottir V, Manolescu A, Jones GT, Rinkel GJ, Blankensteijn JD, Ronkainen A, et al.: The same sequence variant on 9p21 associates with myocardial infarction, abdominal aortic aneurysm and intracranial aneurysm.

    Nat Genet 2008, 40(2):217-224. PubMed Abstract | Publisher Full Text OpenURL

  18. Hu WL, Li SJ, Liu DT, Wang Y, Niu SQ, Yang XC, Zhang Q, Yu SZ, Jin L, Wang XF: Genetic variants on chromosome 9p21 and ischemic stroke in Chinese.

    Brain Res Bull 2009, 79(6):431-435. PubMed Abstract | Publisher Full Text OpenURL

  19. Wahlstrand B, Orho-Melander M, Delling L, Kjeldsen S, Narkiewicz K, Almgren P, Hedner T, Melander O: The myocardial infarction associated CDKN2A/CDKN2B locus on chromosome 9p21 is associated with stroke independently of coronary events in patients with hypertension.

    J Hypertens 2009, 27(4):769-773. PubMed Abstract | Publisher Full Text OpenURL

  20. Shete S, Hosking FJ, Robertson LB, Dobbins SE, Sanson M, Malmer B, Simon M, Marie Y, Boisselier B, Delattre JY, et al.: Genome-wide association study identifies five susceptibility loci for glioma.

    Nat Genet 2009, 41(8):899-904. PubMed Abstract | Publisher Full Text OpenURL

  21. Bishop DT, Demenais F, Iles MM, Harland M, Taylor JC, Corda E, Randerson-Moor J, Aitken JF, Avril MF, Azizi E, et al.: Genome-wide association study identifies three loci associated with melanoma risk.

    Nat Genet 2009, 41(8):920-925. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Silander K, Tang H, Myles S, Jakkula E, Timpson NJ, Cavalli-Sforza L, Peltonen L: Worldwide patterns of haplotype diversity at 9p21.3, a locus associated with type 2 diabetes and coronary heart disease.

    Genome Medicine 2009, 1(5):51. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  23. Scanlon PJ, Faxon DP, Audet AM, Carabello B, Dehmer GJ, Eagle KA, Legako RD, Leon DF, Murray JA, Nissen SE, et al.: ACC/AHA guidelines for coronary angiography. A report of the American College of Cardiology/American Heart Association Task Force on practice guidelines (Committee on Coronary Angiography). Developed in collaboration with the Society for Cardiac Angiography and Interventions.

    J Am Coll Cardiol 1999, 33(6):1756-1824. PubMed Abstract | Publisher Full Text OpenURL

  24. De Luca N, Izzo R, Iaccarino G, Malini PL, Morisco C, Rozza F, Iovino GL, Rao MA, Bodenizza C, Lanni F, et al.: The use of a telematic connection for the follow-up of hypertensive patients improves the cardiovascular prognosis.

    J Hypertens 2005, 23(7):1417-1423. PubMed Abstract | Publisher Full Text OpenURL

  25. Hoppmann P, Erl A, Türk S, Tiroch K, Mehilli J, Schömig A, Kastrati A, Koch W: No Association of Chromosome 9p21.3 Variation With Clinical and Angiographic Outcomes After Placement of Drug-Eluting Stents.

    JACC Cardiovasc Interv 2009, 2(11):1149-1155. PubMed Abstract | Publisher Full Text OpenURL

  26. Stewart AF, Dandona S, Chen L, Assogba O, Belanger M, Ewart G, LaRose R, Doelle H, Williams K, Wells GA, et al.: Kinesin family member 6 variant Trp719Arg does not associate with angiographically defined coronary artery disease in the Ottawa Heart Genomics Study.

    J Am Coll Cardiol 2009, 53(16):1471-1472. PubMed Abstract | Publisher Full Text OpenURL

  27. Dandona S, Roberts R: Creating a genetic risk score for coronary artery disease.

    Curr Atheroscler Rep 2009, 11(3):175-181. PubMed Abstract | Publisher Full Text OpenURL

  28. Gensini GG: A more meaningful scoring system for determining the severity of coronary heart disease.

    Am J Cardiol 1983, 51(3):606-606. PubMed Abstract | Publisher Full Text OpenURL

  29. Montorsi P, Ravagnani PM, Galli S, Salonia A, Briganti A, Werba JP, Montorsi F: Association between erectile dysfunction and coronary artery disease: matching the right target with the right test in the right patient.

    Eur Urol 2006, 50(4):721-731. PubMed Abstract | Publisher Full Text OpenURL

  30. Afonso L, Niraj A, Veeranna V, Fakhry H, Pradhan J: Ethnic and sex differences in disease burden in patients undergoing coronary angiography: the confounding influence of obesity.

    Ethn Dis Winter 2008, 18(1):53-58. OpenURL

  31. Miller JM, Rochitte CE, Dewey M, Arbab-Zadeh A, Niinuma H, Gottlieb I, Paul N, Clouse ME, Shapiro EP, Hoe J, et al.: Diagnostic performance of coronary angiography by 64-Row CT.

    N Engl J Med 2008, 359(22):2324-2336. PubMed Abstract | Publisher Full Text OpenURL

  32. Dahlen GH, Guyton JR, Attar M, Farmer JA, Kautz JA, Gotto AM Jr: Association of levels of lipoprotein Lp(a), plasma lipids, and other lipoproteins with coronary artery disease documented by angiography.

    Circulation 1986, 74(4):758-765. PubMed Abstract | Publisher Full Text OpenURL

  33. Aulchenko Y, Struchalin M: GenABEL: genome-wide SNP association analysis - R package. http://cran.r-project.org/ webcite. 2009

  34. Dandona S, Stewart AF, Chen L, Williams K, So D, O'Brien E, Glover C, Lemay M, Assogba O, Vo L, et al.: Gene dosage of the common variant 9p21 predicts severity of coronary artery disease.

    J Am Coll of Cardiol 2010, 56(6):479-486. Publisher Full Text OpenURL

  35. Harper JW: Cyclin dependent kinase inhibitors.

    Cancer surveys 1997, 29:91-107. PubMed Abstract OpenURL

  36. Fischer PM, Endicott J, Meijer L: Cyclin-dependent kinase inhibitors.

    Progress in cell cycle research 2003, 5:235-248. PubMed Abstract OpenURL

  37. Kalinina N, Agrotis A, Antropova Y, Ilyinskaya O, Smirnov V, Tararak E, Bobik A: Smad expression in human atherosclerotic lesions: evidence for impaired TGF-beta/Smad signaling in smooth muscle cells of fibrofatty lesions.

    Arterioscler Thromb Vasc Biol 2004, 24(8):1391-1396. PubMed Abstract | Publisher Full Text OpenURL

  38. Liu Y, Sanoff HK, Cho H, Burd CE, Torrice C, Mohlke KL, Ibrahim JG, Thomas NE, Sharpless NE: INK4/ARF transcript expression is associated with chromosome 9p21 variants linked to atherosclerosis.

    PLoS One 2009, 4(4):e5027. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2350/14/11/prepub