Next Article in Journal
Resveratrol Inhibition of Cellular Respiration: New Paradigm for an Old Mechanism
Next Article in Special Issue
RORC2 Genetic Variants and Serum Levels in Patients with Rheumatoid Arthritis
Previous Article in Journal
Influence of Long-Distance Bicycle Riding on Serum/Urinary Biomarkers of Prostate Cancer
Previous Article in Special Issue
Sanger Sequencing for BRCA1 c.68_69del, BRCA1 c.5266dup and BRCA2 c.5946del Mutation Screen on Pap Smear Cytology Samples
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 17
Issue 3
10.3390/ijms17030395
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intronic Polymorphisms in the CDKN2B-AS1 Gene Are Strongly Associated with the Risk of Myocardial Infarction and Coronary Artery Disease in the Saudi Population

by
Sayed AbdulAzeez
1,*,
Awatif N. Al-Nafie
2,
Abdullah Al-Shehri
3,
J. Francis Borgio
1,
Ekaterina V. Baranova
4,
Mohammed S. Al-Madan
5,
Rudaynah A. Al-Ali
3,
Fahad Al-Muhanna
3,
Abdullah Al-Ali
6,
Mohammed Al-Mansori
3,
Mohammed Fakhry Ibrahim
3,
Folkert W. Asselbergs
7,
Brendan Keating
8,
Bobby P. C. Koeleman
9 and
Amein K. Al-Ali
10
1
Department of Genetic Research, Institute for Research and Medical Consultation (IRMC), University of Dammam, P.O. Box 1982, 31441 Dammam, Saudi Arabia
2
Department of Pathology, King Fahd Hospital of the University, University of Dammam, 34445 Al-Khobar, Saudi Arabia
3
Department of Internal Medicine, King Fahd Hospital of the University, University of Dammam, 34445 Al-Khobar, Saudi Arabia
4
Department of Pharmaceutical Sciences, Division of Pharmacoepidemiology and Clinical Pharmacology, Utrecht University, 3508 Utrecht, The Netherlands
5
Department of Pediatrics, King Fahd Hospital of the University, University of Dammam, 34445 Al-Khobar, Saudi Arabia
6
Department of Cardiology, King Fahd Hospital, Al-Hafof, 36441 Al-Ahssa, Saudi Arabia
7
Division of Heart & Lungs, Department of Cardiology, University Medical Center Utrecht, 3508 GA Utrecht, The Netherlands
8
Penn Transplant Institute, Hospital of the University of Pennsylvania, Philadelphia, 19104 PA, USA
9
Department of Medical Genetics, University Medical Center Utrecht, 3508 Utrecht, The Netherlands
10
Institute for Research and Medical Consultation (IRMC), University of Dammam, 31441 Dammam, Saudi Arabia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2016, 17(3), 395; https://doi.org/10.3390/ijms17030395
Submission received: 8 February 2016 / Revised: 2 March 2016 / Accepted: 11 March 2016 / Published: 17 March 2016
(This article belongs to the Special Issue Human Single Nucleotide Polymorphisms and Disease Diagnostics)
Browse Figures
Versions Notes

Abstract

:
Recent genome-wide association studies identified single nucleotide polymorphisms (SNPs) on the chromosome 9p21.3 conferring the risk for CAD (coronary artery disease) in individuals of Caucasian ancestry. We performed a genetic association study to investigate the effect of 12 candidate SNPs within 9p21.3 locus on the risk of CAD in the Saudi population of the Eastern Province of Saudi Arabia. A total of 250 Saudi CAD patients who had experienced an myocardial infarction (MI) and 252 Saudi age-matched healthy controls were genotyped using TaqMan assay. Controls with evidenced lack of CAD provided 90% of statistical power at the type I error rate of 0.05. Five percent of the results were rechecked for quality control using Sanger sequencing, the results of which concurred with the TaqMan genotyping results. Association analysis of 12 SNPs indicated a significant difference in the genotype distribution for four SNPs between cases and controls (rs564398 p = 0.0315, χ2 = 4.6, odds ratio (OD) = 1.5; rs4977574 p = 0.0336, χ2 = 4.5, OD = 1.4; rs2891168 p = 1.85 × 10 − 10, χ2 = 40.6, OD = 2.1 and rs1333042 p = 5.14 × 10 − 9, χ2 = 34.1, OD = 2.2). The study identified three protective haplotypes (TAAG p = 1.00 × 10 − 4; AGTA p = 0.022 and GGGCC p = 0.0175) and a risk haplotype (TGGA p = 2.86 × 10 − 10) for the development of CAD. This study is in line with others that indicated that the SNPs located in the intronic region of the CDKN2B-AS1 gene are associated with CAD.

Graphical Abstract

1. Introduction

Coronary artery disease (CAD) is a major cause of morbidity and mortality worldwide and represents a tremendous social and economic burden on society [1]. In Saudi Arabia, the prevalence of CAD is 5.5% in comparison to European Caucasians at a prevalence of 6.4% [2]. Although several independent risk factors for CAD (including smoking, hypercholesterolemia, hypertension, obesity, and diabetes) have been identified, the genetic basis for CAD remains relatively unknown. Recently, genome-wide association studies (GWAS) have reported the increased susceptibility to CAD in carriers of certain single nucleotide polymorphisms (SNPs) within chromosome 9p21.3 locus [3,4,5,6,7]. This locus codes for an anti-sense RNA (CDKN2B-AS1 or ANRIL), which is located nearby the CDKN2A-CDKN2B gene cluster [8,9,10]. The CDKN2A and CDKN2B genes were shown to be significantly associated with an increased risk of CAD (21%) [3,4,11]. Several studies have identified a significant association of specific SNPs in the 9p21.3 locus with CAD [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28]. However, the majority of these studies were primarily conducted on Caucasian, African American, and Chinese populations. Studies conducted on the population of the Eastern Province of Saudi Arabia have reported that 59.2% of sudden death incidences are attributable to CAD [29]. We carried out a case-control study to investigate the association of 12 risk variants at 9p21.3 locus with myocardial infarction (MI) in a Saudi Arabian population for the first time.

