Evolution of Conus peptide genes: duplication and positive selection in the A-superfamily

Chemical Synthesis and Insecticidal Activity Research Based on α-Conotoxins

The escalating resistance of agricultural pests to chemical insecticides necessitates the development of novel, efficient, and safe biological insecticides. Conus quercinus, a vermivorous cone snail, yields a crude venom rich in peptides for marine worm predation. This study screened six α-conotoxins with insecticidal potential from a previously constructed transcriptome database of C. quercinus, characterized by two disulfide bonds. These conotoxins were derived via solid-phase peptide synthesis (SPPS) and folded using two-step iodine oxidation for further insecticidal activity validation, such as CCK-8 assay and insect bioassay. The final results confirmed the insecticidal activities of the six α-conotoxins, with Qc1.15 and Qc1.18 exhibiting high insecticidal activity. In addition, structural analysis via homology modeling and functional insights from molecular docking offer a preliminary look into their potential insecticidal mechanisms. In summary, this study provides essential references and foundations for developing novel insecticides.

The Sobering Sting: Oleoyl Serotonin Is a Novel Stephanoconus Snail Venom-Derived Antagonist of Cannabinoid Receptors That Counteracts Learning and Memory Deficits

Cannabinoid receptors (CB1 and CB2) are promising targets for a better understanding of neurological diseases. Nevertheless, only a few ligands of CB have reached clinical application so far. Venoms are considered as interesting sources of novel biologically active compounds. Here, we describe an endocannabinoid-like molecule, oleoyl serotonin (OS), present in the venom of Stephanoconus snails. Using electrophysiological assays, it was shown that OS inhibits CB1 and CB2. Structure-activity relationship studies using a chimeric CB1/2 revealed that the domain encompassing the transmembrane helix V (TMHV)- intracellular loop 3 (ICL3)-TMHVI of the CB2 is critical for the binding and function of OS. We concluded that OS binds to sites of the CB2 that are different from the binding sites of the non-selective CB agonist WIN55,212-2. Behavioral assays in mice showed that OS counteracted learning and memory deficits caused by WIN55,212-2. Furthermore, a selectivity screening of OS showed high selectivity for CB over various ion channels and receptors. Overall, OS may represent a new approach to the prevention and treatment of learning and memory cognition impairment in neurological diseases.

Proteomic Analysis of the Predatory Venom of Conus striatus Reveals Novel and Population-Specific κA-Conotoxin SIVC

Animal venoms are a rich source of pharmacological compounds with ecological and evolutionary significance, as well as with therapeutic and biotechnological potentials. Among the most promising venomous animals, cone snails produce potent neurotoxic venom to facilitate prey capture and defend against aggressors. Conus striatus, one of the largest piscivorous species, is widely distributed, from east African coasts to remote Polynesian Islands. In this study, we investigated potential intraspecific differences in venom composition between distinct geographical populations from Mayotte Island (Indian Ocean) and Australia (Pacific Ocean). Significant variations were noted among the most abundant components, namely the κA-conotoxins, which contain three disulfide bridges and complex glycosylations. The amino acid sequence of a novel κA-conotoxin SIVC, including its N-terminal acetylated variant, was deciphered using tandem mass spectrometry (MS/MS). In addition, the glycosylation pattern was found to be consisting of two HexNAc and four Hex for the Mayotte population, which diverge from the previously characterized two HexNAc and three Hex combinations for this species, collected elsewhere. Whereas the biological and ecological roles of these modifications remain to be investigated, population-specific glycosylation patterns provide, for the first time, a new level of intraspecific variations in cone snail venoms.

Integrating Venom Peptide Libraries Into a Phylogenetic and Broader Biological Framework

The venomous marine snails are conventionally divided into three groups, the cone snails (family Conidae), the auger snails (family Terebridae) and the turrids (formerly all assigned to a single family, Turridae). In this study, a library of venom peptides from species conventionally assigned to the genus Turris was correlated to a phylogenetic analysis. Nucleotide sequences of multiple genes from transcriptomes were used to assess the phylogenetic relationships across a diverse set of species. The resulting tree shows that as conventionally defined, the conoidean genus Turris, is polyphyletic. We describe a new genus, Purpuraturris gen. nov., that comprises the outlier species. In addition to morphological distinctions, molecular data reveal that this group is divergent from Turris sensu stricto. The correlation between phylogenetic information and a family of peptide sequences was used to highlight those peptides mostly likely to be unique and intimately associated with biological diversity. The plethora of peptide sequences available requires two prioritization decisions: which subset of peptides to initially characterize, and after these are characterized, which to comprehensively investigate for potential biomedical applications such as drug developments. Life Science Identifiers: urn:lsid:zoobank.org; pub: 60D46561-28F0-4C39-BAC4-66DC8B4EAEA4.

A phylogeny-aware approach reveals unexpected venom components in divergent lineages of cone snails

Marine gastropods of the genus Conus are renowned for their remarkable diversity and deadly venoms. While Conus venoms are increasingly well studied for their biomedical applications, we know surprisingly little about venom composition in other lineages of Conidae. We performed comprehensive venom transcriptomic profiling for Conasprella coriolisi and Pygmaeconus traillii, first time for both respective genera. We complemented reference-based transcriptome annotation by a de novo toxin prediction guided by phylogeny, which involved transcriptomic data on two additional 'divergent' cone snail lineages, Profundiconus, and Californiconus. We identified toxin clusters (SSCs) shared among all or some of the four analysed genera based on the identity of the signal region-a molecular tag present in toxins. In total, 116 and 98 putative toxins represent 29 and 28 toxin gene superfamilies in Conasprella and Pygmaeconus, respectively; about quarter of these only found by semi-manual annotation of the SSCs. Two rare gene superfamilies, originally identified from fish-hunting cone snails, were detected outside Conus rather unexpectedly, so we further investigated their distribution across Conidae radiation. We demonstrate that both these, in fact, are ubiquitous in Conidae, sometimes with extremely high expression. Our findings demonstrate how a phylogeny-aware approach circumvents methodological caveats of similarity-based transcriptome annotation.

