Gene transcript abundance profiles distinguish Kawasaki disease from adenovirus infection

HHV-6B detection and host gene expression implicate HHV-6B as pulmonary pathogen after hematopoietic cell transplant

Limited understanding of the immunopathogenesis of human herpesvirus 6B (HHV-6B) has prevented its acceptance as a pulmonary pathogen after hematopoietic cell transplant (HCT). In this prospective multicenter study of patients undergoing bronchoalveolar lavage (BAL) for pneumonia after allogeneic HCT, we test blood and BAL fluid (BALF) for HHV-6B DNA and mRNA transcripts associated with lytic infection and perform RNA-seq on paired blood. Among 116 participants, HHV-6B DNA is detected in 37% of BALs, 49% of which also have HHV-6B mRNA detection. We establish HHV-6B DNA viral load thresholds in BALF that are highly predictive of HHV-6B mRNA detection and associated with increased risk for overall mortality and death from respiratory failure. Participants with HHV-6B DNA in BALF exhibit distinct host gene expression signatures, notable for enriched interferon signaling pathways in participants clinically diagnosed with idiopathic pneumonia. These data implicate HHV-6B as a pulmonary pathogen after allogeneic HCT.

Calculating the fraction of Kawasaki disease potentially attributable to seasonal pathogens: a time series analysis

Background

Kawasaki disease is an acute, febrile, systemic vasculitis of children that primarily affects medium-sized blood vessels with a tropism for the coronary arteries. Although the etiological factors remain unknown, infections have been suggested as the trigger of Kawasaki disease. We sought to calculate the fraction of Kawasaki disease potentially attributable to seasonal infections.

Methods

This cohort study used a population-based time series analysis from the French hospitalisation database (Programme de Médicalisation des Systèmes d'Information), which includes all inpatients admitted to any public or private hospital in France. We included all children aged 0-17 years hospitalised for Kawasaki disease in France over 13 years. The monthly incidence of Kawasaki disease per 10,000 children over time was analysed by a quasi-Poisson regression model. The model accounted for seasonality by using harmonic terms (a pair of sines and cosines with 12-month periods). The circulation of eight common seasonal pathogens (adenovirus, influenza, metapneumovirus, Mycoplasma pneumoniae, norovirus, rhinovirus, rotavirus, respiratory syncytial virus, and Streptococcus pneumonia) over the same period was included in the model to analyse the fraction of Kawasaki disease potentially attributable to each pathogen. Infections were identified on the basis of polymerase chain reaction or rapid antigen testing in hospital laboratories.

Findings

Between Jan 1, 2007, and Dec 31, 2019, we included 10,337 children with Kawasaki disease and 442,762 children with the selected infectious diseases. In the Kawasaki disease cohort, the median age [IQR] was 2 [0-4] years, 6164 [59.6%] were boys. Adenovirus infection was potentially responsible for 24.4% [21.5-27.8] (p < 0.001) of Kawasaki diseases, Norovirus for 6.7% [1.3-11.2] (p = 0.002), and RSV 4.6% [1.2-7.8] (p = 0.022). Sensitivity analyses found similar results.

Interpretation

This cohort study of data from a comprehensive national hospitalisation database indicated that approximately 35% of Kawasaki diseases was potentially attributable to seasonal infections.

Funding

None.

Neutrophil to Lymphocyte Ratio as a Biomarker for Predicting the Coronary Artery Abnormality in Kawasaki Disease: A Meta-Analysis

We conducted a systematic review and meta-analysis on the relationship between the neutrophil to lymphocyte ratio (NLR) and coronary artery abnormalities (CAA) in patients with Kawasaki disease (KD), according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) statements. We searched PubMed, Scopus, Web of Science, Embase, TRIP, Google Scholar, and ProQuest up to the 8th of August 2022. This was done to retrieve eligible studies. No date or language limitations were considered in this study. Methodology quality assessment was conducted according to the Newcastle-Ottawa scale (NOS). Standard mean difference (SMD) and its 95% confidence interval (CI) were used to depict the pooled continuous variables. Finally, 17 articles with 6334 KD patients, of whom 1328 developed CAA, were enrolled in this meta-analysis. NLR level was significantly higher in KD patients with CAA compared to those without CAA (SMD =0.81; 95% CI =0.05-1.57, P = 0.03). In addition, NLR level was significantly higher in patients with coronary artery aneurysms than those without coronary artery aneurysms (SMD =2.29; 95% CI =0.18-4.41, P = 0.03). However, no significant association between NLR and coronary artery dilation was observed in this meta-analysis (SMD =0.56; 95% CI = -0.86-1.99). There was no publication bias for the pooled SMD of NLR for coronary artery abnormality in KD (Egger's test P = 0.82; Begg's test P = 0.32). The NLR may be useful in monitoring CAA development in these patients and may further imply a mechanistic role in potential inflammation that mediates this process.

Innate immune response analysis in COVID-19 and kawasaki disease reveals MIS-C predictors

BACKGROUND/PURPOSE

The association between dysregulated innate immune responses seen in Kawasaki disease (KD) with predisposition to Kawasaki-like multisystem inflammatory syndrome in children (MIS-C) remains unclear. We aimed to compare the innate immunity transcriptome signature between COVID-19 and KD, and to analyze the interactions of these molecules with genes known to predispose to KD.

