Question: Cluster sub-set analysis using Seurat
1
gravatar for singcell
4 weeks ago by
singcell0
singcell0 wrote:

I am trying to dig deeper into my Seurat single-cell data analysis. If I want to further sub-cluster a big cluster then what would be the best way to do it:

1) Decreasing the resolution at FindClusters stage

or

2) extracting the individual cell index and re-clustering and then further analysis.

What would be the major difference in outcome between two approaches? How to perfrom the second approach in Seurat? using Whichcells function to extract and then I am not sure about next step please

Thank you.

seurat single-cell • 159 views
ADD COMMENTlink modified 4 weeks ago by igor11k • written 4 weeks ago by singcell0

It doesn't matter, since the raw count table are the same either way. The real important thing is the algo you use for sub-clustering. Modularity based clustering may not work good. Maybe gaussian mixture will be better?

ADD REPLYlink written 4 weeks ago by shoujun.gu370

the raw count table are the same either way

Clustering is not based on the raw counts table (in the default Seurat workflow, at least). It's based on the reduced dimensions matrix which will change when you remove some of the sources of variation.

ADD REPLYlink modified 4 weeks ago • written 4 weeks ago by igor11k
2
gravatar for swbarnes2
4 weeks ago by
swbarnes28.2k
United States
swbarnes28.2k wrote:

Subsetting and rerunning the umap coords and clustering will separate out different things better.

Personally, I cheat, and use Loupe. It's easy to use the mouse to draw categories, then export your categories into the metadata of your seurat object. But you have to have your data go through the 10xGenomics pipeline for that.

ADD COMMENTlink written 4 weeks ago by swbarnes28.2k

Never perform the analysis geometrically.

ADD REPLYlink written 4 weeks ago by shoujun.gu370

What is a geometric analysis?

ADD REPLYlink written 4 weeks ago by igor11k

Like "use the mouse to draw categories", or make conclusions based on your eyes or based on dimension reduced visualizations, etc. In my opinion, all analysis should be done analytically (mathmatically).

ADD REPLYlink written 4 weeks ago by shoujun.gu370

For scRNA-seq, there has to be substantial manual supervision to guide the analysis and make biologically meaningful conclusions. An exploratory analysis just to help you understand your data does not need to be perfect.

ADD REPLYlink modified 4 weeks ago • written 4 weeks ago by igor11k

In general, I don't quite agree with "manual supervision". All judgement should be based on some mathematical algo. If the results are not "biologically meaningful", it means that algo used is problematic here. Theoretically, researcher should improve or replace algo instead of applying "manual supervision" on the same algo.

ADD REPLYlink modified 29 days ago • written 4 weeks ago by shoujun.gu370
2
gravatar for igor
4 weeks ago by
igor11k
United States
igor11k wrote:

Both approaches are valid. It really depends on your question. You can increase the resolution to get more clusters. You can just select the cluster of interest and analyze it more closely or in different ways. For example, if your population of interest is only 10% of the total cells, it would be difficult to analyze it in the full dataset. If you just keep those cells, that becomes easier.

If you want to extract the cells of a specific cluster, you can use subset(). The Interaction Tips vignette has a few examples such as:

# Can I create a Seurat object of just the NK cells and B cells?
subset(pbmc, idents = c("NK", "B"))
ADD COMMENTlink written 4 weeks ago by igor11k

Thank you for your reply. I have got integrated Seurat object containing KO and WT samples. Can I still use subset() function? Can I still differentiate between WT and KO samples during downstream analysis after subsetting.

ADD REPLYlink written 4 weeks ago by singcell0
1

Subsetting just keeps/removes specific cells. It does not alter their metadata.

ADD REPLYlink written 4 weeks ago by igor11k
Please log in to add an answer.

Help
Access

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