Identification of key signaling pathways and hub genes related to immune infiltration in Kawasaki disease with resistance to intravenous immunoglobulin based on weighted gene co-expression network analysis

Background: Kawasaki disease (KD) is an acute vasculitis, that is, the leading cause of acquired heart disease in children, with approximately 10%-20% of patients with KD suffering intravenous immunoglobulin (IVIG) resistance. Although the underlying mechanism of this phenomenon remains unclear, recent studies have revealed that immune cell infiltration may associate with its occurrence. Methods: In this study, we downloaded the expression profiles from the GSE48498 and GSE16797 datasets in the Gene Expression Omnibus database, analyzed differentially expressed genes (DEGs), and intersected the DEGs with the immune-related genes downloaded from the ImmPort database to obtain differentially expressed immune-related genes (DEIGs). Then CIBERSORT algorithm was used to calculate the immune cell compositions, followed by the WGCNA analysis to identify the module genes associated with immune cell infiltration. Next, we took the intersection of the selected module genes and DEIGs, then performed GO and KEGG enrichment analysis. Moreover, ROC curve validation, Spearman analysis with immune cells, TF, and miRNA regulation network, and potential drug prediction were implemented for the finally obtained hub genes. Results: The CIBERSORT algorithm showed that neutrophil expression was significantly higher in IVIG-resistant patients compared to IVIG-responsive patients. Next, we got differentially expressed neutrophil-related genes by intersecting DEIGs with neutrophil-related module genes obtained by WGCNA, for further analysis. Enrichment analysis revealed that these genes were associated with immune pathways, such as cytokine-cytokine receptor interaction and neutrophil extracellular trap formation. Then we combined the PPI network in the STRING database with the MCODE plugin in Cytoscape and identified 6 hub genes (TLR8, AQP9, CXCR1, FPR2, HCK, and IL1R2), which had good diagnostic performance in IVIG resistance according to ROC analysis. Furthermore, Spearman's correlation analysis confirmed that these genes were closely related to neutrophils. Finally, TFs, miRNAs, and potential drugs targeting the hub genes were predicted, and TF-, miRNA-, and drug-gene networks were constructed. Conclusion: This study found that the 6 hub genes (TLR8, AQP9, CXCR1, FPR2, HCK, and IL1R2) were significantly associated with neutrophil cell infiltration, which played an important role in IVIG resistance. In a word, this work rendered potential diagnostic biomarkers and prospective therapeutic targets for IVIG-resistant patients.

Exploring the relationship between pyroptosis, infiltrating immune cells and Kawasaki disease with resistance to intravenous immunoglobulin (IVIG) via bioinformatic analysis

BACKGROUND

Kawasaki disease (KD) is a kind of vasculitis predominantly afflicting children younger than five. Although intravenous immunoglobulin (IVIG) has been regarded as the first-line therapy, there are some children unresponsive to it, resulting in higher risk of coronary artery aneurysms (CAA), the most severe complication of KD. Pyroptosis is an inflammatory apoptosis, which resembles the traits of IVIG-resistance. Therefore, our research aims to find relationships between KD with IVIG-resistance and pyroptosis, and provide the underlying mechanisms of IVIG-resistance.

METHODS

The transcriptome data of three datasets were downloaded from Gene Expression Omnibus (GEO) database. CIBERSORTx and WGCNA were combined to identify the coexpression gene network correlated with the up-regulated immune cells in KD, using differentially expressed genes (DEGs) overlapped in GSE68004 and GSE73461. The key genes in hub module were intersected with pyroptosis-related genes (PRGs). Then KD patients were divided into subgroups according to the expression of remaining genes, along with the construction of risk score (RS) based on the least absolute shrinkage and selection operator (LASSO) regression analysis. Besides, we explored the clinical value of RS between IVIG-responsive and -resistant KD patients in GSE16797. In addition, the biological pathways between subgroups were evaluated using Gene Set Variation Analysis (GSVA).

RESULTS

A total of 4246 DEGs and three immune cells, including Monocytes, M0 macrophage, and neutrophils, were analyzed with P < 0.05 between KD and healthy controls (HCs). The lightcyan module was the hub module based on WGCNA, and only NLRC4, CASP1, CASP4, GSDMD, IL1B and PYCARD in the hub module were overlapped with PRGs. Then KD patients in GSE68004 were stratified into two clusters on the basis of the expression levels of six genes. RS was built with five out of six genes (exclude PYCARD) according to the LASSO analysis, which could differentiate C1 from C2, IVIG-responsive from -resistant KD patients. Besides, the high-risk group (C1) tended to be with increased levels of inflammation, immune responses and infiltration of neutrophils according to the analysis of GSVA and CIBERSORTx.