METHODS

Transcriptome datasets of COVID-19 and KD cohorts (E-MTAB-9357, GSE-63881, GSE-68004) were downloaded from ArrayExpress for innate immune response analyses. Network analysis was used to determine enriched pathways of interactions.

RESULTS

Upregulations of IRAK4, IFI16, STING, STAT3, PYCARD, CASP1, IFNAR1 and CD14 genes were observed in blood cells of acute SARS-CoV-2 infections with moderate severity. In the same patient group, increased expressions of TLR2, TLR7, IRF3, and CD36 were also noted in blood drawn a few days after COVID-19 diagnosis. Elevated blood PYCARD level was associated with severe COVID-19 in adults. Similar gene expression signature except differences in TLR8, NLRP3, STING and IRF3 levels was detected in KD samples. Network analysis on innate immune genes and genes associated with KD susceptibility identified enriched pathways of interactions. Furthermore, higher expression levels of KD susceptibility genes HLA-DOB, PELI1 and FCGR2A correlated with COVID-19 of different severities.

CONCLUSION

Our findings suggest that most enriched innate immune response pathways were shared between transcriptomes of KD and COVID-19 with moderate severity. Genetic polymorphisms associated with innate immune dysregulation and KD susceptibility, together with variants in STING and STAT3, might predict COVID-19 severity and potentially susceptibility to COVID-19 related MIS-C.

Atlas of circulating immune cells in Kawasaki disease

Increasing evidence shows that the pathogenesis of Kawasaki disease (KD) is caused by abnormal and unbalanced innate and adaptive immune responses. However, the changes in and functions of adaptive immune cells in the peripheral blood of subjects with KD remain controversial. In this study, three different methods, CIBERSORT, Immune Cell Abundance Identifier (ImmuCellAI), and immune cell markers, were used to evaluate the proportions and abundances of immune cells in eight KD datasets (GSE9863, GSE9864, GSE18606, GSE63881, GSE68004, GSE73461, GSE73463, and GSE64486; a total of 1,251 samples). Compared with those in normal controls and convalescent KD samples, the proportions and abundances of innate immune cells such as neutrophils, monocytes, and macrophages in acute KD peripheral blood samples were significantly increased, while those of adaptive immune cells such as B and T cells were significantly decreased. The change tendencies of these immune cells were similar to those observed in other febrile illnesses but were more significant. However, in the coronary artery tissues of patients with convalescent KD, adaptive immune cells, especially B cells and CD8+ T cell subsets, were significantly increased. This result suggests that adaptive immune cells can be selectively recruited from peripheral blood into the coronary arteries. In addition, we found that elevated neutrophils in peripheral blood could be used as a biomarker to assist in the differential diagnosis of KD, but we did not find immune cells that could accurately predict intravenousimmunoglobulin (IVIG) responses in multiple datasets.

Single-cell RNA sequencing of peripheral blood mononuclear cells from acute Kawasaki disease patients

Kawasaki disease (KD) is the most common cause of acquired heart disease in children in developed countries. Although functional and phenotypic changes of immune cells have been reported, a global understanding of immune responses underlying acute KD is unclear. Here, using single-cell RNA sequencing, we profile peripheral blood mononuclear cells from seven patients with acute KD before and after intravenous immunoglobulin therapy and from three age-matched healthy controls. The most differentially expressed genes are identified in monocytes, with high expression of pro-inflammatory mediators, immunoglobulin receptors and low expression of MHC class II genes in acute KD. Single-cell RNA sequencing and flow cytometry analyses, of cells from an additional 16 KD patients, show that although the percentage of total B cells is substantially decreased after therapy, the percentage of plasma cells among the B cells is significantly increased. The percentage of CD8+ T cells is decreased in acute KD, notably effector memory CD8+ T cells compared with healthy controls. Oligoclonal expansions of both B cell receptors and T cell receptors are observed after therapy. We identify biological processes potentially underlying the changes of each cell type. The single-cell landscape of both innate and adaptive immune responses provides insights into pathogenesis and therapy of KD.

A Diagnostic Model for Kawasaki Disease Based on Immune Cell Characterization From Blood Samples

Background: Kawasaki disease (KD) is the leading cause of acquired heart disease in children. However, distinguishing KD from febrile infections early in the disease course remains difficult. Our goal was to estimate the immune cell composition in KD patients and febrile controls (FC), and to develop a tool for KD diagnosis. Methods: We used a machine-learning algorithm, CIBERSORT, to estimate the proportions of 22 immune cell types based on blood samples from children with KD and FC. Using these immune cell compositions, a diagnostic score for predicting KD was then constructed based on LASSO regression for binary outcomes. Results: In the training set (n = 496), a model was fit which consisted of eight types of immune cells. The area under the curve (AUC) values for diagnosing KD in a held-out test set (n = 212) and an external validation set (n = 36) were 0.80 and 0.77, respectively. The most common cell types in KD blood samples were monocytes, neutrophils, CD4⁺-naïve and CD8⁺ T cells, and M0 macrophages. The diagnostic score was highly correlated to genes that had been previously reported as associated with KD, such as interleukins and chemokine receptors, and enriched in reported pathways, such as IL-6/JAK/STAT3 and TNFα signaling pathways. Conclusion: Altogether, the diagnostic score for predicting KD could potentially serve as a biomarker. Prospective studies could evaluate how incorporating the diagnostic score into a clinical algorithm would improve diagnostic accuracy further.

