Genomics Inform Search


Genomics Inform > Volume 19(4); 2021 > Article
Kim, Hong, Yun, Lee, Yoon, Han, Kim, Kil, Song, Lee, Ha, Jun, Choi, Oh, Yu, Jang, Lee, and The Korean Kawasaki Disease Genetics Consortium: Identification of rare coding variants associated with Kawasaki disease by whole exome sequencing


Kawasaki disease (KD) is an acute pediatric vasculitis that affects genetically susceptible infants and children. To identify coding variants that influence susceptibility to KD, we conducted whole exome sequencing of 159 patients with KD and 902 controls, and performed a replication study in an independent 586 cases and 732 controls. We identified five rare coding variants in five genes (FCRLA, PTGER4, IL17F, CARD11, and SIGLEC10) associated with KD (odds ratio [OR], 1.18 to 4.41; p = 0.0027–0.031). We also performed association analysis in 26 KD patients with coronary artery aneurysms (CAAs; diameter > 5 mm) and 124 patients without CAAs (diameter < 3 mm), and identified another five rare coding variants in five genes (FGFR4, IL31RA, FNDC1, MMP8, and FOXN1), which may be associated with CAA (OR, 3.89 to 37.3; p = 0.0058–0.0261). These results provide insights into new candidate genes and genetic variants potentially involved in the development of KD and CAA.


Kawasaki disease (KD) is an acute, self-limited vasculitis that occurs predominantly in infants and young children [1]. KD is often complicated by coronary artery aneurysms (CAAs), which occur in approximately 15%–25% of untreated children and 3%–5% of children treated with high-dose intravenous immunoglobulin (IVIG) [2-4]. Although the pathogenesis of KD remains unknown, it has been postulated that KD is an abnormal immunological reaction to an infection in genetically susceptible individuals [5,6].
Many complex traits are shaped by intricate interplay between common and rare genetic variants. Recently, genome-wide association studies (GWAS) of common single nucleotide polymorphisms (SNPs) have identified several variants associated with susceptibility to KD [7-11]. Although identification of these common variants has contributed somewhat to understanding the biology of KD, many complexities of the disease remain to be unraveled [12]. Furthermore, most GWAS-derived SNPs do not directly influence the disease trait. Rather, they are index markers in linkage disequilibrium with causative variants for the disease; hence, it is necessary to identify rare coding variants that affect susceptibility to KD using other approaches, such as next-generation sequencing technology. Whole exome sequencing (WES) is a powerful technique for identification of rare protein-coding variants. Here, we used WES and association analysis to identify the protein-coding variants that influence susceptibility to KD and CAA.



Children with KD were recruited from 12 tertiary academic hospitals in Korea, all of which participated in the Korean Kawasaki Disease Genetics Consortium. All patients were diagnosed by pediatricians in accordance with the diagnostic criteria of the American Heart Association [13,14]. All patients received a single high-dose of IVIG (2 g/kg). Two-dimensional echocardiography results were interpreted by pediatric cardiologists, and coronary arteries were categorized as normal or abnormal (showing coronary artery lesions (i.e., dilation or aneurysm)). A total of 200 children with KD were initially included in the WES analysis. Our WES data were deposited to the Biobank for Health Sciences, Chungwon, Korea (; accession number 2016-ER7401-00). For association analysis of WES data from KD cases, a total of 902 control WES or whole genome sequencing (WGS) data samples were obtained from Dr. Yeon-Mok Oh (Asan Medical Center, Seoul, Korea; WES data from 147 patients with chronic obstructive pulmonary disease), Dr. Byung-Ok Choi (Samsung Medical Center, Seoul, Korea; WES data from 258 patients with Charcot–Marie–Tooth disease), and the Biobank for Health Sciences (Center for Genome Science, Chungwon, Korea; WES data from 100 healthy subjects and WGS data from 397 healthy adults). For replication study, a total of 586 patients with complete KD were selected from our KD patient pool. The healthy adult control samples (n = 732) were obtained from the Biobank for Health Sciences (Chungwon, Korea). This study was approved by the institutional review board (IRB) at Asan Medical Center (IRB No. 2014-0823), and informed consent was obtained from the parents of patients with KD.

