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

Proteomic analysis of the venom of Conus flavidus from Red Sea reveals potential pharmacological applications

BACKGROUND

Venomous marine cone snails produce unique neurotoxins called conopeptides or conotoxins, which are valuable for research and drug discovery. Characterizing Conus venom is important, especially for poorly studied species, as these tiny and steady molecules have considerable potential as research tools for detecting new pharmacological applications. In this study, a worm-hunting cone snail, Conus flavidus inhabiting the Red Sea coast were collected, dissected and the venom gland extraction was subjected to proteomic analysis to define the venom composition, and confirm the functional structure of conopeptides.

RESULTS

Analysis of C. flavidus venom identified 117 peptide fragments and assorted them to conotoxin precursors and non-conotoxin proteins. In this procedure, 65 conotoxin precursors were classified and identified to 16 conotoxin precursors and hormone superfamilies. In the venom of C. flavidus, the four conotoxin superfamilies T, A, O2, and M were the most abundant peptides, accounting for 75.8% of the total conotoxin diversity. Additionally, 19 non-conotoxin proteins were specified in the venom, as well as several potentially biologically active peptides with putative applications.

CONCLUSION

Our research displayed that the structure of the C. flavidus-derived proteome is similar to other Conus species and includes toxins, ionic channel inhibitors, insulin-like peptides, and hyaluronidase. This study provides a foundation for discovering new conopeptides from C. flavidus venom for pharmaceutical use.

Diversity and Evolutionary Analysis of Venom Insulin Derived from Cone Snails

Cone snails possess a diverse array of novel peptide toxins, which selectively target ion channels and receptors in the nervous and cardiovascular systems. These numerous novel peptide toxins are a valuable resource for future marine drug development. In this review, we compared and analyzed the sequence diversity, three-dimensional structural variations, and evolutionary aspects of venom insulin derived from different cone snail species. The comparative analysis reveals that there are significant variations in the sequences and three-dimensional structures of venom insulins from cone snails with different feeding habits. Notably, the venom insulin of some piscivorous cone snails exhibits a greater similarity to humans and zebrafish insulins. It is important to emphasize that these venom insulins play a crucial role in the predatory strategies of these cone snails. Furthermore, a phylogenetic tree was constructed to trace the lineage of venom insulin sequences, shedding light on the evolutionary interconnections among cone snails with diverse diets.

Biomaterials and Bioactive Natural Products from Marine Invertebrates: From Basic Research to Innovative Applications

Aquatic invertebrates are a major source of biomaterials and bioactive natural products that can find applications as pharmaceutics, nutraceutics, cosmetics, antibiotics, antifouling products and biomaterials. Symbiotic microorganisms are often the real producers of many secondary metabolites initially isolated from marine invertebrates; however, a certain number of them are actually synthesized by the macro-organisms. In this review, we analysed the literature of the years 2010–2019 on natural products (bioactive molecules and biomaterials) from the main phyla of marine invertebrates explored so far, including sponges, cnidarians, molluscs, echinoderms and ascidians, and present relevant examples of natural products of interest to public and private stakeholders. We also describe omics tools that have been more relevant in identifying and understanding mechanisms and processes underlying the biosynthesis of secondary metabolites in marine invertebrates. Since there is increasing attention on finding new solutions for a sustainable large-scale supply of bioactive compounds, we propose that a possible improvement in the biodiscovery pipeline might also come from the study and utilization of aquatic invertebrate stem cells.

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.

Detection of αB-Conotoxin VxXXIVA (αB-CTX) by ic-ELISA Based on an Epitope-Specific Monoclonal Antibody