Severe COVID-19, multisystem inflammatory syndrome in children, and Kawasaki disease: immunological mechanisms, clinical manifestations and management

Multisystem inflammatory syndrome (MIS-C) is a pediatric hyperinflammation disorder caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). It has now been reported from several countries the world over. Some of the clinical manifestations of MIS-C mimic Kawasaki disease (KD) shock syndrome. MIS-C develops 4–6 weeks following SARS-CoV-2 infection, and is presumably initiated by adaptive immune response. Though it has multisystem involvement, it is the cardiovascular manifestations that are most prominent. High titres of anti-SARS-CoV-2 antibodies are seen in these patients. As this is a new disease entity, its immunopathogenesis is not fully elucidated. Whether it has some overlap with KD is still unclear. Current treatment guidelines recommend use of intravenous immunoglobulin and high-dose corticosteroids as first-line treatment. Mortality rates of MIS-C are lower compared to adult forms of severe COVID-19 disease.

Distinguishing Kawasaki Disease from Febrile Infectious Disease Using Gene Pair Signatures

Kawasaki disease (KD) is an acute systemic vasculitis of childhood with prolonged fever, and the diagnosis of KD is mainly based on clinical criteria, which is prone to misdiagnosis with other febrile infectious (FI) diseases. Currently, there remain no effective molecular markers for KD diagnosis. In this study, we aimed to use a relative-expression-based method k-TSP and resampling framework to identify robust gene pair signatures to distinguish KD from bacterial and virus febrile infectious diseases. Our study pool consisted of 808 childhood patients from several studies and assigned to three groups, namely, the discovery set (n = 224), validation set-1 (n = 197), and validation set-2 (n = 387). We had identified 60 biologically relevant gene pairs and developed a top-ranked gene pair classifier (TRGP) using the first seven signatures, with the area under the receiver-operating characteristic curves (AUROC) of 0.947 (95% CI, 0.918-0.976), a sensitivity of 0.936 (95% CI, 0.872-0.987), and a specificity of 0.774 (95% CI, 0.705-0.836) in the discovery set. In the validation set-1, the TRGP classifier distinguished KD from FI with AUROC of 0.955 (95% CI, 0.919-0.991), a sensitivity of 0.959 (95% CI, 0.925-0.986), and a specificity of 0.863 (95% CI, 0.764-0.961). In the validation set-2, the predictive performance of classification was with an AUROC of 0.796 (95% CI, 0.747-0.845), a sensitivity of 0.797 (95% CI, 0.720-0.864), and a specificity of 0.661 (95% CI, 0.606-0.717). Our study reveals that gene pair signatures are robust across diverse studies and can be utilized as objective biomarkers to distinguish KD from FI, helping to develop a fast, simple, and effective molecular approach to improve the diagnosis of KD.

Kawasaki Disease: Global Burden and Genetic Background

Kawasaki disease (KD) is a childhood vasculitides associated with serious coronary artery lesions. It is the most common cause of pediatric acquired heart disease in developed countries, and is increasingly reported from many rapidly industrializing developing countries. The incidence varies widely among different nations and is highest in North-East Asian countries, where almost 1 in 100 children in Japan having the disease by age of 5, where the lowest incidence reported in sub-Saharan Africa. The etiology of KD is still uncertain; interaction between a genetic predisposition and several environmental and immunological factors has been hypothesized. Several susceptibility genes were identified to be associated with the development of KD and increased risk of coronary artery lesions. Gene-gene associations and alteration of deoxyribonucleic acid (DNA) methylation are also found to play key roles in the pathogenesis and prognosis of KD. This article will focus on the global epidemiological patterns of KD, and the currently known genetic predisposition.

The use of host factors in microbial forensics

Advances have been made in the forensic analysis of microbes and toxins. An underdeveloped and underutilized area in microbial forensics is how the host interacts with microorganisms in a way that provides unique signatures for forensic use. For forensic purposes, an immediate goal is to distinguish a potential victim and innocent person from a perpetrator, and to distinguish between a naturally acquired or intentional infection. Principal methods that are sufficiently developed are characterization of the humoral immune response to microbial antigens including vaccine-induced immunity and detection of antibiotics that may be present in a possible perpetrator. This chapter presents central elements of the host response in a simplified fashion and describes a representative example, which, in the appropriate context, has a high potential of providing evidence that may aid an investigation to distinguish a perpetrator from a victim. This chapter also presents information about the immune system so that the interested reader can have a fuller understanding of the immune response in general.

Aetiological Significance of Infectious Stimuli in Kawasaki Disease

Kawasaki disease (KD) is a pediatric vasculitis syndrome that is often involves coronary artery lesions (e. g., coronary artery aneurysms). Although its causal factors and entire pathogenesis remain elusive, the available evidence indicates that the pathogenesis of KD is closely associated with dysregulation of immune responses to various viruses or microbes. In this short review, we address several essential aspects of the etiology of KD with respect to the immune response to infectious stimuli: 1) the role of viral infections, 2) the role of bacterial infections and the superantigen hypothesis, 3) involvement of innate immune response including pathogens/microbe-associated molecular patterns and complement pathways, and 4) the influence of genetic background on the response to infectious stimuli. Based on the clinical and experimental evidence, we discuss the possibility that a wide range of microbes and viruses could cause KD through common and distinct immune processes.

