如何用ggplot2和gggenes同时显示X轴基因标签与刻度条？

核心思路 #

用底部X轴显示连续的碱基对数值刻度，顶部X轴显示离散的基因标签，通过碱基对位置映射把基因标签精准对应到其在基因组上的位置。

步骤与代码示例 #

假设你的基因数据框结构包含start（基因起始碱基）、end（基因终止碱基）、gene（基因名，如geneA-geneF）、strand（链方向）等列。

1. 先准备基因标签位置映射表

# 提取每个基因的中点位置，用于放置标签
gene_labels <- your_gene_data %>% 
  distinct(gene, start, end) %>% 
  mutate(mid = (start + end)/2)

2. 绘制主图（底部X轴为碱基对）

library(ggplot2)

p <- ggplot() +
  # 绘制你的3个基因数据框，示例用geom_rect绘制基因区块
  geom_rect(data = gene_df1, aes(xmin = start, xmax = end, ymin = 0, ymax = 1, fill = "Group1")) +
  geom_rect(data = gene_df2, aes(xmin = start, xmax = end, ymin = 1.2, ymax = 2.2, fill = "Group2")) +
  geom_rect(data = gene_df3, aes(xmin = start, xmax = end, ymin = 2.4, ymax = 3.4, fill = "Group3")) +
  # 底部X轴：连续碱基对刻度
  scale_x_continuous(name = "Base Pairs", expand = c(0,0)) +
  scale_y_continuous(breaks = NULL, name = "") +
  theme_bw()

3. 添加顶部基因标签X轴

p + 
  scale_x_continuous(
    name = "Genes",
    breaks = gene_labels$mid,
    labels = gene_labels$gene,
    position = "top",
    expand = c(0,0)
  ) +
  # 调整主题优化显示
  theme(
    axis.title.x.top = element_text(margin = margin(b = 10)),
    axis.text.x.top = element_text(angle = 45, hjust = 0)
  )

关键说明 #

• 若基因位置不连续，确保gene_labels里的mid值准确对应每个基因的中心，标签会自动对齐。
• 若基因标签重叠，可调整axis.text.x.top的angle和hjust参数优化显示。
• 主图与顶部X轴的scale_x_continuous要保持expand参数一致，避免轴范围错位。

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