阿狗工具

作者，追风少年i

最近一直有人问我空间转录组的饼图该如何绘制，如下图

图片.png

之前呢，说过一些，大家可以借鉴

10X空间转录组绘图之细胞类型百分比饼图和空间密度分布图

其中stlearn提供了封装好的函数，不过美观度上不够好。

相信大多数还是利用Seurat来进行联合分析，如果用了其他软件，可以将联合的矩阵结果进行替换，我们还是以之前的结果为例，HE染色和单细胞空间联合分析的结果部分展示如下

图片.png

图片.png

图片.png

画出来的效果如下图，颜色根据自己的喜好进行搭配

图片.png

好了，我们要开始了，一定要记住，前面的单细胞空间联合要亲自做

library(Seurat)
library(ggplot2)

####读取数据
cortex_sp <- readRDS(spatial.rds)
####单细胞空间联合的结果，这个如果已经导入到了rds，就不用再读取了
decon_mtrx = read.csv(sp.predictions.csv,row.names = 1)

####数据处理
cell_types_all <- colnames(decon_mtrx)[which(colnames(decon_mtrx) != "max")]

decon_df <- decon_mtrx %>%
  data.frame(check.names = F) %>%
  tibble::rownames_to_column("barcodes")
###添加比例到meta
cortex_sp@meta.data <- cortex_sp@meta.data %>%
  tibble::rownames_to_column("barcodes") %>%
  dplyr::left_join(decon_df, by = "barcodes") %>%
  tibble::column_to_rownames("barcodes")

###plot dot
slice <- names(cortex_sp@images)[1]
metadata_ds <- data.frame(cortex_sp@meta.data)
colnames(metadata_ds) <- colnames(cortex_sp@meta.data)
cell_types_interest <- cell_types_all

metadata_ds <- metadata_ds %>% tibble::rownames_to_column("barcodeID") %>%
            dplyr::mutate(rsum = base::rowSums(.[, cell_types_interest,
                drop = FALSE])) %>% dplyr::filter(rsum != 0) %>%
            dplyr::select("barcodeID") %>% dplyr::left_join(metadata_ds %>%
            tibble::rownames_to_column("barcodeID"), by = "barcodeID") %>%
            tibble::column_to_rownames("barcodeID")

spatial_coord <- data.frame(cortex_sp@images[[slice]]@coordinates) %>%
        tibble::rownames_to_column("barcodeID") %>% dplyr::mutate(imagerow_scaled = imagerow *
        cortex_sp@images[[slice]]@scale.factors$lowres, imagecol_scaled = imagecol *
        cortex_sp@images[[slice]]@scale.factors$lowres) %>% dplyr::inner_join(metadata_ds %>%
        tibble::rownames_to_column("barcodeID"), by = "barcodeID")

head(spatial_coord)

图片.png

就这样，得到了每个spot的坐标和细胞类型的比例，接下来就要绘图了，

###读取图片信息，这里读取Spaceranger分析的tissue_lowres_image.png即可
img <- png::readPNG(img)

###网格化
img_grob <- grid::rasterGrob(img, interpolate = FALSE, width = grid::unit(1,
        "npc"), height = grid::unit(1, "npc"))

####百分比饼图
scatterpie_pie <- suppressMessages(ggplot2::ggplot() + ggplot2::annotation_custom(grob = img_grob,
        xmin = 0, xmax = ncol(img), ymin = 0, ymax = -nrow(img)) +
        scatterpie::geom_scatterpie(data = spatial_coord, ggplot2::aes(x = imagecol_scaled,
            y = imagerow_scaled), cols = cell_types_all, color = NA,
            alpha = 1, pie_scale = 0.4) +
        ggplot2::scale_y_reverse() + ggplot2::ylim(nrow(img),
        0) + ggplot2::xlim(0, ncol(img)) + cowplot::theme_half_open(11,
        rel_small = 1) + ggplot2::theme_void() + ggplot2::coord_fixed(ratio = 1,
        xlim = NULL, ylim = NULL, expand = TRUE, clip = "on"))

pdf('sample.pie.pdf')
print(scatterpie_pie)
dev.off()

MK.pie.png

图例和配色有点差，我们来调整一下，首先调节一下图例的大小

