Association of six genetic variants with myocardial infarction
- Authors:
- Published online on: February 27, 2015 https://doi.org/10.3892/ijmm.2015.2115
- Pages: 1451-1459
Abstract
Introduction
Myocardial infarction (MI) is a life-threatening disease, the pathogenesis of which is complex and multifactorial. Epidemiological studies have identified several risk factors, including hypertension, diabetes mellitus, dyslipidemia, chronic kidney disease and smoking (1). In addition to these risk factors, recent studies have demonstrated the importance of genetic factors and interactions between multiple genes and environmental factors in the development of MI (2,3). In fact, twin and family studies have indicated that the heritability of coronary heart disease (CHD) is 40–50% (4).
Recent genome-wide association studies (GWAS) have become a powerful approach to identifying genetic variants that are related to common complex diseases (5,6). To date, meta-analyses of GWAS have identified 46 loci associated with CHD or MI at a genome-wide significance level (P<5×10−8) (5,6), and an additional 104 genetic variants have been found to be associated with CHD with a false discovery rate (FDR) of <0.05 among individuals of European and south Asian ancestry (5). Although various loci and genes that confer susceptibility to MI or CHD have been identified for Caucasian populations by GWAS (5,6), genetic variants related to these conditions in Japanese individuals have not been identified definitively.
Therefore, the aim of this study was to examine a possible association of MI in Japanese individuals with 29 single nucleotide polymorphisms (SNPs) identified as susceptibility loci for MI or CHD in Caucasian populations by meta-analyses of GWAS.
Materials and methods
Study population
The study population comprised 4,153 Japanese individuals, including 1,824 subjects with MI (1,468 males and 356 females) and 2,329 controls (1,029 males and 1,300 females), who either visited outpatient clinics or were admitted to the participating hospitals (Gifu Prefectural General Medical Center, Gifu; Gifu Prefectural Tajimi Hospital, Tajimi; Japanese Red Cross Nagoya First Hospital, Nagoya; Inabe General Hospital, Inabe; Hirosaki University Hospital and Hirosaki Stroke Center, Hirosaki, Japan) between 2002 and 2012, due to various symptoms or for an annual health checkup. The 1,824 subjects with MI all underwent coronary angiography and left ventriculography. The diagnosis of MI was based on typical electrocardiographic changes and on increases both in the serum activity of creatine kinase (MB isoform) and in the serum concentration of troponin T. The diagnosis was confirmed by the presence of stenosis in any of the major coronary arteries or in the left main trunk by coronary angiography. The 2,329 control individuals had no history of CHD or MI, aortic aneurysm or peripheral arterial occlusive disease, or of ischemic or hemorrhagic stroke or any other cerebrovascular diseases, or of any other atherosclerotic, thrombotic, embolic or hemorrhagic disorders. The study protocol complied with the Declaration of Helsinki and was approved by the ethics committees of each participating institution. Written informed consent was obtained from each participant.
Selection and genotyping of polymorphisms
We searched SNPs that were shown to be significantly associated with CHD or MI in Caucasian populations by meta-analyses of GWAS (5,6). We examined these SNPs using the SNP database (dbSNP; National Center for Biotechnology Information, Bethesda, MD, USA) to find SNPs with a minor allele frequency of ≥0.015 in a Japanese population. We finally selected 29 SNPs (Table I) and examined a possible association with MI. The wild-type (ancestral allele) and variant alleles of the SNPs were determined from the original sources.
Venous blood (7 ml) was collected into tubes containing 50 mmol/l ethylenediaminetetraacetic acid (disodium salt), the peripheral blood leukocytes were isolated and genomic DNA was extracted from these cells using a DNA extraction kit (Genomix; Talent Srl, Trieste, Italy). Genotypes of the 29 SNPs were determined at G&G Science Co., Ltd. (Fukushima, Japan) by a method that combines the polymerase chain reaction and sequence-specific oligonucleotide probes with suspension array technology (Luminex, Austin, TX, USA). Primers, probes and other conditions for the genotyping of the SNPs examined in this study are shown in Table II. The overall call rate of the genotypes of the 29 SNPs was 99%. Detailed genotyping methodology was as previously described (7).
Statistical analysis
The χ2 test was used to compare categorical variables. Quantitative data were compared between 2 groups using the Mann-Whitney U test or among 3 groups using the Kruskal-Wallis test, given that data were not normally distributed [P<0.01 by the Lilliefors (Kolmogorov-Smirnov) test]. The allele frequencies of each SNP were compared between the subjects with MI and the controls using the χ2 test. To compensate for multiple comparisons of genotypes, a FDR was calculated from the distribution of P-values for the 29 SNPs, and a FDR of <0.05 was considered as statistically significant for association. The statistical power of each SNP was calculated with an Online Sample Size Estimator (http://osse.bii.a-star.edu.sg/calculation2.php). Multivariable logistic regression analysis was performed with MI as a dependent variable and independent variables, including age, gender (0, female; 1, male), body mass index (BMI), smoking status (0, non-smoker; 1, current or former smoker), the prevalence of hypertension, diabetes mellitus and dyslipidemia (0, no history of these conditions; 1, positive history) and the genotype of each SNP; the P-value, odds ratio and 95% confidence interval were calculated. A genotype was assessed according to a dominant (a combined group of heterozygotes and variant homozygotes vs. wild-type homozygotes), recessive (variant homozygotes vs. a combined group of wild-type homozygotes and heterozygotes) and additive model. Additive models included the additive 1 (heterozygotes vs. wild-type homozygotes) and additive 2 (variant homozygotes vs. wild-type homozygotes) models, which were analyzed simultaneously with a single statistical model. Statistical tests were performed using JMP version 11 and JMP Genomics version 6.0 software (SAS Institute Inc., Cary, NC, USA).
Results
Clinical characteristics of the study subjects
The clinical characteristics of the study subjects are presented in Table III. Age, the frequency of the male gender, BMI, the prevalence of smoking, hypertension, diabetes mellitus and dyslipidemia, the fasting plasma glucose levels, the serum concentrations of triglycerides, low-density lipoprotein (LDL)-cholesterol levels and creatinine levels were all higher, whereas diastolic blood pressure, the serum concentrations of high-density lipoprotein (HDL)-cholesterol and the estimated glomerular filtration rate were all lower in the subjects with MI compared to the controls.
Association of SNPs with MI
The allele frequencies of each SNP were compared between the subjects with MI and the controls using the χ2 test, and the SNPs significantly (FDR <0.05) associated with MI are presented in Table IV. Of the 29 SNPs examined in the present study, rs12413409 of the cyclin and CBS domain divalent metal cation transport mediator 2 gene (CNNM2), rs9319428 of the fms-like tyrosine kinase-1 gene (FLT1), rs264 of the lipoprotein lipase gene (LPL), rs599839 of the proline/serine rich coiled coil 1 gene (PSRC1), rs4977574 of the CDKN2B antisense RNA 1 gene (CDKN2B-AS1) and rs9369640 of the phosphatase and actin regulator 1 gene (PHACTR1) were significantly associated with MI. The statistical power of each SNP was calculated with the number of samples and minor allele frequencies of the cases and the controls, and the significance level (α=0.05) was 49.0–85.5%. The genotype distributions of the 6 SNPs were in Hardy-Weinberg equilibrium (P>0.05) among the control individuals.
Table IVThe 6 single nucleotide polymorphisms (SNPs) significantly (allele frequency FDR of <0.05) associated with myocardial infarction (MI). |
Multivariable logistic regression analysis with adjustment for age, gender, BMI, smoking status and the prevalence of hypertension, diabetes mellitus and dyslipidemia revealed that rs12413409 of CNNM2 (dominant, recessive and additive 1 and 2 models), rs264 of LPL (dominant model), rs599839 of PSRC1 (dominant and additive 1 models) and rs9369640 of PHACTR1 (dominant, recessive and additive 1 and 2 models) were significantly associated with MI, with the minor A, A, G and C alleles, respectively, being protective against this condition (Table V). Similar analysis revealed that rs9319428 of FLT1 (dominant and additive 2 models) and rs4977574 of CDKN2B-AS1 (dominant, recessive and additive 2 models) were also significantly associated with MI, with the A and G alleles, respectively, representing risk factors for this condition (Table V).
Table VMultivariable logistic regression analysis of the 6 single nucleotide polymorphisms (SNPs) associated with myocardial infarction. |
Association of the 6 SNPs to conventional risk factors for CHD
We examined the association of the 6 SNPs associated with MI to conventional risk factors for CHD among genotypes or between 2 groups (dominant or recessive model) (Table VI). Due to the multiple comparisons of genotypes, a P-value of <0.001 (0.05/48) was considered statistically significant. rs599839 of PSRC1 was significantly associated with the serum concentrations of HDL-cholesterol with the minor G allele being related to increased serum HDL-cholesterol. This SNP was also associated (borderline significance) with the serum concentrations of LDL-cholesterol with the G allele being related to decreased serum LDL-cholesterol. The rs264 of LPL was associated (borderline significance) with the serum concentrations of HDL-cholesterol with the minor A allele being associated with increased serum HDL-cholesterol.
Table VIComparisons of coronary risk factors among genotypes or between 2 groups (dominant or recessive model) of the SNPs. |
Discussion
In the present study, we examined the association of MI in Japanese individuals with 29 SNPs identified as susceptibility loci for CHD or MI in Caucasian populations by meta-analyses of GWAS and observed that 6 SNPs in the CNNM2, FLT1, LPL, PSRC1, CDKN2B-AS1 and PHACTR1 genes were significantly associated with MI. The association of the SNPs with MI remained significant following adjustment for covariates, suggesting that these genetic variants were independent risk factors for MI.
PHACTR1 is highly expressed in the heart and brain and regulates protein phosphatase 1, an ubiquitous enzyme regulating essential cellular processes (8). The suppression of PHACTR1 in cultured endothelial cells has been shown to result in attenuated tube formation and cell survival (9), as well as in the impaired synthesis of nitric oxide, an important mediator of cardiovascular function (10). In addition, protein phosphatase 1 is considered to be a critical regulator of Ca2+ cycling and contractility in smooth muscle cells, and thus, is closely related to arterial elasticity and wave reflection which is recognized as an important factor for cardiovascular disease (11). In the present study, rs9369640 of PHACTR1 showed the most significant association with MI with the minor C allele being protective against this condition, although the role of PHACTR1 in the pathogenesis of MI remains unclear.
CDKN2B-AS1 is located within the CDKN2B-CDKN2A gene cluster at chromosome 9p21. The 9p21.3 locus was first discovered by GWAS to be a risk factor for CHD or MI (12–15). This locus was then shown to be associated with ischemic stroke (16,17) and with intracranial and aortic aneurysms (18). The 9p21.3 locus has also been associated with platelet reactivity and the increased platelet reactivity may explain the association with MI and ischemic stroke (19). The 9p21.3 locus includes CDKN2B-AS1, which encodes a functional RNA molecule that interacts with polycomb repressive complex-1 and -2, leading to epigenetic silencing of other genes in this cluster (NCBI Gene) and to the alteration of the expression of several genes related to cellular proliferation (20). Recent studies have also indicated that polymorphisms at 9p21 influence inflammatory signaling (21) and vascular cell proliferation (22). In this study, we demonstrated that rs4977574 of CDKN2B-AS1 was significantly associated with MI in Japanese individuals with the G allele representing a risk factor for this condition. These previous data and the data from the present study indicate that CDKN2B-AS1 at 9p21.3 may be a susceptibility locus for MI, although the underlying molecular mechanisms have not been determined definitively.
The lipoprotein lipase encoded by LPL hydrolyzes triglycerides in circulating lipoprotein into 2 free fatty acids and 1 monoacylglycerol molecules. It also binds simultaneously to both lipoproteins and cell surface receptors and promotes the cellular uptake of lipoproteins and free fatty acids. Genetic variants causing decreased activity of LPL have been shown to result in disorders of lipoprotein metabolism, such as hypertriglyceridemia and hypo-HDL-cholesterolemia (23). In addition, genetic variants of LPL have been implicated in insulin resistance and type 2 diabetes mellitus (24). In the present study, rs264 of LPL was found to be associated with the serum concentrations of HDL-cholesterol with the A allele being related to increased HDL-cholesterol levels. We observed that rs264 was significantly associated with MI with the A allele being protective against MI. This association may be attributable, at least in part, to the effect of rs264 on the metabolism of HDL-cholesterol.
The rs599839 of PSRC1 is located in the cadherin, EGF LAG seven-pass G-type receptor 2 (CELSR2)-PSRC1-sortilin 1 (SORT1) gene cluster on chromosome 1p13.3. The SORT1 protein, the higher expression of which is associated with the minor G allele of rs599839, is a multiligand transmembrane receptor protein that binds to a variety of ligands, including the LDL-receptor associated protein (25), lipoprotein lipase (26) and apolipoprotein A-V (27), and enhances the endocytosis and intracellular degradation of LDL-cholesterol. In the present study, the minor G allele of rs599839 was found to be associated with an increase in serum HDL-cholesterol levels and a decrease in serum LDL-cholesterol levels. We demonstrated that rs599839 of PSRC1 was associated with MI in Japanese individuals with the G allele being protective against MI. The effects of rs599839 on the metabolism of HDL- and LDL-cholesterol may account for its association with MI.
rs12413409 is located in CNNM2 at 10q24.32. CNNM2 was identified as a gene that plays an important role in magnesium homeostasis by mediating the epithelial transport and renal reabsorption of Mg2+ (28,29). We observed the association of rs12413409 of CNNM2 with MI in a Japanese population with the minor A allele being protective, although the role of CNNM2 in the pathogenesis of CHD or MI remains unclear.
FLT1 encodes vascular endothelial growth factor (VEGF) receptor 1 in humans and plays an important role in the pathological mechanisms of atherosclerosis. In addition to VEGF, placental growth factor selectively binds to VEGF receptor 1 and accelerates atherosclerotic processes through the enhancement of intramural angiogenesis and monocyte recruitment (30–32). Recruited monocytes, differentiating into macrophages, produce pro-inflammatory cytokines, such as monocyte chemoattractant protein-1 which cause the augmentation of the inflammatory response and plaque vulnerability (33). In this study, we demonstrated that rs9319428 of FLT1 was significantly associated with MI with the A allele representing a risk factor for MI. The effects of rs9319428 on the development of atherosclerosis may account for its association with MI.
In conclusion, the results from the present study suggest that rs12413409 of CNNM2, rs264 of LPL, rs9369640 of PHACTR1, rs9319428 of FLT1, rs599839 of PSRC1 and rs4977574 of CDKN2B-AS1 are susceptibility loci for MI in Japanese individuals. The determination of genotypes for these SNPs may prove informative for the assessment of the genetic risk for MI. Given that multiple variants, each having a small effect, will ultimately be found to be responsible for a large fraction of the genetic component of MI, further identification of MI susceptibility genes will allow for a more accurate assessment of the genetic component of this condition.
Acknowledgments
The present study was supported by a Collaborative Research Grant from the Gifu Prefectural General Medical Center (H24-26 to Y.Y.) and a Grant-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and Technology of Japan (no. 24590746 to Y.Y.).
References
Yusuf S, Hawken S, Ounpuu S, et al: Effect of potentially modifiable risk factors associated with myocardial infarction In 52 countries (the INTERHEART study): case-control study. Lancet. 364:937–952. 2004. View Article : Google Scholar : PubMed/NCBI | |
Topol EJ, Smith J, Plow EF and Wang QK: Genetic susceptibility to myocardial infarction and coronary artery disease. Hum Mol Genet. 15:R117–R123. 2006. View Article : Google Scholar : PubMed/NCBI | |
Wang Q: Molecular genetics of coronary artery disease. Curr Opin Cardiol. 20:182–188. 2005. View Article : Google Scholar : PubMed/NCBI | |
Peden JF and Farrall M: Thirty-five common variants for coronary artery disease: the fruits of much collaborative labour. Hum Mol Genet. 20:R198–R205. 2011. View Article : Google Scholar : PubMed/NCBI | |
CARDIoGRAMplusC4D Consortium: Large-scale association analysis identifies new risk loci for coronary artery disease. Nat Genet. 45:25–33. 2013. | |
Schunkert H, König IR, Kathiresan S, et al: Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease. Nat Genet. 43:333–338. 2011. View Article : Google Scholar : PubMed/NCBI | |
Itoh Y, Mizuki N, Shimada T, et al: High-throughput DNA typing of HLA-A, -B, -C, and -DRB1 loci by a PCR-SSOP-Luminex method in the Japanese population. Immunogenetics. 57:717–729. 2005. View Article : Google Scholar : PubMed/NCBI | |
Peli W, Nairn AC and Page R: Structural basis for protein phosphatase 1 regulation and specificity. FEBS J. 280:596–611. 2013. View Article : Google Scholar | |
Jarray R, Allain B, Borriello L, et al: Depletion of the novel protein PHACTR-1 from human endothelial cells abolishes tube formation and induces cell death receptor apoptosis. Biochimie. 93:1668–1675. 2011. View Article : Google Scholar : PubMed/NCBI | |
Mount PF, Kemp BE and Power DA: Regulation of endothelial and myocardial NO synthesis by multi-site eNOS phosphorylation. J Mol Cell Cardiol. 42:271–279. 2007. View Article : Google Scholar | |
Patel RS, Morris AA, Ahmed Y, et al: A genetic risk variant for myocardial infarction on chromosome 6p24 is associated with impaired central hemodynamic indexes. Am J Hypertens. 25:797–803. 2012. View Article : Google Scholar : PubMed/NCBI | |
Helgadottir A, Thorleifsson G, Manolescu A, et al: A common variant on chromosome 9p21 affects the risk of myocardial infarction. Science. 316:1491–1493. 2007. View Article : Google Scholar : PubMed/NCBI | |
McPherson R, Pertsemlidis A, Kavaslar N, et al: A common allele on chromosome 9 associated with coronary heart disease. Science. 316:1488–1491. 2007. View Article : Google Scholar : PubMed/NCBI | |
Samani NJ, Erdmann J, Hall AS, et al: Genomewide association analysis of coronary artery disease. N Engl J Med. 357:443–453. 2007. View Article : Google Scholar : PubMed/NCBI | |
Wellcome Trust Case Control Consortium: Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 447:661–678. 2007. View Article : Google Scholar : PubMed/NCBI | |
Matarin M, Brown WM, Singleton A and Hardy JA: Whole genome analyses suggest ischemic stroke and heart disease share an association with polymorphisms on chromosome 9p21. Stroke. 39:1586–1589. 2008. View Article : Google Scholar : PubMed/NCBI | |
Gschwendtner A, Bevan S, Cole JW, et al: Sequence variants on chromosome 9p21.3 confer risk for atherosclerotic stroke. Ann Neurol. 65:531–539. 2009. View Article : Google Scholar : PubMed/NCBI | |
Helgadottir A, Thorleifsson G, Magnusson KP, et al: The same sequence variant on 9p21 associates with myocardial infarction, abdominal aortic aneurysm and intracranial aneurysm. Nat Genet. 40:217–224. 2008. View Article : Google Scholar : PubMed/NCBI | |
Musunuru K, Post WS, Herzog W, et al: Association of single nucleotide polymorphisms on chromosome 9p21.3 with platelet reactivity: a potential mechanism for increased vascular disease. Circ Cardiovasc Genet. 3:445–453. 2010. View Article : Google Scholar : PubMed/NCBI | |
Jarinova O, Stewart AF, Roberts R, et al: Functional analysis of the chromosome 9p21.3 coronary artery disease risk locus. Arterioscler Thromb Vasc Biol. 29:1671–1677. 2009. View Article : Google Scholar : PubMed/NCBI | |
Harismendy O, Notani D, Song X, et al: 9p21 DNA variants associated with coronary artery disease impair interferon-γ signalling response. Nature. 470:264–268. 2011. View Article : Google Scholar : PubMed/NCBI | |
Visel A, Zhu Y, May D, et al: Targeted deletion of the 9p21 non-coding coronary artery disease risk interval in mice. Nature. 464:409–412. 2010. View Article : Google Scholar : PubMed/NCBI | |
Fisher RM, Humphries SE and Talmud PJ: Common variation in the lipoprotein lipase gene: effect on plasma lipids and risk of atherosclerosis. Atherosclerosis. 135:145–159. 1997. View Article : Google Scholar | |
Eriksson JW, Buren J, Svensson M, Olivecrona T and Olivecrona G: Postprandial regulation of blood lipids and adipose tissue lipoprotein lipase in type 2 diabetes patients and healthy control subjects. Atherosclerosis. 166:359–367. 2003. View Article : Google Scholar : PubMed/NCBI | |
Jansen P, Giehl K, Nyengaard JR, et al: Roles for the proneurotrophin receptor sortilin in neuronal development, aging and brain injury. Nat Neurosci. 10:1449–1457. 2007. View Article : Google Scholar : PubMed/NCBI | |
Nielsen MS, Jacobsen C, Olivecrona G, Gliemann J and Petersen CM: Sortilin/neurotensin receptor-3 binds and mediates degradation of lipoprotein lipase. J Biol Chem. 274:8832–8836. 1999. View Article : Google Scholar : PubMed/NCBI | |
Nilsson SK, Christensen S, Raarup MK, Ryan RO, Nielsen MS and Olivecrona G: Endocytosis of apolipoprotein A-V by members of the low density lipoprotein receptor and the VPS10p domain receptor families. J Biol Chem. 283:25920–25927. 2008. View Article : Google Scholar : PubMed/NCBI | |
Stuiver M, Lainez S, Will C, et al: CNNM2, encoding a baso-lateral protein required for renal Mg2+ handling, is mutated in dominant hypomagnesemia. Am J Hum Genet. 88:333–343. 2011. View Article : Google Scholar : PubMed/NCBI | |
De Baaij JHF, Stuiver M, Meij IC, et al: Membrane topology and intracellular processing of cyclin M2 (CNNM2). J Biol Chem. 287:13644–13655. 2012. View Article : Google Scholar : PubMed/NCBI | |
Khurana R, Moons L, Shafi S, et al: Placental growth factor promotes atherosclerotic intimal thickening and macrophage accumulation. Circulation. 111:2828–2836. 2005. View Article : Google Scholar : PubMed/NCBI | |
Pilarczyk K, Sattler KJ, Galili O, et al: Placenta growth factor expression in human atherosclerotic carotid plaques is related to plaque destabilization. Atherosclerosis. 196:333–340. 2008. View Article : Google Scholar | |
Roncal C, Buysschaert I, Gerdes N, et al: Short-term delivery of anti-PlGF antibody delays progression of atherosclerotic plaques to vulnerable lesions. Cardiovasc Res. 86:29–36. 2010. View Article : Google Scholar | |
Carmeliet P, Moons L, Luttun A, et al: Synergism between vascular endothelial growth factor and placental growth factor contributes to angiogenesis and plasma extravasation in pathological conditions. Nat Med. 7:575–583. 2001. View Article : Google Scholar : PubMed/NCBI |