Tardigraph
Genetics and genomics

Genomes, genes and what we found

We joined public tardigrade genomes and expression data on one reference, found the condensin subunits a published survey had reported missing, and measured how much tardigrade a genome AI was trained on. Findings here are preliminary until written up; each says what it rests on.

Everything known about what switches what — every gene-regulatory relationship recorded for tardigrades. Eight arrows, two disconnected pieces, and only two of them rest on anyone perturbing anything.

Two developmental datasets on one genome

We put two measurements of gene activity through H. exemplaris development on one genome: Levin et al. 2016, hour by hour in embryos, and Yoshida et al. 2017, across embryo and adult stages. They used different transcript sets, so their results couldn't be compared gene by gene. We aligned Levin's 12,734 assembled transcripts to the nHd_3.1 genome: 12,724 align (median 99.5% identity) and land on 9,069 genes. We joined them with Yoshida's gene-level counts, and lifted everything onto the chromosome-length Hi-C assembly.

The two datasets agree: across 8,541 shared genes, average expression correlates with Spearman ρ = 0.52, against 0.008 when gene labels are shuffled. That gives an hourly-plus-daily developmental expression atlas on one set of gene models.

Joining them also settled Levin's time axis, which the paper states two ways: the series runs about 98 hours, not 49.5 (see Development).

How this is built
Data
Levin et al. 2016's 12,734 assembled transcripts and their per-embryo counts; Yoshida et al. 2017's gene-level counts and the nHd_3.1 genome with its gene models; DNA Zoo's chromosome-length Hi-C assembly of nHd_3.1.
Built with
We aligned the transcripts to nHd_3.1 with minimap2 and joined everything by gene in our own genome database (build_db.py, layer2.py). hic_lift.py carries the coordinates onto the Hi-C chromosomes using DNA Zoo's map of the pieces, with Liftoff as a cross-check.
Checked
Levin against Yoshida, gene by gene: ρ 0.52 across 8,541 genes, 0.008 with the labels shuffled (pass). Before trusting the lift, the script rebuilds every Hi-C scaffold from nHd_3.1 and compares it base for base with DNA Zoo's file.
Limits
Only transcripts that land on one gene enter the comparison. ρ measures agreement in average expression across genes, not stage by stage. The ~98-hour axis is our inference from matching each embryo to Yoshida's daily samples.
More
How the whole project is built · the claims and tests in the graph

Tardigrades aren't missing condensin

We found every condensin subunit in the H. exemplaris proteome, including two divergent copies of CAP-H2, the one BLAST could not find. Condensin complexes fold chromosomes. In a survey of genome folding across the tree of life, Hoencamp et al. 2021 had found only the two core subunits in the tardigrade, SMC2 and SMC4, searching older genome assemblies, and wrote that this "will need to be investigated further".

In the 2017 nHd_3.1 proteome, with reciprocal best hits against human proteins:

SubunitComplexH. exemplarisAlso in
SMC2, SMC4bothfoundRamazzottius, Paramacrobiotus
CAP-D2, CAP-Gcondensin Ifoundboth (CAP-G not found in Ramazzottius)
CAP-Hcondensin Ifound, via the Paramacrobiotus orthologboth
CAP-D3, CAP-G2condensin IIfoundboth
CAP-H2condensin IItwo divergent copies, found only by profile searchtwo in each

That left condensin II without its kleisin, CAP-H2, which BLAST couldn't find. Kleisins change fast, so we built a profile from CAP-H2 in 33 other animals and searched with it. Each tardigrade has two CAP-H2 genes (in H. exemplaris, BV898_16921 and BV898_15594). Each carries both of CAP-H2's signature domains, matches the CAP-H2 profile and not the CAP-H one, and nothing turns up when the proteome is shuffled. The same profile finds human CAP-H2 with overwhelming confidence, and a CAP-H profile finds exactly the known CAP-H in each tardigrade.

Through development, the two CAP-H2 genes rise and fall with the other condensin II subunits: high in eggs in the first two days and in adults, low in between. So condensin II looks complete. Hoencamp et al.'s type I folding in tardigrades may still hold, but not because condensin II is missing. Credit where it's due: van Hooff et al. 2025 found the two Ramazzottius CAP-H2 genes first, in their survey of these complexes across eukaryotes; their sequences are identical to ours. Their survey didn't include H. exemplaris or Paramacrobiotus, and database annotation had already labelled three of the proteins "CAP-H2 C-terminal domain-containing". What we added: the genes in H. exemplaris, the species behind the 2021 question; the answer to that question; and the tests below.