In view of the toxicological hazard and important applications in analgesics and cancer chemotherapeutics of αB-CTX, it is urgent to develop an accurate, effective and feasible immunoassay for the determination and analysis of αB-CTX in real samples. In this study, MBP-αB-CTX4 tandem fusion protein was used as an immunogen to elicit a strong immune response, and a hybridoma cell 5E4 secreting IgG2b against αB-CTX was successfully screened by hybridoma technology. The affinity of the purified 5E4 monoclonal antibody (mAb) was 1.02 × 108 L/mol, which showed high affinity and specificity to αB-CTX. Epitope 1 of αB-CTX is the major binding region for 5E4 mAb recongnization, and two amino acid residues (14L and 15F) in αB-CTX were critical sites for the interaction between αB-CTX and 5E4 mAb. Indirect competitive ELISA (ic-ELISA) based on 5E4 mAb was developed to detect and analyze αB-CTX in real samples, and the linear range of ic-ELISA to αB-CTX was 117–3798 ng/mL, with a limit of detection (LOD) of 81 ng/mL. All the above results indicated that the developed ic-ELISA had high accuracy and repeatability, and it could be applied for αB-CTX detection and drug analysis in real samples.

Proteomic analysis of Red Sea Conus taeniatus venom reveals potential biological applications

Background:

Diverse and unique bioactive neurotoxins known as conopeptides or conotoxins are produced by venomous marine cone snails. Currently, these small and stable molecules are of great importance as research tools and platforms for discovering new drugs and therapeutics. Therefore, the characterization of Conus venom is of great significance, especially for poorly studied species.

Methods:

In this study, we used bioanalytical techniques to determine the venom profile and emphasize the functional composition of conopeptides in Conus taeniatus, a neglected worm-hunting cone snail.

Results:

The proteomic analysis revealed that 84.0% of the venom proteins were between 500 and 4,000 Da, and 16.0% were > 4,000 Da. In C. taeniatus venom, 234 peptide fragments were identified and classified as conotoxin precursors or non-conotoxin proteins. In this process, 153 conotoxin precursors were identified and matched to 23 conotoxin precursors and hormone superfamilies. Notably, the four conotoxin superfamilies T (22.87%), O1 (17.65%), M (13.1%) and O2 (9.8%) were the most abundant peptides in C. taeniatus venom, accounting for 63.40% of the total conotoxin diversity. On the other hand, 48 non-conotoxin proteins were identified in the venom of C. taeniatus. Moreover, several possibly biologically active peptide matches were identified, and putative applications of the peptides were assigned.

Conclusion:

Our study showed that the composition of the C. taeniatus-derived proteome is comparable to that of other Conus species and contains an effective mix of toxins, ionic channel inhibitors and antimicrobials. Additionally, it provides a guidepost for identifying novel conopeptides from the venom of C. taeniatus and discovering conopeptides of potential pharmaceutical importance.

Diversity of Conopeptides and Conoenzymes from the Venom Duct of the Marine Cone Snail Conus bayani as Determined from Transcriptomic and Proteomic Analyses

Marine cone snails are predatory gastropods characterized by a well-developed venom apparatus and highly evolved hunting strategies that utilize toxins to paralyze prey and defend against predators. The venom of each species of cone snail has a large number of pharmacologically active peptides known as conopeptides or conotoxins that are usually unique in each species. Nevertheless, venoms of only very few species have been characterized so far by transcriptomic approaches. In this study, we used transcriptome sequencing technologies and mass spectrometric methods to describe the diversity of venom components expressed by a worm-hunting species, Conus bayani. A total of 82 conotoxin sequences were retrieved from transcriptomic data that contain 54 validated conotoxin sequences clustered into 21 gene superfamilies including divergent gene family, 17 sequences clustered to 6 different conotoxin classes, and 11 conotoxins classified as unassigned gene family. Seven new conotoxin sequences showed unusual cysteine patterns. We were also able to identify 19 peptide sequences using mass spectrometry that completely overlapped with the conotoxin sequences obtained from transcriptome analysis. Importantly, herein we document the presence of 16 proteins that include five post-translational modifying enzymes obtained from transcriptomic data. Our results revealed diverse and novel conopeptides of an unexplored species that could be used extensively in biomedical research due to their therapeutic potentials.

The first Conus genome assembly reveals a primary genetic central dogma of conopeptides in C. betulinus

