Molecular and temporal profile of monocyte/macrophage subtype heterogeneity acutely after SCI. (A) GO enrichment analysis of the DEGs distinguishing the two macrophage subtypes. GO biological process terms displayed along y axis and –log10 (Bonferroni-corrected P values) displayed along x axis. (B) Proportion of each peripheral myeloid subtype among all peripheral myeloid cells at each time point. (C) Dot plot of marker genes that differentiate peripheral myeloid subtypes. Color of dots represents z-scored expression level, and size of dots represents percentage of cells with at least one UMI detected per gene. (D) Heatmap of the highest DEGs per peripheral myeloid subtype. Color represents z-score expression level. (E) Heatmaps of M1 marker gene (left) and M2 marker gene (right) expression among all myeloid subtypes show that neither inflammatory nor chemotaxis-inducing subtypes exclusively display M1 or M2 gene signatures. A notable exception is higher expression of M2 markers Chil3 and Arg1 in monocytes and chemotaxis-inducing macrophages compared with inflammatory macrophages. Dendrograms along axes display hierarchical clustering results. Color of tiles represents z-scored expression value. Inset stars in C–E indicate statistically significant greater expression compared with all other cells combined (adjusted P value < 10−10; Wilcoxon rank-sum with Bonferroni correction; see Materials and methods). Prop., proportion; Mac, macrophage.