The draft genome of the microscopic Nemertoderma westbladi sheds light on the evolution of Acoelomorpha genomes

Background: Xenacoelomorpha is a marine clade of microscopic worms that is an important model system for understanding the evolution of key bilaterian novelties, such as the excretory system. Nevertheless, Xenacoelomorpha genomics has been restricted to a few species that either can be cultured in the lab or are centimetres long. Thus far, no genomes are available for Nemertodermatida, one of the group's main clades and whose origin has been dated more than 400 million years ago. Methods: DNA was extracted from a single specimen and sequenced with HiFi following the PacBio Ultra-Low DNA Input protocol. After genome assembly, decontamination, and annotation, the genome quality was benchmarked using two acoel genomes and one Illumina genome as reference. The gene content of three cnidarians, three acoelomorphs, four deuterostomes, and eight protostomes was clustered in orthogroups to make inferences of gene content evolution. Finally, we focused on the genes related to the ultrafiltration excretory system to compare patterns of presence/absence and gene architecture among these clades. Results: We present the first nemertodermatid genome sequenced from a single specimen of Nemertoderma westbladi. Although genome contiguity remains challenging (N50: 60 kb), it is very complete (BUSCO: 80.2%, Metazoa; 88.6%, Eukaryota) and the quality of the annotation allows fine-detail analyses of genome evolution. Acoelomorph genomes seem to be relatively conserved in terms of the percentage of repeats, number of genes, number of exons per gene and intron size. In addition, a high fraction of genes present in both protostomes and deuterostomes are absent in Acoelomorpha. Interestingly, we show that all genes related to the excretory system are present in Xenacoelomorpha except Osr, a key element in the development of these organs and whose acquisition seems to be interconnected with the origin of the specialised excretory system. Conclusion: Overall, these analyses highlight the potential of the Ultra-Low Input DNA protocol and HiFi to generate high-quality genomes from single animals, even for relatively large genomes, making it a feasible option for sequencing challenging taxa, which will be an exciting resource for comparative genomics analyses.

Hidden species diversity and mito-nuclear discordance within the Mediterranean cone snail, Lautoconus ventricosus

The Mediterranean cone snail, Lautoconus ventricosus, is currently considered a single species inhabiting the whole Mediterranean basin and the adjacent Atlantic coasts. Yet, no population genetic study has assessed its taxonomic status. Here, we collected 245 individuals from 75 localities throughout the Mediterranean Sea and used cox1 barcodes, complete mitochondrial genomes, and genome skims to test whether L. ventricosus represents a complex of cryptic species. The maximum likelihood phylogeny based on complete mitochondrial genomes recovered six main clades (hereby named blue, brown, green, orange, red, and violet) with sufficient sequence divergence to be considered putative species. On the other hand, phylogenomic analyses based on 437 nuclear genes only recovered four out of the six clades: blue and orange clades were thoroughly mixed and the brown one was not recovered. This mito-nuclear discordance revealed instances of incomplete lineage sorting and introgression, and may have caused important differences in the dating of main cladogenetic events. Species delimitation tests proposed the existence of at least three species: green, violet, and red+blue+orange (i.e., cyan). Green plus cyan (with sympatric distributions) and violet, had West and East Mediterranean distributions, respectively, mostly separated by the Siculo-Tunisian biogeographical barrier. Morphometric analyses of the shell using species hypotheses as factor and shell length as covariate showed that the discrimination power of the studied parameters was only 70.2%, reinforcing the cryptic nature of the uncovered species, and the importance of integrative taxonomic approaches considering morphology, ecology, biogeography, and mitochondrial and nuclear population genetic variation.

Mitogenomics and hidden-trait models reveal the role of phoresy and host shifts in the diversification of parasitoid blister beetles (Coleoptera: Meloidae)

Changes in life history traits are often considered speciation triggers and can have dramatic effects on the evolutionary history of a lineage. Here, we examine the consequences of changes in two life history traits, host‐type and phoresy, in the hypermetamorphic blister beetles, Meloidae. Subfamilies Nemognathinae and Meloinae exhibit a complex life cycle involving multiple metamorphoses and parasitoidism. Most genera and tribes are bee‐parasitoids, and include phoretic or nonphoretic species, while two tribes feed on grasshopper eggs. These different life strategies are coupled with striking differences in species richness among clades. We generated a mitogenomic phylogeny for Nemognathinae and Meloinae, confirming the monophyly of these two clades, and used the dated phylogeny to explore the association between diversification rates and changes in host specificity and phoresy, using state‐dependent speciation and extinction (SSE) models that include the effect of hidden traits. To account for the low taxon sampling, we implemented a phylogenetic‐taxonomic approach based on birth‐death simulations, and used a Bayesian framework to integrate parameter and phylogenetic uncertainty. Results show that the ancestral hypermetamorphic Meloidae was a nonphoretic bee‐parasitoid, and that transitions towards a phoretic bee‐parasitoid and grasshopper parasitoidism occurred multiple times. Nonphoretic bee‐parasitoid lineages exhibit significantly higher relative extinction and lower diversification rates than phoretic bee‐and grasshopper‐parasitoids, but no significant differences were found between the latter two strategies. This suggests that Orthopteran host shifts and phoresy contributed jointly to the evolutionary success of the parasitoid meloidae. We also demonstrate that SSE models can be used to identify hidden traits coevolving with the focal trait in driving a lineage's diversification dynamics.

