Divergent venom effectors correlate with ecological niche in neuropteran predators.

Neuropteran larvae are fierce predators that use venom to attack and feed on arthropod prey. Neuropterans have adapted to diverse and sometimes extreme habitats, suggesting their venom may have evolved accordingly, but the ecology and evolution of venom deployment in different families is poorly understood. We applied spatial transcriptomics, proteomics, morphological analysis, and bioassays to investigate the venom systems in the antlion Euroleon nostras and the lacewing Chrysoperla carnea, which occupy distinct niches. Although the venom system morphology was similar in both species, we observed remarkable differences at the molecular level. E. nostras produces particularly complex venom secreted from three different glands, indicating functional compartmentalization. Furthermore, E. nostras venom and digestive tissues were devoid of bacteria, strongly suggesting that all venom proteins are of insect origin rather than the products of bacterial symbionts. We identified several toxins exclusive to E. nostras venom, including phospholipase A2 and several undescribed proteins with no homologs in the C. carnea genome. The compositional differences have significant ecological implications because only antlion venom conferred insecticidal activity, indicating its use for the immobilization of large prey. Our results indicate that molecular venom evolution plays a role in the adaptation of antlions to their unique ecological niche.

Selection Across the Three-Dimensional Structure of Venom Proteins from North American Scolopendromorph Centipedes.

Gene duplication followed by nucleotide differentiation is one of the simplest mechanisms to develop new functions for genes. However, the evolutionary processes underlying the divergence of multigene families remain controversial. We used multigene families found within the diversity of toxic proteins in centipede venom to test two hypotheses related to venom evolution: the two-speed mode of venom evolution and the rapid accumulation of variation in exposed residues (RAVER) model. The two-speed mode of venom evolution proposes that different types of selection impact ancient and younger venomous lineages with negative selection being the predominant form in ancient lineages and positive selection being the dominant form in younger lineages. The RAVER hypothesis proposes that, instead of different types of selection acting on different ages of venomous lineages, the different types of selection will selectively contribute to amino acid variation based on whether the residue is exposed to the solvent where it can potentially interact directly with toxin targets. This hypothesis parallels the longstanding understanding of protein evolution that suggests that residues found within the structural or active regions of the protein will be under negative or purifying selection, and residues that do not form part of these areas will be more prone to positive selection. To test these two hypotheses, we compared the venom of 26 centipedes from the order Scolopendromorpha from six currently recognized species from across North America using both transcriptomics and proteomics. We first estimated their phylogenetic relationships and uncovered paraphyly among the genus Scolopendra and evidence for cryptic diversity among currently recognized species. Using our phylogeny, we then characterized the diverse venom components from across the identified clades using a combination of transcriptomics and proteomics. We conducted selection-based analyses in the context of predicted three-dimensional properties of the venom proteins and found support for both hypotheses. Consistent with the two-speed hypothesis, we found a prevalence of negative selection across all proteins. Consistent with the RAVER hypothesis, we found evidence of positive selection on solvent-exposed residues, with structural and less-exposed residues showing stronger signal for negative selection. Through the use of phylogenetics, transcriptomics, proteomics, and selection-based analyses, we were able to describe the evolution of venom from an ancient venomous lineage and support principles of protein evolution that directly relate to multigene family evolution.

Phylogeny, envenomation syndrome, and membrane permeabilising venom produced by Australia's electric caterpillar Comana monomorpha.

Zygaenoidea is a superfamily of lepidopterans containing many venomous species, including the Limacodidae (nettle caterpillars) and Megalopygidae (asp caterpillars). Venom proteomes have been recently documented for several species from each of these families, but further data are required to understand the evolution of venom in Zygaenoidea. In this study, we examined the 'electric' caterpillar from North-Eastern Australia, a limacodid caterpillar densely covered in venomous spines. We used DNA barcoding to identify this caterpillar as the larva of the moth Comana monomorpha (Turner, 1904). We report the clinical symptoms of C. monomorpha envenomation, which include acute pain, and erythema and oedema lasting for more than a week. Combining transcriptomics of venom spines with proteomics of venom harvested from the spine tips revealed a venom markedly different in composition from previously examined limacodid venoms that are rich in peptides. In contrast, the venom of C. monomorpha is rich in aerolysin-like proteins similar to those found in venoms of asp caterpillars (Megalopygidae). Consistent with this composition, the venom potently permeabilises sensory neurons and human neuroblastoma cells. This study highlights the diversity of venom composition in Limacodidae.

