Mitogenomic phylogeny of cone snails endemic to Senegal

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.

Mitogenome Characterization of Four Conus Species and Comparative Analysis

Cone snails, as a type of marine organism, have rich species diversity. Traditionally, classifications of cone snails were based mostly on radula, shell, and anatomical characters. Because of these phenotypic features' high population variability and propensity for local adaptation and convergence, identifying species can be difficult and occasionally inaccurate. In addition, mitochondrial genomes contain high phylogenetic information, so complete mitogenomes have been increasingly employed for inferring molecular phylogeny. To enrich the mitogenomic database of cone snails (Caenogastropoda: Conidae), mitogenomes of four Conus species, i.e., C. imperialis (15,505 bp), C. literatus (15,569 bp), C. virgo (15,594 bp), and C. marmoreus (15,579 bp), were characterized and compared. All 4 of these mitogenomes included 13 protein-coding genes, 2 ribosomal RNA genes, 22 tRNA genes, and non-coding regions. All the Protein Codon Genes (PCGs) of both newly sequenced mitogenomes used TAA or TAG as a terminal codon. Most PCGs used conventional start codon ATG, but an alternative initiation codon GTG was detected in a gene (NADH dehydrogenase subunit 4 (nad4)) of C. imperialis. In addition, the phylogenetic relationships were reconstructed among 20 Conus species on the basis of PCGs, COX1, and the complete mitogenome using both Bayesian Inference (BI) and Maximum Likelihood (ML). The phylogenetic results supported that C. litteratus, C. quercinus, and C. virgo were clustered together as a sister group (PP = 1, BS = 99), but they did not support the phylogenetic relation of C. imperialis and C. tribblei (PP = 0.79, BS = 50). In addition, our study established that PCGs and complete mitogenome are the two useful markers for phylogenetic inference of Conus species. These results enriched the data of the cone snail's mitochondrion in the South China Sea and provided a reliable basis for the interpretation of the phylogenetic relationship of the cone snail based on the mitochondrial genome.

Mitogenomic phylogeny of Cypraeidae (Gastropoda: Mesogastropoda)

Cowries (Family Cypraeidae) are widely distributed in tropical and subtropical seas, with the highest diversity of cowries in the Indo-Pacific region. However, the classification of Cypraeidae, especially at the lower taxonomic levels, is still controversial. In the present study, we determined the complete mitochondrial genomes of 10 cowries. All the newly sequenced mtDNA encode 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and several non-coding regions. The reconstructed mitogenomic phylogeny of Cypraeidae recognized two monophyletic clades, with the first clade comprised of Erroneinae, Cypraeinae and Luriinae and the second clade formed by the single subfamily Erosarinae. The congeneric genetic distance values fall within 0.118–0.144, lower than those above genus level ranging from 0.163 to 0.271, consistent with the current division of genera within Cypraeidae. The divergence time estimated here indicated that the ancestor of Indo-Pacific cowries diversified during the Paleocene, and the closure of the Tethys Seaway might lead to the speciation events of several Indo-Pacific species. This study suggests that the complete mtDNA is a promising tool to improve the phylogenetic resolution of Cypraeidae, and mtDNA could also provide important information for future species delimitation especially within the cowries that possess different morphological phenotypes.

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

The Neogastropoda (Mollusca, Gastropoda) encompass more than 15,000 described species of marine predators, including several model organisms in toxinology, embryology and physiology. However, their phylogenetic relationships remain mostly unresolved and their classification unstable. We took advantage of the many mitogenomes published in GenBank to produce a new molecular phylogeny of the neogastropods. We completed the taxon sampling by using an in‐house bioinformatic pipeline to retrieve mitochondrial genes from 13 transcriptomes, corresponding to five families not represented in GenBank, for a final dataset of 113 taxa. Because mitogenomic data are prone to reconstruction artefacts, eight different evolutionary models were applied to reconstruct phylogenetic trees with IQTREE, RAxML and MrBayes. If the over‐parametrization of some models produced trees with aberrant internal long branches, the global topology of the trees remained stable over models and softwares, and several relationships were revealed or found supported here for the first time. However, even if our dataset encompasses 60% of the valid families of neogastropods, some key taxa are missing and should be added in the future before proposing a revision of the classification of the neogastropods. Our study also demonstrates that even complex models struggle to satisfactorily handle the evolutionary history of mitogenomes, still leading to long‐branch attractions in phylogenetic trees. Other approaches, such as reduced‐genome strategies, must be envisaged to fully resolve the neogastropod phylogeny.

Complete mitochondrial genome of Conus lischkeanus Weinkauff, 1875 (Neogastropoda, Conidae) and phylogenetic implications of the evolutionary diversification of dietary types of Conus species