A Combined Transcriptomics and Proteomics Approach Reveals the Differences in the Predatory and Defensive Venoms of the Molluscivorous Cone Snail Cylinder ammiralis (Caenogastropoda: Conidae)

Venoms are complex mixtures of proteins that have evolved repeatedly in the animal kingdom. Cone snail venoms represent one of the best studied venom systems. In nature, this venom can be dynamically adjusted depending on its final purpose, whether to deter predators or hunt prey. Here, the transcriptome of the venom gland and the proteomes of the predation-evoked and defensive venoms of the molluscivorous cone snail Cylinder ammiralis were catalogued. A total of 242 venom-related transcripts were annotated. The conotoxin superfamilies presenting more different peptides were O1, O2, T, and M, which also showed high expression levels (except T). The three precursors of the J superfamily were also highly expressed. The predation-evoked and defensive venoms showed a markedly distinct profile. A total of 217 different peptides were identified, with half of them being unique to one venom. A total of 59 peptides ascribed to 23 different protein families were found to be exclusive to the predatory venom, including the cono-insulin, which was, for the first time, identified in an injected venom. A total of 43 peptides from 20 protein families were exclusive to the defensive venom. Finally, comparisons of the relative abundance (in terms of number of peptides) of the different conotoxin precursor superfamilies showed that most of them present similar abundance regardless of the diet.

The genome of the venomous snail Lautoconus ventricosus sheds light on the origin of conotoxin diversity

BACKGROUND

Venoms are deadly weapons to subdue prey or deter predators that have evolved independently in many animal lineages. The genomes of venomous animals are essential to understand the evolutionary mechanisms involved in the origin and diversification of venoms.

RESULTS

Here, we report the chromosome-level genome of the venomous Mediterranean cone snail, Lautoconus ventricosus (Caenogastropoda: Conidae). The total size of the assembly is 3.59 Gb; it has high contiguity (N50 = 93.53 Mb) and 86.6 Mb of the genome assembled into the 35 largest scaffolds or pseudochromosomes. On the basis of venom gland transcriptomes, we annotated 262 complete genes encoding conotoxin precursors, hormones, and other venom-related proteins. These genes were scattered in the different pseudochromosomes and located within repetitive regions. The genes encoding conotoxin precursors were normally structured into 3 exons, which did not necessarily coincide with the 3 structural domains of the corresponding proteins. Additionally, we found evidence in the L. ventricosus genome for a past whole-genome duplication event by means of conserved gene synteny with the Pomacea canaliculata genome, the only one available at the chromosome level within Caenogastropoda. The whole-genome duplication event was further confirmed by the presence of a duplicated hox gene cluster. Key genes for gastropod biology including those encoding proteins related to development, shell formation, and sex were located in the genome.

CONCLUSIONS

The new high-quality L. ventricosus genome should become a reference for assembling and analyzing new gastropod genomes and will contribute to future evolutionary genomic studies among venomous animals.

Comparative transcriptomics of the venoms of continental and insular radiations of West African cones

The transcriptomes of the venom glands of 13 closely related species of vermivorous cones endemic to West Africa from genera Africonus and Varioconus were sequenced and venom repertoires compared within a phylogenetic framework using one Kalloconus species as outgroup. The total number of conotoxin precursors per species varied between 108 and 221. Individuals of the same species shared about one-fourth of the total conotoxin precursors. The number of common sequences was drastically reduced in the pairwise comparisons between closely related species, and the phylogenetical signal was totally eroded at the inter-generic level (no sequence was identified as shared derived), due to the intrinsic high variability of these secreted peptides. A common set of four conotoxin precursor superfamilies (T, O1, O2 and M) was expanded in all studied cone species, and thus, they are considered the basic venom toolkit for hunting and defense in the West African vermivorous cone snails. Maximum-likelihood ancestral character reconstructions inferred shared conotoxin precursors preferentially at internal nodes close to the tips of the phylogeny (between individuals and between closely related species) as well as in the common ancestor of Varioconus. Besides the common toolkit, the two genera showed significantly distinct catalogues of conotoxin precursors in terms of type of superfamilies present and the abundance of members per superfamily, but had similar relative expression levels indicating functional convergence. Differential expression comparisons between vermivorous and piscivorous cones highlighted the importance of the A and S superfamilies for fish hunting and defense.

Conotoxin Diversity in Chelyconus ermineus (Born, 1778) and the Convergent Origin of Piscivory in the Atlantic and Indo-Pacific Cones

