Systematic interrogation of the Conus marmoreus venom duct transcriptome with ConoSorter reveals 158 novel conotoxins and 13 new gene superfamilies

Coordinated adaptations define the ontogenetic shift from worm- to fish-hunting in a venomous cone snail

Marine cone snails have attracted researchers from all disciplines but early life stages have received limited attention due to difficulties accessing or rearing juvenile specimens. Here, we document the culture of Conus magus from eggs through metamorphosis to reveal dramatic shifts in predatory feeding behaviour between post-metamorphic juveniles and adult specimens. Adult C. magus capture fish using a set of paralytic venom peptides combined with a hooked radular tooth used to tether envenomed fish. In contrast, early juveniles feed exclusively on polychaete worms using a unique "sting-and-stalk" foraging behaviour facilitated by short, unbarbed radular teeth and a distinct venom repertoire that induces hypoactivity in prey. Our results demonstrate how coordinated morphological, behavioural and molecular changes facilitate the shift from worm- to fish-hunting in C. magus, and showcase juvenile cone snails as a rich and unexplored source of novel venom peptides for ecological, evolutionary and biodiscovery studies.

High-Throughput Prediction and Design of Novel Conopeptides for Biomedical Research and Development

Cone snail venoms have been considered a valuable treasure for international scientists and businessmen, mainly due to their pharmacological applications in development of marine drugs for treatment of various human diseases. To date, around 800 Conus species are recorded, and each of them produces over 1,000 venom peptides (termed as conopeptides or conotoxins). This reflects the high diversity and complexity of cone snails, although most of their venoms are still uncharacterized. Advanced multiomics (such as genomics, transcriptomics, and proteomics) approaches have been recently developed to mine diverse Conus venom samples, with the main aim to predict and identify potentially interesting conopeptides in an efficient way. Some bioinformatics techniques have been applied to predict and design novel conopeptide sequences, related targets, and their binding modes. This review provides an overview of current knowledge on the high diversity of conopeptides and multiomics advances in high-throughput prediction of novel conopeptide sequences, as well as molecular modeling and design of potential drugs based on the predicted or validated interactions between these toxins and their molecular targets.

A 4/8 Subtype α-Conotoxin Vt1.27 Inhibits N-Type Calcium Channels With Potent Anti-Allodynic Effect

A novel 4/8 subtype α-conotoxin, Vt1.27 (NCCMFHTCPIDYSRFNC-NH₂), was identified from Conus vitulinus in the South China Sea by RACE methods. The peptide was synthesized and structurally characterized. Similar to other α-conotoxins that target neuronal nicotinic acetylcholine receptor (nAChR) subtypes, Vt1.27 inhibited the rat α3β2 nAChR subtype (IC₅₀ = 1160 nM) and was inactive at voltage-gated sodium and potassium channels in rat sensory neurons. However, Vt1.27 inhibited high voltage-activated N-type (Ca_V2.2) calcium channels expressed in HEK293T cells with an IC₅₀ of 398 nM. An alanine scan of the peptide showed that residues Phe⁵, Pro⁹, Ile¹⁰, and Ser¹³ contribute significantly to the inhibitory activity of Vt1.27. The molecular dockings indicate that Vt1.27 inhibits the transmembrane region of Ca_V2.2, which is different from that of ω-conotoxins. Furthermore, Vt1.27 exhibited potent anti-allodynic effect in rat partial sciatic nerve injury (PNL) and chronic constriction injury (CCI) pain models at 10 nmol/kg level with the intramuscular injection. The pain threshold elevation of Vt1.27 groups was higher than that of α-conotoxin Vc1.1 in CCI rat models. These findings expand our knowledge of targets of α-conotoxins and potentially provide a potent, anti-allodynic peptide for the treatment of neuropathic pain.

Venomics Reveals a Non-Compartmentalised Venom Gland in the Early Diverged Vermivorous Conus distans

The defensive use of cone snail venom is hypothesised to have first arisen in ancestral worm-hunting snails and later repurposed in a compartmentalised venom duct to facilitate the dietary shift to molluscivory and piscivory. Consistent with its placement in a basal lineage, we demonstrate that the C. distans venom gland lacked distinct compartmentalisation. Transcriptomics revealed C. distans expressed a wide range of structural classes, with inhibitory cysteine knot (ICK)-containing peptides dominating. To better understand the evolution of the venom gland compartmentalisation, we compared C. distans to C. planorbis, the earliest diverging species from which a defence-evoked venom has been obtained, and fish-hunting C. geographus from the Gastridium subgenus that injects distinct defensive and predatory venoms. These comparisons support the hypothesis that venom gland compartmentalisation arose in worm-hunting species and enabled repurposing of venom peptides to facilitate the dietary shift from vermivory to molluscivory and piscivory in more recently diverged cone snail lineages.

Comparative Venomics of C. flavidus and C. frigidus and Closely Related Vermivorous Cone Snails

Cone snail venom biodiversity reflects dietary preference and predatory and defensive envenomation strategies across the ≈900 species of Conidae. To better understand the mechanisms of adaptive radiations in closely related species, we investigated the venom of two phylogenetically and spatially related species, C. flavidus and C. frigidus of the Virgiconus clade. Transcriptomic analysis revealed that the major superfamily profiles were conserved between the two species, including 68 shared conotoxin transcripts. These shared transcripts contributed 90% of the conotoxin expression in C. frigidus and only 49% in C. flavidus, which showed greater toxin diversification in the dominant O1, I2, A, O2, O3, and M superfamilies compared to C. frigidus. On the basis of morphology, two additional sub-groups closely resembling C. flavidus were also identified from One Tree Island Reef. Despite the morphological resemblance, the venom duct proteomes of these cryptic sub-groups were distinct from C. flavidus. We suggest rapid conotoxin sequence divergence may have facilitated adaptive radiation and the establishment of new species and the regulatory mechanisms facilitating species-specific venom evolution.