How Should We Classify Kawasaki Disease?

The exact classification of Kawasaki disease (KD) has been debated. Infectious disease specialists have claimed it as an infection with a classic immune responses to an as yet unidentified pathogen that localizes to the coronary arteries. Others have favored an autoreactive hypothesis that KD is triggered by an antigen that shares homology with structures in the vascular wall, and molecular mimicry resulting in an immune response directed to that tissue. Rheumatologists have classified it as a systemic vasculitis, while some immunologists have stressed the robust nature of the innate immune response that causes both systemic inflammation as well as damage to the coronary arterial wall and questioned whether KD falls within the spectrum of autoinflammatory diseases. This review will describe the evidences available up to now regarding these hypotheses.

Diagnosis of Kawasaki Disease Using a Minimal Whole-Blood Gene Expression Signature

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.

Whole blood transcriptional profiles as a prognostic tool in complete and incomplete Kawasaki Disease

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.

Robust classification of bacterial and viral infections via integrated host gene expression diagnostics

Improved diagnostics for acute infections could decrease morbidity and mortality by increasing early antibiotics for patients with bacterial infections and reducing unnecessary antibiotics for patients without bacterial infections. Several groups have used gene expression microarrays to build classifiers for acute infections, but these have been hampered by the size of the gene sets, use of overfit models, or lack of independent validation. We used multicohort analysis to derive a set of seven genes for robust discrimination of bacterial and viral infections, which we then validated in 30 independent cohorts. We next used our previously published 11-gene Sepsis MetaScore together with the new bacterial/viral classifier to build an integrated antibiotics decision model. In a pooled analysis of 1057 samples from 20 cohorts (excluding infants), the integrated antibiotics decision model had a sensitivity and specificity for bacterial infections of 94.0 and 59.8%, respectively (negative likelihood ratio, 0.10). Prospective clinical validation will be needed before these findings are implemented for patient care.

Is Kawasaki disease an infectious disorder?

Although the etiology of Kawasaki disease (KD) is largely unknown, a large body of clinical, epidemiologic, immunologic, pathologic and ultrastructural evidence suggests that an infectious agent triggers a cascade that causes the illness. However, this elusive infectious agent remains unidentified at present. Increasingly sensitive molecular methods for identifying microbial nucleic acids and proteins in tissue samples continue to rapidly emerge, and these methods should be utilized in studies on KD etiology as they become available. Identifying the etiology of this enigmatic disease remains the single most important research goal in the field, and accomplishing this goal is the best means to improve diagnosis, treatment and prevention of this potentially fatal childhood disease.

Pediatric Kawasaki Disease and Adult Human Immunodeficiency Virus Kawasaki-Like Syndrome Are Likely the Same Malady

Background.  Pediatric Kawasaki disease (KD) and human immunodeficiency virus (HIV)⁺ adult Kawasaki-like syndrome (KLS) are dramatic vasculitides with similar physical findings. Both syndromes include unusual arterial histopathology with immunoglobulin (Ig)A⁺ plasma cells, and both impressively respond to pooled Ig therapy. Their distinctive presentations, histopathology, and therapeutic response suggest a common etiology. Because blood is in immediate contact with inflamed arteries, we investigated whether KD and KLS share an inflammatory signature in serum. Methods.  A custom multiplex enzyme-linked immunosorbent assay (ELISA) defined the serum cytokine milieu in 2 adults with KLS during acute and convalescent phases, with asymptomatic HIV⁺ subjects not taking antiretroviral therapy serving as controls. We then prospectively collected serum and plasma samples from children hospitalized with KD, unrelated febrile illnesses, and noninfectious conditions, analyzing them with a custom multiplex ELISA based on the KLS data. Results.  Patients with KLS and KD subjects shared an inflammatory signature including acute-phase reactants reflecting tumor necrosis factor (TNF)-α biologic activity (soluble TNF receptor I/II) and endothelial/smooth muscle chemokines Ccl1 (Th2), Ccl2 (vascular inflammation), and Cxcl11 (plasma cell recruitment). Ccl1 was specifically elevated in KD versus febrile controls, suggesting a unique relationship between Ccl1 and KD/KLS pathogenesis. Conclusions.  This study defines a KD/KLS inflammatory signature mirroring a dysfunctional response likely to a common etiologic agent. The KD/KLS inflammatory signature based on elevated acute-phase reactants and specific endothelial/smooth muscle chemokines was able to identify KD subjects versus febrile controls, and it may serve as a practicable diagnostic test for KD.

The transcriptional profile of coronary arteritis in Kawasaki disease

Background

Kawasaki Disease (KD) can cause potentially life-threatening coronary arteritis in young children, and has a likely infectious etiology. Transcriptome profiling is a powerful approach to investigate gene expression in diseased tissues. RNA sequencing of KD coronary arteries could elucidate the etiology and the host response, with the potential to improve KD diagnosis and/or treatment.

Methods