Venom resistance mechanisms in centipede show tissue specificity.

Venomous animals utilize venom glands to secrete and store powerful toxins for intraspecific and/or interspecific antagonistic interactions, implying that tissue-specific resistance is essential for venom glands to anatomically separate toxins from other tissues. Here, we show the mechanism of tissue-specific resistance in centipedes (Scolopendra subspinipes mutilans), where the splice variant of the receptor repels its own toxin. Unlike the well-known resistance mechanism by mutation in a given exon, we found that the KCNQ1 channel is highly expressed in the venom gland as a unique splice variant in which the pore domain and transmembrane domain six, partially encoded by exon 6 (rather than 7 as found in other tissues), contain eleven mutated residues. Such a splice variant is sufficient to gain resistance to SsTx (a lethal toxin for giant prey capture) in the venom gland due to a partially buried binding site. Therefore, the tissue-specific KCNQ1 modification confers resistance to the toxins, establishing a safe zone in the venom-storing/secreting environment.

Production, composition, and mode of action of the painful defensive venom produced by a limacodid caterpillar, Doratifera vulnerans.

Venoms have evolved independently several times in Lepidoptera. Limacodidae is a family with worldwide distribution, many of which are venomous in the larval stage, but the composition and mode of action of their venom is unknown. Here, we use imaging technologies, transcriptomics, proteomics, and functional assays to provide a holistic picture of the venom system of a limacodid caterpillar, Doratifera vulnerans Contrary to dogma that defensive venoms are simple in composition, D. vulnerans produces a complex venom containing 151 proteinaceous toxins spanning 59 families, most of which are peptides <10 kDa. Three of the most abundant families of venom peptides (vulnericins) are 1) analogs of the adipokinetic hormone/corazonin-related neuropeptide, some of which are picomolar agonists of the endogenous insect receptor; 2) linear cationic peptides derived from cecropin, an insect innate immune peptide that kills bacteria and parasites by disrupting cell membranes; and 3) disulfide-rich knottins similar to those that dominate spider venoms. Using venom fractionation and a suite of synthetic venom peptides, we demonstrate that the cecropin-like peptides are responsible for the dominant pain effect observed in mammalian in vitro and in vivo nociception assays and therefore are likely to cause pain after natural envenomations by D. vulnerans Our data reveal convergent molecular evolution between limacodids, hymenopterans, and arachnids and demonstrate that lepidopteran venoms are an untapped source of novel bioactive peptides.

Phylogenetic analyses suggest centipede venom arsenals were repeatedly stocked by horizontal gene transfer.

Venoms have evolved over a hundred times in animals. Venom toxins are thought to evolve mostly by recruitment of endogenous proteins with physiological functions. Here we report phylogenetic analyses of venom proteome-annotated venom gland transcriptome data, assisted by genomic analyses, to show that centipede venoms have recruited at least five gene families from bacterial and fungal donors, involving at least eight horizontal gene transfer events. These results establish centipedes as currently the only known animals with venoms used in predation and defence that contain multiple gene families derived from horizontal gene transfer. The results also provide the first evidence for the implication of horizontal gene transfer in the evolutionary origin of venom in an animal lineage. Three of the bacterial gene families encode virulence factors, suggesting that horizontal gene transfer can provide a fast track channel for the evolution of novelty by the exaptation of bacterial weapons into animal venoms.

Pick Your Poison: Molecular Evolution of Venom Proteins in Asilidae (Insecta: Diptera).

Robber flies are an understudied family of venomous, predatory Diptera. With the recent characterization of venom from three asilid species, it is possible, for the first time, to study the molecular evolution of venom genes in this unique lineage. To accomplish this, a novel whole-body transcriptome of Eudioctria media was combined with 10 other publicly available asiloid thoracic or salivary gland transcriptomes to identify putative venom gene families and assess evidence of pervasive positive selection. A total of 348 gene families of sufficient size were analyzed, and 33 of these were predicted to contain venom genes. We recovered 151 families containing homologs to previously described venom proteins, and 40 of these were uniquely gained in Asilidae. Our gene family clustering suggests that many asilidin venom gene families are not natural groupings, as delimited by previous authors, but instead form multiple discrete gene families. Additionally, robber fly venoms have relatively few sites under positive selection, consistent with the hypothesis that the venoms of older lineages are dominated by negative selection acting to maintain toxic function.

Weaponisation 'on the fly': Convergent recruitment of knottin and defensin peptide scaffolds into the venom of predatory assassin flies.

