We annotated clusters based on highly expressed genes, then sub-clusters were characterized by performing PCA dimensionality reduction and clustering over those cells alone, and annotated based on highly expressed genes found via one-versus-rest differential expression test (Wilcoxon) within the major cell type.