Although there are various Conus species with publicly available transcriptome and proteome data, no genome assembly has been reported yet. Here, using Chinese tubular cone snail (C. betulinus) as a representative, we sequenced and assembled the first Conus genome with original identification of 133 genome-widely distributed conopeptide genes. After integration of our genomics, transcriptomics, and peptidomics data in the same species, we established a primary genetic central dogma of diverse conopeptides, assuming a rough number ratio of ~1:1:1:10s for the total genes: transcripts: proteins: post-translationally modified peptides. This ratio may be special for this worm-hunting Conus species, due to the high diversity of various Conus genomes and the big number ranges of conopeptide genes, transcripts, and peptides in previous reports of diverse Conus species. Only a fraction (45.9%) of the identified conotopeptide genes from our achieved genome assembly are transcribed with transcriptomic evidence, and few genes individually correspond to multiple transcripts possibly due to intraspecies or mutation-based variances. Variable peptide processing at the proteomic level, generating a big diversity of venom conopeptides with alternative cleavage sites, post-translational modifications, and N-/C-terminal truncations, may explain how the 133 genes and ~123 transcripts can generate thousands of conopeptides in the venom of individual C. betulinus. We also predicted many conopeptides with high stereostructural similarities to the putative analgesic ω-MVIIA, addiction therapy AuIB and insecticide ImI, suggesting that our current genome assembly for C. betulinus is a valuable genetic resource for high-throughput prediction and development of potential pharmaceuticals.

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.

Multi-omics investigations within the Phylum Mollusca, Class Gastropoda: from ecological application to breakthrough phylogenomic studies

Gastropods are the largest and most diverse class of mollusc and include species that are well studied within the areas of taxonomy, aquaculture, biomineralization, ecology, microbiome and health. Gastropod research has been expanding since the mid-2000s, largely due to large-scale data integration from next-generation sequencing and mass spectrometry in which transcripts, proteins and metabolites can be readily explored systematically. Correspondingly, the huge data added a great deal of complexity for data organization, visualization and interpretation. Here, we reviewed the recent advances involving gastropod omics ('gastropodomics') research from hundreds of publications and online genomics databases. By summarizing the current publicly available data, we present an insight for the design of useful data integrating tools and strategies for comparative omics studies in the future. Additionally, we discuss the future of omics applications in aquaculture, natural pharmaceutical biodiscovery and pest management, as well as to monitor the impact of environmental stressors.

Transcriptomic Profiling Reveals Extraordinary Diversity of Venom Peptides in Unexplored Predatory Gastropods of the Genus Clavus

Predatory gastropods of the superfamily Conoidea number over 12,000 living species. The evolutionary success of this lineage can be explained by the ability of conoideans to produce complex venoms for hunting, defense, and competitive interactions. Whereas venoms of cone snails (family Conidae) have become increasingly well studied, the venoms of most other conoidean lineages remain largely uncharacterized. In the present study, we present the venom gland transcriptomes of two species of the genus Clavus that belong to the family Drilliidae. Venom gland transcriptomes of two specimens of Clavus canalicularis and two specimens of Clavus davidgilmouri were analyzed, leading to the identification of a total of 1,176 putative venom peptide toxins (drillipeptides). Based on the combined evidence of secretion signal sequence identity, entire precursor similarity search (BLAST), and the orthology inference, putative Clavus toxins were assigned to 158 different gene families. The majority of identified transcripts comprise signal, pro-, mature peptide, and post-regions, with a typically short (<50 amino acids) and cysteine-rich mature peptide region. Thus, drillipeptides are structurally similar to conotoxins. However, convincing homology with known groups of Conus toxins was only detected for very few toxin families. Among these are Clavus counterparts of Conus venom insulins (drillinsulins), porins (drilliporins), and highly diversified lectins (drillilectins). The short size of most drillipeptides and structural similarity to conotoxins were unexpected, given that most related conoidean gastropod families (Terebridae and Turridae) possess longer mature peptide regions. Our findings indicate that, similar to conotoxins, drillipeptides may represent a valuable resource for future pharmacological exploration.

