Search

Background: Early identification of children with Kawasaki Disease (KD) is key for timely initiation of intravenous immunoglobulin (IVIG) therapy. However, the diagnosis of the disease remains challenging, especially in children with an incomplete presentation (inKD). Moreover, we currently lack objective tools for identification of non-response (NR) to IVIG., Methods: Children with KD were enrolled and samples obtained before IVIG treatment and sequentially at 24 h and 4-6 weeks post-IVIG in a subset of patients. We also enrolled children with other febrile illnesses [adenovirus (AdV); group A streptococcus (GAS)] and healthy controls (HC) for comparative analyses. Blood transcriptional profiles were analyzed to define: a) the cKD and inKD biosignature, b) compare the KD signature with other febrile illnesses and, c) identify biomarkers predictive of clinical outcomes., Results: We identified a cKD biosignature (n = 39; HC, n = 16) that was validated in two additional cohorts of children with cKD (n = 37; HC, n = 20) and inKD (n = 13; HC, n = 8) and was characterized by overexpression of inflammation, platelets, apoptosis and neutrophil genes, and underexpression of T and NK cell genes. Classifier genes discriminated KD from adenovirus with higher sensitivity and specificity (92% and 100%, respectively) than for GAS (75% and 87%, respectively). We identified a genomic score (MDTH) that was higher at baseline in IVIG-NR [median 12,290 vs. 5,572 in responders, p = 0.009] and independently predicted IVIG-NR., Conclusion: A reproducible biosignature from KD patients was identified, and was similar in children with cKD and inKD. A genomic score allowed early identification of children at higher risk for non-response to IVIG., Competing Interests: The authors have declared that no competing interests exist.

Background: The pathogenesis of Kawasaki disease (KD) is commonly ascribed to an exaggerated immunologic response to an unidentified environmental or infectious trigger in susceptible children. A comprehensive framework linking epidemiological data and global distribution of KD has not yet been proposed., Methods and Findings: Patients with KD (n = 81) were enrolled within 6 weeks of diagnosis along with control subjects (n = 87). All completed an extensive epidemiological questionnaire. Geographic localization software characterized the subjects' neighborhood. KD incidence was compared to atmospheric biological particles counts and winds patterns. These data were used to create a comprehensive risk framework for KD, which we tested against published data on the global distribution. Compared to controls, patients with KD were more likely to be of Asian ancestry and were more likely to live in an environment with low exposure to environmental allergens. Higher atmospheric counts of biological particles other than fungus/spores were associated with a temporal reduction in incidence of KD. Finally, westerly winds were associated with increased fungal particles in the atmosphere and increased incidence of KD over the Greater Toronto Area. Our proposed framework was able to explain approximately 80% of the variation in the global distribution of KD. The main limitations of the study are that the majority of data used in this study are limited to the Canadian context and our proposed disease framework is theoretical and circumstantial rather than the result of a single simulation., Conclusions: Our proposed etiologic framework incorporates the 1) proportion of population that are genetically susceptible; 2) modulation of risk, determined by habitual exposure to environmental allergens, seasonal variations of atmospheric biological particles and contact with infectious diseases; and 3) exposure to the putative trigger. Future modelling of individual risk and global distribution will be strengthened by taking into consideration all of these non-traditional elements.

