Analyze, manipulate, compare, annotate, and visualize phylogenetic or other hierarchical trees with ETE 4. Use for Newick/Nexus tree I/O, topology edits and pattern matching, Robinson-Foulds comparisons, gene-tree evolutionary events and reconciliation, NCBI/GTDB taxonomy, SmartView exploration, and publication rendering. Do not use it to infer trees from raw sequences; align sequences and infer a tree first.
npx skills add https://github.com/K-Dense-AI/scientific-agent-skills --skill etetoolkit
Use ETE 4 to work with an existing tree:
Newick trees
TreePatternPhyloTreeETE does not replace sequence alignment or phylogenetic inference software. For
raw sequences, first use MAFFT or another aligner and IQ-TREE 2, FastTree, or
another inference tool; then load the resulting tree into ETE.
This skill targets ETE 4.4.0, released September 3, 2025 and verified as the
current PyPI release on July 23, 2026.
Use https://etetoolkit.github.io/ete/ for ETE 4 documentation. The
etetoolkit.org/docs/latest pages are legacy ETE 3 documentation despite the
URL name.
Do not silently translate these examples back to ETE 3:
ete4, not ete3rely on path-string heuristics retained in ETE 4.4.0
parser=, not format=props, add_prop(), and add_props()leaves(), descendants(), and related methods return iteratorsnode.is_leaf and node.is_root are properties, not methodstree["name"], not tree & "name"For porting older code, load
references/migration-ete3-to-ete4.md.
Install the pinned base package:
uv pip install "ete4==4.4.0"
Add only the visualization extra required by the workflow:
# SmartView static PNG screenshots
uv pip install "ete4[render-sm]==4.4.0"
# Legacy Qt renderer for PNG, PDF, and SVG
uv pip install "ete4[treeview]==4.4.0"
Confirm the active environment:
uv run --with "ete4==4.4.0" python -c "import ete4; print(ete4.__version__)"
No credentials are required. NCBI and GTDB workflows download public taxonomy
data and can consume substantial disk space; see
references/taxonomy.md before the first update.
from pathlib import Path
from ete4 import Tree
# Use an open file object for files; reserve strings for Newick text.
with Path("tree.nw").open(encoding="utf-8") as handle:
tree = Tree(handle, parser=1) # parser 1: internal node names
print(tree.to_str(props=["name", "dist"], compact=True))
print("Leaves:", list(tree.leaf_names()))
# Search and annotate.
focal = tree["species1"]
focal.add_props(host="human", status="focal")
# Keep selected tips while preserving pairwise branch-length distances.
tree.prune(
["species1", "species2", "species3"],
preserve_branch_length=True,
)
# Root and serialize explicitly.
tree.set_midpoint_outgroup()
tree.write(
outfile="processed.nw",
parser=1,
props=["host", "status"],
)
Choose the parser deliberately. A parser mismatch is the most common cause of
NewickError, lost internal labels, or support values being read as names.
See references/api_reference.md.
from ete4 import Tree
tree = Tree("((A:1,B:1)CladeAB:0.4,C:2)Root;", parser=1)
for node in tree.traverse("preorder"):
label = node.name if node.name is not None else node.id
print(label, node.level, node.is_leaf, node.dist)
tree["A"].add_prop("group", "case")
tree["B"].add_prop("group", "control")
mrca = tree.common_ancestor("A", "B")
print(mrca.name)
tree.write(
outfile="annotated.nhx",
parser=1,
props=["group"],
format_root_node=True,
)
Node names need not be unique. tree["A"] returns the first match; use
list(tree.search_nodes(name="A")) and validate the count when duplicates are
possible.
from ete4 import Tree
tree_a = Tree("((A,B),(C,D));")
tree_b = Tree("((A,C),(B,D));")
(
rf,
max_rf,
common_leaves,
edges_a,
edges_b,
discarded_a,
discarded_b,
) = tree_a.robinson_foulds(tree_b)
normalized_rf = rf / max_rf if max_rf else 0.0
print(rf, max_rf, normalized_rf, sorted(common_leaves))
RF comparison uses shared leaf labels and requires meaningful, preferably
unique names. Decide explicitly whether rooted or unrooted comparison is
scientifically appropriate.
from ete4 import PhyloTree
gene_tree = PhyloTree(
"((Hsa|g1,Ptr|g1),(Hsa|g2,Mmu|g1));",
sp_naming_function=lambda name: name.split("|", 1)[0],
)
for event in gene_tree.get_descendant_evol_events(sos_thr=0.0):
relationship = "speciation/orthology" if event.etype == "S" else "duplication/paralogy"
print(relationship, sorted(event.in_seqs), sorted(event.out_seqs))
Species-overlap calls are inferences from the supplied topology and naming
function, not independent evidence of orthology. Pass the naming function
explicitly, and use a rooted, fully bifurcating gene tree. For strict
reconciliation, use a curated species tree and
gene_tree.reconcile(species_tree).
from ete4 import NCBITaxa
ncbi = NCBITaxa()
names = ["Homo sapiens", "Pan troglodytes", "Mus musculus"]
name_to_taxids = ncbi.get_name_translator(names)
missing = [name for name in names if name not in name_to_taxids]
if missing:
raise ValueError(f"Names not resolved by NCBI taxonomy: {missing}")
taxids = [name_to_taxids[name][0] for name in names]
taxonomy_tree = ncbi.get_topology(taxids)
print(taxonomy_tree.to_str(props=["sci_name", "rank"]))
ETE 4 also provides GTDBTaxa for genome-centric bacterial and archaeal
taxonomy. Do not mix NCBI numeric TaxIDs and GTDB string identifiers.
Interactive SmartView:
from ete4 import Tree
tree = Tree("((A:1,B:1)90:0.2,C:1);", parser="support")
tree.explore()
Static SmartView screenshot:
tree.render_sm("tree.png", w=1200, h=800)
render_sm() produces PNG screenshot data; use the Qt treeview renderer when
the deliverable must be vector PDF or SVG. Load
references/visualization.md for layouts,
faces, remote exploration, and renderer selection.
Run from this skill directory. The commands below use a pinned, isolated ETE 4
runtime through uv run --with.
uv run --with "ete4==4.4.0" python scripts/tree_operations.py \
stats tree.nw --parser 1
uv run --with "ete4==4.4.0" python scripts/tree_operations.py \
ascii tree.nw --parser 1 --props name,dist
uv run --with "ete4==4.4.0" python scripts/tree_operations.py \
convert tree.nw output.nw \
--input-parser 1 --output-parser 1
uv run --with "ete4==4.4.0" python scripts/tree_operations.py \
reroot tree.nw rooted.nw \
--parser 1 --midpoint
uv run --with "ete4==4.4.0" python scripts/tree_operations.py \
prune tree.nw pruned.nw \
--parser 1 --keep species1 species2 species3
uv run --with "ete4==4.4.0" python scripts/tree_operations.py \
compare tree_a.nw tree_b.nw
Use --keep-file taxa.txt instead of --keep ... for one taxon per line.
The script refuses ambiguous or missing requested names rather than silently
producing a partial tree.
# Interactive SmartView
uv run --with "ete4==4.4.0" python scripts/quick_visualize.py \
tree.nw --parser 1
# SmartView PNG (requires ete4[render-sm])
uv run --with "ete4[render-sm]==4.4.0" python scripts/quick_visualize.py \
tree.nw tree.png \
--parser support --mode circular --show-support --color-by-support
# Vector output via Qt treeview (requires ete4[treeview])
uv run --with "ete4[treeview]==4.4.0" python scripts/quick_visualize.py \
tree.nw tree.svg \
--parser 1 --engine treeview --title "Species phylogeny"
Before reporting a result:
branch lengths.
comparison.
should remain unchanged.
not evolutionary evidence.
database snapshot in reproducible analyses.
get_cached_content() for repeateddescendant-content queries.
Load only the reference needed for the task:
references/api_reference.md — ETE 4 coreclasses, parsers, properties, traversal, I/O, topology, and comparison
references/workflows.md — complete analysispatterns, validation, reconciliation, batching, and large-tree work
references/visualization.md — SmartView,layouts/faces, PNG screenshots, and Qt vector rendering
references/taxonomy.md — NCBI and GTDB setup,translation, topology, annotation, and reproducibility
references/migration-ete3-to-ete4.md— breaking API changes and porting checklist
Integration with protocols.io API for managing scientific protocols. This skill should be used when working with protocols.io to search, create, update, or publish protocols; manage protocol steps and materials; handle discussions and comments; organize workspaces; upload and manage files; or integrate protocols.io functionality into workflows. Applicable for protocol discovery, collaborative protocol development, experiment tracking, lab protocol management, and scientific documentation.
Analyzes job descriptions and generates tailored resumes that highlight relevant experience, skills, and achievements to maximize interview chances
Generate Excalidraw diagrams from natural language descriptions. Use when asked to "create a diagram", "make a flowchart", "visualize a process", "draw a system architecture", "create a mind map", or "generate an Excalidraw file". Supports flowcharts, relationship diagrams, mind maps, and system architecture diagrams. Outputs .excalidraw JSON files that can be opened directly in Excalidraw.
Build and distribute Expo development clients locally or via TestFlight
Use when you have a written implementation plan to execute in a separate session with review checkpoints
Data structure for annotated matrices in single-cell analysis. Use when working with .h5ad files or integrating with the scverse ecosystem. This is the data format skill—for analysis workflows use scanpy; for probabilistic models use scvi-tools; for population-scale queries use cellxgene-census.
Benchling R&D platform integration. Access registry (DNA, proteins), inventory, ELN entries, workflows via API, build Benchling Apps, query Data Warehouse, for lab data management automation.
Comprehensive molecular biology toolkit. Use for sequence manipulation, file parsing (FASTA/GenBank/PDB), phylogenetics, and programmatic NCBI/PubMed access (Bio.Entrez). Best for batch processing, custom bioinformatics pipelines, BLAST automation. For quick lookups use gget; for multi-service integration use bioservices.
Take k-dense-ai/etetoolkit from the repository into ~/.claude/skills for personal
use, or into .claude/skills inside a project.
The agent identifies a skill by the name field in its header. Two skills with the
same name cannot sit side by side — one of them will be ignored.
The instructions reference pip, uv.
Without those the skill loads but fails at the first command.