This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Identifying DEG's on Integrated samples

Hello,

I am analyzing snRNA-seq data that is composed of a KO and WT sample. I have integrated/anchored and clustered the two datasets to identify similarities and differences across the two samples. I have identified the DEGs for each cluster of each sample independently using FindMarkers... But beyond this we would like to generate a list of DEGs for each cluster that compares the KO and WT in order to identify the differences our KO is causing. Is there a function to specifically look at DEGs between samples in specific clusters?

Thank you for your help.

seurat deg

1 answer

You can add a new column to your meta data and use that column to run FindMarkers again

seuratobj@meta.data$new_column = paste0(seuratobj@meta.data$sample_name,"_",seuratobj@meta.data$cluster,sep="")
Idents(seuratobj) = "new_column"
deg = FindAllMarkers(seuratobj)

Log in to answer this question.