Venom regeneration in the centipede Scolopendra polymorpha: evidence for asynchronous venom component synthesis

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.

The arylsulfatase- and phospholipase-rich venom of the plutoniumid centipede Theatops posticus

Research on centipede venoms has led to the discovery of a diverse array of novel proteins and peptides, including those with homology to previously discovered toxin families (e.g., phospholipase A2s and pM12a metalloproteases) and novel toxin families not previously detected in venoms (e.g., β-pore forming toxins and scoloptoxins). Most of this research has focused on centipedes in the order Scolopendromorpha, particularly those in the families Scolopendridae, Cryptopidae, and Scolopocryptopidae. To generate the first high-throughput venom characterization for a centipede in the scolopendromorph family Plutoniumidae, we performed venom-gland transcriptomics and venom proteomics on two Theatops posticus. We identified a total of 64 venom toxins, 60 of which were detected in both the venom-gland transcriptome and venom proteome and four of which were only detected transcriptomically. We detected a single highly abundant arylsulfatase B (ARSB) toxin, the first ARSB toxin identified from centipede venoms. As ARSBs have been detected in other venomous species (e.g., scorpions), ARSBs in T. posticus highlights a new case of convergent evolution across venoms. Theatops posticus venom also contained a much higher abundance and diversity of phospholipase A2 toxins compared to other characterized centipede venoms. Conversely, we detected other common centipedes toxins, such as CAPs and scoloptoxins, at relatively low abundances and diversities. Our observation of a diverse set of toxins from T. posticus venom, including those from novel toxin families, emphasizes the importance of studying unexplored centipede taxonomic groups and the continued potential of centipede venoms for novel toxin discovery and unraveling the molecular mechanisms underlying trait evolution.

Contrasting patterns of venom regeneration in a centipede (Scolopendra viridis) and a scorpion (Centruroides hentzi)

As biochemical traits with clear fitness consequences, venoms serve a critical ecological role for the animals that produce them. Understanding how venoms are maintained and regenerated after use will, therefore, provide valuable insight into the ecology of venomous animals. Furthermore, most studies on venomous organisms often require removing animals from the wild and waiting extended periods of time between venom extractions. Uncovering the patterns of venom regeneration across different species will likely lead to the development of more efficient venom extraction protocols, reducing both experimental time and the number of animals required. Using reversed-phase high-performance liquid chromatography, we identified asynchronous regeneration of venom protein component abundances in the centipede Scolopendra viridis but found no evidence for asynchronous venom regeneration in the scorpion Centruroides hentzi. We also observed high levels of intraspecific venom variation in C. hentzi, emphasizing the importance of testing for intraspecific venom variation in studies evaluating the synchronicity of venom regeneration. Although the regeneration of relative venom protein component abundances is an asynchronous process in S. viridis, we provide evidence that the presence-absence of major venom components is not an asynchronous process and suggest that studies relying on just the presence/absence of individual proteins (e.g. bioprospecting, drug discovery) could use catch-and-release methods of venom extraction to reduce the number of animals removed from the wild.

Venom Collection from Spiders and Snakes: Voluntary and Involuntary Extractions ("Milking") and Venom Gland Extractions

Venom collection (often called "milking") provides the toxic secretions essential for studying animal venoms and/or generating venom products. Methods of venom collection vary widely, falling into three broad categories: voluntary venom extraction (inducing the animal to willingly release its venom), involuntary venom extraction (glandular massage, electrical stimulation, or administration of induction chemicals to promote venom expulsion), and venom gland extraction (surgical aspiration or trituration of homogenized gland tissue). Choice of method requires consideration of animal species, animal welfare, human safety (avoiding envenomation), venom yield and composition desired, and level of toxin purity required. Here, we summarize the materials and methods used to obtain venom by each of these approaches from spiders and snakes.

Kinetic analysis of effects of temperature and time on the regulation of venom expression in Bungarus multicinctus

Venom gland is a highly efficient venom production system to maintain their predatory arsenal. Venom toxins mRNA has been shown to increase abruptly in snake after venom expenditure, while the dynamics of venom accumulation during synthesis are poorly understood. Here, PacBio long-read sequencing, Illumina RNA sequencing (RNA-seq), and label-free proteome quantification were used to investigate the composition landscape and time- and temperature-dependent dynamics changes of the Bungarus multicinctus venom gland system. Transcriptome data (19.5223 Gb) from six adult B. multicinctus tissues were sequenced using PacBio RS II to generate a reference assembly, and average 7.28 Gb of clean RNA-seq data was obtained from venom glands by Illumina sequencing. Differentially expressed genes (DEGs) mainly were protein processing rather than venom toxins. Post-translational modifications provided the evidence of the significantly different proportions of toxins in the venom proteome with the changing of replenishment time and temperature, but constant of venom toxins mRNA in the venom gland transcriptome. Dynamic of toxins and genes involved in venom gland contraction suggesting the formation of the mature venom gland system would take at least 9 days. In addition, 59 toxin processing genes were identified, peptidylprolyl isomerase B of which underwent positive selection in Toxicofera. These results provide a reference for determining the extraction time of venom, production of polyclonal and monoclonal antibody for precise treatment plans of venomous snakebites, and construction of an in vitro synthesis system for snake venom protein.

The Diversity of Venom: The Importance of Behavior and Venom System Morphology in Understanding Its Ecology and Evolution