Venom Diversity and Evolution in the Most Divergent Cone Snail Genus Profundiconus

Profundiconus is the most divergent cone snail genus and its unique phylogenetic position, sister to the rest of the family Conidae, makes it a key taxon for examining venom evolution and diversity. Venom gland and foot transcriptomes of Profundiconus cf. vaubani and Profundiconus neocaledonicus were de novo assembled, annotated, and analyzed for differential expression. One hundred and thirty-seven venom components were identified from P. cf. vaubani and 82 from P. neocaledonicus, with only four shared by both species. The majority of the transcript diversity was composed of putative peptides, including conotoxins, profunditoxins, turripeptides, insulin, and prohormone-4. However, there were also a significant percentage of other putative venom components such as chymotrypsin and L-rhamnose-binding lectin. The large majority of conotoxins appeared to be from new gene superfamilies, three of which are highly different from previously reported venom peptide toxins. Their low conotoxin diversity and the type of insulin found suggested that these species, for which no ecological information are available, have a worm or molluscan diet associated with a narrow dietary breadth. Our results indicate that Profundiconus venom is highly distinct from that of other cone snails, and therefore important for examining venom evolution in the Conidae family.

Spider Venom: Components, Modes of Action, and Novel Strategies in Transcriptomic and Proteomic Analyses

This review gives an overview on the development of research on spider venoms with a focus on structure and function of venom components and techniques of analysis. Major venom component groups are small molecular mass compounds, antimicrobial (also called cytolytic, or cationic) peptides (only in some spider families), cysteine-rich (neurotoxic) peptides, and enzymes and proteins. Cysteine-rich peptides are reviewed with respect to various structural motifs, their targets (ion channels, membrane receptors), nomenclature, and molecular binding. We further describe the latest findings concerning the maturation of antimicrobial, and cysteine-rich peptides that are in most known cases expressed as propeptide-containing precursors. Today, venom research, increasingly employs transcriptomic and mass spectrometric techniques. Pros and cons of venom gland transcriptome analysis with Sanger, 454, and Illumina sequencing are discussed and an overview on so far published transcriptome studies is given. In this respect, we also discuss the only recently described cross contamination arising from multiplexing in Illumina sequencing and its possible impacts on venom studies. High throughput mass spectrometric analysis of venom proteomes (bottom-up, top-down) are reviewed.

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.

Development of an ic-ELISA and immunochromatographic strip based on IgG antibody for detection of ω-conotoxin MVIIA

ω-conotoxin MVIIA(ω-CTX MVIIA) is a peptide consisting of 25 amino acid residues secreted mainly by Conus magus. In view of the toxin threat to humans and animals and defined application in analgesic therapy, it is necessary to develop a rapid, effective and accuracy method for the quantification and analysis of ω-CTX MVIIA in real samples. In the present study, a hybridoma cell named 2E5 stable secreting IgG antibody against ω-CTX MVIIA was selected successfully, and the subtype of Mab 2E5 was IgG1. The purified monoclonal antibody(Mab) 2E5 has high affinity (about 2.79 × 10⁹ L/mol), and shows high specificity to ω-CTX MVIIA antigen. The linear range of ic-ELISA to detect ω-CTX MVIIA was 0.20˜7.22 μg/mL, with a lower detection limit (LOD) of 0.14 ng/mL. The average recovery of intra- and inter-assay were (85.45 ± 2.28)% and (88.03 ± 4.80)% respectively, with a coefficient of variation from 2.59% to 5.42%. The LOD of colloidal strip by naked eye was 1 μg/mL, and the detection time was less than 10 min without any equipment. The developed ELISA and colloidal test strips based on this IgG antibody could be used to detect ω-CTX MVIIA residue in real Conus samples.

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.

Neurobiology and therapeutic utility of neurotoxins targeting postsynaptic mechanisms of neuromuscular transmission

