313 cells in D54 SST+ clusters were selected for differential expression analysis. The R package “limma” was used to test each gene for variable expression between clusters, and 2304 genes were statistically significant (Benjamini-Hochberg corrected P < 0.05). For the top 50 most significantly variable genes, the expression z-score across the four clusters was calculated and plotted as a heatmap using the R package “pheatmap”.