2. Results

After obtaining written informed consent from all patients and volunteers, blood samples and clinical data were collected. The clinical characteristics of the study population are presented in Table 1. The study included 250 Saudi CAD patients who had experienced an MI and 252 age matched healthy controls with no history of CAD. These numbers have provided 90% of statistical power at the type I error rate of 0.05 with the maximum odd ratio value of 2. The observed genotype frequencies of the tested polymorphisms in this study obeyed the Hardy–Weinberg equilibrium (Table 2). Initially, the statistical analyses of the association of all the SNPs with CAD patients, who had a secondary condition, either hypertension, type 2 diabetes (T2D), or obesity, were evaluated separately. However, as there were statistical differences between the controls and the patient sub-groups, the CAD patients, regardless of their secondary condition, were pooled prior to data analysis.
The genetic association analysis of 12 SNPs indicated a significant difference in the genotypic distribution of four SNPs between cases and controls (rs564398 p = 0.0315, χ2 = 4.6, OD = 1.5; rs4977574 p = 0.0336, χ2 = 4.5, OD = 1.4; rs2891168 p < 0.0001, χ2 = 40.6, OD = 2.1 and rs1333042 p < 0.0001, χ2 = 34.1, OD = 2.2) (Table 2; Figure 1). The test identified a significant association between GG genotype of rs2891168 (OR: 2.7919; 95%CI (Confidence interval): 1.8201–4.2824; χ2: 23.5464 and p < 0.0001) and CAD using a common genetic model, namely the recessive, dominant, and additive model. Permutation association test of single markers conferred the most significant risk effect of rs289116G on CAD (permutation p-value < 0.0001). An extended RR of risk alleles revealed that the rs2891168G allele has the most significant association (p < 0.0001; RR: 1.4593 CI at 95% 1.2848–1.6575). An LD (linkage disequilibrium) test (D prime (D’) = 0.86) showed that the alleles of rs2891168G and rs4977574G are strongly linked with each other among the significant SNPs in the Saudi Arabian population (Figure 2).
LD plot was constructed based on the pairwise correlation between the 12 SNPs on the 9p21 locus, which has three blocks (blocks 1, 2, and 3) of linked variants (Figure 1). Furthermore, a single group of the four significant SNPs (Sig. SNPs) was made as one block and analyzed for its significant association with the CAD (Figure 1). Haplotype analysis of the significant SNPs, rs564398, rs4977574, rs2891168, and rs1333042 indicated that the haplotype TGGA is a risk factor for CAD. The haplotypes TAAG (rs564398, rs4977574, rs2891168, and rs1333042); AGTA (rs523096, rs518394, rs564398, and rs7865618); and GGGCC (rs1333042, rs2383207, rs10757278, rs1333048, and rs1333049) are CAD protective (p < 0.0001, 0.022, and 0.0175, respectively) (Figure 1 and Table 3). The protective haplotypes in Block 1 (AGTA p = 0.022) and in Block 3 (GGGCC p = 0.0175) with more than 30% frequency were observed in controls (Figure 1 and Table 3). Furthermore, the most significant risk haplotype TGG (rs10757272:rs4977574:rs2891168 p < 0.0001) of Block 2 had a frequency of 56.8% in patients and 34.5% in controls.
To evaluate whether the SNPs identified to be associated with CAD (rs564398, rs4977574, rs2891168, and rs1333042) enhance the predictive value of known conventional risk factors such as gender, age, and body mass index (BMI), two multivariate models were built. The first model (the clinical only model) consisted of known CAD risk factors collected in this study: gender, age, and BMI. For the second model (clinical + genetic model), the four CAD-associated SNPs (rs564398, rs4977574, rs2891168, and rs1333042) were entered into the model assuming an additive model of inheritance, in addition to those included in the clinical-only model (Table 4). SNPs rs4977574 and rs564398 were not significantly associated with CAD (p > 0.05) after adjusting for other variables in the model and therefore were removed from the clinical + genetic model. The final clinical + genetic model included age, gender, BMI, rs1333042, and rs2891168. As shown in Table 4, both models were significantly predictive of CAD, with area under curve (AUC) of 0.79 (95% CI: 0.73–0.84) for clinical only model (p < 10−10), and AUC of 0.87 (95% CI: 0.82–0.90) for clinical + genetic model. Most importantly, two CAD-associated SNPs, rs1333042, and rs2891168, improved the predictive power for CAD over the model composing of only conventional known risk factors, with an improvement in AUC of 0.08 (95% CI: 0.04–0.12, p = 0.000236).

3. Discussion

