Neogastropod (Mollusca, Gastropoda) phylogeny: A step forward with mitogenomes

Phylogenomics of Neogastropoda: the backbone hidden in the bush

The molluscan order Neogastropoda encompasses over 15,000 almost exclusively marine species playing important roles in benthic communities and in the economies of coastal countries. Neogastropoda underwent intensive cladogenesis in early stages of diversification, generating a 'bush' at the base of their evolutionary tree, that has been hard to resolve even with high throughput molecular data. In the present study to resolve the bush, we use a variety of phylogenetic inference methods and a comprehensive exon capture dataset of 1,817 loci (79.6% data occupancy) comprising 112 taxa of 48 out of 60 Neogastropoda families. Our results show consistent topologies and high support in all analyses at (super)family level, supporting monophyly of Muricoidea, Mitroidea, Conoidea, and, with some reservations, Olivoidea and Buccinoidea. Volutoidea and Turbinelloidea as currently circumscribed are clearly paraphyletic. Despite our analyses consistently resolving most backbone nodes, three prove problematic: First, uncertain placement of Cancellariidae, as the sister group to either a Ficoidea-Tonnoidea clade, or to the rest of Neogastropoda, leaves monophyly of Neogastropoda unresolved. Second, relationships are contradictory at the base of the major 'core Neogastropoda' grouping. Third, coalescence-based analyses reject monophyly of the Buccinoidea in relation to Vasidae. We analysed phylogenetic signal of targeted loci in relation to potential biases, and we propose most probable resolutions in the latter two recalcitrant nodes. The uncertain placement of Cancellariidae may be explained by orthology violations due to differential paralog loss shortly after the whole genome duplication, which should be resolved with a curated set of longer loci.

Comparative genomic analysis of chemosensory-related gene families in gastropods

Chemoreception is critical for the survival and reproduction of animals. Except for a reduced group of insects and chelicerates, 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 with a lower number of members. Our phylogenetic analysis of the GPCR gene family across protostomes revealed: (i) remarkable gene family expansions in Gastropoda; (ii) clades including members from all protostomes; and (iii) species-specific clades with a substantial 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.

Collaborative Expression: Transcriptomics of Conus virgo Suggests Contribution of Multiple Secretory Glands to Venom Production

Venomous marine gastropods of the family Conidae are among the most diversified predators in marine realm-in large due to their complex venoms. Besides being a valuable source of bioactive neuropeptides conotoxins, cone-snails venoms are an excellent model for molecular evolution studies, addressing origin of key innovations. However, these studies are handicapped by scarce current knowledge on the tissues involved in venom production, as it is generally assumed the sole prerogative of the venom gland (VG). The role of other secretory glands that are present in all Conus species (salivary gland, SG) or only in some species (accessory salivary gland, ASG) remains poorly understood. Here, for the first time, we carry out a detailed analysis of the VG, SG, and ASG transcriptomes in the vermivorous Conus virgo. We detect multiple transcripts clusters in both the SG and ASG, whose annotations imply venom-related functions. Despite the subsets of transcripts highly-expressed in the VG, SG, and ASG being very distinct, SG expresses an L-, and ASG-Cerm08-, and MEFRR- superfamily conotoxins, all previously considered specific for VG. We corroborate our results with the analysis of published SG and VG transcriptomes from unrelated fish-hunting C. geographus, and C. striatus, possibly fish-hunting C. rolani, and worm-hunting Conus quercinus. In spite of low expression levels of conotoxins, some other specific clusters of putative venom-related peptides are present and may be highly expressed in the SG of these species. Further functional studies are necessary to determine the role that these peptides play in envenomation. In the meantime, our results show importance of routine multi-tissue sampling both for accurate interpretation of tissue-specific venom composition in cone-snails, and for better understanding origin and evolution of venom peptides genes.

Mitogenomic Phylogeny of Tonnoidea Suter, 1913 (1825) (Gastropoda: Caenogastropoda)

