还记得之前一直有同学问，如何在一张图上同时展示两个基因的表达情况，方便比较，今天我们来实现一下

library(Seurat)
Seurat_obj = readRDS('rds文件')
FeaturePlot(pbmc, features = c('BPREB3','HLA-DRA'), blend = TRUE)

SCS2502_Cluster_5_co-expression.png

这里我随便挑选了两个基因，但是大家在做分析的时候，这两个基因最好是细胞通讯得到的配受体对，这样就可以在一张UMAP图上看到配体和受体的表达关系了。

空间图上同时展示两个基因的表达情况

library(Seurat)
Seurat_obj = readRDS('rds文件')
coor = data.frame('UMAP_1'= Seurat_obj@images$image@coordinates$imagerow , 'UMAP_2' = Seurat_obj@images$image@coordinates$imagecol)
rownames(coor) = colnames(Seurat_obj)
coor = as.matrix(coor)
Seurat_obj@reductions$umap@cell.embeddings = coor  ###空间位置的坐标替换了UMAP坐标
FeaturePlot(Seurat_obj, features = c('Pecam1','Cd38'), blend = TRUE)  ##这里选择了配受体对

图片.png

效果不错吧

生活很好，等你超越