Venom duct origins of prey capture and defensive conotoxins in piscivorous Conus striatus

The venom duct origins of predatory and defensive venoms has not been studied for hook-and-line fish hunting cone snails despite the pharmacological importance of their venoms. To better understand the biochemistry and evolution of injected predatory and defensive venoms, we compared distal, central and proximal venom duct sections across three specimens of C. striatus (Pionoconus) using proteomic and transcriptomic approaches. A total of 370 conotoxin precursors were identified from the whole venom duct transcriptome. Milked defensive venom was enriched with a potent cocktail of proximally expressed inhibitory α-, ω- and μ-conotoxins compared to milked predatory venom. In contrast, excitatory κA-conotoxins dominated both the predatory and defensive venoms despite their distal expression, suggesting this class of conotoxin can be selectively expressed from the same duct segment in response to either a predatory or defensive stimuli. Given the high abundance of κA-conotoxins in the Pionoconus clade, we hypothesise that the κA-conotoxins have evolved through adaptive evolution following their repurposing from ancestral inhibitory A superfamily conotoxins to facilitate the dietary shift to fish hunting and species radiation in this clade.

The evolutionary dynamics of venom toxins made by insects and other animals

Animal venoms are recognised as unique biological systems in which to study molecular evolution. Venom use has evolved numerous times among the insects, and insects today use venom to capture prey, defend themselves from predators, or to subdue and modulate host responses during parasitism. However, little is known about most insect venom toxins or the mode and tempo by which they evolve. Here, I review the evolutionary dynamics of insect venom toxins, and argue that insects offer many opportunities to examine novel aspects of toxin evolution. The key questions addressed are: How do venomous animals evolve from non-venomous animals, and how does this path effect the composition and pharmacology of the venom? What genetic processes (gene duplication, co-option, neofunctionalisation) are most important in toxin evolution? What kinds of selection pressures are acting on toxin-encoding genes and their cognate targets in envenomated animals? The emerging evidence highlights that venom composition and pharmacology adapts quickly in response to changing selection pressures resulting from new ecological interactions, and that such evolution occurs through a stunning variety of genetic mechanisms. Insects offer many opportunities to investigate the evolutionary dynamics of venom toxins due to their evolutionary history rich in venom-related adaptations, and their quick generation time and suitability for culture in the laboratory.

Non-Peptidic Small Molecule Components from Cone Snail Venoms

Venomous molluscs (Superfamily Conoidea) comprise a substantial fraction of tropical marine biodiversity (>15,000 species). Prior characterization of cone snail venoms established that bioactive venom components used to capture prey, defend against predators and for competitive interactions were relatively small, structured peptides (10–35 amino acids), most with multiple disulfide crosslinks. These venom components (“conotoxins, conopeptides”) have been widely studied in many laboratories, leading to pharmaceutical agents and probes. In this review, we describe how it has recently become clear that to varying degrees, cone snail venoms also contain bioactive non-peptidic small molecule components. Since the initial discovery of genuanine as the first bioactive venom small molecule with an unprecedented structure, a broad set of cone snail venoms have been examined for non-peptidic bioactive components. In particular, a basal clade of cone snails (Stephanoconus) that prey on polychaetes produce genuanine and many other small molecules in their venoms, suggesting that this lineage may be a rich source of non-peptidic cone snail venom natural products. In contrast to standing dogma in the field that peptide and proteins are predominantly used for prey capture in cone snails, these small molecules also contribute to prey capture and push the molecular diversity of cone snails beyond peptides. The compounds so far characterized are active on neurons and thus may potentially serve as leads for neuronal diseases. Thus, in analogy to the incredible pharmacopeia resulting from studying venom peptides, these small molecules may provide a new resource of pharmacological agents.

Curses or Cures: A Review of the Numerous Benefits Versus the Biosecurity Concerns of Conotoxin Research

Conotoxins form a diverse group of peptide toxins found in the venom of predatory marine cone snails. Decades of conotoxin research have provided numerous measurable scientific and societal benefits. These include their use as a drug, diagnostic agent, drug leads, and research tools in neuroscience, pharmacology, biochemistry, structural biology, and molecular evolution. Human envenomations by cone snails are rare but can be fatal. Death by envenomation is likely caused by a small set of toxins that induce muscle paralysis of the diaphragm, resulting in respiratory arrest. The potency of these toxins led to concerns regarding the potential development and use of conotoxins as biological weapons. To address this, various regulatory measures have been introduced that limit the use and access of conotoxins within the research community. Some of these regulations apply to all of the ≈200,000 conotoxins predicted to exist in nature of which less than 0.05% are estimated to have any significant toxicity in humans. In this review we provide an overview of the many benefits of conotoxin research, and contrast these to the perceived biosecurity concerns of conotoxins and research thereof.

Venomic Interrogation Reveals the Complexity of Conus striolatus Venom

