Skip to main content

Advertisement

Genetic variants of glutamate receptor gene family in Taiwanese Kawasaki disease children with coronary artery aneurysms

Abstract

Background

Patients with Kawasaki disease (KD), a pediatric systemic vasculitis, may develop coronary artery aneurysm (CAA) as a complication. To investigate the role of glutamate receptors in KD and its CAA development, we performed genetic association studies.

Methods and results

We examined the whole family of glutamate receptors by genetic association studies in a Taiwanese cohort of 262 KD patients. We identified glutamate receptor ionotropic, kainate 1 (GRIK1) as a novel susceptibility locus associated with CAA formation in KD. Statistically significant differences were noted for factors like fever duration, 1st Intravenous immunoglobulin (IVIG) used time (number of days after the first day of fever) and the GRIK1 (rs466013, rs425507, and rs38700) genetic variants. This significant association persisted even after using multivariate regression analysis (Full model: for rs466013: odds ratio =2.12; 95% CI =1.22-3.65; for rs425507: odds ratio =2.16; 95% CI =1.26-3.76; for rs388700: odds ratio =2.16; 95% CI =1.26-3.76).

Conclusions

We demonstrated that GRIK1 polymorphisms are associated CAA formation in KD, even when adjusted for fever duration and IVIG used time, and may also serve as a genetic marker for the CAA formation in KD.

Background

Patients with Kawasaki disease (KD), an acute systemic vasculitis, may develop coronary artery aneurysm (CAA) as a complication. KD is one of the leading causes of acquired cardiovascular diseases in childhood. Infectious agents, host immune dysregulation, and genetic susceptibility are thought to be responsible for the development of KD and its related complications[13]. However, the pathological mechanisms underlying KD remain to be elucidated.

Numerous genome-wide association studies have been conducted to identify host cellular genes that affect KD susceptibility[414] in the European, Japanese, Korean, and Taiwanese populations. In the European population[11, 13], no common SNPs have been identified as susceptibility loci for European KD. However, a common SNP (rs2233152; MIA gene) was observed in the European, Japanese, and Taiwanese populations[911]. Common gene SNPs among Asians including Japanese, Taiwanese, and Korean populations have also been observed[4, 6, 9, 10, 12, 14, 15]. Six SNPs, namely, rs2736340 (BLK), rs2618479 (BLK), rs6993775 (BLK), rs10401344 (ITPKC), rs2233152 (MIA), and rs4813003 (CD40) have been observed in both Japanese and Taiwanese populations[9, 10] (Additional files1 and2). These studies suggest that genes involved in the immune-regulatory responses and cardiovascular-related pathogenesis may contribute to KD susceptibility.

Glutamate receptors were initially demonstrated to play important roles in excitatory neurotransmission in the brain and interneuronal communication[16]. Based on their different activation mechanisms, glutamate receptors can be divided into 2 groups: ionotropic glutamate receptors (iGluRs) and metabotropic glutamate receptors (mGluRs). The human genome is known to contain at least 16 iGluRs and 8 mGluRs. Based on their agonist binding and electrophysiological properties, iGluRs can be classified to 3 groups: alpha-amino-3-hydroxy-5-methyl-4-isoxazole (AMPA), N-methyl-D-aspartate (NMDA), and kainate (KA) receptors. Genetic mutations in glutamate receptors are associated with a number of human diseases including autism, Huntington’s disease and Parkinson’s disease[17, 18]. In addition, glutamate receptors have been found to influence autoantigen/antibody interactions and multiple sclerosis. GluR3 (GRIA3) is known to act as an autoantigen in Rasmussen’s encephalitis, suggesting a strong link between glutamate receptors and autoimmune interaction in certain degenerative diseases[19]. The regulation of glutamate receptor binding activity can reduce central nervous system (CNS) inflammation, apoptosis, and axonal damage[20]. In addition, glutamate receptors have also been implicated in cardiovascular diseases[21]. Glutamate receptor 1 (GluR1), an AMPA receptor subtype, can mediate the regulation of platelet activation through glutamate and GluR1 knockout mice develop in vivo thrombosis after a prolonged time[22]. The activation of GluR1 may contribute to the development of cardiovascular disease via accelerating thrombus formation.

Endothelial cells are principal targets for ischemic free-radical injury. Glutamate receptors are known to prevent nitric oxide-induced vascular injury[23]. On the other hand, activation of certain glutamate receptors was demonstrated to be a potential strategy for disrupting angiogenesis[24]. Coronary artery damage in KD is strongly associated with endothelial cell dysfunction[25]. Additional evidence suggests that glutamate receptors may influence KD pathogenesis[26, 27].

To explore the role of glutamate receptors in KD development, we investigated the entire family of glutamate receptors by performing genetic association studies on a Taiwanese cohort of 262 KD patients. Our study identified glutamate receptor ionotropic, kainate 1 (GRIK1) as a novel susceptibility locus on 21q21.3. To our knowledge, this is the first instance to screen the glutamate receptor family for the association between genetic variants of glutamate receptors and CAA formation in KD.

