Single-cell mapping of peripheral blood mononuclear cells reveals key transcriptomic changes favoring coronary artery lesion in IVIG-resistant Kawasaki disease

Background: Intravenous immunoglobulin (IVIG)-resistant Kawasaki disease (KD) poses a considerable challenge to patients and their families due to its severe complications. Previous researches have highlighted the critical role of immune disorders in its pathogenesis. However, fragmented studies based on isolated cases hinder a comprehensive understanding of this deadly illness. This study aimed to explore the overall landscape of peripheral blood mononuclear cells (PBMCs) in IVIG-resistant KD patients using single-cell RNA sequencing (scRNA-seq). Methods: The scRNA-seq was used to characterize the transcriptomic profiles of IVIG-resistant KD patients, IVIG-responsive KD patients, and healthy controls. Data quality control (QC) and subsequent analysis were conducted using various R packages. These included DoubletFinder and Harmony for QC, Seurat and SingleR for identifying and annotating major cell types, ggpubr for calculating and visualizing the percentages of each cell type, Seurat for characterizing differentially expressed genes (DEGs) between groups, pheatmap for visualizing the DEGs, clusterProfiler for performing Gene Ontology (GO) enrichment analysis of DEGs, scRepertoire for TCR and BCR data analysis, Monocle for assessing cell differentiation trajectories, and CellChat for intercellular interaction evaluation. Results: High-quality single-cell transcriptome data from 12 participants were analyzed, including five with IVIG-resistant KD, four with IVIG-responsive KD, and three healthy controls. We identified 10 major cell types and observed that the differentiation of CD8+ effector T cells was impeded in IVIG-resistant KD patients with coronary artery lesion (CAL) according to cell differentiation trajectory analysis. Subsequent cell communication analysis demonstrated that myeloid cluster with high expression of LCN2, S100P, and LTF played a key role, potentially signaling through MIF-CD74/CXCR4 and MIF-CD74/CD44 ligand-receptor pairs. Conclusion: Complex immunopathological changes occur during the development of CAL in IVIGresistant KD. Stunted differentiation of CD8+ effector T cells is noted in KD-CAL. Interactions between myeloid cells and T cells activates multiple inflammatory signaling pathways, with ligand-receptor pairs, including MIF-CD74/CXCR4 and MIF-CD74/CD44, potentially playing crucial roles.