Given the complexity of cone snail venoms, high throughput venomics approaches are required to fully investigate venom composition, envenomation strategies, and evolutionary trajectories. This study describes 158 conotoxins in the venom transcriptome of the little studied C. striolatus from the fish hunting clade Pionoconus. Despite similar gene superfamily distributions along the venom duct, only 18 common transcripts were identified between distal, central, and proximal venom duct transcriptomes. Proteomic analysis of the injected predatory venom collected from the same individual revealed an ~18-fold enhanced complexity at the proteomic level, consistent with complex post-translational modifications and variable venom peptide processing occurring in the venom duct. Overall, C. striolatus venom was dominated by M, O1, O2, and A gene superfamily conotoxins and conkunitzins, which are potential modulators of sodium, calcium, and potassium channels. Conkunitzins and gene superfamily A peptides dominated the proximal over the distal duct, the M and O1 gene superfamily peptides were distributed along the full length of the duct, while the O2 gene superfamily peptides dominated the distal duct. Interestingly, the predatory injected venom of C. striolatus was dominated by peptides from gene superfamilies M, O1, O2, A, and conkunitzins, suggesting the predatory venom of C. striolatus may arise at multiple sites along the venom duct.

Conotoxin Diversity in the Venom Gland Transcriptome of the Magician's Cone, Pionoconus magus

The transcriptomes of the venom glands of two individuals of the magician’s cone, Pionoconus magus, from Okinawa (Japan) were sequenced, assembled, and annotated. In addition, RNA-seq raw reads available at the SRA database from one additional specimen of P. magus from the Philippines were also assembled and annotated. The total numbers of identified conotoxin precursors and hormones per specimen were 118, 112, and 93. The three individuals shared only five identical sequences whereas the two specimens from Okinawa had 30 sequences in common. The total number of distinct conotoxin precursors and hormones for P. magus was 275, and were assigned to 53 conotoxin precursor and hormone superfamilies, two of which were new based on their divergent signal region. The superfamilies that had the highest number of precursors were M (42), O1 (34), T (27), A (18), O2 (17), and F (13), accounting for 55% of the total diversity. The D superfamily, previously thought to be exclusive of vermivorous cones was found in P. magus and contained a highly divergent mature region. Similarly, the A superfamily alpha 4/3 was found in P. magus despite the fact that it was previously postulated to be almost exclusive of the genus Rhombiconus. Differential expression analyses of P. magus compared to Chelyconus ermineus, the only fish-hunting cone from the Atlantic Ocean revealed that M and A2 superfamilies appeared to be more expressed in the former whereas the O2 superfamily was more expressed in the latter.

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.

Scorpion Toxins: Positive Selection at a Distal Site Modulates Functional Evolution at a Bioactive Site

The bioactive sites of proteins are those that directly interact with their targets. In many immunity- and predation-related proteins, they frequently experience positive selection for dealing with the changes of their targets from competitors. However, some sites that are far away from the interface between proteins and their targets are also identified to evolve under positive selection. Here, we explore the evolutionary implication of such a site in scorpion α-type toxins affecting sodium (Na⁺) channels (abbreviated as α-ScNaTxs) using a combination of experimental and computational approaches. We found that despite no direct involvement in interaction with Na⁺ channels, mutations at this site by different types of amino acids led to toxicity change on both rats and insects in three α-ScNaTxs, accompanying differential effects on their structures. Molecular dynamics simulations indicated that the mutations changed the conformational dynamics of the positively selected bioactive site-containing functional regions by allosteric communication, suggesting a potential evolutionary correlation between these bioactive sites and the distant nonbioactive site. Our results reveal for the first time the cause of fast evolution at nonbioactive sites of scorpion neurotoxins, which is presumably to adapt to the change of their bioactive sites through coevolution to maintain an active conformation for channel binding. This might aid rational design of scorpion Na⁺ channel toxins with improved phyletic selectivity via modification of a distant nonbioactive site.

High Throughput Identification of Novel Conotoxins from the Vermivorous Oak Cone Snail (Conus quercinus) by Transcriptome Sequencing

The primary objective of this study was to realize the large-scale discovery of conotoxin sequences from different organs (including the venom duct, venom bulb and salivary gland) of the vermivorous Oak cone snail, Conus quercinus. Using high-throughput transcriptome sequencing, we identified 133 putative conotoxins that belong to 34 known superfamilies, of which nine were previously reported while the remaining 124 were novel conotoxins, with 17 in new and unassigned conotoxin groups. A-, O₁-, M-, and I₂- superfamilies were the most abundant, and the cysteine frameworks XIII and VIII were observed for the first time in the A- and I₂-superfamilies. The transcriptome data from the venom duct, venom bulb and salivary gland showed considerable inter-organizational variations. Each organ had many exclusive conotoxins, and only seven of all the inferred mature peptides were common in the three organs. As expected, most of the identified conotoxins were synthesized in the venom duct at relatively high levels; however, a number of conotoxins were also identified in the venom bulb and the salivary gland with very low transcription levels. Therefore, various organs have different conotoxins with high diversity, suggesting greater contributions from several organs to the high-throughput discovery of new conotoxins for future drug development.

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.

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.

Venom Insulins of Cone Snails Diversify Rapidly and Track Prey Taxa