Results

Genetic association study of the glutamate receptor gene family in Taiwanese KD children and controls

To identify KD susceptibility genes, a total of 53 SNPs of 16 genes within the glutamate receptor gene family including GRIK1, GRIK2, GRIK3, GRIK4, GRIK5, GRIA1, GRIA2, GRIA4, GRM1, GRM2, GRM3, GRM4, GRM5, GRM6, GRM7, and GRM8 genes were genotyped in 262 Taiwanese KD children and in 1107 healthy people from the general population of Taiwan who were Han Chinese ethnic background for the SNP association study (Table 1). No significant differences were found between these 2 groups, suggesting that the glutamate receptor family genes may not contribute to KD susceptibility.

Table 1 Genotype distribution of glutamate receptor family gene SNPs in Taiwanese KD patients and controls

GRIK1 genetic polymorphisms may be related to KD-associated CAA complications

To examine the role of glutamate receptors in KD-associated CAA complications, we analyzed the correlation between KD children and the whole glutamate gene family. As shown in Table 2, the genotype distributions (dominant model) of 6 glutamate gene SNPs were statistically different between these 2 groups (p <0.05). These SNPs were rs466013, rs425507, rs388700, rs402280, rs17104835 and rs712723. Among these, 4 SNPs were found to be located in the GRIK1 gene (p =0.007, 0.005, 0.004 and 0.022, respectively) (Additional file1). GRIK1 consists of 18 exons and is located at 21q21.3 as shown in Figure 1. All SNPs were in Hardy-Weinberg equilibrium and had a successful genotyping frequency of >99%. The linkage disequilibrium (LD) structure of this region was also established, with 1 haplotype block determined. Four SNPs were located in that block. To evaluate the relationship among these 4 SNPs, pairwise LD analysis was performed. The D’ statistics were all 1.0. Strong LD was observed in the following 2 groups of SNPs, group1 (rs466013, rs425507, rs388700), with the r2 statistics >0.5 between every 2 SNPs in each group (data not shown). The frequencies of the TT and TC genotypes of GRIK1 (rs466013) were significantly higher in KD patients with CAA than those in patients without CAA (63.2% for KD with CAA and 44.9% for KD without CAA complications; odds ratio =2.11 [95% confidence interval (CI) =1.22-3.65]). Similar results were also observed in rs425507, rs388700 and rs402280. These data suggest that GRIK1 may be a potential susceptibility locus involved in the development of KD with CAA complications.

Table 2 Association of the genetic variants of glutamate receptor family genes in Taiwanese KD children according to the presence or absence of CAA
Figure 1
figure1

Analysis of single nucleotide polymorphisms (SNPs) and the linkage disequilibrium (LD) pattern of the GRIK1 gene. Genomic location of SNPs present on chromosome 21q21.3. Physical position of individual SNPs was based on the NCBI Assembly database: GRCh38 version. Linkage disequilibrium (LD) blocks in the GRIK1 gene, estimated by using HAPLOVIEW software. Pairwise D’ values (%) are indicated in squares; red indicates linkage disequilibrium [D’ =1, logarithm of odds (LOD) ≥2].

Multivariate regression analyses shows that GRIK1 genetic polymorphisms may be related to CAA formation in KD

According to the above results, statistically significant differences in factors associated with CAA formation in KD were noted for the clinical characteristics including fever duration (p <0.0001), first IVIG used time (p <0.0001; number of days after the first day of fever), and the GRIK1 (rs466013, rs425507, rs38700, and rs402280) genetic variants (p =0.007, p =0.005, p =0.004, and p =0.022, respectively) (Tables 1 and3). To further confirm the genetic role of GRIK1, we used multivariate regression analyses to adjust those potential factors (i.e., fever duration and IVIG used time) that may affect the analysis. As shown in Table 3, significant associations between KD with CAA complications and the GRIK1 (rs466013, rs425507, rs38700 and rs402280) genetic variants were observed (Full model: for rs466013: odds ratio = 2.12; 95% CI = 1.22-3.65; for rs425507: odds ratio = 2.16; 95% CI = 1.26-3.76; for rs388700: odds ratio = 2.16; 95% CI = 1.26-3.76; for rs402280: odds ratio = 1.89; 95% CI = 1.09-3.21). Taken together, these data suggest that the significant association observed between CAA complications and the presence of the GRIK1 genotypes persists even after adjusting for the potential factors.

Table 3 Association of GRIK1 genetic polymorphisms with CAA complications in Taiwanese KD children by multivariate regression analysis

Discussion