The neuromuscular junction (NMJ) is the principal site for the translation of motor neurochemical signals to muscle activity. Therefore, the release and sensing machinery of acetylcholine (ACh) along with muscle contraction are two of the main targets of natural toxins and pathogens, causing paralysis. Given pharmacology and medical advances, the active ingredients of toxins that target postsynaptic mechanisms have become of major interest, showing promise as drug leads. Herein, we review key facets of prevalent toxins modulating the mechanisms of ACh sensing and generation of the postsynaptic response, with muscle contraction. We consider the correlation between their outstanding selectivity and potency plus effects on motor function, and discuss emerging data advocating their usage for the development of therapies alleviating neuromuscular dysfunction.

High-Throughput Identification and Analysis of Novel Conotoxins from Three Vermivorous Cone Snails by Transcriptome Sequencing

The venom of each Conus species consists of a diverse array of neurophysiologically active peptides, which are mostly unique to the examined species. In this study, we performed high-throughput transcriptome sequencing to extract and analyze putative conotoxin transcripts from the venom ducts of 3 vermivorous cone snails (C. caracteristicus, C. generalis, and C. quercinus), which are resident in offshore waters of the South China Sea. In total, 118, 61, and 48 putative conotoxins (across 22 superfamilies) were identified from the 3 Conus species, respectively; most of them are novel, and some possess new cysteine patterns. Interestingly, a series of 45 unassigned conotoxins presented with a new framework of C-C-C-C-C-C, and their mature regions were sufficiently distinct from any other known conotoxins, most likely representing a new superfamily. O- and M-superfamily conotoxins were the most abundant in transcript number and transcription level, suggesting their critical roles in the venom functions of these vermivorous cone snails. In addition, we identified numerous functional proteins with potential involvement in the biosynthesis, modification, and delivery process of conotoxins, which may shed light on the fundamental mechanisms for the generation of these important conotoxins within the venom duct of cone snails.

Conodipine-P1-3, the First Phospholipases A2 Characterized from Injected Cone Snail Venom

The phospholipase A₂ (PLA₂s) superfamily are ubiquitous small enzymes that catalyze the hydrolysis of phospholipids at the sn-2 ester bond. PLA₂s in the venom of cone snails (conodipines, Cdpi) are composed of two chains termed as alpha and beta subunits. Conodipines are categorized within the group IX of PLA₂s. Here we describe the purification and biochemical characterization of three conodipines (Cdpi-P1, -P2 and -P3) isolated from the injected venom of Conus purpurascens Using proteomics methods, we determined the full sequences of all three conodipines. Conodipine-P1-3 have conserved consensus catalytic domain residues, including the Asp/His dyad. Additionally, these enzymes are expressed as a mixture of proline hydroxylated isoforms. The activities of the native Conodipine-Ps were evaluated by conventional colorimetric and by MS-based methods, which provide the first detailed cone snail venom conodipine activity monitored by mass spectrometry. Conodipines can have medicinal applications such inhibition of cancer proliferation, bacterial and viral infections among others.

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 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.

Venomix: a simple bioinformatic pipeline for identifying and characterizing toxin gene candidates from transcriptomic data

The advent of next-generation sequencing has resulted in transcriptome-based approaches to investigate functionally significant biological components in a variety of non-model organism. This has resulted in the area of “venomics”: a rapidly growing field using combined transcriptomic and proteomic datasets to characterize toxin diversity in a variety of venomous taxa. Ultimately, the transcriptomic portion of these analyses follows very similar pathways after transcriptome assembly often including candidate toxin identification using BLAST, expression level screening, protein sequence alignment, gene tree reconstruction, and characterization of potential toxin function. Here we describe the Python package Venomix, which streamlines these processes using common bioinformatic tools along with ToxProt, a publicly available annotated database comprised of characterized venom proteins. In this study, we use the Venomix pipeline to characterize candidate venom diversity in four phylogenetically distinct organisms, a cone snail (Conidae; Conus sponsalis), a snake (Viperidae; Echis coloratus), an ant (Formicidae; Tetramorium bicarinatum), and a scorpion (Scorpionidae; Urodacus yaschenkoi). Data on these organisms were sampled from public databases, with each original analysis using different approaches for transcriptome assembly, toxin identification, or gene expression quantification. Venomix recovered numerically more candidate toxin transcripts for three of the four transcriptomes than the original analyses and identified new toxin candidates. In summary, we show that the Venomix package is a useful tool to identify and characterize the diversity of toxin-like transcripts derived from transcriptomic datasets. Venomix is available at: https://bitbucket.org/JasonMacrander/Venomix/.

