如何在gggenes中缩短长基因间区并添加文献截断标记？

我正在使用gggenes可视化一系列基因簇，希望找到内置功能或解决方案，缩短目标基因之间的长无注释区域。相关基因数据如下：

df<-data.frame(start=c(594198,596540,598457,600085,983488,984345),
stop=c(596450,598423,600070,601182,984336,986495),
species=rep("Ferriphaselus amnicola",6),
gene=c("gene1","gene2","gene3","gene4","gene5","gene6"))

使用如下gggenes代码可视化时，效果不佳：

ggplot(df, aes(xmin = start, xmax = stop, y = species, fill = gene)) +
    geom_gene_arrow() +
    facet_wrap(~ species, scales = "free", ncol = 1) +
    scale_fill_brewer(palette = "Set3") +
    theme_genes()

理想效果是当两个基因间的核苷酸数超过阈值x时，用文献常用的双斜线替代该基因组区段，类似PowerPoint修改后的效果。请问在gggenes或其他R包中是否有简便实现方法？

gggenes本身没有内置的“压缩长间隔+添加双斜线”功能，但可以通过数据预处理+自定义绘图元素实现需求，步骤如下：

1. 数据预处理：计算基因间隔并压缩长区段 #

首先计算相邻基因的间隔，设定阈值（示例设为10000bp），对超过阈值的间隔进行压缩，同时记录需要标记双斜线的位置：

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

# 设定间隔阈值
gap_threshold <- 10000
# 按物种排序基因（确保顺序正确）
df_sorted <- df %>% arrange(species, start)

# 计算每个基因与下一个基因的间隔，生成调整参数
df_sorted <- df_sorted %>%
  group_by(species) %>%
  mutate(next_start = lead(start),
         gap_length = next_start - stop,
         is_long_gap = gap_length > gap_threshold,
         # 长间隔只保留500bp可视长度，其余部分转为偏移量
         offset = ifelse(is_long_gap, gap_length - 500, 0),
         cumulative_offset = cumsum(ifelse(is.na(offset), 0, offset))) %>%
  ungroup()

# 调整基因的坐标
df_adj <- df_sorted %>%
  mutate(start_adj = start - cumulative_offset,
         stop_adj = stop - cumulative_offset)

# 提取需要添加双斜线的位置（长间隔的中间点）
gap_markers <- df_sorted %>%
  filter(is_long_gap) %>%
  mutate(x = (stop + next_start)/2 - cumulative_offset,
         y = species)

2. 可视化：绘制基因箭头+添加双斜线标记 #

用调整后的坐标绘制基因，然后在压缩的间隔位置添加双斜线：

ggplot(df_adj, aes(xmin = start_adj, xmax = stop_adj, y = species, fill = gene)) +
  geom_gene_arrow() +
  # 添加双斜线文本标记
  geom_text(data = gap_markers, aes(x = x, y = y, label = "//"), 
            size = 5, color = "gray50", fontface = "bold") +
  facet_wrap(~ species, scales = "free", ncol = 1) +
  scale_fill_brewer(palette = "Set3") +
  theme_genes() +
  labs(x = "Position (compressed)")

可选优化：用线段绘制斜线 #

如果需要更精细的斜线样式，可替换geom_text为两条交叉线段：

ggplot(df_adj, aes(xmin = start_adj, xmax = stop_adj, y = species, fill = gene)) +
  geom_gene_arrow() +
  # 绘制交叉斜线
  geom_segment(data = gap_markers, aes(x = x - 200, xend = x + 200, y = y - 0.1, yend = y + 0.1), 
               color = "gray50", linewidth = 1) +
  geom_segment(data = gap_markers, aes(x = x - 200, xend = x + 200, y = y + 0.1, yend = y - 0.1), 
               color = "gray50", linewidth = 1) +
  facet_wrap(~ species, scales = "free", ncol = 1) +
  scale_fill_brewer(palette = "Set3") +
  theme_genes() +
  labs(x = "Position (compressed)")

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