Previous research from our lab suggests that the NMDA receptor (GRIN3A) from the glutamate receptor family may influence KD pathogenesis[26]. In this study, we screened the entire glutamate receptor family including the iGluRs and mGluRs (GRIK, GRIA and GRM gene families) and identified another member, namely GRIK1, that may be involved in the development of KD-associated CAA complications in Taiwanese children of Han Chinese ethnic background. The most striking finding of this study is that 4 GRIK1 gene variants were found to be strongly associated with the presence of CAA in KD patients, even in the multivariable model.

Our genetic association study showed that none of the genes of the glutamate receptor gene family including GRIK1, GRIK2, GRIK3, GRIK4, GRIK5, GRIA1, GRIA2, GRIA4, GRM1, GRM2, GRM3, GRM4, GRM5, GRM6, GRM7, and GRM8 genes contributed to KD susceptibility. However, genetic variation of the GRIK1 locus may potential induce susceptibility to the development of KD with CAA complications. The significant association observed between KD with CAA complications and the GRIK1 genetic variants (rs466013, rs425507, rs38700, and rs402280) was found to persist even after adjusting for fever duration and first IVIG used time. These results suggest that the GRIK1 gene may be involved in CAA formation of KD. GRIK1 polymorphisms have been investigation for their associations with different diseases including Juvenile absence epilepsy[28, 29], schizophrenia[30, 31], alcohol dependence[32], topiramate’s effects on heavy drinking[33, 34], topiramate-induced side effects[35], and hepatitis B virus (HBV)-related hepatocellular carcinoma[36]. However, these GRIK1 polymorphism data of various studies are also not absolutely consistent and conclusive. These studies show that GRIK1 gene may mainly contribute to neuropsychological diseases. Glutamate is known to signal and is released by nerves, macrophages, lymphocytes, and chondrocytes[37, 38]. These amino acids bind to iGluRs and mGluRs to regulate peripheral pain, release of cytokines and matrix metalloproteinases, and immune responses[3941]. Our studies have firstly showed that glutamate receptors including NMDA[26] and KA receptors are involved in the CAA complications of KD regardless of the fever duration and first IVIG used time. KD is a multi-systemic disorder with a possible underlying pathology of immune-mediated vasculitis[1, 42]. The vascular endothelium forms a functional barrier between the vessel wall and the bloodstream. Recent studies have shown that regulation of certain glutamate receptors may induce the inflammation of endothelial cells, thereby mediating pathogenesis of vascular diseases[43, 44]. Although the current therapy for KD includes high doses of aspirin in conjunction with IVIG treatment[45], reports suggest that this regimen cannot efficiently prevent CAA development.

In this study, we showed that the glutamate receptor GRIK1 is significantly associated with KD with CAA complications in Taiwanese children with Han Chinese ethnic background. Genetic polymorphisms of the GRIK1 gene may play a role in KD pathogenesis and this molecule may serve as a therapeutic target for the KD treatment to prevent CAA development. Children with specific glutamate receptor genotypes related to KD should be careful assessed for CAA status at the time of diagnosis and monitored during the CAA development and related cardiovascular diseases. In addition to aspirin and IVIG therapy, the glutamate receptors may also serve as good targets for the design of novel KD therapeutics.

Methods

Patients

We performed a retrospective study. Individuals fulfilling the diagnostic criteria of KD (n = 262) were identified and enrolled into this study from the Department of Pediatrics at China Medical University Hospital in Taichung, Taiwan[46, 47]. This study population has been previously used for SNP analysis and KD studies[4, 10, 26, 48, 49]. In this study, there were 164 males and 94 females with an average age at diagnosis 1.75 ± 1.61 years. All the patients were diagnosed according to KD criteria[50]. All the patients underwent regular echocardiography examinations at the acute stage, 2 and 6 months after onset and once a year thereafter. CAA was identified when either the right or left coronary artery showed a dilated diameter of 3 mm in children younger than 5 years of age, or 4 mm in the older children[51]. According to the presence or absence of CAA, statistically significant differences between these 2 groups were found with respect to the fever duration and 1st IVIG used time (number of days after the first day of fever)[26, 48]. Only Han Chinese individuals, who account for 98% of the Taiwanese residents, were considered for recruitment. This study was approved by the Human Studies Committee of China Medical University Hospital (CMUH REC No. DMR101-IRB1-313 (CR-1)).

Consent

The written informed consent was obtained from the patient’s guardian/parent/next of kin for the publication of this report and any accompanying images.

SNP genotyping