A tree of 68 kleisins places all six tardigrade genes inside CAP-H2, as one tardigrade branch. The two copies come from one duplication: the same copies keep the same neighbouring genes across species, which dates it to before Hypsibius and Paramacrobiotus split, deep in their order. In single embryos, every condensin II subunit is switched on for the first ~46 hours and off after ~47 hours, while condensin I's kleisin stays on. We also redid Hoencamp et al.'s statistical test with the tardigrade moved to "complete". Their link between condensin II and genome folding still holds (p = 0.019, from 0.007), and the tardigrade becomes one of seven species with a complete condensin II and Rabl-like folding. Preprint: Meow-Meow 2026, doi:10.5281/zenodo.22832510 (all versions; version 2 adds the structure models); code and data on GitHub.

What the clasp looks like

A clasp only folds when it grips the ring, so we predicted each tardigrade CAP-H2 end together with its own species' ring protein, using Boltz-2 on a laptop GPU. The C-terminal end of the clasp is known to grip the head of SMC4. Five of the six tardigrade copies are predicted to do the same with confidence (ipTM 0.76–0.81; human 0.83). The controls behave: offered the other ring protein, SMC2, or scrambled into a random sequence of the same letters, the clasps fall to 0.19–0.59. The clasp's other end holds the neck of SMC2. With a short piece of the neck it gave no confident prediction; given the full head and neck, the b copies of H. exemplaris and R. varieornatus grip it in every sample (ipTM 0.62–0.69, human 0.57–0.61), on the same stretch of the neck as the human clasp, while scrambled clasps score 0.11–0.19. So in these two species copy b is predicted to grip both ring proteins, in separately modelled contacts. H. exemplaris a and Paramacrobiotus b still give no signal at this end.

Grey: the SMC4 head and the start of its coiled coil. Coloured: the CAP-H2 C-terminal region, from blue (confident) to red (uncertain), by the model's per-residue confidence (pLDDT). Drag to turn, scroll to zoom.

Human CAP-H20.83H. exemplaris a0.78H. exemplaris b0.79R. varieornatus a0.68R. varieornatus b0.76Paramacrobiotus a0.81Paramacrobiotus b0.81Condensin I clasp (CAP-H) on SMC40.81Wrong partner: clasp on SMC2 (H. ex. a)0.31Wrong partner: clasp on SMC2 (human)0.59Scrambled clasp (H. ex. b)0.19Scrambled clasp (human)0.250.8
Interface confidence (ipTM, 0–1) of the best of three samples for each prediction. Above about 0.8 a predicted interface is usually right; below 0.5 the model is guessing. In the preprint's version 2, with the models and scripts on GitHub.
How this is built
Data
Public protein sequences: each species' two CAP-H2 genes and its own SMC2 and SMC4, with human as the reference. Sequence alignments for each chain come from the ColabFold MMseqs2 server.
Built with
We cut each clasp end (its Pfam domain plus 40 residues either side) and paired it with a head-and-neck piece of its partner, with boundaries taken from human AlphaFold models (caph2_structure_inputs.py). Boltz-2 predicted each pair on a laptop GPU, and caph2_structure_analysis.py scores the interfaces. The viewer draws the C-terminal models with 3Dmol.js.
Checked
C-terminal end on SMC4: pass for five of six tardigrade copies, with wrong-partner and scrambled controls low. N-terminal end on a short SMC2 neck: inconclusive. Retried on the full head and neck: pass for the b copies of H. exemplaris and R. varieornatus, no support for H. exemplaris a or Paramacrobiotus b.
Limits
These are predictions, not measured structures, and ipTM values of 0.6–0.7 are moderate. R. varieornatus a swings from 0.14 to 0.68 across samples. Version 3 of the preprint includes the N-terminal retry.
More
How the whole project is built · the claims and tests in the graph

We also connect this to an imaging result neither paper links it to: H. exemplaris chromosomes stay individualised through interphase in early embryos, each in its own compartment (Papell et al. 2026). Condensin II is what drives that kind of territory in other animals. Hoencamp's Hi-C came from whole animals, so one idea worth testing is that the architecture changes with stage. This is a hypothesis, not a result.