####百分比饼图
scatterpie_pie <- suppressMessages(ggplot2::ggplot() + ggplot2::annotation_custom(grob = img_grob,
        xmin = 0, xmax = ncol(img), ymin = 0, ymax = -nrow(img)) +
        scatterpie::geom_scatterpie(data = spatial_coord, ggplot2::aes(x = imagecol_scaled,
            y = imagerow_scaled), cols = cell_types_all, color = NA,
            alpha = 1, pie_scale = 0.4) +
        ggplot2::scale_y_reverse() + ggplot2::ylim(nrow(img),
        0) + ggplot2::xlim(0, ncol(img)) + cowplot::theme_half_open(11,
        rel_small = 1) + ggplot2::theme_void() + ggplot2::coord_fixed(ratio = 1,
        xlim = NULL, ylim = NULL, expand = TRUE, clip = "on") + ggplot2::theme(legend.key.size = unit(50,'pt'),legend.title = element_text(size = 60),legend.text = element_text(size = 40)))

png(paste(outdir,paste(sample,'pie.png',sep = '.'),sep = '/'),width = 9 * 200,height = 7 * 200,type = 'cairo')
print(scatterpie_pie)
dev.off()

MK.pie.png

接下来颜色搭配

defined_cols = c('#e6194b', '#3cb44b', '#ffe119', '#4363d8', '#f58231', '#911eb4', '#46f0f0', '#f032e6', '#bcf60c', '#fabebe', '#008080', '#e6beff', '#9a6324', '#fffac8', '#800000', '#aaffc3', '#808000', '#ffd8b1', '#000075', '#808080')

scatterpie_pie <- suppressMessages(ggplot2::ggplot() + ggplot2::annotation_custom(grob = img_grob,
        xmin = 0, xmax = ncol(img), ymin = 0, ymax = -nrow(img)) +
        scatterpie::geom_scatterpie(data = spatial_coord, ggplot2::aes(x = imagecol_scaled,
            y = imagerow_scaled), cols = cell_types_all, color = NA,
            alpha = 1, pie_scale = 0.35) +
        ggplot2::scale_y_reverse() + ggplot2::ylim(nrow(img),
        0) + ggplot2::xlim(0, ncol(img)) + cowplot::theme_half_open(11,
        rel_small = 1) + ggplot2::theme_void() + ggplot2::coord_fixed(ratio = 1,xlim = NULL, ylim = NULL, expand = TRUE, clip = "on") + ggplot2::theme(legend.key.size = unit(50,'pt'),legend.title = element_text(size = 60),legend.text = element_text(size = 40)) + scale_fill_manual(values = defined_cols[1:length(cell_types_all)]))

png(paste(outdir,paste(sample,'pie.png',sep = '.'),sep = '/'),width = 9 * 200,height = 7 * 200,type = 'cairo')
print(scatterpie_pie)
dev.off()

图片.png

好多了，借鉴一下第一张图的搭配

ggplot2::annotation_custom(grob = img_grob,
        xmin = 0, xmax = ncol(img), ymin = 0, ymax = -nrow(img)) +
        scatterpie::geom_scatterpie(data = spatial_coord, ggplot2::aes(x = imagecol_scaled,
            y = imagerow_scaled), cols = cell_types_all, color = NA,
            alpha = 1, pie_scale = 0.35) +
        ggplot2::scale_y_reverse() + ggplot2::ylim(nrow(img),
        0) + ggplot2::xlim(0, ncol(img)) + cowplot::theme_half_open(11,
        rel_small = 1) + ggplot2::theme_void() + ggplot2::coord_fixed(ratio = 1, xlim = NULL, ylim = NULL, expand = TRUE, clip = "on") + ggplot2::theme(legend.key.size = unit(50,'pt'),legend.title = element_text(size = 60),legend.text = element_text(size = 40)) + scale_fill_manual(values = defined_cols[1:length(cell_types_all)]))

图片.png

嗯，还行

来对比一下其他文章的图片形式

图片.png

好了，已经分享给大家了，生活很好，有你更好，百度文库出现了大量抄袭我的文章，对此我深表无奈，我写的文章，别人挂上去赚钱，抄袭可耻，挂到百度文库的人更可耻