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

Jutty Rajan Prashanth, Sebastien Dutertre, Subash Kumar Rai, Richard J Lewis

<figure><img src="image_1.png" /></figure>

## To cite this version:

Jutty Rajan Prashanth, Sebastien Dutertre, Subash Kumar Rai, Richard J Lewis. Venomics Reveals a Non-Compartmentalised Venom Gland in the Early Diverged Vermivorous Conus distans. Toxins, 2022, 14 (3), pp.226. 10.3390/toxins14030226. [https://doi.org/10.3390/toxins14030226](https://doi.org/10.3390/toxins14030226)

HAL Id: hal-04185884

https://hal.science/hal-04185884v1

Submitted on 23 Aug 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

Article

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

Jutty Rajan Prashanth<fnref n="1" />, Sebastien Dutertre<fnref n="2" /> ID, Subash Kumar Rai<fnref n="1" /> ID and Richard J. Lewis<fnref n="1" />,* ID

**Abstract:** 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.

**Keywords:** proteomics; transcriptomics; evolution; defensive venom; conotoxins

**Key Contribution:** This study reveals that the early diverged vermivorous *Conus distans* has an unspecialised venom duct. Our 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. *C. planorbis* represents an intermediate evolutionary stage in the transition from non-compartmentalised venom glands to the highly compartmentalised venom glands found in more recently evolved mollusc and fish-hunting species.

## 1. Introduction

Venoms have evolved multiple times across different lineages to aid predation, defence and competitor deterrence [1]. They are injected into the target species through specialised apparatuses that have evolved under ecological pressures together with the venom peptides they deliver [2,3]. The biological activity of venom peptides and the mechanisms to express and inject venom are influenced by ecological interactions with prey and predators [4–6], with some species of molluscs (cone snails) and insects (assassin bugs) able to deploy separately evolved predatory and defensive venoms [7,8].

Cone snails are circumglobal venomous marine gastropods found in tropical waters that use venoms for predation and defence. More than 800 species of cone snails have been identified mostly in the Indo-Pacific region, including ~700 species of the genus *Conus* [9]. Ancestral and most extant cone snail species are worm-hunters, with mollusc- and fish-hunting species having evolved more recently [10]. Some cone snail species, including the fish-hunting *Conus geographus* and the mollusc hunting *Conus marmoreus*, inject different venoms for defence and predation [7]. These distinct defensive and predatory venoms are injected using the same radula tooth envenomation apparatus but are differentially

<figure><img src="image_4.png" /></figure>

check for
updates

**Citation:** Prashanth, J.R.; Dutertre, S.; Rai, S.K.; Lewis, R.J. Venomics Reveals a Non-Compartmentalised Venom Gland in the Early Diverged Vermivorous *Conus distans*. *Toxins* **2022**, *14*, 226. https://doi.org/10.3390/toxins14030226

Received: 17 February 2022

Accepted: 16 March 2022

Published: 19 March 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_5.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Institute for Molecular Bioscience, University of Queensland, St. Lucia, QLD 4072, Australia; prashanth.rajan185@gmail.com (J.R.P.); subash.rai@uq.edu.au (S.K.R.)</footnote>

<footnote n="2">² IBMM, Université Montpellier, CNRS, ENSCM, 34095 Montpellier, France; sebastien.dutertre@umontpellier.fr</footnote>

<footnote>* Correspondence: r.lewis@imb.uq.edu.au</footnote>
expressed from different sections of the venom duct [11]. Based on the venom duct analysis of *C. geographus*, defence-evoked venoms are typically expressed in the proximal section (near the venom bulb) and those used for predation in the distal section (near the proboscis) [7]. However, in the worm-hunting *Conus planorbis*, the defensive venom peptides were not restricted to the proximal section, despite variable venom production along the duct [12]. In contrast, more recently diverged worm-hunting species, such as *Conus vexillum* [13], produced defensive venoms that dominated only the proximal duct. We hypothesised that the defensive use of venom first arose in ancestral worm-hunting snails and was later repurposed in a compartmentalised venom duct to facilitate the dietary shift to molluscivory and piscivory [7].

To better understand when venom gland compartmentalisation might have evolved in cone snails, we studied the venom of the earliest diverging extant Conidae, *Conus distans*, from the *Fraterconus* subgenus [9,14]. The venom of this species is little studied, with previous assay-guided efforts only identifying proteins between 24 kDa and 25.5 kDa that blocked calcium influx associated with neurotransmitter release in the hippocampus of rats and bovine chromaffin cells [15,16] and a novel excitatory peptide, DiXIXA, possessing a rare triple cysteine motif [17].

Consistent with its placement within a basal lineage, we demonstrate here 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. We compared *C. distans* to *C. planorbis*, the earliest diverging species from which a defence-evoked venom has been obtained [12], and *C. geographus*, a fish-hunting cone snail from the *Gastridium* subgenus that injects highly distinct defensive and predatory venoms [7]. These data support the hypothesis that venom gland compartmentalisation likely arose in an ancestral worm-hunting species and facilitated the dietary shift from vermivory to molluscivory and piscivory in cone snails.

## 2. Results

### 2.1. *The Venom Gland of C. distans Is Not Compartmentalised*

To examine whether venom expression in *C. distans* varied along the venom gland, we performed LC-ESI-MS experiments on venom extracted from the distal (closest to the proboscis), distal central, proximal central and proximal (closest to the venom bulb) venom duct of *C. distans* (Specimen-1). Comparing base peak chromatograms (BPCs), we observed that the elution time of all major masses were similar across the four sections of the gland (Figure 1). Manual reconstruction of the masses within each co-eluting peak confirmed minimal qualitative or quantitative differences between venom peptides secreted from each section, except for several late eluting minor peaks (Figure 1A(i)–D(i)). We also reconstructed the masses using the LC-MS reconstruction and observed that the mass distributions within the four sections were similar, with the venom profile of *C. distans* dominated by 1–4 kDa peptides (Figure 1A(ii)–D(ii)). Each section also had a significant number of masses <1 kDa, which were not characterised further but may include small molecules similar to those discovered in another early diverging species, *Conus imperalis* [18]. We could not detect any masses corresponding to DiXIXA, the only peptide isolated from this species previously [17], but were able to identify two peptide masses (Dis41 and Dis63; Figure 1) that corresponded to two novel conotoxins we identified in the venom gland transcriptome of another *C. distans* specimen.

To confirm whether *C. distans* has a non-compartmentalised venom duct, we analysed a second specimen obtained from a different geographical location (Specimen 2). Again, the BPC revealed that the expression of major masses along the venom gland was uniform (Supplementary Figure S1). We manually reconstructed the masses in each major peak for both specimens and compared the expression levels across the four sections of their respective glands. Despite substantial differences in the specific masses detected between the two specimens, the expression levels of major common masses were similar across the
length of the gland in both species, demonstrating that the venom gland of *C. distans* is indeed non-compartmentalised (Figure 2).

<figure><img src="image_6.png" /><figcaption>**Figure 1.** Proteomics of venom gland sections of *C. distans*. (A(i)–D(i)) show LC-ESI-MS profiles of the four venom duct sections, with major masses labelled. (A(ii)–D(ii)) show mass distributions for the corresponding sections.</figcaption></figure>
Relative Expression Levels of Peak Masses

<figure><img src="image_7.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_8.png" /></figure>

**Figure 2.** Relative expression levels of dominant masses in *C. distans* venom gland sections. Comparison of relative expression levels (% total area under the curve) of masses across the four sections of the venom gland. (A) LC-ESI-MS chromatograms of (A) Specimen 1 and (B) Specimen 2 were run for 65 min and 100 min, respectively.

## 2.2. *C. distans* Venom Gland Transcriptome

We sequenced the venom gland transcriptome of the third *C. distans* (Specimen 3) using 454-pyrosequencing to identify the expressed venom peptide sequences. Specimen 3 produced 209,346 total reads, including 183,763 reads with a sequence quality score cut-off of >30 and an average length of ~332 nucleotides. Trimmed reads (see Materials and Methods) contained 135 conotoxin transcripts belonging to 25 different superfamilies,
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 3.** Venom gland transcriptome of *C. distans*. (A) Expression levels of various superfamilies in the terms of number of transcripts and sequence reads. (B) Highly expressed transcripts (expression level >1% of total reads) in *C. distans*. The superfamily and predicted cysteine framework are indicated. (C) The distribution of cysteine frameworks in the transcriptome dataset.</figcaption></figure>
<table><caption>**Table 1.** Conotoxins in *C. distans*.</caption><thead><tr><th>Name (Superfamily) and Predicted Mass (Da)</th><th>Predicted Mature Peptide</th><th>Framework</th><th>Related Conotoxins</th><th>Likely Activity</th></tr></thead><tbody><tr><td>Dis2 (A)—3582.4</td><td>RIAEPNTEEWNECKDPSCRNNHMDRCAE</td><td>CC-C-C</td><td>Di1.1, Di1.2 (C. distans)</td><td>nAChR antagonists</td></tr><tr><td>Dis4 (B1)—1888.9</td><td>TITAQEAETARERLSTL</td><td>Linear</td><td>Conantokin-E1 (C. eburneus)</td><td>NMDA receptor antagonist</td></tr><tr><td>Dis9 (B1)—1234.7</td><td>LVGEVEIIVHK</td><td>Linear</td><td>Conantokin-V (C. vitulinus)</td><td>NMDA receptor antagonist</td></tr><tr><td>Dis14 (B1)—1947.9</td><td>QTEEVEESQEKLEEL</td><td>Linear</td><td>Conantokin-G (C. geographus)</td><td>NMDA receptor antagonist</td></tr><tr><td>Dis17 (Con-ikot-ikot)—12726.0</td><td>GDPTAECLTLGLCYDTCSHDNTSPDCWGHCKS ESTT GCSLDFSLHYCSEFQDCYGPCVTDKDE SRCFKACRKEAMIDCLDHGSVECCPGFVNCYQN CRRSREDCITHCLKENC</td><td>N/A</td><td>Con-ikot-ikot (C. striatus)</td><td>AMPA receptor antagonist</td></tr><tr><td>Dis20 (SFMi2)—2928.9</td><td>DCGRDCVGCDPANCCCGQTCVNGNKCE</td><td>XXVII</td><td>Mi045, Mi046 (C. miles)</td><td>N/A</td></tr><tr><td>Dis22 (M—L-LLTVA)—6305.5</td><td>DQTASGFLRDDESVFPCNSDRCACLPKEG STTSYQCQSLEASTDDCVNNECITEDEW</td><td>IX</td><td>Cl9.5 (C. californicus), Turripptide IX-01</td><td>N/A</td></tr><tr><td>Dis30 (M—L-LLTVA) (DiXIXA)—4953.8</td><td>GCDPTDGCQTTVCETDTGPC CCKPNFTCQISNSGTKSCSCSGQPSDCPV</td><td>XIX</td><td>N/A</td><td>Hyperexcitability/lethargy (dose-dependent)</td></tr><tr><td>Dis38 (H)—3984.5</td><td>TDVDCGGVSCTFGCCETVNGEKKCKELDCDVTSDTENS</td><td>VI/VII</td><td>Teretoxin Tsu6.5</td><td>N/A</td></tr><tr><td>Dis43 (I1)—3378.3</td><td>CSYSSCKTESCCTGYLCNSVKSCVDPNSGGRF</td><td>VI/VII</td><td>Im11.14 (C. imperialis)</td><td>N/A</td></tr><tr><td>Dis45 (L)—2026.8</td><td>RCPIACKTCPDENTCIPAP</td><td>XIV</td><td>Cl14.2 (C. californicus)</td><td>N/A</td></tr><tr><td>Dis46 (M)—3711.5</td><td>TCDPYYCNDGKVCCPEYPTCGDSTGKLICVRVTD</td><td>VI/VII</td><td>Im6.7 (C. imperialis)</td><td>N/A</td></tr><tr><td>Dis49 (M)—1801.8</td><td>DVKCIGSCDSTVWHRV</td><td>N/A</td><td>M-conotoxin 6 (C. marmoreus)</td><td>N/A</td></tr><tr><td>Dis53 (NSDis1)—5430.2</td><td>SCEGNDSYCRKPDWISDKPCCDPLVC VCTGPMSGGRRRTTCETGKCGPRPSS</td><td>XV</td><td>N/A</td><td>N/A</td></tr><tr><td>Dis58 (NSDis3)—3668.3</td><td>NACELDSSTGDDCTGTQICCTPSGSMSGECREADEC</td><td>VI/VII</td><td>Lp7.1 (C. leopardus)</td><td>N/A</td></tr><tr><td>Dis60 (NSG3)—6759.6</td><td>DESQCPECRCHDLKNAICDISEAC SDEASCPTSPPECNNGNCLCKNFHGGRCVSHSQCNDRHC</td><td>XXVIII</td><td>G125, G126 (C. geographus)</td><td>N/A</td></tr><tr><td>Dis62 (O1)—3509.3</td><td>CDPPGYPCELRENDCCDACKIVSQNPNVCSDE</td><td>VI/VII</td><td>MiK42 (C. miles)</td><td>ω-conotoxin-like</td></tr><tr><td>Dis63 (O1)—3244.3</td><td>CLSIGYACGVAISEKCCHVCDNPTGAGTCVYN</td><td>VI/VII</td><td>Im6.1 (C. imperialis)</td><td>ω-conotoxin-like</td></tr><tr><td>Dis68 (O1)—3205.2</td><td>SCAESGLSCDTRPCCDDKTCVRNGRQSMCS</td><td>VI/VII</td><td>Mr6.11 (C. marmoreus)</td><td>N/A</td></tr></tbody></table>
<table><caption>Table 1. *Cont.*</caption><thead><tr><th>Name (Superfamily) and Predicted Mass (Da)</th><th>Predicted Mature Peptide</th><th>Framework</th><th>Related Conotoxins</th><th>Likely Activity</th></tr></thead><tbody><tr><td>Dis70 (O1)—3078.1</td><td>CTESGLTCWP TNHDCCSGTCNGTMTTGTCT</td><td>VI/VII</td><td>Mi024 (C. miles)</td><td>ω-conotoxin-like</td></tr><tr><td>Dis71 (O1)—3087.0</td><td>CEDEGSPCQFDSECCSGACTPEGVFDFCE</td><td>VI/VII</td><td>ω-like peptide (C. capitaneus)</td><td>ω-conotoxin-like</td></tr><tr><td>Dis72 (O1)—2969.1</td><td>SCANYHESCASDPCCEGLECIGAQGGGVCI</td><td>VI/VII</td><td>ArMKLT2-0322 (C. arenatus)</td><td>N/A</td></tr><tr><td>Dis76 (O1)—2789.9</td><td>CKGTDAPCDDHDECCEHVCDGVCVED</td><td>VI/VII</td><td>MgJr94P (C. magus)</td><td>ω-conotoxin-like</td></tr><tr><td>Dis80 (O1)—2503.8</td><td>CEDPGEPCGSDHSCCSGSCNHNVCA</td><td>VI/VII</td><td>ArMKLT2-0322 (C. arenatus)</td><td>ω-conotoxin-like</td></tr><tr><td>Dis81 (O2)—8000.6</td><td>DEDCVTEEGDHVEEGKASLLTTATRVPAMMENCTVR LCTVN NGVTEVLTNHENSVDDSLITWNLWTPCHCFI</td><td>XIV</td><td>N/A</td><td>N/A</td></tr><tr><td>Dis85 (O2)—4165.4</td><td>QCAPDDFQCDVDEDDCCNDLECKCFTSTDCTSGYKCRN</td><td>VI/VII</td><td>Bt15a (C. betulinus)</td><td>N/A</td></tr><tr><td>Dis91 (O2)— 3759.3</td><td>CKGASALCEEDGECCSGDCKCMHASGCTNDINLRCAA</td><td>VI/VII</td><td>Di6.5 (C. distans)</td><td>N/A</td></tr><tr><td>Dis97 (O2)— 3768.4</td><td>TCKGRLQSCDHDSECCSPYTCYCGMKQG CNLKCI</td><td>XV</td><td>Vt15a (C. vitulinus)</td><td>N/A</td></tr><tr><td>Dis98 (O2)—3264.1</td><td>CDTWRDPCTYEHECCWQYHCGFRTCE</td><td>VI/VII</td><td>Di6.10 (C. distans)</td><td>N/A</td></tr><tr><td>Dis102 (O2—contryphan)—1316.5</td><td>EFDCPWHPWCC</td><td>N/A</td><td>Contryphan-R-like (C. virgo)</td><td>Contryphan-like activity</td></tr><tr><td>Dis104 (O2)—4409.7</td><td>HGGCLNEEGDLVAEDGETIEVE CNRCRCEDGDLACTKMACE</td><td>IX</td><td>Di6.5 (C. distans)</td><td>N/A</td></tr><tr><td>Dis116 (O3)—3535.5</td><td>NVDQECIDACQLEDKNCCGRTDGEPRCAKICL</td><td>VI/VII</td><td>Di6.6 (C. distans)</td><td>N/A</td></tr><tr><td>Dis119 (O3)—3871.5</td><td>LVEGACTSPSN CPTGQECCPNKLDEPEGSCANDCPFY</td><td>VI/VII</td><td>Di6.12 (C. distans)</td><td>N/A</td></tr><tr><td>Dis120 (P)—2780.9</td><td>STCPTSCATHMNCWPECTYCTTSGCT</td><td>IX</td><td>Fla9.1 (C. flavidus), TxIXA (C. textile), GmIXA (C. gloriamaris)</td><td>spasmodic</td></tr><tr><td>Dis121 (S)—3803.1</td><td>SSCSGTCYGSANCDGTCYCREDNCWCTGDSSCACQCA</td><td>VIII</td><td>Di8.1 (C. distans), GVIIIA, RVIIIA</td><td>5-HT3 receptor and nAChR antagonists</td></tr><tr><td>Dis124 (T)—8350.9</td><td>TPSEQNLPGELTPADLEGAETTPEESWYSKIKSGVKHASC KLVGYACDDSETEESLLSKIKGGVEHAACKYVNIACED</td><td>XIV</td><td>Tr5.4 (C. terebra)</td><td>N/A</td></tr><tr><td>Dis134 (T)—6399.6</td><td>APSEPNLQRGLKLGGLKAEPNLQRGLK LGGVKDKLLKVGGNILKGAVQGAVDSLTKEDRKQ</td><td>N/A</td><td>Vr5.4 (C. varius)</td><td>N/A</td></tr></tbody></table>

N/A not ascertained.
<table><caption>**Table 2.** MS/MS identification of transcriptomic sequences.</caption><thead><tr><th></th><th>Sequence</th></tr></thead><tbody><tr><td>Dis1</td><td>MGMRMMFIVFLLVVLATTVVLRSDFARNRKNRRRIAEPNTEEEWNECCKDPSCRNNHLDRCPE</td></tr><tr><td>Dis10</td><td>MELYTYLLVPLVAFHLIQGTGTRSHGGPLTEGRSADVTALKPEPVLLQKSDARSTDNGKDKLTRMKRTLKKGGNMARRQTEEVEESNETLAEAGKR</td></tr><tr><td>Dis102</td><td>MKKTILVLVAAVLLSTQVMVQGDGDQPADRNAVRDANPGGTSGKLMRVLQGREFDCPWHPWCG</td></tr><tr><td>Dis103</td><td>MNKLTMLILVATVLLSIQVMVRGDEDCVNEEGDLVAEDGETVKVECNBTrCRCDDGDLACTKMACE</td></tr><tr><td>Dis11</td><td>MQLTYLLVPLVAFHLIQGTGTRGHGGALTEGRSADVTALKPEPVLLQKSDARSADDNGKDKLTRMRRTLKNKGNMARRQTEEVEESQEKLEELGKR</td></tr><tr><td>Dis118</td><td>MSGLGIMVLTLLLLVPMATSQQDGEKQAMQRDAINAAPGTSITRRETDQECIDTCEQEDKKCCGRTNGEPVCAKICFG</td></tr><tr><td>Dis121</td><td>MMSKMGAMFVLLLLCPLASNQQEGDIKARRTFWKRDLYGDLAGRSSCSGTCYGSANCDGTCYCREDCNCWCTGDSSCACQCA</td></tr><tr><td>Dis125</td><td>MLCLPVFIILLLLASPAVTTPSEQNLPGELTPADLEGAETTPEESWYSKIKGGVKHASCKLVGYACDDSETEESLLSKIKGGVEHAACKYVNIGCED</td></tr><tr><td>Dis126</td><td>MLCLPVFIILLLLASPAVTTLSEQNLPGELTPADLEGAETTPEESWYSKIKSGVKHASCKLVGYACDDSETEESLLSKIKGGVEHAACKYVNIACED</td></tr><tr><td>Dis127</td><td>MLCLPVFIILLLLASPAVTTPSEQNLPGELTPADLEGAETTPEESWYSKIKSGVKHASCKLVGYACDDSETEESLLSKIKGGVEYAACKYVNIACED</td></tr><tr><td>Dis128</td><td>MLCLPVFIILLLLASPAVTTPSEQNLPGELTPADLEGAETTOEESWYSKIKGGVKHASCKLVGYACDDSETEESLLSKIKGSVEHAACKYVNIGCED</td></tr><tr><td>Dis132</td><td>MLCLPVFIILLLLASPAVTTPSEQNLPGELTPADLEGAETTPEESWYSKIKGGVKHASCKLVGYACDDSETEESLLSKIKGVSNMLRANTLI</td></tr><tr><td>Dis133</td><td>MLCLPVFIILLLLAAPAVTAPSEPNLQRGLKGGLKAEPNLQRGLKGGLVKEGLLKVGASAIKGAVNGALNSITKEDRKK</td></tr><tr><td>Dis135</td><td>MLCLPVFIILLLLAAPAVTAPSEPNLQRGLKGGLKAEPNLQRGLKGGLVKDKLLKVGGNIFKGAVQGAVDSLTKEDRKQ</td></tr><tr><td>Dis16</td><td>MLRLIITAVLASACLALPHRRDAAPADMGALKPFEQQMOPMGMPGSMAGMQGMPQQAMPGGMLGNOLMPFGPGMGMGAGYRRAADHNQEKRLDPLT</td></tr><tr><td>Dis18</td><td>MNMWMTPSVLVVVVFTATVVVCSTEDERLTRQRRGDPTAECCLTLGCYDTCSHDNTSPDCWGHCKSESTTGCSLDFSLHYCSEFQDCYGPCVTDKDESRCKACRKEAMIDCLDHGSVECCPGFVNCYQNCRRSREDCITYCLKENC</td></tr><tr><td>Dis2</td><td>MGMRMMFIVFLLVVLATTVVLRSDFARNRKNRRRIAEPNTEEEWNECCKDPSCRNNHMDRCAE</td></tr><tr><td>Dis20</td><td>MNFYLLLTVTLLLASFTGGDARRIQGMDIYRHFVRRDCGRDCVGCDDNOANCCCGGQTCVNGNKCE</td></tr><tr><td>Dis21</td><td>MNFYLLLTVTLLLASFTGGDARRIQGMDIYRHFVRRDCGRDCVGCDDNPANCCCGGQTCVNGNKCE</td></tr><tr><td>Dis26</td><td>MGFRQLVTVGLLLTFFMSTDASHADQTESGFLRDETEVFPČNSDRCACLPKEGSTTSYQCQSLEASTDDCVNNECITEDEWSGRR</td></tr><tr><td>Dis27</td><td>MGFRQLVTVGLLLTFFMSTDASHADQTESGFLRDETEVFPČNSDRCACLPKEGSTTSYQCQSLETSTDGCVNNECVTEDEW</td></tr><tr><td>Dis28</td><td>MRFLLRLTVALFLTWFTETDAAAIGKREVHQVILGEPLTNYATVPPDAFQQKLPEIILGQPLMEYQSTESPEVLS</td></tr><tr><td>Dis29</td><td>MRFLLRLTVALFLTWFTETDAAAIGKREVHQVILGEPLTNYVPPDAFQQKLPEIILGQPLMEYQSTESPEVLS</td></tr><tr><td>Dis3</td><td>MQLTYLLVPLVAFHLIQGTGTLGHGGALTEGRSADATAPKPEPVLLQKSDARSADNSKDKLTQMKRTLKKQGHIARTITAEFAERNRERMSTLGKR</td></tr></tbody></table>
<table><caption>**Table 2. Cont.**</caption><thead><tr><th></th><th>Sequence</th></tr></thead><tbody><tr><td>Dis30</td><td>MSTLGILLIALLPLANPAETGDGQAMPRTRNLRSLSGFRTLRRLEKRGCDOTDGCQTTVCETDTGOCCCCKONFTCQISNSGTKSCSCGQOSDCOV</td></tr><tr><td>Dis43</td><td>MKLSVALLLIVLLLPVVAGEKESGDHVLKKRCSYSSCKTESCCTGYLCNSVKSCVDPNSGGRFGK</td></tr><tr><td>Dis44</td><td>MKLSVALLLIVLLLPVVAGEKESGDHVLKKRCSYSSCKTESCCTGYLCNSVKDCVDPNSGGRFGK</td></tr><tr><td>Dis52</td><td>MMTKLGAVTLLSLVIIPQVLLQQHQDGVDVKSMQRNKGRTAAGSVLSHSLRSTNNEYDAKHERSCEGNNSYCRKPDWVGDKPCCSPLVCVCTGTMSGGRRTTCKRAKCGOHOSSK</td></tr><tr><td>Dis54</td><td>MMTKLGAVTLLSLVIIPQVLLQQHQDSIADVKSMERNKGRTAAGSVLSHSLRSTNNEYDAKHKRSCEGNDSYCRKPDWISDKOCCDPLVCVCTGOMSGGRRTTCETGKCGPRPSSK</td></tr><tr><td>Dis55</td><td>MMTKLGAVTLLSLVIIPQVLLQQHQDSIADVKAMERNKGRTAAGSVLSHSLRSTNNEYDTKHKRSCEGNDSYCRKPDWISDKPCCDPLVCVCTGOMSGGRRTTCETGKCGOROSSK</td></tr><tr><td>Dis56</td><td>MMTKLGAVTLLSLVIIPQVLLQQHQDSIADVKAMERNKGRTAAGSVLSHSLRSTNNEYDAKHKRSCEGNDSYCRKPDWISDKPCCDPLICVCTGPMSGGRRTTCETGKCGPRPSSK</td></tr><tr><td>Dis57</td><td>MIQALASMAWTSMLCSADQVSTSPSVPTFVMVLMATVLLTGIMETEARTLFQMIARRSSDYPCAAGTFADCRGQPDGATCCDTGYCQGNVCHY</td></tr><tr><td>Dis59</td><td>MQLSVILFVLLLTMLPFNGSVLNAINGRKTFERNDRSTDSSQMEKRCPTACKSCSOOGTCQPVR</td></tr><tr><td>Dis60</td><td>MKMYLCLAVVLLLASTIVDSALLDKTETLRNWRRKGRDESQCPECRCHDLKNAICDISEACSDEASCPTSPECNNGNCLCKNFHGGRCVSHSQCNDRHC</td></tr><tr><td>Dis61</td><td>MKMYLCLAVVLLLASTIVDSALLDKTETLRNWRRKGRDESQCPECRCHELKNAICDISEACNDEASCPTSPGCNNGNCLCKNFHGGRCVSHSECNDRHC</td></tr><tr><td>Dis64</td><td>MKLTYALIVAVLFLTACQVITTDDSRDKQDLLAMLFSKKRNSRDSKWLTKRCLSIGYACGVAISEKCCHVCDNOTGAGTCVYN</td></tr><tr><td>Dis65</td><td>MKLTYALIVAVLFLTACQVITTDDSRDKQDLLAMLFNKKRNSRDSKWLTKRclsIGYACGVAISEKCCHVCDNPTGAGTCVYN</td></tr><tr><td>Dis66</td><td>MKLTYALIVAVLFLTACQLITTDSDRDKQDLLAMLFNKKRNSRDSKWLAKRclsIGYACGVAISEKCCHVCDNPTGAGTCVYN</td></tr><tr><td>Dis67</td><td>MKLTYALIVVVVLFLTACQQLTADYSRDQKEYPTMFRDQMRNAKGPKWIRSCAESGKSCDTKVCCDDMγCIGTPGGSMCNG</td></tr><tr><td>Dis72</td><td>MKLTCVLVVAVLFLTACQFNTADDSRNKQEYRAARLRVGMQKSNGFRSCANYHESCASDPCCEGLYCIGAQGGGVCI</td></tr><tr><td>Dis73</td><td>MKLTCVLVVAVLFLTACQFNTADDSRNKQEYRAARLRVGMQSKGFRSCANYHESCASDPCCEGLECIGAQGGGVCI</td></tr><tr><td>Dis87</td><td>MKELMILILVATALLSIQVMVRGDGEKPLMGGIKRNAAAGLSALIRGKRCKGTSAICEEDGECCSDDCKCMIASGCSNHINRRCAA</td></tr><tr><td>Dis88</td><td>MKELMILILVATALLSIQVMVRGDGEKPLMGGIKRNAAAGLSALIRGKRCKGSAICEEDGECCSDDCKCMIASGCSNHINRRCAA</td></tr><tr><td>Dis89</td><td>MKELMILILVATALLSIQVMVRGDGEKPLMGGIKRNAAAGLSALIRGKRCKGSAICEEDGECCSDDCKCMIASGCSNHINRRCAA</td></tr><tr><td>Dis90</td><td>MKELMILILVATALLSIQVMVRGDGEKPLMGGVKRNAAAGLSALIRGKRCKGTSACEEDGECCSDDCKCMIASGCSNHINRRCAA</td></tr><tr><td>Dis93</td><td>MKELMILILVATTLLSIRVMVRGDGEKPLMGGIKRNAAAGLSALIRGKRCKGTSALCEEDDECCSGDCKCMHASGCTNDINLRCAA</td></tr><tr><td>Dis95</td><td>MKELMILILVATALLSIQVMVRGDGEKPLMGRIKRNAAAGLSALIRGKRCKGTSALCEEDDECCSGDCKCMIASGCTNDINLRCAA</td></tr><tr><td>Dis96</td><td>MKELMILILVATALLSIQVMVRGDGEKPLMGRIKRNAAAGLSALIRGKRCKGASALCEEDGECCSGDCKCMIASGCTNDINLRCAA</td></tr></tbody></table>

MS validated sequences (>99% confidence) are bolded. Identified PTMs are O = Hydroxyproline, Btr = Bromotryptophan, γ = γ-Carboxyglutamic acid. Underlined residues indicate oxidation.
## 2.4. *Comparison of C. distans with C. planorbis*

C. planorbis has been placed into the Strategoconus subgenus of Conus based on the molecular phylogenetic classification [9]. C. planorbis remains the earliest diverging species from which a defensive venom has been collected [12]. However, Jin et al. reported that unlike the more recently evolved species [7,13], the expression of defensive venom peptides in C. planorbis was not restricted to the proximal sections [12]. Therefore, we hypothesised that the venom gland of C. planorbis might represent a transitionary state between non-compartmentalised venom glands in C. distans and the compartmentalised glands found in piscivorous species such as C. geographus. To understand this comparison better, we performed additional MS experiments on three specimens of C. planorbis and one specimen of C. geographus, obtained from the northern Great Barrier Reef.

To improve spatial resolution, the venom gland was divided into six sections and MS analysis performed as described for *C. distans*. Principal Component Analyses (PCA) [25] were then used to compare the variance between each section of the venom gland. The venom gland of *C. planorbis* shows differentiation along the gland, as previously reported [12], though the levels of variance differ between individuals (Figure 4). Specimen 3 was least differentiated, Specimen 2 was most differentiated, while Specimen 1 showed intermediate levels of differentiation. Specimen 1 had the highest number of masses and was used as the reference specimen for *C. planorbis*.

<figure><img src="image_12.png" /><figcaption>**Figure 4.** Principal Component Analysis (PCA) analysis of *Conus planorbis*. (A(i)–C(ii)) scores plots indicate the diversity of venom gland segments within the *C. planorbis*, and (A(ii)–C(ii)) loading plots show families of correlated variables. Specimen 1 includes the highest number of mass peaks (430) among three specimens in PCA. P: Proximal, PC: Proximal Central, PMC: Proximal Middle Central, DMC: Distal Middle Central, DC: Distal Central, and D: Distal.</figcaption></figure>
Figure 5 compares the PCA results of *C. distans* Specimen 2 with the *C. geographus* specimen and *C. planorbis* Specimen 2. The ion clustering in the different species (Figure 5 B) suggests that the *C. planorbis* venom gland represents an intermediate (transitionary) phase between the non-compartmentalised venom gland found in *C. distans* and the highly compartmentalised venom gland of *C. geographus*.

<figure><img src="image_13.png" /></figure>

**Figure 5.** PCA analysis of venom gland segments of *C. distans*, *C. planorbis* and *C. geographus*. (A(i)–C(i)) PC1 and PC2 scores plots show the separation of venom gland segments in indicated species and (A(ii)–C(ii)) PC1 and PC2 loadings plots show the families of correlated variables. Percentages of variance used in PCA of each species are indicated in the scores plot. P: proximal, PC: proximal central, PMC: proximal middle central, DMC: distal middle central, DC: distal central, and D: distal.

**3. Discussion**

It has been established that mollusc- and fish-hunting cone snails can inject distinct and separate defensive and predatory venoms [7]. Defensive use of venom has been demonstrated in mollusc-, fish- and worm-hunting species [7,13]. From these observations, it was hypothesised that ancestral defensive venom peptides used by worm-hunting species to deter molluscs and fish predators were later repurposed and compartmentalised to facilitate the transition to mollusc- and fish-hunting [7]. However, it remains unclear when compartmentalisation of the venom gland occurred within the genus *Conus* to enable the secretion and injection of distinct venoms for different ecological roles.
Since *C. distans* belongs to the stem group of the genus *Conus*, we studied the expression of venom peptides across the gland in the species and compared it to another early diverging species, *C. planorbis*, as well as a more recently derived fish-hunter, *C. geographus*. Our proteomic data showed that overall expression of venom remained similar throughout the gland with only minor quantitative variations detected (Figure 2). This contrasts with later diverged species previously studied, supporting the hypothesis that ancestral cone snails had a non-compartmentalised venom gland. Thus, these data suggest *C. distans* diverged prior to the appearance of predatory and defensive venom compartmentalisation and likely uses the same venom for both predation and defence. Given that specialised mollusc- and fish-hunting are derived traits that have originated within the genus *Conus* [14], venom gland compartmentalisation and stimulus-dependent venom deployment could have arisen uniquely within this genus. The *C. distans* venom provides an opportunity to study quasi-ancestral conotoxin structure and pharmacology before predatory and defensive venoms diverged.

We surveyed the venom composition of *C. distans* by sequencing its venom gland transcriptome on the 454-pyrosequencing platform, with RNA-seq data as described earlier [13,26]. The venom gland transcriptome was characteristically diverse in terms of superfamilies, cysteine frameworks and probable pharmacology. A total of 135 transcripts belonging to 25 superfamilies were recovered, with early evolved superfamilies such as O₁, O₂ and T dominating (Figure 3). In addition, several divergent superfamilies originally reported in *C. californicus* were also expressed at high levels as well as four novel superfamilies. Structurally, the transcriptome was dominated by peptides containing frameworks VI/VII, IX or XV, which are characterised by their ICK motif (Figure 3C). Interestingly, a framework VI/VII peptide from the H-superfamily was highly similar to a teretoxin, while some framework IX peptides were similar to *C. californicus* and turrid toxins (Table 1). These toxins likely represent toxins inherited from an ancestral toxoglossate mollusc, given their prevalence across separate genera. Several conotoxins and conopeptides that were similar to others found in mollusc- and fish-hunting cone snails were also identified. Conantokins (NMDA receptor antagonists) [27], con-ikot-ikots (AMPA receptor antagonists) [28] and a contryphan [29] are all expressed by *C. distans*. Interestingly, these toxin families were found to be prominent in the predatory venom of the fish-hunting *C. geographus*, suggesting potential repurposing for predation [7]. Two different peptide classes containing unique triple cysteine motifs were also discovered. One of them, DiXIXA, is the only peptide to have been isolated from this species so far, and elicits excitatory activity when injected in mice [17]. The other peptides belonged to SF-Mi2, a recently discovered superfamily reported originally in *C. miles* [21]. Little is known about the structure and the activity of this group of toxins. Thus, the venom gland transcriptome reveals that the venom of *C. distans* contains a number of novel conotoxins whose pharmacology is yet to be characterised.

We also identified two putative α-conotoxins with a novel 4/8 cysteine arrangement. Though the pharmacological and evolutionary significance of this subtype of α-conotoxins has not been defined, α-conotoxins are widely used in defence across the Conidae [7,12,13]. We also found putative ω-conotoxins among the transcriptomic sequences exhibiting a sequence similarity to ω-conotoxins from fish-hunting species. However, using BLAST we could not find any obvious δ-, μ- or μO-like sequences [30–32]. A third type of μ-conotoxins, framework V, containing peptides belonging to the T-superfamily, has also been identified [33] but could not be detected either. κ-Conotoxin K⁺ channel blockers are presumed to have convergently evolved from a range of different structural folds [34]. While *C. purpurascens* employs PnVIIA, a framework VI/VII peptide from the O1-superfamily [35], a kunitz-containing protein was proposed to act similarly in *C. striatus* [36]. While some sequences in our transcriptome with the framework IX were similar to *C. californicus* sequences that contained a kunitz-motif [37], a search across the Conserved Domain Database retrieved no hits for the kunitz-domain and no other sequences with similarities to κ-conotoxins were identified. While the presence of highly divergent δ-, μ- or κ-conotoxins in *C. distans* cannot be excluded, it appears that these
would be no more than minor components of the venom. Given there are intraspecific variations within cone snail species, sequencing the venom gland transcriptomes of more specimens of *C. distans* would help to further elucidate intraspecific differences in venom expression. Nonetheless, our data suggests early diverging cone snails expressed sequences ancestral to modern α- and ω-conotoxins found in mollusc- and fish-hunting cone snails. Figure 6 provides an overview of how *C. distans'* venom gland transcriptome compares with other cone snail species when placed within its phylogenetic context, as determined by Puillandre et al. [14].

<figure><img src="image_14.png" /><figcaption>**Figure 6.** Overview of the major venom components of Conidae. The phylogenetic reconstruction was adapted from Puillandre et al. [14]. Lineages in green indicate mollusc-hunting species and lineages in red indicate fish-hunting species. All other lineages are predominantly comprised of worm-hunters except for divergent species such as *Conus californicus*, which can prey on fish, worms and molluscs.</figcaption></figure>

Previous screening attempts using *C. distans* venom have identified few conotoxins with potent biological activity. Worm, mollusc and/or fish assays specifically designed to identify molecules targeting ecologically relevant species are expected to identify more bioactive peptides. Interestingly, *C. distans* venom is unique in having a large proportion of low molecular weight molecules <1 kDa. The total ion chromatograms (TICs) from LC-ESI-MS experiments from both *C. distans* specimens did not show any major peaks despite numerous attempts to optimise the sample runs with different LC-MS conditions. Unfortunately, we were unable to elicit a defensive response from *C. distans* when stimulated by a predator such as *C. marmoreus* or *C. textile*, as previously demonstrated in another early diverging species, *C. imperialis*. Given the relatively large size and shell strength of *C. distans*, the shell alone might offer sufficient protection from most predators for this species and *C. distans* may potentially lack a defensive venom. We were also unable to find suitable prey for *C. distans* to collect the predation-evoked venom to establish if the full complement of venom peptides, or only a specific subset, are expressed during predation. The low expression levels of masses that correspond to the typical conotoxin mass range
raises the possibility that *C. distans* may have evolved defensive and predatory strategies in
the absence of obvious venom gland compartmentalisation.

To further understand the evolution of the venom gland within the genus *Conus*, we compared *C. distans* with *C. planorbis* and *C. geographus*. *C. planorbis* is the earliest diverging species from which a defence-evoked venom was obtained, notably in the presence of *C. marmoreus*, a molluscivore cone snail [12]. Uniquely among species from which defensive stings have been collected, expression of defensive venom peptides in *C. planorbis* was not restricted to the proximal section [12]. We used *C. geographus* as a reference species, as the use of defensive and predatory venoms in this species has been well established [7,38]. Our comparison here using PCA analysis showed that the expression of peptides across the venom gland of *C. planorbis* is more compartmentalised than *C. distans* but less compartmentalised than *C. geographus*. This suggests that the venom gland of *C. planorbis* represents an evolutionary stage that is intermediate between the non-compartmentalised *C. distans* venom gland and the highly compartmentalised glands observed in more recently evolved worm and fish-hunting species [13].

**4. Materials and Methods**

**4.1. Sample Collection, RNA Extraction and Sequencing**

Specimens of *C. distans* were collected from Gould Reef in the central Great Barrier Reef (Queensland, Australia) and maintained alive in marine aquaria before use. One specimen was carefully dissected on ice and the venom gland placed in an Eppendorf tube containing 1 mL of TRIZOL (Invitrogen, Carlsbad, CA, USA) and total RNA extracted following the manufacturer's instructions. Next, mRNA was purified using the Oligotex mRNA kit and ~200 ng was sequenced by the Australian Genome Research Facility (AGRF, Brisbane) on a Roche 454-GS FLX plus titanium sequencer. Sequencing yielded 209,346 reads with an average read length of ~350 bp after filtering poor quality reads.

4.2. *Transcriptomic Analysis*

Transcriptome analysis was performed as previously described [26]. Briefly, raw reads were filtered based on quality scores (Quality >30; 1 in 1000 incorrect base call frequency) and sequences possessing an archetypal conotoxin precursor structure [39] were sorted into superfamilies using Conosorter. Unclassified sequences with >50 A.A., a read count ≥2, a class score >1, and superfamily score >0 were also included. Reads with incomplete signal sequences, unrecognised amino acids, frameshifts or truncations, and duplicated sequences were discarded. After discarding housekeeping proteins, the remaining sequences were clustered by their signal region and classified using BLASTp. ConoPrec [40] was then used to place conotoxin-like sequences into known superfamilies or to designate novel Superfamilies based on a signal sequence similarity cut-off of 53.3% [26]. Sequences were then manually inspected, aligned and visualised using Ugene [41].

4.3. *Sample Collection for Mass Spectrometry*

*C. distans* specimen 1 was collected from Lady Musgrave Island on the southern Great Barrier Reef and *C. distans* specimen 2 was collected off Cairns in the northern Great Barrier Reef. Both were maintained in an aquarium with a 12 h day/night cycle until use. All three *C. planorbis* specimens and *C. geographus* were collected off Cairns in the northern Great Barrier Reef.

4.4. *Venom Extraction*

Specimens were dissected on ice and the venom gland removed intact and unravelled
carefully to avoid breakage. The gland was then divided into four equal lengths of distal
(near the proboscis), distal central, proximal central and proximal (near the venom bulb)
venom duct. Venom from each section was stripped into 500 µL of 30% acetonitrile
(ACN)/0.1% formic acid (FA), vortexed for 1 min, sonicated for 30 s and centrifuged on a
bench-top machine at 12,000× g for 10 min. The supernatant was collected, and the pellet
was subjected to two additional rounds of extraction, as described above. Supernatants
were pooled, lyophilised overnight, and stored at -20 °C prior to use.

4.5. *Reduction, Alkylation and Enzymatic Digest of Venom Samples*

Powdered venoms were resuspended in 1 M (NH₄)₂HCO₃ at pH 11 (1 µg/µL) and reduction and alkylation solution (97.5% ACN, 2% iodoethanol, 0.5% triethylphospine) added in a 1:1 ratio, followed by incubation at 37 °C for 2 h [42]. MS-grade trypsin and endoproteinase-GluC (Promega, Madison, WI, USA) dissolved in 100 µL of 40 mM NH₄HCO₃ were then added to speed-vac-dried sample at a protease-to-protein ratio of 1:100 and incubated overnight at 37 °C. Digested samples were dried on a speed-vac prior to use.

4.6. *Venom Gland Mass Spectrometry (LC-ESI-MS and MS/MS)*

LC-ESI-MS and MS/MS were performed on native and reduced-alkylated-digested
samples using a 5600 tripleTOF mass spectrometer with a quadruple TOF system equipped
with a TurboV DuoSpray (ABSCIEX, Ontario, Canada) source set to 5300 V and 500 °C.
Native and processed venoms were dried to remove buffer, resuspended in 1% formic acid
(FA), vortexed and centrifuged at 12,000 × g for 10 min to remove particulates. Samples were
separated on a ZORBAX 300SB-C18 (2.1mm × 100 mm × 1.8 µm) column using a Shimadzu
30 series HPLC system (Shimadzu, Kyoto, Japan). For LC-ESI-MS experiments on C. *distans*
specimen 2, C. *planorbis* and C. *geographus* samples, samples were eluted at 0.2 mL/min
with a gradient of 1% Solvent A (0.1% FA) to 40% Solvent B (40% acetonitrile/0.1% FA) over
80 min. Full-scan MS data was accumulated over 500 ms with cycle time set to 0.525 s and
a mass range of 350–2400 Da. For LC-ESI-MS/MS for C. *distans* specimen 1, samples were
eluted at 180 µL/min with a linear gradient of 1–40% over 45 min. The columns were then
flushed over 15 min with a linear gradient increasing to 98% buffer B. Separated samples
were injected using the TurboV ion spray source with the ionspray voltage of the source
set to 5300 V and the temperature to 450 °C. Full-scanning MS data was accumulated over
100 ms followed by full-scan product ion data in the high-resolution information dependent
acquisition (IDA) mode. The total duration of the full-scan MS was 60 min, consisting of
1645 cycles of 2.15 s per cycle of full-scan MS. A rolling collision energy based on precursor
ion *m/z* was used to generate product ions for the TOF MS/MS. The mass range was set
to 300–1800 (*m/z*) for TOF MS mode and 80–1400 (*m/z*) for full-scan TOF MS/MS mode.
Product ions with charge states between 2 and 5, with a minimum intensity of 150 cps,
were recorded, and isotopes within 4 Da were binned, though former target ions were
retained. To capture the mass range of undigested venom samples more thoroughly, with
a particular focus on higher molecular weight components in the venom, we performed
LC-ESI-MS experiments on native C. *distans* specimen 1. The LC gradient was set to 2–40%
B over 30 min with a flow rate of 0.5 µL/min. Full-scan MS was performed over 54 min,
comprising of 3161 cycles of 1.002 s per cycle. Full-scan TOF MS data was obtained over a
mass range of 400–2000 Da and accumulated over 1 s.

4.7. *MS and MS/MS Data Analysis*

MS data analysis was performed using Analyst 1.6 and Peakview 2.1 (ABSCIEX,
Ontario, Canada). Masses were reconstructed using the LC-MS reconstruct option with
a mass range of 400–12,000 Da, a tolerance of 0.2 Da and the S/N ratio threshold of 10.
Masses of major peaks were confirmed manually. MS/MS data analysis was performed
using ProteinPilot 4.1. Briefly, theoretical spectra generated from transcriptomic sequences
were matched to MS/MS spectra from proteomics using the paragon algorithm. Only
sequences assigned to spectra at a confidence value >99% were analysed.

4.8. *Principal Component Analysis (PCA) of LC-ESI-MS Data*

Unprocessed LC-ESI-MS data were imported into MarkerView™ Software (v1.3.1) to
generate a list of unique masses for each venom gland section. Peaks were detected using a
peak width of <100 scans, a noise threshold of 10 and retention time between 10–65 min, with a maximum of 8000 peaks generated for each dataset. Peaks from different samples were aligned and filtered using a retention time tolerance of 0.5 min and mass tolerance of 25 ppm. Unsupervised PCA using Pareto scaling was performed on monoisotopic masses >400 Da. The generated PC1 and PC2 score and loading values were plotted on GraphPad Prism (v7).

**5. Conclusions**

In conclusion, our study of the basal lineage cone snail species *C. distans* revealed it has an undifferentiated venom duct and likely uses a single venom for defence and prey capture, as previously hypothesised for ancestral species [7]. Transcriptomic analysis identified novel α- and ω-conotoxin-like sequences similar to conotoxins associated with predation in piscivores, suggesting that these sequences might represent ancestral defensive venom in this worm-hunting species. As the venom gland began to compartmentalise, as evidenced by the proto-compartmentalised gland seen in *C. planorbis*, the expression of these sequences was likely localised to discrete duct sections, facilitating the separate evolution of predatory and defensive venoms. Our study supports a pivotal role of venom gland compartmentalisation in guiding the evolution of venom peptides that allowed mollusc- and fish-hunting to evolve in cone snails.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/toxins14030226/s1, Figure S1: LC-ESI-MS analysis of the four duct sections from *C. distans* specimen 2.; Figure S2: O1 and O2 superfamilies in *C. distans*; Figure S3: O1 and O2 superfamilies in *C. distans*; Figure S4: Conopeptides (A) and divergent superfamilies (B) in *C. distans*; Figure S5: Novel superfamilies in *C. distans*. (A) and (B) indicate novel superfamilies reported only in *C. distans* and novel superfamilies also reported in other species.

**Author Contributions:** Conceptualization, R.J.L. and S.D.; methodology, J.R.P. and S.K.R.; investigation, J.R.P. and S.K.R.; resources, R.J.L.; writing—original draft preparation, J.R.P.; writing—review and editing, R.J.L. and S.D.; visualization, J.R.P., S.K.R. and R.J.L.; supervision, R.J.L. and S.D.; project administration, R.J.L.; funding acquisition, R.J.L. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by an Australian Research Council grant DP170104792 (to RJL and SD), and a National Health and Medical Research Council grant APP1119056 and research fellowship APP1119056 (to RJL).

**Data Availability Statement:** All available data are published in the manuscript and Supplementary Data.

**Acknowledgments:** We thank Alun Jones for mass spectrometry support and AGRF for sequencing.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Casewell, N.R.; Wüster, W.; Vonk, F.J.; Harrison, R.A.; Fry, B.G. Complex cocktails: The evolutionary novelty of venoms. *Trends Ecol. Evol.* **2013**, *28*, 219–229. [CrossRef]

2. Savitzky, A.H. The role of venom delivery strategies in snake evolution. *Evolution* **1980**, *34*, 1194–1204. [CrossRef]

3. Undheim, E.A.; Hamilton, B.R.; Kurniawan, N.D.; Bowlay, G.; Cribb, B.W.; Merritt, D.J.; Fry, B.G.; King, G.F.; Venter, D.J. Production and packaging of a biological arsenal: Evolution of centipede venoms under morphological constraint. *Proc. Natl. Acad. Sci. USA* **2015**, *112*, 4026–4031. [CrossRef]

4. Fry, B.G.; Scheib, H.; van der Weerd, L.; Young, B.; McNaughtan, J.; Ramjan, S.F.R.; Vidal, N.; Poelmann, R.E.; Norman, J.A. Evolution of an arsenal: Structural and functional diversification of the venom system in the advanced snakes (Caenophidia). *Mol. Cell. Proteom.* **2008**, *7*, 215–246. [CrossRef]

5. Li, M.; Fry, B.G.; Kini, R.M. Eggs-only diet: Its implications for the toxin profile changes and ecology of the marbled sea snake (*Aipysurus eydouxii*). *J. Mol. Evol.* **2005**, *60*, 81–89. [CrossRef]

6. Wright, J.J. Adaptive significance of venom glands in the tadpole madtom *Noturus gyrinus* (Siluriformes: *Ictaluridae*). *J. Exp. Biol.* **2012**, *215*, 1816–1823. [CrossRef]

7. Dutertre, S.; Jin, A.-H.; Vetter, I.; Hamilton, B.; Sunagar, K.; Lavergne, V.; Dutertre, V.; Fry, B.G.; Antunes, A.; Venter, D.J. Evolution of separate predation-and defence-evoked venoms in carnivorous cone snails. *Nat. Commun.* **2014**, *5*, 3521. [CrossRef]
8. Walker, A.A.; Mayhew, M.L.; Jin, J.; Herzig, V.; Undheim, E.A.B.; Sombke, A.; Fry, B.G.; Meritt, D.J.; King, G.F. The assassin bug *Pristhesancus plagipennis* produces two distinct venoms in separate gland lumens. *Nat. Commun.* **2018**, *9*, 755. [CrossRef]

9. Puillandre, N.; Duda, T.F.; Meyer, C.; Olivera, B.M.; Bouchet, P. One, four or 100 genera? A new classification of the cone snails. *J. Molluscan Stud.* **2015**, *81*, 1–23. [CrossRef]

10. Julita, S.I.; Natalie, S.; Olivera, B.M.; Bandyopadhyay, P.K.; Annett, S.; Ferber, M.; Heinrich, T. Using chemistry to reconstruct evolution: On the origins of fish-hunting in venomous cone snails. *Proc. Am. Philos. Soc.* **2007**, *151*, 185–200.

11. Dutertre, S.; Griffin, J.; Lewis, R.J. Phyla Molluska: The venom apparatus of cone snails. In *Marine and Freshwater Toxins*; Gopalakrishnakone, P., Haddad, J.V., Tubaro, A., Kim, E., Kem, W.R., Eds.; Springer: Dordrecht, The Netherlands, 2016; pp. 327–340. [CrossRef]

12. Jin, A.H.; Vetter, I.; Himaya, S.W.; Alewood, P.F.; Lewis, R.J.; Dutertre, S. Transcriptome and proteome of *Conus planorbis* identify the nicotinic receptors as primary target for the defensive venom. *Proteomics* **2015**, *15*, 4030–4040. [CrossRef] [PubMed]

13. Prashanth, J.R.; Dutertre, S.; Jin, A.H.; Lavergne, V.; Hamilton, B.; Cardoso, F.C.; Griffin, J.; Venter, D.J.; Alewood, P.F.; Lewis, R.J. The role of defensive ecological interactions in the evolution of conotoxins. *Mol. Ecol.* **2016**, *25*, 598–615. [CrossRef]

14. Puillandre, N.; Bouchet, P.; Duda, T.; Kauferstein, S.; Kohn, A.; Olivera, B.; Watkins, M.; Meyer, C. Molecular phylogeny and evolution of the cone snails (Gastropoda, Conoidea). *Mol. Phylogenet. Evol.* **2014**, *78*, 290–303. [CrossRef]

15. Saxena, V.; Partoens, P.; De Block, J.; Coen, E.; Vauquelin, G.; De Potter, W. Inhibition of evoked neurotransmitter release from rat hippocampus by a polypeptide toxin isolated from the marine snail *Conus distans*. *Neurochem. Int.* **1992**, *20*, 69–74. [CrossRef]

16. Partoens, P.; Wang, J.; Coen, E.; Vauquelin, G.; De Potter, W. Two polypeptide toxins with opposite effects on calcium uptake in bovine chromaffin cells: Isolation from the venom of the marine snail *Conus distans*. *Neurochem. Int.* **1996**, *28*, 619–624. [CrossRef]

17. Chen, P.; Garrett, J.E.; Watkins, M.; Olivera, B.M. Purification and characterization of a novel excitatory peptide from *Conus distans* venom that defines a novel gene superfamily of conotoxins. *Toxicon* **2008**, *52*, 139–145. [CrossRef]

18. Neves, J.L.B.; Lin, Z.; Imperial, J.S.; Antunes, A.; Vasconcelos, V.; Olivera, B.M.; Schmidt, E.W. Small molecules in the cone snail arsenal. *Org. Lett.* **2015**, *17*, 4933–4935. [CrossRef] [PubMed]

19. Biggs, J.S.; Watkins, M.; Puillandre, N.; Ownby, J.-P.; Lopez-Vera, E.; Christensen, S.; Moreno, K.J.; Bernaldez, J.; Licea-Navarro, A.; Corneli, P.S.; et al. Evolution of *Conus* peptide toxins: Analysis of *Conus californicus* Reeve, 1844. *Mol. Phylogenet. Evol.* **2010**, *56*, 1–12. [CrossRef] [PubMed]

20. Miles, L.A.; Dy, C.Y.; Nielsen, J.; Barnham, K.J.; Hinds, M.G.; Olivera, B.M.; Bulaj, G.; Norton, R.S. Structure of a novel P-superfamily spasmodic conotoxin reveals an inhibitory cystine knot motif. *J. Biol. Chem.* **2002**, *277*, 43033–43040. [CrossRef]

21. Jin, A.-H.; Dutertre, S.; Kaas, Q.; Lavergne, V.; Kubala, P.; Lewis, R.J.; Alewood, P.F. Transcriptomic messiness in the venom duct of *Conus miles* contributes to conotoxin diversity. *Mol. Cell. Proteom.* **2013**, *12*, 3824–3833. [CrossRef]

22. Lirazan, M.B.; Hooper, D.; Corpuz, G.P.; Ramilo, C.A.; Bandyopadhyay, P.; Cruz, L.J.; Olivera, B.M. The spasmodic peptide defines a new conotoxin superfamily. *Biochemistry* **2000**, *39*, 1583–1588. [CrossRef] [PubMed]

23. England, L.J.; Imperial, J.; Jacobsen, R.; Craig, A.G.; Gulyas, J.; Akhtar, M.; Rivier, J.; Julius, D.; Olivera, B.M. Inactivation of a serotonin-gated ion channel by a polypeptide toxin from marine snails. *Science* **1998**, *281*, 575–578. [CrossRef]

24. Teichert, R.W.; Jimenez, E.C.; Olivera, B.M. αS-conotoxin RVIII: A structurally unique conotoxin that broadly targets nicotinic acetylcholine receptors. *Biochemistry* **2005**, *44*, 7897–7902. [CrossRef]

25. Himaya, S.W.A.; Mari, F.; Lewis, R.J. Accelerated proteomic visualization of individual predatory venoms of *Conus purpurascens* reveals separately evolved predation-evoked venom cabals. *Sci. Rep.* **2018**, *8*, 330. [CrossRef]

26. Prashanth, J.R.; Lewis, R.J. An efficient transcriptome analysis pipeline to accelerate venom peptide discovery and characterisation. *Toxicon* **2015**, *107*, 282–289. [CrossRef]

27. Mena, E.E.; Gullak, M.F.; Pagnozzi, M.J.; Richter, K.E.; Rivier, J.; Cruz, L.J.; Olivera, B.M. Conantokin-G: A novel peptide antagonist to the N-methyl-D-aspartic acid (NMDA) receptor. *Neurosci. Lett.* **1990**, *118*, 241–244. [CrossRef]

28. Walker, C.S.; Jensen, S.; Ellison, M.; Matta, J.A.; Lee, W.Y.; Imperial, J.S.; Duclos, N.; Brockie, P.J.; Madsen, D.M.; Isaac, J.T. A Novel *Conus* snail polypeptide causes excitotoxicity by blocking desensitization of AMPA receptors. *Curr. Biol.* **2009**, *19*, 900–908. [CrossRef]

29. Jimenéz, E.C.; Olivera, B.M.; Gray, W.R.; Cruz, L.J. Contryphan is a D-tryptophan-containing *Conus* peptide. *J. Biol. Chem.* **1996**, *271*, 28002–28005. [CrossRef]

30. Lewis, R.J.; Dutertre, S.; Vetter, I.; Christie, M.D.J. *Conus* venom peptide pharmacology. *Pharmacol. Rev.* **2012**, *64*, 259–298. [CrossRef]

31. Cruz, L.; Gray, W.; Olivera, B.; Zeikus, R.; Kerr, L.; Yoshikami, D.; Moczydlowski, E. *Conus geographus* toxins that discriminate between neuronal and muscle sodium channels. *J. Biol. Chem.* **1985**, *260*, 9280–9288. [CrossRef]

32. McIntosh, J.M.; Hasson, A.; Spira, M.E.; Gray, W.R.; Li, W.; Marsh, M.; Hillyard, D.R.; Olivera, B.M. A new family of conotoxins that blocks voltage-gated sodium channels. *J. Biol. Chem.* **1995**, *270*, 16796–16802. [CrossRef] [PubMed]

33. Liu, J.; Wu, Q.; Pi, C.; Zhao, Y.; Zhou, M.; Wang, L.; Chen, S.; Xu, A. Isolation and characterization of a T-superfamily conotoxin from *Conus litteratus* with targeting tetrodotoxin-sensitive sodium channels. *Peptides* **2007**, *28*, 2313–2319. [CrossRef] [PubMed]

34. Zakon, H.H. Convergent evolution on the molecular level. *Brain Behav. Evol.* **2002**, *59*, 250–261. [CrossRef]

35. Terlau, H.; Shon, K.-J.; Grilley, M.; Stocker, M.; Stuehmer, W.; Olivera, B.M. Strategy for rapid immobilization of prey by a fish-hunting marine snail. *Nature* **1996**, *381*, 148–151. [CrossRef]
36. Bayrhuber, M.; Vijayan, V.; Ferber, M.; Graf, R.; Korukottu, J.; Imperial, J.; Garrett, J.E.; Olivera, B.M.; Terlau, H.; Zweckstetter, M.; et al. Conkunitzin-S1is the first member of a new Kunitz-type neurotoxin family: Structural and functional characterization. *J. Biol. Chem.* **2005**, *280*, 23766–23770. [CrossRef]
37. Elliger, C.A.; Richmond, T.A.; Lebaric, Z.N.; Pierce, N.T.; Sweedler, J.V.; Gilly, W.F. Diversity of conotoxin types from Conus californicus reflects a diversity of prey types and a novel evolutionary history. *Toxicon* **2011**, *57*, 311–322. [CrossRef] [PubMed]
38. Dutertre, S.; Jin, A.-H.; Alewood, P.F.; Lewis, R.J. Intraspecific variations in *Conus geographus* defence-evoked venom and estimation of the human lethal dose. *Toxicon* **2014**, *91*, 135–144. [CrossRef]
39. Lavergne, V.; Dutertre, S.; Jin, A.-H.; Lewis, R.J.; Taft, R.J.; Alewood, P.F. Systematic interrogation of the *Conus marmoreus* venom duct transcriptome with ConoSorter reveals 158 novel conotoxins and 13 new gene superfamilies. *BMC Genom.* **2013**, *14*, 708. [CrossRef]
40. Kaas, Q.; Yu, R.; Jin, A.-H.; Dutertre, S.; Craik, D.J. ConoServer: Updated content, knowledge, and discovery tools in the conopeptide database. *Nucleic Acids Res.* **2012**, *40*, D325–D330. [CrossRef] [PubMed]
41. Okonechnikov, K.; Golosova, O.; Fursov, M. Unipro UGENE: A unified bioinformatics toolkit. *Bioinformatics* **2012**, *28*, 1166–1167. [CrossRef]
42. Hale, J.E.; Butler, J.P.; Gelfanova, V.; You, J.-S.; Knierman, M.D. A simplified procedure for the reduction and alkylation of cysteine residues in proteins prior to proteolytic digestion and mass spectral analysis. *Anal. Biochem.* **2004**, *333*, 174–181. [CrossRef] [PubMed]