Many arthropod venom peptides have potential as bioinsecticides, drug leads, and pharmacological tools due to their specific neuromodulatory functions. Assassin flies (Asilidae) are a family of predaceous dipterans that produce a unique and complex peptide-rich venom for killing insect prey and deterring predators. However, very little is known about the structure and function of their venom peptides. We therefore used an E. coli periplasmic expression system to express four disulfide-rich peptides that we previously reported to exist in venom of the giant assassin fly Dolopus genitalis. After purification, each recombinant peptide eluted from a C18 column at a position closely matching its natural counterpart, strongly suggesting adoption of the native tertiary fold. Injection of purified recombinant peptides into blowflies (Lucilia cuprina) and crickets (Acheta domestica) revealed that two of the four recombinant peptides, named rDg3b and rDg12, inhibited escape behaviour in a manner that was rapid in onset (<1 min) and reversible. Homonuclear NMR solution structures revealed that rDg3b and rDg12 adopt cystine-stabilised α/ß defensin and inhibitor cystine knot folds, respectively. Although the closest known homologues of rDg3b at the level of primary structure are dipteran antimicrobial peptides such as sapecin and lucifensin, a DALI search showed that the tertiary structure of rDg3b most closely resembles the KV11.1-specific α-potassium channel toxin CnErg1 from venom of the scorpion Centruroides noxius. This is mainly due to the deletion of a large, unstructured loop between the first and second cysteine residues present in Dg3b homologues from non-asiloid, but not existing in asiloid, species. Patch-clamp electrophysiology experiments revealed that rDg3b shifts the voltage-dependence of KV11.1 channel activation to more depolarised potentials, but has no effect on KV1.3, KV2.1, KV10.1, KCa1.1, or the Drosophila Shaker channel. Although rDg12 shares the inhibitor cystine knot structure of many gating modifier toxins, rDg12 did not affect any of these KV channel subtypes. Our results demonstrate that multiple disulfide-rich peptide scaffolds have been convergently recruited into asilid and other animal venoms, and they provide insight into the molecular evolution accompanying their weaponisation.

Bioinformatic analysis suggests potential mechanisms underlying parasitoid venom evolution and function.

Hymenopteran parasitoid wasps are a diverse collection of species that infect arthropod hosts and use factors found in their venoms to manipulate host immune responses, physiology, and behaviour. Whole parasitoid venoms have been profiled using proteomic approaches, and here we present a bioinformatic characterization of the venom protein content from Ganaspis sp. 1, a parasitoid that infects flies of the genus Drosophila. We find evidence that diverse evolutionary processes including multifunctionalization, co-option, gene duplication, and horizontal gene transfer may be acting in concert to drive venom gene evolution in Ganaspis sp.1. One major role of parasitoid wasp venom is host immune evasion. We previously demonstrated that Ganaspis sp. 1 venom inhibits immune cell activation in infected Drosophila melanogaster hosts, and our current analysis has uncovered additional predicted virulence functions. Overall, this analysis represents an important step towards understanding the composition and activity of parasitoid wasp venoms.

The Latoia consocia Caterpillar Induces Pain by Targeting Nociceptive Ion Channel TRPV1.

Accidental contact with caterpillar bristles causes local symptoms such as severe pain, intense heat, edema, erythema, and pruritus. However, there is little functional evidence to indicate a potential mechanism. In this study, we analyzed the biological characteristics of the crude venom from the larval stage of Latoia consocia living in South-West China. Intraplantar injection of the venom into the hind paws of mice induced severe acute pain behaviors in wild type (WT) mice; the responses were much reduced in TRPV1-deficit (TRPV1 KO) mice. The TRPV1-specific inhibitor, capsazepine, significantly attenuated the pain behaviors. Furthermore, the crude venom evoked strong calcium signals in the dorsal root ganglion (DRG) neurons of WT mice but not those of TRPV1 KO mice. Among the pain-related ion channels we tested, the crude venom only activated the TRPV1 channel. To better understand the venom components, we analyzed the transcriptome of the L. consocia sebaceous gland region. Our study suggests that TRPV1 serves as a primary nociceptor in caterpillar-induced pain and forms the foundation for elucidating the pain-producing mechanism.

Missiles of Mass Disruption: Composition and Glandular Origin of Venom Used as a Projectile Defensive Weapon by the Assassin Bug Platymeris rhadamanthus.

