Question: scRNAseq Seurat conserved markers among clusters
gravatar for michelafrascoli
21 months ago by
michelafrascoli10 wrote:

Hello, I am trying to find the most highly expressed genes that are conserved in all my scRNAseq clusters with Seurat.

Here is what I tried to do:

YSa.markers <- FindMarkers(object = YSa, ident.1 = 0, ident.2 = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14), min.pct = 0.8)
print (x = head(x = YSa.markers, n = 10))

Where 0 to 14 are my clusters.

But when I make Violin plots of the genes that it gives me it’s obvious they are not expressed in all the clusters. I am clearly doing something wrong and I was wondering if anyone could help me with this.

Thank you very much for your time, Michela

rna-seq • 1.3k views
ADD COMMENTlink modified 20 months ago by geek_y11k • written 21 months ago by michelafrascoli10
gravatar for geek_y
20 months ago by
geek_y11k wrote:

Here you are trying to find markers per cluster but not genes expressed across all clusters.

You can get a list of genes that do not vary across clusters by using rowVars() function on normalised expression matrix. If the matrix is too sparse, you can average the data across similar cells by using KNN or just by pooling counts from each cluster.

ADD COMMENTlink written 20 months ago by geek_y11k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1827 users visited in the last hour