During the last five decades, Saudi Arabia has undergone tremendous socio-economic changes, which have resulted in an increase in the prevalence of common diseases, such as diabetes, cardiovascular diseases, and obesity [1,2]. However, the high prevalence of cardiovascular diseases in the Kingdom cannot be solely attributed to these socio-economic changes. It is believed that genetics and ion channels are contributory factors in the increased prevalence rate of these diseases [30,31]. Recent genome-wide association studies have indicated that there is an association between the increased susceptibility to CAD and specific single nucleotide polymorphisms (SNPs) within the genome [3,4,5,6,7] which play a role in conjunction with other known traditional CAD risk factors. One such locus has been identified within chromosome 9p21.3 which codes for an anti-sense RNA (CDKN2B-AS1 or ANRIL) and is located near the CDKN2A-CDKN2B gene cluster [3,7,8,9,10]. The variants within this locus were shown to be significantly associated with an increased risk of CAD and diabetes in populations of European origin and some other populations [3,4,11]. It is of great importance that these studies are replicated, especially in populations that have distinctive genetic backgrounds, such as that in the Eastern Province of Saudi Arabia, in order to understand the pathophysiology of CAD. This is the first study undertaken on an Arab population to replicate previous studies conducted on the association of these variants with CAD in other populations [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]. Type 2 diabetes mellitus, dyslipidemia, and hypertension were considered separately for the significance analysis. These three parameters did not show any significant association during the non-random association of 12 alleles.
Four of the 12 SNPs tested in the present study were significantly associated with CAD (p < 0.05). One of the four SNPs (rs564398) has been shown to be moderately associated with T2D in Israeli population [32]. The high OR of the risk alleles harboring these four SNPs may be due to the high percentage of consanguinity in the Eastern Province of Saudi Arabia, which increases the possibility of inheriting a recessive allele [33,34,35,36,37].
Haplotype analysis of the rs564398, rs4977574, rs2891168, and rs1333042 indicated that the haplotype TGGA is a risk factor for CAD, while the haplotypes TAAG (rs564398, rs4977574, rs2891168, and rs1333042); AGTA (rs523096, rs518394, rs564398, and rs7865618); and GGGCC (rs1333042, rs2383207, rs10757278, rs1333048, and rs1333049) are CAD protective (p < 0.0001, 0.022, and 0.0175, respectively). These haplotypes are reported here for the first time in an Arab population. However, the SNPs of haplotypes TAAG and TGGA (rs564398, rs4977574, rs2891168, and rs1333042) have been found to be associated with CAD in other populations, including European Hispanic, Chinese, white American, and Turkish populations [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27].
In comparing two multivariate models built with and without CAD-associated SNPs, we found these risk SNPs adding predictive value of CAD on top of known risk factors such as age, gender, and BMI. Due to limited collection of clinical variables, traditional risk factors commonly included in the CAD risk prediction model, such as systolic blood pressure and blood lipid levels HDL (High-density lipoprotein), LDL (Low-density lipoprotein), and TC (total cholesterol), were not evaluated in the model building process. Further studies that collect these clinical risk factors in order to evaluate the improvement in prediction value from risk SNPs are needed.

4. Materials and Methods

4.1. Study Population

The study, which was conducted from 2012–2014, included 250 Saudi CAD patients who had experienced an MI, (STEMI and NSTEMI) as diagnosed by the medical history such as, signs and symptoms (like fatigue), chest discomfort (with physical examination), in addition to electrocardiogram and high troponin level. The patients were selected randomly from major hospitals in Al-Ahssa Qatif and Al-Khobar in the Eastern Province of Saudi Arabia. Two-hundred and fifty-two age-matched healthy volunteers attending the blood banks of the same hospitals were included in the study as controls. The majority of patients attending these hospitals either live in these areas or originate from these areas, with very little population admixture due to the high rate of consanguinity. This has been confirmed by our recent research results [33,34]. Stringent inclusion criteria for both cases and controls were followed as a result of the high incident rate of CAD (11.7%) in this age group (>40 years) in the Saudi population [38,39,40,41,42]. CAD diagnosis was based on electrocardiogram, echocardiography, blood tests and coronary angiography and cardiac catheterization. Exclusion criteria for the control group included a present or past or family history of CAD. Peripheral blood samples (5 mL) were collected for genetic analysis. SNPs selected for inclusion in the study have been previously reported to be strongly associated with CAD [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28].
This study was approved by the Institutional Review Board and Committee for Biological and Medical Ethics, University of Dammam. Written informed consent was obtained from all study participants.

4.2. TaqMan® Single Nucleotide Polymorphisms (SNPs) Genotyping

The genomic DNA was extracted from peripheral blood samples of patients and controls using a standard extraction kit (Wizard Genomic DNA purification kit, Promega, Madison, WI, USA). The Rotor-Gene Q (Qiagen, Hilden, Germany) real time polymerase chain reaction (PCR) system was used to perform SNP genotyping using TaqMan® assay (Applied Biosystems, Foster City, CA, USA). Genotyping of 12 SNPs (rs1333042, rs2891168, rs4977574, rs10757278, rs2383207, rs523096, rs518394, rs564398, rs7865618, rs1333049, rs10757272, and rs1333048) was carried out separately in a 15 μL volume containing 20 ng of genomic DNA, 6.25 μL of 2× TaqMan master mix, 1.25 μL of 20× assay working stock, and nuclease free water to 15 μL. An allele calling report was generated with Rotor-Gene Q Series Software 2.0.2 (Build 4) (Corbett Life Science, Qiagen).

4.3. Quality for SNP Genotyping—Sanger Sequencing

The initial TaqMan genotyping results were verified through Sanger sequencing on 5% of the samples. Oligonucleotide primers were designed for the PCR amplification product of ~150 bp up and downstream to the position of each SNP (primers are available on request) and were amplified separately. The PCR mixture contained 100 ng of DNA, top Taq buffer (10×) (5 μL), 25 mM MgCl2 (3 µL), 25 mM dNTP (0.4 µL), top Taq DNA polymerase 5 U/µL (0.2 µL), 10 µM forward and reverse (respective primers, accordingly) oligonucleotides (2 µL), and water (to 50 µL). The temperature profile was as follows: 5 min at 95 °C, 30 PCR cycles of 95 °C/30 s, 56 °C/30 s, and 72 °C/30 s, final extension of 7 min at 72 °C. The amplicons were purified using a standard PCR purification kit (QIAquick PCR Purification Kit, Qiagen). Cycle sequencing was performed in a BIO-RAD MyCycler™ (Bio-Rad, Hercules, CA, USA) with a total volume of 20 μL. Purified products (10 μL) were separated in a Genetic Analyzer 3500 DNA sequencer (Applied Biosystems). Data were analyzed using Sequencing Analysis Software V5.4 (Applied Biosystems). MAFFT version 7 was used for multiple sequence alignment. The quality for SNP genotyping through TaqMan® assay was assured by direct DNA sequence analysis of 5% of the samples. Furthermore, the DNA samples for controls and cases were analyzed in the same batches.

