Question: Skipping Clustering in Seurat
gravatar for suvratha
8 months ago by
suvratha60 wrote:


I was wondering if there is a possibility to skip the clustering step in Seurat and just do the differential gene expression as I know what are the cell populations I've used as my input data. I don't want Seurat to find any novel subpopulations of the cells but I just want it to find the differential gene expression between 2 different populations of the cells in my input data.

Thanks Suvi

ADD COMMENTlink modified 8 months ago by jared.andrews077.9k • written 8 months ago by suvratha60

I think you can perform differential expression analysis without performing the initial PCA, but you have group cells based on some gene expression pattern.

ADD REPLYlink modified 8 months ago • written 8 months ago by Arup Ghosh2.7k
gravatar for jared.andrews07
8 months ago by
Memphis, TN
jared.andrews077.9k wrote:

Sure, if you can assign whatever labels necessary to each cell (by barcode) to group them as necessary in the metadata of the Seurat object. Then you can just run FindMarkers between the two groups after setting the Ident of each cell to that metadata variable.

So in short, the commands you want are AddMetaData, StashIdent (to save the original barcodes to a new metadata column for posterity), Idents (to set the cell identities to whatever group variable you added with AddMetaData), and FindMarkers (to actually perform the differential expression between the groups).

ADD COMMENTlink written 8 months ago by jared.andrews077.9k

my input is a read count matrix into Seurat. I'm making this into a Seurat object initially, so to this I do AddMetaData to say which cell type this is, then add Idents and use them in Findmarkers. Is this correct?

ADD REPLYlink written 8 months ago by suvratha60

Yes. Make your Seurat object, then perform those steps.

ADD REPLYlink written 8 months ago by jared.andrews077.9k
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: 2119 users visited in the last hour