Isolation and characterization of five novel disulfide-poor conopeptides from Conus marmoreus venom

Background:

Conopeptides from cone snail venom have aroused great interest related to the discovery of novel bioactive candidates, due to their excellent prospects for the treatment of various health problems such as pain, addiction, psychosis and epilepsy. In order to explore novel biopeptides, we investigated the structure and function of five novel conopeptides isolated from the venom of Conus marmoreus from South China Sea.

Methods:

C. marmoreus crude venom was prepared, fractionated and purified by HPLC system. The primary sequences of the five novel disulfide-poor conopeptides Mr-1 to Mr-5 were identified by comprehensive analysis of de novo MALDI-TOF tandem mass spectrometry and Edman degradation data. In order to investigate their function, these five conopeptides were synthesized by Fmoc-SPPS chemistry, and their biological effects at several heterologous rat nicotinic acetylcholine receptor (nAChR) subtypes (α1β1δε, α3β2, α3β4, α4β2) were determined by electrophysiological technique.

Results:

Five novel disulfide-poor conopeptides were identified and named as follows: Mr-1 (DWEYHAHPKPNSFWT), Mr-2 (YPTRAYPSNKFG), Mr-3 (NVIQAPAQSVAPP NTST), Mr-4 [KENVLNKLKSK(L/I)] and Mr-5 [NAVAAAN(L/I)PG(L/I)V]. None of them contains a disulfide bond. The sequences of conopeptides Mr-2 to Mr-5 do not belong to any category of the known disulfide-poor conopeptides. No significant activity against the above nAChR subtypes were observed for the five conopeptides at 100 µM.

Conclusion:

We purified and structurally characterized five novel disulfide-poor conopeptides from C. marmoreus crude venom and first investigated their nAChR inhibitory effects. This work expanded our knowledge on the structure and function of disulfide-poor conopeptides from this cone snail venom.

Identification of Novel Conotoxin Precursors from the Cone Snail Conus spurius by High-Throughput RNA Sequencing