Assassin bugs (Reduviidae) produce venoms that are insecticidal, and which induce pain in predators, but the composition and function of their individual venom components is poorly understood. We report findings on the venom system of the red-spotted assassin bug Platymeris rhadamanthus, a large species of African origin that is unique in propelling venom as a projectile weapon when threatened. We performed RNA sequencing experiments on venom glands (separate transcriptomes of the posterior main gland, PMG, and the anterior main gland, AMG), and proteomic experiments on venom that was either defensively propelled or collected from the proboscis in response to electrostimulation. We resolved a venom proteome comprising 166 polypeptides. Both defensively propelled venom and most venom samples collected in response to electrostimulation show a protein profile similar to the predicted secretory products of the PMG, with a smaller contribution from the AMG. Pooled venom samples induce calcium influx via membrane lysis when applied to mammalian neuronal cells, consistent with their ability to cause pain when propelled into the eyes or mucus membranes of potential predators. The same venom induces rapid paralysis and death when injected into fruit flies. These data suggest that the cytolytic, insecticidal venom used by reduviids to capture prey is also a highly effective defensive weapon when propelled at predators.

A family of orthologous proteins from centipede venoms inhibit the hKir6.2 channel.

Inhibitors targeting ion channels are useful tools for studying their functions. Given the selectivity of any inhibitor for a channel is relative, more than one inhibitor of different affinities may be used to help identify the channel in a biological preparation. Here, we describe a family of small proteins in centipede venoms that inhibit the pore (hKir6.2) of a human ATP-sensitive K+ channel (hKATP). While the traditional peptide-sequencing service gradually vanishes from academic institutions, we tried to identify the sequences of inhibitory proteins purified from venoms by searching the sequences of the corresponding transcriptomes, a search guided by the key features of a known hKir6.2 inhibitor (SpTx1). The candidate sequences were cross-checked against the masses of purified proteins, and validated by testing the activity of recombinant proteins against hKir6.2. The four identified proteins (SsdTx1-3 and SsTx) inhibit hKATP channels with a Kd of <300 nM, compared to 15 nM for SpTx1. SsTx has previously been discovered to block human voltage-gated KCNQ K+ channels with a 2.5 µM Kd. Given that SsTx inhibits hKir6.2 with >10-fold lower Kd than it inhibits hKCNQ, SsTx may not be suitable for probing KCNQ channels in a biological preparation that also contains more-SsTx-sensitive KATP channels.

Parallel Evolution of Complex Centipede Venoms Revealed by Comparative Proteotranscriptomic Analyses.

Centipedes are among the most ancient groups of venomous predatory arthropods. Extant species belong to five orders, but our understanding of the composition and evolution of centipede venoms is based almost exclusively on one order, Scolopendromorpha. To gain a broader and less biased understanding we performed a comparative proteotranscriptomic analysis of centipede venoms from all five orders, including the first venom profiles for the orders Lithobiomorpha, Craterostigmomorpha, and Geophilomorpha. Our results reveal an astonishing structural diversity of venom components, with 93 phylogenetically distinct protein and peptide families. Proteomically-annotated gene trees of these putative toxin families show that centipede venom composition is highly dynamic across macroevolutionary timescales, with numerous gene duplications as well as functional recruitments and losses of toxin gene families. Strikingly, not a single family is found in the venoms of representatives of all five orders, with 67 families being unique for single orders. Ancestral state reconstructions reveal that centipede venom originated as a simple cocktail comprising just four toxin families, with very little compositional evolution happening during the approximately 50 My before the living orders had diverged. Venom complexity then increased in parallel within the orders, with scolopendromorphs evolving particularly complex venoms. Our results show that even venoms composed of toxins evolving under the strong constraint of negative selection can have striking evolutionary plasticity on the compositional level. We show that the functional recruitments and losses of toxin families that shape centipede venom arsenals are not concentrated early in their evolutionary history, but happen frequently throughout.

Sex-based venom variation in the eastern bark centipede (Hemiscolopendra marginata).