4.4. Statistical Analysis

Statistical power analysis to ascertain association study was carried out to verify the adequate power of the selected sample size using the Sampsize online tool [42]. Clinical characteristics of the study participants were recorded upon inclusion. Differences in clinical characteristics between cases and controls were calculated by the two-sample t-test or the χ2 test if appropriate, odds ratios (OR) and relative risks (RR) for the association of risk alleles with CAD were calculated and 95% confidence intervals (CI) were constructed [43]. The allelic association of SNPs on cases and controls were assessed using the χ2 test. The Hardy-Weinberg equilibrium test was applied to confirm the independent segregation of the alleles. Non-random association of 12 alleles at 9p21 was tested by LD test. Haplotype blocks were constructed from the genotype data. To understand how genotypes of each SNP alleles are associated with CAD, we used a genetic model (Recessive, Dominant, and Additive model) in the association analysis. SNP haplotype analysis was extended to reveal the risk of SNP patterns on the development of CAD by treating each haplotype as a single variant and the other haplotypes as an alternative allele. Statistical significance was set at p ≤ 0.05. A permutation p-value for a single marker only was also calculated. All genetic analyses were performed using Haploview, version 4.2 [44].
Multivariate models to discriminate CAD cases and healthy controls were built using logistic regression. Demographic and patient characteristic variables (gender, age, and BMI) that were significantly associated with CAD were included in the model building (clinical only model). Individual SNPs that were found to be associated with CAD were then added to clinical model. The variables that were not significant after adjusting for other variables in the model were then removed, resulting in the final most parsimonious model including both demographic/clinical variables and SNPs (clinical + genetic model). The improvement in model performance of the clinical + genetic model compared to clinical only model was evaluated as outlined by Delong [45].

5. Conclusions

Four SNPs tested in the present study were significantly associated with CAD (p < 0.05). Haplotype analysis indicated that the haplotypes TGGA is a risk factor for CAD, while the haplotypes TAAG, AGTA, and GGGCC are CAD protective in the Saudi population of the Eastern Province. The study clearly replicates the findings of previous studies conducted on the association of these polymorphisms and CAD.

Acknowledgments

We thank Ranilo Moring Tumbaga, Horace Tumulak Pacifico, Florentino Jr Catuguran Mata, and Mohammed Hassan Al Shamlan for their technical expertise. The study was partially funded by the Deanship of Scientific Research, University of Dammam (To Sayed AbdulAzeez Ref No.: 2013261).

Author Contributions

