The decision to move forward with three clinical trials of IL-1 blockade for treatment of acute Kawasaki disease is a case study in translational science. These trials were born on the one hand from transcriptome studies of host response during the acute disease coupled with animal model investigations of key immune signaling pathways and, on the other hand, out of clinical desperation to intervene in patients with severe inflammation in the setting of acute Kawasaki disease. The convergence of laboratory science and clinical observations led to the clinical trials described here and serves as a model for how such observations can be translated into new therapies.

Kawasaki disease (KD) is a pediatric vasculitis that damages the coronary arteries in 25% of untreated and approximately 5% of treated children. Epidemiologic data suggest that KD is triggered by unidentified infection(s) in genetically susceptible children. To investigate genetic determinants of KD susceptibility, we performed a genome-wide association study (GWAS) in 119 Caucasian KD cases and 135 matched controls with stringent correction for possible admixture, followed by replication in an independent cohort and subsequent fine-mapping, for a total of 893 KD cases plus population and family controls. Significant associations of 40 SNPs and six haplotypes, identifying 31 genes, were replicated in an independent cohort of 583 predominantly Caucasian KD families, with NAALADL2 (rs17531088, p(combined) = 1.13 x 10(-6)) and ZFHX3 (rs7199343, p(combined) = 2.37 x 10(-6)) most significantly associated. Sixteen associated variants with a minor allele frequency of >0.05 that lay within or close to known genes were fine-mapped with HapMap tagging SNPs in 781 KD cases, including 590 from the discovery and replication stages. Original or tagging SNPs in eight of these genes replicated the original findings, with seven genes having further significant markers in adjacent regions. In four genes (ZFHX3, NAALADL2, PPP1R14C, and TCP1), the neighboring markers were more significantly associated than the originally associated variants. Investigation of functional relationships between the eight fine-mapped genes using Ingenuity Pathway Analysis identified a single functional network (p = 10(-13)) containing five fine-mapped genes-LNX1, CAMK2D, ZFHX3, CSMD1, and TCP1-with functional relationships potentially related to inflammation, apoptosis, and cardiovascular pathology. Pair-wise blood transcript levels were measured during acute and convalescent KD for all fine-mapped genes, revealing a consistent trend of significantly reduced transcript levels prior to treatment. This is one of the first GWAS in an infectious disease. We have identified novel, plausible, and functionally related variants associated with KD susceptibility that may also be relevant to other cardiovascular diseases.

Bloodstream infection is a hallmark of sepsis, a medically emergent condition requiring rapid treatment. However, upregulation of host defense proteins through toll-like receptors and NFκB requires hours after endotoxin detection. Using confocal pulmonary intravital microscopy, we identified that the lung provides a TLR4-Myd88-and abl tyrosine kinase-dependent niche for immediate CD11b-dependent neutrophil responses to endotoxin and Gram-negative bloodstream pathogens. In an in vivo model of bacteremia, neutrophils crawled to and rapidly phagocytosed Escherichia coli sequestered to the lung endothelium. Therefore, the lung capillaries provide a vascular defensive niche whereby endothelium and neutrophils cooperate for immediate detection and capture of disseminating pathogens.

Importance

Because clinical features do not reliably distinguish bacterial from viral infection, many children worldwide receive unnecessary antibiotic treatment, while bacterial infection is missed in others.

Objective

To identify a blood RNA expression signature that distinguishes bacterial from viral infection in febrile children.

Design, setting, and participants

Febrile children presenting to participating hospitals in the United Kingdom, Spain, the Netherlands, and the United States between 2009-2013 were prospectively recruited, comprising a discovery group and validation group. Each group was classified after microbiological investigation as having definite bacterial infection, definite viral infection, or indeterminate infection. RNA expression signatures distinguishing definite bacterial from viral infection were identified in the discovery group and diagnostic performance assessed in the validation group. Additional validation was undertaken in separate studies of children with meningococcal disease (n = 24) and inflammatory diseases (n = 48) and on published gene expression datasets.

Exposures

A 2-transcript RNA expression signature distinguishing bacterial infection from viral infection was evaluated against clinical and microbiological diagnosis.

Main outcomes and measures

Definite bacterial and viral infection was confirmed by culture or molecular detection of the pathogens. Performance of the RNA signature was evaluated in the definite bacterial and viral group and in the indeterminate infection group.

Results

