The clustering of cells was performed using the Seurat v2.3.4 R package59. In brief, initial bulk samples containing single cell data matrices were column-normalized and log-transformed. Quality control parameters were used to filter cells with 200–5000 genes with a mitochondrial percentage ...