Marine gastropods of the genus Conus, comprising more than 800 species, have the characteristic of injecting worms and other prey with venom. These conopeptide toxins, highly diverse in structure and action, are highly potent and specific for their molecular targets (ion channels, receptors, and transporters of the prey's nervous system), and thus are important research tools and source for drug discovery. Next-generation sequencing technologies are speeding up the discovery of novel conopeptides in many of these species, but only limited information is available for Conus spurius, which inhabits sandy mud. To search for new precursor conopeptides, we analyzed the transcriptome of the venous ducts of C. spurius and identified 55 putative conotoxins. Seven were selected for further study and confirmed by Sanger sequencing to belong to the M-superfamily (Sr3.M01 and Sr3.M02), A-superfamily (Sr1.A01 and Sr1.A02), O-superfamily (Sr15.O01), and Con-ikot-ikot (Sr21.CII01 and Sr22.CII02). Six of these have never been reported. To our knowledge, this report is the first to use high-throughput RNA sequencing for the study of the diversity of C. spurius conotoxins.

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

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

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.

Improved prediction of conopeptide superfamilies with ConoDictor 2.0

Motivation

Cone snails are among the richest sources of natural peptides with promising pharmacological and therapeutic applications. With the reduced costs of RNAseq, scientists now heavily rely on venom gland transcriptomes for the mining of novel bioactive conopeptides, but the bioinformatic analyses often hamper the discovery process.

Results

Here, we present ConoDictor 2.0 as a standalone and user-friendly command-line program. We have updated the program originally published as a web server 10 years ago using novel and updated tools and algorithms and improved our classification models with new and higher quality sequences. ConoDictor 2.0 is now more accurate, faster, multiplatform and able to deal with a whole cone snail venom gland transcriptome (raw reads or contigs) in a very short time. The new version of Conodictor also improves the identification and subsequent classification for entirely novel or relatively distant conopeptides. We conducted various tests on known conopeptides from public databases and on the published venom duct transcriptome of Conus geographus, and compared previous results with the output of ConoDictor 2.0, ConoSorter and BLAST. Overall, ConoDictor 2.0 is 4 to 8 times faster for the analysis of a whole transcriptome on a single core computer and performed better at predicting gene superfamily.

Availability and implementation

ConoDictor 2.0 is available as a python 3 git folder at https://github.com/koualab/conodictor.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Combined Proteotranscriptomic-Based Strategy to Discover Novel Antimicrobial Peptides from Cone Snails

Despite their impressive diversity and already broad therapeutic applications, cone snail venoms have received less attention as a natural source in the investigation of antimicrobial peptides than other venomous animals such as scorpions, spiders, or snakes. Cone snails are among the largest genera (Conus sp.) of marine invertebrates, with more than seven hundred species described to date. These predatory mollusks use their sophisticated venom apparatus to capture prey or defend themselves. In-depth studies of these venoms have unraveled many biologically active peptides with pharmacological properties of interest in the field of pain management, the treatment of epilepsy, neurodegenerative diseases, and cardiac ischemia. Considering sequencing efficiency and affordability, cone snail venom gland transcriptome analyses could allow the discovery of new, promising antimicrobial peptides. We first present here the need for novel compounds like antimicrobial peptides as a viable alternative to conventional antibiotics. Secondly, we review the current knowledge on cone snails as a source of antimicrobial peptides. Then, we present the current state of the art in analytical methods applied to crude or milked venom followed by how antibacterial activity assay can be implemented for fostering cone snail antimicrobial peptides studies. We also propose a new innovative profile Hidden Markov model-based approach to annotate full venom gland transcriptomes and speed up the discovery of potentially active peptides from cone snails.

Insertions and Deletions Play an Important Role in the Diversity of Conotoxins

Previous studies have indicated that each conotoxin precursor has a hyperconserved signal region, a rather conserved pro region and a hypervariable mature region, and nucleotide mutations are the main driven factor. However, in this study, we made an in-depth analysis on the M-superfamily conotoxin precursors and found that the diversity of the signal, pro and mature regions are more complicated than previous findings. Different conotoxin precursors can have same signal, pro and/or mature regions, especially different conotoxin precursors with same mature region but different signal and pro regions. In addition, insertions and deletions (indels) were detected in conotoxin precursors. Indels are infrequent in the signal region but frequent in the pro and mature regions. In contrast to deletions that dominate in the pro region, insertions dominate in the mature region. The number of amino acids is crucial for the physiological functions of mature conotoxins, therefore indels, especially insertions in the mature region, play an important role in the sequence and function diversity of conotoxins.

Neuropeptide signalling systems - An underexplored target for venom drug discovery

Neuropeptides are signalling molecules mainly secreted from neurons that act as neurotransmitters or peptide hormones to affect physiological processes and modulate behaviours. In humans, neuropeptides are implicated in numerous diseases and understanding their role in physiological processes and pathologies is important for therapeutic development. Teasing apart the (patho)physiology of neuropeptides remains difficult due to ligand and receptor promiscuity and the complexity of the signalling pathways. The current approach relies on a pharmacological toolbox of agonists and antagonists displaying high selectivity for independent receptor subtypes, with the caveat that only few selective ligands have been discovered or developed. Animal venoms represent an underexplored source for novel receptor subtype-selective ligands that could aid in dissecting human neuropeptide signalling systems. Multiple endogenous-like neuropeptides as well as peptides acting on neuropeptide receptors are present in venoms. In this review, we summarise current knowledge on neuropeptides and discuss venoms as a source for ligands targeting neuropeptide signalling systems.

Diversity of Conopeptides and Their Precursor Genes of Conus Litteratus

The venom of various Conus species is composed of a rich variety of unique bioactive peptides, commonly referred to as conotoxins (conopeptides). Most conopeptides have specific receptors or ion channels as physiologically relevant targets. In this paper, high-throughput transcriptome sequencing was performed to analyze putative conotoxin transcripts from the venom duct of a vermivorous cone snail species, Conus litteratus native to the South China Sea. A total of 128 putative conotoxins were identified, most of them belonging to 22 known superfamilies, with 43 conotoxins being regarded as belonging to new superfamilies. Notably, the M superfamily was the most abundant in conotoxins among the known superfamilies. A total of 15 known cysteine frameworks were also described. The largest proportion of cysteine frameworks were VI/VII (C-C-CC-C-C), IX (C-C-C-C-C-C) and XIV (C-C-C-C). In addition, five novel cysteine patterns were also discovered. Simple sequence repeat detection results showed that di-nucleotide was the major type of repetition, and the codon usage bias results indicated that the codon usage bias of the conotoxin genes was weak, but the M, O1, O2 superfamilies differed in codon preference. Gene cloning indicated that there was no intron in conotoxins of the B1- or J superfamily, one intron with 1273-1339 bp existed in a mature region of the F superfamily, which is different from the previously reported gene structure of conotoxins from other superfamilies. This study will enhance our understanding of conotoxin diversity, and the new conotoxins discovered in this paper will provide more potential candidates for the development of pharmacological probes and marine peptide drugs.

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.

Conotoxins: Chemistry and Biology

The venom of the marine predatory cone snails (genus Conus) has evolved for prey capture and defense, providing the basis for survival and rapid diversification of the now estimated 750+ species. A typical Conus venom contains hundreds to thousands of bioactive peptides known as conotoxins. These mostly disulfide-rich and well-structured peptides act on a wide range of targets such as ion channels, G protein-coupled receptors, transporters, and enzymes. Conotoxins are of interest to neuroscientists as well as drug developers due to their exquisite potency and selectivity, not just against prey but also mammalian targets, thereby providing a rich source of molecular probes and therapeutic leads. The rise of integrated venomics has accelerated conotoxin discovery with now well over 10,000 conotoxin sequences published. However, their structural and pharmacological characterization lags considerably behind. In this review, we highlight the diversity of new conotoxins uncovered since 2014, their three-dimensional structures and folds, novel chemical approaches to their syntheses, and their value as pharmacological tools to unravel complex biology. Additionally, we discuss challenges and future directions for the field.

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.

Transcriptomic Analysis of Marine Gastropod Hemifusus tuba Provides Novel Insights into Conotoxin Genes

The marine gastropod Hemifusus tuba is served as a luxury food in Asian countries and used in traditional Chinese medicine to treat lumbago and deafness. The lack of genomic data on H. tuba is a barrier to aquaculture development and functional characteristics of potential bioactive molecules are poorly understood. In the present study, we used high-throughput sequencing technologies to generate the first transcriptomic database of H. tuba. A total of 41 unique conopeptides were retrieved from 44 unigenes, containing 6-cysteine frameworks belonging to four superfamilies. Duplication of mature regions and alternative splicing were also found in some of the conopeptides, and the de novo assembly identified a total of 76,306 transcripts with an average length of 824.6 nt, of which including 75,620 (99.1%) were annotated. In addition, simple sequence repeats (SSRs) detection identified 14,000 unigenes containing 20,735 SSRs, among which, 23 polymorphic SSRs were screened. Thirteen of these markers could be amplified in Hemifusus ternatanus and seven in Rapana venosa. This study provides reports of conopeptide genes in Buccinidae for the first time as well as genomic resources for further drug development, gene discovery and population resource studies of this species.

Transcriptomic-Proteomic Correlation in the Predation-Evoked Venom of the Cone Snail, Conus imperialis

Individual variation in animal venom has been linked to geographical location, feeding habit, season, size, and gender. Uniquely, cone snails possess the remarkable ability to change venom composition in response to predatory or defensive stimuli. To date, correlations between the venom gland transcriptome and proteome within and between individual cone snails have not been reported. In this study, we use 454 pyrosequencing and mass spectrometry to decipher the transcriptomes and proteomes of the venom gland and corresponding predation-evoked venom of two specimens of Conus imperialis. Transcriptomic analyses revealed 17 conotoxin gene superfamilies common to both animals, including 5 novel superfamilies and two novel cysteine frameworks. While highly expressed transcripts were common to both specimens, variation of moderately and weakly expressed precursor sequences was surprisingly diverse, with one specimen expressing two unique gene superfamilies and consistently producing more paralogs within each conotoxin gene superfamily. Using a quantitative labelling method, conotoxin variability was compared quantitatively, with highly expressed peptides showing a strong correlation between transcription and translation, whereas peptides expressed at lower levels showed a poor correlation. These results suggest that major transcripts are subject to stabilizing selection, while minor transcripts are subject to diversifying selection.

'Messy' Processing of χ-conotoxin MrIA Generates Homologues with Reduced hNET Potency

Integrated venomics techniques have shown that variable processing of conotoxins from Conus marmoreus resulted in a dramatic expansion in the number of expressed conotoxins. One conotoxin from C. marmoreus, the χ-conotoxin MrIA, is a selective inhibitor of human norepinephrine transporters (hNET) and therefore a drug candidate for attenuating chronic neuropathic pain. It has been found that “messy” processing of the MrIA transcripts results in the expression of MrIA analogs with different truncations of the pro-peptide that contains portions of the MrIA molecule. The aim of this study was to investigate if variable processing of the expressed peptides results in modulation of the existing hNET pharmacology or creates new pharmacologies. To this end, a number of MrIA analogs found in C. marmoreus venom were synthesized and evaluated for their activity at hNET receptors. While several of the analogs exhibited norepinephrine transporter inhibitory activity comparable to that of MrIA, none significantly improved on the potency of conotoxin MrIA, and those analogs with disrupted pharmacophores produced greatly reduced NET inhibition, confirming previous structure-activity relationships seen on χ-class conopeptides. Additionally, analogs were screened for new activities on ion channels using calcium influx assays, although no major new pharmacology was revealed.

Venomics Reveals Venom Complexity of the Piscivorous Cone Snail, Conus tulipa

The piscivorous cone snail Conus tulipa has evolved a net-hunting strategy, akin to the deadly Conus geographus, and is considered the second most dangerous cone snail to humans. Here, we present the first venomics study of C. tulipa venom using integrated transcriptomic and proteomic approaches. Parallel transcriptomic analysis of two C. tulipa specimens revealed striking differences in conopeptide expression levels (2.5-fold) between individuals, identifying 522 and 328 conotoxin precursors from 18 known gene superfamilies. Despite broad overlap at the superfamily level, only 86 precursors (11%) were common to both specimens. Conantokins (NMDA antagonists) from the superfamily B1 dominated the transcriptome and proteome of C. tulipa venom, along with superfamilies B2, A, O1, O3, con-ikot-ikot and conopressins, plus novel putative conotoxins precursors T1.3, T6.2, T6.3, T6.4 and T8.1. Thus, C. tulipa venom comprised both paralytic (putative ion channel modulating α-, ω-, μ-, δ-) and non-paralytic (conantokins, con-ikot-ikots, conopressins) conotoxins. This venomic study confirms the potential for non-paralytic conotoxins to contribute to the net-hunting strategy of C. tulipa.

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.

Discovery of Novel Conotoxin Candidates Using Machine Learning

Cone snails (genus Conus) are venomous marine snails that inject prey with a lethal cocktail of conotoxins, small, secreted, and cysteine-rich peptides. Given the diversity and often high affinity for their molecular targets, consisting of ion channels, receptors or transporters, many conotoxins have become invaluable pharmacological probes, drug leads, and therapeutics. Transcriptome sequencing of Conus venom glands followed by de novo assembly and homology-based toxin identification and annotation is currently the state-of-the-art for discovery of new conotoxins. However, homology-based search techniques, by definition, can only detect novel toxins that are homologous to previously reported conotoxins. To overcome these obstacles for discovery, we have created ConusPipe, a machine learning tool that utilizes prominent chemical characters of conotoxins to predict whether a certain transcript in a Conus transcriptome, which has no otherwise detectable homologs in current reference databases, is a putative conotoxin. By using ConusPipe on RNASeq data of 10 species, we report 5148 new putative conotoxin transcripts that have no homologues in current reference databases. 896 of these were identified by at least three out of four models used. These data significantly expand current publicly available conotoxin datasets and our approach provides a new computational avenue for the discovery of novel toxin families.

Discovery Methodology of Novel Conotoxins from Conus Species

Cone snail venoms provide an ideal resource for neuropharmacological tools and drug candidates discovery, which have become a research hotspot in neuroscience and new drug development. More than 1,000,000 natural peptides are produced by cone snails, but less than 0.1% of the estimated conotoxins has been characterized to date. Hence, the discovery of novel conotoxins from the huge conotoxin resources with high-throughput and sensitive methods becomes a crucial key for the conotoxin-based drug development. In this review, we introduce the discovery methodology of new conotoxins from various Conus species. It focuses on obtaining full N- to C-terminal sequences, regardless of disulfide bond connectivity through crude venom purification, conotoxin precusor gene cloning, venom duct transcriptomics, venom proteomics and multi-omic methods. The protocols, advantages, disadvantages, and developments of different approaches during the last decade are summarized and the promising prospects are discussed as well.

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.

Identifying novel conopepetides from the venom ducts of Conus litteratus through integrating transcriptomics and proteomics

The venom ducts of marine cone snails secrete highly complex mixtures of cysteine-rich active peptides, which are generally known as conotoxins or conopeptides and provide a potential fertile resource for pharmacological neuroscience research and the discovery of new drugs. Previous studies have devoted substantial effort to the identification of novel conopeptides, and the 109 cone snail species have yielded 7000 known conopeptides to date. Here, we used de novo deep transcriptome sequencing analyses combined with traditional Sanger sequencing and liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) to identify 30 distinct conopeptide precursors. Twenty of these were previously reported and the other 10 were novel conopeptide precursors. The study provides the first identification of the Con-ikot-ikot, NSF-bt05, O3 and I1 gene superfamilies in C. litteratus. A new putative superfamily was identified. In addition, the following cysteine frameworks were first identified in this study: CC-C-C-C-C-C-C-C-C-C-C-C-CC-C-C-C-C-C and C-C-C-C-C-CC-C. Several isomerases involved in post-translational modification of conopeptides were identified as well. The discovery of new conopeptides in C. litteratus will enhance our understanding of the conopeptide diversity in this particular clade of cone snails. We also found the existence of intraspecific variations in vermivorous species. Finally, the analysis strategy offers a relatively reliable workflow for screening for peptide drug candidates. SIGNIFICANCE: These novel conopeptides provide a potential resource for the development of new channel-targeting drugs. The intraspecific variation in C. litteratus enhance our understanding of the conopeptide diversity in this particular clade of cone snails. The identified three cysteine residues, which might participate in the formation of disulfide bonds, provide a clue to get the connectivity of cysteine frameworks. Finally, the analysis strategy offers a relatively reliable workflow for screening for peptide drug candidates.

