The 10× single-cell RNA-seq was used in the experiments. (A) The results were analysed using Loupe Cell Browser with t-SNE method. The percentages of STAT3⁺ cells in NC and RRV groups is presented. (B) The two subclusters of Gr-1⁺ cells was showed using Ly6G and Ly6C as markers and the percentage was indicated. (C) The expression of STAT3 in each subclusters and the percentage was indicated. (D) Five up- and down-regulated genes each in BA compared with the control group. *P<0.1, **P<0.05, ***P<0.01 and ****P<0.001, the P-values adjusted using the Benjamini–Hochberg correction for multiple testes was showed in the graph. (E) Five up- and down-regulated genes expressed in the sub clusters compared BA and control group with Ly6G⁺ and Ly6C⁺indicted. The significant P-value label was the same as previous. (F) Metascape analysis showing the enrichment of gene expression in different signalling pathway and interaction of molecules in the context of genes that are highly differentially expressed between NC and RRV groups. (G) The analysis of the most highly expressed genes and comparison of the each subcluster to determine functional differences.