Lack of Signal for the Impact of Conotoxin Gene Diversity on Speciation Rates in Cone Snails.

Understanding why some groups of organisms are more diverse than others is a central goal in macroevolution. Evolvability, or the intrinsic capacity of lineages for evolutionary change, is thought to influence disparities in species diversity across taxa. Over macroevolutionary time scales, clades that exhibit high evolvability are expected to have higher speciation rates. Cone snails (family: Conidae, $>$900 spp.) provide a unique opportunity to test this prediction because their toxin genes can be used to characterize differences in evolvability between clades. Cone snails are carnivorous, use prey-specific venom (conotoxins) to capture prey, and the genes that encode venom are known and diversify through gene duplication. Theory predicts that higher gene diversity confers a greater potential to generate novel phenotypes for specialization and adaptation. Therefore, if conotoxin gene diversity gives rise to varying levels of evolvability, conotoxin gene diversity should be coupled with macroevolutionary speciation rates. We applied exon capture techniques to recover phylogenetic markers and conotoxin loci across 314 species, the largest venom discovery effort in a single study. We paired a reconstructed timetree using 12 fossil calibrations with species-specific estimates of conotoxin gene diversity and used trait-dependent diversification methods to test the impact of evolvability on diversification patterns. Surprisingly, we did not detect any signal for the relationship between conotoxin gene diversity and speciation rates, suggesting that venom evolution may not be the rate-limiting factor controlling diversification dynamics in Conidae. Comparative analyses showed some signal for the impact of diet and larval dispersal strategy on diversification patterns, though detection of a signal depended on the dataset and the method. If our results remain true with increased taxonomic sampling in future studies, they suggest that the rapid evolution of conid venom may cause other factors to become more critical to diversification, such as ecological opportunity or traits that promote isolation among lineages.

Targeted Sequencing of Venom Genes from Cone Snail Genomes Improves Understanding of Conotoxin Molecular Evolution.

To expand our capacity to discover venom sequences from the genomes of venomous organisms, we applied targeted sequencing techniques to selectively recover venom gene superfamilies and nontoxin loci from the genomes of 32 cone snail species (family, Conidae), a diverse group of marine gastropods that capture their prey using a cocktail of neurotoxic peptides (conotoxins). We were able to successfully recover conotoxin gene superfamilies across all species with high confidence (> 100× coverage) and used these data to provide new insights into conotoxin evolution. First, we found that conotoxin gene superfamilies are composed of one to six exons and are typically short in length (mean = ∼85 bp). Second, we expanded our understanding of the following genetic features of conotoxin evolution: 1) positive selection, where exons coding the mature toxin region were often three times more divergent than their adjacent noncoding regions, 2) expression regulation, with comparisons to transcriptome data showing that cone snails only express a fraction of the genes available in their genome (24-63%), and 3) extensive gene turnover, where Conidae species varied from 120 to 859 conotoxin gene copies. Finally, using comparative phylogenetic methods, we found that while diet specificity did not predict patterns of conotoxin evolution, dietary breadth was positively correlated with total conotoxin gene diversity. Overall, the targeted sequencing technique demonstrated here has the potential to radically increase the pace at which venom gene families are sequenced and studied, reshaping our ability to understand the impact of genetic changes on ecologically relevant phenotypes and subsequent diversification.

Exon-Capture-Based Phylogeny and Diversification of the Venomous Gastropods (Neogastropoda, Conoidea).

Transcriptome-based exon capture methods provide an approach to recover several hundred markers from genomic DNA, allowing for robust phylogenetic estimation at deep timescales. We applied this method to a highly diverse group of venomous marine snails, Conoidea, for which published phylogenetic trees remain mostly unresolved for the deeper nodes. We targeted 850 protein coding genes (678,322 bp) in ca. 120 samples, spanning all (except one) known families of Conoidea and a broad selection of non-Conoidea neogastropods. The capture was successful for most samples, although capture efficiency decreased when DNA libraries were of insufficient quality and/or quantity (dried samples or low starting DNA concentration) and when targeting the most divergent lineages. An average of 75.4% of proteins was recovered, and the resulting tree, reconstructed using both supermatrix (IQ-tree) and supertree (Astral-II, combined with the Weighted Statistical Binning method) approaches, are almost fully supported. A reconstructed fossil-calibrated tree dates the origin of Conoidea to the Lower Cretaceous. We provide descriptions for two new families. The phylogeny revealed in this study provides a robust framework to reinterpret changes in Conoidea anatomy through time. Finally, we used the phylogeny to test the impact of the venom gland and radular type on diversification rates. Our analyses revealed that repeated losses of the venom gland had no effect on diversification rates, while families with a breadth of radula types showed increases in diversification rates, thus suggesting that trophic ecology may have an impact on the evolution of Conoidea.

Range instability leads to cytonuclear discordance in a morphologically cryptic ground squirrel species complex.