For this research samples were collected by Awatif N. Al-Nafie, Abdullah Al-Shehri, Fahad Al-Muhanna; Abdullah Al-Ali, Mohammed Al-Mansori, Mohammed S. Al-Madan, Mohammed Fakhry Ibrahim; the experiments were designed by Amein K. Al-Ali, Sayed AbdulAzeez, J. Francis Borgio; the experiments were performed by Sayed AbdulAzeez, J. Francis Borgio; analysis was done by Sayed AbdulAzeez, J. Francis Borgio, Ekaterina V. Baranova, Folkert W. Asselbergs, Rudaynah A. Al-Ali, Brendan Keating, Bobby P.C. Koeleman. partial reagents/materials provided by Amein K. Al-Ali.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zou, J.G.; Ma, Y.T.; Xie, X.; Yang, Y.N.; Pan, S.; Adi, D.; Liu, F.; Chen, B.D. The association between CYP1A1 genetic polymorphisms and coronary artery disease in the Uygur and Han of China. Lipids Health Dis. 2014, 5. [Google Scholar] [CrossRef] [PubMed]
  2. Go, A.S.; Mozaffarian, D.; Roger, V.L.; Benjamin, E.J.; Berry, J.D.; Borden, W.B.; Bravata, D.M.; Dai, S.; Ford, E.S.; Fox, C.S.; et al. Executive summary. Circulation 2013, 127, 143–152. [Google Scholar] [CrossRef] [PubMed]
  3. Wellcome Trust Case Control Consort. Genome-wide association study of 14,000 cases of seven common diseases and 3000 shared controls. Nature 2007, 447, 661–678. [Google Scholar]
  4. Samani, N.J.; Erdmann, J.; Hall, A.S.; Hengstenberg, C.; Mangino, M.; Mayer, B.; Dixon, R.J.; Meitinger, T.; Braund, P.; Wichmann, H.E.; et al. Genomewide association analysis of coronary artery disease. N. Engl. J. Med. 2007, 357, 443–453. [Google Scholar] [CrossRef] [PubMed]
  5. Abdullah, K.G.; Li, L.; Shen, G.Q.; Hu, Y.; Yang, Y.; MacKinlay, K.G.; Topol, E.J.; Wang, Q.K. Four SNPS on chromosome 9p21 confer risk to premature, familial CAD and MI in an American Caucasian population (GeneQuest). Ann. Hum. Genet. 2008, 72, 654–657. [Google Scholar] [CrossRef] [PubMed]
  6. Durda, P.; Sabourin, J.; Lange, E.M.; Nalls, M.A.; Mychaleckyj, J.C.; Jenny, N.S.; Li, J.; Walston, J.; Harris, T.B.; Psaty, B.M.; et al. Plasma levels of soluble interleukin-2 receptor α associations with clinical cardiovascular events and genome-wide association scan. Arterioscler. Thromb. Vasc. Biol. 2015, 35, 2246–2253. [Google Scholar] [CrossRef] [PubMed]
  7. Roberts, R.; Stewart, A.F. 9p21 and the genetic revolution for coronary artery disease. Clin. Chem. 2012, 58, 104–112. [Google Scholar] [CrossRef] [PubMed]
  8. McPherson, R.; Pertsemlidis, A.; Kavaslar, N.; Stewart, A.; Roberts, R.; Cox, D.R.; Hinds, D.A.; Pennacchio, L.A.; Tybjaerg-Hansen, A.; Folsom, A.R.; et al. A common allele on chromosome 9 associated with coronary heart disease. Science 2007, 16, 1488–1491. [Google Scholar] [CrossRef] [PubMed]
  9. Broadbent, H.M.; Peden, J.F.; Lorkowski, S.; Goel, A.; Ongen, H.; Green, F.; Clarke, R.; Collins, R.; Franzosi, M.G.; 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. Gen. 2008, 17, 806–814. [Google Scholar] [CrossRef] [PubMed]
  10. Szpakowicz, A.; Pepinski, W.; Waszkiewicz, E.; Maciorkowska, D.; Skawronska, M.; Niemcunowicz-Janica, A.; Milewski, R.; Dobrzycki, S.; Musial, W.J.; Kaminski, K.A. Polymorphism of 9p21. 3 locus is associated with 5-year survival in high-risk patients with myocardial infarction. PLoS ONE 2013, 8, e72333. [Google Scholar] [CrossRef] [PubMed]
  11. 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, 8, 1491–1493. [Google Scholar] [CrossRef] [PubMed]
  12. Björck, H.M.; 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, 373–381. [Google Scholar] [CrossRef] [PubMed]
  13. Gori, F.; Specchia, C.; Pietri, S.; Crociati, L.; Barlera, S.; Franciosi, M.; Nicolucci, A.; Signorini, S.; Brambilla, P.; Franzosi, M.G. Common genetic variants on chromosome 9p21 are associated with myocardial infarction and type 2 diabetes in an Italian population. BMC Med. Genet. 2010, 11. [Google Scholar] [CrossRef] [PubMed]
  14. Qi, L.; Ma, J.; Qi, Q.; Hartiala, J.; Allayee, H.; Campos, H. Genetic risk score and risk of myocardial infarction in Hispanics. Circulation 2011, 123, 374–380. [Google Scholar] [CrossRef] [PubMed]
  15. Saade, S.; Cazier, J.B.; Ghassibe-Sabbagh, M.; Youhanna, S.; Badro, D.A.; Kamatani, Y.; Hager, J.; Yeretzian, J.S.; El-Khazen, G.; Haber, M.; et al. Large scale association analysis identifies three susceptibility loci for coronary artery disease. PLoS ONE 2011, 6, e29427. [Google Scholar] [CrossRef] [PubMed]
  16. Qi, L.; Li, J.M.; Sun, H.; Huang, X.Q.; Lin, K.Q.; Chu, J.Y.; Yang, Z.Q. Association between gene polymorphisms and myocardial infarction in Han Chinese of Yunnan province. Chin. J. Med. Genet. 2012, 29, 413–419. [Google Scholar]
  17. Gränsbo, K.; Almgren, P.; Sjögren, M.; Smith, J.G.; Engström, G.; Hedblad, B.; Melander, O. Chromosome 9p21 genetic variation explains 13% of cardiovascular disease incidence but does not improve risk prediction. J. Intern. Med. 2013, 274, 233–240. [Google Scholar] [CrossRef] [PubMed]
  18. Lee, J.Y.; Lee, B.S.; Shin, D.J.; Park, K.W.; Shin, Y.A.; Kim, K.J.; Heo, L.; Lee, J.Y.; Kim, Y.K.; Kim, Y.J.; et al. A genome-wide association study of a coronary artery disease risk variant. J. Hum. Genet. 2013, 58, 120–126. [Google Scholar] [CrossRef] [PubMed]
  19. Lövkvist, H.; Sjögren, M.; Höglund, P.; Engström, G.; Jern, C.; Olsson, S.; Smith, J.G.; Hedblad, B.; Andsberg, G.; Delavaran, H.; et al. Are 25 SNPs from the CARDIoGRAM study associated with ischaemic stroke? Eur. J. Neurol. 2013, 20, 1284–1291. [Google Scholar] [CrossRef] [PubMed]
  20. Sakalar, C.; Gurbuz, E.; Kalay, N.; Kaya, M.G. Higher frequency of rs4977574 (the G Allele) on chromosome 9p21. 3 in patients with myocardial infarction as revealed by PCR-RFLP analysis. Tohoku J. Exp. Med. 2013, 230, 171–176. [Google Scholar] [CrossRef] [PubMed]
  21. Tragante, V.; Doevendans, P.A.; Nathoe, H.M.; van der Graaf, Y.; Spiering, W.; Algra, A.; de Borst, G.J.; de Bakker, P.I.; Asselbergs, F.W. The impact of susceptibility loci for coronary artery disease on other vascular domains and recurrence risk. Eur. Heart J. 2013, 34, 2896–2904. [Google Scholar] [CrossRef] [PubMed]
  22. van Setten, J.; Isgum, I.; Smolonska, J.; Ripke, S.; de Jong, P.A.; Oudkerk, M.; de Koning, H.; Lammers, J.W.; Zanen, P.; Groen, H.J.; et al. Genome-wide association study of coronary and aortic calcification implicates risk loci for coronary artery disease and myocardial infarction. Atherosclerosis 2013, 228, 400–405. [Google Scholar] [CrossRef] [PubMed]
  23. Hamrefors, V.; Hedblad, B.; Hindy, G.; Smith, J.G.; Almgren, P.; Engström, G.; Sjögren, M.; Gränsbo, K.; Orho-Melander, M.; Melander, O. Smoking modifies the associated increased risk of future cardiovascular disease by genetic variation on chromosome 9p21. PLoS ONE 2014, 9, e85893. [Google Scholar] [CrossRef] [PubMed]
  24. Huang, Y.; Ye, H.; Hong, Q.; Xu, X.; Jiang, D.; Xu, L.; Dai, D.; Sun, J.; Gao, X.; Duan, S. Association of CDKN2BAS polymorphism rs4977574 with coronary heart disease: A case-control study and a meta-analysis. Int. J. Mol. Sci. 2014, 15, 17478–17492. [Google Scholar] [CrossRef] [PubMed]
  25. Lee, I.T.; Goodarzi, M.O.; Lee, W.J.; Rotter, J.I.; Chen, Y.D.; Liang, K.W.; Lee, W.L.; Sheu, W.H. The chromosome 9p21 variant not predicting long-term cardiovascular mortality in Chinese with established coronary artery disease: An eleven-year follow-up study. Biomed. Res. Int. 2014, 2014. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, Y.; Wang, L.; Liu, X.; Zhang, Y.; Yu, L.; Zhang, F.; Liu, L.; Cai, J.; Yang, X.; Wang, X. Genetic variants associated with myocardial infarction and the risk factors in Chinese population. PLoS ONE 2014, 9, e86332. [Google Scholar] [CrossRef] [PubMed]
  27. Yinko, S.L.; Thanassoulis, G.; Stark, K.D.; Tsadok, M.A.; Engert, J.C.; Pilote, L. ω-3 Fatty acids and the genetic risk of early onset acute coronary syndrome. Nutr. Metab. Cardiovasc. Dis. 2014, 24, 1234–1239. [Google Scholar] [CrossRef] [PubMed]
  28. Yayla, Ç.; Okyay, K.; Yılmaz, A.; Sahinarslan, A.; Seda, A.; Eyiol, A.; Bolayır, A.; Sezenöz, B.; Menevse, S.; Çengel, A. OP-195 the relationship of coronary artery disease with alleles of chromosome 9p21 in Turkish population. Am. J. Cardiol. 2015, 115, S88–S89. [Google Scholar] [CrossRef]
  29. Nofal, H.; Abdulmohsen, M.; Khamis, A. Incidence and causes of sudden death in a university hospital in eastern Saudi Arabia. East. Mediterr. Health 2011, 17, 665–670. [Google Scholar]
  30. Fedele, F.; Mancone, M.; Chilian, W.M.; Severino, P.; Canali, E.; Logan, S.; De Marchis, M.L.; Volterrani, M.; Palmirotta, R.; et al. Role of genetic polymorphisms of ion channels in the pathophysiology of coronary microvascular dysfunction and ischemic heart disease. Basic Res. Cardiol. 2013, 108, 1–12. [Google Scholar] [CrossRef] [PubMed]
  31. Fedele, F.; Severino, P.; Bruno, N.; Stio, R.; Caira, C.; D’Ambrosi, A.; Brasolin, B.; Ohanyan, V.; Mancone, M. Role of ion channels in coronary microcirculation: A review of the literature. Future Cardiol. 2013, 9, 897–905. [Google Scholar] [CrossRef] [PubMed]
  32. Cauchi, S.; Meyre, D.; Durand, E.; Proença, C.; Marre, M.; Hadjadj, S.; Choquet, H.; De Graeve, F.; Gaget, S.; Allegaert, F.; et al. Post genome-wide association studies of novel genes associated with type 2 diabetes show gene-gene interaction and high predictive value. PLoS ONE 2008, 3, e2031. [Google Scholar] [CrossRef] [PubMed]
  33. Akhtar, M.S.; Qaw, F.; Borgio, J.F.; Albuali, W.; Suliman, A.; Nasserullah, Z.; Al-Jarrash, S.; Al-Ali, A. Spectrum of alpha-thalassemia mutations in transfusion-dependent beta-thalassemia patients from the Eastern Province of Saudi Arabia. Hemoglobin 2013, 37, 65–73. [Google Scholar] [CrossRef] [PubMed]
  34. Borgio, J.F.; AbdulAzeez, S.; Al-Nafie, A.N.; Naserullah, Z.A.; Al-Jarrash, S.; Al-Madan, M.S.; Al-Muhanna, F.; Steinberg, M.H.; Al-Ali, A.K. A novel HBA2 gene conversion in cis or trans: “alpha12 allele” in a Saudi population. Blood Cells Mol. Dis. 2014, 53, 199–203. [Google Scholar] [CrossRef] [PubMed]
  35. Becker, S.M.; Al Halees, Z.; Molina, C.; Paterson, R.M. Consanguinity and congenital heart disease in Saudi Arabia. Am. J. Med. Genet. 2001, 99, 8–13. [Google Scholar] [CrossRef]
  36. Al-Nafie, A.N.; Borgio, J.F.; AbdulAzeez, S.; Al-Suliman, A.M.; Qaw, F.S.; Naserullah, Z.A.; Al-Jarrash, S.; Al-Madan, M.S.; Al-Ali, R.A.; AlKhalifah, M.A.; et al. Co-inheritance of novel ATRX gene mutation and globin (α & β) gene mutations in transfusion dependent beta-thalassemia patients. Blood Cells Mol. Dis. 2015, 55, 27–29. [Google Scholar] [PubMed]
  37. Borgio, J.F. Molecular nature of alpha-globin genes in the Saudi population. Saudi Med. J. 2015, 36, 1271–1276. [Google Scholar] [CrossRef] [PubMed]
  38. Al-Nozha, M.M.; Arafah, M.R.; Al-Mazrou, Y.Y.; Al-Maatouq, M.A.; Khan, N.B.; Khalil, M.Z.; Al-Khadra, A.H.; Al-Marzouki, K.; Abdullah, M.A.; Al-Harthi, S.S.; et al. Coronary artery disease in Saudi Arabia. Saudi Med. J. 2004, 25, 1165–1171. [Google Scholar] [PubMed]
  39. Al-Sheikh, S.O.; Aljabri, B.A.; Al-Ansary, L.A.; Al-Khayal, L.A.; Al-Salman, M.M.; Al-Omran, M.A. Prevalence of and risk factors for peripheral arterial disease in Saudi Arabia. A pilot cross-sectional study. Saudi Med. J. 2007, 28, 412–414. [Google Scholar] [PubMed]
  40. Al-Omran, M. Atherosclerotic disease and risk factor modification in Saudi Arabia: A call to action. Vasc. Health Risk Manag. 2012, 8, 349. [Google Scholar] [CrossRef] [PubMed]
  41. Al-Nozha, M.M.; Al-Mazrou, Y.Y.; Arafah, M.R.; Al-Maatouq, M.A.; Khalil, M.Z.; Khan, N.B.; Al-Khadra, A.; Al-Marzouki, K.; Al-Harthi, S.S.; Abdullah, M.; Al-Shahid, M.S. Smoking in Saudi Arabia and its relation to coronary artery disease. J Saudi Heart Assoc. 2009, 2, 169–176. [Google Scholar] [CrossRef] [PubMed]
  42. Sampsize. Available online: http://sampsize. sourceforge.net/iface/s3.html (accessed on 5 July 2015).
  43. MedCalc. Available online: http://www.medcalc.org (accessed on 1 August 2015).
  44. Barrett, J.C.; Fry, B.; Maller, J.D.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 15, 263–265. [Google Scholar] [CrossRef] [PubMed]
  45. De Long, E.R.; De Long, D.M.; Clarke-Pearson, D.L. Comparing the areas under two or more correlated receiver operating characteristic curves: A nonparametric approach. Biometrics 1988, 44, 837–845. [Google Scholar] [CrossRef]
