When calling a gene and something from the metadata along with a gene, the plots remain in alphabetical order (I think), but the names are in the order listed in the code. library(SeuratData) library(Seurat) InstallData("pbmc3k") data('pbmc3k') pbmc3k <- UpdateSeuratObject(pbmc3k) ...