Deep RNA sequencing was performed on KD (n = 8) and childhood control (n = 7) coronary artery tissues, revealing 1074 differentially expressed mRNAs. Non-human RNA sequences were subjected to a microbial discovery bioinformatics platform, and microbial sequences were analyzed by Metastats for association with KD.

Results

T lymphocyte activation, antigen presentation, immunoglobulin production, and type I interferon response were significantly upregulated in KD arteritis, while the tumor necrosis factor α pathway was not differentially expressed. Transcripts from known infectious agents were not specifically associated with KD coronary arteritis.

Conclusions

The immune transcriptional profile in KD coronary artery tissues has features of an antiviral immune response such as activated cytotoxic T lymphocyte and type I interferon-induced gene upregulation. These results provide new insights into the pathogenesis of KD arteritis that can guide selection of new immunomodulatory therapies for high-risk KD patients, and provide direction for future etiologic studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2323-5) contains supplementary material, which is available to authorized users.

Post-infectious group A streptococcal autoimmune syndromes and the heart

There is a pressing need to reduce the high global disease burden of rheumatic heart disease (RHD) and its harbinger, acute rheumatic fever (ARF). ARF is a classical example of an autoimmune syndrome and is of particular immunological interest because it follows a known antecedent infection with group A streptococcus (GAS). However, the poorly understood immunopathology of these post-infectious diseases means that, compared to much progress in other immune-mediated diseases, we still lack useful biomarkers, new therapies or an effective vaccine in ARF and RHD. Here, we summarise recent literature on the complex interaction between GAS and the human host that culminates in ARF and the subsequent development of RHD. We contrast ARF with other post-infectious streptococcal immune syndromes - post-streptococcal glomerulonephritis (PSGN) and the still controversial paediatric autoimmune neuropsychiatric disorders associated with streptococcal infections (PANDAS), in order to highlight the potential significance of variations in the host immune response to GAS. We discuss a model for the pathogenesis of ARF and RHD in terms of current immunological concepts and the potential for application of in depth "omics" technologies to these ancient scourges.

CXCL10/IP-10 is a biomarker and mediator for Kawasaki disease

RATIONALE

Kawasaki disease (KD), an acute febrile vasculitis, is the most common cause of acquired heart disease in childhood; however, diagnosing KD can be difficult.

OBJECTIVE

To identify unique proteomic biomarkers that can be used to facilitate earlier diagnosis of KD.

METHODS AND RESULTS

We enrolled 214 children with fever and clinical features suggestive of KD. Of those, only 100 were diagnosed with KD. Their plasma samples were globally analyzed for cytokines, chemokines, and cell adhesion molecules using an unbiased, large-scale, quantitative protein array. This study was conducted in 3 stages: discovery, replication, and blinded validation. During the discovery phase (n [KD]=37; n [control]=20), the expression of interleukin-17F, sCD40L, E-selectin, CCL23 (myeloid progenitor inhibitory factor 1), and CXCL10 (IFN-γ-inducible protein 10 [IP-10]) were upregulated during the acute phase in patients with KD when compared with that in the controls. A notable increase was observed in the IP-10 levels (KD, 3037 ± 226.7 pg/mL; control, 672 ± 130.4 pg/mL; P=4.1 × 10(-11)). Receiver-operating characteristic analysis of the combined discovery and replication data (n [KD]=77; n [control]=77) showed that the IP-10 level had high area under the curve values (0.94 [95% confidence interval, 0.9055-0.9778]; sensitivity, 100%; and specificity, 77%). With 1318 pg/mL as the optimal cutoff, the blinded validation study confirmed that the IP-10 levels were a good predictor of KD. With intravenous immunoglobulin treatment, the IP-10 levels returned to normal. The downstream receptor of IP-10, CXCR3, was activated in the T cells of patients with acute KD.

CONCLUSIONS

IP-10 may be used as a biomarker to facilitate KD diagnosis, and it may provide clues about the pathogenesis of KD.

The application of transcriptional blood signatures to enhance our understanding of the host response to infection: the example of tuberculosis

Despite advances in antimicrobials, vaccination and public health measures, infectious diseases remain a leading cause of morbidity and mortality worldwide. With the increase in antimicrobial resistance and the emergence of new pathogens, there remains a need for new and more accurate diagnostics, the ability to monitor adequate treatment response as well as the ability to predict prognosis for an individual. Transcriptional approaches using blood signatures have enabled a better understanding of the host response to diseases, leading not only to new avenues of basic research, but also to the identification of potential biomarkers for use in diagnosis, prognosis and treatment monitoring.

Evaluation of gene expression classification studies: factors associated with classification performance

Classification methods used in microarray studies for gene expression are diverse in the way they deal with the underlying complexity of the data, as well as in the technique used to build the classification model. The MAQC II study on cancer classification problems has found that performance was affected by factors such as the classification algorithm, cross validation method, number of genes, and gene selection method. In this paper, we study the hypothesis that the disease under study significantly determines which method is optimal, and that additionally sample size, class imbalance, type of medical question (diagnostic, prognostic or treatment response), and microarray platform are potentially influential. A systematic literature review was used to extract the information from 48 published articles on non-cancer microarray classification studies. The impact of the various factors on the reported classification accuracy was analyzed through random-intercept logistic regression. The type of medical question and method of cross validation dominated the explained variation in accuracy among studies, followed by disease category and microarray platform. In total, 42% of the between study variation was explained by all the study specific and problem specific factors that we studied together.

