When the exposed group was analyzed in the time course (starting at 4dpe), a statistically significant increasing trend (before Bonferroni correction) was found for IgM+ B cells over the course of the infection (F([1],[36]) = 4.89, p = .033, adjusted R2 = .10, coeff = 0.20%, SE = 0.09%) and a statistically significant decreasing trend (before Bonferroni correction) for myeloid cells (F([1],[43]) = 4.66, p = .036, adjusted R2 = .08, coeff = −0.22%, SE = 0.10%), while no statistically significant trend was observed for CD8+ T cells (F([1],[50]) = 1.74, p = .193) (Fig 3D). The gene discussed is CD8A; the disease is infection.