The Tonnoidea Suter, 1913 (1825) is a moderately diverse group of large predatory gastropods, the systematics of which remain unclear. In the present study, the complete mitochondrial genomes of nine Tonnoidean species were sequenced. All newly sequenced mitogenomes contain 13 protein-coding genes (PCGs), 22 transfer RNA genes and two ribosomal RNA genes, showing similar patterns in genome size, gene order and nucleotide composition. The ratio of nonsynonymous to synonymous of PCGs indicated that NADH complex genes of Tonnoideans were experiencing a more relaxed purifying selection compared with the COX genes. The reconstructed phylogeny based on the combined amino acid sequences of 13 protein-coding genes and the nucleotide sequences of two rRNA genes supported that Ficidae Meek, 1864 (1840) is a sister to Tonnoidea. The monophylies of all Tonnoidean families were recovered and the internal phylogenetic relationships were consistent with the current classification. The phylogeny also revealed that Tutufa rebuta (Linnaeus, 1758) is composed of at least two different species, indicating that the species diversity within Bursidae Thiele, 1925 might be underestimated. The present study contributes to the understanding of the Tonnoidean systematics, and it could provide important information for the revision of Tonnoidean systematics in the future.

Whole Genome Duplication and Gene Evolution in the Hyperdiverse Venomous Gastropods

The diversity of venomous organisms and the toxins they produce have been increasingly investigated, but taxonomic bias remains important. Neogastropods, a group of marine predators representing almost 22% of the known gastropod diversity, evolved a wide range of feeding strategies, including the production of toxins to subdue their preys. However, whether the diversity of these compounds is at the origin of the hyperdiversification of the group and how genome evolution may correlate with both the compounds and species diversities remain understudied. Among the available gastropods genomes, only eight, with uneven quality assemblies, belong to neogastropods. Here, we generated chromosome-level assemblies of two species belonging to the Tonnoidea and Muricoidea superfamilies (Monoplex corrugatus and Stramonita haemastoma). The two obtained high-quality genomes had 3 and 2.2 Gb, respectively, and 92-89% of the total assembly conformed 35 pseudochromosomes in each species. Through the analysis of syntenic blocks, Hox gene cluster duplication, and synonymous substitutions distribution pattern, we inferred the occurrence of a whole genome duplication event in both genomes. As these species are known to release venom, toxins were annotated in both genomes, but few of them were found in homologous chromosomes. A comparison of the expression of ohnolog genes (using transcriptomes from osphradium and salivary glands in S. haemastoma), where both copies were differentially expressed, showed that most of them had similar expression profiles. The high quality of these genomes makes them valuable reference in their respective taxa, facilitating the identification of genome-level processes at the origin of their evolutionary success.

Chromosome-level genome of the venomous snail Kalloconus canariensis: a valuable model for venomics and comparative genomics

BACKGROUND

Genomes are powerful resources to understand the evolutionary mechanisms underpinning the origin and diversification of the venoms of cone snails (Conidae: Caenogastropoda) and could aid in the development of novel drugs.

FINDINGS

Here, we used PacBio continuous long reads and Omni-C data to assemble the chromosome-level genome of Kalloconus canariensis, a vermivorous cone endemic to the Canary Islands. The final genome size was 2.87 Gb, with a N50 of 79.75 Mb and 91% of the reads located into the 35 largest scaffolds. Up to 55.80% of the genome was annotated as repetitive regions, being class I of transposable elements (16.65%) predominant. The annotation estimated 34,287 gene models. Comparative analysis of this genome with the 2 cone snail genomes released to date (Dendroconus betulinus and Lautoconus ventricosus) revealed similar genome sizes and organization, although chromosome sizes tended to be shorter in K. canariensis. Phylogenetic relationships within subclass Caenogastropoda were recovered with strong statistical support. The family Conidae was recovered as a clade, with K. canariensis plus L. ventricosus sister to D. betulinus.

CONCLUSIONS

Despite the great diversity of cone snails (>900 species) and their venoms (hundreds of peptides per species), only 2 recently reported genomes are available for the group. The high-quality chromosome-level assembly of K. canariensis will be a valuable reference for studying the origin and evolution of conotoxin genes as well as whole-genome duplication events during gastropod evolution.