The discovery group of 240 children (median age, 19 months; 62% male) included 52 with definite bacterial infection, of whom 36 (69%) required intensive care, and 92 with definite viral infection, of whom 32 (35%) required intensive care. Ninety-six children had indeterminate infection. Analysis of RNA expression data identified a 38-transcript signature distinguishing bacterial from viral infection. A smaller (2-transcript) signature (FAM89A and IFI44L) was identified by removing highly correlated transcripts. When this 2-transcript signature was implemented as a disease risk score in the validation group (130 children, with 23 definite bacterial, 28 definite viral, and 79 indeterminate infections; median age, 17 months; 57% male), all 23 patients with microbiologically confirmed definite bacterial infection were classified as bacterial (sensitivity, 100% [95% CI, 100%-100%]) and 27 of 28 patients with definite viral infection were classified as viral (specificity, 96.4% [95% CI, 89.3%-100%]). When applied to additional validation datasets from patients with meningococcal and inflammatory diseases, bacterial infection was identified with a sensitivity of 91.7% (95% CI, 79.2%-100%) and 90.0% (95% CI, 70.0%-100%), respectively, and with specificity of 96.0% (95% CI, 88.0%-100%) and 95.8% (95% CI, 89.6%-100%). Of the children in the indeterminate groups, 46.3% (63/136) were classified as having bacterial infection, although 94.9% (129/136) received antibiotic treatment.

Conclusions and relevance

This study provides preliminary data regarding test accuracy of a 2-transcript host RNA signature discriminating bacterial from viral infection in febrile children. Further studies are needed in diverse groups of patients to assess accuracy and clinical utility of this test in different clinical settings.

The aetiology of Kawasaki disease (KD), an acute inflammatory disorder of childhood, remains unknown despite various triggers of KD having been proposed. Host 'omic profiles offer insights into the host response to infection and inflammation, with the interrogation of multiple 'omic levels in parallel providing a more comprehensive picture. We used differential abundance analysis, pathway analysis, clustering, and classification techniques to explore whether the host response in KD is more similar to the response to bacterial or viral infections at the transcriptomic and proteomic levels through comparison of 'omic profiles from children with KD to those with bacterial and viral infections. Pathways activated in patients with KD included those involved in anti-viral and anti-bacterial responses. Unsupervised clustering showed that the majority of KD patients clustered with bacterial patients on both 'omic levels, whilst application of diagnostic signatures specific for bacterial and viral infections revealed that many transcriptomic KD samples had low probabilities of having bacterial or viral infections, suggesting that KD may be triggered by a different process not typical of either common bacterial or viral infections. Clustering based on the transcriptomic and proteomic responses during KD revealed three clusters of KD patients on both 'omic levels, suggesting heterogeneity within the inflammatory response during KD. The observed heterogeneity may reflect differences in the host response to a common trigger, or variation dependent on different triggers of the condition.

Introduction

Kawasaki disease (KD) is a pediatric vasculitis that can result in coronary artery aneurysm (CAA) formation, which is a dangerous complication. Treatment with intravenous immunoglobulin (IVIg) significantly decreases the risk of CAA, possibly through competitive binding to Fc-gamma receptors (FcγRs), which reduces the binding of pathological immune complexes. However, ~20% of children have recrudescence of fever and have an increased risk of CAA. Therefore, we aimed to identify genetic markers at the FCGR2/3 locus associated with susceptibility to KD, IVIg resistance, or CAA.

Materials and methods

We investigated the association of single-nucleotide polymorphisms (SNPs) and copy number variations (CNVs) at the FCGR2/3 locus with KD susceptibility, IVIg resistance, and CAA risk using a family-based test (KD susceptibility) and case-control analyses (IVIg resistance and CAA risk) in different cohorts, adding up to a total of 1,167 KD cases. We performed a meta-analysis on IVIg resistance and CAA risk including all cohorts supplemented by previous studies identified through a systematic search.

Results

FCGR2A-p.166His was confirmed to be strongly associated with KD susceptibility (Z = 3.17, p = 0.0015). In case-control analyses, all of the investigated genetic variations at the FCGR2/3 locus were generally not associated with IVIg resistance or with CAA risk, apart from a possible association in a Polish cohort for the FCGR3B-NA2 haplotype (OR = 2.15, 95% CI = 1.15-4.01, p = 0.02). Meta-analyses of all available cohorts revealed no significant associations of the FCGR2/3 locus with IVIg resistance or CAA risk.

Discussion

FCGR2/3 polymorphisms are associated with susceptibility to KD but not with IVIg resistance and CAA formation. Currently known genetic variations at the FCGR2/3 locus are not useful in prediction models for IVIg resistance or CAA risk.

The human Fc-gamma receptors (FcγRs) link adaptive and innate immunity by binding immunoglobulin G (IgG). All human low-affinity FcγRs are encoded by the FCGR2/3 locus containing functional single nucleotide polymorphisms (SNPs) and gene copy number variants. This locus is notoriously difficult to genotype and high-throughput methods commonly used focus on only a few SNPs. We performed multiplex ligation-dependent probe amplification for all relevant genetic variations at the FCGR2/3 locus in >4,000 individuals to define linkage disequilibrium (LD) and allele frequencies in different populations. Strong LD and extensive ethnic variation in allele frequencies was found across the locus. LD was strongest for the FCGR2C-ORF haplotype (rs759550223+rs76277413), which leads to expression of FcγRIIc. In Europeans, the FCGR2C-ORF haplotype showed strong LD with, among others, rs201218628 (FCGR2A-Q27W, r ² = 0.63). LD between these two variants was weaker (r ² = 0.17) in Africans, whereas the FCGR2C-ORF haplotype was nearly absent in Asians (minor allele frequency <0.005%). The FCGR2C-ORF haplotype and rs1801274 (FCGR2A-H131R) were in weak LD (r ² = 0.08) in Europeans. We evaluated the importance of ethnic variation and LD in Kawasaki Disease (KD), an acute vasculitis in children with increased incidence in Asians. An association of rs1801274 with KD was previously shown in ethnically diverse genome-wide association studies. Now, we show in 1,028 European KD patients that the FCGR2C-ORF haplotype, although nearly absent in Asians, was more strongly associated with susceptibility to KD than rs1801274 in Europeans. Our data illustrate the importance of interpreting findings of association studies concerning the FCGR2/3 locus with knowledge of LD and ethnic variation.