A specialized insulin was recently found in the venom of a fish-hunting cone snail, Conus geographus Here we show that many worm-hunting and snail-hunting cones also express venom insulins, and that this novel gene family has diversified explosively. Cone snails express a highly conserved insulin in their nerve ring; presumably this conventional signaling insulin is finely tuned to the Conus insulin receptor, which also evolves very slowly. By contrast, the venom insulins diverge rapidly, apparently in response to biotic interactions with prey and also possibly the cones' own predators and competitors. Thus, the inwardly directed signaling insulins appear to experience predominantly purifying sele\ction to target an internal receptor that seldom changes, while the outwardly directed venom insulins frequently experience directional selection to target heterospecific insulin receptors in a changing mix of prey, predators and competitors. Prey insulin receptors may often be constrained in ways that prevent their evolutionary escape from targeted venom insulins, if amino-acid substitutions that result in escape also degrade the receptor's signaling functions.

Target-Driven Positive Selection at Hot Spots of Scorpion Toxins Uncovers Their Potential in Design of Insecticides

Positive selection sites (PSSs), a class of amino acid sites with an excess of nonsynonymous to synonymous substitutions, are indicators of adaptive molecular evolution and have been detected in many protein families involved in a diversity of biological processes by statistical approaches. However, few studies are conducted to evaluate their functional significance and the driving force behind the evolution (i.e., agent of selection). Scorpion α-toxins are a class of multigene family of peptide neurotoxins affecting voltage-gated Na(+ )(Nav) channels, whose members exhibit differential potency and preference for insect and mammalian Nav channels. In this study, we undertook a systematical molecular dissection of nearly all the PSSs newly characterized in the Mesobuthus α-toxin family and a two-residue insertion ((19)AlaPhe(20)) located within a positively selected loop via mutational analysis of α-like MeuNaTxα-5, one member affecting both insect and mammalian Nav channels. This allows to identify hot-spot residues on its functional face involved in interaction with the receptor site of Nav channels, which comprises two PSSs (Ile(40) and Leu(41)) and the small insertion, both located on two spatially separated functional loops. Mutations at these hot-spots resulted in a remarkably decreased anti-mammalian activity in MeuNaTxα-5 with partially impaired or enhanced insecticide activity, suggesting the potential of PSSs in designing promising candidate insecticides from scorpion α-like toxins. Based on an experiment-guided toxin-channel complex model and high evolutionary variability in the receptor site of predators and prey of scorpions, we provide new evidence for target-driven adaptive evolution of scorpion toxins to deal with their targets' diversity.

Rapid expansion of the protein disulfide isomerase gene family facilitates the folding of venom peptides

Formation of correct disulfide bonds in the endoplasmic reticulum is a crucial step for folding proteins destined for secretion. Protein disulfide isomerases (PDIs) play a central role in this process. We report a previously unidentified, hypervariable family of PDIs that represents the most diverse gene family of oxidoreductases described in a single genus to date. These enzymes are highly expressed specifically in the venom glands of predatory cone snails, animals that synthesize a remarkably diverse set of cysteine-rich peptide toxins (conotoxins). Enzymes in this PDI family, termed conotoxin-specific PDIs, significantly and differentially accelerate the kinetics of disulfide-bond formation of several conotoxins. Our results are consistent with a unique biological scenario associated with protein folding: The diversification of a family of foldases can be correlated with the rapid evolution of an unprecedented diversity of disulfide-rich structural domains expressed by venomous marine snails in the superfamily Conoidea.

Structural features of conopeptide genes inferred from partial sequences of the Conus tribblei genome

The evolvability of venom components (in particular, the gene-encoded peptide toxins) in venomous species serves as an adaptive strategy allowing them to target new prey types or respond to changes in the prey field. The structure, organization, and expression of the venom peptide genes may provide insights into the molecular mechanisms that drive the evolution of such genes. Conus is a particularly interesting group given the high chemical diversity of their venom peptides, and the rapid evolution of the conopeptide-encoding genes. Conus genomes, however, are large and characterized by a high proportion of repetitive sequences. As a result, the structure and organization of conopeptide genes have remained poorly known. In this study, a survey of the genome of Conus tribblei was undertaken to address this gap. A partial assembly of C. tribblei genome was generated; the assembly, though consisting of a large number of fragments, accounted for 2160.5 Mb of sequence. A large number of repetitive genomic elements consisting of 642.6 Mb of retrotransposable elements, simple repeats, and novel interspersed repeats were observed. We characterized the structural organization and distribution of conotoxin genes in the genome. A significant number of conopeptide genes (estimated to be between 148 and 193) belonging to different superfamilies with complete or nearly complete exon regions were observed, ~60 % of which were expressed. The unexpressed conopeptide genes represent hidden but significant conotoxin diversity. The conotoxin genes also differed in the frequency and length of the introns. The interruption of exons by long introns in the conopeptide genes and the presence of repeats in the introns may indicate the importance of introns in facilitating recombination, evolution and diversification of conotoxins. These findings advance our understanding of the structural framework that promotes the gene-level molecular evolution of venom peptides.

Comparative Venomics Reveals the Complex Prey Capture Strategy of the Piscivorous Cone Snail Conus catus

