Fig. 5

Functional-shifting inflammatory myeloid cells are associated with complications emerge in KD. UMAP projection of immune cell profiles from KD PBMCs, (A) group by cell types (B) and predicted groups. C Boxplots showing CytoTRACE values for different predicted groups of KD PBMCs. D Alluvial diagram showed the different shifting rates from KD PBMCs mapping to TAIA. Axis 1, 2, and 3 represented real age, predicted groups, and cell types. Dotchart showed GO-BP (E) and KEGG (F) based enrichment results of shift- versus preserve- myeloid upregulated genes. Color indicates log₁₀Pvalue, and size indicates gene counts. (G) lipid-metabolic signature pathways in the main immune cell types of different age groups. Both color and size indicate the VISION score (H) UMAP projection of immune cell profiles from KD PBMCs after IVIG treatment, group by cell types. I Alluvial diagram showed the different shifting rates from KD PBMCs after IVIG treatment mapping to TAIA. Axis 1, 2, 3 represented real age, predicted groups and cell types. J Percentages of different myeloid and NK cell shift- subsets between before and after IVIG treatment, color-coded by cell subsets. K Violin plot showing shift- myeloid (upper) and NK (lower) signature score before and after IVIG treatment (GSE168732). The P values were calculated with wilcox.test. L Box plot showing shift- myeloid signature scores before and after IVIG treatment (left) or IFX treatment (right) based on ssGSEA algorithm (GSE48498). M Box plot showing shift- myeloid signature scores before and after IVIG treatment (left) or differential response groups (right) based on ssGSEA algorithm (GSE16797). N Spearman correlation analysis between shift- myeloid and T cell signature score (GSE16797). O Box plot showing shift- myeloid signature scores in different periods of onset based on ssGSEA algorithm (GSE73463). In figure L, M and O, the P values are calculated with t-tests (two-sided). The P values in Figure C is calculated with kruskal.test