Venomics: A Mini-Review

Venomics is the integration of proteomic, genomic and transcriptomic approaches to study venoms. Advances in these approaches have enabled increasingly more comprehensive analyses of venoms to be carried out, overcoming to some extent the limitations imposed by the complexity of the venoms and the small quantities that are often available. Advances in bioinformatics and high-throughput functional assay screening approaches have also had a significant impact on venomics. A combination of all these techniques is critical for enhancing our knowledge on the complexity of venoms and their potential therapeutic and agricultural applications. Here we highlight recent advances in these fields and their impact on venom analyses.

Venomics-Accelerated Cone Snail Venom Peptide Discovery

Cone snail venoms are considered a treasure trove of bioactive peptides. Despite over 800 species of cone snails being known, each producing over 1000 venom peptides, only about 150 unique venom peptides are structurally and functionally characterized. To overcome the limitations of the traditional low-throughput bio-discovery approaches, multi-omics systems approaches have been introduced to accelerate venom peptide discovery and characterisation. This “venomic” approach is starting to unravel the full complexity of cone snail venoms and to provide new insights into their biology and evolution. The main challenge for venomics is the effective integration of transcriptomics, proteomics, and pharmacological data and the efficient analysis of big datasets. Novel database search tools and visualisation techniques are now being introduced that facilitate data exploration, with ongoing advances in related omics fields being expected to further enhance venomics studies. Despite these challenges and future opportunities, cone snail venomics has already exponentially expanded the number of novel venom peptide sequences identified from the species investigated, although most novel conotoxins remain to be pharmacologically characterised. Therefore, efficient high-throughput peptide production systems and/or banks of miniaturized discovery assays are required to overcome this bottleneck and thus enhance cone snail venom bioprospecting and accelerate the identification of novel drug leads.