Gene expression profiles in febrile children with defined viral and bacterial infection

Viral infections are common causes of fever without an apparent source in young children. Despite absence of bacterial infection, many febrile children are treated with antibiotics. Virus and bacteria interact with different pattern recognition receptors in circulating blood leukocytes, triggering specific host transcriptional programs mediating immune response. Therefore, unique transcriptional signatures may be defined that discriminate viral from bacterial causes of fever without an apparent source. Gene expression microarray analyses were conducted on blood samples from 30 febrile children positive for adenovirus, human herpesvirus 6, or enterovirus infection or with acute bacterial infection and 22 afebrile controls. Blood leukocyte transcriptional profiles clearly distinguished virus-positive febrile children from both virus-negative afebrile controls and afebrile children with the same viruses present in the febrile children. Virus-specific gene expression profiles could be defined. The IFN signaling pathway was uniquely activated in febrile children with viral infection, whereas the integrin signaling pathway was uniquely activated in children with bacterial infection. Transcriptional profiles classified febrile children with viral or bacterial infection with better accuracy than white blood cell count in the blood. Similarly accurate classification was shown with data from an independent study using different microarray platforms. Our results support the paradigm of using host response to define the etiology of childhood infections. This approach could be an important supplement to highly sensitive tests that detect the presence of a possible pathogen but do not address its pathogenic role in the patient being evaluated.

Surface-enhanced Raman scattering molecular sentinel nanoprobes for viral infection diagnostics

In this paper, we describe a surface-enhanced Raman scattering (SERS)-based detection approach, referred to as "molecular sentinel" (MS) plasmonic nanoprobes, to detect an RNA target related to viral infection. The MS method is essentially a label-free technique incorporating the SERS effect modulation scheme associated with silver nanoparticles and Raman dye-labeled DNA hairpin probes. Hybridization with target sequences opens the hairpin and spatially separates the Raman label from the silver surface thus reducing the SERS signal of the label. Herein, we have developed a MS nanoprobe to detect the human radical S-adenosyl methionine domain containing 2 (RSAD2) RNA target as a model system for method demonstration. The human RSAD2 gene has recently emerged as a novel host-response biomarker for diagnosis of respiratory infections. Our results showed that the RSAD2 MS nanoprobes exhibits high specificity and can detect as low as 1 nM target sequences. With the use of a portable Raman spectrometer and total RNA samples, we have also demonstrated for the first time the potential of the MS nanoprobe technology for detection of host-response RNA biomarkers for infectious disease diagnostics.

Can a systems biology approach unlock the mysteries of Kawasaki disease?

Kawasaki disease (KD) is a systemic inflammatory illness of childhood that particularly affects the coronary arteries. It can lead to coronary artery aneurysms, myocardial infarction, and sudden death. Clinical and epidemiologic data support an infectious cause, and the etiology remains unknown, but recent data support infection with a 'new' virus. Genetic factors influence KD susceptibility; the incidence is 10-fold higher in children of Asian when compared with Caucasian ethnicity. Recent research has identified genes affecting immune response that are associated with KD susceptibility and outcome. A re-examination of the pathologic features of KD has yielded a three process model of KD vasculopathy, providing a framework for understanding the KD arterial immune response and the damage it inflicts and for identifying new therapeutic targets for KD patients with coronary artery abnormalities. The researcher is faced with many challenges in determining the pathogenesis of KD. A systems biology approach incorporating genomics, proteomics, transcriptomics, and microbial bioinformatics analysis of high-throughput sequence data from KD tissues could provide the keys to unlocking the mysteries of this potentially fatal illness of childhood.

Temporal dynamics of the transcriptional response to dengue virus infection in Nicaraguan children

Dengue is the most prevalent mosquito-borne human illness worldwide. The ability to predict disease severity during the earliest days of the illness is a long-sought, but unachieved goal.

We examined human genome-wide transcript abundance patterns in daily peripheral blood mononuclear cell (PBMC) samples from 41 children hospitalized with dengue virus (DENV) infection in Nicaragua, as well as 8 healthy control subjects. Nine patients had primary dengue fever (DF1), 11 had dengue fever with serologic evidence of prior DENV infection, i.e., secondary dengue fever (DF2), 12 had dengue hemorrhagic fever (DHF), and 9 had dengue shock syndrome (DSS). We identified 2,092 genes for which transcript abundance differed significantly between patients on days 3–6 of fever and healthy subjects (FDR<1%). Prior DENV infection explained the greatest amount of variation in gene expression among patients. The number of differentially expressed genes was greatest on fever day 3 in patients with DF1, while the number in patients with DF2 or DHF/DSS was greatest on day 5. Genes associated with the mitotic cell cycle and B cell differentiation were expressed at higher levels, and genes associated with signal transduction and cell adhesion were expressed at lower levels, in patients versus healthy controls. On fever day 3, a set of interferon-stimulated gene transcripts was less abundant in patients who subsequently developed DSS than in other patient groups (p<0.05, ranksum). Patients who later developed DSS also had higher levels of transcripts on day 3 associated with mitochondrial function (p<0.01, ranksum). These day 3 transcript abundance findings were not evident on subsequent fever days.

