如何用R绘制等长基因邻域图？解决gggenes绘制问题

要实现所有基因矩形等长、按坐标顺序排列的需求，核心是放弃使用原始的start/end坐标，改为基于每个物种内基因的排序位置生成等长区间。具体步骤如下：

1. 预处理数据集 #

使用dplyr对数据按物种分组，再按原始start坐标排序，为每个基因分配顺序号并生成等长的xmin和xmax：

library(ggplot2)
library(gggenes)
library(dplyr)

# 按物种分组，基于start坐标排序，生成等长区间
example_genes <- example_genes %>%
  group_by(molecule) %>%
  arrange(start, .by_group = TRUE) %>% # 保证基因按原始坐标顺序排列
  mutate(
    gene_order = row_number(),
    xmin = gene_order - 0.5,  # 每个基因占据1个单位长度
    xmax = gene_order + 0.5
  ) %>%
  ungroup()

2. 绘制符合要求的基因邻域图 #

修改ggplot代码，使用新生成的区间，并满足其他需求：

ggplot(example_genes, aes(
  xmin = xmin, xmax = xmax, 
  y = molecule, 
  fill = label,  # 按label唯一着色
  forward = orientation  # 箭头方向由orientation决定
)) +
  geom_gene_arrow() +
  geom_gene_label(aes(label = gene)) +  # 在矩形中心标注基因名称
  facet_wrap(~ molecule, scales = "free_y", ncol = 1) +  # 物种纵向排列，y轴自适应
  scale_fill_discrete(palette = "Set3") +  # 可选：用Set3调色板，也可替换为其他离散配色
  theme_genes() +
  theme(axis.text.y = element_text(hjust = 1))  # 物种名称左对齐显示在左侧

关键修改说明 #

• 替换原始坐标：用gene_order生成的等长区间替代start/end，确保所有基因矩形大小一致、无间隙
• 着色逻辑：将fill = gene改为fill = label，实现按label唯一着色
• 标注优化：添加geom_gene_label自动将基因名放在矩形中心
• 布局调整：scales = "free_y"保证不同物种的y轴显示正常，theme调整让物种名称靠左显示

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