Mitochondrial genome sequencing of a vermivorous cone snail Conus quercinus supports the correlative analysis between phylogenetic relationships and dietary types of Conus species

Complete mitochondrial genome (mitogenome) sequence of a worm-hunting cone snail, Conus quercinus, was reported in this study. Its mitogenome, the longest one (16,460 bp) among reported Conus specie, is composed of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and one D-loop region. The mitochondrial gene arrangement is highly-conserved and identical to other reported. However, the D-loop region of C. quercinus is the longest (943 bp) with the higher A+T content (71.3%) and a long AT tandem repeat stretch (68 bp). Subsequent phylogenetic analysis demonstrated that three different dietary types (vermivorous, molluscivorous and piscivorous) of cone snails are clustered separately, suggesting that the phylogenetics of cone snails is related to their dietary types. In conclusion, our current work improves our understanding of the mitogenomic structure and evolutionary status of the vermivorous C. quercinus, which support the putative hypothesis that the Conus ancestor was vermivorous.

Evaluating the Performance of De Novo Assembly Methods for Venom-Gland Transcriptomics

Venom-gland transcriptomics is a key tool in the study of the evolution, ecology, function, and pharmacology of animal venoms. In particular, gene-expression variation and coding sequences gained through transcriptomics provide key information for explaining functional venom variation over both ecological and evolutionary timescales. The accuracy and usefulness of inferences made through transcriptomics, however, is limited by the accuracy of the transcriptome assembly, which is a bioinformatic problem with several possible solutions. Several methods have been employed to assemble venom-gland transcriptomes, with the Trinity assembler being the most commonly applied among them. Although previous evidence of variation in performance among assembly software exists, particularly regarding recovery of difficult-to-assemble multigene families such as snake venom metalloproteinases, much work to date still employs a single assembly method. We evaluated the performance of several commonly used de novo assembly methods for the recovery of both nontoxin transcripts and complete, high-quality venom-gene transcripts across eleven snake and four scorpion transcriptomes. We varied k-mer sizes used by some assemblers to evaluate the impact of k-mer length on transcript recovery. We showed that the recovery of nontoxin transcripts and toxin transcripts is best accomplished through different assembly software, with SDT at smaller k-mer lengths and Trinity being best for nontoxin recovery and a combination of SeqMan NGen and a seed-and-extend approach implemented in Extender as the best means of recovering a complete set of toxin transcripts. In particular, Extender was the only means tested capable of assembling multiple isoforms of the diverse snake venom metalloproteinase family, while traditional approaches such as Trinity recovered at most one metalloproteinase transcript. Our work demonstrated that traditional metrics of assembly performance are not predictive of performance in the recovery of complete and high quality toxin genes. Instead, effective venom-gland transcriptomic studies should combine and quality-filter the results of several assemblers with varying algorithmic strategies.

Mass spectrometric identification of bromotryptophan containing conotoxin sequences from the venom of C. amadis

Four 30 residue conotoxin have been identified from the venom of C. amadis. MS/MS analysis of crude venom subjected to global reduction/alkylation yielded fragmentation patterns, which permitted searching and matching with a database of putative mature toxin sequences obtained from transcriptomic analysis. Of the four sequences identified, Am3408(Am6.1b), Am3452(Am6.1c), Am3136(Am6.2a) and Am3214(Am6.2b), three contain bromotryptophan residues, while an additional post translational modification, gamma carboxylation of glutamic acid, is present in Am3408(Am6.1b)/3452(Am6.1c). The conotoxins belong to the O1/O2 gene superfamily and possess cysteine framework VI/VII. While, the cysteine patterns show a similarity to omega conotoxins, the three C. amadis peptides are highly negatively charged and possess a significant content of hydrophobic residues.

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.

