Mouse skin wound infection atlas.

(A) Gene expression of healing markers Pdgfa, Tgfb1 and Egf four days post-infection (4 dpi) for uninfected and E. faecalis-infected 6–7-week-old C57BL/6J mouse skin wounds, normalized to intact skin (n = 6 for skin; n = 8 for wounds; one-way ANOVA). (B) Single-cell RNA sequencing workflow of full-thickness mouse wounds. (C) Integrated dataset of ∼23,000 scRNA-seq libraries from uninfected and infected wounds identifies 5 mega cell classes indicated in the UMAP. (D) UMAP colored by the uninfected and infected conditions. (E) Schematic describing the color-matched cell types with that of the clusters in Figure S1F. (F) Dot plot depicting the top two cell type-specific markers in the integrated data. Legend indicates average expression and dot size represents percent expression. (G) Density plots depict cell types described in C. (H) Heat map of weighted gene co-expression network analysis for annotated cell populations, colored by bars matching mega clusters and annotated cell types in C.

Sub-clustering of keratinocyte populations reveals infection-specific cell types.

(A) UMAP of integrated keratinocyte (basal, suprabasal, hair follicle, bulge, and sebaceous gland) population reveals 19 clusters. (B) Infected keratinocytes (green) show unique and shared clusters with uninfected keratinocytes (purple). (C) Spatial dispersion of Krt5 and Krt10 abundance in keratinocytes. (D) Expression of Lgr5, Ivl, Lor, and Krt15 in Louvain clusters shown in A. (E) Heat map of top 5 differentially expressed marker genes for each cluster in keratinocytes. Rectangle boxes indicate infection-specific Louvain clusters. (F-G) The bar plots show the top 15 Gene Ontology terms for infection-specific (F) Zeb2hi (cluster 0) and (G) Gjb6hi (cluster 6) keratinocyte populations. (H-I) Dynamic RNA velocity estimation of uninfected (H) and infected (I) keratinocytes. (J) The top lineage driver gene, Rgs1, in infected keratinocytes, was ubiquitously expressed in infection-specific Louvain clusters. (K) Gene expression dynamics resolved along latent time in the top 50 likelihood-ranked genes of infected keratinocytes. The colored bar at the top indicates Louvain clusters in I. Legend describes scaled gene expression. (L) Inferred ligand-receptor pairs outgoing from infected keratinocytes. (M) The hierarchy tree depicts the trajectory of differentiating and terminally differentiated keratinocyte cells originating from basal keratinocytes.

E. faecalis delays immune response in fibroblasts.

(A) UMAP of integrated fibroblasts reveals 12 clusters. (B) Infected fibroblasts (green) show unique and shared clusters with uninfected fibroblasts (purple). (C) Spatial dispersion of Col1a1 and Col1a2 abundance in fibroblasts. (D) Expression of Eln, Aspn, Pdgfrb and Fap in Louvain clusters shown in A. (E) Heat map of top 5 differentially expressed marker genes for each cluster in fibroblasts. Rectangle boxes indicate infection-specific Louvain clusters. (F-G) The bar plots show the top 15 Gene Ontology terms for infection-specific (F) Lyz2hi/Taglnhi (cluster 0) and (G) Timp1hi (cluster 2) fibroblast populations. (H-I) Dynamic RNA velocity estimation of uninfected (H) and infected (I) fibroblasts. (J) The top lineage driver gene, Csgalnac1, expression in infected fibroblasts. (K) Gene expression dynamics resolved along latent time in the top 50 likelihood-ranked genes of infected fibroblasts. The colored bar at the top indicates Louvain clusters in I. Legend describes scaled gene expression. (L) While uninfected fibroblasts show healing phenotypes, infected fibroblasts undergo two transitioning phases: (i) contractile and (ii) pathologic.

Macrophages display M2-like polarization.

