nature级别图表：一个注释气泡热图函数（适用于单细胞及普通数据）

之前我们写过一个函数，主要是展示单细胞比例图（nature级别图表：单细胞转录组细胞比例统计可视化函数），很多小伙伴购买函数，在使用过程却出现问题，主要的原因是我们在帖子里写的不够清楚，小伙伴也没有理解代码的意思，鉴于一个一个解释太费时间，我们决定写函数的帖子一律录制视频解说(视频在B站，搜索：KS科研分享与服务)，方便大家使用。所以，本次函数的使用我们已经制作好视频了，B站链接： 气泡热图注释函数_哔哩哔哩_bilibili需要购买函数的小伙伴可添加作者微信联系。微信VIP群已提前发布过，没有看到的小伙伴可查看记录。QQ群成员还是半价。

正式内容

我们写了一个作图函数Dotplot_anno()。首先写的初衷是为了展示单细胞marker基因，并对基因进行注释。但是后来我们将这个函数的功能扩大了，不仅仅使用在单细胞中，而且可以使用在普通基因表达气泡热图或者方块热图的使用上，并对需要的基因进行注释。或者在多组富集气泡图展示上，并对通路绘制相关的注释。而且，对legend给加上了边框。理解函数内容，学习更多。

函数的相关参数在开头已经进行详细了的注释。接下来就使用下看看效果！这个函数输入的数据单细胞是seurat对象，其他的数据则是ggplot可识别的长数据，如果是宽数据需要自己转化，对照我们的示例数据。

image.png

1、单细胞marker基因展示

首先我们运行单细胞数据看看效果：这里的features基因的顺序可以自定义，至于基因的分组要和group结合一致。

library(Seurat)
library(ggplot2)
library(dittoSeq)
library(viridis)
setwd('D:/KS项目/公众号文章/seurat中Dotplot注释函数')
source('./Dotplot_anno.R')
#单细胞的效果
DefaultAssay(uterus) <- "RNA"
markers <- c("ACTA2", "RGS5", #smooth muscle cells---7, 16
             "MS4A6A", "CD68","LYZ",#macrophages---13
             "CCL5", "STK17B","PTPRC",#lymphocytes---0,3,4,5,6,14,15,17,23,18,19
             "DCN", "COL6A3", "LUM",#stromal fibroblasts---2,20
             "PECAM1","PCDH17", "VWF",#endothelial cells---8,11,22
             "EPCAM", "CDH1",#(unciliated)epithelial cells---1,9,21
             "FOXJ1","CDHR3","DYDC2")#(ciliated)epithelial cells---10,12

Dotplot_anno(uterus, features = markers, celltype_color = dittoColors(),
             group = c(rep('SMC',2), rep('MAC',3),rep('Ly',3),
                       rep('SF',3), rep('EC',3),rep("UEC",2), rep('CEC',3)),
             color = colorRampPalette(c("navy","white","firebrick3"))(100),
             order = T)

image.png

2、基因表达热图 接下来看看bulk的效果，对于bulk表达数据，bulk基因表达量那一列列名要修改命名成exp，对于非单细胞的数据， 首先要进行设置的一个参数就是single_type=F，表明我们的对象不是单细胞对象。

df <- read.csv('Bulk.csv', header = T)
df$exp <- log10(df$exp)
#气泡图效果
Dotplot_anno(df, single_type = F, group = c(rep('inflammation',7),
                                            rep("chemotaxis",6),
                                            rep("cytotoxicity",5)),
             celltype_color = dittoColors(), 
             color = colorRampPalette(c("navy","white","firebrick3"))(100),
             order = F,heatmap = F,
             bulk_feature = 'features.plot',
             bulk_samples = 'id')

image.png

#热图效果
Dotplot_anno(df, single_type = F, group = c(rep('inflammation',7),
                                            rep("chemotaxis",6),
                                            rep("cytotoxicity",5)),
             celltype_color = dittoColors(), 
             color = colorRampPalette(c("navy","white","firebrick3"))(100),
             order = F,heatmap = T,
             bulk_feature = 'features.plot',
             bulk_samples = 'id')

image.png

3、富集气泡图

富集分析则是需要将log（P）那一列列名修改为exp，而且gene ratio那一列的列名必须修改为ratio。当然，首先要进行设置的一个参数就是single_type=F。

#看一下富集效果
df2 <- read.csv('df_enrich.csv',header = T)
df2$exp <- -log10(df2$exp)
Dotplot_anno(df2, single_type = F, group = c(rep('immune',5),
                                            rep("Metabolism",5),
                                            rep("stress",5),
                                            rep("antigen",5)),
             celltype_color = dittoColors(), 
             color = viridis(256),
             order = F,
             bulk_feature = 'features.plot',
             bulk_samples = 'id',)

image.png

修改下顺序，效果会好点。

#自定义顺序
Dotplot_anno(df2, single_type = F, group = c(rep('immune',5),
                                             rep("Metabolism",5),
                                             rep("stress",5),
                                             rep("antigen",5)),
             celltype_color = dittoColors(), 
             color = viridis(256),
             order = T,
             bulk_feature = 'features.plot',
             bulk_samples = 'id', level=c('NK',"CD14+ Mono","DC","FCGR3A+ Mono"))

image.png

这就是这个函数的作用了，使用起来还是可以的，希望对您的学习有用。如果觉得分享有用点个赞再走呗！