如何用R基于Galaxy的DESeq2数据绘制甲虫雌雄基因表达堆叠条形图

现有代码问题 #

你现有代码的核心问题有两个：

1. geom_bar默认使用stat="count"统计观测数量，没有将log2fc映射到Y轴，无法展示表达量数值
2. 手动拼接的gender-log2fc是离散分组变量，无法实现连续的性别偏向颜色渐变效果

实现步骤 #

1. 数据准备 #

确认数据包含三个核心字段即可，无需额外拼接分组字段：

• chromosomes：染色体编号
• sex：性别分组（雄性/雌性）
• log2fc：DESeq2输出的差异表达倍数log2转换值

2. 核心可视化代码 #

根据你的展示需求选择对应实现方式：

方案一：展示每个染色体上所有基因的log2FC分布（推荐用于表达范围展示）

# 加载依赖包
library(ggplot2)

# 读取数据，你原有剪贴板读取逻辑可保留
sexlinked <- read.table("clipboard", header = TRUE)
df <- sexlinked

# 绘图
ggplot(df, aes(x = chromosomes, y = log2fc)) +
  # 箱线图展示表达分布，fill映射log2fc实现渐变
  geom_boxplot(aes(fill = log2fc), colour = "grey30", outlier.size = 0.8) +
  # 配置渐变配色：蓝色对应雄性高表达，橙色对应雌性高表达
  scale_fill_gradient2(low = "#1f77b4", 
                       mid = "lightgrey", 
                       high = "#ff7f0e", 
                       midpoint = 0, 
                       name = "表达偏向性\n蓝=雄高 橙=雌高") +
  # 坐标轴与主题设置
  labs(x = "染色体", y = "表达量 log2(FC)") +
  theme_bw()

如果需要同时展示单个基因的表达点位，可以在geom_boxplot后追加一行：
geom_jitter(size = 0.6, alpha = 0.2, width = 0.2)

方案二：展示每个染色体雌雄的平均log2FC条形图

# 加载依赖包
library(ggplot2)
library(dplyr)

# 读取数据
sexlinked <- read.table("clipboard", header = TRUE)
df <- sexlinked

# 按染色体+性别分组计算平均表达量，若你的数据已聚合为单条记录可跳过该步
df_summary <- df %>%
  group_by(chromosomes, sex) %>%
  summarise(mean_log2fc = mean(log2fc), .groups = "drop")

# 绘图
ggplot(df_summary, aes(x = chromosomes, y = mean_log2fc, fill = sex)) +
  geom_col(position = position_dodge(0.8), colour = "grey30") +
  # 配置性别对应渐变配色
  scale_fill_manual(values = c("雄性" = "#1f77b4", "雌性" = "#ff7f0e"),
                    name = "性别") +
  labs(x = "染色体", y = "平均表达量 log2(FC)") +
  theme_bw()

3. 自定义调整说明 #

• 可修改scale_fill_gradient2的low/high参数替换渐变两端的颜色
• 若你的log2FC是雄性比雌性的计算结果，将midpoint对应颜色逻辑反过来即可
• 可通过调整geom_boxplot/geom_col的width参数修改图形宽度

内容的提问来源于stack exchange，提问作者datanerd90