The curious case of a Chilean copepod (Tigriopus aff. angulatus) genome assembly

Abstract

Tigriopus copepods are found in splash pools on all seven continents from the equator to Arctic and Antarctic regions. Given their geographic distribution, frequent exposure to extreme environmental conditions, and strong signatures of local adaptation, these copepods have become models for exploring patterns of adaptation to stressful environments. However, most studies focus on a small subset of Tigriopus species, and there are few genome resources representing the diversity of Tigriopus species and populations. Here, we combine long-read, Pacific Biosciences HiFi data with short-read, Illumina HiC and RNA-seq data to assemble and annotate a genome representing a Tigriopus population from the coast of central Chile. Based on the level of divergence that we observed in mitochondrial genes, we also performed a comparison of morphological characteristics between this population and members of the T. angulatus complex. The assembly that we generated (qhTigAngs1.1.pri) includes 12 major scaffolds (N50 19Mbp, L50 7), equivalent to the number of chromosomes in other Tigriopus species. BUSCO and k-mer analyses of the assembly and BUSCO analyses of gene models are relatively complete (89-99%) with respect to gene or k-mer content. The level of divergence that we observed, morphological differences we recorded, and the spatial distance of this population from the type locality suggest that this Chilean population of Tigriopus may represent a novel species that we call Tigriopus aff. angulatus. These genomic resources will help us understand the diversity and structure of Tigriopus species and populations as well as facilitate future comparisons of adaptation across parallel environmental gradients.