At last, we run the t-SNE dimensionality reduction algorithm on a combined data sample using the Seurat package based on harmony embedding (https://github.com/immunogenomics/harmony, version 1.0.0).