单细胞分析这个cluster可以去除吗

Question

在跑单细胞时碰到这个cluster（图中的cluster29）nCount_RNA和nFeature_RNA和别的cluster都不一样，并且在umap图中位置比较偏，细胞数也少，请问这个是由于doublets未去除干净导致的吗？这个cluster可以删除吗

Answer 1

下午好🌅🌅🌅
本答案参考通义千问

在单细胞分析中，遇到一个与其它cluster差异显著、位置偏移、细胞数量少的cluster（如你提到的cluster29），确实需要仔细评估其是否为双细胞（doublets）或异常细胞。以下是详细的分析和解决方案：

---

一、为什么这个cluster可能存在问题？

1. nCount_RNA 和 nFeature_RNA 显著不同

   • 这是判断是否为双细胞的一个重要指标。通常，双细胞的UMAP图位置会偏离正常细胞群，且表达量较高（nCount_RNA 和 nFeature_RNA 偏高）。
   • 如果cluster29的nCount_RNA和nFeature_RNA明显高于其他cluster，可能是由于两个细胞被错误地合并成一个。

2. UMAP位置偏移

   • UMAP图中如果某个cluster的位置与其他cluster有明显分离，尤其是孤立存在，则很可能是双细胞或技术噪声。

3. 细胞数量少

   • 单个cluster中只有少数细胞，也可能是双细胞或数据异常。

---

二、这个cluster是否可以删除？

可以删除，但需谨慎操作，并确保不影响整体分析结果。

---

三、如何判断该cluster是否为双细胞？

方法1：使用 doubletScores 或 scDblFinder 检测双细胞

# 使用 scDblFinder 检测双细胞
library(scDblFinder)
dblf <- scDblFinder(seu, use.dimred = "UMAP")
plot(dblf)

# 查看哪些细胞被标记为双细胞
[e-mail]$doublet <- dblf$doublet

方法2：检查 nCount_RNA 和 nFeature_RNA 分布

# 绘制每个cluster的nCount_RNA分布
plot_data <- as.data.frame([e-mail])
plot_data$nCount_RNA <- as.numeric(plot_data$nCount_RNA)
ggplot(plot_data, aes(x = cluster, y = nCount_RNA, fill = cluster)) +
  geom_boxplot() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

---

四、解决方案（有序列表）

1. 确认cluster29是否为双细胞

   • 使用 scDblFinder 或 doubletScores 工具进行检测。
   • 观察其在UMAP中的位置和表达特征是否异常。

2. 检查 nCount_RNA 和 nFeature_RNA 的分布

   • 如果该cluster的数值明显高于其他cluster，可能是双细胞。

3. 若确定为双细胞或异常细胞，可将其从数据集中删除

   • 使用以下代码删除cluster29：

   # 删除cluster29
   seu <- subset(seu, subset = !is.in.categorical("cluster", "29"))
   

4. 重新运行UMAP或tSNE等降维分析

   • 确保删除后其他cluster的分布没有异常。

5. 重新进行聚类和注释

   • 确保后续分析不受影响，避免误判。

---

五、注意事项

• 不要盲目删除，建议先通过多种方法验证cluster29的性质。
• 保留足够多的细胞，以避免因删除过多细胞导致后续分析不可靠。
• 记录删除过程，便于后续复现和溯源。

---

六、总结

如果你发现某个cluster在UMAP中位置偏移、表达量异常且细胞数少，它很可能是双细胞或异常细胞。在这种情况下，你可以安全地将它从数据集中删除，但需先进行充分验证。

---

如有需要，我可以帮助你进一步分析你的数据集或提供具体的代码示例。欢迎继续提问！