FindConservedMarkers vs FindMarkers vs FindAllMarkers Seurat

Question

[Payal]

Can someone please explain the difference between FindConservedMarkers and FindMarkers.

As per definition:

FindConservedMarkers- Finds markers that are conserved between the groups

1. But does that mean that the genes are similarly expressed between groups/conditions or genes are differentially expressed between groups/conditions?
2. Also is this supposed to be for all clusters or a single cluster? If all clusters, according to the function we are identifying them for one cluster (ident.1) or max two clusters (ident.2) at a time right? Then, what ? Do I look for similar genes for all the clusters?

FindConservedMarkers(object, ident.1, ident.2 = NULL, grouping.var,
    assay = "RNA", slot = "data", meta.method = minimump,
    verbose = TRUE, ...)

FindMarkers - Finds markers (differentially expressed genes) for identity classes

I don't understand what's the significance of this?

FindAllMarkers - Finds markers (differentially expressed genes) for each of the identity classes in a dataset

What is this doing and how is this different from FindConserved Markers?

It would be very helpful if someone can explain this with a biological use cases.

Thanks in advance!
Payal

Answer 1

[jared.andrews07]

FindMarkers will find markers between two different identity groups - you have to specify both identity groups. This is useful for comparing the differences between two specific groups.

FindAllMarkers will find markers differentially expressed in each identity group by comparing it to all of the others - you don't have to manually define anything. Note that markers may bleed over between closely-related groups - they are not forced to be specific to only one group. This is what most people use (and likely what you want).

FindConservedMarkers will find markers that are conserved between two groups - this can be useful if you want to find markers that are conserved between a treated and untreated condition for a specific cell type or group of cells. It means they are differentially expressed compared to other groups, but have similar expression between the two groups you're actually comparing.

[Payal]

Thank you! This is helpful.

I have a combined dataset (control vs treatment ). I was following this - https://satijalab.org/seurat/v3.1/immune_alignment.html tutorial. I got confused after the FindConservedMarkers step.

In general on what basis should we select the marker genes to draw the feature plot or identify the clusters? In my case, should I be looking at the genes with low p-value and high log fold change difference between control and treatment for say cluster 1?

[jared.andrews07]

This depends on your experimental setup, research question, and cells. It's difficult to answer without more information. If it seems like cell types are clustering together, despite the treatment, then I think your approach would be valid. If not, you may want to assign cell types to your cells and then compare between treatment and control for each given cell type.

[m.sadman.sakib]

Hi Jared! Thanks a lot for this explanation! I am trying to identify differential expression in each cell cluster for control vs treatment. For this, which function will be appropriate? Thanks a lot in advance!

[jared.andrews07]

See their DE vignette. Generally, you'll want to use FindMarkers, though in your case you'll probably want to loop through each cluster and subset to your cell of interest before comparing the conditions. Or combine the cluster/treatment info into a single column in the metadata and use that to specify the two groups.

Answer 2

[Payal]

I found this link. Its very well explained there too - https://hbctraining.github.io/scRNA-seq/lessons/sc_exercises_integ_marker_identification.html