Fifty-three single nucleotide polymorphisms (SNPs) of 16 genes within the glutamate receptor gene family including GRIK1, GRIK2, GRIK3, GRIK4, GRIK5, GRIA1, GRIA2, GRIA4, GRM1, GRM2, GRM3, GRM4, GRM5, GRM6, GRM7, and GRM8 were selected from the NCBI SNP database and HAPMAP website (Tables 2 and3)[52]. Selection criteria for including SNPs in the analysis were a minimum allele frequency of >0.05 in the Han Chinese population and a Hardy-Weinberg equilibrium (p >0.05). A summary of information on the SNPs in the glutamate receptor genes (location, position, rs number, and genotype) is presented in Table 1. Briefly, genomic DNA was extracted from peripheral blood leukocytes according to the standard protocols (Genomic DNA kit; Qiagen). SNPs were genotyped using a custom-designed VeraCode GoldenGate Genotyping Assay System (Illumina); genotyping was performed as outlined inhttp://www.illumina.com/.

Primers and probes were designed and created using Custom VeraCode GoldenGate Genotyping Assay System software. Genotype calls were automatically generated using GenCall software version 3.1.3. We assessed the 8 VeraCode runs individually for intra-plate inconsistencies (e.g., variation in fluorescence intensities). Genotype cluster plots generated by individual VeraCode and SAM assays were visually inspected for call quality. Plots that appeared to be "unusually" clustered (i.e., unlike the predicted spread in terms of software-generated HWE or distance between clusters [θ]) were investigated further by selecting samples via direct Sanger sequencing for genotype confirmation. Samples were sequenced using Big Dye Terminator v3.1 (AB, Foster City, CA, USA) according to the manufacturer’s guidelines, and sequenced with an AB 3730 genetic analyzer.

Analysis of haplotype blocks

Based on the HAPLOVIEW software, we used the Lewontin D’ measure to estimate the intermarker coefficient of LD of patients[53]. The confidence interval (CI) of LD was estimated using a resampling procedure and then used to construct the haplotype blocks.

Statistical analyses

Data are expressed as means ± standard deviation for continuous variables. Genotypes were obtained by direct count, followed by allele frequency calculations (Table 1). χ2 tests were performed to determine the differences in categorical variables, and the odds ratio and 95% CI were calculated for the factors under consideration. Forward stepwise multivariate regression analyses were also performed to identify factors that contribute independently to CAA formation in KD. All statistical analyses were performed using SPSS (v12.0) for Windows.