The processes responsible for cytonuclear discordance frequently remain unclear. Here, we employed an exon capture data set and demographic methods to test hypotheses generated by species distribution models to examine how contrasting histories of range stability vs. fluctuation have caused cytonuclear concordance and discordance in ground squirrel lineages from the Otospermophilus beecheyi species complex. Previous studies in O. beecheyi revealed three morphologically cryptic and highly divergent mitochondrial DNA lineages (named the Northern, Central and Southern lineages based on geography) with only the Northern lineage exhibiting concordant divergence for nuclear genes. Here, we showed that these mtDNA lineages likely formed in allopatry during the Pleistocene, but responded differentially to climatic changes that occurred since the last interglacial (~120,000 years ago). We find that the Northern lineage maintained a stable range throughout this period, correlating with genetic distinctiveness among all genetic markers and low migration rates with the other lineages. In contrast, our results suggested that the Southern lineage expanded from Baja California Sur during the Late Pleistocene to overlap and potentially swamp a contracting Central lineage. High rates of intraspecific gene flow between Southern lineage individuals among expansion origin and expansion edge populations largely eroded Central ancestry from autosomal markers. However, male-biased dispersal in this system preserved signals of this past hybridization and introgression event in matrilineal-biased X-chromosome and mtDNA markers. Our results highlight the importance of range stability in maintaining the persistence of phylogeographic lineages, whereas unstable range dynamics can increase the tendency for lineages to merge upon secondary contact.

Dietary breadth is positively correlated with venom complexity in cone snails.

BACKGROUND: Although diet is believed to be a major factor underlying the evolution of venom, few comparative studies examine both venom composition and diet across a radiation of venomous species. Cone snails within the family, Conidae, comprise more than 700 species of carnivorous marine snails that capture their prey by using a cocktail of venomous neurotoxins (conotoxins or conopeptides). Venom composition across species has been previously hypothesized to be shaped by (a) prey taxonomic class (i.e., worms, molluscs, or fish) and (b) dietary breadth. We tested these hypotheses under a comparative phylogenetic framework using ecological data from past studies in conjunction with venom duct transcriptomes sequenced from 12 phylogenetically disparate cone snail species, including 10 vermivores (worm-eating), one molluscivore, and one generalist. RESULTS: We discovered 2223 unique conotoxin precursor peptides that encoded 1864 unique mature toxins across all species, >90 % of which are new to this study. In addition, we identified two novel gene superfamilies and 16 novel cysteine frameworks. Each species exhibited unique venom profiles, with venom composition and expression patterns among species dominated by a restricted set of gene superfamilies and mature toxins. In contrast with the dominant paradigm for interpreting Conidae venom evolution, prey taxonomic class did not predict venom composition patterns among species. We also found a significant positive relationship between dietary breadth and measures of conotoxin complexity. CONCLUSIONS: The poor performance of prey taxonomic class in predicting venom components suggests that cone snails have either evolved species-specific expression patterns likely as a consequence of the rapid evolution of conotoxin genes, or that traditional means of categorizing prey type (i.e., worms, mollusc, or fish) and conotoxins (i.e., by gene superfamily) do not accurately encapsulate evolutionary dynamics between diet and venom composition. We also show that species with more generalized diets tend to have more complex venoms and utilize a greater number of venom genes for prey capture. Whether this increased gene diversity confers an increased capacity for evolutionary change remains to be tested. Overall, our results corroborate the key role of diet in influencing patterns of venom evolution in cone snails and other venomous radiations.

Nitric oxide signaling differentially affects habitat choice by two larval morphs of the sea slug Alderia willowi: mechanistic insight into evolutionary transitions in dispersal strategies.

In many marine animals, adult habitat is selected by lecithotrophic (non-feeding) larvae with a limited lifespan. In generalist species, larvae may increasingly accept sub-optimal habitat over time as energy stores are depleted ('desperate larva' hypothesis). If the fitness cost of suboptimal habitat is too high, larvae of specialists may prolong the searching phase until they encounter a high-quality patch or die ('death before dishonor' hypothesis). In generalists, starvation is hypothesized to lead to a decline in inhibitory nitric oxide (NO) signaling, thereby triggering metamorphosis. Here, we document alternative functions for identified signaling pathways in larvae having 'desperate' versus 'death before dishonor' strategies in lecithotrophic clutches of a habitat specialist, the sea slug Alderia willowi. In an unusual dimorphism, each clutch of A. willowi hatches both non-selective larvae that settle soon after hatching and siblings that delay settlement in the absence of cues from the alga Vaucheria, the sole adult food. Pharmacological manipulation of NO signaling induced metamorphosis in non-selective but not selective stages. However, decreased NO signaling in selective larvae lowered the threshold for response to habitat cues, mimicking the effect of declining energy levels. Manipulation of cGMP or dopamine production induced metamorphosis in selective and non-selective larvae alike, highlighting a distinct role for the NO pathway in the two larval morphs. We propose a model in which NO production (1) links nitrogen metabolism with sensory receptor signaling, and (2) shifts from a regulatory role in 'desperate larva' strategies to a modulatory role in 'death before dishonor' strategies. This study provides new mechanistic insight into how the function of conserved signaling pathways may change in response to selection on larval habitat choice behaviors.

Selecting Potential Neuronal Drug Leads from Conotoxins of Various Venomous Marine Cone Snails in Bali, Indonesia.

Many conotoxins, natural peptides of marine cone snails, have been identified to target neurons. Here, we provide data on pharmacological families of the conotoxins of 11 species of cone snails collected in Bali. The identified definitive pharmacological families possibly targeting neuronal tissues were α (alpha), ι (iota), κ (kappa), and ρ (rho). These classes shall target nicotinic acetylcholine receptors, voltage-gated Na channels, voltage-gated K channels, and α1-adrenoceptors, respectively. The VI/VII-O3 conotoxins might be prospected as an inhibitor of N-methyl-d-aspartate. Con-ikot-ikot could be applied as an α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor blocker medicine. The definitive pharmacology classes of conotoxins as well as those yet to be elucidated need to be further established and verified.