使用gggenes绘制长基因组区域：能否拆分为多行展示？

可以实现，核心思路是预处理数据，将长基因组拆分为多个子区域，给每个子区域分配独立的y轴分组，再调整基因的x坐标为子区域内的相对位置，最后用gggenes绘制。

具体步骤与代码示例 #

假设你的genes数据包含start（基因起始位置）、end（基因终止位置）、gene（基因名称）、molecule（分子名称）列，以下是实现代码：

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

# 定义每个分段的碱基长度，根据你的基因组实际情况调整
chunk_size <- 10000

# 预处理数据：拆分基因组为多个分段，计算相对位置
genes_processed <- genes %>%
  # 计算每个基因所属的分段编号
  mutate(segment = floor((start - min(start)) / chunk_size) + 1) %>%
  # 生成新的y轴分组（原分子名+分段编号）
  mutate(molecule_segment = paste(molecule, "part", segment, sep = "_")) %>%
  # 计算每个分段内的相对起始/终止位置，让每个分段的起点对齐
  group_by(molecule_segment) %>%
  mutate(rel_start = start - min(start),
         rel_end = end - min(start)) %>%
  ungroup()

# 绘制拆分后的基因图谱
ggplot(genes_processed, aes(xmin = rel_start, xmax = rel_end, y = molecule_segment, label = gene)) +
  geom_gene_arrow() +
  # 可选：添加每个分段的实际位置范围标签
  geom_text(aes(x = (max(rel_end) + min(rel_start))/2, y = molecule_segment, 
                label = paste("区间:", min(start), "-", max(start))), 
            vjust = 1.5, size = 3) +
  theme_gggenes()

说明 #

• 你可以根据需求调整chunk_size（分段长度），也可以改为按基因数量拆分（比如每20个基因一段）。
• 如果需要保留原分子的分组样式，可以在molecule_segment的命名上做调整，或者用facet分面配合，但直接生成新y轴分组的方式更直观。

内容的提问来源于stack exchange，提问作者Antonio J. Pérez