English: Maximum likelihood tree showing kinorhynch evolutionary relationships. Topology and branch lengths are based on the concatenated 70% occupancy matrix (≥20 taxa, 382 genes, 212,043 amino acids). Identical topology was also found with all 15,976 gene trees summarized with ASTRAL (Supplementary figure 2). Annotated nodes do not have full support at least in one of the analyses. Non-annotated nodes have full support across analyses. Support values are: bootstrap 70% matrix/ASTRAL local posterior probabilities. Outgroups in grayscale, Kinorhynch in-groups in colours. Asterisks mark aberrant species. SEM images represent the species included in each clade except for Pycnophyes sp. and Zelinkaderes yong. A Campyloderes vanhoeffeni, B Centroderes spinosus, C Tubulideres seminoli, D Antygomonas paulae, E Semnoderes armiger, F Sphenoderes poseidon, G Zelinkaderes brightae, H Meristoderes macracanthus, I Echinoderes dujardinii, J Echinoderes ohtsukai, K Echinoderes rex, L Echinoderes kohni, M Pycnophyes ilyocryptus, N Pycnophyes giganteus, O Cristaphyes yushini, P Cateria styx, Q Franciscideres kalenesos, R Dracoderes abei, S Paracentrophyes quadridentatus.