The spider genus Dysdera has undergone a remarkable diversification in the oceanic archipelago of the Canary Islands, ~60 endemic species originated during the 20 million years since the origin of the archipelago. This evolutionary radiation has been accompanied by substantial dietary shifts, often characterized by phenotypic modifications encompassing morphological, metabolic and behavioral changes. Hence, these endemic spiders represent an excellent model for understanding the evolutionary drivers and to pinpoint the genomic determinants underlying adaptive radiations. Recently, we achieved the first chromosome-level genome assembly of one of the endemic species, D. silvatica, providing a high-quality reference sequence for evolutionary genomics studies. Here, we conducted a low-coverage based resequencing study of a natural population of D. silvatica from La Gomera island. Taking advantage of the new high-quality genome, we characterized genome-wide levels of nucleotide polymorphism, divergence, and linkage disequilibrium, and inferred the demographic history of this population. We also performed comprehensive genome-wide scans for recent positive selection. Our findings uncovered exceptionally high levels of nucleotide diversity and recombination in this geographically restricted endemic species, indicative of large historical effective population sizes. Furthermore, we identified genomic regions potentially under positive selection, shedding light on relevant biological processes, such as vision and nitrogen extraction as possible targets of adaptation and eventually, as drivers of the species diversification. This pioneering study in spiders endemic of an oceanic archipelago lays the groundwork for broader population genomics investigations aimed at understanding the genetic mechanisms driven adaptive radiations in island ecosystems.

Chemoreception is critical for the survival and reproduction of animals. Except for a reduced group of insects and spiders, the molecular identity of chemosensory proteins is poorly understood in invertebrates. Gastropoda is the extant mollusk class with the greatest species richness, including marine, freshwater, and terrestrial lineages, and likely, highly diverse chemoreception systems. Here, we performed a comprehensive comparative genome analysis taking advantage of the chromosome-level information of two Gastropoda species, one of which belongs to a lineage that underwent a whole genome duplication event. We identified thousands of previously uncharacterized chemosensory-related genes, the majority of them encoding G protein-coupled receptors (GPCR), mostly organized into clusters distributed across all chromosomes. We also detected gene families encoding degenerin epithelial sodium channels (DEG-ENaC), ionotropic receptors (IR), sensory neuron membrane proteins (SNMP), Niemann–Pick type C2 (NPC2) proteins, and lipocalins, although much smaller in size. Our phylogenetic analysis of the GPCR gene family across protostomes revealed: (i) large gene family expansions in Gastropoda; (ii) clades including members from all protostomes; and (iii) species-specific clades with a huge number of receptors. For the first time, we provide new and valuable knowledge into the evolution of the chemosensory gene families in invertebrates other than arthropods.

We present the chromosome-level genome assembly of Dysdera silvatica Schmidt, 1981, a nocturnal ground-dwelling spider endemic from the Canary Islands. The genus Dysdera has undergone a remarkable diversification in this archipelago mostly associated with shifts in the level of trophic specialization, becoming an excellent model to study the genomic drivers of adaptive radiations. The new assembly (1.37 Gb; and scaffold N50 of 174.2 Mb), was performed using the chromosome conformation capture scaffolding technique, represents a continuity improvement of more than 4,500 times with respect to the previous version. The seven largest scaffolds or pseudochromosomes cover 87% of the total assembly size and match consistently with the seven chromosomes of the karyotype of this species, including the characteristic large X chromosome. To illustrate the value of this new resource we performed a comprehensive analysis of the two major arthropod chemoreceptor gene families (i.e., gustatory and ionotropic receptors). We identified 545 chemoreceptor sequences distributed across all pseudochromosomes, with a notable underrepresentation in the X chromosome. At least 54% of them localize in 83 genomic clusters with a significantly lower evolutionary distances between them than the average of the family, suggesting a recent origin of many of them. This chromosome-level assembly is the first high-quality genome representative of the Synspermiata clade, and just the third among spiders, representing a new valuable resource to gain insights into the structure and organization of chelicerate genomes, including the role that structural variants, repetitive elements and large gene families played in the extraordinary biology of spiders.

Gene annotation is a critical bottleneck in genomic research, especially for the comprehensive study of very large gene families in the genomes of non-model organisms. Despite the recent progress in automatic methods, state-of-the-art tools used for this task often produce inaccurate annotations, such as fused, chimeric, partial or even completely absent gene models for many family copies, errors that require considerable extra efforts to be corrected. Here we present BITACORA, a bioinformatics solution that integrates popular sequence similarity-based search tools and Perl scripts to facilitate both the curation of these inaccurate annotations and the identification of previously undetected gene family copies directly in genomic DNA sequences. We tested the performance of BITACORA in annotating the members of two chemosensory gene families with different repertoire size in seven available genome sequences, and compared its performance with that of Augustus-PPX, a tool also designed to improve automatic annotations using a sequence similarity-based approach. Despite the relatively high fragmentation of some of these drafts, BITACORA was able to improve the annotation of many members of these families and detected thousands of new chemoreceptors encoded in genome sequences. The program creates general feature format (GFF) files, with both curated and newly identified gene models, and FASTA files with the predicted proteins. These outputs can be easily integrated in genomic annotation editors, greatly facilitating subsequent manual annotation and downstream evolutionary analyses.