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.
Keywords: 454 sequencing; Conus imperialis; cone snail venom; conopeptide; conotoxin; iTRAQ; mass spectrometry; proteomics; transcriptomics; venom transcriptome.