Importance

To date, there is no diagnostic test for Kawasaki disease (KD). Diagnosis is based on clinical features shared with other febrile conditions, frequently resulting in delayed or missed treatment and an increased risk of coronary artery aneurysms.

Objective

To identify a whole-blood gene expression signature that distinguishes children with KD in the first week of illness from other febrile conditions.

Design, setting, and participants

The case-control study comprised a discovery group that included a training and test set and a validation group of children with KD or comparator febrile illness. The setting was pediatric centers in the United Kingdom, Spain, the Netherlands, and the United States. The training and test discovery group comprised 404 children with infectious and inflammatory conditions (78 KD, 84 other inflammatory diseases, and 242 bacterial or viral infections) and 55 healthy controls. The independent validation group comprised 102 patients with KD, including 72 in the first 7 days of illness, and 130 febrile controls. The study dates were March 1, 2009, to November 14, 2013, and data analysis took place from January 1, 2015, to December 31, 2017.

Main outcomes and measures

Whole-blood gene expression was evaluated using microarrays, and minimal transcript sets distinguishing KD were identified using a novel variable selection method (parallel regularized regression model search). The ability of transcript signatures (implemented as disease risk scores) to discriminate KD cases from controls was assessed by area under the curve (AUC), sensitivity, and specificity at the optimal cut point according to the Youden index.

Results

Among 404 patients in the discovery set, there were 78 with KD (median age, 27 months; 55.1% male) and 326 febrile controls (median age, 37 months; 56.4% male). Among 202 patients in the validation set, there were 72 with KD (median age, 34 months; 62.5% male) and 130 febrile controls (median age, 17 months; 56.9% male). A 13-transcript signature identified in the discovery training set distinguished KD from other infectious and inflammatory conditions in the discovery test set, with AUC of 96.2% (95% CI, 92.5%-99.9%), sensitivity of 81.7% (95% CI, 60.0%-94.8%), and specificity of 92.1% (95% CI, 84.0%-97.0%). In the validation set, the signature distinguished KD from febrile controls, with AUC of 94.6% (95% CI, 91.3%-98.0%), sensitivity of 85.9% (95% CI, 76.8%-92.6%), and specificity of 89.1% (95% CI, 83.0%-93.7%). The signature was applied to clinically defined categories of definite, highly probable, and possible KD, resulting in AUCs of 98.1% (95% CI, 94.5%-100%), 96.3% (95% CI, 93.3%-99.4%), and 70.0% (95% CI, 53.4%-86.6%), respectively, mirroring certainty of clinical diagnosis.

Conclusions and relevance

In this study, a 13-transcript blood gene expression signature distinguished KD from other febrile conditions. Diagnostic accuracy increased with certainty of clinical diagnosis. A test incorporating the 13-transcript disease risk score may enable earlier diagnosis and treatment of KD and reduce inappropriate treatment in those with other diagnoses.

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.

Kawasaki disease (KD) is a paediatric vasculitis associated with coronary artery aneurysms (CAA). Genetic variants influencing susceptibility to KD have been previously identified, but no risk alleles have been validated that influence CAA formation. We conducted a genome-wide association study (GWAS) for CAA in KD patients of European descent with 200 cases and 276 controls. A second GWAS for susceptibility pooled KD cases with healthy paediatric controls from vaccine trials in the UK (n = 1609). Logistic regression mixed models were used for both GWASs. The susceptibility GWAS was meta-analysed with 400 KD cases and 6101 controls from a previous European GWAS, these results were further meta-analysed with Japanese GWASs at two putative loci. The CAA GWAS identified an intergenic region of chromosome 20q13 with multiple SNVs showing genome-wide significance. The risk allele of the most associated SNV (rs6017006) was present in 13% of cases and 4% of controls; in East Asian 1000 Genomes data, the allele was absent or rare. Susceptibility GWAS with meta-analysis with previously published European data identified two previously associated loci (ITPKC and FCGR2A). Further meta-analysis with Japanese GWAS summary data from the CASP3 and FAM167A genomic regions validated these loci in Europeans showing consistent effects of the top SNVs in both populations. We identified a novel locus for CAA in KD patients of European descent. The results suggest that different genes determine susceptibility to KD and development of CAA and future work should focus on the function of the intergenic region on chromosome 20q13.