Divergence of the Venom Exogene Repertoire in Two Sister Species of Turriconus

The genus Conus comprises approximately 700 species of venomous marine cone snails that are highly efficient predators of worms, snails, and fish. In evolutionary terms, cone snails are relatively young with the earliest fossil records occurring in the Lower Eocene, 55 Ma. The rapid radiation of cone snail species has been accompanied by remarkably high rates of toxin diversification. To shed light on the molecular mechanisms that accompany speciation, we investigated the toxin repertoire of two sister species, Conus andremenezi and Conus praecellens, that were until recently considered a single variable species. A total of 196 and 250 toxin sequences were identified in the venom gland transcriptomes of C. andremenezi and C. praecellens belonging to 25 and 29 putative toxin gene superfamilies, respectively. Comparative analysis with closely (Conus tribblei and Conus lenavati) and more distantly related species (Conus geographus) suggests that speciation is associated with significant diversification of individual toxin genes (exogenes) whereas the expression pattern of toxin gene superfamilies within lineages remains largely conserved. Thus, changes within individual toxin sequences can serve as a sensitive indicator for recent speciation whereas changes in the expression pattern of gene superfamilies are likely to reflect more dramatic differences in a species' interaction with its prey, predators, and competitors.

Venomics of Remipede Crustaceans Reveals Novel Peptide Diversity and Illuminates the Venom's Biological Role

We report the first integrated proteomic and transcriptomic investigation of a crustacean venom. Remipede crustaceans are the venomous sister group of hexapods, and the venom glands of the remipede Xibalbanus tulumensis express a considerably more complex cocktail of proteins and peptides than previously thought. We identified 32 venom protein families, including 13 novel peptide families that we name xibalbins, four of which lack similarities to any known structural class. Our proteomic data confirm the presence in the venom of 19 of the 32 families. The most highly expressed venom components are serine peptidases, chitinase and six of the xibalbins. The xibalbins represent Inhibitory Cystine Knot peptides (ICK), a double ICK peptide, peptides with a putative Cystine-stabilized α-helix/β-sheet motif, a peptide similar to hairpin-like β-sheet forming antimicrobial peptides, two peptides related to different hormone families, and four peptides with unique structural motifs. Remipede venom components represent the full range of evolutionary recruitment frequencies, from families that have been recruited into many animal venoms (serine peptidases, ICKs), to those having a very narrow taxonomic range (double ICKs), to those unique for remipedes. We discuss the most highly expressed venom components to shed light on their possible functional significance in the predatory and defensive use of remipede venom, and to provide testable ideas for any future bioactivity studies.

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.

Screening and Validation of Highly-Efficient Insecticidal Conotoxins from a Transcriptome-Based Dataset of Chinese Tubular Cone Snail

Most previous studies have focused on analgesic and anti-cancer activities for the conotoxins identified from piscivorous and molluscivorous cone snails, but little attention has been devoted to insecticidal activity of conotoxins from the dominant vermivorous species. As a representative vermivorous cone snail, the Chinese tubular cone snail (Conus betulinus) is the dominant Conus species inhabiting the South China Sea. We sequenced related venom transcriptomes from C. betulinus using both the next-generation sequencing and traditional Sanger sequencing technologies, and a comprehensive library of 215 conotoxin transcripts was constructed. In our current study, six conotoxins with potential insecticidal activity were screened out from our conotoxin library by homologous search with a reported positive control (alpha-conotoxin ImI from C. imperialis) as the query. Subsequently, these conotoxins were synthesized by chemical solid-phase and oxidative folding for further insecticidal activity validation, such as MTT assay, insect bioassay and homology modeling. The final results proved insecticidal activities of our achieved six conotoxins from the transcriptome-based dataset. Interestingly, two of them presented a lot of high insecticidal activity, which supports their usefulness for a trial as insecticides in field investigations. In summary, our present work provides a good example for high throughput development of biological insecticides on basis of the accumulated genomic resources.