The family Conidae, commonly known as cone snails, is one of the most intriguing gastropod groups owing to their diverse array of feeding behaviors (diets) and toxin peptides (conotoxins). Conuslischkeanus Weinkauff, 1875 is a worm-hunting species widely distributed from Africa to the Northwest Pacific. In this study, we report the mitochondrial genome sequence of C.lischkeanus and inferred its phylogenetic relationship with other Conus species. Its mitochondrial genome is a circular DNA molecule (16,120 bp in size) composed of 37 genes: 13 protein-coding genes (PCGs), 22 transfer RNA genes, and two ribosomal RNA genes. Phylogenetic analyses of concatenated nucleotide sequences of 13 PCGs and two ribosomal RNA genes showed that C.lischkeanus belongs to the subgenus Lividoconus group, which is grouped with species of the subgenus Virgiconus, and a member of the largest assemblage of worm-hunting (vermivorous) species at the most basal position in this group. Mitochondrial genome phylogeny supports the previous hypothesis that the ancestral diet of cone snails was worm-hunting, and that other dietary types (molluscivous or piscivorous) have secondarily evolved multiple times from different origins. This new, complete mitochondrial genome information provides valuable insights into the mitochondrial genome diversity and molecular phylogeny of Conus species.

The mitochondrial genome and phylogenetic position of a marine snail Nerita (Heminerita) japonica (Gastropoda: Neritimorpha: Neritidae)

The mitochondrial genome of the neritid snail Nerita (Heminerita) japonica (Mollusca: Neritimorpha) from Kumamoto, Japan was determined by whole-genome sequencing. This mitogenome is comprised of 13 protein-coding genes, 2 ribosomal RNA (12S and 16S) genes, and 22 transfer RNA genes, with the same gene order as in the other species of the family Neritidae. A likelihood-based phylogenetic reconstruction recovered the subgenus Heminerita (including N. japonica as its type and N. yoldii from China) as monophyletic and sister to a clade with four species of the subgenera Nerita and Theliostyla.

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.

Complete Mitochondrial Genomes of Two Toxin-Accumulated Nassariids (Neogastropoda: Nassariidae: Nassarius) and Their Implication for Phylogeny

The Indo-Pacific nassariids (genus Nassarius) possesses the highest diversity within the family Nassariidae. However, the previous shell or radula-based classification of Nassarius is quite confusing due to the homoplasy of certain morphological characteristics. The toxin accumulators Nassarius glans and Nassarius siquijorensis are widely distributed in the subtidal regions of the Indo-Pacific Ocean. In spite of their biological significance, the phylogenetic positions of N. glans and N. siquijorensis are still undetermined. In the present study, the complete mitochondrial genomes of N. glans and N. siquijorensis were sequenced. The present mitochondrial genomes were 15,296 and 15,337 bp in length, respectively, showing negative AT skews and positive GC skews as well as a bias of AT rich on the heavy strand. They contained 13 protein coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and several noncoding regions, and their gene order was identical to most caenogastropods. Based on the nucleotide sequences combining 13 protein coding genes and two rRNA genes, a well-supported phylogeny of Nassarius was reconstructed, and several morphological synapomorphies were observed corresponding to the phylogenetic framework. In addition, the sister group relationship between N. variciferus and the remaining toxin-accumulated nassariids was determined, suggesting that the phylogeny might be related to their diet. The divergence time estimation analysis revealed a correlation between speciation events of nassariids and glacial cycles during the Pliocene-Pleistocene epoch.

Comparative Characterization of the Complete Mitochondrial Genomes of the Three Apple Snails (Gastropoda: Ampullariidae) and the Phylogenetic Analyses

The apple snails Pomacea canaliculata, Pomacea diffusa and Pomacea maculate (Gastropoda: Caenogastropoda: Ampullariidae) are invasive pests causing massive economic losses and ecological damage. We sequenced and characterized the complete mitochondrial genomes of these snails to conduct phylogenetic analyses based on comparisons with the mitochondrial protein coding sequences of 47 Caenogastropoda species. The gene arrangements, distribution and content were canonically identical and consistent with typical Mollusca except for the tRNA-Gln absent in P. diffusa. An identifiable control region (d-loop) was absent. Bayesian phylogenetic analysis indicated that all the Ampullariidae species clustered on the same branch. The genus Pomacea clustered together and then with the genus Marisa. The orders Architaenioglossa and Sorbeoconcha clustered together and then with the order Hypsogastropoda. Furthermore, the intergenic and interspecific taxonomic positions were defined. Unexpectedly, Ceraesignum maximum, Dendropoma gregarium, Eualetes tulipa and Thylacodes squamigerus, traditionally classified in order Hypsogastropoda, were isolated from the order Hypsogastropoda in the most external branch of the Bayesian inference tree. The divergence times of the Caenogastropoda indicated that their evolutionary process covered four geological epochs that included the Quaternary, Neogene, Paleogene and Cretaceous periods. This study will facilitate further investigation of species identification to aid in the implementation of effective management and control strategies of these invasive species.

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.