The transcriptome of the venom duct of the Atlantic piscivorous cone species Chelyconus ermineus (Born, 1778) was determined. The venom repertoire of this species includes at least 378 conotoxin precursors, which could be ascribed to 33 known and 22 new (unassigned) protein superfamilies, respectively. Most abundant superfamilies were T, W, O1, M, O2, and Z, accounting for 57% of all detected diversity. A total of three individuals were sequenced showing considerable intraspecific variation: each individual had many exclusive conotoxin precursors, and only 20% of all inferred mature peptides were common to all individuals. Three different regions (distal, medium, and proximal with respect to the venom bulb) of the venom duct were analyzed independently. Diversity (in terms of number of distinct members) of conotoxin precursor superfamilies increased toward the distal region whereas transcripts detected toward the proximal region showed higher expression levels. Only the superfamilies A and I3 showed statistically significant differential expression across regions of the venom duct. Sequences belonging to the alpha (motor cabal) and kappa (lightning-strike cabal) subfamilies of the superfamily A were mainly detected in the proximal region of the venom duct. The mature peptides of the alpha subfamily had the α4/4 cysteine spacing pattern, which has been shown to selectively target muscle nicotinic-acetylcholine receptors, ultimately producing paralysis. This function is performed by mature peptides having a α3/5 cysteine spacing pattern in piscivorous cone species from the Indo-Pacific region, thereby supporting a convergent evolution of piscivory in cones.

Phylogenetic relationships of cone snails endemic to Cabo Verde based on mitochondrial genomes

BACKGROUND

Due to their great species and ecological diversity as well as their capacity to produce hundreds of different toxins, cone snails are of interest to evolutionary biologists, pharmacologists and amateur naturalists alike. Taxonomic identification of cone snails still relies mostly on the shape, color, and banding patterns of the shell. However, these phenotypic traits are prone to homoplasy. Therefore, the consistent use of genetic data for species delimitation and phylogenetic inference in this apparently hyperdiverse group is largely wanting. Here, we reconstruct the phylogeny of the cones endemic to Cabo Verde archipelago, a well-known radiation of the group, using mitochondrial (mt) genomes.

RESULTS

The reconstructed phylogeny grouped the analyzed species into two main clades, one including Kalloconus from West Africa sister to Trovaoconus from Cabo Verde and the other with a paraphyletic Lautoconus due to the sister group relationship of Africonus from Cabo Verde and Lautoconus ventricosus from Mediterranean Sea and neighboring Atlantic Ocean to the exclusion of Lautoconus endemic to Senegal (plus Lautoconus guanche from Mauritania, Morocco, and Canary Islands). Within Trovaoconus, up to three main lineages could be distinguished. The clade of Africonus included four main lineages (named I to IV), each further subdivided into two monophyletic groups. The reconstructed phylogeny allowed inferring the evolution of the radula in the studied lineages as well as biogeographic patterns. The number of cone species endemic to Cabo Verde was revised under the light of sequence divergence data and the inferred phylogenetic relationships.

CONCLUSIONS

The sequence divergence between continental members of the genus Kalloconus and island endemics ascribed to the genus Trovaoconus is low, prompting for synonymization of the latter. The genus Lautoconus is paraphyletic. Lautoconus ventricosus is the closest living sister group of genus Africonus. Diversification of Africonus was in allopatry due to the direct development nature of their larvae and mainly triggered by eustatic sea level changes during the Miocene-Pliocene. Our study confirms the diversity of cone endemic to Cabo Verde but significantly reduces the number of valid species. Applying a sequence divergence threshold, the number of valid species within the sampled Africonus is reduced to half.

Mitogenomic phylogeny of cone snails endemic to Senegal

Cone snails attain in Senegal one of their highest peaks of species diversity throughout the continental coast of Western Africa. A total of 15 endemic species have been described, all placed in the genus Lautoconus. While there is ample data regarding the morphology of the shell and the radular tooth of these species, virtually nothing is known regarding the genetic diversity and phylogenetic relationships of one of the most endangered groups of cones. In this work, we determined the complete or near-complete (only lacking the control region) mitochondrial (mt) genomes of 17 specimens representing 11 endemic species (Lautoconus belairensis, Lautoconus bruguieresi, Lautoconus cacao, Lautoconus cloveri, Lautoconus cf. echinophilus, Lautoconus guinaicus, Lautoconus hybridus, Lautoconus senegalensis, Lautoconus mercator, Lautoconus taslei, and Lautoconus unifasciatus). We also sequenced the complete mt genome of Lautoconus guanche from the Canary Islands, which has been related to the cones endemic to Senegal. All mt genomes share the same gene arrangement, which conforms to the consensus reported for Conidae, Neogastropoda and Caenogastropoda. Phylogenetic analyses using probabilistic methods recovered three major lineages, whose divergence coincided in time with sea level and ocean current changes as well as temperature fluctuations during the Messinian salinity crisis and the Plio-Pleistocene transition. Furthermore, the three lineages corresponded to distinct types of radular tooth (robust, small, and elongated), suggesting that dietary specialization could be an additional evolutionary driver in the diversification of the cones endemic to Senegal. The reconstructed phylogeny showed several cases of phenotypic convergence (cryptic species) and questions the validity of some species (ecotypes or phenotypic plasticity), both results having important taxonomic and conservation consequences.