In conclusion, we identified differences in the timing and magnitude of human gene transcript abundance changes in DENV patients that were associated with serologic evidence of prior infection and with disease severity. Some of these differential features may predict the outcome of DENV infection.

Infection with dengue virus (DENV) causes dengue fever, the most prevalent mosquito-borne illness of humans worldwide. Tens of millions of cases occur annually; up to 500,000 patients develop additional life-threatening complications, including hemorrhage and shock. The clinical course of the disease evolves rapidly, making it difficult to identify patients at risk for severe disease and suggesting that biological events associated with the development of severe disease may be short-lived. We examined gene expression patterns in the blood of children hospitalized with DENV infection, and found that patients with differences in disease severity and history of previous DENV infection shared a common set of gene expression features, but the timing and magnitude of these features differed. In our study, prior DENV infection explained the greatest amount of variation in gene expression among patients. We discovered features of gene expression on day 3 that were associated with subsequent disease severity—features that were not apparent on subsequent days, emphasizing the importance of looking at temporal patterns of gene expression in acute infection.

A diagnostic algorithm combining clinical and molecular data distinguishes Kawasaki disease from other febrile illnesses

BACKGROUND

Kawasaki disease is an acute vasculitis of infants and young children that is recognized through a constellation of clinical signs that can mimic other benign conditions of childhood. The etiology remains unknown and there is no specific laboratory-based test to identify patients with Kawasaki disease. Treatment to prevent the complication of coronary artery aneurysms is most effective if administered early in the course of the illness. We sought to develop a diagnostic algorithm to help clinicians distinguish Kawasaki disease patients from febrile controls to allow timely initiation of treatment.

METHODS

Urine peptidome profiling and whole blood cell type-specific gene expression analyses were integrated with clinical multivariate analysis to improve differentiation of Kawasaki disease subjects from febrile controls.

RESULTS

Comparative analyses of multidimensional protein identification using 23 pooled Kawasaki disease and 23 pooled febrile control urine peptide samples revealed 139 candidate markers, of which 13 were confirmed (area under the receiver operating characteristic curve (ROC AUC 0.919)) in an independent cohort of 30 Kawasaki disease and 30 febrile control urine peptidomes. Cell type-specific analysis of microarrays (csSAM) on 26 Kawasaki disease and 13 febrile control whole blood samples revealed a 32-lymphocyte-specific-gene panel (ROC AUC 0.969). The integration of the urine/blood based biomarker panels and a multivariate analysis of 7 clinical parameters (ROC AUC 0.803) effectively stratified 441 Kawasaki disease and 342 febrile control subjects to diagnose Kawasaki disease.

CONCLUSIONS

A hybrid approach using a multi-step diagnostic algorithm integrating both clinical and molecular findings was successful in differentiating children with acute Kawasaki disease from febrile controls.

Modulation of the host interferon response and ISGylation pathway by B. pertussis filamentous hemagglutinin

Bordetella pertussis filamentous hemagglutinin (FHA) is a surface-associated and secreted protein that serves as a crucial adherence factor, and displays immunomodulatory activity in human peripheral blood mononuclear cells (PBMCs). In order to appreciate more fully the role of secreted FHA in pathogenesis, we analyzed FHA-induced changes in genome-wide transcript abundance in human PBMCs. Among the 683 known unique genes with greater than 3-fold change in transcript abundance following FHA treatment, 125 (18.3%) were identified as interferon (IFN)-regulated. Among the latter group were genes encoding several members of the IFN type I response, as well as 3 key components of the ISGylation pathway. Using real-time RT-PCR, we confirmed FHA-associated increases in transcript abundance for the genes encoding ubiquitin-like protein, ISG15, and its specific protease USP18. Western-blot analysis demonstrated the presence of both, free ISG15 and several ISGylated conjugates in FHA-stimulated PBMC lysates, but not in unstimulated cells. Intracellular FACS analysis provided evidence that monocytes and a natural killer-enriched cell population were the primary producers of ISG15 in PBMCs after FHA stimulation. Our data reveal previously-unrecognized effects of B. pertussis FHA on host IFN and ISGylation responses, and suggest previously-unsuspected mechanisms by which FHA may alter the outcome of the host-pathogen interaction.

Assessing the human immune system through blood transcriptomics

Blood is the pipeline of the immune system. Assessing changes in transcript abundance in blood on a genome-wide scale affords a comprehensive view of the status of the immune system in health and disease. This review summarizes the work that has used this approach to identify therapeutic targets and biomarker signatures in the field of autoimmunity and infectious disease. Recent technological and methodological advances that will carry the blood transcriptome research field forward are also discussed.

Transcript abundance patterns in Kawasaki disease patients with intravenous immunoglobulin resistance

Intravenous immunoglobulin (IVIG)-resistant Kawasaki disease (KD) patients comprise at least 20% of treated patients and are at high risk for coronary artery abnormalities. If identified early in the course of the disease, such patients may benefit from additional anti-inflammatory therapy. The aim of this study was to compare the transcript abundance between IVIG resistant and -responsive KD patients, to identify biomarkers that might differentiate between these two groups and to generate new targets for therapies in IVIG resistant KD patients. We compared the transcript abundance profiles of whole-blood RNA on Agilent arrays from acute and convalescent KD subjects and age-similar, healthy controls. KD subjects were stratified as IVIG resistant or -responsive based on response to initial IVIG therapy. Transcript abundance was higher for IL-1 pathway genes (IL-1 receptor, interleukin receptor associated kinase, p38 mitogen-activated protein kinase), and MMP-8. These findings point to candidate biomarkers that may predict IVIG resistance in acute KD patients. The results also underscore the importance of the IL-1 pathway as a mediator of inflammation in KD and suggest that IL-1 or its receptor may be reasonable targets for therapy, particularly for IVIG resistant patients.