Venoms are one of the most convergent of animal traits known, and encompass a much greater taxonomic and functional diversity than is commonly appreciated. This knowledge gap limits the potential of venom as a model trait in evolutionary biology. Here, we summarize the taxonomic and functional diversity of animal venoms and relate this to what is known about venom system morphology, venom modulation, and venom pharmacology, with the aim of drawing attention to the importance of these largely neglected aspects of venom research. We find that animals have evolved venoms at least 101 independent times and that venoms play at least 11 distinct ecological roles in addition to predation, defense, and feeding. Comparisons of different venom systems suggest that morphology strongly influences how venoms achieve these functions, and hence is an important consideration for understanding the molecular evolution of venoms and their toxins. Our findings also highlight the need for more holistic studies of venom systems and the toxins they contain. Greater knowledge of behavior, morphology, and ecologically relevant toxin pharmacology will improve our understanding of the evolution of venoms and their toxins, and likely facilitate exploration of their potential as sources of molecular tools and therapeutic and agrochemical lead compounds.

Venom content and toxicity regeneration after venom gland depletion by electrostimulation in the scorpion Centruroides limpidus

The scorpion venom is a cocktail of many components. Its composition can exhibit a level of plasticity in response to different behavioral and environmental factors, leading to intraspecific variation. The toxicity and specificity of scorpion venoms appear to be taxon-dependent, due to a co-evolutionary interaction with prey and predators, which shaped the composition at the molecular level. The venom regeneration by the venom glands is an asynchronous process, in which particular components are expressed at different stages and at different rates. According to this, it can be reasonably assumed that the regeneration of toxicity in the venom is also asynchronous. In this work, we studied the toxicity regeneration dynamics by the scorpion Centruroides limpidus after full venom depletion by electrical stimulation. For this, we evaluated the toxicity of venom samples extracted at different days post depletion, against insects (crickets) and mammals (humans, by assessing the venom activity on the human voltage-dependent Na⁺ channel Nav1.6). The regeneration of toxicity against humans lagged behind that against crickets (13 vs 10 days, respectively). Thirteen days after depletion the venom seems to be replenished. Our results show asynchrony in the regeneration of species-specific toxic activity in the venom of Centruroides limpidus. The understanding of the venom regeneration kinetics for the different scorpion species will help to design venom extraction protocols that could maximize the yield and quality of the collected venoms.

Studying Smaller and Neglected Organisms in Modern Evolutionary Venomics Implementing RNASeq (Transcriptomics)-A Critical Guide

Venoms are evolutionary key adaptations that species employ for defense, predation or competition. However, the processes and forces that drive the evolution of venoms and their toxin components remain in many aspects understudied. In particular, the venoms of many smaller, neglected (mostly invertebrate) organisms are not characterized in detail, especially with modern methods. For the majority of these taxa, even their biology is only vaguely known. Modern evolutionary venomics addresses the question of how venoms evolve by applying a plethora of -omics methods. These recently became so sensitive and enhanced that smaller, neglected organisms are now more easily accessible to comparatively study their venoms. More knowledge about these taxa is essential to better understand venom evolution in general. The methodological core pillars of integrative evolutionary venomics are genomics, transcriptomics and proteomics, which are complemented by functional morphology and the field of protein synthesis and activity tests. This manuscript focuses on transcriptomics (or RNASeq) as one toolbox to describe venom evolution in smaller, neglected taxa. It provides a hands-on guide that discusses a generalized RNASeq workflow, which can be adapted, accordingly, to respective projects. For neglected and small taxa, generalized recommendations are difficult to give and conclusions need to be made individually from case to case. In the context of evolutionary venomics, this overview highlights critical points, but also promises of RNASeq analyses. Methodologically, these concern the impact of read processing, possible improvements by perfoming multiple and merged assemblies, and adequate quantification of expressed transcripts. Readers are guided to reappraise their hypotheses on venom evolution in smaller organisms and how robustly these are testable with the current transcriptomics toolbox. The complementary approach that combines particular proteomics but also genomics with transcriptomics is discussed as well. As recently shown, comparative proteomics is, for example, most important in preventing false positive identifications of possible toxin transcripts. Finally, future directions in transcriptomics, such as applying 3rd generation sequencing strategies to overcome difficulties by short read assemblies, are briefly addressed.

Ultrastructure of the salivary glands in Lithobius forficatus (Myriapoda, Chilopoda, Lithobiidae) according to seasonal and circadian rhythms

The salivary glands (mandibular epidermal glands) of adult males and females of Lithobius forficatus (Myriapoda, Chilopoda) were isolated during spring, summer and autumn. In addition, the organs were isolated at different times of the day - at about 12:00 (noon) and about 00:00 (midnight). The ultrastructure of these organs depending on seasonal and circadian rhythms was analyzed using transmission and scanning electron microscopy and histochemical methods. The paired salivary glands of L. forficatus are situated in the vicinity of the foregut and they are formed by numerous acini that are surrounded by the fat body, hemocytes and tracheolae. The salivary glands are composed of a terminal acinar component and a system of tubular ducts that are lined with a cuticle. The glandular part is composed of secretory epithelial cells that are at various stages of their secretory activity. The saliva that is produced by the secretory cells of the acini is secreted into the salivary ducts, which are lined with a simple epithelium that is based on the non-cellular basal lamina. The ultrastructural variations suggest that salivary glands function differently depending on seasonal rhythms and prepare the animal for overwintering. Therefore, the salivary glands of the centipedes that were analyzed participate in the accumulation of proteins, lipids and polysaccharides during the spring, summer and autumn. Subtle differences in the ultrastructure of the secretory cells of the salivary glands during the circadian cycle must be related to the physiological reactions of the organism. The salivary ducts showed no differences in the specimens that were analyzed during the day/night cycle or during the seasonal cycle.

Advances in venomics

The term “venomics” was coined to describe the global study of venom and venom glands, targeting comprehensive characterization of the whole toxin profile of a venomous animal by means of proteomics, transcriptomics, genomics and bioinformatics studies.