Kawasaki disease (KD) is an acute systemic vasculitis that can potentially cause coronary artery aneurysms in some children. KD occurs approximately 1.5 times more frequently in males than in females. To identify sex-specific genetic variants that are involved in KD pathogenesis in children, we performed a sex-stratified genome-wide association study (GWAS), using the Illumina HumanOmni1-Quad BeadChip data (249 cases and 1,000 controls) and a replication study for the 34 sex-specific candidate SNPs in an independent sample set (671 cases and 3,553 controls). Male-specific associations were detected in three common variants: rs1801274 in FCGR2A [odds ratio (OR) = 1.40, P = 9.31 × 10-5], rs12516652 in SEMA6A (OR = 1.87, P = 3.12 × 10-4), and rs5771303 near IL17REL (OR = 1.57, P = 2.53 × 10-5). The male-specific association of FCGR2A, but not SEMA6A and IL17REL, was also replicated in a Japanese population (OR = 1.74, P = 1.04 × 10-4 in males vs. OR = 1.22, P = 0.191 in females). In a meta-analysis with 1,461 cases and 5,302 controls, a very strong association of KD with the nonsynonymous SNP rs1801274 (p.His167Arg, previously assigned as p.His131Arg) in FCGR2A was confirmed in males (OR = 1.48, P = 1.43 × 10-7), but not in the females (OR = 1.17, P = 0.055). The present study demonstrates that p.His167Arg, a KD-associated FCGR2A variant, acts as a susceptibility gene in males only. Overall, the gender differences associated with FCGR2A in KD provide a new insight into KD susceptibility.

Kawasaki disease (KD) is the most common acquired pediatric heart disease. We analyzed Whole Genome Sequences (WGS) from a 6-member African American family in which KD affected two of four children. We sought rare, potentially causative genotypes by sequentially applying the following WGS filters: sequence quality scores, inheritance model (recessive homozygous and compound heterozygous), predicted deleteriousness, allele frequency, genes in KD-associated pathways or with significant associations in published KD genome-wide association studies (GWAS), and with differential expression in KD blood transcriptomes. Biologically plausible genotypes were identified in twelve variants in six genes in the two affected children. The affected siblings were compound heterozygous for the rare variants p.Leu194Pro and p.Arg247Lys in Toll-like receptor 6 (TLR6), which affect TLR6 signaling. The affected children were also homozygous for three common, linked (r2 = 1) intronic single nucleotide variants (SNVs) in TLR6 (rs56245262, rs56083757 and rs7669329), that have previously shown association with KD in cohorts of European descent. Using transcriptome data from pre-treatment whole blood of KD subjects (n = 146), expression quantitative trait loci (eQTL) analyses were performed. Subjects homozygous for the intronic risk allele (A allele of TLR6 rs56245262) had differential expression of Interleukin-6 (IL-6) as a function of genotype (p = 0.0007) and a higher erythrocyte sedimentation rate at diagnosis. TLR6 plays an important role in pathogen-associated molecular pattern recognition, and sequence variations may affect binding affinities that in turn influence KD susceptibility. This integrative genomic approach illustrates how the analysis of WGS in multiplex families with a complex genetic disease allows examination of both the common disease-common variant and common disease-rare variant hypotheses., Competing Interests: Whole genome sequencing was donated by Jay Flatley who was the CEO of Illumina at that time. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Kawasaki disease (KD) or Kawasaki syndrome is known as a vasculitis of small to medium-sized vessels, and coronary arteries are predominantly involved in childhood. Generally, 20-25% of untreated with IVIG and 3-5% of treated KD patients have been developed coronary artery lesions (CALs), such as dilatation and aneurysm. Understanding how coronary artery aneurysms (CAAs) are established and maintained in KD patients is therefore of great importance. Upon our previous genotyping data of 157 valid KD subjects, a genome-wide association study (GWAS) has been conducted among 11 (7%) CAA-developed KD patients to reveal five significant genetic variants passed pre-defined thresholds and resulted in two novel susceptibility protein-coding genes, which are NEBL (rs16921209 (P = 7.44 × 10(-9); OR = 32.22) and rs7922552 (P = 8.43 × 10(-9); OR = 32.0)) and TUBA3C (rs17076896 (P = 8.04 × 10(-9); OR = 21.03)). Their known functions have been reported to associate with cardiac muscle and tubulin, respectively. As a result, this might imply their putative roles of establishing CAAs during KD progression. Additionally, various model analyses have been utilized to determine dominant and recessive inheritance patterns of identified susceptibility mutations. Finally, all susceptibility genes hit by significant genetic variants were further investigated and the top three representative gene-ontology (GO) clusters were regulation of cell projection organization, neuron recognition, and peptidyl-threonine phosphorylation. Our results help to depict the potential routes of the pathogenesis of CAAs in KD patients and will facilitate researchers to improve the diagnosis and prognosis of KD in personalized medicine.

