如何缩短ggtree中由APE生成的phylo树的过长分支？

嘿，这个问题我在做种内种群遗传分析的时候也碰到过！外类群的超长分支确实会把种内的精细结构挤得没地方展示，下面分享几个我亲测有效的解决办法，都是用R里的ape和ggtree就能实现的：

方法1：直接修改树对象的分支长度（适合需要保留调整后树用于后续分析的场景） #

你可以手动定位外类群对应的分支，然后把它们的长度按比例缩小。步骤很简单：

library(ape)
library(ggtree)

# 假设你的系统发育树对象是tree，外类群物种名为"species_out1"和"species_out2"
# 找到外类群对应的末端分支索引
outgroup_tips <- which(tree$tip.label %in% c("species_out1", "species_out2"))
# 按比例缩小分支长度，这里用0.1系数，你可以根据实际情况调整
tree$edge.length[outgroup_tips] <- tree$edge.length[outgroup_tips] * 0.1

# 如果外类群内部还有长分支，还可以找到它们的共同祖先节点调整内部分支
outgroup_mrca <- getMRCA(tree, c("species_out1", "species_out2"))
outgroup_internal_edges <- which(tree$edge[,2] == outgroup_mrca)
tree$edge.length[outgroup_internal_edges] <- tree$edge.length[outgroup_internal_edges] * 0.1

# 绘图查看效果
ggtree(tree) + geom_tiplab()

方法2：绘图层面限制X轴范围（适合仅需优化可视化、不想修改原始树的场景） #

这个方法不需要改动树本身，只是通过调整绘图的视图范围，把超长的外类群分支“移出”可视区域，同时手动调整外类群的标签位置避免被截断：

# 先计算种内分支的最大长度，作为X轴上限的参考
intra_branch_max <- max(tree$edge.length[-outgroup_tips])

ggtree(tree) +
  # 给外类群标签设置右对齐，避免被X轴截断
  geom_tiplab(aes(hjust = ifelse(label %in% c("species_out1", "species_out2"), 1.1, 0))) +
  # 限制X轴范围，聚焦种内区域
  coord_cartesian(xlim = c(0, intra_branch_max * 1.5))

方法3：折叠外类群分支（最简洁的可视化方案） #

用ggtree的折叠功能把外类群的整个分支变成一个紧凑的三角块，既能保留根的定位作用，又不会占用大量空间：

# 找到外类群的最近共同祖先节点
outgroup_mrca <- getMRCA(tree, c("species_out1", "species_out2"))

ggtree(tree) +
  # 折叠外类群分支，设置填充色区分
  geom_collapse(node = outgroup_mrca, fill = "#a6bddb") +
  geom_tiplab()

方法4：精准调整特定分支长度（依赖phytools包） #

如果你想给外类群分支设置一个固定的长度，而不是按比例缩放，可以用phytools包的rescaleBranch函数：

library(phytools)

# 将外类群的共同祖先分支长度设为0.5（数值可自定义）
tree_rescaled <- rescaleBranch(tree, node = outgroup_mrca, length = 0.5)

ggtree(tree_rescaled) + geom_tiplab()

你可以根据自己的需求选择合适的方法：如果只是做图展示，方法2和3最方便；如果后续还要用调整后的树做分析，方法1或4更合适。

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