Structural plasticity of mini-M conotoxins - expression of all mini-M subtypes by Conus regius

The mini-M conotoxins are peptidic scaffolds found in the venom of cones snails. These scaffolds are tightly folded structures held together by three disulfide bonds with a CC-C-C-CC arrangement (conotoxin framework III) and belong to the M Superfamily of conotoxins. Here, we describe mini-M conotoxins from the venom of Conus regius, a Western Atlantic worm-hunting cone snail species using transcriptomic and peptidomic analyses. These C. regius conotoxins belong to three different subtypes: M1, M2, and M3. The subtypes show little sequence homology, and their loop sizes (intercysteine amino acid chains) vary significantly. The mini-Ms isolated from dissected venom contains preferentially hydroxylated proline residues, thus augmenting the structural reach of this conotoxin class. Using 2D-NMR methods, we have determined the 3D structure of reg3b, an M2 subtype conotoxin, which shows a constrained multi-turn scaffold. The structural diversity found within mini-M conotoxin scaffolds of C. regius is indicative of structural hypervariability of the conotoxin M superfamily that is not seen in other superfamilies. These stable minimalistic scaffolds may be investigated for the development of engineered peptides for therapeutic applications.

DATABASES

Sequences are available in GenBank under accession numbers MF588935-MF588952. Structural data are available in the RCSB protein database under the accession code 6BX9.

