Transcriptomic messiness in the venom duct of Conus miles contributes to conotoxin diversity

Marine cone snails have developed sophisticated chemical strategies to capture prey and defend themselves against predators. Among the vast array of bioactive molecules in their venom, peptide components called conotoxins or conopeptides dominate, with many binding with high affinity and selectivity to a broad range of cellular targets, including receptors and transporters of the nervous system. Whereas the conopeptide gene precursor organization has a conserved topology, the peptides in the venom duct are highly processed. Indeed, deep sequencing transcriptomics has uncovered on average fewer than 100 toxin gene precursors per species, whereas advanced proteomics has revealed >10-fold greater diversity at the peptide level. In the present study, second-generation sequencing technologies coupled to highly sensitive mass spectrometry methods were applied to rapidly uncover the conopeptide diversity in the venom of a worm-hunting species, Conus miles. A total of 662 putative conopeptide encoded sequences were retrieved from transcriptomic data, comprising 48 validated conotoxin sequences that clustered into 10 gene superfamilies, including 3 novel superfamilies and a novel cysteine framework (C-C-C-CCC-C-C) identified at both transcript and peptide levels. A surprisingly large number of conopeptide gene sequences were expressed at low levels, including a series of single amino acid variants, as well as sequences containing deletions and frame and stop codon shifts. Some of the toxin variants generate alternative cleavage sites, interrupted or elongated cysteine frameworks, and highly variable isoforms within families that could be identified at the peptide level. Together with the variable peptide processing identified previously, background genetic and phenotypic levels of biological messiness in venoms contribute to the hypervariability of venom peptides and their ability to evolve rapidly.