References

  1. 1.

    Burns JC, Glode MP: Kawasaki syndrome. Lancet. 2004, 364: 533-544. 10.1016/S0140-6736(04)16814-1

  2. 2.

    Chang LY, Chang IS, Lu CY, Chiang BL, Lee CY, Chen PJ, Wang JT, Ho HN, Chen DS, Huang LM: Epidemiologic features of Kawasaki disease in Taiwan, 1996–2002. Pediatrics. 2004, 114: e678-682. 10.1542/peds.2004-0726

  3. 3.

    Lin W, Liu HP, Chang JS, Lin YJ: Genetic variations within the PSORS1 region affect Kawasaki disease development and coronary artery aneurysm formation. Biomed. 2013, 3: 73-81. 10.1016/j.biomed.2013.03.001. 10.1016/j.biomed.2013.03.001

  4. 4.

    Chang CJ, Kuo HC, Chang JS, Lee JK, Tsai FJ, Khor CC, Chang LC, Chen SP, Ko TM, Liu YM, Chen YJ, Hong YM, Jang GY, Hibberd ML, Kuijpers T, Burgner D, Levin M, Burns JC, Davila S, Chen YT, Chen CH, Wu JY, Lee YC: Replication and meta-analysis of GWAS identified susceptibility loci in Kawasaki disease confirm the importance of B lymphoid tyrosine kinase (BLK) in disease susceptibility. PLoS One. 2013, 8: e72037. 10.1371/journal.pone.0072037

  5. 5.

    Kim JJ, Park YM, Yoon D, Lee KY, Seob Song M, Doo Lee H, Kim KJ, Park IS, Nam HK, Weon Yun S, Ki Han M, Mi Hong Y, Young Jang G, Lee JK: Identification of KCNN2 as a susceptibility locus for coronary artery aneurysms in Kawasaki disease using genome-wide association analysis. J Hum Genet. 2013, 58: 521-525. 10.1038/jhg.2013.43

  6. 6.

    Yan Y, Ma Y, Liu Y, Hu H, Shen Y, Zhang S, Tao D, Wu Q, Peng Q, Yang Y: Combined analysis of genome-wide-linked susceptibility loci to Kawasaki disease in Han Chinese. Hum Genet. 2013, 132: 669-680. 10.1007/s00439-013-1279-2

  7. 7.

    Lin MT, Hsu CL, Chen PL, Yang WS, Wang JK, Fann CS, Wu MH: A genome-wide association analysis identifies novel susceptibility loci for coronary arterial lesions in patients with Kawasaki disease. Transl Res. 2013, 161: 513-515. 10.1016/j.trsl.2013.02.002

  8. 8.

    Onouchi Y: Genetics of Kawasaki disease: what we know and don’t know. Circ J. 2012, 76: 1581-1586. 10.1253/circj.CJ-12-0568

  9. 9.

    Onouchi Y, Ozaki K, Burns JC, Shimizu C, Terai M, Hamada H, Honda T, Suzuki H, Suenaga T, Takeuchi T, Yoshikawa N, Suzuki Y, Yasukawa K, Ebata R, Higashi K, Saji T, Kemmotsu Y, Takatsuki S, Ouchi K, Kishi F, Yoshikawa T, Nagai T, Hamamoto K, Sato Y, Honda A, Kobayashi H, Sato J, Shibuta S, Miyawaki M, Oishi K: A genome-wide association study identifies three new risk loci for Kawasaki disease. Nat Genet. 2012, 44: 517-521. 10.1038/ng.2220

  10. 10.

    Lee YC, Kuo HC, Chang JS, Chang LY, Huang LM, Chen MR, Liang CD, Chi H, Huang FY, Lee ML, Huang YC, Hwang B, Chiu NC, Hwang KP, Lee PC, Chang LC, Liu YM, Chen YJ, Chen CH, Chen YT, Tsai FJ, Wu JY: Two new susceptibility loci for Kawasaki disease identified through genome-wide association analysis. Nat Genet. 2012, 44: 522-525. 10.1038/ng.2227

  11. 11.

    Khor CC, Davila S, Breunis WB, Lee YC, Shimizu C, Wright VJ, Yeung RS, Tan DE, Sim KS, Wang JJ, Wong TY, Pang J, Mitchell P, Cimaz R, Dahdah N, Cheung YF, Huang GY, Yang W, Park IS, Lee JK, Wu JY, Levin M, Burns JC, Burgner D, Kuijpers TW, Hibberd ML: Genome-wide association study identifies FCGR2A as a susceptibility locus for Kawasaki disease. Nat Genet. 2011, 43: 1241-1246. 10.1038/ng.981

  12. 12.

    Kim JJ, Hong YM, Sohn S, Jang GY, Ha KS, Yun SW, Han MK, Lee KY, Song MS, Lee HD, Kim DS, Lee JE, Shin ES, Jang JH, Lee YS, Kim SY, Lee JY, Han BG, Wu JY, Kim KJ, Park YM, Seo EJ, Park IS, Lee JK: A genome-wide association analysis reveals 1p31 and 2p13.3 as susceptibility loci for Kawasaki disease. Hum Genet. 2011, 129: 487-495. 10.1007/s00439-010-0937-x

  13. 13.

    Burgner D, Davila S, Breunis WB, Ng SB, Li Y, Bonnard C, Ling L, Wright VJ, Thalamuthu A, Odam M, Shimizu C, Burns JC, Levin M, Kuijpers TW, Hibberd ML: A genome-wide association study identifies novel and functionally related susceptibility Loci for Kawasaki disease. PLoS Genet. 2009, 5: e1000319. 10.1371/journal.pgen.1000319

  14. 14.

    Tsai FJ, Lee YC, Chang JS, Huang LM, Huang FY, Chiu NC, Chen MR, Chi H, Lee YJ, Chang LC, Liu YM, Wang HH, Chen CH, Chen YT, Wu JY: Identification of novel susceptibility Loci for kawasaki disease in a Han chinese population by a genome-wide association study. PLoS One. 2011, 6: e16853. 10.1371/journal.pone.0016853

  15. 15.

    Peng Q, Chen C, Zhang Y, He H, Wu Q, Liao J, Li B, Luo C, Hu X, Zheng Z, Yang Y: Single-nucleotide polymorphism rs2290692 in the 3′UTR of ITPKC associated with susceptibility to Kawasaki disease in a Han Chinese population. Pediatr Cardiol. 2012, 33: 1046-1053. 10.1007/s00246-012-0223-x

  16. 16.

    Debanne D, Daoudal G, Sourdet V, Russier M: Brain plasticity and ion channels. J Physiol Paris. 2003, 97: 403-414. 10.1016/j.jphysparis.2004.01.004

  17. 17.

    Diguet E, Fernagut PO, Normand E, Centelles L, Mulle C, Tison F: Experimental basis for the putative role of GluR6/kainate glutamate receptor subunit in Huntington’s disease natural history. Neurobiol Dis. 2004, 15: 667-675. 10.1016/j.nbd.2003.12.010

  18. 18.

    Meldrum B: Amino acids as dietary excitotoxins: a contribution to understanding neurodegenerative disorders. Brain Res Brain Res Rev. 1993, 18: 293-314. 10.1016/0165-0173(93)90014-Q

  19. 19.

    Rogers SW, Andrews PI, Gahring LC, Whisenand T, Cauley K, Crain B, Hughes TE, Heinemann SF, McNamara JO: Autoantibodies to glutamate receptor GluR3 in Rasmussen’s encephalitis. Science. 1994, 265: 648-651. 10.1126/science.8036512

  20. 20.

    Bolton C, Paul C: Glutamate receptors in neuroinflammatory demyelinating disease. Mediators Inflamm. 2006, 2006: 93684.

  21. 21.

    Chen H: Possible Role of Platelet GluR1 Receptors in Comorbid Depression and Cardiovascular Disease. Cardiovasc Psychiatry Neurol. 2009, 2009: 424728.

  22. 22.

    Morrell CN, Sun H, Ikeda M, Beique JC, Swaim AM, Mason E, Martin TV, Thompson LE, Gozen O, Ampagoomian D, Sprengel R, Rothstein J, Faraday N, Huganir R, Lowenstein CJ: Glutamate mediates platelet activation through the AMPA receptor. J Exp Med. 2008, 205: 575-584. 10.1084/jem.20071474

  23. 23.

    Lin SH, Maiese K: Group I metabotropic glutamate receptors prevent endothelial programmed cell death independent from MAP kinase p38 activation in rat. Neurosci Lett. 2001, 298: 207-211. 10.1016/S0304-3940(00)01766-3

  24. 24.

    Chen CH, Beard RS, Bearden SE: Homocysteine impairs endothelial wound healing by activating metabotropic glutamate receptor 5. Microcirculation. 2012, 19: 285-295. 10.1111/j.1549-8719.2012.00159.x

  25. 25.

    Chen Z, Du ZD, Liu JF, Lu DX, Li L, Guan YQ, Wan SG: Endothelial progenitor cell transplantation ameliorates elastin breakdown in a Kawasaki disease mouse model. Chin Med J (Engl). 2012, 125: 2295-2301.

  26. 26.

    Lin YJ, Chang JS, Liu X, Hung CH, Lin TH, Huang SM, Jeang KT, Chen CY, Liao CC, Lin CW, Lai CH, Tien N, Lan YC, Ho MW, Chien WK, Chen JH, Huang YC, Tsang H, Wu JY, Chen CH, Chang LC, Tsai FJ: Association between GRIN3A Gene Polymorphism in Kawasaki Disease and Coronary Artery Aneurysms in Taiwanese Children. PLoS One. 2013, 8: e81384. 10.1371/journal.pone.0081384

  27. 27.

    Burgner D, Curtis N: Kawasaki disease as a cause of encephalitis. Arch Dis Child. 2011, 96: 988-989.

  28. 28.

    Sander T, Hildmann T, Kretz R, Furst R, Sailer U, Bauer G, Schmitz B, Beck-Mannagetta G, Wienker TF, Janz D: Allelic association of juvenile absence epilepsy with a GluR5 kainate receptor gene (GRIK1) polymorphism. Am J Med Genet. 1997, 74: 416-421. 10.1002/(SICI)1096-8628(19970725)74:4<416::AID-AJMG13>3.0.CO;2-L

  29. 29.

    Izzi C, Barbon A, Kretz R, Sander T, Barlati S: Sequencing of the GRIK1 gene in patients with juvenile absence epilepsy does not reveal mutations affecting receptor structure. Am J Med Genet. 2002, 114: 354-359. 10.1002/ajmg.10254

  30. 30.

    Shibata H, Joo A, Fujii Y, Tani A, Makino C, Hirata N, Kikuta R, Ninomiya H, Tashiro N, Fukumaki Y: Association study of polymorphisms in the GluR5 kainate receptor gene (GRIK1) with schizophrenia. Psychiatr Genet. 2001, 11: 139-144. 10.1097/00041444-200109000-00005

  31. 31.

    Hirata Y, Zai CC, Souza RP, Lieberman JA, Meltzer HY, Kennedy JL: Association study of GRIK1 gene polymorphisms in schizophrenia: case–control and family-based studies. Hum Psychopharmacol. 2012, 27: 345-351. 10.1002/hup.2233

  32. 32.

    Kranzler HR, Gelernter J, Anton RF, Arias AJ, Herman A, Zhao H, Burian L, Covault J: Association of markers in the 3′ region of the GluR5 kainate receptor subunit gene to alcohol dependence. Alcohol Clin Exp Res. 2009, 33: 925-930. 10.1111/j.1530-0277.2009.00913.x

  33. 33.

    Kranzler HR, Covault J, Feinn R, Armeli S, Tennen H, Arias AJ, Gelernter J, Pond T, Oncken C, Kampman KM: Topiramate treatment for heavy drinkers: moderation by a GRIK1 polymorphism. Am J Psychiatry. 2014, 171: 445-452. 10.1176/appi.ajp.2013.13081014

  34. 34.

    Kranzler HR, Armeli S, Feinn R, Tennen H, Gelernter J, Covault J: GRIK1 Genotype moderates topiramate’s effects on daily drinking level, expectations of alcohol’s positive effects and desire to drink. Int J Neuropsychopharmacol. 2014, 17: 1549-1556. 10.1017/S1461145714000510

  35. 35.

    Ray LA, Miranda R, MacKillop J, McGeary J, Tidey JW, Rohsenow DJ, Gwaltney C, Swift RW, Monti PM: A preliminary pharmacogenetic investigation of adverse events from topiramate in heavy drinkers. Exp Clin Psychopharmacol. 2009, 17: 122-129.

  36. 36.

    Li S, Qian J, Yang Y, Zhao W, Dai J, Bei JX, Foo JN, McLaren PJ, Li Z, Yang J, Shen F, Liu L, Pan S, Wang Y, Li W, Zhai X, Zhou B, Shi L, Chen X, Chu M, Yan Y, Wang J, Cheng S, Shen J, Jia W, Liu J, Wen Z, Li A, Zhang Y, Zhang G: GWAS identifies novel susceptibility loci on 6p21.32 and 21q21.3 for hepatocellular carcinoma in chronic hepatitis B virus carriers. PLoS Genet. 2012, 8: e1002791. 10.1371/journal.pgen.1002791

  37. 37.

    Lawand NB, McNearney T, Westlund KN: Amino acid release into the knee joint: key role in nociception and inflammation. Pain. 2000, 86: 69-74. 10.1016/S0304-3959(99)00311-5

  38. 38.

    Piepoli T, Mennuni L, Zerbi S, Lanza M, Rovati LC, Caselli G: Glutamate signaling in chondrocytes and the potential involvement of NMDA receptors in cell proliferation and inflammatory gene expression. Osteoarthritis Cartilage. 2009, 17: 1076-1083. 10.1016/j.joca.2009.02.002

  39. 39.

    Flood S, Parri R, Williams A, Duance V, Mason D: Modulation of interleukin-6 and matrix metalloproteinase 2 expression in human fibroblast-like synoviocytes by functional ionotropic glutamate receptors. Arthritis Rheum. 2007, 56: 2523-2534. 10.1002/art.22829

  40. 40.

    Lindblad SS, Mydel P, Hellvard A, Jonsson IM, Bokarewa MI: The N-methyl-d-aspartic acid receptor antagonist memantine ameliorates and delays the development of arthritis by enhancing regulatory T cells. Neurosignals. 2012, 20: 61-71. 10.1159/000329551

  41. 41.

    Miller KE, Hoffman EM, Sutharshan M, Schechter R: Glutamate pharmacology and metabolism in peripheral primary afferents: physiological and pathophysiological mechanisms. Pharmacol Ther. 2011, 130: 283-309. 10.1016/j.pharmthera.2011.01.005

  42. 42.

    Burns JC: Commentary: translation of Dr. Tomisaku Kawasaki’s original report of fifty patients in 1967. Pediatr Infect Dis J. 2002, 21: 993-995. 10.1097/00006454-200211000-00002

  43. 43.

    Sharp CD, Houghton J, Elrod JW, Warren A, Jackson TH, Jawahar A, Nanda A, Minagar A, Alexander JS: N-methyl-D-aspartate receptor activation in human cerebral endothelium promotes intracellular oxidant stress. Am J Physiol Heart Circ Physiol. 2005, 288: H1893-1899.

  44. 44.

    Yoshio T, Okamoto H, Hirohata S, Minota S: IgG anti-NR2 glutamate receptor autoantibodies from patients with systemic lupus erythematosus activate endothelial cells. Arthritis Rheum. 2013, 65: 457-463. 10.1002/art.37745

  45. 45.

    Weng KP, Ou SF, Lin CC, Hsieh KS: Recent advances in the treatment of Kawasaki disease. J Chin Med Assoc. 2011, 74: 481-484. 10.1016/j.jcma.2011.09.001

  46. 46.

    Newburger JW, Takahashi M, Gerber MA, Gewitz MH, Tani LY, Burns JC, Shulman ST, Bolger AF, Ferrieri P, Baltimore RS, Wilson WR, Baddour LM, Levison ME, Pallasch TJ, Falace DA, Taubert KA: Diagnosis, treatment, and long-term management of Kawasaki disease: a statement for health professionals from the Committee on Rheumatic Fever, Endocarditis, and Kawasaki Disease, Council on Cardiovascular Disease in the Young, American Heart Association. Pediatrics. 2004, 114: 1708-1733. 10.1542/peds.2004-2182

  47. 47.

    Falcini F: Kawasaki disease. Curr Opin Rheumatol. 2006, 18: 33-38. 10.1097/01.bor.0000197998.50450.f6

  48. 48.

    Lin YJ, Chang JS, Liu X, Lin TH, Huang SM, Liao CC, Lin CW, Chien WK, Chen JH, Wu JY, Chen CH, Chang LC, Tsang H, Jeang KT, Chen CY, Tsai FJ: Sorting nexin 24 genetic variation associates with coronary artery aneurysm severity in Kawasaki disease patients. Cell Biosci. 2013, 3: 44. 10.1186/2045-3701-3-44

  49. 49.

    Lin YJ, Lan YC, Lai CH, Lin TH, Huang SM, Liao CC, Lin CW, Hung CH, Tien N, Liu X, Chien WK, Chen JH, Tsai FJ: Association of Promoter Genetic Variants in Interleukin-10 and Kawasaki Disease With Coronary Artery Aneurysms. J Clin Lab Anal. 2014, 28: 461-464. 10.1002/jcla.21710

  50. 50.

    Kim S, Dedeoglu F: Update on pediatric vasculitis. Curr Opin Pediatr. 2005, 17: 695-702. 10.1097/01.mop.0000187190.97166.83

  51. 51.

    Matsubara T, Furukawa S, Yabuta K: Serum levels of tumor necrosis factor, interleukin 2 receptor, and interferon-gamma in Kawasaki disease involved coronary-artery lesions. Clin Immunol Immunopathol. 1990, 56: 29-36. 10.1016/0090-1229(90)90166-N

  52. 52.

    Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K: dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001, 29: 308-311. 10.1093/nar/29.1.308

  53. 53.

    Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457

