Whole-exome sequencing analysis identifies novel variants associated with Kawasaki disease susceptibility
Pediatric Rheumatology volume 21, Article number: 78 (2023)
Kawasaki disease (KD) is an acute pediatric vasculitis affecting genetically susceptible infants and children. Although the pathogenesis of KD remains unclear, growing evidence links genetic susceptibility to the disease.
To explore the genes associated with susceptibility in KD, we applied whole-exome sequencing to KD and control subjects from Yunnan province, China. We conducted association study analysis on the two groups.
In this study, we successfully identified 11 significant rare variants in two genes (MYH14 and RBP3) through the genotype/allele frequency analysis. A heterozygous variant (c.2650G > A, p.V884M) of the RBP3 gene was identified in 12 KD cases, while eight heterozygous variants (c.566G > A, p.R189H; c.1109 C > T, p.S370L; c.3917T > G, p.L1306R; c.4301G > A, p.R1434Q; c.5026 C > T, p.R1676W; c.5329 C > T, p.R1777C; c.5393 C > A, p.A1798D and c.5476 C > T, p.R1826C) of the MYH14 gene were identified in 8 KD cases respectively.
This study suggested that nine variants in MYH14 and RBP3 gene may be associated with KD susceptibility in the population from Yunnan province.
Kawasaki disease (KD; OMIM611775) is an acute, self-limiting systemic vasculitis syndrome with the main clinical manifestations of fever, oral mucosal changes, rash, cervical lymphadenopathy, conjunctivitisin, and extremity changes, known as mucocutaneous lymph node syndrome (MCLS) [1, 2]. It was first described by Japanese pediatrician Tomisaku Kawasaki and particularly affects children under five years of age. With an almost worldwide increase in incidence, KD becomes now the leading cause of acquired heart disease in children, as it may cause coronary artery lesions in 15–25% of untreated patients or in 5–10% of patients treated with intravenous immunoglobulin (IVIG) . The incidence of incomplete KD, which accounts for 15% and 47% of all KD cases , has also been reported to be increasing, posing a threat to the health of children’s coronary arteries [5,6,7].
Although KD can be diagnosed based on its typical features ( such as fever, conjunctivitis, skin rashes, increased fibrinogen, etc.), while the immunopathogenic mechanisms of this disease remain unclear. Common and rare genetic variants could form many complex traits with complex interactions [9,10,11]. Domestic and foreign studies have found that inflammation-related genes IL-18 and IL-1B [1, 12, 13], inositol 1, 4, 5-trisphosphate 3-kinase C (ITPKC) [14, 15], and other gene polymorphisms are associated with KD . Meanwhile, the family aggregation of KD patients indicated that genetic factors play an important role in the occurrence of KD [17,18,19]. However, the susceptibility loci obtained by the candidate gene method have been controversial because the results of various studies cannot lead to more accurate and consistent conclusions due to differences in race, environment, and sample content. Finding susceptibility genes associated with complex diseases at the genome-wide level is an effective approach to investigating polygenic diseases. Commonly used methods include genome-wide linkage studies, genome-wide association studies (GWAS), and whole-exome sequencing (WES) . Most GWAS-derived Single nucleotide polymorphisms (SNPs) do not directly affect disease characteristics, but are an index marker linking disease-specific imbalances and pathogenic variants [16, 21, 22]. Therefore, it is necessary to use other methods to identify rare coding variants that affect KD susceptibility. WES is one of the efficient sequencing techniques to identify rare protein-coding variants. In this study, we determined to identify the KD-associated protein-coding variants through WES that may provide new insights into diagnosis and treatment of KD.
Materials and methods
Patients and samples
The case-control sample set used in this study included 93 KD patients and 91 non-KD control cases. All cases were obtained from the Kunming Children’s Hospital, Yunnan province, southern China, and unrelated to each other. Inclusion criteria for KD cases included patients met the criteria for KD (fever together with principal symptoms such as conjunctivitis, skin rashes, increased fibrinogen, etc.)[8, 23], non-infectious and no previous cancer or metastases. The control cases were non-KD patients without fever. This study was approved by the Ethical Review Board of Kunming Children’s Hospital, and informed consent was obtained from all of KD and control subjects.
Whole-exome sequencing and Sanger sequencing
WES was conducted using genomic DNA samples obtained from 93 children with KD. The exome sequences were efficiently enriched from 1 µg genomic DNA extracted from the peripheral blood using Agilent liquid capture system (Agilent Sure Select Human All Exon V6 kit, Agilent Technologies, Santa Clara, CA, USA) according to the manufacturer’s protocol. Finally, Illumina Novaseq 6000 platform (Illumina, San Diego, CA, USA), with 150 bp pair fragments sequencing mode, was used for sequencing the genomic DNA for shotgun library construction. The overall genotyping success rate was 99.5%. Raw image files were processed using CASAVA v1.82 for base calling and generating raw data.
Sanger sequencing was performed to confirm the variants identified by WES. PCR was conducted with TaKaRa Taq (Takara, Osaka, Japan) under the following conditions: 95 °C for 5 min; followed by 34 cycles at 95 °C for 30 s, 59 °C for 30 s and 72 °C for 30 s; 72 °C for 5 min. PCR products were purified by gel electrophoresis and sequenced using ABI 3730xl DNA Analyzer with the BigDye™ Terminator Cycle Sequencing Kit (Applied Biosystems, Foster, CA, USA).
Variation frequencies were described as proportions, and SNP allele frequency comparisons between cases and controls were analyzed by Fisher’s exact tests and odds ratios , and 95% confidence intervals (CIs) calculated by unconditional logistic regression were used to analyze the association between SNPs and KD susceptibility. Two-tailed p-value < 0.05 was considered statistically significant.
In the KD group, there were 87 classic KD patients (93.55%), 8 incomplete KD patients (8.6%) and 2 unresponsive to IVIG (2.2%). The mean age of KD group was 2.40 ± 1.84 years (ranged from 7 days to 12 years), and the male to female ratio was 1.7:1 (59/34). It was slightly higher than the ratio of 1.5: 1 of KD patients generally, p-value < 0.05 according to Fisher’s exact test. The proportion of KD patients with coronary artery aneurysms was 8.6% (8, 93). The mean age of control group was 2.50 ± 1.76 years (ranged from 3 months to 12 years), and the male to female ratio was 1.6:1 (56/35). In KD and control groups, the proportions of Han were 80.65% (75, 93) and 84.62% (77, 91) respectively. The proportions of ethnicity were 19.35% (18, 93) and 15.38% (14, 91) respectively. The detail information of the patients was in supplement file 1.
Filtering of candidate variants
A total of 349,054 variants were identified from exome sequence data of 93 KD cases and 91 controls. The filtering steps including: filtering by MAF (MAF < 0.01), Picking of damage variants, case-control analysis. MAF filtering according to the population frequency databases include including 1000G(1000 Genome), ExAC(Exome Aggregation Consortium), and gnomAD (Genome Aggregation Database). Frameshift mutations, terminator mutations, splice mutations, and missense mutation with a Combined Annotation Dependent Depletion (CADD)  scores > 25 were retained. In the step, 8413 protein-altering variants remained, among which 6466 were missense SNV, 851 were frameshift, 394 were in splice-acceptor/donor sites, 687 were nonsense, and 15 were start lost.
Associations with KD
ORs were valued by R package ggplot2, which performed deleterious when ORs > 1. Through association analysis of these 8413 variants, we successfully identified 11 variants in three genes (TTI1, MYH14, and RBP3) from 46 cases (Fig. 1; Table 1), which showed nominal significances (ORs = 2.3177 to 13.1963; p = 0.0025 to 0.0346) (Table 1). However, two variants in TTI1 gene (appeared in 26 cases) were excluded as the high allele frequency in control group. A heterozygous variant (c.2650G > A, p.V884M) of the RBP3 gene was identified in 12 KD cases, while eight heterozygous variants (c.566G > A, p.R189H; c.1109 C > T, p.S370L; c.3917T > G, p.L1306R; c.4301G > A, p.R1434Q; c.5026 C > T, p.R1676W; c.5329 C > T, p.R1777C; c.5393 C > A, p.A1798D and c.5476 C > T, p.R1826C) of the MYH14 gene were identified in 8 KD cases respectively. All of the variants were confirmed by Sanger sequencing.
KD is one of the most common systemic vasculitic illness of children under the age of five years, leads to coronary artery aneurysms in 25% of untreated patients [8, 23]. It’s a multisystem inflammatory process, presumably, the etiology is an excessive immune response to possible infection or environmental triggers in genetically susceptible individuals . People with KD may be inherently prone to other diseases, especially children younger than five years. Previous studies indicated that the incidence of KD in Asia was higher than that in the United States and Europe [25,26,27,28,29], and a higher incidence of males than females [2, 30]. Incidence within families is higher than in sporadic cases [17, 29]. KD could be regarded as a multifactorial and polygenic (complex) disorder [31, 32].
GWAS has identified some well-defined KD-associated loci and part of the genetic background successfully in recent studies, while it does not contribute significantly to exploring the pathogenesis of KD [33, 34]. Different from the GWAS, WES technology can explore global genetic mutations of many other complex diseases. It could discover rare mutations in the encoding sequence, which may cause its protein-coding variants that contribute to KD susceptibility. Jae-Jung Kim et al. explore the impact of coding variation on KD using WES for the first time , while no studies in the Chinese population.
In this study, we performed WES to identify rare protein-coding variants responsible for KD susceptibility. Nine variants in RBP3 and MYH14 gene were significantly enriched in KD cases. c.2650G > A (p.V884M) in exon 1 of RBP3 gene were identified in 12 KD cases, appears to be present at a high rate (12.9%) in the KD group. Eight variants of the MYH14 gene were identified in 8 KD cases respectively. All the allele frequencies were lower than 0.0275%, which indicated both of them were rare variants of genes. RBP3 and MYH14 gene are the first time reported to be associated with KD.
RBP3 gene encode interphotoreceptor retinol-binding protein, transport retinoids between the retinal pigment epithelium and the photoreceptors . In 2015, Arno et al. first described retinal dystrophy in children caused by homozygous nonsense RBP3 mutations, highlighting the requirement for IRBP in normal eye development and visual function . Yokomizo et al.  found that elevated expression of photoreceptor-secreted RBP3 may play a role in protection against the progression of diabetic retinopathy. To date, sixteen RBP3 gene variants have been recorded in the HGMD database (https://www.hgmd.cf.ac.uk/ac/index.php), including eleven missense variants, three nonsense variants, one frameshift variant, and one fragment deletion variant (Fig. 2). In the study, c.2650G > A (p.V884M) locates in the third of four tandem homology modules [27, 38, 39], causes amino acid change from Valine to methionine. Chen P et al  indicated c.2650G > A (p.V884M) was associated with corneal curvature in Asian populations. Although previous researchers thought RBP3 was associated with retinal retinoid transport and corneal changes, our research suggests an association with KD also. Conjunctivitis and subconjunctival hemorrhage are common phenotypes in KD [8, 23]. While in this research, c.2650G > A (p.V884M) maybe related to the ocular phenotypes in KD such as conjunctivitis and subconjunctival hemorrhage, which needs further research.
MYH14 is a member of the nonmuscle myosin II family of ATP-dependent molecular motors, which interact with cytoskeletal actin and regulate cytokinesis, cell motility, and cell polarity . Sixty gene variants of MYH14 gene have been recorded in the HGMD database, within 54 missense variants, three nonsense variants, and three frameshift variants. Missense variant is the most common pathogenic variant, scattered across the whole gene. Eight missense variants of MYH14 gene were discovered in this study (Fig. 3), only c.5393 C > A( p.A1798D) has been published before. The homozygous variant (c.5393 C > A) may cause perineal fistulas in Anorectal malformations, based on the genetic and computer analyses, related to normal cloaca development by nonmuscle myosin heavy chain IIC (NMHC IIC) localization analysis . It seems no significantly relevant to the result of this study. Wang M et al  indicated that c.5417 C > A (p.A1806D) in MYH14 gene led to sensorineural hearing loss (SNHL). SNHL was reported in approximately 36% patients with Kawasaki , rare variants in MYH14 gene maybe potentially associate with the symptom.
Above all, RBP3 and MYH14 were first reported to be associated with KD susceptibility in chinses population. This study’s limitation is a relatively small sample size, so the samples of coronary artery aneurysms and ethnic minorities were limited, unable to conduct more diversified data analysis. Studies on large sample sizes are needed in future to further reveal the relationship between candidate genes and KD.
WES was performed on the KD and control group to identify susceptibility genes in patients from southwest China, and two protein-coding gene (RBP3 and MYH14) were identified in the case-control analysis (ORs, 8.3945 to 13.1963; p-value, 0.0346 to 0.0025). These results provide insights into novel candidate genes and genetic variants that may be involved in KD and related KD complications. Further association studies with expanded KD samples from southwestern China or different ethnic groups are needed to confirm these results.
All data generated or analyzed during this study are included in this published article and tables, and whole-exome sequencing analysis data has been uploaded to NCBI database (No. PRJNA869779).
Mucocutaneous lymph node syndrome
Inositol 1, 4, 5-trisphosphate 3-kinase C
Genome-wide association studies
Myosin Heavy Chain 14
Retinoid-binding protein 3
Burns JC, Glodé MP. Kawasaki syndrome. The Lancet. 2004;364(9433):533–44.
Kuo HC, Liang CD, Wang CL, Yu HR, Hwang KP, Yang KD. Serum albumin level predicts initial intravenous immunoglobulin treatment failure in Kawasaki disease. Acta Paediatr. 2010;99(10):1578–83.
Sheu JJ, Lin YJ, Chang JS, Wan L, Chen SY, Huang YC, Chan C, Chiu IW, Tsai FJ. Association of COL11A2 polymorphism with susceptibility to Kawasaki disease and development of coronary artery lesions. Int J Immunogenet. 2010;37(6):487–92.
Kwon YC, Kim JJ, Yun SW, Yu JJ, Yoon KL, Lee KY, Kil HR, Kim GB, Han MK, Song MS, et al. BCL2L11 is Associated with Kawasaki Disease in Intravenous Immunoglobulin Responder Patients. Circ Genom Precis Med. 2018;11(2):e002020.
Marrani E, Burns JC, Cimaz R. How should we classify Kawasaki Disease? Front Immunol. 2018;9:2974.
Lin MT, Wang JK, Yeh JI, Sun LC, Chen PL, Wu JF, Chang CC, Lee WL, Shen CT, Wang NK, et al. Clinical implication of the C allele of the ITPKC Gene SNP rs28493229 in Kawasaki Disease: Association with Disease susceptibility and BCG scar reactivation. Pediatr Infect Dis J. 2011;30(2):148–52.
Huang YH, Chen KD, Lo MH, Cai XY, Chang LS, Kuo YH, Huang WD, Kuo HC. Decreased DNA methyltransferases expression is associated with coronary artery lesion formation in Kawasaki disease. Int J Med Sci. 2019;16(4):576–82.
McCrindle BW, Rowley AH, Newburger JW, Burns JC, Bolger AF, Gewitz M, Baker AL, Jackson MA, Takahashi M, Shah PB, 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(17):e927–99.
Kuo HC, Chang WC. Genetic polymorphisms in Kawasaki disease. Acta Pharmacol Sin. 2011;32(10):1193–8.
Onouchi Y. The genetics of Kawasaki disease. Int J Rheum Dis. 2018;21(1):26–30.
Chen MR, Kuo HC, Lee YJ, Chi H, Li SC, Lee HC, Yang KD. Phenotype, susceptibility, autoimmunity, and Immunotherapy between Kawasaki Disease and Coronavirus Disease-19 Associated Multisystem Inflammatory Syndrome in Children. Front Immunol. 2021;12:632890.
Chen SY, Wan L, Huang YC, Sheu JJ, Lan YC, Lai CH, Lin CW, Chang JS, Tsai Y, Liu SP, et al. Interleukin-18 gene 105A/C genetic polymorphism is associated with the susceptibility of Kawasaki disease. J Clin Lab Anal. 2009;23(2):71–6.
Fu LY, Qiu X, Deng QL, Huang P, Pi L, Xu Y, Che D, Zhou H, Lu Z, Tan Y, et al. The IL-1B gene polymorphisms rs16944 and rs1143627 contribute to an increased risk of coronary artery lesions in Southern Chinese Children with Kawasaki Disease. J Immunol Res. 2019;2019:4730507.
Kim KY, Bae YS, Ji W, Shin D, Kim HS, Kim DS. ITPKC and SLC11A1 gene polymorphisms and gene-gene interactions in korean patients with Kawasaki Disease. Yonsei Med J. 2018;59(1):119–27.
Peng Q, Chen C, Zhang Y, He H, Wu Q, Liao J, Li B, Luo C, Hu X, Zheng Z, et al. 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(7):1046–53.
Khor CC, Davila S, Breunis WB, Lee YC, Shimizu C, Wright VJ, Yeung RS, Tan DE, Sim KS, Wang JJ, et al. Genome-wide association study identifies FCGR2A as a susceptibility locus for Kawasaki disease. Nat Genet. 2011;43(12):1241–6.
Fujita Y, Nakamura Y, Sakata K, Hara N, Kobayashi M, Nagai M, Yanagawa H, Kawasaki T. Kawasaki disease in families. Pediatrics. 1989;84(4):666–9.
Hata A, Onouchi Y. Susceptibility genes for Kawasaki disease: toward implementation of personalized medicine. J Hum Genet. 2009;54(2):67–73.
Hsu JSJ, So M, Tang CSM, Karim A, Porsch RM, Wong C, Yu M, Yeung F, Xia H, Zhang R, et al. De novo mutations in Caudal Type Homeo Box transcription factor 2 (CDX2) in patients with persistent cloaca. Hum Mol Genet. 2018;27(2):351–8.
Kim JJ, Hong YM, Yun SW, Lee KY, Yoon KL, Han MK, Kim GB, Kil HR, Song MS, Lee HD, et al. Identification of rare coding variants associated with Kawasaki disease by whole exome sequencing. Genomics Inf. 2021;19(4):e38.
Lo MS. A framework for understanding Kawasaki disease pathogenesis. Clin Immunol. 2020;214:108385.
Kircher M, Witten DM, Jain P, O’Roak BJ, Cooper GM, Shendure J. A general framework for estimating the relative pathogenicity of human genetic variants. Nat Genet. 2014;46(3):310–5.
de Graeff N, Groot N, Ozen S, Eleftheriou D, Avcin T, Bader-Meunier B, Dolezalova P, Feldman BM, Kone-Paut I, Lahdenne P, et al. European consensus-based recommendations for the diagnosis and treatment of Kawasaki disease - the SHARE initiative. Rheumatology (Oxford). 2019;58(4):672–82.
Rentzsch P, Witten D, Cooper GM, Shendure J, Kircher M. CADD: predicting the deleteriousness of variants throughout the human genome. Nucleic Acids Res. 2019;47(D1):D886–94.
Wang CL, Wu YT, Liu CA, Kuo HC, Yang KD. Kawasaki disease: infection, immunity and genetics. Pediatr Infect Dis J. 2005;24(11):998–1004.
Uehara R, Belay ED. Epidemiology of Kawasaki disease in Asia, Europe, and the United States. J Epidemiol. 2012;22(2):79–85.
Chen P, Miyake M, Fan Q, Liao J, Yamashiro K, Ikram MK, Chew M, Vithana EN, Khor CC, Aung T, et al. CMPK1 and RBP3 are associated with corneal curvature in asian populations. Hum Mol Genet. 2014;23(22):6129–36.
Bae Y, Shin D, Nam J, Lee HR, Kim JS, Kim KY, Kim DS, Chung YJ. Variants in the gene EBF2 are Associated with Kawasaki Disease in a korean Population. Yonsei Med J. 2018;59(4):519–23.
Dergun M, Kao A, Hauger SB, Newburger JW, Burns JC. Familial occurrence of Kawasaki syndrome in North America. Arch Pediatr Adolesc Med. 2005;159(9):876–81.
Kwon YC, Kim JJ, Yun SW, Yu JJ, Yoon KL, Lee KY, Kil HR, Kim GB, Han MK, Song MS, et al. Male-specific association of the FCGR2A His167Arg polymorphism with Kawasaki disease. PLoS ONE. 2017;12(9):e0184248.
Buda P, Chyb M, Smorczewska-Kiljan A, Wieteska-Klimczak A, Paczesna A, Kowalczyk-Domagala M, Okarska-Napierala M, Sobalska-Kwapis M, Grochowalski L, Slomka M, et al. Association between rs12037447, rs146732504, rs151078858, rs55723436, and rs6094136 polymorphisms and Kawasaki Disease in the Population of Polish Children. Front Pead. 2021;9:624798.
Chen MR, Chang TY, Chiu NC, Chi H, Yang KD, Chang L, Huang DT, Huang FY, Lien YP, Lin WS, et al. Validation of genome-wide associated variants for Kawasaki disease in a taiwanese case-control sample. Sci Rep. 2020;10(1):11756.
Lee YC, Kuo HC, Chang JS, Chang LY, Huang LM, Chen MR, Liang CD, Chi H, Huang FY, Lee ML, et al. Two new susceptibility loci for Kawasaki disease identified through genome-wide association analysis. Nat Genet. 2012;44(5):522–5.
Kim JJ, Yun SW, Yu JJ, Yoon KL, Lee KY, Kil HR, Kim GB, Han MK, Song MS, Lee HD, et al. A genome-wide association analysis identifies NMNAT2 and HCP5 as susceptibility loci for Kawasaki disease. J Hum Genet. 2017;62(12):1023–9.
Zeng S, Zhang T, Madigan MC, Fernando N, Aggio-Bruce R, Zhou F, Pierce M, Chen Y, Huang L, Natoli R, et al. Interphotoreceptor Retinoid-Binding protein (IRBP) in Retinal Health and Disease. Front Cell Neurosci. 2020;14:577935.
Arno G, Hull S, Robson AG, Holder GE, Cheetham ME, Webster AR, Plagnol V, Moore AT. Lack of Interphotoreceptor Retinoid binding protein caused by homozygous mutation of RBP3 is Associated with high myopia and retinal dystrophy. Invest Ophthalmol Vis Sci. 2015;56(4):2358–65.
Yokomizo H, Maeda Y, Park K, Clermont AC, Hernandez SL, Fickweiler W, Li Q, Wang CH, Paniagua SM, Simao F, et al. Retinol binding protein 3 is increased in the retina of patients with diabetes resistant to diabetic retinopathy. Sci Transl Med. 2019;11(499):eaau6627.
Wu Q, Blakeley LR, Cornwall MC, Crouch RK, Wiggert BN, Koutalos Y. Interphotoreceptor retinoid-binding protein is the physiologically relevant carrier that removes retinol from rod photoreceptor outer segments. Biochemistry. 2007;46(29):8669–79.
Gonzalez-Fernandez F. Interphotoreceptor retinoid-binding protein–an old gene for new eyes. Vis Res. 2003;43(28):3021–36.
Golomb E, Ma X, Jana SS, Preston YA, Kawamoto S, Shoham NG, Goldin E, Conti MA, Sellers JR, Adelstein RS. Identification and characterization of nonmuscle myosin II-C, a new member of the myosin II family. J Biol Chem. 2004;279(4):2800–8.
Zhu Z, Peng L, Chen G, Jiang W, Shen Z, Du C, Zang R, Su Y, Xie H, Li H, et al. Mutations of MYH14 are associated to anorectal malformations with recto-perineal fistulas in a small subset of chinese population. Clin Genet. 2017;92(5):503–9.
Wang M, Zhou Y, Zhang F, Fan Z, Bai X, Wang H. A novel MYH14 mutation in a chinese family with autosomal dominant nonsyndromic hearing loss. BMC Med Genet. 2020;21(1):154.
Aggarwal V, Etinger V, Orjuela AF. Sensorineural hearing loss in Kawasaki disease. Ann Pediatr Cardiol. 2016;9(1):87–9.
We thank all of our patients and their families for participating in this study. The authors thank Huifeng Han for her contribution to this study.
This work was supported by a grant from Kunming “Spring City Plan” High-level talent introduced by the Engineering Young Talents Special Project and a fund project of Yunnan Province Clinical Research Center for Children’s Health and Disease(2022-ETYY-YJ-16).
Ethics approval and consent to participate
This study was approved by the Ethical Committee, Kunming Children’s Hospital, Yunnan Province under the number 2022/722.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
About this article
Cite this article
Zhang, X., Sun, Y., Meng, L. et al. Whole-exome sequencing analysis identifies novel variants associated with Kawasaki disease susceptibility. Pediatr Rheumatol 21, 78 (2023). https://doi.org/10.1186/s12969-023-00857-0