单细胞marker基因表达密度图-(还有一个包装函数）

有小伙伴说想要做单细胞marker基因表达密度图，我一想，好像之前是做过的（单细胞marker基因可视化的补充---密度图与等高线图）。但是他又说没有文献中的效果。后来我一看，是因为着色的问题。其实用Nebulosa包（GitHub - powellgenomicslab/Nebulosa: R package to visualize gene expression data based on weighted kernel density estimation）就可以完成。

首先我们读入数据，分析下marker基因，用于作图：

setwd('D:\\KS项目\\公众号文章\\单细胞marker基因密度图')library(Seurat)
library(Nebulosa)
library(ggnetwork)
library(dplyr)markers <- FindAllMarkers(mouse_data, logfc.threshold = 0.5, min.pct = 0.5)
markers_plot <- markers %>% group_by(cluster) %>%slice(1:2)

需要注意，这里我算是投机取巧了，用了一个和这个作图毫无关系的包ggnetwork，主要的目的是用他的theme_blank()主题，一次性将所有清空，当然也可以自己慢慢写，不过有现成的，为啥不用呢？单个图看看：

plot_density(mouse_data, features = c("Ltf"),pal = 'magma', raster = T, size = 0.8) &theme_blank()&#这个是ggnetwork里面的主题，可以去掉所有ggplot背景theme(legend.frame = element_rect(colour = "black"),legend.ticks = element_line(colour = "black", linewidth  = 0),legend.key.width = unit(0.3, "cm"),legend.key.height = unit(0.8, "cm"),legend.title = element_text(color = 'black', face = "bold", size=8))

想要黑夜效果，背景设置呈黑色即可：

plot_density(mouse_data, features = c("Ltf"),pal = 'magma', raster = T, size = 0.8) &theme_blank()&theme(panel.background = element_rect(fill = "black"))&theme(legend.frame = element_rect(colour = "black"),legend.ticks = element_line(colour = "black", linewidth  = 0),legend.key.width = unit(0.3, "cm"),legend.key.height = unit(0.8, "cm"),legend.title = element_text(color = 'black', face = "bold", size=8))

我们借助Nebulosa，将上面的过程包装为一个函数，还是那句话，受累麻烦的事我们来解决，您轻松即可，先看看函数参数：需要注意的是，如果你只想框选某一类celltype，理想模式是这群celltype是单独聚类的，没有散落在其他类群，否则将会全部框选，效果不好，建议使用AI/PS添加。

很多是默认参数，我们演示一下：

#1黑夜模式，框选所有celltype轮廓
single_CM_density(object = mouse_data,features = "Ltf",night = T)

#2黑夜模式，只选择需要的celltype，我这里选择PMN(7），聚类好
#设置这个模式，需要设置参数idents
single_CM_density(object = mouse_data,features = "Ltf",night = T,idents = "celltype",label_celltype = "PMN(7)")

批量出一下黑夜模式的图：完美！这质感不就上去了嘛。

#批量出图
plist2 <- list()
for (i in 1:nrow(markers_plot)) {p = single_CM_density(object = mouse_data,features = markers_plot$gene[i],night = T)plist2[[i]] <- p
}#拼图
library(cowplot)
plot_grid(plotlist = plist2, ncol = 4)