Transcription of proteinase 3 and related myelopoiesis genes in peripheral blood mononuclear cells of patients with active Wegener's granulomatosis

OBJECTIVE

Wegener's granulomatosis (WG) is a systemic inflammatory disease that is associated with substantial morbidity. The aim of this study was to understand the biology underlying WG and to discover markers of disease activity that would be useful for prognosis and treatment guidance.

METHODS

Gene expression profiling was performed using total RNA from peripheral blood mononuclear cells (PBMCs) and granulocyte fractions from 41 patients with WG and 23 healthy control subjects. Gene set enrichment analysis (GSEA) was performed to search for candidate WG-associated molecular pathways and disease activity biomarkers. Principal components analysis was used to visualize relationships between subgroups of WG patients and controls. Longitudinal changes in proteinase 3 (PR3) gene expression were evaluated using reverse transcription-polymerase chain reaction, and clinical outcomes, including remission status and disease activity, were determined using the Birmingham Vasculitis Activity Score for WG (BVAS-WG).

RESULTS

Eighty-six genes in WG PBMCs and 40 in WG polymorphonuclear neutrophils (PMNs) were significantly up-regulated relative to controls. Genes up-regulated in WG PBMCs were involved in myeloid differentiation, and these included the WG autoantigen PR3. The coordinated regulation of myeloid differentiation genes was confirmed by GSEA. The median expression values of the 86 up-regulated genes in WG PBMCs were associated with disease activity (P = 1.3 x 10(-4)), and WG patients with low-level expression of the WG signature genes showed expression profiles that were only modestly different from that in healthy controls (P = 0.07). PR3 transcription was significantly up-regulated in WG PBMCs (P = 1.3 x 10(-5), false discovery rate [FDR] 0.002), but not in WG PMNs (P = 0.03, FDR 0.28), and a preliminary longitudinal analysis showed that the fold change in PR3 RNA levels in WG PBMCs corresponded to changes in the BVAS-WG score over time.

CONCLUSION

Transcription of PR3 and related myeloid differentiation genes in PBMCs may represent novel markers of disease activity in WG.

A genomic approach to human autoimmune diseases

The past decade has seen an explosion in the use of DNA-based microarrays. These techniques permit assessment of RNA abundance on a genome-wide scale. Medical applications emerged in the field of cancer, with studies of both solid tumors and hematological malignancies leading to the development of tests that are now used to personalize therapeutic options. Microarrays have also been used to analyze the blood transcriptome in a wide range of diseases. In human autoimmune diseases, these studies are showing potential for identifying therapeutic targets as well as biomarkers for diagnosis, assessment of disease activity, and response to treatment. More quantitative and sensitive high-throughput RNA profiling methods are starting to be available and will be necessary for transcriptome analyses to become routine tests in the clinical setting. We expect this to crystallize within the coming decade, as these methods become part of the personalized medicine armamentarium.

Dissecting interferon-induced transcriptional programs in human peripheral blood cells

Interferons are key modulators of the immune system, and are central to the control of many diseases. The response of immune cells to stimuli in complex populations is the product of direct and indirect effects, and of homotypic and heterotypic cell interactions. Dissecting the global transcriptional profiles of immune cell populations may provide insights into this regulatory interplay. The host transcriptional response may also be useful in discriminating between disease states, and in understanding pathophysiology. The transcriptional programs of cell populations in health therefore provide a paradigm for deconvoluting disease-associated gene expression profiles.We used human cDNA microarrays to (1) compare the gene expression programs in human peripheral blood mononuclear cells (PBMCs) elicited by 6 major mediators of the immune response: interferons alpha, beta, omega and gamma, IL12 and TNFalpha; and (2) characterize the transcriptional responses of purified immune cell populations (CD4+ and CD8+ T cells, B cells, NK cells and monocytes) to IFNgamma stimulation. We defined a highly stereotyped response to type I interferons, while responses to IFNgamma and IL12 were largely restricted to a subset of type I interferon-inducible genes. TNFalpha stimulation resulted in a distinct pattern of gene expression. Cell type-specific transcriptional programs were identified, highlighting the pronounced response of monocytes to IFNgamma, and emergent properties associated with IFN-mediated activation of mixed cell populations. This information provides a detailed view of cellular activation by immune mediators, and contributes an interpretive framework for the definition of host immune responses in a variety of disease settings.

Shifting the paradigm: host gene signatures for diagnosis of infectious diseases

Instead of focusing on the pathogen, in a paradigm shift, Zaas et al. (2009) identified host gene profiles as a strategy for diagnosis of respiratory infections. Application of host gene profiles offers tremendous possibilities for identification of diagnostic signatures, markers of disease severity, and eventually, prognostic indicators in the clinical setting.