Differential gene expression at 20 DPI. (A) Venn diagram shows the number of overlapping and unique DE genes for macrophages, CD4 T cells, CD8 T cells, and B cells at 20 DPI. Plasma cells were included in the B cells. (B) Bubble plots presenting all enriched GO terms shared across B cell, CD4 T cell, and CD8 T cell populations at 20 DPI, based on all genes with an adjusted P-value of less than 0.05. GO terms annotated with fewer than 3 DE genes were omitted. Bubble size represents the ratio of the number of DE genes associated with a specific GO term to the total number of DE genes identified in that cluster. Bubble color represents the adjusted P-value. Terms are ranked according to gene ratio values within the 3 cell populations. Macrophage populations did not show significantly enriched GO terms. Clusters with no shared enriched GO terms are not shown. (C) Bubble plots presenting immune-related enriched KEGG pathways appeared across B cell, CD4 T cell, CD8 T cell, and macrophage populations at 20 DPI, based on all genes with an adjusted P-value of less than 0.05. Bubble size represents the ratio of the number of DE genes associated with a specific KEGG pathways to the total number of DE genes identified in that cluster. Significantly enriched KEGG pathways were manually filtered to concentrate on immune-related pathways. Bubble color represents the adjusted P-value. Terms are ranked according to gene ratio values within the 4 cell populations. Clusters with no shared enriched KEGG pathways are not shown. Macrophages (M).