Sexually dimorphic traits are widespread across metazoans and are often the result of sex-specific inheritance or sex-based differences in gene expression. Intersexual differences have even been observed in invertebrate venoms, although the identification of these differences has been limited to the more well-studied groups, such as scorpions and spiders, where sex-based differences in morphology and behavior are apparent. Recent studies on centipede venom have identified evidence of intraspecific variation, but intersexual differences have not been reported. To investigate the potential for sex-based differences in centipede venom composition, we performed reversed-phase high performance liquid chromatography (RP-HPLC) analyses on five male and 15 female eastern bark centipedes (Hemiscolopendra marginata) from the Apalachicola National Forest in northern Florida. After detecting a significant sex-based difference in H. marginata venom composition, we completed a high-throughput venom-gland transcriptomic and venom proteomic analysis of one male and one female to determine the genetic basis for differences in venom composition. We identified 47 proteomically confirmed toxins and 717 nontoxin transcripts in H. marginata venom-glands. Of these proteomically confirmed toxins, the most abundantly expressed in the male venom included ion channel-modulating toxins and toxins so divergent from any characterized homologs that they could not be given a functional classification, whereas the most abundantly expressed in the female venom were γ-glutamyl transferases and CAPs (cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 proteins). These differences were then confirmed by performing replicate LC-MS/MS analyses on the venom from an additional three male and three female H. marginata. Our RP-HPLC and high-throughput transcriptomic and proteomic approach resulted in not only an in-depth characterization of H. marginata venom, but represents the first example of sex-based variation in centipede venoms.

Toxins from scratch? Diverse, multimodal gene origins in the predatory robber fly Dasypogon diadema indicate a dynamic venom evolution in dipteran insects.

BACKGROUND: Venoms and the toxins they contain represent molecular adaptations that have evolved on numerous occasions throughout the animal kingdom. However, the processes that shape venom protein evolution are poorly understood because of the scarcity of whole-genome data available for comparative analyses of venomous species. RESULTS: We performed a broad comparative toxicogenomic analysis to gain insight into the genomic mechanisms of venom evolution in robber flies (Asilidae). We first sequenced a high-quality draft genome of the hymenopteran hunting robber fly Dasypogon diadema, analysed its venom by a combined proteotranscriptomic approach, and compared our results with recently described robber fly venoms to assess the general composition and major components of asilid venom. We then applied a comparative genomics approach, based on 1 additional asilid genome, 10 high-quality dipteran genomes, and 2 lepidopteran outgroup genomes, to reveal the evolutionary mechanisms and origins of identified venom proteins in robber flies. CONCLUSIONS: While homologues were identified for 15 of 30 predominant venom protein in the non-asilid genomes, the remaining 15 highly expressed venom proteins appear to be unique to robber flies. Our results reveal that the venom of D. diadema likely evolves in a multimodal fashion comprising (i) neofunctionalization after gene duplication, (ii) expression-dependent co-option of proteins, and (iii) asilid lineage-specific orphan genes with enigmatic origin. The role of such orphan genes is currently being disputed in evolutionary genomics but has not been discussed in the context of toxin evolution. Our results display an unexpected dynamic venom evolution in asilid insects, which contrasts the findings of the only other insect toxicogenomic evolutionary analysis, in parasitoid wasps (Hymenoptera), where toxin evolution is dominated by single gene co-option. These findings underpin the significance of further genomic studies to cover more neglected lineages of venomous taxa and to understand the importance of orphan genes as possible drivers for venom evolution.

Convergent recruitment of adamalysin-like metalloproteases in the venom of the red bark centipede (Scolopocryptops sexspinosus).

Many venom proteins have presumably been convergently recruited by taxa from diverse venomous lineages. These toxic proteins have characteristics that allow them to remain stable in solution and have a high propensity for toxic effects on prey and/or potential predators. Despite this well-established convergent toxin recruitment, some toxins seem to be lineage specific. To further investigate the toxic proteins found throughout venomous lineages, venom proteomics and venom-gland transcriptomics were performed on two individual red bark centipedes (Scolopocryptops sexspinosus). Combining the protein phenotype with the transcript genotype resulted in the first in-depth venom characterization of S. sexspinosus, including 72 venom components that were identified in both the transcriptome and proteome and 1468 nontoxin transcripts identified in the transcriptome. Ten different toxin families were represented in the venom and venom gland with the majority of the toxins belonging to metalloproteases, CAPS (cysteine-rich secretory protein, antigen 5, and pathogenesis-related 1 proteins), and ß-pore-forming toxins. Nine of these toxin families shared a similar proteomic structure to venom proteins previously identified from other centipedes. However, the most highly expressed toxin family, the adamalysin-like metalloproteases, has until now only been observed in the venom of snakes. We confirmed adamalysin-like metalloprotease activity by means of in vivo functional assays. The recruitment of an adamalysin-like metalloprotease into centipede venom represents a striking case of convergent evolution.

Venom-gland transcriptomics and venom proteomics of the giant Florida blue centipede, Scolopendra viridis.