Venomous marine cone snails produce a unique and remarkably diverse range of venom peptides (conotoxins and conopeptides) that have proven to be invaluable as pharmacological probes and leads to new therapies. Conus catus is a hook-and-line fish hunter from clade I, with ∼20 conotoxins identified, including the analgesic ω-conotoxin CVID (AM336). The current study unravels the venom composition of C. catus with tandem mass spectrometry and 454 sequencing data. From the venom gland transcriptome, 104 precursors were recovered from 11 superfamilies, with superfamily A (especially κA-) conotoxins dominating (77%) their venom. Proteomic analysis confirmed that κA-conotoxins dominated the predation-evoked milked venom of each of six C. catus analyzed and revealed remarkable intraspecific variation in both the intensity and type of conotoxins. High-throughput FLIPR assays revealed that the predation-evoked venom contained a range of conotoxins targeting the nAChR, Cav, and Nav ion channels, consistent with α- and ω-conotoxins being used for predation by C. catus. However, the κA-conotoxins did not act at these targets but induced potent and rapid immobilization followed by bursts of activity and finally paralysis when injected intramuscularly in zebrafish. Our venomics approach revealed the complexity of the envenomation strategy used by C. catus, which contains a mix of both excitatory and inhibitory venom peptides.

Comparison of the Venom Peptides and Their Expression in Closely Related Conus Species: Insights into Adaptive Post-speciation Evolution of Conus Exogenomes

Genes that encode products with exogenous targets, which comprise an organism's "exogenome," typically exhibit high rates of evolution. The genes encoding the venom peptides (conotoxins or conopeptides) in Conus sensu lato exemplify this class of genes. Their rapid diversification has been established and is believed to be linked to the high speciation rate in this genus. However, the molecular mechanisms that underlie venom peptide diversification and ultimately emergence of new species remain poorly understood. In this study, the sequences and expression levels of conotoxins from several specimens of two closely related worm-hunting species, Conus tribblei and Conus lenavati, were compared through transcriptome analysis. Majority of the identified putative conopeptides were novel, and their diversity, even in each specimen, was remarkably high suggesting a wide range of prey targets for these species. Comparison of the interspecific expression patterns of conopeptides at the superfamily level resulted in the discovery of both conserved as well as species-specific expression patterns, indicating divergence in the regulatory network affecting conotoxin gene expression. Comparison of the transcriptomes of the individual snails revealed that each specimen produces a distinct set of highly expressed conopeptides, reflecting the capability of individual snails to fine-tune the composition of their venoms. These observations reflect the role of sequence divergence and divergence in the control of expression for specific conopeptides in the evolution of the exogenome and hence venom composition in Conus.

One, four or 100 genera? A new classification of the cone snails

We present a new classification for the genus Conus sensu lato (family Conidae), based on molecular phylogenetic analyses of 329 species. This classification departs from both the traditional classification in only one genus and from a recently proposed shell- and radula-based classification scheme that separates members of this group into five families and 115 genera. Roughly 140 genus-group names are available for Recent cone snails. We propose to place all cone snails within a single family (Conidae) containing four genera-Conus, Conasprella, Profundiconus and Californiconus (with Conus alone encompassing about 85% of known species)-based on the clear separation of cone snails into four distinct and well-supported groups/lineages in molecular phylogenetic analyses. Within Conus and Conasprella, we recognize 57 and 11 subgenera, respectively, that represent well-supported subgroupings within these genera, which we interpret as evidence of intrageneric distinctiveness. We allocate the 803 Recent species of Conidae listed as valid in the World Register of Marine Species into these four genera and 71 subgenera, with an estimate of the confidence for placement of species in these taxonomic categories based on whether molecular or radula and/or shell data were used in these determinations. Our proposed classification effectively departs from previous schemes by (1) limiting the number of accepted genera, (2) retaining the majority of species within the genus Conus and (3) assigning members of these genera to species groups/subgenera to enable the effective communication of these groups, all of which we hope will encourage acceptance of this scheme.

Conotoxin gene superfamilies

Conotoxins are the peptidic components of the venoms of marine cone snails (genus Conus). They are remarkably diverse in terms of structure and function. Unique potency and selectivity profiles for a range of neuronal targets have made several conotoxins valuable as research tools, drug leads and even therapeutics, and has resulted in a concerted and increasing drive to identify and characterise new conotoxins. Conotoxins are translated from mRNA as peptide precursors, and cDNA sequencing is now the primary method for identification of new conotoxin sequences. As a result, gene superfamily, a classification based on precursor signal peptide identity, has become the most convenient method of conotoxin classification. Here we review each of the described conotoxin gene superfamilies, with a focus on the structural and functional diversity present in each. This review is intended to serve as a practical guide to conotoxin superfamilies and to facilitate interpretation of the increasing number of conotoxin precursor sequences being identified by targeted-cDNA sequencing and more recently high-throughput transcriptome sequencing.

Molecular phylogeny and evolution of the cone snails (Gastropoda, Conoidea)

We present a large-scale molecular phylogeny that includes 320 of the 761 recognized valid species of the cone snails (Conus), one of the most diverse groups of marine molluscs, based on three mitochondrial genes (COI, 16S rDNA and 12S rDNA). This is the first phylogeny of the taxon to employ concatenated sequences of several genes, and it includes more than twice as many species as the last published molecular phylogeny of the entire group nearly a decade ago. Most of the numerous molecular phylogenies published during the last 15years are limited to rather small fractions of its species diversity. Bayesian and maximum likelihood analyses are mostly congruent and confirm the presence of three previously reported highly divergent lineages among cone snails, and one identified here using molecular data. About 85% of the species cluster in the single Large Major Clade; the others are divided between the Small Major Clade (∼12%), the Conus californicus lineage (one species), and a newly defined clade (∼3%). We also define several subclades within the Large and Small major clades, but most of their relationships remain poorly supported. To illustrate the usefulness of molecular phylogenies in addressing specific evolutionary questions, we analyse the evolution of the diet, the biogeography and the toxins of cone snails. All cone snails whose feeding biology is known inject venom into large prey animals and swallow them whole. Predation on polychaete worms is inferred as the ancestral state, and diet shifts to molluscs and fishes occurred rarely. The ancestor of cone snails probably originated from the Indo-Pacific; rather few colonisations of other biogeographic provinces have probably occurred. A new classification of the Conidae, based on the molecular phylogeny, is published in an accompanying paper.

