Varying Modes of Selection Among Toxin Families in the Venoms of the Giant Desert Hairy Scorpions (Hadrurus).

Venoms are primarily believed to evolve under strong diversifying selection resulting from persistent coevolution between predator and prey. Recent research has challenged this hypothesis, proposing that venoms from younger venomous lineages (e.g., snakes and cone snails) are governed predominantly by diversifying selection, while venoms from older venomous lineages (e.g., centipedes, scorpions, and spiders) are under stronger purifying selection. However, most research in older lineages has tested selection at more diverse phylogenetic scales. Although these tests are important for evaluating broad macroevolutionary trends underlying venom evolution, they are less equipped to detect species-level evolutionary trends, which likely have large impacts on venom variation seen at more diverse phylogenetic scales. To test for selection among closely related species from an older venomous lineage, we generated high-throughput venom-gland transcriptomes and venom proteomes for four populations of Giant Desert Hairy Scorpions (Hadrurus), including three Hadrurus arizonensis populations and one Hadrurus spadix population. We detected significant episodic and pervasive diversifying selection across a highly abundant toxin family that likely has a major role in venom function ([Formula: see text]KTxs), providing a contrast to the stronger purifying selection identified from other studies on scorpion venoms. Conversely, we detected weak episodic diversifying and/or stronger purifying selection in four toxin families (non-disulfide bridged peptides, phospholipase A2s, scorpine-like antimicrobial peptides, and serine proteases), most of which were less abundant and likely have ancillary functional roles. Finally, although we detected several major toxin families at disproportionate transcriptomic and/or proteomic abundances, we did not identify significant sex-based variation in Hadrurus venoms.

Identification and Characterization of Novel Proteins from Arizona Bark Scorpion Venom That Inhibit Nav1.8, a Voltage-Gated Sodium Channel Regulator of Pain Signaling.

The voltage-gated sodium channel Nav1.8 is linked to neuropathic and inflammatory pain, highlighting the potential to serve as a drug target. However, the biophysical mechanisms that regulate Nav1.8 activation and inactivation gating are not completely understood. Progress has been hindered by a lack of biochemical tools for examining Nav1.8 gating mechanisms. Arizona bark scorpion (Centruroides sculpturatus) venom proteins inhibit Nav1.8 and block pain in grasshopper mice (Onychomys torridus). These proteins provide tools for examining Nav1.8 structure-activity relationships. To identify proteins that inhibit Nav1.8 activity, venom samples were fractioned using liquid chromatography (reversed-phase and ion exchange). A recombinant Nav1.8 clone expressed in ND7/23 cells was used to identify subfractions that inhibited Nav1.8 Na+ current. Mass-spectrometry-based bottom-up proteomic analyses identified unique peptides from inhibitory subfractions. A search of the peptides against the AZ bark scorpion venom gland transcriptome revealed four novel proteins between 40 and 60% conserved with venom proteins from scorpions in four genera (Centruroides, Parabuthus, Androctonus, and Tityus). Ranging from 63 to 82 amino acids, each primary structure includes eight cysteines and a "CXCE" motif, where X = an aromatic residue (tryptophan, tyrosine, or phenylalanine). Electrophysiology data demonstrated that the inhibitory effects of bioactive subfractions can be removed by hyperpolarizing the channels, suggesting that proteins may function as gating modifiers as opposed to pore blockers.

The Chemosensory Repertoire of the Eastern Diamondback Rattlesnake (Crotalus adamanteus) Reveals Complementary Genetics of Olfactory and Vomeronasal-Type Receptors.

