Examples gallery ZOMBI2 0.38.0
Click any figure for the code that produces it. Every example simulates with ZOMBI2 and plots with Phylustrator, a separate package — so to run one you need both:
pip install zombi2 phylustrator
Species trees
8Forward birth-death trees. The run keeps the whole history, survivors and extinctions, and the diversification model shows on the tree.
Sp1Yule tree
Pure birth, no extinction — a forward tree of 100 lineages.
birth
Sp2Extinct lineages
The full history behind 50 survivors — their branches solid, extinct lineages dashed.
birth · death
Sp3Rate shifts
Speciation runs slow, then fast, then slow — the burst packs branches between the two dashed regime lines.
birth · changing_at
Sp4Diversity-dependent
Speciation slows as diversity fills up; the skyline rises and plateaus at the cap of 100.
birth · TotalDiversity
Sp5Inherited rates
Each lineage starts at its parent's rate and is nudged from it, so close relatives run at similar speeds and one clade radiates while its sister stays sparse.
birth · Drift
Sp6Independent rates
The same distribution without the inheritance: every lineage draws for itself, so a fast lineage tells you nothing about its neighbours.
birth · LogNormal
Sp7Mass extinction
A pulse at t = 3 culls 75% of lineages — the skyline drops sharply at the dashed line, then recovers.
mass_extinctions
Sp8Shape statistics over many trees
Two thousand trees under each of two processes. zombi2 tools tree --gamma separates them almost perfectly.
Genomes
11Genes on chromosomes. A genome draws as a ring, and two genomes show their synteny. Gene-family events and copy number read against the species tree.
Ge1Gene-family events on the tree
One family's history on the species tree: duplications (squares), losses (crosses), transfers (arrows, donor→recipient).
duplication · transfer · loss
Ge2Profile copy-number
A family × genome copy-number heatmap, its rows locked to the tips. beside(tree, heatmap(profiles)).
Ge3A transfer highway between clades
Transfers steered to run between two clades, by topology rather than by a trait. The barplot counts them by clade pair, so A↔B towers over within-clade.
transfer_to · Clades
Ge4Gene family heterogeneity
Two runs at the same mean rates. Every family alike gives no core at all; letting families differ gives 28 core families and a U-shaped spectrum.
varying_among('families')
Ge5A clade slows down
The shaded clade loses genes twenty times faster after the dashed line, and nothing about it differs before. It ends with about 140 genes against 270 outside — one factor, scoped to a group and to a time.
Clade · changing_at
Ge6Circular genome (ordered)
A genome as a ring — genes evenly spaced by rank, coloured by family, arrows by strand. plot(g, layout="circular") + genes().
Ge7Synteny between two genomes
Two genomes, one per row; ribbons link same-family genes and cross where the order was rearranged. stack([a,b]) + synteny().
Ge8Synteny across a whole clade
Every tip's gene order beside the tree. Genes are coloured by their ancestral position, so each rearrangement is a break in the gradient.
inversion
Ge9An inversion, before → after
One inversion on a circular genome: the segment is reversed and its strands flip. The band marks it in both rings.
inversion · inversion_extent
Ge10Karyotypes
Seven lineages from the three circular chromosomes drawn above the stem. Fission and fusion are marked on the branch where they happened, so the karyotypes are what those marks did.
fission · fusion
Ge11Real genome (Mycoplasma)
A real bacterium: Mycoplasma genitalium, 546 genes at their true base positions. The forward/reverse switch marks the replication origin.
nucleotide · --gffSequences
9The dated tree the sequences evolve down, and an alignment lined up row-for-row with its tips.
Sq1Alignment beside the tree
A single-copy family across 20 species, residues coloured (with a nucleotide key), each row locked to its tip. beside(tree, alignment(aln)).
Sq2A clade's own model
A rate says how fast a lineage evolves; the model says what the change looks like. One clade is pulled toward A and T, and its tips arrive at 0.20 GC against 0.50 outside — each on its own model's equilibrium.
model · Clade
Sq3A protein alignment
The same twenty species evolving amino acids under lg, coloured by chemical class rather than one hue per residue, so a conserved column reads as conserved chemistry. Also poisson, jtt, dayhoff, wag.
Sq4Uncorrelated lognormal clock
Every lineage draws its own rate, with no memory of its parent, so the colour is salt-and-pepper. substitution = PerSite().varying_among('lineages', LogNormal(0.0, 0.55)).
Sq5Uncorrelated gamma clock
The same independent draw with a gamma instead of a lognormal. varying_among('lineages', Gamma(shape=3.31, scale=0.302)).
Sq6Autocorrelated clock
A daughter starts at its parent's rate and is nudged, so the colour moves in clades rather than branch to branch. substitution = PerSite().varying_among('lineages', Drift(LogNormal(0.0, 0.4))).
Sq7Discrete-bin clock
The same inherited drift in steps: the rate takes one of a few values and a daughter moves to a neighbouring one. varying_among('lineages', Drift(LogNormal(0.0, 0.45), bins=6)).
Sq8Ancestral sequences at the nodes
A small tree with its internal nodes numbered, and beside it the sequence at each. The rows are one per node, not aligned to the tips.
substitution · .ancestral
Sq9Indels in a real gene
A real gene down twenty species. The pale band at 158–170 is one insertion, carried by exactly one of the two clades at the root — an indel is a shared derived character.
insertion · deletionTraits
5A trait evolving down the tree. Branches take the colour of its value. Some examples add a companion panel.
Tr1Brownian motion
Free diffusion — sister lineages drift apart with time.
rate
Tr2Ornstein–Uhlenbeck
Pulled to an optimum: a high start (yellow) converges to blue.
rate · pull · reverts_to
Tr3Discrete states
A two-state trait hops between habitats; each branch is painted by its state history.
switch
Tr4Dependent continuous traits
Two traits evolve together (r = 0.9) — two trees, coloured by each trait, and the tip scatter.
rate · correlated
Tr5Dependent discrete traits
Two binary characters where one's flip rate depends on the other's state. X in green, Y in purple, so you can see Y is present where X is; the 2×2 chain is the model.
switch · dependentConditioning
18Two runs, in order. The first run grows the driver on the tree and holds it fixed. The second run reads it. A driver is a trait, a gene family or a whole module. It drives a rate, or which lineage receives a transfer.
Co1Genome reduction
A lifestyle trait drives gene loss. Endosymbionts shed genes faster and gain fewer, so their genomes end up a fifth the size.
trait → loss
Co2Genome expansion
A selection trait drives duplication. Under relaxed selection duplicates accumulate and the genomes grow fivefold.
trait → duplication
Co3HGT uptake by competence
A competence trait drives who receives a transfer rather than a rate. Competent lineages take up DNA more often.
trait → transfer uptake
Co4A continuous driver
A diffusing trait drives gene gain. A Curve turns each value into a factor, so genome size follows the trait.
Co5A saturating curve
The same run with a different Curve. Gene gain rises with the trait and then levels off, so the factor is bounded.
Co6A curve of your own
The response is any function you write. Here it peaks at an intermediate value, which no table of per-state multipliers can express.
continuous trait → origination
Co7A trait drives rate and extent
Hot lineages invert fifteen times as often, and each inversion takes a longer run of genes. Rate and extent are separate targets and multiply.
trait → inversion
Co8A trait sets the substitution rate
The same target reached from the other end of the map. A driver is read wherever it changes, so a lineage that switches halfway down accumulates at two rates.
trait → substitution
Co9A trait driving a trait
A habitat trait sets how fast body size diffuses: twenty times faster where the habitat fluctuates. rate = PerLineage(0.25).scaled_by(habitat, {…}).
Co10One trait drives another
A temperature trait is grown first; body size then diffuses at a rate that reads it. The scale is centred on where it started, so white means it has not moved.
trait → trait
Co11A mobile element spreads its genome
One level conditioning itself. Carrying the element makes a lineage donate twenty-five times more often — a driven transfer is about giving, not receiving.
gene → transfer
Co12Lose the repair gene, evolve faster
A gene family drives the substitution rate. The same tree in time, then in substitutions: branches that have lost the mismatch-repair family run four times longer.
gene → substitution
Co13A gene drives a trait
Carrying a toxin family makes a lineage become pathogenic forty times faster. 74% of the tips with the gene end up pathogenic, against none of those without.
gene → trait
Co14How much of an operon is left
The driver is a fraction, not a yes/no. At the ordered resolution a loss takes a run of neighbours, so an operon goes in blocks and completion drops in steps.
module → substitution
Co15A module, through a step
How much of a four-gene module a lineage keeps decides its metabolism, through a step: lambda f: 20.0 if f > 0.5 else 1.0.
Co16The same driver, another level
The repair operon again, now setting how fast a trait diffuses instead of how fast sequences evolve. One driver reaches every level on the tree.
module → trait rate
Co17One gene’s composition, another’s rate
A sequence drives a sequence. It cannot drive the gene it grows inside — that would read a run’s own output — but it can drive a different one, run after it.
GC → substitution
Co18Composition drives a character
The only pair whose driver is grown at the sequence level and whose target is not. gc() is a number, so the connection carries a Curve.
Joining
4One run makes both. The trait sets the speciation or extinction rate of the lineage carrying it. The trait and the tree therefore come out together.
Jo1A gene drives the splitting
Gene content drives speciation, so the tree is an output. Loss runs at twice origination, yet the family reaches 83% of tips — its carriers split five times as often.
gene content → speciation
Jo2BiSSE
A two-state trait drives speciation — the fast state's clades take over; the inset is the state Markov chain.
trait → speciation
Jo3State-dependent extinction
One state dies far faster; the doomed lineages (dashed) drop out.
trait → extinction
Jo4MuSSE
Three graded speciation rates with constant death — the fastest state fills the tree, extinct lineages dashed.
trait → speciation