Downsampling one of the sample on the UMAP clustering to match the number of cells of the other group

- Individual seurat object for WT and KO - Merged WT_KO - Filtered_WT_KO - Split_cond_WT_KO (Normalize, cell cycle scoring and SCTransform) - SelectIntegrationFeatures - PrepSCTIntegration - FindIntegrationAnchors and SCT Normalization - IntegrateData - Run PCA - FindNeighbors - FindClusters - Run UMAP

10 months ago by

jared.andrews07 ♦ 8.3k

Memphis, TN

jared.andrews07 ♦ 8.3k wrote:

Approach 1 is a poor idea. Downsampling removes information, lowering the power of your differential expression analysis. This could result in marker genes being lost. If they just want to visualize them equally, I would just grab equal numbers of cells from each condition for visualization and feed them to the cells parameter for FeaturePlot or cells.use parameter for DimPlot.

Alternatively, you could plot differences using methods that show summary statistics for each group, like boxplots, violin plots, or ridgeplots. While the Seurat functions do okay with these, I prefer using dittoSeq, which allows for much greater customization and generally just looks better by default.

You don't need to re-run your entire analysis. Approach 2 is fine.

ADD COMMENT • link written 10 months ago by jared.andrews07 ♦ 8.3k

Yeah I tried the approach2.

WT_cells <-  which(WT_KO_integrated_seurat$sample == 'WT')
KO_cells <- which(WT_KO_integrated_seurat$sample == 'KO')
downsampled_KO_cells <- sample(KO_cells, 895)
WT_KO_integrated_downsampled <- WT_KO_integrated_seurat[,c(WT_cells, downsampled_KO_cells)]
DimPlot(WT_KO_integrated_downsampled, reduction = "umap", split.by ="orig.ident", ncol=2)

I have downsampled 4 times and recorded the number of cells in each cluster for all 4 versions (V1, V2, V3 and V4).

1) The downsampled percentage of cells in WT and KO is more over same compared to the actual % of cells in WT and KO

2) In each versions, I have highlighted the KO cells for cluster 1, 4, 5, 6 and 7 where the downsampled number is less than the WT cells. But before downsampling, if you see KO cells are higher compared to WT cells. However for cluster 0, 2 and 3 the trend is preserved.

Is it normal?

ADD REPLY • link modified 10 months ago • written 10 months ago by bioinforesearchquestions • 280

I mean, you're downsampling to a specific number of cells, so yes? Clusters 1, 4, 5, 6, 7 had closer to equal ratios of KO to WT cells, so losing cells from those groups means you're going to end up with fewer KO cells than WT cells.

ADD REPLY • link written 10 months ago by jared.andrews07 ♦ 8.3k

Yes, the collaborator requested equal number of cells on both the WT and KO samples. This is what I expected initially equal number of cells in each conditions under each cluster as well. The below image is based on the dummy data I manually entered.

But in the above image actual KO value vs downsampled KO value is fluctuating highly instead of downsampled versions V1, V2, V3, and V4.

From these different version of downsampled analyses, I am thinking of downsampling only on selected clusters like cluster 0, 1, 2, and 3 alone. Is there a way to do it?

You mentioned in your reply "If they just want to visualize them equally, I would just grab equal numbers of cells from each condition for visualization and feed them to the cells parameter for FeaturePlot or cells.use parameter for DimPlot."

Does the code which I provided above carry out that taks? Do you have any other way to do it?

ADD REPLY • link modified 10 months ago • written 10 months ago by bioinforesearchquestions • 280

So your collaborator wants equal numbers of cells per cluster? Seems silly to me.

Presumably you've already done what I suggested by performing your manual subsetting - your UMAP DimPlot should only contain the subsetted cells. My suggestion was just to use the parameters from the plotting functions to indicate which cells should be used from your full data set. Two approaches, should result in the same outcome.

This still feels like a job for violin or box plots. Particularly with dittoSeq, where you can hide the individual points and just show the box or violin plot. Then you could use your full dataset without your collaborator complaining about unequal cell numbers (why exactly do they want this?).

ADD REPLY • link written 10 months ago by jared.andrews07 ♦ 8.3k

Similar posts • Search »