This is a test version of Biostars. For the public version, visit https://www.biostars.org.
scRNA integration problem with cluster not being separated

For single cell analysis, seuratobojects were created for 6 samples and integration was performed after going through basic procedures. However, we confirmed that the clusters were not clearly divided. Can you tell me the cause and solution for this? The overall process followed the Seurat V5 process. However, even when following the Harmony or V4 process, there was no clear division. In the figure, the first figure is before integration, and the second and third figures are after integration.

Seurat <- merge(AS1, c(AS2, AS3, CO1, CO2, CO3))
Seurat <- NormalizeData(object = Seurat)
Seurat <- FindVariableFeatures(object = Seurat)
Seurat <- ScaleData(object = Seurat)
Seurat <- RunPCA(object = Seurat)
Seurat <- FindNeighbors(object = Seurat, dims = 1:15)
Seurat <- FindClusters(object = Seurat, resolution = 0.2)
Seurat <- RunUMAP(object = Seurat, dims = 1:15)
S1 <- DimPlot(Seurat, reduction = "umap", label = TRUE, group.by = "orig.ident")
Seurat_I <- IntegrateLayers(object = Seurat, method = CCAIntegration, orig.reduction = "pca", new.reduction = "integrated.cca",
                           verbose = FALSE)
Seurat_I[["RNA"]] <- JoinLayers(Seurat_I[["RNA"]])
Seurat_I <- FindNeighbors(Seurat_I, reduction = "integrated.cca", dims = 1:15)
Seurat_I <- FindClusters(Seurat_I, resolution = 0.3)
Seurat_I <- RunUMAP(Seurat_I, dims = 1:15, reduction = "integrated.cca")
S2 <- DimPlot(Seurat_I, reduction = "umap", group.by = "orig.ident")
S3 <- DimPlot(Seurat_I, reduction = "umap", label=T)
S1+S2+S3

enter image description here

seurat scrna

is that integrated or not depends on batch effect and biological variance in the data itseft? here you performed Harmony, meaning that you were trying to remove batch effect, then 'similar' clusters are now integrated.

So can we say that this data is clustered well?

are AS2, AS3, CO1, CO2, CO3 in same batch? if then harmony is not needed

0 answers

No answers yet.

Log in to answer this question.