Cone Snails: A Big Store of Conotoxins for Novel Drug Discovery

Marine drugs have developed rapidly in recent decades. Cone snails, a group of more than 700 species, have always been one of the focuses for new drug discovery. These venomous snails capture prey using a diverse array of unique bioactive neurotoxins, usually named as conotoxins or conopeptides. These conotoxins have proven to be valuable pharmacological probes and potential drugs due to their high specificity and affinity to ion channels, receptors, and transporters in the nervous systems of target prey and humans. Several research groups, including ours, have examined the venom gland of cone snails using a combination of transcriptomic and proteomic sequencing, and revealed the existence of hundreds of conotoxin transcripts and thousands of conopeptides in each Conus species. Over 2000 nucleotide and 8000 peptide sequences of conotoxins have been published, and the number is still increasing quickly. However, more than 98% of these sequences still lack 3D structural and functional information. With the rapid development of genomics and bioinformatics in recent years, functional predictions and investigations on conotoxins are making great progress in promoting the discovery of novel drugs. For example, ω-MVIIA was approved by the U.S. Food and Drug Administration in 2004 to treat chronic pain, and nine more conotoxins are at various stages of preclinical or clinical evaluation. In short, the genus Conus, the big family of cone snails, has become an important genetic resource for conotoxin identification and drug development.

A Transcriptomic Survey of Ion Channel-Based Conotoxins in the Chinese Tubular Cone Snail (Conus betulinus)

Conotoxins in the venom of cone snails (Conus spp.) are a mixture of active peptides that work as blockers, agonists, antagonists, or inactivators of various ion channels. Recently we reported a high-throughput method to identify 215 conotoxin transcripts from the Chinese tubular cone snail, C. betulinus. Here, based on the previous datasets of four transcriptomes from three venom ducts and one venom bulb, we explored ion channel-based conotoxins and predicted their related ion channel receptors. Homologous analysis was also performed for the most abundant ion channel protein, voltage-gated potassium (Kv; with Kv1.1 as the representative), and the most studied ion channel receptor, nicotinic acetylcholine receptor (nAChR; with α2-nAChR as the representative), in different animals. Our transcriptomic survey demonstrated that ion channel-based conotoxins and related ion channel proteins/receptors transcribe differentially between the venom duct and the venom bulb. In addition, we observed that putative κ-conotoxins were the most common conotoxins with the highest transcription levels in the examined C. betulinus. Furthermore, Kv1.1 and α2-nAChR were conserved in their functional domains of deduced protein sequences, suggesting similar effects of conotoxins via the ion channels in various species, including human beings. In a word, our present work suggests a high-throughput way to develop conotoxins as potential drugs for treatment of ion channel-associated human diseases.

Pharmacological screening technologies for venom peptide discovery

Venomous animals occupy one of the most successful evolutionary niches and occur on nearly every continent. They deliver venoms via biting and stinging apparatuses with the aim to rapidly incapacitate prey and deter predators. This has led to the evolution of venom components that act at a number of biological targets - including ion channels, G-protein coupled receptors, transporters and enzymes - with exquisite selectivity and potency, making venom-derived components attractive pharmacological tool compounds and drug leads. In recent years, plate-based pharmacological screening approaches have been introduced to accelerate venom-derived drug discovery. A range of assays are amenable to this purpose, including high-throughput electrophysiology, fluorescence-based functional and binding assays. However, despite these technological advances, the traditional activity-guided fractionation approach is time-consuming and resource-intensive. The combination of screening techniques suitable for miniaturization with sequence-based discovery approaches - supported by advanced proteomics, mass spectrometry, chromatography as well as synthesis and expression techniques - promises to further improve venom peptide discovery. Here, we discuss practical aspects of establishing a pipeline for venom peptide drug discovery with a particular emphasis on pharmacology and pharmacological screening approaches. This article is part of the Special Issue entitled 'Venom-derived Peptides as Pharmacological Tools.'

From Marine Venoms to Drugs: Efficiently Supported by a Combination of Transcriptomics and Proteomics

The potential of marine natural products to become new drugs is vast; however, research is still in its infancy. The chemical and biological diversity of marine toxins is immeasurable and as such an extraordinary resource for the discovery of new drugs. With the rapid development of next-generation sequencing (NGS) and liquid chromatography-tandem mass spectrometry (LC-MS/MS), it has been much easier and faster to identify more toxins and predict their functions with bioinformatics pipelines, which pave the way for novel drug developments. Here we provide an overview of related bioinformatics pipelines that have been supported by a combination of transcriptomics and proteomics for identification and function prediction of novel marine toxins.

Is Mutation Random or Targeted?: No Evidence for Hypermutability in Snail Toxin Genes