Download references

Acknowledgments

The authors wish to thank the Department of Pediatrics, China Medical University Hospital (CMUH) for administrative assistance and China Medical University (CMU) under the Aim for Top University Plan of the Ministry of Education, Taiwan. We also thank Drs. Kuan-Teh Jeang, Chia-Yen Chen, and Willy W. L. Hong for technical help and suggestions.

Funding

Financial support for this research was provided by CMU (CMU100-S-01), CMUH (DMR-103-039), and the Republic of China National Science Council (NSC100-2320-B-039-012-MY3).

Author information

Correspondence to Fuu-Jen Tsai.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

YJL, JSC, XL, and FJT conceived and designed the experiments. THL, CCL, SMH, CWL and HT performed the experiments. WKC and JHC analyzed the data. JSC, XL, JYW, CHC, LCC, and TJH contributed reagents/materials/analysis tools. YJL and XL wrote the manuscript. All the authors have read and approved the final manuscript.

Electronic supplementary material

13578_2014_193_MOESM1_ESM.ppt

Additional file 1: Figure S1.: Search results of single nucleotide polymorphism (SNP) of rs466013 of the GRIK1 gene used in this study (http://www.ncbi.nlm.nih.gov/projects/SNP/snp_ref.cgi?rs=466013). Above: Genomic location of rs466013 (pointed by red arrows; the NCBI Assembly database: GRCh37.p10 version). Down: Genomic location of rs466013 for 6 versions of the NCBI Assembly database (pointed by red arrows). Figure S2. Search results of single nucleotide polymorphisms (SNPs) of rs466013, rs425507, rs388700 and rs402280 of the GRIK1 gene used in this study (http://genome.ucsc.edu/cgi-bin/hgTracks?db=hg18&position=chr21%3A30120300-30129700&hgsid=370279953_3haDCdtlwLEpPqkcmUFYdaAFYNhx). Above: Genomic location of the GRIK1 gene. Down: Genomic location of rs466013, rs425507, rs388700 and rs402280 (pointed by red arrows; the NCBI Assembly database: NCBI36/hg18 version). Figure S3. GRIK1 mRNA expression levels in peripheral blood mononuclear cells according to the GRIK1 SNPs (rs388700 and rs402280) genotypes. The relative GRIK1 expression was detected by quantitative real-time RT-PCR, and expression from individuals with TT + TA genotypes was compared to that from individuals with AA genotype. The GRIK1 (NM_000830.3) primer sequences were 5′-gcggttagagatggatcaaca-3′ (located at nucleotide 2559–2579 of the transcript (NM_000830.3) and 5′-tcatgaaagcccacatcttct-3′ (located at nucleotide 2617–2637 of the transcript (NM_000830.3)). The relative expression levels were expressed as GRIK1 mRNA/ HPRT mRNA ratio. Figure S4. Venn diagram of 4 GWAS studies. Gene SNPs from 4 GWAS studies were used for searching for common gene SNPs by using Venny website (http://bioinfogp.cnb.csic.es/tools/venny/). Figure S5. Venn diagram of 4 GWAS studies. Gene SNPs from 4 GWAS studies were used for searching for common gene SNPs by using Venny website (http://bioinfogp.cnb.csic.es/tools/venny/). Figure S5. Venn diagram of 4 GWAS studies. Gene SNPs from 4 GWAS studies were used for searching for common gene SNPs by using Venny website (http://bioinfogp.cnb.csic.es/tools/venny/). (PPT 1 MB)

Additional file 2: Table S1.: Characteristics of GWAS studies for KD susceptibility included in this meta-analysis. Table S2. Meta-analysis for previous reported GWAS studies for KD susceptibility. (DOC 54 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • KD
  • GRIK1
  • Single nucleotide polymorphism
  • CAA