Application of community phylogenetic approaches to understand gene expression: differential exploration of venom gene space in predatory marine gastropods

Background

Predatory marine gastropods of the genus Conus exhibit substantial variation in venom composition both within and among species. Apart from mechanisms associated with extensive turnover of gene families and rapid evolution of genes that encode venom components (‘conotoxins’), the evolution of distinct conotoxin expression patterns is an additional source of variation that may drive interspecific differences in the utilization of species’ ‘venom gene space’. To determine the evolution of expression patterns of venom genes of Conus species, we evaluated the expression of A-superfamily conotoxin genes of a set of closely related Conus species by comparing recovered transcripts of A-superfamily genes that were previously identified from the genomes of these species. We modified community phylogenetics approaches to incorporate phylogenetic history and disparity of genes and their expression profiles to determine patterns of venom gene space utilization.

Results

Less than half of the A-superfamily gene repertoire of these species is expressed, and only a few orthologous genes are coexpressed among species. Species exhibit substantially distinct expression strategies, with some expressing sets of closely related loci (‘under-dispersed’ expression of available genes) while others express sets of more disparate genes (‘over-dispersed’ expression). In addition, expressed genes show higher d_N/d_S values than either unexpressed or ancestral genes; this implies that expression exposes genes to selection and facilitates rapid evolution of these genes. Few recent lineage-specific gene duplicates are expressed simultaneously, suggesting that expression divergence among redundant gene copies may be established shortly after gene duplication.

Conclusions

Our study demonstrates that venom gene space is explored differentially by Conus species, a process that effectively permits the independent and rapid evolution of venoms in these species.

Evolutionary origins of a bioactive peptide buried within Preproalbumin

The de novo evolution of proteins is now considered a frequented route for biological innovation, but the genetic and biochemical processes that lead to each newly created protein are often poorly documented. The common sunflower (Helianthus annuus) contains the unusual gene PawS1 (Preproalbumin with SFTI-1) that encodes a precursor for seed storage albumin; however, in a region usually discarded during albumin maturation, its sequence is matured into SFTI-1, a protease-inhibiting cyclic peptide with a motif homologous to unrelated inhibitors from legumes, cereals, and frogs. To understand how PawS1 acquired this additional peptide with novel biochemical functionality, we cloned PawS1 genes and showed that this dual destiny is over 18 million years old. This new family of mostly backbone-cyclic peptides is structurally diverse, but the protease-inhibitory motif was restricted to peptides from sunflower and close relatives from its subtribe. We describe a widely distributed, potential evolutionary intermediate PawS-Like1 (PawL1), which is matured into storage albumin, but makes no stable peptide despite possessing residues essential for processing and cyclization from within PawS1. Using sequences we cloned, we retrodict the likely stepwise creation of PawS1's additional destiny within a simple albumin precursor. We propose that relaxed selection enabled SFTI-1 to evolve its inhibitor function by converging upon a successful sequence and structure.

Evolution stings: the origin and diversification of scorpion toxin peptide scaffolds

The episodic nature of natural selection and the accumulation of extreme sequence divergence in venom-encoding genes over long periods of evolutionary time can obscure the signature of positive Darwinian selection. Recognition of the true biocomplexity is further hampered by the limited taxon selection, with easy to obtain or medically important species typically being the subject of intense venom research, relative to the actual taxonomical diversity in nature. This holds true for scorpions, which are one of the most ancient terrestrial venomous animal lineages. The family Buthidae that includes all the medically significant species has been intensely investigated around the globe, while almost completely ignoring the remaining non-buthid families. Australian scorpion lineages, for instance, have been completely neglected, with only a single scorpion species (Urodacus yaschenkoi) having its venom transcriptome sequenced. Hence, the lack of venom composition and toxin sequence information from an entire continent’s worth of scorpions has impeded our understanding of the molecular evolution of scorpion venom. The molecular origin, phylogenetic relationships and evolutionary histories of most scorpion toxin scaffolds remain enigmatic. In this study, we have sequenced venom gland transcriptomes of a wide taxonomical diversity of scorpions from Australia, including buthid and non-buthid representatives. Using state-of-art molecular evolutionary analyses, we show that a majority of CSα/β toxin scaffolds have experienced episodic influence of positive selection, while most non-CSα/β linear toxins evolve under the extreme influence of negative selection. For the first time, we have unraveled the molecular origin of the major scorpion toxin scaffolds, such as scorpion venom single von Willebrand factor C-domain peptides (SV-SVC), inhibitor cystine knot (ICK), disulphide-directed beta-hairpin (DDH), bradykinin potentiating peptides (BPP), linear non-disulphide bridged peptides and antimicrobial peptides (AMP). We have thus demonstrated that even neglected lineages of scorpions are a rich pool of novel biochemical components, which have evolved over millions of years to target specific ion channels in prey animals, and as a result, possess tremendous implications in therapeutics.

Molecular evolution and diversity of Conus peptide toxins, as revealed by gene structure and intron sequence analyses

