How to get top 10 genes expressed in a cluster ( for cell type identification)?
1
0
Entering edit mode
15 months ago
Chris ▴ 340

Hi all,

I try to label clusters but singleR result is not consistent with the marker gene for a cell type. Findmarker() or Findallmarker() seems could not get genes express in a cluster but they are used for comparison. Would you please suggest the function or code I should use in this case? I appreciate it!

seurat • 1.6k views
ADD COMMENT
0
Entering edit mode

These two functions- Findmarker() or Findallmarker() had been discussed before. Please check it out in the link below:-

FindConservedMarkers vs FindMarkers vs FindAllMarkers Seurat

ADD REPLY
0
Entering edit mode

Thank you for the information. I have read it before, very helpful. Unfortunately, I can't use it for this question.

ADD REPLY
0
Entering edit mode
15 months ago
yhdist ▴ 70

You want to find the genes with the highest expression or the highest rates of differential expression? Surely you can find the 10 "most epxressed" genes on your vst data before drawing your Heatmap, if you have any.

There can be a million reasons why your marker does not reflect your expectations depending on your process.

ADD COMMENT
0
Entering edit mode

Just a modification that I want to find all genes expressed in a cluster. I think Findallmarker() will do that.

ADD REPLY
0
Entering edit mode

The Seurat object already tells you what the clusters are, if that's what you want.

ADD REPLY
0
Entering edit mode

Thank you. Seurat object will tell the cluster number but not which genes or cell type of that cluster. I think I figured it out.

ADD REPLY
0
Entering edit mode

It will tell you which cells are involved, but it won't be able to tell you the cell type unless you've already labeled them. I hope this worked out in the end.

ADD REPLY

Login before adding your answer.

Traffic: 1553 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6