复现《Cell》文章图表：气泡热图展示基因表达+分组添加

今天我们来复现一下一篇《cell》文章中的图标，图如下：

（Liver Immune Profiling Reveals Pathogenesis and Therapeutics for Biliary Atresia）

我们重点复现左边部分，右侧的作图在单细胞系列提到过：跟着Cell学单细胞转录组分析(十一):单细胞基因评分|AUCell评分，这里的复现不仅适用于单细胞，其他的数据也是可以做这样的图的。

这个图的重点在于将每个细胞群的各个样本分开展示基因表达，添加样本分组的注释。第二点在于x轴坐标标签设置。我复现的结果如下，颜色是看到网上一个比较好的色彩搭配所以使用了，基本要素都实现了，后期需要组图！

首先我们使用单细胞数据构建做图数据。其他数据只要做成ggplot气泡图读入的作图格式即可。

DefaultAssay(immune) <- "RNA"
markers <- FindAllMarkers(immune, logfc.threshold = 0.25, min.pct = 0.25, only.pos = T)
top5 <- markers %>% group_by(cluster) %>% top_n(5, avg_log2FC)
library(ggplot2)
p1 <- DotPlot(immune, features = top5$gene,
             cols = my36colors, group.by = "celltype", split.by = "orig.ident")+coord_flip()
exp <- p1$data
library(forcats)
exp$features.plot <- as.factor(exp$features.plot)
exp$features.plot <- fct_inorder(exp$features.plot)

更多内容请至我的公众号《KS科研分享与服务》，感谢支持！