plotting two genes in Seurat ScRNAseq

Question

[rob.costa1234]

I have a clarification about about using Seurat package (link is https://satijalab.org/seurat/pbmc3k_tutorial.html)-

FeaturePlot(object = pbmc, features.plot = c("MS4A1", "GNLY", "CD3E", "CD14", 
    "FCER1A", "FCGR3A", "LYZ", "PPBP", "CD8A"), cols.use = c("grey", "blue"), 
    reduction.use = "tsne")

At this time only one gene is plotted on the cluster. Can I plot two genes simultaneously in different colors on the clusters? if so how can I do that.

Thanks

[dppb05]

Can I plot two genes simultaneously in different colors on the clusters?

Let me see if I understood correctly: you wish to plot the expression of several genes, each with a different colour, in a single tSNE plot?

[rob.costa1234]

Yes that is exactly what I want to do

Answer 1

[dppb05]

Can I plot two genes simultaneously in different colors on the clusters?

If it is just two genes, the overlay option in Seurat's FeaturePlot can help you with that.

FeaturePlot(object = pbmc,
    features.plot = c("LYZ", "CD8A"),
    cols.use = c("grey", "blue", "red", "pink"), 
    reduction.use = "tsne",
    overlay = TRUE)

The cols.use parameter indicate, respectively, low expression of both genes, high expression of LYZ, high epxression of CD8A, and high expression of both genes. Of course, you can change the colours and genes as you prefer, but keep in mind that you can only use two genes and four colours. If you pass less than four colours the other colours are going to be defined automatically by Seurat.

[Za]

Sorry, how I can change the color of some cells in feature plot? I have a list of cells, I would like to map them in feature plot

[dppb05]

@Za There are a few different ways you can do that. I suggest you create a new post with your question detailing a bit how your list of cells is stored.

[biologo]

Hi! Za i wanna know do you figure how to plot several genes simultaneously in Featureplot? because i am also find how to do with it.

[Za]

For mapping several cells this is the solution

https://bioinformatics.stackexchange.com/questions/4470/mapping-a-list-of-cells-in-seurat-featureplot

for mapping several genes if I am not wring by a vector of those genes you can do that by seurat

FeaturePlot(object = seurat,features.plot = c("list of your genes"))