The limited number of centipede venom characterizations have revealed a rich diversity of toxins, and recent work has suggested centipede toxins may be more rapidly diversifying than previously considered. Additionally, many identified challenges in venomics research, including assembly and annotation methods, toxin quantification, and the ability to provide biological or technical replicates, have yet to be addressed in centipede venom characterizations. We performed high-throughput, quantifiable transcriptomic and proteomic methods on two individual Scolopendra viridis centipedes from North Florida. We identified 39 toxins that were proteomically confirmed, and 481 nontoxins that were expressed in the venom gland of S. viridis. The most abundant toxins expressed in the venom of S. viridis belonged to calcium and potassium ion-channel toxins, venom allergens, metalloproteases, and ß-pore forming toxins. We compared our results to the previously characterized S. viridis from Morelos, Mexico, and found only five proteomically confirmed toxins in common to both localities, suggesting either extreme toxin divergence within S. viridis, or that these populations may represent entirely different species. By using multiple assembly and annotation methods, we generated a comprehensive and quantitative reference transcriptome and proteome of a Scolopendromorpha centipede species, while overcoming some of the challenges present in venomics research.

True Lies: Using Proteomics to Assess the Accuracy of Transcriptome-Based Venomics in Centipedes Uncovers False Positives and Reveals Startling Intraspecific Variation in Scolopendra Subspinipes.

Centipede venoms have emerged as a rich source of novel bioactive compounds. However, most centipede species are commonly considered too small for venom extraction and transcriptomics is likely to be an attractive way of probing the molecular diversity of these venoms. Examining the venom composition of Scolopendra subspinipes, we test the accuracy of this approach. We compared the proteomically determined venom profile with four common toxin transcriptomic toxin annotation approaches: BLAST search against toxins in UniProt, lineage-specific toxins, or species-specific toxins and comparative expression analyses of venom and non-venom producing tissues. This demonstrated that even toxin annotation based on lineage-specific homology searches is prone to substantial errors compared to a proteomic approach. However, combined comparative transcriptomics and phylogenetic analysis of putative toxin families substantially improves annotation accuracy. Furthermore, comparison of the venom composition of S. subspinipes with the closely related S. subspinipes mutilans revealed a surprising lack of overlap. This first insight into the intraspecific venom variability of centipedes contrasts the sequence conservation expected from previous findings that centipede toxins evolve under strong negative selection. Our results highlight the importance of proteomic data in studies of even comparably well-characterized venoms and warrants caution when sourcing venom from centipedes of unknown origin.

Giant fish-killing water bug reveals ancient and dynamic venom evolution in Heteroptera.

True Bugs (Insecta: Heteroptera) produce venom or saliva with diverse bioactivities depending on their feeding strategies. However, little is known about the molecular evolution of the venom toxins underlying these biological activities. We examined venom of the giant fish-killing water bug Lethocerus distinctifemur (Insecta: Belostomatidae) using infrared spectroscopy, transcriptomics, and proteomics. We report 132 venom proteins including putative enzymes, cytolytic toxins, and antimicrobial peptides. Over 73% (96 proteins) showed homology to venom proteins from assassin bugs (Reduviidae), including 21% (28 proteins from seven families) not known from other sources. These data suggest that numerous protein families were recruited into venom and diversified rapidly following the switch from phytophagy to predation by ancestral heteropterans, and then were retained over > 200 my of evolution. In contrast, trophic switches to blood-feeding (e.g. in Triatominae and Cimicidae) or reversions to plant-feeding (e.g., in Pentatomomorpha) were accompanied by rapid changes in the composition of venom/saliva, including the loss of many protein families.

The assassin bug Pristhesancus plagipennis produces two distinct venoms in separate gland lumens.

The assassin bug venom system plays diverse roles in prey capture, defence and extra-oral digestion, but it is poorly characterised, partly due to its anatomical complexity. Here we demonstrate that this complexity results from numerous adaptations that enable assassin bugs to modulate the composition of their venom in a context-dependent manner. Gland reconstructions from multimodal imaging reveal three distinct venom gland lumens: the anterior main gland (AMG); posterior main gland (PMG); and accessory gland (AG). Transcriptomic and proteomic experiments demonstrate that the AMG and PMG produce and accumulate distinct sets of venom proteins and peptides. PMG venom, which can be elicited by electrostimulation, potently paralyses and kills prey insects. In contrast, AMG venom elicited by harassment does not paralyse prey insects, suggesting a defensive role. Our data suggest that assassin bugs produce offensive and defensive venoms in anatomically distinct glands, an evolutionary adaptation that, to our knowledge, has not been described for any other venomous animal.