candidate_markers <- c("GAPDH","NANOG","ACTB")

for(this_gene in candidate_markers){

   ## Violin Plot
  vln_plot_avgdiff <- VlnPlot(expression_seurat.updated,
            this_gene,
            point.size.use = -1,
            x.lab.rot=TRUE,
            remove.legend = TRUE,
            size.title.use = 0 ,
            y.max = 6,
            size.y.use = 0,
            do.return = TRUE
            )
  
  vln_plot_avgdiff <- vln_plot_avgdiff +
    theme(
  axis.text.x = element_blank(),
  axis.ticks = element_blank(),
  axis.line.x = element_blank(),
  axis.line.y = element_blank()) +
    coord_flip()
  
  vln_plot_avgdiff