Targeted capture and exome sequencing

WES was conducted using genomic DNA samples obtained from 200 children with KD. For exome sequencing, 3 μg genomic DNA was initially sheared into 150-base-pair (bp) fragments for shotgun library construction. Fragmented DNA from each sample was captured using the SureSelect Human All Exon V5 Kit (Agilent Technologies, Santa Clara, CA, USA), designed to cover 50 Mb of human genomic sequences corresponding to the NCBI Consensus Coding Sequence database. Each captured library was sequenced at Macrogen (Seoul, Korea) using an Illumina HiSeq 2500 instrument (Illumina, San Diego, CA, USA), according to the manufacturer’s instructions.

Read mapping and variant analysis

Sequence reads in FASTQ format were aligned to a human genome reference sequence (GRCh37) using Burrows–Wheeler Aligner (BWA) software [15]. Duplicate reads were marked in BAM files using Picard MarkDuplicates ( Local realignment, base quality recalibration, variant calling, joint genotyping, and variant quality score recalibration and filtration were applied using Genome Analysis ToolKit (GATK) v3.7 ( Variant annotation was performed using SnpEff [16].


Genotyping for the replication study was performed on the high-throughput Fluidigm EP1 system (Fluidigm Corp., South San Francisco, CA, USA), using the Fluidigm SNPType assay platform, and nanofluidic 96.96 and 48.48 Dynamic Array Integrated Fluidic Circuits (Fluidigm). Genotype calls were made using the Fluidigm SNP Genotyping Analysis program (Fluidigm). All but one of the selected SNPs were genotyped for replication. The overall genotyping success rate was 99.5%.

SNP filtering

After WES, a total of 870,062 SNPs were identified from 200 KD cases and 902 controls. To filter the SNP markers, we excluded 337,282 SNP markers with missing call rates > 3%, 138,329 monomorphic markers, and 3,404 multiallelic markers in patients with KD. As coding variants that cause amino acid changes may be hypothesized to confer larger phenotypic effects, we also excluded 313,470 intergenic, intronic, untranslated-region, and synonymous SNPs. After filtering, 77,577 protein-altering variants remained, among which 75,365 were missense, 1,448 were nonsense, and 764 were in splice-acceptor/donor sites. Next, we excluded 3,209 SNPs with missing call rates > 5% (slightly looser threshold due to heterogeneous control data compared to case data) in control data. Finally, a total of 74,368 SNPs were included in the association analysis.

Statistical analysis

All statistical analyses were performed using PLINK software v1.07 ( [17]. The chi-square or Fisher’s exact tests (when the expected cell count was less than five) were used to compare SNP allele frequencies between cases and controls to test for associations with KD and CAA.


Associations with KD

WES was initially performed using samples from 200 children with KD, of which 159 had complete KD with a fever lasting ≥ 5 days and > 4 of the principal symptoms of KD. Therefore, we performed association analysis of 159 complete KD cases and 902 controls. For the replication study, we selected 12 coding SNPs with significant associations with KD susceptibility (chi-square test p < 0.001) and immune-related functions. Extremely rare variants are more likely to be deleterious; thus we also selected an additional 51 immune-related rare coding SNPs detected in case samples only, or with odds ratios (ORs) ≥ 3 in case-control comparisons. We tested these 63 SNPs in 586 patients with complete KD and 732 healthy controls using the Fluidigm SNPType assay platform; however, none showed statistical significance (p > 0.05) due to their very low frequencies in both case and control samples; however, among them, five coding SNPs (rs2275603 in FCRLA, rs755244149 in PTGER4, rs117796773 in IL17F, rs41493047 in CARD11, and rs201376644 in SIGLEC10) showed nominal significances (OR = 1.18 to 4.41, p = 0.0027–0.031) in combined analysis of data generated by WES and Fluidigm SNPType assay (Table 1).

Associations with CAA

Of the 200 children with KD, 74 had coronary artery lesions. To investigate genetic loci that affect the development of CAAs in KD, we initially stratified our KD samples by CAA size, according to the American Heart Association statement [13], as follows: small, 3 to <5 mm; medium, 5 to 8 mm, or giant, >8 mm. Next, we performed association analysis in 124 KD patients without CAAs (<3 mm) and 26 KD patients with medium or giant aneurysms that showed more extreme phenotypes (internal lumen diameter of coronary arteries > 5 mm). For the replication study, we selected five SNPs with significant associations with CAAs (chi-square test p < 0.001) and CAA-related functions, and tested them in 420 patients with KD without CAAs and 19 with CAAs (diameter > 5 mm). Although no SNP showed significant results in the replication study, five coding variants (rs148721785 in IL31RA, rs201812753 in FGFR4, rs374967242 in FNDC1, rs61754773 in MMP8, and rs188424977 in FOXN1) showed significant associations (OR, 3.89–37.3; p = 0.0058 to 0.0261) in combined analysis of the WES and replication samples (Table 2).


Recent GWAS identified some clear KD-associated loci and successfully identified part of the genetic background of KD; however, knowledge of these genetic factors has not made a significant contribution to understanding KD pathogenesis [12]. In this study, we performed WES with the aim of identifying rare protein-coding variants that contribute to KD susceptibility. We identified five rare coding variants in five genes (FCRLA, PTGER4, IL17F, CARD11, and SIGLEC10) associated with KD and another five coding variants in five genes (FGFR4, IL31RA, FNDC1, MMP8, and FOXN1) associated with CAA.
Among the five KD-associated genes, FCRLA has previously been reported to be highly significantly associated with low levels of affinity immunoglobulin gamma Fc region receptor II-a/b (products of FCGR2A and FCGR2B) in human blood (rs6668534: 0.434-unit decrease, p = 6 × 10-17) [18], as well as being associated with monocyte and lymphocyte counts [19]. FCGR2A is an established KD susceptibility gene [7,10]. These observations suggest that the effect of FCGR2A gene variants on KD susceptibility may be regulated by FCRLA. PTGER4 is one of four identified receptors for prostaglandin E2 (PGE2). Previous studies investigated the function of PGE2 in relation to KD [20-23]. PGE2 levels are markedly elevated during the acute stage of KD and then decrease during the recovery stage [20]. PGE2 level was also correlated with prevention of IVIG resistance and CAA formation through CD40L [23]. Thus, PTGER4, a PGE2 receptor, may function in KD pathogenesis. The role of the CARD11 gene in KD remains unknown; however, mutations of CARD11 can cause a condition referred to as ‘B-cell expansion with NFκB and T cell anergy’, and immunodeficiency 11, due to defective B-cell receptor signaling [24,25]. CARD11 is also associated with various immune-related diseases and traits, including multiple sclerosis, eczema, and platelet count [19,26,27]. Furthermore, defective B-cell development and B-cell functions are well known to be involved in KD pathogenesis [28,29]. Therefore, CARD11 may play a role in susceptibility to KD. The roles of IL17F and SIGLEC10 in KD susceptibility are unknown.
To date, several GWAS have been performed to identify susceptibility loci for coronary artery formation in KD, using a very limited number of patients with KD and CAA [30-34]. No overlapping genetic loci were detected from these GWAS. In this study, we identified five rare coding variants in five genes associated with CAAs in KD. None of these five genes, except for MMP8, have been reported to be associated with susceptibility to KD or CAAs in KD. MMP8 encodes a member of the matrix metalloproteinase (MMP) family of proteins, which may be involved in maintaining the structure and function of coronary vascular walls [35]. An increasing number of studies have indicated elevated expression, activity, or protein levels of MMPs in KD [36,37]. MMP8 expression levels are significantly upregulated in KD before IVIG treatment and decrease after undergoing IVIG treatment [37]. MMP8 is also highly abundant in IVIG-resistant KD, suggesting that this is a potential biomarker for identifying those KD patients at highest risk for CAAs, and who may benefit from additional anti-inflammatory therapy [38,39].
To our knowledge, this is the first WES study of KD, and we identified several rare coding risk variants potentially associated with KD and CAA through WES and association analysis; however, we failed to confirm the significance of the selected SNPs in the replication study, mainly due to their extremely low frequency. Another reason for lack of significance in a replication study can be partly due to sample selection bias. To enrich for possible causative variants in the WES stage, we selected more severe KD cases than in the replication stage, such as patients with medium-sized or giant aneurysms (6.3% in WES samples vs. 1.7% in replication samples), IVIG resistance (50.3% vs. 22.5%), recurrence (14.5% vs. 5.5%), and family history (6.3% vs. 1.9%). Therefore, more variants that result in these severe phenotypes may have been discovered in the WES stage; however, in combined analysis of WES and replication data, statistically significant associations were detected. WES is a powerful technique for identification of protein-coding variants, which are more straightforward to annotate for biological functions and help to pinpoint causative genes; however, functional variants are more likely to be rare. Therefore, large sample sizes are usually required to detect loci with exome-wide significance, particularly for highly polygenic traits; however, we were only able to sequence a small number of samples, particularly in the case of children with CAA. Both common and rare variant association approaches are drastically underpowered to detect associations in small samples. Thus, our conclusions are tentative and preliminary, and further studies with large sample size are necessary to confirm our findings, as well as the functional role of the newly identified genes and variants in KD pathogenesis.
In summary, we sequenced the whole exomes of patients with KD and identified five candidate genes (FCRLA, PTGER4, IL17F, CARD11, and SIGLEC10) that may influence KD susceptibility, and another five candidate genes (FGFR4, IL31RA, FNDC1, MMP8, and FOXN1) with potential influence on the development of coronary aneurysms in KD. These results provide insights into new candidate genes and genetic variants potentially involved in the development of KD and CAA. Further association studies with expanded KD sample numbers from the Korean population or from cohorts with different ethnicities are required to confirm these results.


Authors’ Contribution

Conceptualization: JJY, GYJ, JKL. Data curation: JJK, YMH, SWY, KYL, KLY, MKH, GBK, HRK, MSS, HDL, KSH, HOJ, BOC, YMO, JJY, GYJ. Formal analysis: JJK. Funding acquisition: JKL. Methodology: JJK. Writing - original draft: JJK. Writing - review & editing: JJY, GYJ, JKL.

Conflicts of Interest

No potential conflict of interest relevant to this article was reported.


We thank all of our patients and their families for participating in this study. This work was supported by a grant from the Ministry of Health & Welfare of the Republic of Korea (HI15C1575) and a grant from the Korea Center for Disease Control and Prevention (2016-ER7401-00).

Table 1.
SNPs associated with KD
CHR SNP Gene AA change Allelea
WES (KD = 159 vs. Control = 902)
Replication (KD = 586 vs. Control = 732)
Combined (KD = 745 vs. Control = 1,634)
Genotype 11/12/22 (RAF)
OR p Genotype (RAF)
OR p Genotype (RAF)
OR p
1 2 KD Control KD Control KD Control
1 rs2275603 FCRLA p.Ser209Gly G A 58/67/34 (0.575) 202/432/268 (0.463) 1.56 0.00023b 151/281/151 (0.500) 166/382/181 (0.490) 1.04 0.600 209/348/185 (0.516) 368/814/449 (0.475) 1.18 0.0088b
5 rs755244149 PTGER4 p.Arg358His A G 0/3/156 (0.009) 0/2/900 (0.001) 8.58 0.026b,c 0/3/580 (0.003) 0/1/728 (0.001) 3.76 0.329c 0/6/736 (0.004) 0/3/1628 (0.001) 4.41 0.031b,c
6 rs117796773 IL17F Splice donor: c.254+1G>T A C 0/13/146 (0.041) 0/20/882 (0.011) 3.80 7.52E-05b 0/22/561 (0.019) 0/19/710 (0.013) 1.46 0.231 0/35/707 (0.024) 0/39/1592 (0.012) 2.00 0.0027b
7 rs41493047 CARD11 p.Ala687Val A G 0/9/150 (0.028) 0/13/889 (0.007) 4.01 0.00062b 0/12/571 (0.010) 0/11/718 (0.008) 1.37 0.453 0/21/721 (0.014) 0/24/1607 (0.007) 1.94 0.025b
19 rs201376644 SIGLEC10 p.Asn292Asp C T 0/18/141 (0.057) 0/22/880 (0.012) 4.86 7.91E-08b 0/4/579 (0.003) 0/1/728 (0.001) 5.02 0.109 0/22/720 (0.015) 0/23/1608 (0.007) 2.12 0.010b

Chi-square or Fisher’s exact tests were used, as appropriate.

SNP, single nucleotide polymorphism; KD, Kawasaki disease; CHR, chromosome; AA, amino acid; OR, odds ratio; RAF, risk allele frequency.

a Allele 1 refers to a risk allele.

b Significant p-values.

c p-values of the exact test.

Table 2.
SNPs associated with CAA (≥5 mm)
CHR SNP Gene AA change Allelea
WES (CAA = 26 vs. Non-CAA = 124)
Replication (CAA = 19 vs. Non-CAA = 420)
Combined (CAA = 45 vs. Non-CAA = 544)
Genotype 11/12/22 (RAF)
OR P Genotype (RAF)
OR P Genotype (RAF)
1 2 CAAb Non-CAAc CAAb Non-CAAc CAAb Non-CAAc
5 rs148721785 IL31RA p.Val672Leu C G 0/6/20 (0.115) 0/5/119 (0.020) 6.34 0.0048d 0/1/18 (0.026) 0/18/400 (0.022) 1.23 0.574 0/7/38 (0.078) 0/23/519 (0.021) 3.89 0.0058d
5 rs201812753 FGFR4 p.Pro374Leu T C 0/4/22 (0.077) 0/1/123 (0.004) 20.6 0.0036d 0/0/19 (0.000) 0/8/410 (0.010) NA 1 0/4/41 (0.044) 0/9/533 (0.008) 5.56 0.0135d
6 rs374967242 FNDC1 p.Trp59Arg C T 0/3/23 (0.058) 0/0/124 (0.000) NA 0.0050d 0/0/19 (0.000) 0/1/417 (0.001) NA 1 0/3/42 (0.033) 0/1/541 (0.001) 37.3 0.0017d
11 rs61754773 MMP8 p.Asp300Asn T C 0/3/23 (0.058) 0/1/123 (0.004) 15.1 0.0174d 0/1/18 (0.026) 0/10/408 (0.012) 2.23 0.389 0/4/41 (0.044) 0/11/531 (0.010) 4.54 0.0229d
17 rs188424977 FOXN1 p.Gly238Asp A G 0/3/23 (0.058) 0/0/124 (0.000) NA 0.0050d 0/0/19 (0.000) 0/6/412 (0.007) NA 1 0/3/42 (0.033) 0/6/536 (0.006) 6.20 0.0261d

Fisher’s exacts test was used.

SNP, single nucleotide polymorphism; CAA, coronary artery aneurysms; CAA, coronary artery aneurysms; CHR, chromosome; AA, amino acid; WES, whole exome sequencing; RAF, risk allele frequency; OR, odds ratio; NA, not available.

a Allele 1 refers to a risk allele.

b Kawasaki disease (KD) patients with CAAs (internal diameter ≥ 5 mm)

c KD patients without CAAs.

d Significant p-values.


1. Burns JC, Glode MP. Kawasaki syndrome. Lancet 2004;364:533–544.
crossref pmid
2. Kato H, Sugimura T, Akagi T, Sato N, Hashino K, Maeno Y, et al. Long-term consequences of Kawasaki disease: a 10- to 21-year follow-up study of 594 patients. Circulation 1996;94:1379–1385.
3. Newburger JW, Takahashi M, Burns JC, Beiser AS, Chung KJ, Duffy CE, et al. The treatment of Kawasaki syndrome with intravenous gamma globulin. N Engl J Med 1986;315:341–347.
crossref pmid
4. Durongpisitkul K, Gururaj VJ, Park JM, Martin CF. The prevention of coronary artery aneurysm in Kawasaki disease: a meta-analysis on the efficacy of aspirin and immunoglobulin treatment. Pediatrics 1995;96:1057–1061.
crossref pmid
5. Newburger JW, Fulton DR. Kawasaki disease. Curr Opin Pediatr 2004;16:508–514.
crossref pmid
6. Lee KY, Han JW, Lee JS. Kawasaki disease may be a hyperimmune reaction of genetically susceptible children to variants of normal environmental flora. Med Hypotheses 2007;69:642–651.
crossref pmid pmc
7. Khor CC, Davila S, Breunis WB, Lee YC, Shimizu C, Wright VJ, et al. Genome-wide association study identifies FCGR2A as a susceptibility locus for Kawasaki disease. Nat Genet 2011;43:1241–1246.
crossref pmid
8. Lee YC, Kuo HC, Chang JS, Chang LY, Huang LM, Chen MR, et al. Two new susceptibility loci for Kawasaki disease identified through genome-wide association analysis. Nat Genet 2012;44:522–525.
crossref pmid
9. Onouchi Y, Ozaki K, Burns JC, Shimizu C, Terai M, Hamada H, et al. A genome-wide association study identifies three new risk loci for Kawasaki disease. Nat Genet 2012;44:517–521.
crossref pmid
10. Kim JJ, Yun SW, Yu JJ, Yoon KL, Lee KY, Kil HR, et al. A genome-wide association analysis identifies NMNAT2 and HCP5 as susceptibility loci for Kawasaki disease. J Hum Genet 2017;62:1023–1029.
crossref pmid
11. Kwon YC, Kim JJ, Yun SW, Yu JJ, Yoon KL, Lee KY, et al. BCL2L11 is associated with Kawasaki disease in intravenous immunoglobulin responder patients. Circ Genom Precis Med 2018;11:e002020.
12. Onouchi Y. The genetics of Kawasaki disease. Int J Rheum Dis 2018;21:26–30.
crossref pmid
13. Newburger JW, Takahashi M, Gerber MA, Gewitz MH, Tani LY, Burns JC, et al. 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. Circulation 2004;110:2747–2771.
crossref pmid
14. McCrindle BW, Rowley AH, Newburger JW, Burns JC, Bolger AF, Gewitz M, et al. Diagnosis, treatment, and long-term management of Kawasaki disease: a scientific statement for health professionals from the American Heart Association. Circulation 2017;135:e927–e999.
crossref pmid
15. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009;25:1754–1760.
crossref pmid pmc
16. Cingolani P, Platts A, Wang le L, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin) 2012;6:80–92.
crossref pmid pmc
17. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007;81:559–575.
crossref pmid pmc
18. Suhre K, Arnold M, Bhagwat AM, Cotton RJ, Engelke R, Raffler J, et al. Connecting genetic risk to disease end points through the human blood plasma proteome. Nat Commun 2017;8:14357.
crossref pmid pmc
19. Vuckovic D, Bao EL, Akbari P, Lareau CA, Mousas A, Jiang T, et al. The polygenic and monogenic basis of blood traits and diseases. Cell 2020;182:1214–1231.
crossref pmid pmc
20. Lee T, Furukawa S, Fukuda Y, Yabuta K, Kato H. Plasma prostaglandin E2 level in Kawasaki disease. Prostaglandins Leukot Essent Fatty Acids 1988;31:53–57.
crossref pmid
21. Kajimoto M, Ichiyama T, Ueno Y, Shiraishi M, Hasegawa M, Furukawa S. Enhancement of activated beta1-integrin expression by prostaglandin E2 via EP receptors in isolated human coronary arterial endothelial cells: implication for the treatment of Kawasaki disease. Inflamm Res 2009;58:224–228.
crossref pmid
22. Sasai K. Plasma PGE2, TXB2 and 6-keto PGF1 alpha levels in patients with Kawasaki disease. Arerugi 1988;37:952–958.
23. Kuo HC, Wang CL, Yang KD, Lo MH, Hsieh KS, Li SC, et al. Plasma prostaglandin E2 levels correlated with the prevention of intravenous immunoglobulin resistance and coronary artery lesions formation via CD40L in Kawasaki disease. PLoS One 2016;11:e0161265.
crossref pmid pmc
24. Snow AL, Xiao W, Stinson JR, Lu W, Chaigne-Delalande B, Zheng L, et al. Congenital B cell lymphocytosis explained by novel germline CARD11 mutations. J Exp Med 2012;209:2247–2261.
crossref pmid pmc
25. Stepensky P, Keller B, Buchta M, Kienzler AK, Elpeleg O, Somech R, et al. Deficiency of caspase recruitment domain family, member 11 (CARD11), causes profound combined immunodeficiency in human subjects. J Allergy Clin Immunol 2013;131:477–485.
26. International Multiple Sclerosis Genetics Consortium, Beecham AH, Patsopoulos NA, Xifara DK, Davis MF, Kemppinen A, et al. Analysis of immune-related loci identifies 48 new susceptibility variants for multiple sclerosis. Nat Genet 2013;45:1353–1360.
crossref pmid pmc
27. Kichaev G, Bhatia G, Loh PR, Gazal S, Burch K, Freund MK, et al. Leveraging polygenic functional enrichment to improve GWAS power. Am J Hum Genet 2019;104:65–75.
crossref pmid
28. Yamazaki-Nakashimada MA, Gamez-Gonzalez LB, Murata C, Honda T, Yasukawa K, Hamada H. IgG levels in Kawasaki disease and its association with clinical outcomes. Clin Rheumatol 2019;38:749–754.
crossref pmid
29. Lindquist ME, Hicar MD. B cells and antibodies in Kawasaki disease. Int J Mol Sci 2019;20:1834.
crossref pmc
30. Kim JJ, Park YM, Yoon D, Lee KY, Seob Song M, Doo Lee H, et al. 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.
crossref pmid
31. Kuo HC, Li SC, Guo MM, Huang YH, Yu HR, Huang FC, et al. Genome-wide association study identifies novel susceptibility genes associated with coronary artery aneurysm formation in Kawasaki disease. PLoS One 2016;11:e0154943.
crossref pmid pmc
32. Lin YJ, Chang JS, Liu X, Tsang H, Chien WK, Chen JH, et al. Genetic variants in PLCB4/PLCB1 as susceptibility loci for coronary artery aneurysm formation in Kawasaki disease in Han Chinese in Taiwan. Sci Rep 2015;5:14762.
crossref pmid pmc
33. Kwon YC, Kim JJ, Yu JJ, Yun SW, Yoon KL, Lee KY, et al. Identification of the TIFAB gene as a susceptibility locus for coronary artery aneurysm in patients with Kawasaki disease. Pediatr Cardiol 2019;40:483–488.
crossref pmid
34. Hoggart C, Shimizu C, Galassini R, Wright VJ, Shailes H, Bellos E, et al. Identification of novel locus associated with coronary artery aneurysms and validation of loci for susceptibility to Kawasaki disease. Eur J Hum Genet 2021;29:1734–1744.
crossref pmid
35. Hopps E, Caimi G. Matrix metalloproteases as a pharmacological target in cardiovascular diseases. Eur Rev Med Pharmacol Sci 2015;19:2583–2589.
36. Senzaki H. The pathophysiology of coronary artery aneurysms in Kawasaki disease: role of matrix metalloproteinases. Arch Dis Child 2006;91:847–851.
crossref pmid pmc
37. Kuo HC, Li SC, Huang LH, Huang YH. Epigenetic hypomethylation and upregulation of matrix metalloproteinase 9 in Kawasaki disease. Oncotarget 2017;8:60875–60891.
crossref pmid pmc
38. Fury W, Tremoulet AH, Watson VE, Best BM, Shimizu C, Hamilton J, et al. Transcript abundance patterns in Kawasaki disease patients with intravenous immunoglobulin resistance. Hum Immunol 2010;71:865–873.
crossref pmid pmc
39. Pan Y, Fan Q. Identification of potential core genes in immunoglobulin-resistant Kawasaki disease using bioinformatics analysis. Crit Rev Eukaryot Gene Expr 2020;30:85–91.
crossref pmid


Browse all articles >

Editorial Office
Room No. 806, 193 Mallijae-ro, Jung-gu, Seoul 04501, Korea
Tel: +82-2-558-9394    Fax: +82-2-558-9434    E-mail:                

Copyright © 2022 by Korea Genome Organization.

Developed in M2PI

Close layer
prev next