Pitviper sensory perception incorporates diverse stimuli through the integration of trichromatic color vision, bifocal heat-sensing, and dual-system chemoperception. Chemoperception, or olfaction, is mediated by chemoreceptors in the olfactory bulb and the vomeronasal organ, but the true genomic complexity of the gene families and their relative contributions is unknown. A full genomic accounting of pitviper chemoperception directly complements our current understanding of their venoms by generating a more complete polyphenic representation of their predatory arsenal. To characterize the genetic repertoire of pitviper chemoperception, we analyzed a full-genome assembly for Crotalus adamanteus, the eastern diamondback rattlesnake. We identified hundreds of genes encoding both olfactory receptors (ORs; 362 full-length genes) and type-2 vomeronasal receptors (V2Rs; 430 full-length genes). Many chemoreceptor genes are organized into large tandem repeat arrays. Comparative analysis of V2R orthologs across squamates demonstrates how gene array expansion and contraction underlies the evolution of the chemoreceptor repertoire, which likely reflects shifts in life history traits. Chromosomal assignments of chemosensory genes identified sex chromosome specific chemoreceptor genes, providing gene candidates underlying observed sex-specific chemosensory-based behaviors. We detected widespread episodic evolution in the extracellular, ligand-binding domains of both ORs and V2Rs, suggesting the diversification of chemoreceptors is driven by transient periods of positive selection. We provide a robust genetic framework for studying pitviper chemosensory ecology and evolution.

Venom Peptides with Dual Modulatory Activity on the Voltage-Gated Sodium Channel NaV1.1 Provide Novel Leads for Development of Antiepileptic Drugs.

Voltage-gated sodium (NaV) channels play a fundamental role in normal neurological function, especially via the initiation and propagation of action potentials. The NaV1.1 subtype is found in inhibitory interneurons of the brain and it is essential for maintaining a balance between excitation and inhibition in neuronal networks. Heterozygous loss-of-function mutations of SCN1A, the gene encoding NaV1.1, underlie Dravet syndrome (DS), a severe pediatric epilepsy. We recently demonstrated that selective inhibition of NaV1.1 inactivation prevents seizures and premature death in a mouse model of DS. Thus, selective modulators of NaV1.1 might be useful therapeutics for treatment of DS as they target the underlying molecular deficit. Numerous scorpion-venom peptides have been shown to modulate the activity of NaV channels, but little is known about their activity at NaV1.1. Here we report the isolation, sequence, three-dimensional structure, recombinant production, and functional characterization of two peptidic modulators of NaV1.1 from venom of the buthid scorpion Hottentotta jayakari. These peptides, Hj1a and Hj2a, are potent agonists of NaV1.1 (EC50 of 17 and 32 nM, respectively), and they present dual α/ß activity by modifying both the activation and inactivation properties of the channel. NMR studies of rHj1a indicate that it adopts a cystine-stabilized αß fold similar to known scorpion toxins. Although Hj1a and Hj2a have only limited selectivity for NaV1.1, their unusual dual mode of action provides an alternative approach to the development of selective NaV1.1 modulators for the treatment of DS.

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.

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.

Tipping the Scales: The Migration-Selection Balance Leans toward Selection in Snake Venoms.

The migration-selection interaction is the strongest determinant of whether a beneficial allele increases in frequency within a population. Results of empirical studies examining the role of gene flow in an adaptive context, however, have largely been equivocal, with examples of opposing outcomes being repeatedly documented (e.g., local adaptation with high levels of gene flow vs. gene swamping). We compared neutral genomic and venom expression divergence for three sympatric pit vipers with differing ecologies to determine if and how migration-selection disequilibria result in local adaptation. We specifically tested whether neutral differentiation predicted phenotypic differentiation within an isolation-by-distance framework. The decoupling of neutral and phenotypic differentiation would indicate selection led to adaptive divergence irrespective of migration, whereas a significant relationship between neutral and venom expression differentiation would provide evidence in favor of the constraining force of gene flow. Neutral differentiation and geographic distance predicted phenotypic differentiation only in the generalist species, indicating that selection was the predominant mechanism in the migration-selection balance underlying adaptive venom evolution in both specialists. Dispersal is thought to be a stronger influence on genetic differentiation than specialization, but our results suggest the opposite. Greater specialization may lead to greater diversification rates, and extreme spatial and temporal variation in selective pressures can favor generalist phenotypes evolving under strong stabilizing selection. Our results are consistent with these expectations, suggesting that the equivocal findings of studies examining the role of gene flow in an adaptive context may be explained by ecological specialization theory.

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.

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.