Cone snails, which are predatory marine gastropods, produce a cocktail of venoms used for predation, defense and competition. The major venom component, conotoxin, has received significant attention because it is useful in neuroscience research, drug development and molecular diversity studies. In this study, we report the genomic characterization of nine conotoxin gene superfamilies from 18 Conus species and investigate the relationships among conotoxin gene structure, molecular evolution and diversity. The I1, I2, M, O2, O3, P, S, and T superfamily precursors all contain three exons and two introns, while A superfamily members contain two exons and one intron. The introns are conserved within a certain gene superfamily, and also conserved across different Conus species, but divergent among different superfamilies. The intronic sequences contain many simple repeat sequences and regulatory elements that may influence conotoxin gene expression. Furthermore, due to the unique gene structure of conotoxins, the base substitution rates and the number of positively selected sites vary greatly among exons. Many more point mutations and trinucleotide indels were observed in the mature peptide exon than in the other exons. In addition, the first example of alternative splicing in conotoxin genes was found. These results suggest that the diversity of conotoxin genes has been shaped by point mutations and indels, as well as rare gene recombination or alternative splicing events, and that the unique gene structures could have made a contribution to the evolution of conotoxin genes.

On the importance of oxidative folding in the evolution of conotoxins: cysteine codon preservation through gene duplication and adaptation

Conotoxin genes are among the most rapidly evolving genes currently known; however, despite the well-established hypervariability of the intercysteine loops, the cysteines demonstrate significant conservation, with a site-specific codon bias for each cysteine in a family of conotoxins. Herein we present a novel rationale behind the codon-level conservation of the cysteines that comprise the disulfide scaffold. We analyze cysteine codon conservation using an internal reference and phylogenetic tools; our results suggest that the established codon conservation can be explained as the result of selective pressures linked to the production efficiency and folding of conotoxins, driving the conservation of cysteine at the amino-acid level. The preservation of cysteine has resulted in maintenance of the ancestral codon in most of the daughter lineages, despite the hypervariability of adjacent residues. We propose that the selective pressures acting on the venom components of cone snails involve an interplay of biosynthetic efficiency, activity at the target receptor and the importance of that activity to effective prey immobilization. Functional redundancy in the venom can thus serve as a buffer for the energy expenditure of venom production.

Characterizing the evolution and functions of the M-superfamily conotoxins

Conotoxins from cone snails are valuable in physiology research and therapeutic applications. Evolutionary mechanisms of conotoxins have been investigated in several superfamilies, but there is no phylogenetic analysis on M-superfamily conotoxins. In this study, we characterized identical sequences, gene structure, novel cysteine frameworks, functions and evolutionary mechanisms of M-superfamily conotoxins. Identical M-superfamily conotoxins can be found in different Conus species from the analysis of novel 467 M-superfamily conotoxin sequences and other published M-superfamily conotoxins sequences. M-superfamily conotoxin genes consist of two introns and three exons from the results of genome walking. Eighteen cysteine frameworks were identified from the M-superfamily conotoxins, and 10 of the 18 may be generated from framework III. An analysis between diet types and phylogeny of the M-superfamily conotoxins indicate that M-superfamily conotoxins might not evolve in a concerted manner but were subject to birth-and-death evolution. Codon usage analysis shows that position-specific codon conservation is not restricted to cysteines, but also to other conserved residues. By analysing primary structures and physiological functions of M-superfamily conotoxins, we proposed a hypothesis that insertions and deletions, especially insertions in the third cysteine loop, are involved in the creation of new functions and structures of the M-superfamily conotoxins.

Structure of the O-glycosylated conopeptide CcTx from Conus consors venom

The glycopeptide CcTx, isolated from the venom of the piscivorous cone snail Conus consors, belongs to the κA-family of conopeptides. These toxins elicit excitotoxic responses in the prey by acting on voltage-gated sodium channels. The structure of CcTx, a first in the κA-family, has been determined by high-resolution NMR spectroscopy together with the analysis of its O-glycan at Ser7. A new type of glycopeptide O-glycan core structure, here registered as core type 9, containing two terminal L-galactose units {α-L-Galp-(1→4)-α-D-GlcpNAc-(1→6)-[α-L-Galp-(1→2)-β-D-Galp-(1→3)-]α-D-GalpNAc-(1→O)}, is highlighted. A sequence comparison to other putative members of the κA-family suggests that O-linked glycosylation might be more common than previously thought. This observation alone underlines the requirement for more careful and in-depth investigations into this type of post-translational modification in conotoxins.

Adaptive radiation of venomous marine snail lineages and the accelerated evolution of venom peptide genes

An impressive biodiversity (>10,000 species) of marine snails (suborder Toxoglossa or superfamily Conoidea) have complex venoms, each containing approximately 100 biologically active, disulfide-rich peptides. In the genus Conus, the most intensively investigated toxoglossan lineage (∼500 species), a small set of venom gene superfamilies undergo rapid sequence hyperdiversification within their mature toxin regions. Each major lineage of Toxoglossa has its own distinct set of venom gene superfamilies. Two recently identified venom gene superfamilies are expressed in the large Turridae clade, but not in Conus. Thus, as major venomous molluscan clades expand, a small set of lineage-specific venom gene superfamilies undergo accelerated evolution. The juxtaposition of extremely conserved signal sequences with hypervariable mature peptide regions is unprecedented and raises the possibility that in these gene superfamilies, the signal sequences are conserved as a result of an essential role they play in enabling rapid sequence evolution of the region of the gene that encodes the active toxin.

Molecular phylogeny, classification and evolution of conopeptides