How this is built
Data
The proteomes of H. exemplaris (nHd_3.1, 2017), R. varieornatus and Paramacrobiotus; human condensin proteins; CAP-H2 from 33 other animals; expression from Yoshida et al. 2017 and Levin et al. 2016; Hoencamp et al. 2021's Table S6.
Built with
Reciprocal BLAST against human proteins (condensin_check.py); a profile HMM built from other animals' CAP-H2 (caph2_profile.py); an IQ-TREE tree of 68 kleisins aligned with MAFFT (kleisin_tree.py); shared neighbouring genes found with DIAMOND (caph2_synteny.py); and Hoencamp's Fisher test redone (hoencamp_fisher.py).
Checked
Profile search: pass (human CAP-H2 found, the CAP-H profile finds only CAP-H, a shuffled proteome gives nothing). Kleisin tree: pass. Shared neighbours against 20,000 random gene pairs: pass. Hoencamp's test with the tardigrade moved to "complete": partial (the link holds, p 0.007 becomes 0.019).
Limits
Co-expression supports the pairing but can't prove it, since dividing cells switch on many genes together. R. varieornatus b keeps none of its neighbours, so its place in the duplication rests on the tree. Whether condensin II is active in the animal is not shown, and the embryo link is a hypothesis.
More
Preprint · code and data · how the whole project is built · the graph

How much tardigrade does Evo 2 know?

We checked which tardigrade sequence actually went into Evo 2, a large DNA language model trained on OpenGenome2, a public collection of genomes across life. Its data card lists five tardigrade genomes.

99.6%97.3%74.8%35.2%13.7% P. metropolitanusR. varieornatusH. dujardini (Sciento)P. richtersiI. dastychi 0%50%100%
Share of each listed tardigrade genome present in Evo 2's whole-genome training files. Every contig under 10 kb was dropped, which hollows out the two fragmented assemblies.

So what Evo 2 "knows" about tardigrades is mostly one species, and what counted was annotation format, not biology. Small, fragmented and unannotated genomes of other groups are likely hollowed out the same way.

How this is built
Data
OpenGenome2, Evo 2's public training data from the Arc Institute on Hugging Face: its species list, the whole-genome batches, the mRNA files and the genic-window file. NCBI copies of the five listed tardigrade assemblies.
Built with
We indexed all 90 genome batches and streamed the mRNA and window files for tardigrade records (tar_index.py, stream_filter.sh). Then we matched every record to the NCBI assemblies by its sequence, not its name (compare_ncbi.py), and counted the bases present.
Checked
Every tardigrade interval in the training files is identical to NCBI at its coordinates, so the shares are exact counts, not estimates. Not yet tested against the model itself.
Limits
This is a census of what went in, not a measure of what Evo 2 learned. Whether the 10 kb floor hollows out other groups is argued here, not yet measured.
More
How the whole project is built · the claims in the graph

DNA methylation, from data already public

We're reading DNA methylation straight from public long-read (PacBio HiFi) files, whose methylation calls nobody has looked at, to test a 2022 prediction from gene content that tardigrades lost it. Methylated and unmethylated animals from the same pipeline serve as controls. Work in progress: early results agree with the prediction, and the main catch is that many tardigrade libraries were amplified, which erases methylation.

How this is built
Data
Public PacBio HiFi read files from the Wellcome Sanger Institute's Tree of Life meiofauna project, whose reads carry 5mC calls at CpG sites. Controls from the same project's pipeline: C. elegans (unmethylated), a cockle, a tubeworm and a fish (methylated). The prediction is Engelhardt et al. 2022's.
Built with
A reference-free reader we wrote (hifi_5mc.py) decodes each read's methylation tags, checks that every call sits at a CpG, and summarises how confident the calls are. It adds the reads' CpG depletion as a second, independent signal.
Checked
The reader separates C. elegans from the three methylated controls. Native (unamplified) heterotardigrade runs from the older Sequel II instrument look like the unmethylated control. This is checked by hand, not yet a test in the graph.
Limits
Amplified libraries read as unmethylated whatever the animal. Calls from the newer Revio instrument aren't calibrated yet, no native eutardigrade run has been found, and runs were sampled from their first 25 MB.
More
How the whole project is built · status on the findings page