“…To identify cell types and states at chosen clustering resolution (0.6 for Control only, 0.5 for Control and Infected integrated), we performed differential expression analysis using Wilcoxon rank sum-based method (log2 fold change > 0.25, % of expressing cells > 0.25). Known markers were used to identify T cells (cd247l, lck, il7r, and cxcr4a) (53)(54)(55), dendritic cell-like populations (ctsbb, tlr7) (56,57), B cells (cd37, pax5) (58,59), macrophages (mpeg1.1, grn1) (60,61), neutrophils (mpx, il6r) (11,62), erythrocytes (hbba2, hemgn) (63), thrombocytes (thbs1b, fn1b) (57), superficial epithelial cells (krt1-19d, cldne) (64), ionocytes (trpv6, foxi3b) (65), intermediate epithelial cells (cldna, tp63) (64), basal epithelial cells (cldn1, cldni) (64), mesenchymal cells (vcana, clu) (66), lateral line-like cells (prox1a, prr15la) (67), and epidermal mucous cells (agr2, cldnh) (64). We found one population with some cells locating closer to the T cell populations but without clear enrichment for known signatures and included this population for a subclustering analysis with other T cells.…”