Mainly because of the problems with species delineation, the biogeography of microscopic organisms is notoriously difficult to elucidate. In this contribution, variable nuclear and mitochondrial DNA markers were sequenced from individual specimens representing the Echiniscus virginicus complex that are morphologically indistinguishable under light microscopy (five populations from the temperate Eastern Nearctic and 13 populations from the subtropical and tropical zone). A range of methods was used to dissect components of variability within the complex (Bayesian inference, haplotype networks, Poisson tree processes, automatic barcode gap discovery delineations, principal components analysis and ANOVA). We found deep divergence between the temperate Eastern Nearctic E. virginicus and pantropical Echiniscus lineatus in all three genetic markers. In contrast, intraspecific genetic variation was very low, regardless of the geographical distance between the populations. Moreover, for the first time, statistical predictions of tardigrade geographical distributions were modelled. The factor determining the allopatric geographical ranges of deceptively similar species analysed in this study is most likely to be the type of climate. Our study shows that widespread tardigrade species exist, and both geographical distribution modelling and the genetic structure of populations of the pantropical E. lineatus suggest wind-mediated (aeolian) passive long-distance dispersal.

Homeotherms and many poikilotherms display a positive relationship between body size and latitude, but this has rarely been investigated in microscopic animals. We analysed all published records of marine Tardigrada to address whether microscopic marine invertebrates have similar ecogeographical patterns to macroscopic animals. The data were analysed using spatially explicit generalized least squares models and linear models. We looked for latitudinal patterns in body size and species richness, testing for sampling bias and phylogenetic constraints. No latitudinal pattern was detected for species richness, and sampling bias was the strongest correlate of species richness. A hump-shaped increase in median body size with latitude was found, and the effect remained significant for the Northern Hemisphere but not for the Southern. The most significant effect supporting the latitudinal gradient was on minimum body size, with smaller species disappearing at higher latitudes. Our results suggest that biogeographical signals were observed for body size, albeit difficult to detect in poorly studied groups because of swamping from biased sampling effort and from low sample size. We did not find a significant correlation with the latitudinal pattern of body size and ecologically relevant net primary productivity.