一行代码绘制高分SCI火山图

一、概述

在近半年中，我读了很多的高分SCI文章，很多文章中都有多种不同的火山图，包括「普通的火山图、渐变火山图、以及包含GO通路信息的火山图」！

经过一段时间的文献阅读和资料查询，终于找到了一个好用而且简单的包——ggVolcano， 它是一个基于R语言和ggplot2绘图包开发的生物信息学数据可视化工具。它可以用于绘制火山图（Volcano plot），帮助研究者分析高通量实验数据，如基因表达谱或蛋白质组学数据，以识别差异表达或差异富集的基因或蛋白质。

二、数据集

1. 安装及使用

由于无法直接安装ggVolcano，我们需要下载devtools,通过devtools直接从github下载，代码如下：

# 下载器
install.packages("devtools")
library(devtools)
# 火山图
devtools::install_github("BioSenior/ggvolcano")
library(ggvolcano)
# 拼图
install.packages("patchwork")
library(patchwork)
# 渐变色色彩修改
install.packages("RColorBrewer")
library(RColorBrewer)

2. 读取数据

首先，我们需要读取aSAH数据集，查看数据集信息

data(deg_data)
head(deg_data,3)

数据集展示

        row  baseMean log2FoldChange     lfcSE      stat       pvalue         padj
GCR1   GCR1 7201.5782       2.244064 0.2004959  11.19256 4.434241e-29 2.153711e-25
OPI10 OPI10 1009.4171      -2.257454 0.2096469 -10.76789 4.880607e-27 1.185255e-23
AGA2   AGA2  249.1173       3.829474 0.3623263  10.56913 4.143136e-26 6.707736e-23

三、基础用法

ggVolcano绘制的火山图可以根据指定的差异表达或统计显著性阈值，将基因或蛋白质标记为显著差异的或富集的（例如p值小于0.05或fold change大于2），并且根据其差异表达或显著性水平值的高低，用颜色或大小进行编码。此外，它还支持在图中显示用户自定义的标注、注释或基因集富集分析结果，以帮助观察者更好地理解和解释结果。该包主要配置了三个函数：

1. 普通火山图：ggvolcano

• 方法参数解释：

  result: DEG结果数据。
  x: 对应于x轴的列名，默认为“log2FoldChange”。
  y: 对应于y轴的列名，默认为“padj”。
  pointSize: 点的大小。
  pointShape: 点的形状。
  fills: 包含点的填充颜色的向量。
  colors: 包含点的画笔颜色的向量。
  x_lab: x轴的标签。
  y_lab: y轴的标签。
  legend_title: 图例的标题。
  legend_position: 图例的位置。您可以从“UL”–左上，“UR”–右上，“DL”–左下和“DR”–右下中选择一个位置。
  log2FC_cut: log2FC的截止值。
  FDR_cut: FDR的截止值。
  add_line: 一个逻辑值，表示是否添加虚线，默认值为TRUE。
  add_label: 一个逻辑值，表示是否添加基因标签，默认值为TRUE。
  label: 对应于标签的列名。
  label_number: 您想在图中显示多少基因标签。
  custom_label: 包含您感兴趣的基因名称的向量，您想要添加到图中。
  output: 一个逻辑值，表示是否保存图片，默认值为TRUE。
  filename: 如果output=TRUE，请设置一个文件名。

• 代码演示

  data <- add_regulate(deg_data, log2FC_name="log2FoldChange",
                    fdr_name = "padj",log2FC = 1, fdr = 0.05)
  ggvolcano(data, x = "log2FoldChange", y = "padj",
          label = "row", label_number = 10, output = FALSE)

  
• 火山图样式修改

p1 <- ggvolcano(data, x = "log2FoldChange", y = "padj",
          fills = c("red","green","blue"),
          colors = c("red","green","blue"),
          label = "row", label_number = 8, output = FALSE)

p2 <- ggvolcano(data, x = "log2FoldChange", y = "padj",
          label = "row", label_number = 8, output = FALSE,legend_position="DR")+
  ggsci::scale_color_aaas()+
  ggsci::scale_fill_aaas()
p1|p2

2. 渐变色火山图：gradual_volcano

• 方法参数解释：

  result: DEG结果数据。
  x: 对应于x轴的列名，默认为“log2FoldChange”。
  y: 对应于y轴的列名，默认为“padj”。
  pointSize: 点的大小。
  pointShape: 点的形状。
  fills: 包含点的填充颜色的向量。
  colors: 包含点的画笔颜色的向量。
  x_lab: x轴的标签。
  y_lab: y轴的标签。
  legend_title: 图例的标题。
  legend_position: 图例的位置。您可以从“UL”–左上，“UR”–右上，“DL”–左下和“DR”–右下中选择一个位置。
  log2FC_cut: log2FC的截止值。
  FDR_cut: FDR的截止值。
  add_line: 一个逻辑值，表示是否添加虚线，默认值为TRUE。
  add_label: 一个逻辑值，表示是否添加基因标签，默认值为TRUE。
  label: 对应于标签的列名。
  label_number: 您想在图中显示多少基因标签。
  custom_label: 包含您感兴趣的基因名称的向量，您想要添加到图中。
  output: 一个逻辑值，表示是否保存图片，默认值为TRUE。
  filename: 如果output=TRUE，请设置一个文件名。

• 代码演示

gradual_volcano(deg_data, x = "log2FoldChange", y = "padj",
                label = "row", label_number = 10, output = FALSE,legend_title="")

• 火山图样式修改

这里不能直接修改颜色了，需要引入RColorBrewer包，用来调整色彩。

p1 <- gradual_volcano(data, x = "log2FoldChange", y = "padj",
          fills = brewer.pal(5, "RdYlBu"),
          colors = brewer.pal(10, "RdYlBu"),
          label = "row", label_number = 10, output = FALSE,legend_title="")

p2 <- gradual_volcano(data, x = "log2FoldChange", y = "padj",
          label = "row", label_number = 10, output = FALSE,legend_title="")+
  ggsci::scale_color_gsea()+
  ggsci::scale_fill_gsea()
p1|p2

3. GO通路火山图：term_volcano

参数和样式修改方法和上述两个方法类似，不再赘述，接下来展示基础使用方法：

data("term_data")

term_volcano(deg_data, term_data,
             x = "log2FoldChange", y = "padj",
             label = "row", label_number = 10, output = FALSE,
             legend_background_fill = "green")

本文由博客一文多发平台 OpenWrite 发布！