Ever since Luria and Delbruck, the notion that mutation is random with respect to fitness has been foundational to modern biology. However, various studies have claimed striking exceptions to this rule. One influential case involves toxin-encoding genes in snails of the genus Conus, termed conotoxins, a large gene family that undergoes rapid diversification of their protein-coding sequences by positive selection. Previous reconstructions of the sequence evolution of conotoxin genes claimed striking patterns: (1) elevated synonymous change, interpreted as being due to targeted "hypermutation" in this region; (2) elevated transversion-to-transition ratios, interpreted as reflective of the particular mechanism of hypermutation; and (3) much lower rates of synonymous change in the codons encoding several highly conserved cysteine residues, interpreted as strong position-specific codon bias. This work has spawned a variety of studies on the potential mechanisms of hypermutation and on causes for cysteine codon bias, and has inspired hypermutation hypotheses for various other fast-evolving genes. Here, I show that all three findings are likely to be artifacts of statistical reconstruction. First, by simulating nonsynonymous change I show that high rates of dN can lead to overestimation of dS. Second, I show that there is no evidence for any of these three patterns in comparisons of closely related conotoxin sequences, suggesting that the reported findings are due to breakdown of statistical methods at high levels of sequence divergence. The current findings suggest that mutation and codon bias in conotoxin genes may not be atypical, and that random mutation and selection can explain the evolution of even these exceptional loci.

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.

δ-Conotoxin SuVIA suggests an evolutionary link between ancestral predator defence and the origin of fish-hunting behaviour in carnivorous cone snails

Some venomous cone snails feed on small fishes using an immobilizing combination of synergistic venom peptides that target Kv and Nav channels. As part of this envenomation strategy, δ-conotoxins are potent ichtyotoxins that enhance Nav channel function. δ-Conotoxins belong to an ancient and widely distributed gene superfamily, but any evolutionary link from ancestral worm-eating cone snails to modern piscivorous species has not been elucidated. Here, we report the discovery of SuVIA, a potent vertebrate-active δ-conotoxin characterized from a vermivorous cone snail (Conus suturatus). SuVIA is equipotent at hNaV1.3, hNaV1.4 and hNaV1.6 with EC50s in the low nanomolar range. SuVIA also increased peak hNaV1.7 current by approximately 75% and shifted the voltage-dependence of activation to more hyperpolarized potentials from -15 mV to -25 mV, with little effect on the voltage-dependence of inactivation. Interestingly, the proximal venom gland expression and pain-inducing effect of SuVIA in mammals suggest that δ-conotoxins in vermivorous cone snails play a defensive role against higher order vertebrates. We propose that δ-conotoxins originally evolved in ancestral vermivorous cones to defend against larger predators including fishes have been repurposed to facilitate a shift to piscivorous behaviour, suggesting an unexpected underlying mechanism for this remarkable evolutionary transition.

High-throughput identification of novel conotoxins from the Chinese tubular cone snail (Conus betulinus) by multi-transcriptome sequencing

Background

The venom of predatory marine cone snails mainly contains a diverse array of unique bioactive peptides commonly referred to as conopeptides or conotoxins. These peptides have proven to be valuable pharmacological probes and potential drugs because of their high specificity and affinity to important ion channels, receptors and transporters of the nervous system. Most previous studies have focused specifically on the conopeptides from piscivorous and molluscivorous cone snails, but little attention has been devoted to the dominant vermivorous species.

Results

The vermivorous Chinese tubular cone snail, Conus betulinus, is the dominant Conus species inhabiting the South China Sea. The transcriptomes of venom ducts and venom bulbs from a variety of specimens of this species were sequenced using both next-generation sequencing and traditional Sanger sequencing technologies, resulting in the identification of a total of 215 distinct conopeptides. Among these, 183 were novel conopeptides, including nine new superfamilies. It appeared that most of the identified conopeptides were synthesized in the venom duct, while a handful of conopeptides were identified only in the venom bulb and at very low levels.

Conclusions

We identified 215 unique putative conopeptide transcripts from the combination of five transcriptomes and one EST sequencing dataset. Variation in conopeptides from different specimens of C. betulinus was observed, which suggested the presence of intraspecific variability in toxin production at the genetic level. These novel conopeptides provide a potentially fertile resource for the development of new pharmaceuticals, and a pathway for the discovery of new conotoxins.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-016-0122-9) contains supplementary material, which is available to authorized users.

The role of defensive ecological interactions in the evolution of conotoxins

Venoms comprise of complex mixtures of peptides evolved for predation and defensive purposes. Remarkably, some carnivorous cone snails can inject two distinct venoms in response to predatory or defensive stimuli, providing a unique opportunity to study separately how different ecological pressures contribute to toxin diversification. Here, we report the extraordinary defensive strategy of the Rhizoconus subgenus of cone snails. The defensive venom from this worm-hunting subgenus is unusually simple, almost exclusively composed of αD-conotoxins instead of the ubiquitous αA-conotoxins found in the more complex defensive venom of mollusc- and fish-hunting cone snails. A similarly compartmentalized venom gland as those observed in the other dietary groups facilitates the deployment of this defensive venom. Transcriptomic analysis of a Conus vexillum venom gland revealed the αD-conotoxins as the major transcripts, with lower amounts of 15 known and four new conotoxin superfamilies also detected with likely roles in prey capture. Our phylogenetic and molecular evolution analysis of the αD-conotoxins from five subgenera of cone snails suggests they evolved episodically as part of a defensive strategy in the Rhizoconus subgenus. Thus, our results demonstrate an important role for defence in the evolution of conotoxins.

Transcriptome and proteome of Conus planorbis identify the nicotinic receptors as primary target for the defensive venom