Ijms 17 00395 g001 1024
Figure 1. Upper half: Haploview LD (linkage disequilibrium) plot of the 12 SNPs in 9p21 locus. The pairwise correlation between the single nucleotide polymorphisms (SNPs) were measured as r2 and shown (×100) in each diamond. Lower half: Haplotypes of patient group and controls. Protective haplotypes are in purple boxes and risk haplotypes are in red boxes. Sig. SNPs: Block and haplotypes of significant SNPs.
Figure 1. Upper half: Haploview LD (linkage disequilibrium) plot of the 12 SNPs in 9p21 locus. The pairwise correlation between the single nucleotide polymorphisms (SNPs) were measured as r2 and shown (×100) in each diamond. Lower half: Haplotypes of patient group and controls. Protective haplotypes are in purple boxes and risk haplotypes are in red boxes. Sig. SNPs: Block and haplotypes of significant SNPs.
Ijms 17 00395 g001
Ijms 17 00395 g002 1024
Figure 2. The chromosomal positions of the 12 SNPs analyzed in 9p21 locus and their LD. Upper table: Indicating the coordinates according to the reference sequence NT_008413.18. Middle row: indicates the distance (bp) between two adjacent SNPs. Lower half: Indicating the D’ value in each diamond among 12 SNPs in 9p21.
Figure 2. The chromosomal positions of the 12 SNPs analyzed in 9p21 locus and their LD. Upper table: Indicating the coordinates according to the reference sequence NT_008413.18. Middle row: indicates the distance (bp) between two adjacent SNPs. Lower half: Indicating the D’ value in each diamond among 12 SNPs in 9p21.
Ijms 17 00395 g002
Table
Table 1. Clinical characteristics of coronary artery disease (CAD) cases and controls in a population from Eastern Province, Saudi Arabia.
Table 1. Clinical characteristics of coronary artery disease (CAD) cases and controls in a population from Eastern Province, Saudi Arabia.
Characteristics *Cases (n = 250)Controls (n = 252)
Males, n (%)163 (65.2)200 (79.4)
Females, n (%)87 (34.8)52 (20.6)
Age, (years) #52.3 ± 13.848.1 ± 7.9
BMI, (kg/m2)33.6 ± 11.027.9 ± 4.0
Hypertension, n (%)174 (69.6)0 (0)
Type 2 diabetes n, (%)146 (58.4)1 (0.4)
Dyslipidemia29 (11.6)0 (0)
Unstable angina4 (1.6)0 (0)
* Data are shown as means ± standard deviations (SD) or percentages. # For cases age-at-diagnosis and for controls age-at-collection. Individuals with systolic blood pressure >140 mm Hg and diastolic blood pressure >90 mm Hg were defined as hypertension. BMI: Body mass index.
Table
Table 2. Association of CAD with risk alleles of single nucleotide polymorphisms (SNPs) in 9p21.3 locus.
Table 2. Association of CAD with risk alleles of single nucleotide polymorphisms (SNPs) in 9p21.3 locus.
SNP IDAssoc Allelep-HWObsHETPredHETOdds Ratio (95% CI)Case; Control Ratioχ2p Value
rs523096G0.20090.3120.3141.0369 (0.7466–1.4399)99:401, 80:3360.0470.8287
rs518394G0.51370.3160.3141.0112 (0.7294–1.4018)403:97, 341:830.0050.9465
rs564398C0.19350.2570.2611.4917 (1.0345–2.1511)89:411, 54:3724.6250.0315 *
rs7865618G0.18220.2730.2721.1249 (0.7906–1.6004)85:415, 65:3570.4290.5127
rs10757272T0.20090.4530.4841.1612 (0.8882–1.5180)301:197, 225:1691.0110.3146
rs4977574G0.17450.4730.4861.3515 (1.0462–1.7459)308:192, 233:1934.5150.0336 *
rs2891168G0.06980.5290.5002.1908 (1.6920–2.8368)301:197, 167:25740.6191.85 × 10−10 **
rs1333042A0.96680.7210.4852.2012 (1.6932–2.8616)251:249, 133:29334.1375.14 × 10−9 **
rs2383207A0.08270.2840.3121.0005 (0.7165–1.3971)96:400, 77:32100.9976
rs10757278G0.35870.4450.4991.1438 (0.8824–1.4827)249:251, 196:2261.0310.3099
rs1333048C0.15950.4410.4741.2095 (0.9225–1.5859)317:183, 232:1621.8970.1684
rs1333049C0. 43100.4440.4991.1034 (0.8523–1.4283)246:254, 200:2280.5640.4526
p-HW, p value for Hardy–Weinberg equilibrium analysis; ObsHET: Observed heterozygosity; PredHET: Expected or predicted heterozygosity. * Significant at p < 0.05; ** Significant at p < 0.0001. CI: Confidence interval.
Table
Table 3. Significant haplotypes associated with CAD in Saudi Arabian population.
Table 3. Significant haplotypes associated with CAD in Saudi Arabian population.
BlockHaplotypeFrequencyCase, Control FrequencyChi Squarep Value
Block 1AGTA *0.7320.701, 0.7685.2450.0220
Block 2TGG **0.4650.568, 0.34546.389.74 × 10−12
Block 3GGGCC *0.3210.288, 0.3605.6470.0175
Sig. SNPsTAAG *0.1470.106, 0.19514.4731.00 × 10−4
TGGA **0.1470.214, 0.06739.7692.86 × 10−10
* Protective haplotype; ** Risk haplotype. Sig. SNPs: Block of significant SNPs (Order of Significant SNPs: rs564398, rs4977574, rs2891168, rs1333042).
Table
Table 4. Comparison between two multivariate models with and without CAD associated SNPs. AUC: area under curve.
Table 4. Comparison between two multivariate models with and without CAD associated SNPs. AUC: area under curve.
ModelVariablesAUC (95% CI)p ValueDifference in ACUs (95%CI): Model 2–Model 1p Value
Clinical variables only (Model 1)Age, Gender, BMI0.79 (0.73–0.84)1× 10−100.08 (0.04–0.12)0.000236
Clinical variables + SNPs (Model 2)Age, Gender, BMI, rs1333042, rs28911680.87 (0.82–0.90)1× 10−10

