GSVA was performed according to the published protocol (https://bioconductor.org/packages/devel/bioc/vignettes/GSVA/inst/doc/GSVA.html).
library("Seurat")
library("GSVA")
# load your seurat object and the signature genes you are interested in, here: "your.seurat.object" and "signature.genes.1"
signature.genes.1<-c("Adam19", "Amica1", "Ap1s3", "Ass1", "Bcl11a", "Btla", "Ccr7", "Flt3", "Gpr114",
"Gpr132", "Gpr68", "Gpr82", "H2-Eb2", "Hmgn3", "Kit", "Klri1", "Kmo", "P2ry10", "Pvrl1",
"Rab30", "Septin6", "Slamf7", "Traf1", "Zbtb46")
gs <- list(sig1=signature.genes.1)
Y<-GetAssayData(object=your.seurat.object, assay = "RNA", slot="data")
gsva.es <- gsva(Y, gs, verbose=FALSE)
your.seurat.object[["sig1"]]<-gsva.es["sig1",]
FeaturePlot(your.seurat.object, "sig1", cols=c("blue", "lightgrey", "red") )