Based on ggtree::theme_tree2(), customized for publication quality.
Includes coord_cartesian(clip = "off") to prevent collapsed clade
triangles from being clipped at the plot panel boundary — triangle vertices
(especially the MRCA node apex) often extend beyond the tip-based y-axis range.