To classify each cell subpopulation in an unbiased manner, we separately reclustered the cells of each lineage. By analyzing the most significantly upregulated genes in each cluster in scRNA-seq analysis, we identified five distinct subsets of CD3+ TCs (Fig. S3A), five distinct subsets of CD4+ TCs (Fig. S3B), four distinct subsets of CD8+ TCs (Fig. S3C), three distinct subsets of NKs (Fig. S3D), four distinct subsets of BCs (Fig. S3E), three distinct subsets of MCs (Fig. S3F) and four distinct subsets of DCs (Fig. S3G, see Table S3A for the details).