Single cell sequencing analysis

Sequence data analysis was performed using Cell Ranger (v6.0.2)⁹ multi and aggr functions with default parameters. Downstream analysis was performed using R (R Core Team, 2020) and Seurat¹⁰ using filtered UMI count matrices and default parameters. Harmony was applied to integrate individual data sets; harmony reduction was used to generate uniform manifold approximation and projection (UMAP) embeddings and graph-based clustering¹¹.