CONCLUSION

We built a pyroptosis-related RS to evaluate the degree of pyroptosis and infiltrating immune cells in subgroups of KD, and associated it with the responsiveness to IVIG, which might help us to further understand the pathological process during IVIG-nonresponse.

A low-frequency IL4R locus variant in Japanese patients with intravenous immunoglobulin therapy-unresponsive Kawasaki disease

BACKGROUND

Kawasaki disease (KD) is a systemic vasculitis which may be associated with coronary artery aneurysms. A notable risk factor for the development of coronary artery aneurysms is resistance to intravenous immunoglobulin (IVIG) therapy, which comprises standard treatment for the acute phase of KD. The cause of IVIG resistance in KD is largely unknown; however, the contribution of genetic factors, especially variants in immune-related genes, has been suspected.

METHODS

To explore genetic variants related to IVIG-unresponsiveness, we designated KD patients who did not respond to both first and second courses of IVIG therapy as IVIG-unresponsive patients. Using genomic DNA from 30 IVIG-unresponsive KD patients, we performed pooled genome sequencing targeting 39 immune-related cytokine receptor genes.

RESULTS

The single nucleotide variant (SNV), rs563535954 (located in the IL4R locus), was concentrated in IVIG-unresponsive KD patients. Individual genotyping showed that the minor allele of rs563535954 was present in 4/33 patients with IVIG-unresponsive KD, compared with 20/1063 individuals in the Japanese genome variation database (odds ratio = 7.19, 95% confidence interval 2.43-21.47). Furthermore, the minor allele of rs563535954 was absent in 42 KD patients who responded to IVIG treatment (P = 0.0337), indicating that a low-frequency variant, rs563535954, is associated with IVIG-unresponsiveness in KD patients. Although rs563535954 is located in the 3'-untranslated region of IL4R, there was no alternation in IL4R expression associated with the mior allele of rs563535954. However, IVIG-unresponsive patients that exhibited the minor allele of rs563535954 tended to be classified into the low-risk group (based on previously reported risk scores) for prediction of IVIG-resistance. Therefore, IVIG-unresponsiveness associated with the minor allele of rs563535954 might differ from IVIG-unresponsiveness associated with previous risk factors used to evaluate IVIG-unresponsiveness in KD.

CONCLUSION

These findings suggest that the SNV rs563535954 could serve as a predictive indicator of IVIG-unresponsiveness, thereby improving the sensitivity of risk scoring systems, and may aid in prevention of coronary artery lesions in KD patients.

The Kawasaki Disease Comparative Effectiveness (KIDCARE) trial: A phase III, randomized trial of second intravenous immunoglobulin versus infliximab for resistant Kawasaki disease

BACKGROUND

Although intravenous immunoglobulin (IVIG) is effective therapy for Kawasaki disease (KD), the most common cause of acquired heart disease in children, 10-20% of patients are IVIG-resistant and require additional therapy. This group has an increased risk of coronary artery aneurysms (CAA) and there has been no adequately powered, randomized clinical trial in a multi-ethnic population to determine the optimal therapy for IVIG-resistant patients.

OBJECTIVES

The primary outcome is duration of fever in IVIG-resistant patients randomized to treatment with either infliximab or a second IVIG infusion. Secondary outcomes include comparison of inflammatory markers, duration of hospitalization, and coronary artery outcome. An exploratory aim records parent-reported outcomes including signs, symptoms and treatment experience.

METHODS

The KIDCARE trial is a 30-site randomized Phase III comparative effectiveness trial in KD patients with fever ≥36 h after the completion of their first IVIG treatment. Eligible patients will be randomized to receive either a second dose of IVIG (2 g/kg) or infliximab (10 mg/kg). Subjects with persistent or recrudescent fever at 24 h following completion of the first study treatment will cross-over to the other treatment arm. Subjects will exit the study after their first outpatient visit (5-18 days following last study treatment). The parent-reported outcomes, collected daily during hospitalization and at home, will be compared by study arm.

CONCLUSION

This trial will contribute to the management of IVIG-resistant patients by establishing the relative efficacy of a second dose of IVIG compared to infliximab and will provide data regarding the patient/parent experience of these treatments.