Female-biased population divergence in the venom of the Hentz striped scorpion (Centruroides hentzi).

Sex-biased genes are expressed at higher levels in one sex and contribute to phenotypic differences between males and females, as well as overall phenotypic variation within and among populations. Venom has evolved primarily for predation and defense, making venom expression a highly variable phenotype as a result of local adaptation. Several scorpion species have shown both intraspecific and intersexual venom variation, and males have been observed using venom in courtship and mating, suggesting the existence of venom-specific, sex-biased genes that may contribute to population divergence. We used reversed-phase high-performance liquid chromatography (RP-HPLC), Agilent protein bioanalyzer chips, nano-liquid chromatography mass spectrometry (nLC/MS/MS), and median lethal dose (LD50) assays in fruit flies (Drosophila melanogaster) and banded crickets (Gryllodes sigillatus) to investigate proteomic and functional venom variation within and among three Florida populations of the Hentz striped scorpion (Centruroides hentzi). We found significant venom variation among populations, with females, not males, being responsible for this divergence. We also found significant variation in venom expression within populations, with males contributing more to within population variation than females. Our results provide evidence that male and female scorpions experience different natural and sexual selective pressures that have led to the expression of sex-biased venom genes and that these genes may be consequential in population divergence.

A global accounting of medically significant scorpions: Epidemiology, major toxins, and comparative resources in harmless counterparts.

Scorpions are an ancient and diverse venomous lineage, with over 2200 currently recognized species. Only a small fraction of scorpion species are considered harmful to humans, but the often life-threatening symptoms caused by a single sting are significant enough to recognize scorpionism as a global health problem. The continued discovery and classification of new species has led to a steady increase in the number of both harmful and harmless scorpion species. The purpose of this review is to update the global record of medically significant scorpion species, assigning each to a recognized sting class based on reported symptoms, and provide the major toxin classes identified in their venoms. We also aim to shed light on the harmless species that, although not a threat to human health, should still be considered medically relevant for their potential in therapeutic development. Included in our review is discussion of the many contributing factors that may cause error in epidemiological estimations and in the determination of medically significant scorpion species, and we provide suggestions for future scorpion research that will aid in overcoming these errors.

Venom-gland transcriptomics and venom proteomics of the Hentz striped scorpion (Centruroides hentzi; Buthidae) reveal high toxin diversity in a harmless member of a lethal family.

Of the 14 extant scorpion families, Buthidae has the most thoroughly characterized venoms. Most of this characterization, however, has been limited to species with medically significant stings, including members of the Centruroides genus, which have caused human deaths (e.g., Centruroides sculpturatus). To understand the origin and evolution of highly toxic venoms, we should also characterize the more harmless venoms of close relatives. We used Illumina sequencing to separately characterize the venom-gland transcriptomes of a male and female Hentz striped scorpion (Centruroides hentzi) and performed independent quantitative mass-spectrometry analysis of the venom from each individual, providing the first full venom characterization of a Centruroides species that poses no serious threat to humans. We identified 59 venom proteins that were proteomically confirmed, 63 additional transcripts that were identified on the basis of homology to known toxins, and 355 nontoxins expressed in the venom-glands. The most abundant toxins belonged to the Na+ and K+-channel toxin classes. Antimicrobial peptides and peptidases were also identified, along with a large group of venom proteins that could not be classified based on homology, suggesting C. hentzi is a source of previously untapped toxin diversity.

Quantity, Not Quality: Rapid Adaptation in a Polygenic Trait Proceeded Exclusively through Expression Differentiation.

