Screening Key Differentially Expressed Genes in Kawasaki Disease via Integrated Analysis

authors:

avatar Keliang Li 1 , avatar Yan Jiao 2 , avatar Jingjing Liang 1 , * , avatar Pingping Pan 3 , avatar Yanji Zhu 1 , avatar Yilei Li 4

Department of Pediatrics, People's Hospital of Rizhao, Rizhao, China
Department of Neonatology, People's Hospital of Rizhao, Rizhao, China
Department of Pediatrics, Wulian County Rehabilitation Hospital, Rizhao, China
Department of Clinical Laboratory, People's Hospital of Rizhao, Rizhao, China
Iranian Journal of Pediatrics: Vol.31, issue 6; e114730
published online: November 29, 2021
article type: Research Article
received: April 22, 2021
revised: October 18, 2021
accepted: October 20, 2021
DOI: https://doi.org/10.5812/ijp.114730

how to cite: Li K , Jiao Y , Liang J, Pan P , Zhu Y , et al. Screening Key Differentially Expressed Genes in Kawasaki Disease via Integrated Analysis. Iran J Pediatr. 2021;31(6):e114730. https://doi.org/10.5812/ijp.114730.

Abstract

Background:

The current study was done to identify key genes associated with Kawasaki disease (KD).

Methods:

Three datasets were collected from Gene Expression Omnibus (GEO) database. Then, differentially expressed genes (DEGs) analysis, gene ontology (GO) annotation, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment were performed. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to determine the expression levels of DEGs in KD. Receiver operating characteristic (ROC) analysis was performed to assess the diagnostic value of DEGs.

Results:

In total, 2923 DEGs (1239 up- and 1684 down-regulated genes) were detected in KD. Ribosome, Leishmaniasis, and Tuberculosis significantly enriched KEGG pathways of DEGs. Six DEGs, including ADM, S100A12, ZNF438, MYD88, FCGR2A, and FCGR3B, were selected for qRT-PCR validation. Except for MYD88, the qRT-PCR results displayed similar expression patterns with that in our integrated analysis. ROC analysis revealed the diagnostic value of the six DEGs.

Conclusions:

Our study was expected to provide clues toward understanding the pathophysiology of KD inflammation.

1. Background

Kawasaki disease (KD) is an acute multisystem vasculitis, preferentially affecting coronary arteries that mainly occurs in infants and young children (1). In asset-rich countries, it is the leading cause of pediatric acquired heart disease (2). The Japanese pediatrician Tomisaku Kawasaki first described KD in 1967; however, both the cause and pathogenesis of which remain unclear (3). Intravenous immunoglobulin (IVIG) is used to treat a wide variety of rheumatologic and other autoimmune disorders, including KD; however, the mechanism of action is not fully understood (4). Although timely treatment with high-dose IVIG for patients with acute stage KD shortens the duration of fever and reduces the incidence of coronary artery lesions, approximately 10% of the patients fail to respond to IVIG treatment (5, 6). The self-limited nature is unique in KD, and the most patients recover well (7). However, KD carries significant clinical importance and needs special attention because it can lead to serious cardiovascular complications.

Advance development in the investigation on KD includes searching for genetic susceptibility correlated with KD and a study on immunopathogenesis based on the review, in which it was reported that innate immunity performs a more prominent role than the acquired immunity. With the advent of genome-wide association studies (GWAS), a number of genes, including ITPKC, CASP3, BLK, CD40, and HLA were identified as candidates of KD-susceptibility genes (7). With emerging of various high-throughput technologies, a huge amount of RNA-seq data have been available, and numerous genes have been identified to be associated with KD. The comprehensive analysis of multiple gene expression profile data from different platforms will help to explore the pathogenesis of KD with a larger sample size more accurately than a single gene expression profile.

In the attempt to provide a more comprehensive analysis, we performed an integrated analysis based on multiple microarray datasets of KD collected from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) in KD were acquired. In addition, functional annotation of DEGs was conducted. Our study is expected to represent a new avenue to explore the specific mechanism for KD.

2. Methods

2.1. Datasets Collection

The mRNA expression profiles of KD were acquired from the GEO database. From 2016 to 2018, datasets meeting the following criteria were included in our study: (1) selected datasets should be whole-genome mRNA expression profile by array; (2) data derived from blood samples of patients with KD and normal controls; and (3) normalized or original datasets. After filtering, three datasets, including GSE73464, GSE109351, and GSE68004, were included in this study (Appendix 1).

2.2. Screening of DEGs Between Patients with KD and Normal Controls and Functional Annotation

MetaMA was utilized to acquire the DEGs. The threshold for the significance of DEGs was the false discovery rate (FDR) < 0.01 and |Combined.ES| > 1. Hierarchical clustering was visualized using the R package “pheatmap”. CPDB (http://cpdb.molgen.mpg.de/CPDB) was used to perform gene ontology (GO) classification and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The threshold was set as FDR < 0.01.

2.3. Quantitative Real-time Polymerase Chain Reaction (qRT-PCR) Validation

Twelve blood samples were collected from five healthy individuals and seven patients with KD (Table 1). All patients were first diagnosed with KD and during the acute phase. Patients with infection of Epstein-Barr virus or other detectable pathogens and organ dysfunction caused by other reasons were excluded. The healthy control group underwent regular health examinations, had no fever, and no immune-related diseases, two weeks before and after sampling. Children with immunosuppressive therapy or fever in the past two weeks were excluded. All participants or their legally authorized representatives provided written informed consent to participate in the study, and the ethics committee of our hospital approved this study protocol.

Table 1.

Clinical Features of Patients with KD and Control Patients

SexAge (mo)WBC (109/L)Neutrophils (%)Platelets (109/L)CRP (mg/L)Duration of Fever (days)
Case 1Female2612.8670.322044.37
Case 2Male2521.7472.340852.47
Case 3Male712.7956.619853.87
Case 4Male817.8886.343361.96
Case 5Male613.6359.2435142.97
Case 6Male1618.6481.730145.65
Case 7Female7324.8295.9177157.910
Control 1Male15
Control 2Male16
Control 3Female60
Control 4Female24
Control 5Male10

Abbreviations: KD, Kawasaki disease; WBC, white blood cells; CRP, C-reactive protein.

Total RNA was extracted with the Trizol reagent (Invitrogen, USA). The qRT-PCR was performed on SuperReal PreMix Plus (Invitrogen, USA) in ABI 7500 real-time PCR detection system. Using the 2-ΔΔCt method, relative gene expression was determined. The human ACTB was served as a control in analysis. The PCR primers are shown in Appendix 2.

2.4. Receiver Operating Characteristic Analysis

Receiver operating characteristic (ROC) analysis is an effective method to evaluate the performance of diagnostic tests, which was used to study the diagnostic capacity of DEGs. The area under the curve (AUC) was further calculated.

2.5. Validation in GEO Database

GSE18606, consisting of 48 whole blood samples from nine healthy age-appropriate controls, eight IVIG non-responding, and 12 IVIG-responding patients at acute and convalescent stages, and GSE63881, consisting of the acute and convalescent whole blood samples from 171 KD subjects, were downloaded from the GEO database. We used GSE18606 in two steps: in the first step to validate the expression of DEGs regarding eight IVIG non-responding and 12 IVIG-responding patients at acute and convalescent stages as 20 KD patients, and in the second step, to explore the association between the expression of DEGs and the stage of KD.

3. Results

3.1. Identification of DEGs in KD

Compared with normal controls, a total of 2923 DEGs (1239 up- and 1684 down-regulated genes) were detected in KD, of which H2AFJ and KLRG1 were the most up- and down-regulated genes, respectively (Table 2). Hierarchical clustering analysis of the DEGs is shown in Figure 1.

Table 2.

Top 20 Up- and Down-regulation DEGs in KD

IDSymbolCombined.ESP-ValueFDRRegulation
55766H2AFJ3.70001600up
199675C19orf593.6844200up
306ANXA33.61134600up
723790HIST2H2AA43.58803300up
6583SLC22A43.39151400up
7100TLR53.38952500up
8291DYSF3.34603100up
133ADM3.32455600up
249ALPL3.28028200up
4001LMNB13.21987900up
5836PYGL3.20407100up
6283S100A123.20283600up
8349HIST2H2BE3.18324500up
11213IRAK33.17950300up
9334B4GALT53.16061300up
2992GYG13.080600up
100132417FCGR1C3.06723500up
220929ZNF4383.06354700up
3101HK33.05312100up
604BCL63.03520700up
10219KLRG1-2.6413700down
3003GZMK-2.5586100down
10425ARIH2-2.5311200down
2205FCER1A-2.5071700down
3122HLA-DRA-2.5039500down
9404LPXN-2.4483800down
257101ZNF683-2.4242200down
83988NCALD-2.4110200down
4242MFNG-2.404500down
914CD2-2.3453100down
445347TARP-2.3451800down
6136RPL12-2.3426300down
3113HLA-DPA1-2.3034200down
353APRT-2.2953800down
3587IL10RA-2.2915800down
9806SPOCK2-2.2837900down
917CD3G-2.2674600down
6138RPL15-2.2628400down
11184MAP4K1-2.2590100down
925CD8A-2.2441800down

Abbreviations: DEGs, differentially expressed genes; KD, Kawasaki disease; ES, effect size; FDR, false discovery rate.

Hierarchical clustering of top 100 up- and down-regulated differentially expressed genes (DEGs) in Kawasaki disease (KD).
Hierarchical clustering of top 100 up- and down-regulated differentially expressed genes (DEGs) in Kawasaki disease (KD).
Figure 1.

Hierarchical clustering of top 100 up- and down-regulated differentially expressed genes (DEGs) in Kawasaki disease (KD).

3.2. Functional Annotation of DEGs

GO analysis indicated that these DEGs were mainly associated with leukocyte activation (FDR = 1.03E-54) and cell activation (FDR = 1.91E-54) involved in immune response, cytoplasmic part (FDR = 8.78E-58), and protein binding (FDR = 3.52E-32) (Figure 2A-C). The most significantly enriched KEGG pathway were ribosome (FDR = 4.62E-13), leishmaniasis (FDR = 1.44E-07) and tuberculosis (FDR = 1.83E-07) (Figure 2D).

Significantly enriched gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of differentially expressed genes (DEGs) in Kawasaki disease (KD).
Significantly enriched gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of differentially expressed genes (DEGs) in Kawasaki disease (KD).
Figure 2.

Significantly enriched gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of differentially expressed genes (DEGs) in Kawasaki disease (KD).

3.3. QRT-PCR Validation

The expression levels of six DEGs, including ADM, S100A12, ZNF438, MYD88, FCGR2A, and FCGR3B, were determined by qRT-PCR. In this analysis, all these six DEGs were up-regulated in KD. The expression levels of these six DEGs in GEO are displayed in Appendix 3. As shown in Figure 3, the qRT-PCR results for all the DEGs except for MYD88 were generally consistent with that in our integrated analysis (Figure 3).

The quantitative real-time polymerase chain reaction (qRT-PCR) results of differentially expressed genes (DEGs) in Kawasaki disease (KD).
The quantitative real-time polymerase chain reaction (qRT-PCR) results of differentially expressed genes (DEGs) in Kawasaki disease (KD).
Figure 3.

The quantitative real-time polymerase chain reaction (qRT-PCR) results of differentially expressed genes (DEGs) in Kawasaki disease (KD).

3.4. ROC Analysis

The AUC of these six DEGs, including ADM (0.937), S100A12 (0.956), ZNF438 (0.922), MYD88 (0.882), FCGR2A (0.877), and FCGR3B (0.861), was more than 0.85, indicating that these six DEGs may be potential diagnostic biomarkers for KD (Figure 4).

The receiver operating characteristic (ROC) curves of differentially expressed genes (DEGs) in Kawasaki disease (KD). A) ADM, B) S100A12, C) ZNF438, D) MYD88, E) FCGR2A, and F) FCGR3B.
The receiver operating characteristic (ROC) curves of differentially expressed genes (DEGs) in Kawasaki disease (KD). A) ADM, B) S100A12, C) ZNF438, D) MYD88, E) FCGR2A, and F) FCGR3B.
Figure 4.

The receiver operating characteristic (ROC) curves of differentially expressed genes (DEGs) in Kawasaki disease (KD). A) ADM, B) S100A12, C) ZNF438, D) MYD88, E) FCGR2A, and F) FCGR3B.

3.5. Validation in GEO Database

The expression levels of six DEGs, including ADM, S100A12, ZNF438, MYD88, FCGR2A, and FCGR3B, were verified using the GSE18606 dataset. As shown in Figure 5, these six DEGs exhibited the same patterns with our integrated analysis, and except ADM (0.684) and FCGR2A (0.695), the AUC of these DEGs, including S100A12 (0.811), ZNF438 (0.758), MYD88 (0.843), and FCGR3B (0.775) was more than 0.7. The expression levels of these six DEGs in different stages of KD based on GSE18606 and GSE63881 are displayed in Appendix 4. The expression levels of these DEGs exhibited a significant increase in the acute phase compared to healthy controls in GSE18606 and exhibited a significant decrease in the convalescent phase compared to the acute phase in GSE18606 and GSE63881. In other words, the expression levels of these DEGs decreased with the remission of KD symptoms.

Validation in GSE18606. A-F, box-plots; A) ADM, B) S100A12, C) ZNF438, D) MYD88, E) FCGR2A, F) FCGR3B; A'-F', ROC curves; A') ADM, B') S100A12, C') ZNF438, D') MYD88, E') FCGR2A, and F') FCGR3B.
Validation in GSE18606. A-F, box-plots; A) ADM, B) S100A12, C) ZNF438, D) MYD88, E) FCGR2A, F) FCGR3B; A'-F', ROC curves; A') ADM, B') S100A12, C') ZNF438, D') MYD88, E') FCGR2A, and F') FCGR3B.
Figure 5.

Validation in GSE18606. A-F, box-plots; A) ADM, B) S100A12, C) ZNF438, D) MYD88, E) FCGR2A, F) FCGR3B; A'-F', ROC curves; A') ADM, B') S100A12, C') ZNF438, D') MYD88, E') FCGR2A, and F') FCGR3B.

4. Discussion

KD is a systemic vasculitis syndrome of unclear etiology. Although the clinical features, diagnosis, and treatment of KD are well determined, its pathogenesis remains unclear (8). In this study, a series of bioinformatics analyses were performed to investigate crucial genes correlated with KD.

Adrenomedullin (ADM) is a potent vasodilator peptide and is involved in the dilation of coronary arteries in acute KD (9, 10). It has been reported that ADM can significantly reduce cardiac contractility and induce the expression of adhesion molecules on human vascular endothelial cells and during endothelial cell proliferation, which may be one of the causes of acute KD-induced vascular inflammation (11-13). In an immunoradiometric assay, remarkably higher plasma ADM levels were detected in acute KD (14). Based on oligonucleotide microarray, RT-PCR, and ELISA, Nomura et al. indicated an increase in ADM levels in acute KD, as well (10). In addition, down-regulated ADM was observed in KD patients after intravenous infusion of high-dose IVIG therapy (15). In agreement with previous studies, ADM was significantly up-regulated in our analysis, which confirmed the critical role of ADM in KD.

S100 calcium-binding protein A12 (S100A12), a calcium-binding protein with proinflammatory properties, directly activates endothelial cells, macrophages, and lymphocytes by interacting with the multiligand receptor for advanced glycation end products (RAGE) (16). Many studies have demonstrated that S100A12 is closely related to KD. S100A12 was reported to function as a mediator in murine inflammatory disorders (17). Foell et al. observed that the serum concentrations of S100A12 increased in the acute phase of KD and subsequently decreased in response to high-dose IVIG therapy, indicating that S100A12 could be used as an additional serum marker for KD (18). The significantly up-regulated S100A12 in peripheral blood mononuclear cells (PBMCs) was observed at the acute phase of KD (19). Consistent with previous reports, significantly up-regulated S100A12 was detected in the current analysis, which indicated the important role of S100A12 in the development of KD.

To the best of our knowledge, there were no previous studies linking ZNF438 to KD. ZNF438 is a C2H2 type zinc finger gene situated on human chromosome 10p11.2. It is widely expressed in all human adult tissues and reported to possess transcription inhibition activity and may function as a transcription factor (20). ZNF438 showed an increasing trend in KD in our study, which may suggest that elevated ZNF438 may play an important role in KD.

MYD88 encodes a cytosolic adapter protein that exerts a critical role in the innate and adaptive immune response, which serves as an essential signal transducer for inflammatory signaling pathways, downstream of members of the Toll-like receptor and interleukin-1 receptor families. It was reported that MYD88 contributed to the activation of proinflammatory innate immune mechanisms in the mouse model of coronary arteritis (21). MYD88 was elevated in the PBMCs of untreated KD patients at the acute phase and significantly decreased after IVIG therapy (22). In this integrated analysis, MYD88 was up-regulated in patients with KD, while it was down-regulated in our qRT-PCR results. Hence, more studies are required to investigate the precise roles of MYD88 in KD.

The binding of Fc gamma receptors (FCGRs) to the Fc region of IgG is an important immune response mechanism, and then, this process may induce biological responses, such as phagocytosis and immune complexes (ICs) processing (23). There are five low-affinity genes encoding FCGRs in humans, including FCGR2A, FCGR2B, FCGR2C, FCGR3A, and FCGR3B (24). FCGR2A is usually expressed on immune-responsive cells, such as macrophages, monocytes, neutrophils, and dendritic cells to promote phagocytosis and the production of inflammatory mediators (25, 26). A meta-analysis on the Asians linked the H131R polymorphism in the FCGR2A gene to KD susceptibility (27). A genome-wide association study (GWAS) detected FCGR2A as a susceptibility locus for KD (28). Besides, FCGR3B has been associated with non-response to IVIG and risk of coronary artery aneurysms (CAA) (29). Both FCGR2A and FCGR3B were identified as significantly up-regulated DEGs in this study, which may indicate that these two genes play momentous roles in the pathophysiology of KD.

In addition, the expression levels of these six DEGs in GSE18606 and GSE63881 revealed that elevated levels of ADM, S100A12, ZNF438, MYD88, FCGR2A, and FCGR3B in the acute phase of KD significantly decreased in the convalescent phase. This may further indicate the important roles of these DEGs in KD, which can be served as potential therapeutic targets of KD. One of the limitations of the present study was the small sample size for validation. Hence, further experimental studies on a larger sample size are needed to confirm these results.

4.1. Conclusions

In conclusion, 2923 DEGs (1239 up- and 1684 down-regulated genes) were detected in KD, and the importance and diagnostic value of several DEGs, which may implicate in KD were highlighted. This study is expected to provide clues toward understanding the pathophysiology of KD inflammation. In addition, the closed relationship between the expressions of key DEGs and the stage of KD needs to be further explored in a large sample containing different stages, including follow-up studies. Also, it would be interesting to investigate the effect of some factors, such as gender, age, ethnicity, and outcome, on the DEGs in the future.

References

  • 1.

    Wittkowski H, Hirono K, Ichida F, Vogl T, Ye F, Yanlin X, et al. Acute Kawasaki disease is associated with reverse regulation of soluble receptor for advance glycation end products and its proinflammatory ligand S100A12. Arthritis Rheum. 2007;56(12):4174-81. [PubMed ID: 18050248]. https://doi.org/10.1002/art.23042.

  • 2.

    Armaroli G, Verweyen E, Pretzer C, Kessel K, Hirono K, Ichida F, et al. Monocyte-Derived Interleukin-1beta As the Driver of S100A12-Induced Sterile Inflammatory Activation of Human Coronary Artery Endothelial Cells: Implications for the Pathogenesis of Kawasaki Disease. Arthritis Rheumatol. 2019;71(5):792-804. [PubMed ID: 30447136]. https://doi.org/10.1002/art.40784.

  • 3.

    Kawasaki T. [Acute febrile mucocutaneous syndrome with lymphoid involvement with specific desquamation of the fingers and toes in children]. Arerugi. 1967;16(3):178-222. [PubMed ID: 6062087].

  • 4.

    Rodriguez MM, Wagner-Weiner L. Intravenous Immunoglobulin in Pediatric Rheumatology: When to Use It and What Is the Evidence. Pediatr Ann. 2017;46(1):e19-24. [PubMed ID: 28079914]. https://doi.org/10.3928/19382359-20161214-01.

  • 5.

    Leung DY, Schlievert PM, Meissner HC. The immunopathogenesis and management of Kawasaki syndrome. Arthritis Rheum. 1998;41(9):1538-47. [PubMed ID: 9751085]. https://doi.org/10.1002/1529-0131(199809)41:9<1538::AID-ART3>3.0.CO;2-M.

  • 6.

    Nakamura Y, Yanagawa H, Harada K, Kato H, Kawasaki T. Mortality among persons with a history of Kawasaki disease in Japan: existence of cardiac sequelae elevated the mortality. J Epidemiol. 2000;10(6):372-5. [PubMed ID: 11210105]. https://doi.org/10.2188/jea.10.372.

  • 7.

    Kim KY, Kim DS. Recent Advances in Kawasaki Disease. Yonsei Med J. 2016;57(1):15-21. [PubMed ID: 26632378]. [PubMed Central ID: PMC4696947]. https://doi.org/10.3349/ymj.2016.57.1.15.

  • 8.

    Chen Y, Ding YY, Ren Y, Cao L, Xu QQ, Sun L, et al. Identification of differentially expressed microRNAs in acute Kawasaki disease. Mol Med Rep. 2018;17(1):932-8. [PubMed ID: 29115644]. [PubMed Central ID: PMC5780174]. https://doi.org/10.3892/mmr.2017.8016.

  • 9.

    Kitamura K, Kangawa K, Kawamoto M, Ichiki Y, Nakamura S, Matsuo H, et al. Adrenomedullin: a novel hypotensive peptide isolated from human pheochromocytoma. 1993. Biochem Biophys Res Commun. 2012;425(3):548-55. [PubMed ID: 22925672]. https://doi.org/10.1016/j.bbrc.2012.08.022.

  • 10.

    Nomura I, Abe J, Noma S, Saito H, Gao B, Wheeler G, et al. Adrenomedullin is highly expressed in blood monocytes associated with acute Kawasaki disease: a microarray gene expression study. Pediatr Res. 2005;57(1):49-55. [PubMed ID: 15531734]. https://doi.org/10.1203/01.PDR.0000147745.52711.DD.

  • 11.

    Kaygisiz Z, Erksap N, Uyar R, Kabadere S, Kabadere TE, Dernek S. The effect of adrenomedullin, amylin fragment 8-37 and calcitonin gene-related peptide on contractile force, heart rate and coronary perfusion pressure in isolated rat hearts. Acta Physiol Hung. 2003;90(2):133-46. [PubMed ID: 12903912]. https://doi.org/10.1556/APhysiol.90.2003.2.6.

  • 12.

    Hagi-Pavli E, Farthing PM, Kapas S. Stimulation of adhesion molecule expression in human endothelial cells (HUVEC) by adrenomedullin and corticotrophin. Am J Physiol Cell Physiol. 2004;286(2):C239-46. [PubMed ID: 14534081]. https://doi.org/10.1152/ajpcell.00036.2003.

  • 13.

    Miyashita K, Itoh H, Sawada N, Fukunaga Y, Sone M, Yamahara K, et al. Adrenomedullin promotes proliferation and migration of cultured endothelial cells. Hypertens Res. 2003;26 Suppl:S93-8. [PubMed ID: 12630817]. https://doi.org/10.1291/hypres.26.s93.

  • 14.

    Nishida K, Watanabe K, Echigo S, Mayumi M, Nishikimi T. Increased plasma adrenomedullin levels in Kawasaki disease with coronary artery involvement. Am J Med. 2001;111(2):165-6. [PubMed ID: 11501552]. https://doi.org/10.1016/s0002-9343(01)00781-1.

  • 15.

    Abe J, Jibiki T, Noma S, Nakajima T, Saito H, Terai M. Gene expression profiling of the effect of high-dose intravenous Ig in patients with Kawasaki disease. J Immunol. 2005;174(9):5837-45. [PubMed ID: 15843588]. https://doi.org/10.4049/jimmunol.174.9.5837.

  • 16.

    Hofmann MA, Drury S, Fu C, Qu W, Taguchi A, Lu Y, et al. RAGE mediates a novel proinflammatory axis: a central cell surface receptor for S100/calgranulin polypeptides. Cell. 1999;97(7):889-901. [PubMed ID: 10399917]. https://doi.org/10.1016/s0092-8674(00)80801-6.

  • 17.

    Schmidt AM, Yan SD, Yan SF, Stern DM. The multiligand receptor RAGE as a progression factor amplifying immune and inflammatory responses. J Clin Invest. 2001;108(7):949-55. [PubMed ID: 11581294]. [PubMed Central ID: PMC200958]. https://doi.org/10.1172/JCI14002.

  • 18.

    Foell D, Ichida F, Vogl T, Yu X, Chen R, Miyawaki T, et al. S100A12 (EN-RAGE) in monitoring Kawasaki disease. Lancet. 2003;361(9365):1270-2. [PubMed ID: 12699958]. https://doi.org/10.1016/S0140-6736(03)12986-8.

  • 19.

    Ikeda K, Yamaguchi K, Tanaka T, Mizuno Y, Hijikata A, Ohara O, et al. Unique activation status of peripheral blood mononuclear cells at acute phase of Kawasaki disease. Clin Exp Immunol. 2010;160(2):246-55. [PubMed ID: 20015095]. [PubMed Central ID: PMC2857948]. https://doi.org/10.1111/j.1365-2249.2009.04073.x.

  • 20.

    Zhong Z, Wan B, Qiu Y, Ni J, Tang W, Chen X, et al. Identification of a novel human zinc finger gene, ZNF438, with transcription inhibition activity. J Biochem Mol Biol. 2007;40(4):517-24. [PubMed ID: 17669267]. https://doi.org/10.5483/bmbrep.2007.40.4.517.

  • 21.

    Rosenkranz ME, Schulte DJ, Agle LM, Wong MH, Zhang W, Ivashkiv L, et al. TLR2 and MyD88 contribute to Lactobacillus casei extract-induced focal coronary arteritis in a mouse model of Kawasaki disease. Circulation. 2005;112(19):2966-73. [PubMed ID: 16275884]. https://doi.org/10.1161/CIRCULATIONAHA.105.537530.

  • 22.

    Mortazavi SH, Amin R, Alyasin S, Kashef S, Karimi MH, Babaei M, et al. Down-regulation of TLR2, 3, 9 and Signaling Mediators, MyD88 and TRIF, Gene Transcript Levels in Patients with Kawasaki Disease Treated with IVIG. Iran J Allergy Asthma Immunol. 2015;14(2):188-97. [PubMed ID: 25780885].

  • 23.

    Ravetch JV, Bolland S. IgG Fc receptors. Annu Rev Immunol. 2001;19:275-90. [PubMed ID: 11244038]. https://doi.org/10.1146/annurev.immunol.19.1.275.

  • 24.

    Brown EE, Edberg JC, Kimberly RP. Fc receptor genes and the systemic lupus erythematosus diathesis. Autoimmunity. 2007;40(8):567-81. [PubMed ID: 18075791]. https://doi.org/10.1080/08916930701763710.

  • 25.

    Pleass RJ, Woof JM. Fc receptors and immunity to parasites. Trends Parasitol. 2001;17(11):545-51. [PubMed ID: 11872400]. https://doi.org/10.1016/s1471-4922(01)02086-4.

  • 26.

    Gerber JS, Mosser DM. Stimulatory and inhibitory signals originating from the macrophage Fcgamma receptors. Microbes Infect. 2001;3(2):131-9. [PubMed ID: 11251299]. https://doi.org/10.1016/s1286-4579(00)01360-5.

  • 27.

    Lin L, Wang SY, Yang SB, Xiao FC. Association between the FCGR2A gene H131R polymorphism and risk of Kawasaki disease: a meta-analysis. Genet Mol Res. 2015;14(2):6256-64. [PubMed ID: 26125827]. https://doi.org/10.4238/2015.June.9.12.

  • 28.

    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(12):1241-6. [PubMed ID: 22081228]. https://doi.org/10.1038/ng.981.

  • 29.

    Eleftheriou D, Levin M, Shingadia D, Tulloh R, Klein NJ, Brogan PA. Management of Kawasaki disease. Arch Dis Child. 2014;99(1):74-83. [PubMed ID: 24162006]. [PubMed Central ID: PMC3888612]. https://doi.org/10.1136/archdischild-2012-302841.