Conopeptides are toxins expressed in the venom duct of cone snails (Conoidea, Conus). These are mostly well-structured peptides and mini-proteins with high potency and selectivity for a broad range of cellular targets. In view of these properties, they are widely used as pharmacological tools and many are candidates for innovative drugs. The conopeptides are primarily classified into superfamilies according to their peptide signal sequence, a classification that is thought to reflect the evolution of the multigenic system. However, this hypothesis has never been thoroughly tested. Here we present a phylogenetic analysis of 1,364 conopeptide signal sequences extracted from GenBank. The results validate the current conopeptide superfamily classification, but also reveal several important new features. The so-called "cysteine-poor" conopeptides are revealed to be closely related to "cysteine-rich" conopeptides; with some of them sharing very similar signal sequences, suggesting that a distinction based on cysteine content and configuration is not phylogenetically relevant and does not reflect the evolutionary history of conopeptides. A given cysteine pattern or pharmacological activity can be found across different superfamilies. Furthermore, a few conopeptides from GenBank do not cluster in any of the known superfamilies, and could represent yet-undefined superfamilies. A clear phylogenetically based classification should help to disentangle the diversity of conopeptides, and could also serve as a rationale to understand the evolution of the toxins in the numerous other species of conoideans and venomous animals at large.

A diverse family of novel peptide toxins from an unusual cone snail, Conus californicus

Diversity among Conus toxins mirrors the high species diversity in the Indo-Pacific region, and evolution of both is thought to stem from feeding-niche specialization derived from intra-generic competition. This study focuses on Conus californicus, a phylogenetic outlier endemic to the temperate northeast Pacific. Essentially free of congeneric competitors, it preys on a wider variety of organisms than any other cone snail. Using molecular cloning of cDNAs and mass spectrometry, we examined peptides isolated from venom ducts to elucidate the sequences and post-translational modifications of two eight-cysteine toxins (cal12a and cal12b of type 12 framework) that block voltage-gated Na(+) channels. Based on homology of leader sequence and mode of action, these toxins are related to the O-superfamily, but differ significantly from other members of that group. Six of the eight cysteine residues constitute the canonical framework of O-members, but two additional cysteine residues in the N-terminal region define an O+2 classification within the O-superfamily. Fifteen putative variants of Cal12.1 toxins have been identified by mRNAs that differ primarily in two short hypervariable regions and have been grouped into three subtypes (Cal12.1.1-3). This unique modular variation has not been described for other Conus toxins and suggests recombination as a diversity-generating mechanism. We propose that these toxin isoforms show specificity for similar molecular targets (Na(+) channels) in the many species preyed on by C. californicus and that individualistic utilization of specific toxin isoforms may involve control of gene expression.

Characterization of the Conus bullatus genome and its venom-duct transcriptome

BACKGROUND

The venomous marine gastropods, cone snails (genus Conus), inject prey with a lethal cocktail of conopeptides, small cysteine-rich peptides, each with a high affinity for its molecular target, generally an ion channel, receptor or transporter. Over the last decade, conopeptides have proven indispensable reagents for the study of vertebrate neurotransmission. Conus bullatus belongs to a clade of Conus species called Textilia, whose pharmacology is still poorly characterized. Thus the genomics analyses presented here provide the first step toward a better understanding the enigmatic Textilia clade.

RESULTS

We have carried out a sequencing survey of the Conus bullatus genome and venom-duct transcriptome. We find that conopeptides are highly expressed within the venom-duct, and describe an in silico pipeline for their discovery and characterization using RNA-seq data. We have also carried out low-coverage shotgun sequencing of the genome, and have used these data to determine its size, genome-wide base composition, simple repeat, and mobile element densities.

CONCLUSIONS

Our results provide the first global view of venom-duct transcription in any cone snail. A notable feature of Conus bullatus venoms is the breadth of A-superfamily peptides expressed in the venom duct, which are unprecedented in their structural diversity. We also find SNP rates within conopeptides are higher compared to the remainder of C. bullatus transcriptome, consistent with the hypothesis that conopeptides are under diversifying selection.

The Terebridae and teretoxins: Combining phylogeny and anatomy for concerted discovery of bioactive compounds

The Conoidea superfamily, comprised of cone snails, terebrids, and turrids, is an exceptionally promising group for the discovery of natural peptide toxins. The potential of conoidean toxins has been realized with the distribution of the first Conus (cone snail) drug, Prialt (ziconotide), an analgesic used to alleviate chronic pain in HIV and cancer patients. Cone snail toxins (conotoxins) are highly variable, a consequence of a high mutation rate associated to duplication events and positive selection. As Conus and terebrids diverged in the early Paleocene, the toxins from terebrids (teretoxins) may demonstrate highly divergent and unique functionalities. Recent analyses of the Terebridae, a largely distributed family with more than 300 described species, indicate they have evolutionary and pharmacological potential. Based on a three gene (COI, 12S and 16S) molecular phylogeny, including ~50 species from the West-Pacific, five main terebrid lineages were discriminated: two of these lineages independently lost their venom apparatus, and one venomous lineage was previously unknown. Knowing the phylogenetic relationships within the Terebridae aids in effectively targeting divergent lineages with novel peptide toxins. Preliminary results indicate that teretoxins are similar in structure and composition to conotoxins, suggesting teretoxins are an attractive line of research to discover and develop new therapeutics that target ion channels and receptors. Using conotoxins as a guideline, and innovative natural products discovery strategies, such as the Concerted Discovery Strategy, the potential of the Terebridae and their toxins are explored as a pioneering pharmacological resource.