A trait's genomic architecture can affect the rate and mechanism of adaptation, and although many ecologically-important traits are polygenic, most studies connecting genotype, phenotype, and fitness in natural populations have focused on traits with relatively simple genetic bases. To understand the genetic basis of polygenic adaptation, we must integrate genomics, phenotypic data, ecology, and fitness effects for a genetically tractable, polygenic trait; snake venoms provide such a system for studying polygenic adaptation because of their genetic tractability and vital ecological role in feeding and defense. We used a venom transcriptome-proteome map, quantitative proteomics, genomics, and fitness assays in sympatric prey to construct a genotype-phenotype-fitness map for the venoms of an island-mainland pair of rattlesnake populations. Reciprocal fitness experiments demonstrated that each population was locally adapted to sympatric prey. We identified significant expression differentiation with little to no coding-sequence variation across populations, demonstrating that expression differentiation was exclusively the genetic basis of polygenic adaptation. Previous research on the genetics of adaptation, however, has largely been biased toward investigating protein-coding regions because of the complexity of gene regulation. Our results showed that biases at the molecular level can be in the opposite direction, highlighting the need for more systematic comparisons of different molecular mechanisms underlying rapid, adaptive evolution in polygenic traits.

The genetics of venom ontogeny in the eastern diamondback rattlesnake (Crotalus adamanteus).

The same selective forces that give rise to rapid inter- and intraspecific divergence in snake venoms can also favor differences in venoms across life-history stages. Ontogenetic changes in venom composition are well known and widespread in snakes but have not been investigated to the level of unambiguously identifying the specific loci involved. The eastern diamondback rattlesnake was previously shown to undergo an ontogenetic shift in venom composition at sexual maturity, and this shift accounted for more venom variation than geography. To characterize the genetics underlying the ontogenetic venom compositional change in C. adamanteus, we sequenced adult/juvenile pairs of venom-gland transcriptomes from five populations previously shown to have different adult venom compositions. We identified a total of 59 putative toxin transcripts for C. adamanteus, and 12 of these were involved in the ontogenetic change. Three toxins were downregulated, and nine were upregulated in adults relative to juveniles. Adults and juveniles expressed similar total levels of snake-venom metalloproteinases but differed substantially in their featured paralogs, and adults expressed higher levels of Bradykinin-potentiating and C-type natriuretic peptides, nerve growth factor, and specific paralogs of phospholipases A2 and snake venom serine proteinases. Juvenile venom was more toxic to mice, indicating that the expression differences resulted in a phenotypically, and therefore potentially ecologically, significant difference in venom function. We also showed that adult and juvenile venom-gland transcriptomes for a species with known ontogenetic venom variation were equally effective at individually providing a full characterization of the venom genes of a species but that any particular individual was likely to lack several toxins in their transcriptome. A full characterization of a species' venom-gene complement therefore requires sequencing more than one individual, although the ages of the individuals are unimportant.

Venom-gland transcriptomics and venom proteomics of the black-back scorpion (Hadrurus spadix) reveal detectability challenges and an unexplored realm of animal toxin diversity.

The order Scorpiones is one of the most ancient and diverse lineages of venomous animals, having originated approximately 430 million years ago and diversified into 14 extant families. Although partial venom characterizations have been described for numerous scorpion species, we provided the first quantitative transcriptome/proteome comparison for a scorpion species using single-animal approaches. We sequenced the venom-gland transcriptomes of a male and female black-back scorpion (Hadrurus spadix) from the family Caraboctonidae using the Illumina sequencing platform and conducted independent quantitative mass-spectrometry analyses of their venoms. We identified 79 proteomically confirmed venom proteins, an additional 69 transcripts with homology to toxins from other species, and 596 nontoxin proteins expressed at high levels in the venom glands. The venom of H. spadix was rich in antimicrobial peptides, K+-channel toxins, and several classes of peptidases. However, the most diverse and one of the most abundant classes of putative toxins could not be assigned even a tentative functional role on the basis of homology, indicating that this venom contained a wealth of previously unexplored animal toxin diversity. We found good agreement between both transcriptomic and proteomic abundances across individuals, but transcriptomic and proteomic abundandances differed substantially within each individual. Small peptide toxins such as K+-channel toxins and antimicrobial peptides proved challenging to detect proteomically, at least in part due to the significant proteolytic processing involved in their maturation. In addition, we found a significant tendency for our proteomic approach to overestimate the abundances of large putative toxins and underestimate the abundances of smaller toxins.