Most venomous predators have evolved complex venom primarily to immobilize their prey and secondarily to defend against predators. In a new paradigm, carnivorous marine gastropods of the genus Conus were shown to rapidly and reversibly switch between two types of venoms in response to predatory or defensive stimulus, suggesting that the defensive use of venom may have a more important role in venom evolution and specialization than previously thought. To further investigate this phenomenon, the defensive repertoire of a vermivorous species, Conus planorbis, was deciphered using second-generation sequencing coupled to high-throughput proteomics. The venom gland transcriptome of C. planorbis revealed 182 unique conotoxin precursors from 25 gene superfamilies, with superfamily T dominating in terms of read and paralog numbers. Analysis of the defense-evoked venom revealed that this vermivorous species uses a similarly complex arsenal to deter aggressors as more recently evolved fish- and mollusk-hunting species, with MS/MS validating 23 conotoxin sequences from six superfamilies. Pharmacological characterization of the defensive venom on human receptors identified the nicotinic acetylcholine receptors as a primary target. This work provides the first insights into the composition and biological activity of specifically evolved defensive venoms in vermivorous cone snails.

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.

Cone snail venomics: from novel biology to novel therapeutics

Peptide neurotoxins from cone snails called conotoxins are renowned for their therapeutic potential to treat pain and several neurodegenerative diseases. Inefficient assay-guided discovery methods have been replaced by high-throughput bioassays integrated with advanced MS and next-generation sequencing, ushering in the era of 'venomics'. In this review, we focus on the impact of venomics on the understanding of cone snail biology as well as the application of venomics to accelerate the discovery of new conotoxins. We also discuss the continued importance of medicinal chemistry approaches to optimize conotoxins for clinical use, with a descriptive case study of MrIA featured.

Optimized deep-targeted proteotranscriptomic profiling reveals unexplored Conus toxin diversity and novel cysteine frameworks

Cone snails are predatory marine gastropods characterized by a sophisticated venom apparatus responsible for the biosynthesis and delivery of complex mixtures of cysteine-rich toxin peptides. These conotoxins fold into small highly structured frameworks, allowing them to potently and selectively interact with heterologous ion channels and receptors. Approximately 2,000 toxins from an estimated number of >70,000 bioactive peptides have been identified in the genus Conus to date. Here, we describe a high-resolution interrogation of the transcriptomes (available at www.ddbj.nig.ac.jp) and proteomes of the diverse compartments of the Conus episcopatus venom apparatus. Using biochemical and bioinformatic tools, we found the highest number of conopeptides yet discovered in a single Conus specimen, with 3,305 novel precursor toxin sequences classified into 9 known superfamilies (A, I1, I2, M, O1, O2, S, T, Z), and identified 16 new superfamilies showing unique signal peptide signatures. We were also able to depict the largest population of venom peptides containing the pharmacologically active C-C-CC-C-C inhibitor cystine knot and CC-C-C motifs (168 and 44 toxins, respectively), as well as 208 new conotoxins displaying odd numbers of cysteine residues derived from known conotoxin motifs. Importantly, six novel cysteine-rich frameworks were revealed which may have novel pharmacology. Finally, analyses of codon usage bias and RNA-editing processes of the conotoxin transcripts demonstrate a specific conservation of the cysteine skeleton at the nucleic acid level and provide new insights about the origin of sequence hypervariablity in mature toxin regions.

Bioinformatics-Aided Venomics

Venomics is a modern approach that combines transcriptomics and proteomics to explore the toxin content of venoms. This review will give an overview of computational approaches that have been created to classify and consolidate venomics data, as well as algorithms that have helped discovery and analysis of toxin nucleic acid and protein sequences, toxin three-dimensional structures and toxin functions. Bioinformatics is used to tackle specific challenges associated with the identification and annotations of toxins. Recognizing toxin transcript sequences among second generation sequencing data cannot rely only on basic sequence similarity because toxins are highly divergent. Mass spectrometry sequencing of mature toxins is challenging because toxins can display a large number of post-translational modifications. Identifying the mature toxin region in toxin precursor sequences requires the prediction of the cleavage sites of proprotein convertases, most of which are unknown or not well characterized. Tracing the evolutionary relationships between toxins should consider specific mechanisms of rapid evolution as well as interactions between predatory animals and prey. Rapidly determining the activity of toxins is the main bottleneck in venomics discovery, but some recent bioinformatics and molecular modeling approaches give hope that accurate predictions of toxin specificity could be made in the near future.

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.

High conopeptide diversity in Conus tribblei revealed through analysis of venom duct transcriptome using two high-throughput sequencing platforms

The venom of each species of Conus contains different kinds of pharmacologically active peptides which are mostly unique to that species. Collectively, the ~500-700 species of Conus produce a large number of these peptides, perhaps exceeding 140,000 different types in total. To date, however, only a small fraction of this diversity has been characterized via transcriptome sequencing. In addition, the sampling of this chemical diversity has not been uniform across the different lineages in the genus. In this study, we used high-throughput transcriptome sequencing approach to further investigate the diversity of Conus venom peptides. We chose a species, Conus tribblei, as a representative of a poorly studied clade of Conus. Using the Roche 454 and Illumina platforms, we discovered 136 unique and novel putative conopeptides belonging to 30 known gene superfamilies and 6 new conopeptide groups, the greatest diversity so far observed from a transcriptome. Most of the identified peptides exhibited divergence from the known conopeptides, and some contained cysteine frameworks observed for the first time in cone snails. In addition, several enzymes involved in posttranslational modification of conopeptides and also some proteins involved in efficient delivery of the conopeptides to prey were identified as well. Interestingly, a number of conopeptides highly similar to the conopeptides identified in a phylogenetically distant species, the generalist feeder Conus californicus, were observed. The high diversity of conopeptides and the presence of conopeptides similar to those in C. californicus suggest that C. tribblei may have a broad range of prey preferences.