De novo assembly and annotation of the European abalone Haliotis tuberculata transcriptome

The European abalone Haliotis tuberculata is a delicacy and consequently a commercially valuable gastropod species. Aquaculture production and wild populations are subjected to multiple climate-associated stressors and anthropogenic pressures, including rising sea-surface temperatures, ocean acidification and an emerging pathogenic Vibrio infection. Transcript expression data provides a valuable resource for understanding abalone responses to variation in the biotic and abiotic environment. To generate an extensive transcriptome, we performed next-generation sequencing of RNA on larvae exposed to temperature and pH variation and on haemolymph of adults from two wild populations after experimental infection with Vibrio harveyi. We obtained more than 1.5 billion raw paired-end reads, which were assembled into 328,519 contigs. Filtration and clustering produced a transcriptome of 41,099 transcripts, of which 10,626 (25.85%) were annotated with Blast hits, and 7380 of these were annotated with Gene Ontology (GO) terms in Blast2Go. A differential expression analysis comparing all samples from the two life stages identified 5690 and 10,759 transcripts with significantly higher expression in larvae and adult haemolymph respectively. This is the greatest sequencing effort yet in the Haliotis genus, and provides the first high-throughput transcriptomic resource for H. tuberculata.

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.

From Mollusks to Medicine: A Venomics Approach for the Discovery and Characterization of Therapeutics from Terebridae Peptide Toxins

Animal venoms comprise a diversity of peptide toxins that manipulate molecular targets such as ion channels and receptors, making venom peptides attractive candidates for the development of therapeutics to benefit human health. However, identifying bioactive venom peptides remains a significant challenge. In this review we describe our particular venomics strategy for the discovery, characterization, and optimization of Terebridae venom peptides, teretoxins. Our strategy reflects the scientific path from mollusks to medicine in an integrative sequential approach with the following steps: (1) delimitation of venomous Terebridae lineages through taxonomic and phylogenetic analyses; (2) identification and classification of putative teretoxins through omics methodologies, including genomics, transcriptomics, and proteomics; (3) chemical and recombinant synthesis of promising peptide toxins; (4) structural characterization through experimental and computational methods; (5) determination of teretoxin bioactivity and molecular function through biological assays and computational modeling; (6) optimization of peptide toxin affinity and selectivity to molecular target; and (7) development of strategies for effective delivery of venom peptide therapeutics. While our research focuses on terebrids, the venomics approach outlined here can be applied to the discovery and characterization of peptide toxins from any venomous taxa.

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.

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.

TRUFA: A User-Friendly Web Server for de novo RNA-seq Analysis Using Cluster Computing

Application of next-generation sequencing (NGS) methods for transcriptome analysis (RNA-seq) has become increasingly accessible in recent years and are of great interest to many biological disciplines including, eg, evolutionary biology, ecology, biomedicine, and computational biology. Although virtually any research group can now obtain RNA-seq data, only a few have the bioinformatics knowledge and computation facilities required for transcriptome analysis. Here, we present TRUFA (TRanscriptome User-Friendly Analysis), an open informatics platform offering a web-based interface that generates the outputs commonly used in de novo RNA-seq analysis and comparative transcriptomics. TRUFA provides a comprehensive service that allows performing dynamically raw read cleaning, transcript assembly, annotation, and expression quantification. Due to the computationally intensive nature of such analyses, TRUFA is highly parallelized and benefits from accessing high-performance computing resources. The complete TRUFA pipeline was validated using four previously published transcriptomic data sets. TRUFA’s results for the example datasets showed globally similar results when comparing with the original studies, and performed particularly better when analyzing the green tea dataset. The platform permits analyzing RNA-seq data in a fast, robust, and user-friendly manner. Accounts on TRUFA are provided freely upon request at https://trufa.ifca.es.