Share and Cite

MDPI and ACS Style

AbdulAzeez, S.; Al-Nafie, A.N.; Al-Shehri, A.; Borgio, J.F.; Baranova, E.V.; Al-Madan, M.S.; Al-Ali, R.A.; Al-Muhanna, F.; Al-Ali, A.; Al-Mansori, M.; et al. Intronic Polymorphisms in the CDKN2B-AS1 Gene Are Strongly Associated with the Risk of Myocardial Infarction and Coronary Artery Disease in the Saudi Population. Int. J. Mol. Sci. 2016, 17, 395. https://doi.org/10.3390/ijms17030395

AMA Style

AbdulAzeez S, Al-Nafie AN, Al-Shehri A, Borgio JF, Baranova EV, Al-Madan MS, Al-Ali RA, Al-Muhanna F, Al-Ali A, Al-Mansori M, et al. Intronic Polymorphisms in the CDKN2B-AS1 Gene Are Strongly Associated with the Risk of Myocardial Infarction and Coronary Artery Disease in the Saudi Population. International Journal of Molecular Sciences. 2016; 17(3):395. https://doi.org/10.3390/ijms17030395

Chicago/Turabian Style

AbdulAzeez, Sayed, Awatif N. Al-Nafie, Abdullah Al-Shehri, J. Francis Borgio, Ekaterina V. Baranova, Mohammed S. Al-Madan, Rudaynah A. Al-Ali, Fahad Al-Muhanna, Abdullah Al-Ali, Mohammed Al-Mansori, and et al. 2016. "Intronic Polymorphisms in the CDKN2B-AS1 Gene Are Strongly Associated with the Risk of Myocardial Infarction and Coronary Artery Disease in the Saudi Population" International Journal of Molecular Sciences 17, no. 3: 395. https://doi.org/10.3390/ijms17030395

APA Style

AbdulAzeez, S., Al-Nafie, A. N., Al-Shehri, A., Borgio, J. F., Baranova, E. V., Al-Madan, M. S., Al-Ali, R. A., Al-Muhanna, F., Al-Ali, A., Al-Mansori, M., Ibrahim, M. F., Asselbergs, F. W., Keating, B., Koeleman, B. P. C., & Al-Ali, A. K. (2016). Intronic Polymorphisms in the CDKN2B-AS1 Gene Are Strongly Associated with the Risk of Myocardial Infarction and Coronary Artery Disease in the Saudi Population. International Journal of Molecular Sciences, 17(3), 395. https://doi.org/10.3390/ijms17030395

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop