Translational Bioinformatics and Computational Systems MedicineView this Special Issue
Understanding the Pathogenesis of Kawasaki Disease by Network and Pathway Analysis
Kawasaki disease (KD) is a complex disease, leading to the damage of multisystems. The pathogen that triggers this sophisticated disease is still unknown since it was first reported in 1967. To increase our knowledge on the effects of genes in KD, we extracted statistically significant genes so far associated with this mysterious illness from candidate gene studies and genome-wide association studies. These genes contributed to susceptibility to KD, coronary artery lesions, resistance to initial IVIG treatment, incomplete KD, and so on. Gene ontology category and pathways were analyzed for relationships among these statistically significant genes. These genes were represented in a variety of functional categories, including immune response, inflammatory response, and cellular calcium ion homeostasis. They were mainly enriched in the pathway of immune response. We further highlighted the compelling immune pathway of NF-AT signal and leukocyte interactions combined with another transcription factor NF-κB in the pathogenesis of KD. STRING analysis, a network analysis focusing on protein interactions, validated close contact between these genes and implied the importance of this pathway. This data will contribute to understanding pathogenesis of KD.
KD is a systemic vascular disease preferentially occurring in infants and children [1, 2]. It is characterized by the development of coronary artery aneurysms (CAA) which may result in fatal thrombosis and sudden cardiac failure. Clinical manifestations of KD include prolonged fever (1-2 weeks, mean 10-11 days), conjunctival infection, oral lesions, polymorphous skin rashes, extremity changes, and cervical lymphadenopathy, all of which comprise diagnostic criteria . However, great majority of children failed to manifest typical characteristics. In addition to the diagnostic criteria, there are a broad range of nonspecific clinical features, including irritability, uveitis, aseptic meningitis, cough, vomiting, diarrhea, abdominal pain, gallbladder hydrops, urethritis, arthralgia, arthritis, hypoalbuminemia , liver function impairment, and heart failure [5, 6]. The peaked incidence at 9 to 11 months of age coincides with fading of maternal immunity, and symptoms partly similar to other infectious disorders suggest that some microorganisms may trigger this disease. Despite great efforts to identify the cause for nearly a half a century, the etiology of KD still remains unknown . However, the role of genetic susceptibility to KD has long been evident through its striking predilection for children of Japanese ethnicity regardless of their country of residence; compared with Caucasian children, Japanese children have a relative risk of KD that is 10 to 15 times higher [8–10]. Siblings of KD children have a relative risk that is 6 to 10 times greater than that of children without a family history, and the parents of Japanese children with KD are twice as likely to have had KD themselves as children than other adults in the general Japanese population [11–14].
Candidate gene studies and genome-wide studies have been successively applied to explore the association between genetic effect and this mysterious disease [15, 16]. Many suspicious genes related to innate and acquired immune functions or to vascular remodeling have been studied [15, 17–19]. Genetic studies of KD were conducted not only to clarify the genetic background but also in the hope of providing clues about its etiology and pathogenesis. However, none of these studies have analyzed the internal association between these significant association genes and explored the possible pathogenic process in KD from overall level.
In this paper, we aim to extract statistically significant genes associated with KD (Up to September 2012, from all English databases) to explore their association and analyze their function in the pathogenesis of KD. This study is a systematic summary of previous research. Further studies on clinical validation will be summarized in our next study.
2. Methods and Materials
2.1. Extracting Genes with Statistical Significance
We performed a computerized search of Ovid, Google Scholar, and PubMed databases up to September 2012 and reviewed cited references to identify the relevant studies. Citations were screened at the title/abstract level and retrieved as full reports. Search keywords were “Kawasaki disease,” “Kawasaki syndrome,” “lymph node syndrome,” “mucocutaneous lymph node syndrome” combined with “polymorphism,” “gene,” “genetic,” “allele,” and “genotype”. The inclusion criteria of genes were those who have significant association with KD contributed to susceptibility, vascular lesions, resistance to initial IVIG treatment, late diagnosis of KD, and incomplete KD.
2.2. Data Analysis
DAVID (http://david.abcc.ncifcrf.gov/, version: 6.7) was used to process the bioinformatics analysis of these candidate gene markers, including gene classification (based on Biological Process Ontology and Molecular Function Ontology, resp.), enrichment analysis for significant gene ontology categories, KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway mapping, and significant pathway computing. GeneGo MetaCore (http://www.genego.com/, version: 6.5) was used to analyze the pathways of these significant genes. The association between these statistically significant genes were analyzed using STRING (http://string-db.org/), a database of known and predicted protein interactions.
3.1. Extracting Genes with Statistical Significance
3.2. Gene Ontology Analysis
Genes with statistical significance were submitted to functional analysis using DAVID software. Defense response, response to wounding, and inflammatory response were identified as significantly enriched (Enrichment Score = 15.91). Furthermore, DAVID analysis identified clusters of genes with annotations related to cellular calcium ion homeostasis, cell chemotaxis (enrichment Score: 3.75), and positive regulation of immune system process (Enrichment Score: 3.58) which is involved in autoimmune thyroid disease (hsa05320), asthma (hsa05310), type I diabetes mellitus (hsa04940), and allograft rejection (hsa04672). The functional annotation table can be available in supplementary material available online at http://dx.doi.org/10.1155/2013/989307.
3.3. Enrichment Analysis
Enrichment analysis consists of matching genes in functional ontologies by GeneGo MetaCore (Figure 1). The probability of a random intersection between a set of gene list with ontology entities was estimated with the “” value of the hypergeometric intersection. A lower “” value means higher relevance of the entity to the dataset, which appears in higher rating for the entity. All maps were drawn by GeneGo. The height of the histogram corresponded to the relative expression value for a particular gene.
The most significant GeneGo Pathway Maps were (1) immune response: HSP60 and HSP70/TLR signaling pathway; (2) immune response: Inflammasome in inflammatory response; (3) cell adhesion: plasmin signaling; (4) immune response: NF-AT signaling and leukocyte interactions. In addition, there are other pathways including Role of HMGB1 in dendritic cell maturation and migration; histamine signaling in dendritic cells: plasmin signaling in cell adhesion; cross-talk between VEGF and angiopoietin 1 signaling pathways; regulation of epithelial-to-mesenchymal transition (EMT); TGF-beta-dependent induction of EMT via SMADs in Development; role of IAP-proteins in apoptosis pathway in apoptosis and survival, and so forth. Meanwhile, immune system process, defense response, and response to stress were the most significantly enriched GO processes of these genes. With the disease folders, representing over 112 human diseases annotated by GeneGo, these 76 genes were mainly related to autoimmune diseases and some kinds of vascular inflammatory diseases.
The abstracted genes involved in significant pathways are summarized in Table 3.
3.4. STRING Analysis
Now specifically, we are interested in finding functional associations among these genes. We broadcast our data to STRING (a database of known and predicted protein interactions), which responds by displaying a network of nodes (proteins) connected by colored edges representing functional relationships.
Figure 2 summarizes the network of predicted associations between proteins encoded by these genes. The results indicate that CASP3, IL18, BLK, FCGR2B, FCGR2A, CRP, CCR5, CCL5, CCR3, CCL3L1, TNFRSF1A, TNF, IL4, ERAP1, LTA, CD40, NOD1, CTLA4, NLRP1, TGFBR2, SMAD3, TGFB2, VEGFA, KDR, and CCR2 are associated according to experimental evidence, with involvement in many signaling pathways; TNF was the key of nodes, linking to CRP, IL-4, CD40, CD40LG, IL-18, IL-10, and so on. They linked to many immune and inflammatory responses. All of these proteins (encoded by genes) are interrelated, forming a large network. However, many proteins are not linked to others, indicating that their functions are unrelated or unknown.
4.1. Immune Response in the Pathogenesis of KD
KD has long been considered as an abnormal immune disease. The activation of immune system and the cascade release of inflammatory factors are the important features in KD. A large number of T cells (increased activated CD4 T cells, depressed CD8 T cells and CD4+ CD25+ regulatory T cells), large mononuclear cells, macrophages and plasma cells, with a smaller number of neutrophils, are observed in various organ tissues of fatal cases of acute KD [102–106]. Additionally, various inflammatory cytokines and chemokines [107, 108], matrix metalloproteinases, nitric oxide production , autoantibody production [110, 111], and adhesive molecule expression [112, 113] are also overactivated in the acute stage of KD which are considered to facilitate vascular endothelial inflammation and then participate in the pathogenesis of KD and CAL formation. Go processes and DAVID analysis revealed that these genes are significantly enriched in immune responses which have the parallel results with clinical and laboratory findings. In addition, these genes are widely involved in other immune systemic and inflammatory diseases, for example, autoimmune thyroid disease, asthma, type I diabetes mellitus, allograft rejection, inflammatory bowel disease, vasculitis, arthritis, and rheumatic disease. Furthermore, the signal pathway produced in GeneGo contains many immune response pathways that participate in inflammation, apoptosis, injury, and remodeling process, which have been listed in Table 3.
4.2. ECM-Remodeling and Plasmin Signaling Pathway in the Pathogenesis of KD
In addition to the signal pathway of the immune response, ECM-remodeling and plasmin signaling pathway associated with cell adhesion were enriched in GeneGo MetaCore software ().
Numerous studies suggest that they participated in the pathophysiological process of KD. Activation of the fibrinolytic system, vascular injury, and remodeling were the prominent outcome in these pathways. Activated plasmin in the plasmin signaling pathway which is a major fibrinolytic protease can directly degrade fibrinogen, laminin, and fibronectin . On the cell surface, plasmin can activate a number of matrix metalloproteinases (MMPs) MMP1, MMP13 . Other MMPs (MMP-9 and so on) were subsequently activated. Moreover, IL-1 β, IL-6, TNF-α, and IFN-γ can stimulate the endothelial cells to produce more MMP-9. These MMPs degrade extracellular matrix proteins and components of basal membranes leading to the disruption of the internal elastic lamina and the trilaminar structure of the vascular wall [116–118]. Many examinations have showed that many MMPs were highly expressed in the acute stage of KD. MMPs are prominent during the remodeling process, contributing to the formation of coronary artery lesions , and consequently the intima proliferates and thickens, while in rare cases the vessel wall becomes stenotic or occluded by either stenosis or thrombosis. Endogenous tissue inhibitors of metalloproteinases (TIMPs) such as TIMP1, TIMP2, and TIMP3 can reduce excessive proteolytic ECM degradation by MMPs. The balance between MMPs and TIMPs controls the extent of ECM remodeling [120, 121]. One study indicated that MMPs and TIMPs were in a state of imbalance in KD patients . Therefore, ECM-remodeling and plasmin signaling pathway may have played a certain role in the vascular damage in KD.
4.3. NF-AT Signaling and Leukocyte Interactions
NF-AT signaling and leukocyte interactions ( value = 2.28 × 10−5) in the immune response cause our great concern. In this pathway, the activation of NFAT proteins is induced by the engagement of receptors that are coupled to the calcium/calcineurin signals, such as the antigen receptors that are expressed by T cells (TCR) and B cells (BCR), the Fc-epsilon receptors (e.g, Fc epsilon R1) that are expressed by mast cells and basophil cells or receptors coupled to heterotrimeric G-proteins (e.g., CCR3 on eosinophils) [123, 124] (Figure 3).
The NFAT signal is activated in T cell and can promote the expression of the immune-related genes. Antigen presenting cells present antigenic peptides to the T helper cell via major histocompatibility complex, class (II) (MHC class II). MHC class II can upregulate the expression of CD4+T cells and downregulate the expression of CD8+T cells which has been confirmed in acute phase of KD. Then, MHC class II peptides activate the T-cell receptor (TCR alpha/beta-CD3 complex) that starts a signal leading to the increase in cytosolic Ca(II) through both the transient release of calcium from intracellular stores and the influx of calcium through Ca(II) channels. That leads to activation of the calcium-regulated phosphatase, Calcineurin A. The activated Calcineurin A cleaves an inhibitory phosphate residue from the transcription fator NF-AT (e.g., NF-AT1 and NF-AT2). Consequently, NF-AT is transported into the nucleus, where it cooperates with other transcription factors for promoter binding and thereby induces the expression of cytokines and many other T-cell-activation-induced proteins. NF-AT in T cells is critical for the expression of a number of immunologically important genes, including IL-2, IL-4, IL-5, and IL-13, as well as several related membrane-bound proteins such as CD40 Ligand (CD40L) and Fas Ligand (Fasl) [125–127].
IL-4 plays an important role in cell-to-cell activation to activate NFAT signal to release leukotrienes and prostaglandins. Activated by NFAT signal in T cell, IL-4 activates nearby B cells that express corresponding receptor, IL-4R. In conjunction with BCR, IL-4 signaling pathway leads to the activation of several transcription factors, including nuclear factor kappa-B(NF-κB), signal transducer, and activator of transcription 6 (STAT6), that regulate immunoglobulin class switching and the production of immunoglobulin E (IgE) by some B cells [128–130]. IgE in turn activates NF-AT1 translocation and function in mast cells and basophils through the IgE receptor (Fc epsilon R1) leading to production of an array of cytokines, including IL-4, IL-5, and IL-13 [131, 132]. Fc epsilon R1 pathway also leads to activation of the cytosolic phospholipase A2 (Cpla2) that contributes to the secretion of leukotrienes and prostaglandins, the main mediators of inflammatory response . IL-4 and IL-13, in turn, activate epithelial cells and/or fibroblasts to release eosinophil-activating cytokines, such as chemokine ligand 11 (Eotaxin). These cytokines recruit eosinophils to the inflammatory focus in the tissue and induce intracellular signaling, mainly via chemokine receptor 3 (CCR3) activation, which leads to the leukotrienes and prostaglandins synthesis and also can use NF-AT1 transcription complex to activate cytokines and chemokines. IL-4 plays an important role in the interaction between the leukocytes and induces the release of variety of inflammatory mediators.
Additionally, CD40L activates nearby B cells that express corresponding receptor CD40. IL-2 binds to IL-2 receptors at the T Cells surface to drive clonal expansion of the activated cell that induces autocrine proliferation . Fasl activates the adjacent T Cells via binding to its receptors; FasR (CD95)  mediates apoptosis through the FAS signaling cascades (apoptosis). Fas-Fas ligand system has been considered to be involved in inducing apoptosis in KD resulting in marked decrease of peripheral blood lymphocytes .
4.3.1. What Is the NFATs?
NFATs are nuclear factors of activated T cells. The NFAT family consists of five members: NFAT1, NFAT2, NFAT3, NFAT4, NFAT4, and NFAT5. Four (except NFAT5) of these proteins are regulated by calcium signaling and four (except NFAT3) are expressed in the immune system . They are initially identified as Ca2+-sensitive transcription factors that regulate gene transcription in response to intracellular Ca2+signals. NFAT family members are expressed by almost every cell type, including the immune system and nonimmune cells, contributed to the regulation of immune response, as well as development and differentiation. In the immune system, NFATs have pivotal roles in the development and function of immune organs and regulate numerous physiological processes. With the best described effects on T cell activation and phenotype, NFATs also regulate gene expression in other immune cells such as B cells , mast cells [137, 138], eosinophils , basophils  and NK cells , macrophage , and dendritic cells . They can regulate the release of various cytokines in immune cells. In nonimmune cells, they regulate development and differentiation in a variety of organ systems . It has been examined that they control gene expression during remodeling and are activated by growth factors [144, 145] or histamine  in the endothelium, contributing to cell growth, remodeling of smooth muscle cells [147–149], and vascular development and angiogenesis [150–152] (including the isoforms c1 and c3) and are activated in response to inflammatory processes  and high intravascular pressure  in the vascular system. The isoforms NFATc3 and NFATc4 are active during pathophysiological conditions that affect the cardiovascular system, including atrial fibrillation [155, 156] and hypertrophy . Loss of specific NFAT isoforms has been found to result in cardiovascular, skeletal muscle, cartilage, neuronal, or immune system defects [158–162]. Therefore, we can conclude that the Ca2+/NFAT pathway plays a wide range role in inflammatory processes, immune responses, and the remodeling of vascular tissues. All of these physiological processes occur in KD. It is suggested that the Ca2+/NFAT pathway may involve in the pathological processes of KD.
4.3.2. The Upstream Adjustment Signals of NFAT Signal
NFATs are mainly Ca2+-sensitive transcription factors that regulate gene transcription in response to intracellular Ca2+ signals. Four (except NFAT5) of these proteins are regulated by calcium signaling. Activity of NFATs is regulated by phosphorylation. Inactive NFATs are highly phosphorylated and localized in the cytoplasm. Intracellular Ca2+ signals activate the calmodulin-dependent serine/threonine phosphatase calcineurin (CaN), which dephosphorylates NFATs and induces translocation to the nucleus.
Inositol-trisphosphate 3-kinase C (ITPKC) is a negative regulator of the Ca2+/NFAT pathway. NF-AT signaling was first mentioned to be associated with regulation of ITPKC in the KD. ITPKC is a kinase of inositol 1,4,5-triphosphate (IP3) which is a second messenger molecule that releases calcium from the endoplasmic and sarcoplasmic reticulum. First identified by genome-wide study and following confirmation by candidate genetic studies in both Japanese, Taiwanese and US children, ITPKC was considered to be associated with KD which confers both susceptibility to KD and the risk for CAL and IVIG resistance [26, 94, 95, 163], which has been thought to be involved in the Ca2+-dependent NFAT signaling pathways in T cells. It has been considered that C allele of rs28493229 in ITPKC can reduce the splicing efficiency of the ITPKC mRNA, inducing the hyperactivation of Ca2+-dependent NFAT signal in T cells, leading to a reduction in the phosphorylation of IP3 to IP4, resulting in the increase of IP3 levels. This would result in an increase of calcium levels and excessive activation of the NFAT signal, thus leading to immune dysregulation.
Caspase-3 (CASP3) is a key molecule of activation-induced cell death (AICD) ; it is profoundly related to the apoptosis of immune cells. It has also been reported to cleave the inositol 1,4,5-triphosphate receptor, type 1 (ITPR1) in apoptotic T cells (ITPR1 is a receptor for inositol 1,4,5-trisphosphate (IP3), a substrate for ITPKC in T cells ). Thereby, it is a positive regulatory factor of NFAT signal. Additionally, the mutation of CASP3 (rs113420705) can reduce the binding of NFAT to the DNA surrounding the SNP. Its gene variant (4q34-35, rs113420705) has been identified contributing to KD susceptibility in Euro-American triads and Taiwanese [35, 166]. Other studies [167, 168] also stated that CASP3 plays an important role in the execution phase of apoptosis of immune cells in KD.
Calcineurin inhibitors (e.g., CsA, FK506) have been extensively used as immunosuppressive agents to improve graft survival and to treat autoimmune diseases . They act by blocking calcineurin enzymatic activity. CsA has been an effective [169–171] therapeutic drug in the treatment of IVIG resistance patients in KD.
4.3.3. The Downstream of Adjustment Signals of NFAT Signal: NF-κB (Nuclear Factor Kappa-B)
NF-κB is another transcription factor of eukaryotes, which is evolutionarily related to the NF-AT family of transcription factors. It is activated in response to signals that lead to cell growth, differentiation, apoptosis, and other events. It takes part in expression of numerous cytokines and adhesion molecules which are critical elements involved in the regulation of immune responses.
NF-κB plays pivotal roles in the immune and inflammatory responses by regulating the interaction between CD40 and CD40L in T cells and B cells. NF-κB can be activated by IL-4 signaling pathway in B cells to induce the expression of CD40 which has been illustrated above. CD40 plays a crucial role as a costimulatory molecule in the cooperation between T and B cells. It is important in the pathogenesis of autoimmune diseases in humans and animal models such as autoimmune thyroiditis, inflammatory bowel disease, psoriasis, systemic lupus erythematosus, allergic encephalomyelitis, multiple sclerosis, rheumatoid arthritis, collagen-induced arthritis, and autoimmune type of diabetes mellitus [172–174]. CD40 signaling leads to isotype switching and autoantibody production in B cells and in T-cell priming, altering TCR expression through the expression and nuclear translocation of recombinases, which increases the risk of developing autoimmunity . CD40 engagement in both T or B cells leads to the production of cytokines, such as IL-12, IL-2, TNF-α, IFN-α, and CD80, developing an environment which is conducive to autoimmune diseases [172–174].
Additionally, the interaction between CD40 and CD40L regulated by NF-κB can regulate the expression of numerous biomolecules in other cells. They can enhance the expression of cytokines (such as IL-2, IL-6, IL-10, TNF-α, lymphotoxin-α, and transforming growth factor-β by B cells; the synthesis of granulocyte macrophage colony-stimulating factor (GM-CSF) by dendritic cells and eosinophils and the synthesis of TNF-α, IL-1, IL-6, and IL-8 by peripheral blood mononuclear cells), chemokines (monocyte chemotactic protein-1 (MCP-1), IL-8, MCP-1, matrix metalloproteinases (MMP-1,-2,-3,-9,-11, and -13) by peripheral blood mononuclear cells, macrophages, endothelial and smooth muscle cells endothelial), adhesion molecules (E-selectin, vascular cell adhesion molecule-1 (VCAM-1) and intercellular adhesion molecule-1 (ICAM-1) in endothelial cells and fibroblasts), platelet-activating factor , prostaglandin E2 , vascular endothelial growth factor [177, 178], and NO . which are involved in the pathophysiology of inflammatory and autoimmune diseases.
NF-κB may participate in the pathogenesis of vasculitis of KD in acute stage. Some studies have indicated that NF-κB is excessively activated in the acute phase of KD and the inhibition of NF-κB can reduce the generation of inflammatory cytokines which plays important roles in vascular damage of KD [179, 180]. NF-κB signaling pathway is a complex system; it perhaps involves in immune damage of KD in different levels. Activation of NF-κB can be used as the trigger of key links of the inflammatory response and induce the cascade release of inflammatory response factor, eventually leading to inflammatory pathological damage.
4.4. The NF-AT Signaling and Leukocyte Interactions and NF-κB Signaling Together May Be Involved in the Pathogenetic Process of KD
Given the important role of NFAT signaling and NF-κB signaling in the activation of immune system and the regulating of vascular remodeling, we speculate that the interaction between NFAT signaling and NF-κB signaling together may also be involved in the pathogenesis of KD.
Initially due to exposure to some inflammatory stimuli or certain pathogens, antigen presenting cells present antigenic peptides to the T-cell receptors via MHC class II leading to the stimulation of PLC-gamma 1 and hydrolyzation of PIP2. The second messengers IP3 in the T cells start a signal leading to the increase in cytosolic Ca(II) through both the transient release of calcium from intracellular stores and influx of calcium through Ca(II) channels. The high calcium levels lead to activation of the calcium-regulated phosphatase, Calcineurin A. The activated Calcineurin A cleaves an inhibitory phosphate residue from the transcription fator NF-AT (e.g., NF-AT1 and NF-AT2). Consequently, NF-AT is transported into the nucleus, where it cooperates with other transcription factors for promoter binding and activates T cells inducing the expression of a number of immunologically important genes including IL-2, IL-4, IL-5, IL-13, CD40 Ligand (CD40L), and Fas Ligand (Fasl). Through the leukocyte interactions, other immune cells were activated and release other inflammatory cytokines, such as leukotrienes and prostaglandins. In B cells and T-cell, CD40 signaling leads to isotype switching, autoantibody production, and altering TCR expression. CD40 signaling can also enhance the expression of cytokines, chemokines, matrix metalloproteinases, adhesion molecules, platelet-activating factors, prostaglandin E2, vascular endothelial growth factor, and NO. in other cells. The combined effect of these factors causes the vascular damage and formation of coronary artery lesions in KD. The process of NFAT signal in regulating development and differentiation was also excessively induced by the pathological damage of vasculature and then contributed to the remodeling of vascular system.
IL-4, CD40, and CD40L, which are enriched in the pathway of NF-AT signaling and leukocyte interactions and play a crucial role in the immune response and remodeling process, are located in the center position of the network (analysed by STRING) and are closely linked with the other factors. It further demonstrate the importance of this pathway.
KD is a complex disease. Many studies have shown that it is associated with a variety of gene polymorphism. Through GeneGo and DAVID analysis, we speculated that NF-AT signaling and leukocyte interactions combined with another transcription factor NF-κB may play an important role in pathological damage of KD. Their importance needs our follow-up clinical validation.
This work was financially supported by the Jiangsu Province Natural Science Foundation (no. BK2010032), Jiangsu province health department (H201127), and Suzhou Science and Technology Bureau (SYS201144). Additionally, the authors would like to thank the Systems Biology Center of Soochow University of China for their technical support.
The Supplementary Material was made by the software of DAVID (The Database for Annotation, Visualization and Integrated Discovery). It provides a comprehensive set of functional annotation tools for investigators to understand biological meaning behind large list of genes. The Supplementary Material contains the gene ontology analysis and the pathway analysis (KEGG pathway) of my genes which has been presented in Table 1 (candidate genetic studies), Table 2 (genome-wide studies). Gene ontology analysed the gene function from molecular function, biological process, cellular component. KEGG pathway showed the biological pathways and diseases these genes involved in.
T. Kawasaki, “Acute febrile mucocutaneous syndrome with lymphoid involvement with specific desquamation of the fingers and toes in children,” Japanese Journal of Allergology, vol. 16, no. 3, pp. 178–222, 1967.View at: Google Scholar
Y. C. Liu, C. P. Hou, C. M. Kuo, C. D. Liang, and H. C. Kuo, “Atypical Kawasaki disease: literature review and clinical nursing,” Journal of Nursing, vol. 57, no. 6, pp. 104–110, 2010.View at: Google Scholar
H. Yanagawa, Y. Nakamura, M. Yashiro et al., “Results of the nationwide epidemiologic survey of Kawasaki disease in 1995 and 1996 in Japan,” Pediatrics, vol. 102, no. 6, article E65, 1998.View at: Google Scholar
M. T. Lin, J. K. Wang, J. I. Yeh 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,” Pediatric Infectious Disease Journal, vol. 30, no. 2, pp. 148–152, 2011.View at: Publisher Site | Google Scholar
K. Yamamura, K. Ihara, K. Ikeda, H. Nagata, Y. Mizuno, and T. Hara, “Histo-blood group gene polymorphisms as potential genetic modifiers of the development of coronary artery lesions in patients with Kawasaki disease,” International Journal of Immunogenetics, vol. 39, no. 2, pp. 119–125, 2012.View at: Publisher Site | Google Scholar
Y. Y. Hsieh, C. C. Chang, C. M. Hsu, S. Y. Chen, W. H. Lin, and F. J. Tsai, “Major histocompatibility complex class I chain-related gene polymorphisms: associated with susceptibility to kawasaki disease and coronary artery aneurysms,” Genetic Testing and Molecular Biomarkers, vol. 15, no. 11, pp. 755–763, 2011.View at: Publisher Site | Google Scholar
F. Y. Huang, Y. J. Lee, M. R. Chen et al., “Polymorphism of transmembrane region of MICA gene and Kawasaki disease,” Experimental and Clinical Immunogenetics, vol. 17, no. 3, pp. 130–137, 2000.View at: Google Scholar
Y. Y. Hsieh, Y. J. Lin, C. C. Chang et al., “Human lymphocyte antigen B-associated transcript 2, 3, and 5 polymorphisms and haplotypes are associated with susceptibility of Kawasaki disease and coronary artery aneurysm,” Journal of Clinical Laboratory Analysis, vol. 24, no. 4, pp. 262–268, 2010.View at: Publisher Site | Google Scholar
J. Yang, C. R. Li, Y. B. Li et al., “The correlation between Kawasaki disease and polymorphisms of tumor necrosis factor α and interleukin-10 gene promoter,” Zhonghua Er Ke Za Zhi, vol. 41, no. 8, pp. 598–602, 2003.View at: Google Scholar
K. Chaudhuri, T. Singh Ahluwalia, S. Singh, G. Binepal, and M. Khullar, “Polymorphism in the promoter of the CCL5 gene (CCL5 G-403A) in a cohort of North Indian children with Kawasaki disease. A preliminary study,” Clinical and Experimental Rheumatology, vol. 29, no. 1, supplement 64, pp. S126–S130, 2011.View at: Google Scholar
J. A. Park, K. S. Shin, and W. K. Youn, “Polymorphism of matrix metalloproteinase-3 promoter gene as a risk factor for coronary artery lesions in Kawasaki disease,” Journal of Korean Medical Science, vol. 20, no. 4, pp. 607–611, 2005.View at: Google Scholar
J. Gao, H. Y. Wang, N. J. Wu, and S. H. Zhang, “Relationship between fibrinogen Bβ-148C/T polymorphism and coronary artery lesions in children with Kawasaki disease,” Zhongguo Dang Dai Er Ke Za Zhi, vol. 12, no. 7, pp. 518–520, 2010.View at: Google Scholar
Y. C. Huang, Y. J. Lin, J. S. Chang et al., “Single nucleotide polymorphism rs2229634 in the ITPR3 gene is associated with the risk of developing coronary artery aneurysm in children with Kawasaki disease,” International Journal of Immunogenetics, vol. 37, no. 6, pp. 439–443, 2010.View at: Publisher Site | Google Scholar
J. Yang, C. R. Li, Y. B. Li, H. J. Huang, R. X. Li, and G. B. Wang, “Correlation between mannose-binding lectin gene codon 54 polymorphism and susceptibility of Kawasaki disease,” Zhonghua Er Ke Za Zhi, vol. 42, no. 3, pp. 176–179, 2004.View at: Google Scholar
G. Simonini, E. Corinaldesi, C. Massai et al., “Macrophage migration inhibitory factor-173 polymorphism and risk of coronary alterations in children with Kawasaki disease,” Clinical and Experimental Rheumatology, vol. 27, no. 6, pp. 1026–1030, 2009.View at: Google Scholar
K. C. Hsueh, Y. J. Lin, J. S. Chang, L. Wan, and F. J. Tsai, “BTNL2 gene polymorphisms may be associated with susceptibility to Kawasaki disease and formation of coronary artery lesions in Taiwanese children,” European Journal of Pediatrics, vol. 169, no. 6, pp. 713–719, 2010.View at: Publisher Site | Google Scholar
S. W. Park, J. Y. Ban, K. L. Yoon et al., “Involvement of tryptophan hydroxylase 2 (TPH2) gene polymorphisms in susceptibility to coronary artery lesions in Korean children with Kawasaki disease,” European Journal of Pediatrics, vol. 169, no. 4, pp. 457–461, 2010.View at: Publisher Site | Google Scholar
K. C. Hsueh, Y. J. Lin, J. S. Chang et al., “Influence of interleukin 18 promoter polymorphisms in susceptibility to Kawasaki disease in Taiwan,” Journal of Rheumatology, vol. 35, no. 7, pp. 1408–1413, 2008.View at: Google Scholar
J. Y. Ban, K. L. Yoon, S. K. Kim, S. Kang, and J. H. Chung, “Promoter polymorphism (rs3755724, -55C/T) of tissue inhibitor of metalloproteinase 4 (TIMP4) as a risk factor for Kawasaki disease with coronary artery lesions in a Korean population,” Pediatric Cardiology, vol. 30, no. 3, pp. 331–335, 2009.View at: Publisher Site | Google Scholar
M. W. Quasney, D. E. Bronstein, R. M. Cantor et al., “Increased frequency of alleles associated with elevated tumor necrosis factor-α levels in children with Kawasaki disease,” Pediatric Research, vol. 49, no. 5, pp. 686–690, 2001.View at: Google Scholar
H. S. Yoon, S. K. Hae, S. Sohn, and M. H. Young, “Insertion/deletion polymorphism of angiotensin converting enzyme gene in Kawasaki disease,” Journal of Korean Medical Science, vol. 21, no. 2, pp. 208–211, 2006.View at: Google Scholar
R. Fukazawa, T. Sonobe, K. Hamamoto et al., “Possible synergic effect of angiotensin-I converting enzyme gene insertion/deletion polymorphism and angiotensin-II type-1 receptor 1166A/C gene polymorphism on ischemic heart disease in patients with Kawasaki disease,” Pediatric Research, vol. 56, no. 4, pp. 597–601, 2004.View at: Publisher Site | Google Scholar
T. Minami, H. Suzuki, T. Takeuchi, S. Uemura, J. Sugatani, and N. Yoshikawa, “A polymorphism in plasma platelet-activating factor acetylhydrolase is involved in resistance to immunoglobulin treatment in Kawasaki disease,” Journal of Pediatrics, vol. 147, no. 1, pp. 78–83, 2005.View at: Publisher Site | Google Scholar
I. H. Choi, Y. J. Chwae, W. S. Shim et al., “Clonal expansion of CD8+ T Cells in Kawasaki disease,” Journal of Immunology, vol. 159, no. 1, pp. 481–486, 1997.View at: Google Scholar
H. Fujiwara and Y. Hamashima, “Pathology of the heart in Kawasaki disease,” Pediatrics, vol. 61, no. 1, pp. 100–107, 1978.View at: Google Scholar
S. Amano, F. Hazama, and H. Kubagawa, “General pathology of Kawasaki disease. On the morphological alterations corresponding to the clinical manifestations,” Acta Pathologica Japonica, vol. 30, no. 5, pp. 681–694, 1980.View at: Google Scholar
Y. Mitani, H. Sawada, H. Hayakawa et al., “Elevated levels of high-sensitivity C-reactive protein and serum amyloid-A late after Kawasaki disease: association between inflammation and late coronary sequelae in Kawasaki disease,” Circulation, vol. 111, no. 1, pp. 38–43, 2005.View at: Publisher Site | Google Scholar
J. Verstappen and J. W. von den Hoff, “Tissue inhibitors of metalloproteinases (TIMPs): their biological functions and involvement in oral disease,” Journal of Dental Research, vol. 85, no. 12, pp. 1075–1084, 2006.View at: Google Scholar
H. Senzaki, S. Masutani, J. Kobayashi et al., “Circulating matrix metalloproteinases and their inhibitors in patients with Kawasaki disease,” Circulation, vol. 104, no. 8, pp. 860–863, 2001.View at: Google Scholar
F. Rusnak and P. Mertz, “Calcineurin: form and function,” Physiological Reviews, vol. 80, no. 4, pp. 1483–1521, 2000.View at: Google Scholar
M. Lee and J. Park, “Regulation of NFAT activation: a potential therapeutic target for immunosuppression,” Molecules and Cells, vol. 22, no. 1, pp. 1–7, 2006.View at: Google Scholar
K. Silver and R. J. Cornall, “Isotype control of B cell signaling,” Science's STKE, vol. 2003, no. 184, article pe21, 2003.View at: Google Scholar
T. Mizuno and T. L. Rothstein, “B cell receptor (BCR) cross-talk: CD40 engagement enhances BCR-induced ERK activation,” Journal of Immunology, vol. 174, no. 6, pp. 3369–3376, 2005.View at: Google Scholar
A. Lorentz, I. Klopp, T. Gebhardt, M. P. Manns, and S. C. Bischoff, “Role of activator protein 1, nuclear factor-κB, and nuclear factor of activated T cells in IgE receptor-mediated cytokine expression in mature human mast cells,” Journal of Allergy and Clinical Immunology, vol. 111, no. 5, pp. 1062–1068, 2003.View at: Publisher Site | Google Scholar
H. Y. Kim, H. G. Lee, and D. S. Kim, “Apoptosis of peripheral blood mononuclear cells in Kawasaki disease,” Journal of Rheumatology, vol. 27, no. 3, pp. 801–806, 2000.View at: Google Scholar
J. T. Schroeder, K. Miura, H. H. Kim, A. Sin, A. Cianferoni, and V. Casolaro, “Selective expression of nuclear factor of activated T cells 2/c1 in human basophils: evidence for involvement in IgE-mediated IL-4 generation,” Journal of Allergy and Clinical Immunology, vol. 109, no. 3, pp. 507–513, 2002.View at: Publisher Site | Google Scholar
J. Aramburu, L. Azzoni, A. Rao, and B. Perussia, “Activation and expression of the nuclear factors of activated T cells, NFATp and NFATc, in human natural killer cells: regulation upon CD16 ligand binding,” Journal of Experimental Medicine, vol. 182, no. 3, pp. 801–810, 1995.View at: Publisher Site | Google Scholar
A. Yarilina, K. Xu, J. Chen, and L. B. Ivashkiv, “TNF activates calcium-nuclear factor of activated T cells (NFAT)c1 signaling pathways in human macrophages,” Proceedings of the National Academy of Sciences of the United States of America, vol. 108, no. 4, pp. 1573–1578, 2011.View at: Publisher Site | Google Scholar
L. Hadri, C. Pavoine, L. Lipskaia, S. Yacoubi, and A. M. Lompré, “Transcription of the sarcoplasmic/endoplasmic reticulum Ca2+-ATPase type 3 gene, ATP2A3, is regulated by the calcineurin/NFAT pathway in endothelial cells,” Biochemical Journal, vol. 394, no. 1, pp. 27–33, 2006.View at: Publisher Site | Google Scholar
V. Boss, X. Wang, L. F. Koppelman, K. Xu, and T. J. Murphy, “Histamine induces nuclear factor of activated T cell-mediated transcription and cyclosporin A-sensitive interleukin-8 mRNA expression in human umbilical vein endothelial cells,” Molecular Pharmacology, vol. 54, no. 2, pp. 264–272, 1998.View at: Google Scholar
V. Boss, K. L. Abbott, X. F. Wang, G. K. Pavlath, and T. J. Murphy, “The cyclosporin A-sensitive nuclear factor of activated T cells (NFAT) proteins are expressed in vascular smooth muscle cells. Differential localization of NFAT isoforms and induction of NFAT-mediated transcription by phospholipase C-coupled cell surface receptors,” The Journal of Biological Chemistry, vol. 273, no. 31, pp. 19664–19671, 1998.View at: Publisher Site | Google Scholar
L. V. Gonzalez Bosc, M. K. Wilkerson, K. N. Bradley, D. M. Eckman, D. C. Hill-Eubanks, and M. T. Nelson, “Intraluminal pressure is a stimulus for NFATc3 nuclear accumulation: role of calcium, endothelium-derived nitric oxide, and cGMP-dependent protein kinase,” The Journal of Biological Chemistry, vol. 279, no. 11, pp. 10702–10709, 2004.View at: Publisher Site | Google Scholar
J. Hirota, T. Furuichi, and K. Mikoshiba, “Inositol 1,4,5-trisphosphate receptor type I is a substrate for caspase-3 and is cleaved during apoptosis in a caspase-3-dependent manner,” The Journal of Biological Chemistry, vol. 274, no. 48, pp. 34433–34437, 1999.View at: Publisher Site | Google Scholar
Q. J. Yi, C. R. Li, and X. Q. Yang, “Effect of intravenous immunoglobulin on inhibiting peripheral blood lymphocyte apoptosis in acute Kawasaki disease,” Acta Paediatrica, vol. 90, no. 6, pp. 623–627, 2001.View at: Google Scholar
U. Schönbeck and P. Libby, “The CD40/CD154 receptor/ligand dyad,” Cellular and Molecular Life Sciences, vol. 58, no. 1, pp. 4–43, 2001.View at: Google Scholar
M. E. Munroe and G. A. Bishop, “A costimulatory function for T cell CD40,” Journal of Immunology, vol. 178, no. 2, pp. 671–682, 2007.View at: Google Scholar
J. D. Laman, M. de Boer, and B. A. Hart, “CD40 in clinical inflammation: from multiple sclerosis to atherosclerosis,” Developmental Immunology, vol. 6, no. 3-4, pp. 215–222, 1998.View at: Google Scholar
S. Russo, B. Bussolati, I. Deambrosis, F. Mariano, and G. Camussi, “Platelet-activating factor mediates CD40-dependent angiogenesis and endothelial-smooth muscle cell interaction,” Journal of Immunology, vol. 171, no. 10, pp. 5489–5497, 2003.View at: Google Scholar
Y. Inoue, T. Otsuka, H. Niiro et al., “Novel regulatory mechanisms of CD40-induced prostanoid synthesis by IL-4 and IL-10 in human monocytes,” Journal of Immunology, vol. 172, no. 4, pp. 2147–2154, 2004.View at: Google Scholar
M. Melter, M. E. J. Reinders, M. Sho et al., “Ligation of CD40 induces the expression of vascular endothelial growth factor by endothelial cells and monocytes and promotes angiogenesis in vivo,” Blood, vol. 96, no. 12, pp. 3801–3808, 2000.View at: Google Scholar