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.pycarries 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:
| Subunit | Complex | H. exemplaris | Also in |
|---|---|---|---|
| SMC2, SMC4 | both | found | Ramazzottius, Paramacrobiotus |
| CAP-D2, CAP-G | condensin I | found | both (CAP-G not found in Ramazzottius) |
| CAP-H | condensin I | found, via the Paramacrobiotus ortholog | both |
| CAP-D3, CAP-G2 | condensin II | found | both |
| CAP-H2 | condensin II | two divergent copies, found only by profile search | two 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.
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, andcaph2_structure_analysis.pyscores 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.
- The reference genome is absent. The H. exemplaris reference, nHd_3.1, isn't in the training data. The only Hypsibius genome is an older assembly from a commercial culture (Sciento), known for bacterial contamination.
- Fragmented genomes lose most of their sequence. Every contig under 10 kb was dropped, leaving 35% of P. richtersi and 14% of I. dastychi.
- Gene-focused data see one tardigrade. The mRNA set holds 19,905 tardigrade sequences, all from P. metropolitanus, the only one with a RefSeq annotation. The genic-window set has 2,481 windows from it and 37 from Ramazzottius, whose annotation lists genes but no mRNAs, and none from the others.
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