Kawasaki disease (KD; MIM#61175) is a systemic vasculitis syndrome with unknown etiology which predominantly affects infants and children. Recent findings of susceptibility genes for KD suggest possible involvement of the Ca(2+)/NFAT pathway in the pathogenesis of KD. ORAI1 is a Ca(2+) release activated Ca(2+) (CRAC) channel mediating store-operated Ca(2+) entry (SOCE) on the plasma membrane. The gene for ORAI1 is located in chromosome 12q24 where a positive linkage signal was observed in our previous affected sib-pair study of KD. A common non-synonymous single nucleotide polymorphism located within exon 2 of ORAI1 (rs3741596) was significantly associated with KD (P = 0.028 in the discovery sample set (729 KD cases and 1,315 controls), P = 0.0056 in the replication sample set (1,813 KD cases vs. 1,097 controls) and P = 0.00041 in a meta-analysis by the Mantel-Haenszel method). Interestingly, frequency of the risk allele of rs3741596 is more than 20 times higher in Japanese compared to Europeans. We also found a rare 6 base-pair in-frame insertion variant associated with KD (rs141919534; 2,544 KD cases vs. 2,414 controls, P = 0.012). These data indicate that ORAI1 gene variations are associated with KD and may suggest the potential importance of the Ca(2+)/NFAT pathway in the pathogenesis of this disorder.

Kawasaki disease (KD) is a systemic vasculitis primarily affecting children < 5 years old. Genes significantly associated with KD mostly involve cardiovascular, immune, and inflammatory responses. Recent studies have observed stronger associations for KD risk with multiple genes compared to individual genes. Therefore, we investigated whether gene combinations influenced KD susceptibility or coronary artery lesion (CAL) formation. We examined 384 single-nucleotide polymorphisms (SNPs) for 159 immune-related candidate genes in DNA samples from KD patients with CAL (n = 73), KD patients without CAL (n = 153), and cohort controls (n = 575). Individual SNPs were first assessed by univariate analysis (UVA) and multivariate analysis (MVA). We used multifactor dimensionality reduction (MDR) to examine individual SNPs in one-, two-, and three-locus best fit models. UVA identified 53 individual SNPs that were significantly associated with KD risk or CAL formation (p < 0.10), while 35 individual SNPs were significantly associated using MVA (p ≤ 0.05). Significant associations in MDR analysis were only observed for the two-locus models after permutation testing (p ≤ 0.05). In logistic regression, combined possession of PDE2A (rs341058) and CYFIP2 (rs767007) significantly increased KD susceptibility (OR = 3.54; p = 4.14 x 10(-7)), while combinations of LOC100133214 (rs2517892) and IL2RA (rs3118470) significantly increased the risk of CAL in KD patients (OR = 5.35; p = 7.46 x 10(-5)). Our results suggest varying gene-gene associations respectively predispose individuals to KD risk or its complications of CAL.

Background: Kawasaki disease (KD) is a systemic vasculitis with unknown etiology mainly affecting children in Asian countries. Dendritic cell-specific intercellular adhesion molecule-3 grabbing non-integrin (DC-SIGN, CD209) in humans was showed to trigger an anti-inflammatory cascade and associated with KD susceptibility. This study was conducted to investigate the association between genetic polymorphisms of CD209 and the risk KD., Methods: A total of 948 subjects (381 KD and 567 controls) were recruited. Nine tagging SNPs (rs8112310, rs4804800, rs11465421, rs1544766, rs4804801, rs2287886, rs735239, rs735240, rs4804804) were selected for TaqMan allelic discrimination assay. Clinical phenotypes, coronary artery lesions (CAL) and intravenous immunoglobulin (IVIG) treatment outcomes were collected for analysis., Results: Significant associations were found between CD209 polymorphisms (rs4804800, rs2287886, rs735240) and the risk of KD. Haplotype analysis for CD209 polymorphisms showed that A/A/G haplotype (P = 0.0002, OR = 1.61) and G/A/G haplotype (P = 0.0365, OR = 1.52) had higher risk of KD as compared with G/G/A haplotype in rs2287886/rs735239/rs735240 pairwise allele analysis. There were no significant association in KD with regards to CAL formation and IVIG treatment responses., Conclusion: CD209 polymorphisms were responsible for the susceptibility of KD, but not CAL formation and IVIG treatment responsiveness.

Objectives: Recent genome-wide association study found rs1801274, a functional single nucleotide polymorphism (SNP) in IgG receptor gene FCGR2A, was associated with increased risk of Kawasaki disease (KD). However, subsequent studies on the role of this SNP were limited and controversial., Methods: A case-control study was conducted in a Chinese Han population including 428 KD patients and 493 controls to examine the association between rs1801274 and KD susceptibility. A meta-analysis was performed in combination with the relevant published studies to further clarify such an association., Results: Our case-control study found that rs1801274 was significantly associated with increased risk of KD in the Chinese Han population, with an odds ratio (OR) of 1.58 (95% CI = 0.96-2.62) for the GA genotype and 1.93 (95% CI = 1.16-3.19) for the AA genotype compared with the GG genotype. The result of meta-analysis further demonstrated that the A allele of rs1801274 was significantly correlated with KD risk under the allelic model (OR = 1.35, 95% CI = 1.27-1.44) without heterogeneity by fixed-effects model analysis (Q = 17.30, p = 0.139). Moreover, sensitivity analysis supported the robustness of this meta-analysis., Conclusion: These results further confirm that rs1801274 in the FCGR2A gene is significantly associated with increased risk of KD.

Kawasaki disease (KD) is a multi-systemic vasculitis that preferentially affects children. A single nucleotide polymorphism (SNP) in inositol 1,4,5-trisphosphate 3-kinase C (ITPKC) has been identified to be an important polymorphism in the risk of KD. This study was conducted to comprehensively investigate the associations between all tagging SNPs of ITPKC in the risk of KD in a Taiwanese population. A total of 950 subjects (381 KD patients and 569 controls) were recruited. Seven tagging SNPs (rs11673492, rs7257602, rs7251246, rs890934, rs10420685, rs2607420, rs2290692) were selected for TaqMan allelic discrimination assay. Clinical data of coronary artery lesions (CAL) and aneurysms were collected for analysis. A significant association was found between rs7251246 in ITPKC and CAL formation. Haplotype analysis for ITPKC polymorphisms also confirmed this association in the patients with CAL and aneurysm formation. This is the first study to identify that SNP rs7251246 in ITPKC is associated with the severity of KD.

Kawasaki disease (KD) is pediatric systemic vasculitis with the classic complication of coronary artery aneurysm (CAA). It is the leading cause of acquired cardiovascular diseases in children. Some severe cases present with multi-organ involvement or neurological dysfunction. To identify the role of the glutamate receptor, ionotropic, N-methyl-d-aspartate 3A (GRIN3A) in KD, we investigated genetic variations in GRIN3A in a Taiwanese cohort of 262 KD patients (76 with and 186 without CAA complications). We used univariate and multivariate regression analyses to identify the associations between clinical characteristics and GRIN3A genetic variations in KD. According to univariate regression analysis, CAA formation in KD was significantly associated with fever duration (p < 0.0001), first Intravenous immunoglobulin (IVIG) used (days after day one of fever) (p < 0.0001), and the GRIN3A (rs7849782) genetic variant (p < 0.001). KD patients with GG+GC genotype showed a lower rate of developing CAA (GG+GC genotype: odds ratio = 0.26; 95% CI = 0.14-0.46). Significant associations were identified between KD with CAA complication and the GRIN3A (rs7849782) genetic variant by using multivariate regression analysis. Specifically, significant correlations were observed between KD with CAA complications and the presence of GG+GC genotypes for the GRIN3A rs7849782 single-nucleotide polymorphism (full model: odds ratio = 0.25; 95% CI = 0.14-0.46). Our results suggest that a polymorphism of the GRIN3A gene may play a role in KD pathogenesis.

The BLK and CD40 loci have been associated with Kawasaki disease (KD) in two genome-wide association studies (GWAS) conducted in a Taiwanese population of Han Chinese ancestry (Taiwanese) and in Japanese cohorts. Here we build on these findings with replication studies of the BLK and CD40 loci in populations of Korean and European descent. The BLK region was significantly associated with KD susceptibility in both populations. Within the BLK gene the rs2736340-located linkage disequilibrium (LD ) comprising the promoter and first intron was strongly associated with KD, with the combined results of Asian studies including Taiwanese, Japanese, and Korean populations (2,539 KD patients and 7,021 controls) providing very compelling evidence of association (rs2736340, OR = 1.498, 1.354-1.657; P = 4.74×10(-31)). We determined the percentage of B cells present in the peripheral blood mononuclear cell (PBMC) population and the expression of BLK in the peripheral blood leukocytes (leukocytes) of KD patients during the acute and convalescent stages. The percentage of B cells in the PBMC population and the expression of BLK in leukocytes were induced in patients in the acute stage of KD. In B cell lines derived from KD patients, and in purified B cells from KD patients obtained during the acute stage, those with the risk allele of rs2736340 expressed significantly lower levels of BLK. These results suggest that peripheral B cells play a pathogenic role during the acute stage of KD. Decreased BLK expression in peripheral blood B cells may alter B cell function and predispose individuals to KD. These associative data suggest a role for B cells during acute KD. Understanding the functional implications may facilitate the development of B cell-mediated therapy for KD.

Single-nucleotide polymorphisms (SNPs) in inositol 1,4,5-trisphosphate 3-kinase C (ITPKC, rs28493229) and caspase-3 (CASP3, rs113420705) are associated with susceptibility to KD in Japanese and Taiwanese populations. This study was conducted to investigate the involvement of these 2 SNPs in the risk for intravenous immunoglobulin (IVIG) resistance and coronary artery lesion (CAL) in Taiwanese population. A total of 340 KD patients were subjected to assess by the identification of 2-locus genes model. A combinatorial association between ITPKC (rs28493229) and CASP3 (rs113420705) was found in CAL formation (P = 0.0227, OR: 3.06). KD patients with high-risk genotype had a trend of overrepresentation in IVIG resistance compared with individual SNPs. Our findings suggest the existence of genetic factors affecting patients' risk for CAL formation and IVIG responsiveness in a Taiwanese population.

N-glycolylneuraminic acid (Neu5Gc) is an immunogenic sugar of dietary origin that metabolically incorporates into diverse native glycoconjugates in humans. Anti-Neu5Gc antibodies are detected in all human sera, though with variable levels and epitope-recognition profiles. These antibodies likely play a role in several inflammation-mediated pathologies including cardiovascular diseases and cancer. In cancer, they have dualistic and opposing roles, either stimulating or repressing disease, as a function of their dose, and some of these antibodies serve as carcinoma biomarkers. Thus, anti-Neu5Gc antibodies may signify risk of inflammation-mediated diseases, and changes in their levels could potentially be used to monitor disease progression and/or response to therapy. Currently, it is difficult to determine levels of anti-Neu5Gc antibodies in individual human samples because these antibodies recognize multiple Neu5Gc-epitopes. Here we describe a simple and specific method for detection and overall estimation of human anti-Neu5Gc antibodies. We exploit the difference between two mouse models that differ only by Neu5Gc-presence (wild-type) or Neu5Gc-absence (Cmah(-/-) knockout). We characterize mouse serum from both strains by HPLC, lectin and mass-spectrometry analysis and show the target Neu5Gc-epitopes. We then use Cmah(-/-) knockout sera to inhibit all non-Neu5Gc-reactivity followed by binding to wild-type sera to detect overall anti-Neu5Gc response in a single assay. We applied this methodology to characterize and quantify anti-Neu5Gc IgG and IgA in sera of patients with Kawasaki disease (KD) at various stages compared to controls. KD is an acute childhood febrile disease characterized by inflammation of coronary arteries that untreated may lead to coronary artery aneurysms with risk of thrombosis and myocardial infarction. This estimated response is comparable to the average of detailed anti-Neu5Gc IgG profile analyzed by a sialoglycan microarray. Both assays revealed an elevated response in acute KD patients with normal coronaries compared to patients with aneurysm or dilated coronaries. Implications of these findings are discussed.

Background: Kawasaki disease is an acute, self-limited vasculitis of childhood that can result in structural damage to the coronary arteries. Previous studies have implicated the TGF-β pathway in disease pathogenesis and generation of myofibroblasts in the arterial wall. microRNAs are small non-coding RNAs that modulate gene expression at the post-transcriptional level and can be transported between cells in extracellular vesicles. To understand the role that microRNAs play in modifying gene expression in Kawasaki disease, we studied microRNAs from whole blood during the acute and convalescent stages of the illness., Methodology/principal Findings: RNA isolated from the matched whole blood of 12 patients with acute and convalescent Kawasaki disease were analyzed by sequencing of small RNA. This analysis revealed six microRNAs (miRs-143, -199b-5p, -618, -223, -145 and -145* (complementary strand)) whose levels were significantly elevated during the acute phase of Kawasaki disease. The result was validated using targeted qRT-PCR using an independent cohort (n = 16). miR-145, which plays a critical role in the differentiation of neutrophils and vascular smooth muscle cells, was expressed at high levels in blood samples from acute Kawasaki disease but not adenovirus-infected control patients (p = 0.005). miR-145 was also detected in small extracellular vesicles isolated from acute Kawasaki disease plasma samples. Pathway analysis of the predicted targets of the 6 differentially expressed microRNAs identified the TGF-β pathway as the top pathway regulated by microRNAs in Kawasaki disease., Conclusion: Sequencing of small RNA species allowed discovery of microRNAs that may participate in Kawasaki disease pathogenesis. miR-145 may participate, along with other differentially expressed microRNAs, in regulating expression of genes in the TGF-β pathway during the acute illness. If the predicted target genes are confirmed, our findings suggest a model of Kawasaki disease pathogenesis whereby miR-145 modulates TGF-β signaling in the arterial wall.

Background: Kawasaki disease (KD) of unknown immunopathogenesis is an acute febrile systemic vasculitis and the leading cause of acquired heart diseases in childhood. To search for a better strategy for the prevention and treatment of KD, this study compared and validated human KD immunopathogenesis in a mouse model of Lactobacillus casei cell wall extract (LCWE)-induced coronary arteritis., Methods: Recruited subjects fulfilled the criteria of KD and were admitted for intravenous gamma globulin (IVIG) treatment at the Kaohsiung Chang Gung Memorial Hospital from 2001 to 2009. Blood samples from KD patients were collected before and after IVIG treatment, and cardiovascular abnormalities were examined by transthoracic echocardiography. Wild-type male BALB/c mice (4-week-old) were intraperitoneally injected with LCWE (1 mg/mL) to induce coronary arteritis. The induced immune response in mice was examined on days 1, 3, 7, and 14 post injections, and histopathology studies were performed on days 7 and 14., Results: Both human KD patients and LCWE-treated mice developed coronary arteritis, myocarditis, valvulitis, and pericarditis, as well as elevated plasma levels of interleukin (IL)-2, IL-6, IL-10, monocyte chemoattractant protein (MCP)-1, and tumor necrosis factor (TNF)-α in acute phase. Most of these proinflammatory cytokines declined to normal levels in mice, whereas normal levels were achieved in patients only after IVIG treatment, with a few exceptions. Toll-like receptor (TLR)-2, but not TLR4 surface enhancement on circulating CD14+ monocytes, was augmented in KD patients before IVIG treatment and in LCWE-treated mice, which declined in patients after IVIG treatment., Conclusion: This result suggests that that not only TLR2 augmentation on CD14+ monocytes might be an inflammatory marker for both human KD patients and LCWE-induced CAL mouse model but also this model is feasible for studying therapeutic strategies of coronary arteritis in human KD by modulating TLR2-mediated immune activation on CD14+ monocytes.

Kawasaki disease (KD) is characterized by systemic vasculitis with unknown etiology. Previous studies from Japan indicated that a gene polymorphism of ITPKC (rs28493229) is responsible for susceptibility to KD. We collected DNA samples from 1,531 Taiwanese subjects (341 KD patients and 1,190 controls) for genotyping ITPKC. In this study, no significant association was noted for the ITPKC polymorphism (rs28493229) between the controls and KD patients, although the CC genotype was overrepresented. We further combined our data with previously published case/control KD studies in the Taiwanese population and performed a meta-analysis. A significant association between rs28493229 and KD was found (Odds Ratio:1.36, 95% Confidence Interval 1.12-1.66). Importantly, a significant association was obtained between rs28493229 and KD patients with aneurysm formation (P = 0.001, under the recessive model). Taken together, our results indicated that C-allele of ITPKC SNP rs28493229 is associated with the susceptibility and aneurysm formation in KD patients in a Taiwanese population.

Kawasaki disease (KD) is an acute systemic vasculitis syndrome that primarily affects infants and young children. Its etiology is unknown; however, epidemiological findings suggest that genetic predisposition underlies disease susceptibility. Taiwan has the third-highest incidence of KD in the world, after Japan and Korea. To investigate novel mechanisms that might predispose individuals to KD, we conducted a genome-wide association study (GWAS) in 250 KD patients and 446 controls in a Han Chinese population residing in Taiwan, and further validated our findings in an independent Han Chinese cohort of 208 cases and 366 controls. The most strongly associated single-nucleotide polymorphisms (SNPs) detected in the joint analysis corresponded to three novel loci. Among these KD-associated SNPs three were close to the COPB2 (coatomer protein complex beta-2 subunit) gene: rs1873668 (p = 9.52×10⁻⁵), rs4243399 (p = 9.93×10⁻⁵), and rs16849083 (p = 9.93×10⁻⁵). We also identified a SNP in the intronic region of the ERAP1 (endoplasmic reticulum amino peptidase 1) gene (rs149481, p(best) = 4.61×10⁻⁵). Six SNPs (rs17113284, rs8005468, rs10129255, rs2007467, rs10150241, and rs12590667) clustered in an area containing immunoglobulin heavy chain variable regions genes, with p(best)-values between 2.08×10⁻⁵ and 8.93×10⁻⁶, were also identified. This is the first KD GWAS performed in a Han Chinese population. The novel KD candidates we identified have been implicated in T cell receptor signaling, regulation of proinflammatory cytokines, as well as antibody-mediated immune responses. These findings may lead to a better understanding of the underlying molecular pathogenesis of KD.

Background: The etiology of Kawasaki Disease (KD) is enigmatic, although an infectious cause is suspected. Polymorphisms in CC chemokine receptor 5 (CCR5) and/or its potent ligand CCL3L1 influence KD susceptibility in US, European and Korean populations. However, the influence of these variations on KD susceptibility, coronary artery lesions (CAL) and response to intravenous immunoglobulin (IVIG) in Japanese children, who have the highest incidence of KD, is unknown., Methodology/principal Findings: We used unconditional logistic regression analyses to determine the associations of the copy number of the CCL3L1 gene-containing duplication and CCR2-CCR5 haplotypes in 133 Japanese KD cases [33 with CAL and 25 with resistance to IVIG] and 312 Japanese controls without a history of KD. We observed that the deviation from the population average of four CCL3L1 copies (i.e., four copies) was associated with an increased risk of KD and IVIG resistance (adjusted odds ratio (OR)=2.25, p=0.004 and OR=6.26, p=0.089, respectively). Heterozygosity for the CCR5 HHF*2 haplotype was associated with a reduced risk of both IVIG resistance (OR=0.21, p=0.026) and CAL development (OR=0.44, p=0.071)., Conclusions/significance: The CCL3L1-CCR5 axis may play an important role in KD pathogenesis. In addition to clinical and laboratory parameters, genetic markers may also predict risk of CAL and resistance to IVIG.