(A) UMAP of the integrated myeloid population reveals 9 macrophages, 2 dendritic cells, and one Langerhans cell cluster. (B) Infected myeloid (green) population show unique and shared clusters with the uninfected myeloid (purple) population. (C) Spatial distribution of Lyz2 (macrophage) and Itgax (DC and LC) abundance in fibroblasts. (D) Gene expression Mrc1, Nos2 and Arg1 in mouse wounds at 4dpi, normalized to homeostatic skin (n = 6 for skin; n = 8 for wounds; one-way ANOVA). (E) In vitro infection of unpolarized (M0) bone marrow-derived murine macrophages (BMDMs) resulted in a down-regulation of TRM-associated marker Mrc1 and an upregulation of M2-like markers Nos2 and Arg1. Data are pooled from 2 biological replicates (shown in light and dark circles respectively) of 3 technical replicates each. (F) Spatial distribution of TRM and M2-like macrophage markers, Mrc1 and Arg1, respectively. (G) Expression of Mrc1, Adgre1, Arg1, Itgax, Cd207 and Nos2 in the integrated myeloid dataset. (H-I) Dynamic RNA velocity estimation of uninfected (H) and infected (I) myeloid cells. (J) The top lineage driver gene, Fth1, was ubiquitously expressed in terminal macrophage populations. (K) Putative driver genes of infected macrophages. (L) The proposed model describes macrophage characteristics, where neutrophil-attracting and wound repair-associated macrophages were involved in uninfected wound healing. In contrast, bacteria-infected wounds are enriched in efferocytotic macrophages and matrix-producing macrophages.

Crosstalk between neutrophils and anti-inflammatory macrophages regulates the CCL signaling pathway.

(A) UMAP of the integrated myeloid population reveals 6 Louvain clusters. (B) The infected neutrophil (green) population shows unique and shared clusters with the uninfected neutrophil (purple) population. (B) Spatial organization of Csf3r and Itgam abundance in neutrophils. (D) Expression of Cxcr2, Fcgr3, Fth1 and Camk1d in Louvain clusters. (E) Heat map of top 5 differentially expressed marker genes in neutrophils. Rectangle boxes indicate infection-specific Louvain clusters. (F-G) The bar plots show the top 15 Gene Ontology terms for infection-specific (F) Lrg1hi (cluster 0) and (G) Csf1hi (cluster 2) populations. (H) Dynamic RNA velocity estimation of infected neutrophils. (I) The top lineage driver gene, Entpd1, was ubiquitously expressed in Lrg1hi neutrophils (cluster 0). (J) Putative driver genes of infected neutrophil clusters. (K) Cytokine signaling pathway (CCL) cell-cell interaction map. (L) Gene expression of cytokines Ccl3, Ccl4, and Ccl6 in neutrophils. (M) Ccl3:Ccr1 ligand-receptor interaction mediate neutrophil-macrophage crosstalk.

Macrophage-EC interactions display an anti-inflammatory niche.

(A) UMAP of integrated endothelial cells reveals 13 clusters. (B) Infected endothelial cells (green) show unique and shared clusters with uninfected endothelial cells (purple). (C-D) Gene Ontology analysis of infection-specific clusters 0 (C) and 8 (D). (E) Dynamic RNA velocity estimation of infected endothelial cells. (F) The top lineage driver genes, Malat1, Tcf4, Plcb1, Diaph2, Bmpr2, and Adamts9 in infected endothelial cells. (G) NicheNet interaction heat map between keratinocytes, fibroblasts, macrophages, and endothelial cells. Note the macrophage Il1b-specific Bgn induction in infected ECs. (H) The dot plot depicts interactions between endothelial cells (receptors) and keratinocytes, fibroblasts, and macrophages (ligands). Rows demonstrate a ligand-receptor pair for the indicated cell-cell interactions (column). (I) Differential interaction strengths of a cellular interactome for TNF, SPP1 and CXCL signaling pathways between all cell types in infection. (J) Spp1, Cxcl2, and Cxcl3 gene expression in keratinocytes, M2-like macrophages, and fibroblasts (Wilcoxon Rank Sum test, ****p < 0.0001).

Key resources table