Draw and export phylogenetic trees with Bio.Phylo plus matplotlib, and route rich figures to ggtree, ETE4, or iTOL. Covers why a tree figure is an argument not a neutral picture, the cladogram-vs-phylogram-vs-chronogram choice that hides or reveals rate and time, how ladderization manufactures a false arrow of progress, why an unlabeled support number always flatters the result (bootstrap vs posterior vs SH-aLRT vs UFBoot are different scales), why Bio.Phylo silently drops BEAST HPD bars so annotated Bayesian trees must go through treeio plus ggtree, and the tip-count and raster-vs-vector thresholds for legible publication figures. Use when drawing a tree, choosing a layout, coloring branches, showing or labeling support, exporting vector figures, or deciding a drawing tool. Routes annotation-preserving reads to tree-io, and rooting and ladderizing to tree-manipulation.
Draw and export phylogenetic trees with Bio.Phylo plus matplotlib, and route rich figures to ggtree, ETE4, or iTOL. Covers why a tree figure is an argument not a neutral picture, the cladogram-vs-phylogram-vs-chronogram choice that hides or reveals rate and time, how ladderization manufactures a false arrow of progress, why an unlabeled support number always flatters the result (bootstrap vs posterior vs SH-aLRT vs UFBoot are different scales), why Bio.Phylo silently drops BEAST HPD bars so annotated Bayesian trees must go through treeio plus ggtree, and the tip-count and raster-vs-vector thresholds for legible publication figures. Use when drawing a tree, choosing a layout, coloring branches, showing or labeling support, exporting vector figures, or deciding a drawing tool. Routes annotation-preserving reads to tree-io, and rooting and ladderizing to tree-manipulation.
Before using code patterns, verify installed versions match. If versions differ:
Python: then to check signatures
pip show biopython
help(Phylo.draw)
R: packageVersion('ggtree') then ?ggtree to verify layout and geom arguments
If code throws ImportError, AttributeError, or TypeError, introspect the installed
package and adapt the example to match the actual API rather than retrying.
Bio.Phylo draw() renders rectangular layouts only and parses no BEAST/MrBayes annotation block; circular/fan/unrooted layouts and HPD bars require ggtree, ETE4, or iTOL. The ETE3->ETE4 migration changed APIs; introspect rather than assume ETE3 signatures.
Tree Visualization -- A Figure Is an Argument
"Draw and export my tree figure" -> Render a topology under a chosen layout and branch-length encoding, decide which support to show and how to label it, and export in a format that survives print.
Publication with metadata, rich support, or BEAST HPD bars: ggtree + treeio (R)
Scope: drawing, styling, annotating, and exporting tree figures, and choosing the drawing tool. Reading/converting files and preserving [&...] annotations before plotting -> tree-io. Rooting, ladderizing, pruning, collapsing low-support nodes -> tree-manipulation. ggplot2 grammar, themes, and ggsave underlying ggtree -> data-visualization/ggplot2-fundamentals. Composing a tree plus aligned panels into one figure -> data-visualization/multipanel-figures.
The Single Most Important Modern Insight
A phylogenetic tree figure is a rhetorical and modeling artifact, not a neutral picture of the data. A Newick or Nexus file encodes a topology, optionally branch lengths, optionally one node label; it does NOT encode which child is drawn on top, whether the horizontal axis means evolutionary distance or just nesting depth, where the root sits, the aspect ratio, or which of several support measures a bare "95" refers to. Every one of those is supplied by the person drawing the tree, so the figure is the data filtered through a dozen interpretive choices the reader cannot see. Four argument-level decisions dominate, and the same file can be drawn to tell mutually contradictory stories that are all technically correct renderings of one topology:
The geometry is a claim about what the lengths mean. A cladogram asserts only nesting; a phylogram asserts nesting plus amount of change (substitutions per site); a chronogram asserts timing and demands a clock model plus calibrations. Drawing a phylogram as a cladogram erases rate variation, the very signal that flags long-branch attraction; drawing untrusted branch lengths as a phylogram fabricates a quantitative claim.
Ladderization implies a directionality that is not there. Sibling subclades are unordered, so any node rotation is the same tree, but the eye reads top-to-bottom ordering as a march of progress toward the bottom-most tip. Rotation can also bury non-monophyly: a paraphyletic group can be rotated to look contiguous, or a clean clade made to look scattered.
Which support is shown, and whether it is labeled, changes the conclusion. Support is not accuracy, and the measures are not interchangeable: a posterior of 0.95 is generally weaker evidence than a bootstrap of 95 for the same data, and UFBoot 95, SH-aLRT 80, and TBE live on different scales. A bare integer with no legend is the quietest lie because the reader assumes bootstrap and the ambiguity always flatters the result.
The drawing tool is decided by the I/O layer, not the surrounding language. If the tree carries BEAST/MrBayes HPD intervals and posteriors, only a tool whose data model can carry those annotations (treeio + ggtree) can plot them; Bio.Phylo flattens the tree to topology + length + one label and silently drops the uncertainty that was the result.
The skill's job is to make these choices conscious and DECLARE them in the caption, not to default into them. Default ladderization plus cladogram-ish rendering plus unlabeled nodes is the path of least resistance, and it is an argument the user did not intend to make.
Tool Taxonomy
The central decision is which drawing tool to reach for, and it is coupled to the I/O decision: a tool that flattens the tree to topology plus one label cannot draw BEAST HPD bars or posteriors it cannot represent. The plotting path through treeio + ggtree is the same I/O choice made in tree-io.
Tool (lang)
Citation
Layouts and annotation power
When to choose
Bio.Phylo + matplotlib (Py)
Talevich 2012
rectangular phylogram/cladogram only; one label per node via label_func/branch_labels; no circular/unrooted, no BEAST-block parse
a fast scripted look inside a Python pipeline; headless/CI rendering of simple trees; NOT publication figures with metadata or HPD bars
ggtree + treeio + ggtreeExtra (R)
Yu 2017; Wang 2020; Xu 2021
rectangular, slanted, circular, fan, radial/unrooted; %<+% attaches a data.frame, geom_tiplab/geom_nodelab/geom_cladelab, gheatmap, geom_facet (Cartesian) or geom_fruit rings (circular); treeio reads BEAST/MrBayes/IQ-TREE so HPD bars and posterior plot directly
the publication standard whenever a figure integrates metadata, rich/dual support, HPD bars, or aligned panels; the only clean path for BEAST HPD bars and posterior
ETE4 (Py; ETE3 paper)
Huerta-Cepas 2016
rectangular and circular; programmatic per-node NodeStyle/faces and layout functions; integrated NCBI taxonomy
styling thousands of nodes by rule, taxonomy-driven annotation, or rendering many trees headlessly in Python
iTOL v6 (web)
Letunic 2024
rectangular, circular, unrooted, very large trees; tab-delimited dataset templates for color strips, heatmaps, bars, binary symbols, clade collapse; SVG/PDF/PNG/EPS export
large trees and template-based metadata with no code; collaborative web-driven polished figures
FigTree (desktop GUI)
Rambaut (software)
rectangular, polar, radial; reads BEAST NEXUS node annotations to display HPD bars and posterior on screen; manual rotate/collapse/color
interactive inspection of a single BEAST/MrBayes tree before scripting the final figure; not reproducible, use for exploration not the pipeline
Decision rule: large tree or no-code web annotation -> iTOL v6. Publication figure with metadata, dual support, HPD bars, or aligned heatmaps -> ggtree + treeio (+ ggtreeExtra for circular rings). Programmatic styling of thousands of nodes or NCBI-taxonomy annotation -> ETE4. Fast topology look from a Python script -> Bio.Phylo. Interactive rooting/inspection of one BEAST tree -> FigTree, then move the reproducible figure to ggtree. The recurring mistake is defaulting to Bio.Phylo because the pipeline is Python, then discovering it cannot draw a circular layout, place a heatmap, or show the BEAST HPD bars that are the whole point.
What Each Layout Reveals and Hides
Layout
Reveals
Hides or distorts
Good for
Bad for
Rectangular phylogram
branch lengths read directly off a linear axis; honest distance comparison
tip labels overplot past a few hundred tips; uses vertical space
the default when length matters and tip count is moderate (<~150)
very large trees
Slanted / triangular
compactness; quick topology read
the diagonal implies a ladder/progression; imprecise length reading
quick schematic topology views
precise branch-length comparison
Circular / fan
hundreds to thousands of tips in one panel; metadata rings
radial distortion -- the same length subtends a larger arc near the rim, so distance is perceptually compressed near the root; near-root structure cramped
big trees where the message is broad structure plus aligned metadata
quantitative branch-length comparison
Unrooted / radial
honestly shows NO assumed root or time direction; overall shape and long-branch outliers
no time direction; clade membership harder to trace; clusters over-read as clades
exploratory views, "no root committed"
directional stories; any "basal/early-diverging" claim
Chronogram (time-tree)
node timing on a time axis
requires a clock model and calibrations the reader must trust; without HPD bars timing looks falsely precise
dated BEAST/treePL/MCMCtree trees where timing is the message
any tree lacking a clock; never draw without age uncertainty bars
Circular is still a ROOTED tree bent into a ring (center = root = past); unrooted/radial explicitly refuses a root. Presenting an unrooted radial tree and narrating "X is basal" is a contradiction.
Bio.Phylo + matplotlib Recipes
Quick text and ASCII inspection, no figure needed:
from Bio import Phylo
tree = Phylo.read('tree.nwk', 'newick')
print(tree) # indented text summary
Phylo.draw_ascii(tree) # ASCII-art diagram, useful in a terminal or log
Draw to a vector file (always pass an axes and do_show=False for headless/scripted use):
from Bio import Phylo
import matplotlib.pyplot as plt
tree = Phylo.read('tree.nwk', 'newick')
tree.ladderize() # legibility only; ordering carries NO phylogenetic meaning -- say so in the caption
fig, ax = plt.subplots(figsize=(10, 8))
Phylo.draw(tree, axes=ax, do_show=False)
ax.set_title('Phylogenetic tree (phylogram, branch length = subs/site)')
fig.savefig('tree.pdf', bbox_inches='tight') # vector: text and lines stay sharp at any size
plt.close(fig)
Label tips, and show support with its measure named (never a bare integer):
deftip_only(clade):
return clade.name if clade.is_terminal() else''defsupport_label(clade):
# the measure MUST be stated in the legend/caption; here values are bootstrap percentagesifnot clade.is_terminal() and clade.confidence isnotNone:
returnf'{clade.confidence:.0f}'return''
fig, ax = plt.subplots(figsize=(12, 10))
Phylo.draw(tree, axes=ax, do_show=False, label_func=tip_only, branch_labels=support_label)
ax.set_title('Bootstrap support shown at internal nodes')
fig.savefig('supported_tree.svg', bbox_inches='tight')
plt.close(fig)
Color branches by group (convert to phyloXML for native color support):
from Bio.Phylo.PhyloXML import BranchColor
xtree = tree.as_phyloxml() # phyloXML carries branch color through draw()for clade in xtree.find_clades():
if clade.name and clade.name.startswith('Homo'):
clade.color = BranchColor.from_name('red')
fig, ax = plt.subplots(figsize=(10, 8))
Phylo.draw(xtree, axes=ax, do_show=False)
fig.savefig('colored_tree.pdf', bbox_inches='tight')
plt.close(fig)
Scale the panel to tip count so labels stay legible, and drop the axis frame:
For circular/fan/unrooted layouts, metadata heatmaps, dual support, or BEAST HPD bars, Bio.Phylo cannot help -- route to ggtree + treeio (R), ETE4, or iTOL.
Per-Method Failure Modes
Drawing Meaningless Branch Lengths as a Phylogram
Trigger: Rendering a constraint tree, a supertree, or a tree whose lengths are non-comparable with the horizontal axis proportional to length.
Mechanism: The phylogram geometry asserts a quantitative claim about evolutionary distance that the lengths do not support.
Symptom: A topology-only or arbitrary-length tree appears to make precise distance statements.
Fix: Draw as a cladogram and SAY so (branch.length='none' in ggtree); or, if lengths are real, keep them and state the unit (subs/site vs time) plus a scale bar.
Unlabeled or Mislabeled Support
Trigger: A node shows "98" with no legend, or two measures are printed without saying which is which.
Mechanism: The reader defaults to assuming bootstrap, but it may be a posterior (much weaker for the same number), an SH-aLRT (cutoff 80), a UFBoot (cutoff 95, not the BP-70 scale), or a TBE.
Symptom: A weakly resolved bush reads as a confident comb because the displayed number is over-read.
Fix: Always state the measure(s) and their order (e.g. "SH-aLRT/UFBoot" at each node); collapse nodes below threshold into polytomies rather than drawing fake resolution, since bootstrap, posterior, SH-aLRT, and UFBoot sit on different scales.
Tip-Label Overplotting on Large Trees
Trigger: A rectangular layout past a few hundred tips.
Mechanism: Horizontal labels collide into an unreadable black band; authors then shrink the font to illegibility or silently drop labels.
Symptom: Tip labels are unreadable or missing.
Fix: Rotate labels, switch to circular/fan with radial labels, collapse uninformative clades, annotate by colored strips/rings (ggtreeExtra, iTOL datasets) instead of text, or move to iTOL which is built for large trees.
Non-Monophyly Hidden by Ladderization or Rotation
Trigger: Rotating nodes to make a paraphyletic or polyphyletic group look contiguous (or a clade look scattered).
Mechanism: Node rotation is information-free with respect to the tree, but the eye reads contiguity as a clade and ordering as direction.
Symptom: A group that is not monophyletic appears unified, or a directional narrative is implied.
Fix: Color by group and let the topology speak; never narrate ordering as meaning; state "tips ladderized for legibility, ordering carries no phylogenetic meaning."
Chronogram Without Age Uncertainty
Trigger: Drawing a BEAST/MrBayes time-tree with point-estimate node ages and no HPD bars.
Mechanism: The 95% HPD intervals on node heights ARE the result; omitting them asserts false precision.
Symptom: Node ages look known to the day; a reviewer asks where the credible intervals went.
Fix: Draw the HPD bars via treeio read.beast() -> ggtree geom_range('height_0.95_HPD') (introspect the column name; it varies by source program and treeio version), or enable node bars in FigTree; never draw an annotated Bayesian tree with Bio.Phylo, which drops the annotations silently.
Raster Export for Publication
Trigger: Saving a tree as a 150-dpi PNG for a paper.
Mechanism: Raster pixelates thin branches and small tip labels at print size.
Symptom: Fuzzy labels and pixelated branches; journal rejection or blurry print.
Fix: Export SVG/PDF/EPS (vector) so text and lines stay sharp at any scale; only if forced to rasterize, do so at final size with >=600 dpi for line art.
Quantitative Thresholds
These are operational defaults for a standard portrait panel with ~6-8 pt labels; adapt to font, page, and journal specs.
Quantity
Threshold
Rationale / source
Tips, rectangular horizontal labels comfortable
<=~50
labels do not collide at legible font
Tips, start rotating/shrinking labels
~50-150
label height x tip count approaches panel height
Tips, switch to circular/fan or collapse clades
~150-500
rectangular labels collide; radial labels recover space
Tips, per-tip text impractical; use strips/rings
>~500-1000
annotate by color/groups (iTOL, ggtreeExtra); iTOL is built for large trees
Publication export
vector (SVG/PDF/EPS)
text and lines stay sharp at any scale; the default
Forced raster, line-art/tree
>=600 dpi (many journals 600-1200)
~300 dpi is the photo floor but too coarse for thin branches and small labels
Scale bar on any phylogram
mandatory (subs/site or a time axis)
without it the reader cannot recover true distances under aspect-ratio distortion
Bootstrap strong
>=95 (>=70 moderate)
Hillis-Bull heuristic; state the measure (Felsenstein 1985)
UFBoot strong
>=95 (not the BP-70 scale)
recalibrated; IQ-TREE joint rule SH-aLRT>=80 AND UFBoot>=95
SH-aLRT strong
>=80
IQ-TREE recommendation
Posterior probability strong
>=0.95, but WEAKER than bootstrap 95
PP systematically higher for the same data; do not equate
Lock the branch-length scale; do not let the figure engine non-uniformly stretch a phylogram to fill a fixed panel, which distorts the lengths the figure exists to communicate.
Common Errors
Error / symptom
Cause
Solution
BEAST HPD bars absent from a Python figure
drew an annotated tree with Bio.Phylo
route through treeio read.beast + ggtree geom_range
Figure not saving / blank
do_show=True opens a window instead of writing
pass do_show=False, then fig.savefig(...)
Branch colors not appearing
plain Newick tree has no color slot
convert with tree.as_phyloxml() and set clade.color
Labels overlap into a black band
too many tips for rectangular layout
increase panel height, rotate labels, or switch to circular/iTOL
"Basal" claim on a radial tree
narrated an unrooted layout as if rooted
root explicitly (tree-manipulation) and show the root before any directional claim
Support number misread
bare integer with no measure stated
label the measure(s) and order; collapse sub-threshold nodes to polytomies
Phylo.draw has no circular option
Bio.Phylo is rectangular-only
use ggtree layout='circular', ETE4, or iTOL
References
Talevich E, Invergo BM, Cock PJA, Chapman BA. 2012. Bio.Phylo: a unified toolkit for processing, analyzing and visualizing phylogenetic trees in Biopython. BMC Bioinformatics 13:209.
Yu G, Smith DK, Zhu H, Guan Y, Lam TT-Y. 2017. ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods in Ecology and Evolution 8(1):28-36.
Yu G. 2020. Using ggtree to visualize data on tree-like structures. Current Protocols in Bioinformatics 69(1):e96.
Wang L-G, Lam TT-Y, Xu S, Dai Z, Zhou L, Feng T, Guo P, Dunn CW, Jones BR, Bradley T, Zhu H, Guan Y, Jiang Y, Yu G. 2020. treeio: an R package for phylogenetic tree input and output with richly annotated and associated data. Molecular Biology and Evolution 37(2):599-603.
Xu S, Dai Z, Guo P, Fu X, Liu S, Zhou L, Tang W, Feng T, Chen M, Zhan L, Wu T, Hu E, Jiang Y, Bo X, Yu G. 2021. ggtreeExtra: compact visualization of richly annotated phylogenetic data. Molecular Biology and Evolution 38(9):4039-4042.
Huerta-Cepas J, Serra F, Bork P. 2016. ETE 3: reconstruction, analysis, and visualization of phylogenomic data. Molecular Biology and Evolution 33(6):1635-1638.
Letunic I, Bork P. 2024. Interactive Tree of Life (iTOL) v6: recent updates to the phylogenetic tree display and annotation tool. Nucleic Acids Research 52(W1):W78-W82.
Felsenstein J. 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39(4):783-791.
Related Skills
tree-io - parsing and preserving BEAST/MrBayes/IQ-TREE annotations so they survive into the figure; the tightest coupling
tree-manipulation - rooting, ladderizing, pruning, and collapsing low-support nodes before drawing
data-visualization/ggplot2-fundamentals - the grammar, themes, and ggsave vector export underlying ggtree
data-visualization/multipanel-figures - composing a tree plus aligned metadata panels into one figure
