RNA-seq and high-definition mass spectrometry reveal the complex and divergent venoms of two rear-fanged colubrid snakes

Dietary Breadth Predicts Toxin Expression Complexity in the Venoms of North American Gartersnakes

Selection on heritable phenotypic variation has played a prominent role in shaping the remarkable adaptations found across the Tree of Life. Complex ecological traits, such as snake venoms, are thought to be the products of selection because they directly link to fitness and survival. Snake venom increases the efficiency of prey capture and processing and is thus likely under intense selection. While many studies of snake venom have investigated the relationship between venom and diet, they have primarily focused on medically relevant front-fanged snakes. However, recent work has suggested that many non-front fanged snakes also rely on venom for subduing prey, despite having reduced toxicity toward humans. Here, we set out to uncover variation in toxin-producing genes, along with the ecological and evolutionary pressures impacting snake venom characteristics in the North American gartersnakes (Squamata: Natricidae: Thamnophis), a model group of non-front-fanged snakes. We annotated and analyzed Duvernoy's venom gland transcriptomes from 16 species representing all the major lineages within Thamnophis. We then generated measures of complexity of both toxins and dietary breadth. We found strong correlations between the complexity of toxin gene expression and phylogenetic diversity of diet, but no relationship between the complexity of the genetic makeup of the transcriptomes (allelic or sequence variation) and diet complexities. We also found phylogenetic signal associated with venom complexity, suggesting some influence of ancestry on venom characteristics. We suggest that, in non-front-fanged snakes, expression of toxins rather than sequence complexity is under strong selection by dietary diversity. These findings contradict similar studies from front-fanged snakes where increased transcriptomic complexity varies positively with dietary diversity, exposing a potential novel relationship between a complex phenotype-toxin expression-and its selective pressures-diet.

Venom proteomics and Duvernoy's venom gland histology of Pseudoboa neuwiedii (Neuwied's false boa; Dipsadidae, Pseudoboini)

The venom of Colombian specimens of the rear-fanged snake Pseudoboa neuwiedii contains proteolytic and phospholipase A2 (PLA2) activities, but is devoid of esterases. Mass spectrometric analysis of electrophoretic bands indicated that this venom contains C-type lectins (CTL), cysteine-rich secretory proteins (CRiSP), PLA2, snake venom metalloproteinases (SVMP), and snake venom matrix metalloproteinases (svMMP). In this investigation, we extended our characterization of P. neuwiedii by undertaking a shotgun proteomic analysis of the venom and comparing the results with a transcriptomic database for Brazilian P. neuwiedii; proteomic data previously obtained by in-gel digestion of electrophoretic bands coupled with mass spectrometry were also reanalyzed by comparing them with the transcriptomic results. The histology of the Duvernoy's venom gland was also examined. Histological analysis revealed a structural organization similar to that of other colubrids that consisted of a serous venom gland and a mucous supralabial gland. When the shotgun proteomic data were run against a general UniProt database for serpents, only metalloproteinases were identified (99% SVMPs, 1% snake endogenous matrix metalloproteinases-9 or seMMP-9). In contrast, when run against a transcriptomic database derived from the venom gland of Brazilian P. neuwiedii that contains predominantly SVMP, CRiSP, type IIE PLA2 (PLA2-IIE), CTL and seMMP-9, the main components identified were seMMP-9 (49%), SVMP (47%), CRiSP (3%) and minor components that included CTL and PLA2-IIE. These findings confirmed the previously reported general composition of P. neuwiedii venom, with metalloproteinases (SVMP and seMMP-9) being the major components, and refined the identification of certain components, e.g., type IIA PLA2 now identified as PLA2-IIE and the detection of seMMP-9 rather than svMMP. The data also indicate compositional similarity between Brazilian and Colombian P. neuwiedii venoms, and stress the need for specific databases for non-front-fanged colubroid snakes to allow accurate and more comprehensive identification of the venom components of these snakes.

Thamnophis sirtalis and their toxic relationship: Testing for intraspecific venom variation in Common Garter Snakes

Intraspecific phenotypic variation can be used as a window into the ecological differences among individuals of a species and lead to a better understanding of adaptive evolution. Adaptive traits, such as venom, that play an important ecological role for a species are useful models for understanding the sources of intraspecific variation. Intraspecific studies on front-fanged venomous snakes have offered deeper insights into the diverse mechanisms and adaptations that support the effectiveness of venom across species. Despite the extensive research on front-fanged venomous snakes, rear-fanged snakes, representing two-thirds of all snake species, have been largely overlooked. To test for sex and age-based intraspecific venom variation, we sequenced the messenger RNA from the Duvernoy's gland of 9 male and 10 female Common Garter Snakes, Thamnophis sirtalis, of different sizes from a single location. Our data represent the most venom gland transcriptomes of any venomous snake species from a single location and represent the first Duvernoy's venom gland transcriptomes for Thamnophis sirtalis. We found four toxin families dominate the Thamnophis sirtalis transcriptome: Snake Venom Metalloproteinases (SVMPs), Three-finger toxins (3FTxs), Cysteine-Rich Secretory Proteins (CRISPs), and C-type lectins (CTLs). Thamnophis sirtalis exhibits a unique balance in toxin expression, with approximately 30% each of neurotoxic (3FTx-dominated) and enzymatic (SVMP-dominated) components. No other published RFS Duvernoy's gland transcriptome displays this ratio, rather they are dominated by one or the other. Additionally, venom expression varies with sex and size, with differences in toxin gene expression between males and females as they grow. Our study provides new insights on venom composition in a RFS species and highlights the amount of intraspecific variation possible among individuals from a single population.

Highly conserved and extremely variable: The paradoxical pattern of toxin expression revealed by comparative venom-gland transcriptomics of Phalotris (Serpentes: Dipsadidae)

Although non-front fanged snakes account for almost two-thirds of snake diversity, most studies on venom composition and evolution focus exclusively on front-fanged species, which comprise most of the clinically relevant accidents. Comprehensive reports on venom composition of non-front fanged snakes are still scarce for several groups. In this study, we address such shortage of knowledge by providing new insights about the venom composition among species of Phalotris, a poorly studied Neotropical dipsadid genus. Phalotris are known for their specialized venom delivery system and toxic venoms, which can cause life-threatening accidents in humans. We evaluate the venom-gland transcriptome of Phalotris, comparing the following three South American species: P. reticulatus for the Araucaria Pine forests, P. lemniscatus for the Pampa grasslands, and P. mertensi for the Brazilian Cerrado. Our results indicate similar venom profiles, in which they share a high expression level of Kunitz-type inhibitors (KUNZ). On the other hand, comparative analyses revealed substantial differences in the expression levels of C-type lectins (CTL) and snake venom metalloproteinases (SVMP). The diverse set of SVMP and CTL isoforms shows signals of positive selection, and we also identified truncated forms of type III SVMPs, which resemble type II and type I SVMPs of viperids. Additionally, we identified a CNP precursor hosting a proline-rich region containing a BPP motif resembling those commonly detected in viperid venoms with hypotensive activity. Altogether, our results suggest an evolutionary history favoring high expression levels of few KUNZ isoforms in Phalotris venoms, contrasting with a highly diverse set of SVMP and CTL isoforms. Such diversity can be comparable with the venom variability observed in some viperids. Our findings highlight the extreme phenotypic diversity of non-front fanged snakes and the importance to allocate greater effort to study neglected groups of Colubroidea.

Assessing Target Specificity of the Small Molecule Inhibitor MARIMASTAT to Snake Venom Toxins: A Novel Application of Thermal Proteome Profiling

New treatments that circumvent the pitfalls of traditional antivenom therapies are critical to address the problem of snakebite globally. Numerous snake venom toxin inhibitors have shown promising cross-species neutralization of medically significant venom toxins in vivo and in vitro. The development of high-throughput approaches for the screening of such inhibitors could accelerate their identification, testing, and implementation and thus holds exciting potential for improving the treatments and outcomes of snakebite envenomation worldwide. Energetics-based proteomic approaches, including thermal proteome profiling and proteome integral solubility alteration (PISA) assays, represent "deep proteomics" methods for high throughput, proteome-wide identification of drug targets and ligands. In the following study, we apply thermal proteome profiling and PISA methods to characterize the interactions between venom toxin proteoforms in Crotalus atrox (Western Diamondback Rattlesnake) and the snake venom metalloprotease (SVMP) inhibitor marimastat. We investigate its venom proteome-wide effects and characterize its interactions with specific SVMP proteoforms, as well as its potential targeting of non-SVMP venom toxin families. We also compare the performance of PISA thermal window and soluble supernatant with insoluble precipitate using two inhibitor concentrations, providing the first demonstration of the utility of a sensitive high-throughput PISA-based approach to assess the direct targets of small molecule inhibitors for snake venom.

Evolution of Three-Finger Toxin Genes in Neotropical Colubrine Snakes (Colubridae)

Snake venom research has historically focused on front-fanged species (Viperidae and Elapidae), limiting our knowledge of venom evolution in rear-fanged snakes across their ecologically diverse phylogeny. Three-finger toxins (3FTxs) are a known neurotoxic component in the venoms of some rear-fanged snakes (Colubridae: Colubrinae), but it is unclear how prevalent 3FTxs are both in expression within venom glands and more broadly among colubrine species. Here, we used a transcriptomic approach to characterize the venom expression profiles of four species of colubrine snakes from the Neotropics that were dominated by 3FTx expression (in the genera Chironius, Oxybelis, Rhinobothryum, and Spilotes). By reconstructing the gene trees of 3FTxs, we found evidence of putative novel heterodimers in the sequences of Chironius multiventris and Oxybelis aeneus, revealing an instance of parallel evolution of this structural change in 3FTxs among rear-fanged colubrine snakes. We also found positive selection at sites within structural loops or "fingers" of 3FTxs, indicating these areas may be key binding sites that interact with prey target molecules. Overall, our results highlight the importance of exploring the venoms of understudied species in reconstructing the full evolutionary history of toxins across the tree of life.

The Venom Composition of the Snake Tribe Philodryadini: 'Omic' Techniques Reveal Intergeneric Variability among South American Racers

Snakes of the Philodryadini tribe are included in the Dipsadidae family, which is a diverse group of rear-fanged snakes widespread in different ecological conditions, including habitats and diet. However, little is known about the composition and effects of their venoms despite their relevance for understanding the evolution of these snakes or even their impact on the occasional cases of human envenoming. In this study, we integrated venom gland transcriptomics, venom proteomics and functional assays to characterize the venoms from eight species of the Philodryadini tribe, which includes the genus Philodryas, Chlorosoma and Xenoxybelis. The most abundant components identified in the venoms were snake venom metalloproteinases (SVMPs), cysteine-rich secretory proteins (CRISPs), C-type lectins (CTLs), snake endogenous matrix metalloproteinases type 9 (seMMP-9) and snake venom serinoproteinases (SVSPs). These protein families showed a variable expression profile in each genus. SVMPs were the most abundant components in Philodryas, while seMMP-9 and CRISPs were the most expressed in Chlorosoma and Xenoxybelis, respectively. Lineage-specific differences in venom composition were also observed among Philodryas species, whereas P. olfersii presented the highest amount of SVSPs and P. agassizii was the only species to express significant amounts of 3FTx. The variability observed in venom composition was confirmed by the venom functional assays. Philodryas species presented the highest SVMP activity, whereas Chlorosoma species showed higher levels of gelatin activity, which may correlate to the seMMP-9 enzymes. The variability observed in the composition of these venoms may be related to the tribe phylogeny and influenced by their diets. In the presented study, we expanded the set of venomics studies of the Philodryadini tribe, which paves new roads for further studies on the evolution and ecology of Dipsadidae snakes.

Divergent Specialization of Simple Venom Gene Profiles among Rear-Fanged Snake Genera (Helicops and Leptodeira, Dipsadinae, Colubridae)

Many venomous animals express toxins that show extraordinary levels of variation both within and among species. In snakes, most studies of venom variation focus on front-fanged species in the families Viperidae and Elapidae, even though rear-fanged snakes in other families vary along the same ecological axes important to venom evolution. Here we characterized venom gland transcriptomes from 19 snakes across two dipsadine rear-fanged genera (Leptodeira and Helicops, Colubridae) and two front-fanged genera (Bothrops, Viperidae; Micrurus, Elapidae). We compared patterns of composition, variation, and diversity in venom transcripts within and among all four genera. Venom gland transcriptomes of rear-fanged Helicops and Leptodeira and front-fanged Micrurus are each dominated by expression of single toxin families (C-type lectins, snake venom metalloproteinase, and phospholipase A2, respectively), unlike highly diverse front-fanged Bothrops venoms. In addition, expression patterns of congeners are much more similar to each other than they are to species from other genera. These results illustrate the repeatability of simple venom profiles in rear-fanged snakes and the potential for relatively constrained venom composition within genera.

Genetic characterization of potential venom resistance proteins in California ground squirrels ( Otospermophilus beecheyi ) using transcriptome analyses

Understanding the molecular basis of adaptations in coevolving species requires identifying the genes that underlie reciprocally selected phenotypes, such as those involved in venom in snakes and resistance to the venom in their prey. In this regard, California ground squirrels (CGS; Otospermophilus beecheyi) are eaten by northern Pacific rattlesnakes (Crotalus oreganus oreganus), but individual squirrels may still show substantial resistance to venom and survive bites. A recent study using proteomics identified venom interactive proteins (VIPs) in the blood serum of CGS. These VIPs represent possible resistance proteins, but the sequences of genes encoding them are unknown despite the value of such data to molecular studies of coevolution. To address this issue, we analyzed a de novo assembled transcriptome from CGS liver tissue-where many plasma proteins are synthesized-and other tissues from this species. We then examined VIP sequences in terms of three characteristics that identify them as possible resistance proteins: evidence for positive selection, high liver expression, and nonsynonymous variation across CGS populations. Based on these characteristics, we identified five VIPs (i.e., α-2-macroglobulin, α-1-antitrypsin-like protein GS55-LT, apolipoprotein A-II, hibernation-associated plasma protein HP-20, and hibernation-associated plasma protein HP-27) as the most likely candidates for resistance proteins among VIPs identified to date. Four of these proteins have been previously implicated in conferring resistance to the venom in mammals, validating our approach. When combined with the detailed information available for rattlesnake venom proteins, these results set the stage for future work focused on understanding coevolutionary interactions at the molecular level between these species.

Venom production and secretion in reptiles

The venom glands of reptiles, particularly those of front-fanged advanced snakes, must satisfy conflicting biological demands: rapid synthesis of potentially labile and highly toxic proteins, storage in the gland lumen for long periods, stabilization of the stored secretions, immediate activation of toxins upon deployment and protection of the animal from the toxic effects of its own venom. This dynamic system could serve as a model for the study of a variety of different phenomena involving exocrine gland activation, protein synthesis, stabilization of protein products and secretory mechanisms. However, these studies have been hampered by a lack of a long-term model that can be propagated in the lab (as opposed to whole-animal studies). Numerous attempts have been made to extend the lifetime of venom gland secretory cells, but only recently has an organoid model been shown to have the requisite qualities of recapitulation of the native system, self-propagation and long-term viability (>1 year). A tractable model is now available for myriad cell- and molecular-level studies of venom glands, protein synthesis and secretion. However, venom glands of reptiles are not identical, and many differ very extensively in overall architecture, microanatomy and protein products produced. This Review summarizes the similarities among and differences between venom glands of helodermatid lizards and of rear-fanged and front-fanged snakes, highlighting those areas that are well understood and identifying areas where future studies can fill in significant gaps in knowledge of these ancient, yet fascinating systems.

An integrative view of the toxic potential of Conophis lineatus (Dipsadidae: Xenodontinae), a medically relevant rear-fanged snake

Most traditional research on snake venoms has focused on front-fanged snake families (Viperidae, Elapidae, and Atractaspididae). However, venom is now generally accepted as being a much more broadly possessed trait within snakes, including species traditionally considered harmless. Unfortunately, due to historical inertia and methodological challenges, the toxin repertoires of non-front-fanged snake families (e.g., Colubridae, Dipsadidae, and Natricidae) have been heavily neglected despite the knowledge of numerous species capable of inflicting medically relevant envenomations. Integrating proteomic data for validation, we perform a de novo assembly and analysis of the Duvernoy's venom gland transcriptome of the Central American Road Guarder (Dipsadidae: Xenodontinae: Conophis lineatus), a species known for its potent bite. We identified 28 putative toxin transcripts from 13 toxin families in the Duvernoy's venom gland transcriptome, comprising 63.7% of total transcriptome expression. In addition to ubiquitous snake toxin families, we proteomically confirmed several atypical venom components. The most highly expressed toxins (55.6% of total toxin expression) were recently described snake venom matrix metalloproteases (svMMPs), with 48.0% of svMMP expression contributable to a novel svMMP isoform. We investigate the evolution of the new svMMP isoform in the context of rear-fanged snakes using phylogenetics. Finally, we examine the morphology of the venom apparatus using μCT and explore how the venom relates to autecology and the highly hemorrhagic effects seen in human envenomations. Importantly, we provide the most complete venom characterization of this medically relevant snake species to date, producing insights into the effects and evolution of its venom, and point to future research directions to better understand the venoms of 'harmless' non-front-fanged snakes.

Venom of the Annulated Sea Snake Hydrophis cyanocinctus: A Biochemically Simple but Genetically Complex Weapon

Given that the venom system in sea snakes has a role in enhancing their secondary adaption to the marine environment, it follows that elucidating the diversity and function of venom toxins will help to understand the adaptive radiation of sea snakes. We performed proteomic and de novo NGS analyses to explore the diversity of venom toxins in the annulated sea snake (Hydrophis cyanocinctus) and estimated the adaptive molecular evolution of the toxin-coding unigenes and the toxicity of the major components. We found three-finger toxins (3-FTxs), phospholipase A2 (PLA2) and cysteine-rich secretory protein (CRISP) in the venom proteome and 59 toxin-coding unigenes belonging to 24 protein families in the venom-gland transcriptome; 3-FTx and PLA2 were the most abundant families. Nearly half of the toxin-coding unigenes had undergone positive selection. The short- (i.p. 0.09 μg/g) and long-chain neurotoxin (i.p. 0.14 μg/g) presented fairly high toxicity, whereas both basic and acidic PLA2s expressed low toxicity. The toxicity of H. cyanocinctus venom was largely determined by the 3-FTxs. Our data show the venom is used by H. cyanocinctus as a biochemically simple but genetically complex weapon and venom evolution in H. cyanocinctus is presumably driven by natural selection to deal with fast-moving prey and enemies in the marine environment.

Venom-gland transcriptomic, venomic, and antivenomic profiles of the spine-bellied sea snake (Hydrophis curtus) from the South China Sea

Background

A comprehensive evaluation of the -omic profiles of venom is important for understanding the potential function and evolution of snake venom. Here, we conducted an integrated multi-omics-analysis to unveil the venom-transcriptomic and venomic profiles in a same group of spine-bellied sea snakes (Hydrophis curtus) from the South China Sea, where the snake is a widespread species and might generate regionally-specific venom potentially harmful to human activities. The capacity of two heterologous antivenoms to immunocapture the H. curtus venom was determined for an in-depth evaluation of their rationality in treatment of H. curtus envenomation. In addition, a phylogenetic analysis by maximum likelihood was used to detect the adaptive molecular evolution of full-length toxin-coding unigenes.

Results

A total of 90,909,384 pairs of clean reads were generated via Illumina sequencing from a pooled cDNA library of six specimens, and yielding 148,121 unigenes through de novo assembly. Sequence similarity searching harvested 63,845 valid annotations, including 63,789 non-toxin-coding and 56 toxin-coding unigenes belonging to 22 protein families. Three protein families, three-finger toxins (3-FTx), phospholipase A₂ (PLA₂), and cysteine-rich secretory protein, were detected in the venom proteome. 3-FTx (27.15% in the transcriptome/41.94% in the proteome) and PLA₂ (59.71%/49.36%) were identified as the most abundant families in the venom-gland transcriptome and venom proteome. In addition, 24 unigenes from 11 protein families were shown to have experienced positive selection in their evolutionary history, whereas four were relatively conserved throughout evolution. Commercial Naja atra antivenom exhibited a stronger capacity than Bungarus multicinctus antivenom to immunocapture H. curtus venom components, especially short neurotoxins, with the capacity of both antivenoms to immunocapture short neurotoxins being weaker than that for PLA₂s.

Conclusions

Our study clarified the venom-gland transcriptomic and venomic profiles along with the within-group divergence of a H. curtus population from the South China Sea. Adaptive evolution of most venom components driven by natural selection appeared to occur rapidly during evolutionary history. Notably, the utility of commercial N. atra and B. multicinctus antivenoms against H. curtus toxins was not comprehensive; thus, the development of species-specific antivenom is urgently needed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07824-7.

Duvernoy's Gland Transcriptomics of the Plains Black-Headed Snake, Tantilla nigriceps (Squamata, Colubridae): Unearthing the Venom of Small Rear-Fanged Snakes

The venoms of small rear-fanged snakes (RFS) remain largely unexplored, despite increased recognition of their importance in understanding venom evolution more broadly. Sequencing the transcriptome of venom-producing glands has greatly increased the ability of researchers to examine and characterize the toxin repertoire of small taxa with low venom yields. Here, we use RNA-seq to characterize the Duvernoy's gland transcriptome of the Plains Black-headed Snake, Tantilla nigriceps, a small, semi-fossorial colubrid that feeds on a variety of potentially dangerous arthropods including centipedes and spiders. We generated transcriptomes of six individuals from three localities in order to both characterize the toxin expression of this species for the first time, and to look for initial evidence of venom variation in the species. Three toxin families-three-finger neurotoxins (3FTxs), cysteine-rich secretory proteins (CRISPs), and snake venom metalloproteinases (SVMPIIIs)-dominated the transcriptome of T. nigriceps; 3FTx themselves were the dominant toxin family in most individuals, accounting for as much as 86.4% of an individual's toxin expression. Variation in toxin expression between individuals was also noted, with two specimens exhibiting higher relative expression of c-type lectins than any other sample (8.7-11.9% compared to <1%), and another expressed CRISPs higher than any other toxin. This study provides the first Duvernoy's gland transcriptomes of any species of Tantilla, and one of the few transcriptomic studies of RFS not predicated on a single individual. This initial characterization demonstrates the need for further study of toxin expression variation in this species, as well as the need for further exploration of small RFS venoms.

ToxCodAn: a new toxin annotator and guide to venom gland transcriptomics

MOTIVATION

Next-generation sequencing has become exceedingly common and has transformed our ability to explore nonmodel systems. In particular, transcriptomics has facilitated the study of venom and evolution of toxins in venomous lineages; however, many challenges remain. Primarily, annotation of toxins in the transcriptome is a laborious and time-consuming task. Current annotation software often fails to predict the correct coding sequence and overestimates the number of toxins present in the transcriptome. Here, we present ToxCodAn, a python script designed to perform precise annotation of snake venom gland transcriptomes. We test ToxCodAn with a set of previously curated transcriptomes and compare the results to other annotators. In addition, we provide a guide for venom gland transcriptomics to facilitate future research and use Bothrops alternatus as a case study for ToxCodAn and our guide.

RESULTS

Our analysis reveals that ToxCodAn provides precise annotation of toxins present in the transcriptome of venom glands of snakes. Comparison with other annotators demonstrates that ToxCodAn has better performance with regard to run time ($>20x$ faster), coding sequence prediction ($>3x$ more accurate) and the number of toxins predicted (generating $>4x$ less false positives). In this sense, ToxCodAn is a valuable resource for toxin annotation. The ToxCodAn framework can be expanded in the future to work with other venomous lineages and detect novel toxins.

First insights into the biochemical and toxicological characterization of venom from the Banded Cat-eyed Snake Leptodeira annulata pulchriceps

With the aim to widen the current knowledge of toxinological implications of bites from rear-fanged snakes and biological roles of their venoms, this study focuses on the biochemical composition and toxic effects of the venom of Leptodeira annulata pulchriceps from Argentina. We analyzed the protein composition by electrophoresis and mass spectrometry, and enzymatic properties by quantitative assays on different substrates. Additionally, we evaluated local and systemic toxicity in mice, and tested its cross-reactivity with elapid and viperid antivenoms used in Argentina. This venom showed features reminiscent of venoms from snakes of Bothrops genus, containing components ranging from ~17 to 75 kDa, which are mainly tissue-damaging toxins such as proteinases. Although showing low lethality to mice (LD₅₀ = 20 μg/g body weight), prominent hemorrhage developed locally in mice intramuscularly and intradermally injected with the venom, and the minimum hemorrhagic dose was found to be 12.7 μg/mouse. This study is the first comprehensive investigation of the venom of L. a. pulchriceps, and sheds new light on differences between this and those of the other two subspecies of L. annulata. Additionally, the study provides new insights into the venom components of "colubrid" snakes, advocating for considering bites from this rich diversity of snakes as a public health problem that needs to be addressed worldwide.

Size Matters: An Evaluation of the Molecular Basis of Ontogenetic Modifications in the Composition of Bothrops jararacussu Snake Venom

Ontogenetic changes in venom composition have been described in Bothrops snakes, but only a few studies have attempted to identify the targeted paralogues or the molecular mechanisms involved in modifications of gene expression during ontogeny. In this study, we decoded B. jararacussu venom gland transcripts from six specimens of varying sizes and analyzed the variability in the composition of independent venom proteomes from 19 individuals. We identified 125 distinct putative toxin transcripts, and of these, 73 were detected in venom proteomes and only 10 were involved in the ontogenetic changes. Ontogenetic variability was linearly related to snake size and did not correspond to the maturation of the reproductive stage. Changes in the transcriptome were highly predictive of changes in the venom proteome. The basic myotoxic phospholipases A₂ (PLA₂s) were the most abundant components in larger snakes, while in venoms from smaller snakes, PIII-class SVMPs were the major components. The snake venom metalloproteinases (SVMPs) identified corresponded to novel sequences and conferred higher pro-coagulant and hemorrhagic functions to the venom of small snakes. The mechanisms modulating venom variability are predominantly related to transcriptional events and may consist of an advantage of higher hematotoxicity and more efficient predatory function in the venom from small snakes.

Venoms of New World Vinesnakes (Oxybelis aeneus and O. fulgidus)

Species of Oxybelis are extremely elongate arboreal snakes that are broadly distributed in the Americas, from extreme southeastern Arizona (USA) to central South America. Primarily feeding on lizards and birds, Oxybelis venoms are poorly known in general, but a prominent taxon-specific three-finger toxin (fulgimotoxin) was isolated from and is a prominent component of O. fulgidus venom; a homolog is also present in O. aeneus venom. As part of ongoing characterization of venoms from rear-fanged snakes, we describe here the composition of two broadly distributed species, O. aeneus and O. fulgidus. Venom proteomes were of very low complexity, and four protein families (LAAO, PIII SVMP, CRiSP and 3FTx) account for more than 90% of total protein composition. Venoms from both species are moderately toxic to mice and to Hemidactylus geckos, but they are nearly an order of magnitude more toxic to Anolis lizards (a native prey species). These results reflect a trend in colubrid venom composition that is becoming increasingly more common: the presence of taxon-specific toxins, specifically three-finger toxins, preferentially targeting lizards and/or birds.

Spontaneous Tumor Regression in Tasmanian Devils Associated with RASL11A Activation

Spontaneous tumor regression has been documented in a small proportion of human cancer patients, but the specific mechanisms underlying tumor regression without treatment are not well understood. Tasmanian devils are threatened with extinction from a transmissible cancer due to universal susceptibility and a near 100% case fatality rate. In over 10,000 cases, <20 instances of natural tumor regression have been detected. Previous work in this system has focused on Tasmanian devil genetic variation associated with the regression phenotype. Here, we used comparative and functional genomics to identify tumor genetic variation associated with tumor regression. We show that a single point mutation in the 5' untranslated region of the putative tumor suppressor RASL11A significantly contributes to tumor regression. RASL11A was expressed in regressed tumors but silenced in wild-type, nonregressed tumors, consistent with RASL11A downregulation in human cancers. Induced RASL11A expression significantly reduced tumor cell proliferation in vitro The RAS pathway is frequently altered in human cancers, and RASL11A activation may provide a therapeutic treatment option for Tasmanian devils as well as a general mechanism for tumor inhibition.

Proteomics of Select Neglected Tropical Diseases

Technological advances in mass spectrometry have enabled the extensive identification, characterization, and quantification of proteins in any biological system. In disease processes proteins are often altered in response to external stimuli; therefore, proteomics, the large-scale study of proteins and their functions, represents an invaluable tool for understanding the molecular basis of disease. This review highlights the use of mass spectrometry-based proteomics to study the pathogenesis, etiology, and pathology of several neglected tropical diseases (NTDs), a diverse group of disabling diseases primarily associated with poverty in tropical and subtropical regions of the world. While numerous NTDs have been the subject of proteomic studies, this review focuses on Buruli ulcer, dengue, leishmaniasis, and snakebite envenoming. The proteomic studies highlighted provide substantial information on the pathogenic mechanisms driving these diseases; they also identify molecular targets for drug discovery and development and uncover promising biomarkers that can assist in early diagnosis.

Many Options, Few Solutions: Over 60 My Snakes Converged on a Few Optimal Venom Formulations

Gene expression changes contribute to complex trait variations in both individuals and populations. However, the evolution of gene expression underlying complex traits over macroevolutionary timescales remains poorly understood. Snake venoms are proteinaceous cocktails where the expression of each toxin can be quantified and mapped to a distinct genomic locus and traced for millions of years. Using a phylogenetic generalized linear mixed model, we analyzed expression data of toxin genes from 52 snake species spanning the 3 venomous snake families and estimated phylogenetic covariance, which acts as a measure of evolutionary constraint. We find that evolution of toxin combinations is not constrained. However, although all combinations are in principle possible, the actual dimensionality of phylomorphic space is low, with envenomation strategies focused around only four major toxin families: metalloproteases, three-finger toxins, serine proteases, and phospholipases A2. Although most extant snakes prioritize either a single or a combination of major toxin families, they are repeatedly recruited and lost. We find that over macroevolutionary timescales, the venom phenotypes were not shaped by phylogenetic constraints, which include important microevolutionary constraints such as epistasis and pleiotropy, but more likely by ecological filtering that permits a small number of optimal solutions. As a result, phenotypic optima were repeatedly attained by distantly related species. These results indicate that venoms evolve by selection on biochemistry of prey envenomation, which permit diversity through parallelism, and impose strong limits, since only a few of the theoretically possible strategies seem to work well and are observed in extant snakes.

Antimicrobial Activity of Protein Fraction from Naja ashei Venom Against Staphylococcus epidermidis

One of the key problems of modern infectious disease medicine is the growing number of drug-resistant and multi-drug-resistant bacterial strains. For this reason, many studies are devoted to the search for highly active antimicrobial substances that could be used in therapy against bacterial infections. As it turns out, snake venoms are a rich source of proteins that exert a strong antibacterial effect, and therefore they have become an interesting research material. We analyzed Naja ashei venom for such antibacterial properties, and we found that a specific composition of proteins can act to eliminate individual bacterial cells, as well as the entire biofilm of Staphylococcus epidermidis. In general, we used ion exchange chromatography (IEX) to obtain 10 protein fractions with different levels of complexity, which were then tested against certified and clinical strains of S. epidermidis. One of the fractions (F2) showed exceptional antimicrobial effects both alone and in combination with antibiotics. The protein composition of the obtained fractions was determined using mass spectrometry techniques, indicating a high proportion of phospholipases A₂, three-finger toxins, and L-amino acids oxidases in F2 fraction, which are most likely responsible for the unique properties of this fraction. Moreover, we were able to identify a new group of low abundant proteins containing the Ig-like domain that have not been previously described in snake venoms.

Omics Technologies for Profiling Toxin Diversity and Evolution in Snake Venom: Impacts on the Discovery of Therapeutic and Diagnostic Agents

Snake venoms are primarily composed of proteins and peptides, and these toxins have developed high selectivity to their biological targets. This makes venoms interesting for exploration into protein evolution and structure-function relationships. A single venom protein superfamily can exhibit a variety of pharmacological effects; these variations in activity originate from differences in functional sites, domains, posttranslational modifications, and the formations of toxin complexes. In this review, we discuss examples of how the major venom protein superfamilies have diversified, as well as how newer technologies in the omics fields, such as genomics, transcriptomics, and proteomics, can be used to characterize both known and unknown toxins.Because toxins are bioactive molecules with a rich diversity of activities, they can be useful as therapeutic and diagnostic agents, and successful examples of toxin applications in these areas are also reviewed. With the current rapid pace of technology, snake venom research and its applications will only continue to expand.

Transcriptomic-Proteomic Correlation in the Predation-Evoked Venom of the Cone Snail, Conus imperialis

Individual variation in animal venom has been linked to geographical location, feeding habit, season, size, and gender. Uniquely, cone snails possess the remarkable ability to change venom composition in response to predatory or defensive stimuli. To date, correlations between the venom gland transcriptome and proteome within and between individual cone snails have not been reported. In this study, we use 454 pyrosequencing and mass spectrometry to decipher the transcriptomes and proteomes of the venom gland and corresponding predation-evoked venom of two specimens of Conus imperialis. Transcriptomic analyses revealed 17 conotoxin gene superfamilies common to both animals, including 5 novel superfamilies and two novel cysteine frameworks. While highly expressed transcripts were common to both specimens, variation of moderately and weakly expressed precursor sequences was surprisingly diverse, with one specimen expressing two unique gene superfamilies and consistently producing more paralogs within each conotoxin gene superfamily. Using a quantitative labelling method, conotoxin variability was compared quantitatively, with highly expressed peptides showing a strong correlation between transcription and translation, whereas peptides expressed at lower levels showed a poor correlation. These results suggest that major transcripts are subject to stabilizing selection, while minor transcripts are subject to diversifying selection.

A Single Mutation Unlocks Cascading Exaptations in the Origin of a Potent Pitviper Neurotoxin

Evolutionary innovations and complex phenotypes seemingly require an improbable amount of genetic change to evolve. Rattlesnakes display two dramatically different venom phenotypes. Type I venoms are hemorrhagic with low systemic toxicity and high expression of tissue-destroying snake venom metalloproteinases. Type II venoms are highly neurotoxic and lack snake venom metalloproteinase expression and associated hemorrhagic activity. This dichotomy hinges on Mojave toxin (MTx), a phospholipase A2 (PLA2) based β-neurotoxin expressed in Type II venoms. MTx is comprised of a nontoxic acidic subunit that undergoes extensive proteolytic processing and allosterically regulates activity of a neurotoxic basic subunit. Evolution of the acidic subunit presents an evolutionary challenge because the need for high expression of a nontoxic venom component and the proteolytic machinery required for processing suggests genetic changes of seemingly little immediate benefit to fitness. We showed that MTx evolved through a cascading series of exaptations unlocked by a single nucleotide change. The evolution of one new cleavage site in the acidic subunit unmasked buried cleavage sites already present in ancestral PLA2s, enabling proteolytic processing. Snake venom serine proteases, already present in the venom to disrupt prey hemostasis, possess the requisite specificities for MTx acidic subunit proteolysis. The dimerization interface between MTx subunits evolved by exploiting a latent, but masked, hydrophobic interaction between ancestral PLA2s. The evolution of MTx through exaptation of existing functional and structural features suggests complex phenotypes that depend on evolutionary innovations can arise from minimal genetic change enabled by prior evolution.

The Harderian gland transcriptomes of Caraiba andreae, Cubophis cantherigerus and Tretanorhinus variabilis, three colubroid snakes from Cuba

The Harderian gland is a cephalic structure, widely distributed among vertebrates. In snakes, the Harderian gland is anatomically connected to the vomeronasal organ via the nasolacrimal duct, and in some species can be larger than the eyes. The function of the Harderian gland remains elusive, but it has been proposed to play a role in the production of saliva, pheromones, thermoregulatory lipids and growth factors, among others. Here, we have profiled the transcriptomes of the Harderian glands of three non-front-fanged colubroid snakes from Cuba: Caraiba andreae (Cuban Lesser Racer); Cubophis cantherigerus (Cuban Racer); and Tretanorhinus variabilis (Caribbean Water Snake), using Illumina HiSeq2000 100 bp paired-end. In addition to ribosomal and non-characterized proteins, the most abundant transcripts encode putative transport/binding, lipocalin/lipocalin-like, and bactericidal/permeability-increasing-like proteins. Transcripts coding for putative canonical toxins described in venomous snakes were also identified. This transcriptional profile suggests a more complex function than previously recognized for this enigmatic organ.

Defining the pathogenic threat of envenoming by South African shield-nosed and coral snakes (genus Aspidelaps), and revealing the likely efficacy of available antivenom

While envenoming by the southern African shield-nosed or coral snakes (genus Aspidelaps) has caused fatalities, bites are uncommon. Consequently, this venom is not used in the mixture of snake venoms used to immunise horses for the manufacture of regional SAIMR (South African Institute for Medical Research) polyvalent antivenom. Aspidelaps species are even excluded from the manufacturer's list of venomous snakes that can be treated by this highly effective product. This leaves clinicians, albeit rarely, in a therapeutic vacuum when treating envenoming by these snakes. This is a significantly understudied small group of nocturnal snakes and little is known about their venom compositions and toxicities. Using a murine preclinical model, this study determined that the paralysing toxicity of venoms from Aspidelaps scutatus intermedius, A. lubricus cowlesi and A. l. lubricus approached that of venoms from highly neurotoxic African cobras and mambas. This finding was consistent with the cross-genus dominance of venom three-finger toxins, including numerous isoforms which showed extensive interspecific variation. Our comprehensive analysis of venom proteomes showed that the three Aspidelaps species possess highly similar venom proteomic compositions. We also revealed that the SAIMR polyvalent antivenom cross-reacted extensively in vitro with venom proteins of the three Aspidelaps. Importantly, this cross-genus venom-IgG binding translated to preclinical (in a murine model) neutralisation of A. s. intermedius venom-induced lethality by the SAIMR polyvalent antivenom, at doses comparable with those that neutralise venom from the cape cobra (Naja nivea), which the antivenom is directed against. Our results suggest a wider than anticipated clinical utility of the SAIMR polyvalent antivenom, and here we seek to inform southern African clinicians that this readily available antivenom is likely to prove effective for victims of Aspidelaps envenoming. BIOLOGICAL SIGNIFICANCE: Coral and shield-nosed snakes (genus Aspidelaps) comprise two species and several subspecies of potentially medically important venomous snakes distributed in Namibia, Botswana, Zimbabwe, Mozambique and South Africa. Documented human fatalities, although rare, have occurred from both A. lubricus and A. scutatus. However, their venom proteomes and the pathological effects of envenomings by this understudied group of nocturnal snakes remain uncharacterised. Furthermore, no commercial antivenom is made using venom from species of the genus Aspidelaps. To fill this gap, we have conducted a transcriptomics-guided comparative proteomics analysis of the venoms of the intermediate shield-nose snake (A. s. intermedius), southern coral snake (A. l. lubricus), and Cowle's shield snake (A. l. cowlesi); investigated the mechanism of action underpinning lethality by A. s. intermedius in the murine model; and assessed the in vitro immunoreactivity of the SAIMR polyvalent antivenom towards the venom toxins of A. l. lubricus and A. l. cowlesi, and the in vivo capability of this antivenom at neutralising the lethal effect of A. s. intermedius venom. Our data revealed a high degree of conservation of the global composition of the three Aspidelaps venom proteomes, all characterised by the overwhelming predominance of neurotoxic 3FTxs, which induced classical signs of systemic neurotoxicity in mice. The SAIMR polyvalent antivenom extensively binds to Aspidelaps venom toxins and neutralised, with a potency of 0.235 mg venom/mL antivenom, the lethal effect of A. s. intermedius venom. Our data suggest that the SAIMR antivenom could be a useful therapeutic tool for treating human envenomings by Aspidelaps species.

Proteomic Deep Mining the Venom of the Red-Headed Krait, Bungarus flaviceps

The use of -omics technologies allows for the characterization of snake venom composition at a fast rate and at high levels of detail. In the present study, we investigated the protein content of Red-headed Krait (Bungarus flaviceps) venom. This analysis revealed a high diversity of snake venom protein families, as evidenced by high-throughput mass spectrometric analysis. We found all six venom protein families previously reported in a transcriptome study of the venom gland of B. flaviceps, including phospholipases A₂ (PLA₂s), Kunitz-type serine proteinase inhibitors (KSPIs), three-finger toxins (3FTxs), cysteine-rich secretory proteins (CRISPs), snaclecs, and natriuretic peptides. A combined approach of automated database searches and de novo sequencing of tandem mass spectra, followed by sequence similarity searches, revealed the presence of 12 additional toxin families. De novo sequencing alone was able to identify 58 additional peptides, and this approach contributed significantly to the comprehensive description of the venom. Abundant protein families comprise 3FTxs (22.3%), KSPIs (19%), acetylcholinesterases (12.6%), PLA₂s (11.9%), venom endothelial growth factors (VEGFs, 8.4%), nucleotidases (4.3%), and C-type lectin-like proteins (snaclecs, 3.3%); an additional 11 toxin families are present at significantly lower concentrations, including complement depleting factors, a family not previously detected in Bungarus venoms. The utility of a multifaceted approach toward unraveling the proteome of snake venoms, employed here, allowed detection of even minor venom components. This more in-depth knowledge of the composition of B. flaviceps venom facilitates a better understanding of snake venom molecular evolution, in turn contributing to more effective treatment of krait bites.

Three-Finger Toxin Diversification in the Venoms of Cat-Eye Snakes (Colubridae: Boiga)

The Asian genus Boiga (Colubridae) is among the better studied non-front-fanged snake lineages, because their bites have minor, but noticeable, effects on humans. Furthermore, B. irregularis has gained worldwide notoriety for successfully invading Guam and other nearby islands with drastic impacts on the local bird populations. One of the factors thought to allow B. irregularis to become such a noxious pest is irditoxin, a dimeric neurotoxin composed of two three-finger toxins (3FTx) joined by a covalent bond between two newly evolved cysteines. Irditoxin is highly toxic to diapsid (birds and reptiles) prey, but roughly 1000 × less potent to synapsids (mammals). Venom plays an important role in the ecology of all species of Boiga, but it remains unknown if any species besides B. irregularis produce irditoxin-like dimeric toxins. In this study, we use transcriptomic analyses of venom glands from five species [B. cynodon, B. dendrophila dendrophila, B. d. gemmicincta, B. irregularis (Brisbane population), B. irregularis (Sulawesi population), B. nigriceps, B. trigonata] and proteomic analyses of B. d. dendrophila and a representative of the sister genus Toxicodryas blandingii to investigate the evolutionary history of 3FTx within Boiga and its close relative. We found that 92.5% of Boiga 3FTx belong to a single clade which we refer to as denmotoxin-like because of the close relation between these toxins and the monomeric denmotoxin according to phylogenetic, sequence clustering, and protein similarity network analyses. We show for the first time that species beyond B. irregularis secrete 3FTx with additional cysteines in the same position as both the A and B subunits of irditoxin. Transcripts with the characteristic mutations are found in B. d. dendrophila, B. d. gemmicincta, B. irregularis (Brisbane population), B. irregularis (Sulawesi population), and B. nigriceps. These results are confirmed by proteomic analyses that show direct evidence of dimerization within the venom of B. d. dendrophila, but not T. blandingii. Our results also suggest the possibility of novel dimeric toxins in other genera such as Telescopus and Trimorphodon. All together, this suggests that the origin of these peculiar 3FTx is far earlier than was appreciated and their evolutionary history has been complex.

Transcriptome-facilitated proteomic characterization of rear-fanged snake venoms reveal abundant metalloproteinases with enhanced activity

High-throughput technologies were used to identify venom gland toxin expression and to characterize the venom proteomes of two rear-fanged snakes, Ahaetulla prasina (Asian Green Vine Snake) and Borikenophis portoricensis (Puerto Rican Racer). Sixty-nine complete toxin-coding transcripts from 12 venom protein superfamilies (A. prasina) and 50 complete coding transcripts from 11 venom protein superfamilies (B. portoricensis) were identified in the venom glands. However, only 18% (A. prasina) and 32% (B. portoricensis) of the translated protein isoforms were detected in the proteome of these venoms. Both venom gland transcriptomes and venom proteomes were dominated by P-III metalloproteinases. Three-finger toxins, cysteine-rich secretory proteins, and C-type lectins were present in moderate amounts, but other protein superfamilies showed very low abundances. Venoms contained metalloproteinase activity comparable to viperid snake venom levels, but other common venom enzymes were absent or present at negligible levels. Western blot analysis showed metalloproteinase and cysteine-rich secretory protein epitopes shared with the highly venomous Boomslang (Dispholidus typus). The abundance of metalloproteinases emphasizes the important trophic role of these toxins. Comprehensive, transcriptome-informed definition of proteomes and functional characterization of venom proteins in rear-fanged snake families help to elucidate toxin evolution and provide models for protein structure-function analyses.

Adaptive evolution of distinct prey-specific toxin genes in rear-fanged snake venom

Venom proteins evolve rapidly, and as a trophic adaptation are excellent models for predator-prey evolutionary studies. The key to a deeper understanding of venom evolution is an integrated approach, combining prey assays with analysis of venom gene expression and venom phenotype. Here, we use such an approach to study venom evolution in the Amazon puffing snake, Spilotes sulphureus, a generalist feeder. We identify two novel three-finger toxins: sulditoxin and sulmotoxin 1. These new toxins are not only two of the most abundant venom proteins, but are also functionally intriguing, displaying distinct prey-specific toxicities. Sulditoxin is highly toxic towards lizard prey, but is non-toxic towards mammalian prey, even at greater than 22-fold higher dosage. By contrast, sulmotoxin 1 exhibits the reverse trend. Furthermore, evolutionary analysis and structural modelling show highest sequence variability in the central loop of these proteins, probably driving taxon-specific toxicity. This is, to our knowledge, the first case in which a bimodal and contrasting pattern of toxicity has been shown for proteins in the venom of a single snake in relation to diet. Our study is an example of how toxin gene neofunctionalization can result in a venom system dominated by one protein superfamily and still exhibit flexibility in prey capture efficacy.

Transcriptomics-guided bottom-up and top-down venomics of neonate and adult specimens of the arboreal rear-fanged Brown Treesnake, Boiga irregularis, from Guam

The Brown Treesnake (Boiga irregularis) is an arboreal, nocturnal, rear-fanged venomous snake native to northern and eastern regions of Australia, Papua New Guinea and the Solomon Islands. It was inadvertently introduced onto the island of Guam during the late 1940's to early 1950's, and it has caused massive declines and extirpations of the native bird, lizard, and mammal populations. In the current study, we report the characterization of the venom proteome of an adult and a neonate B. irregularis specimens from Guam by a combination of venom gland transcriptomic and venomic analyses. Venom gland transcriptomic analysis of an adult individual identified toxins belonging to 18 protein families, with three-finger toxin isoforms being the most abundantly expressed transcripts, comprising 94% of all venom protein transcript reads. Transcripts for PIII-metalloproteinases, C-type lectins, cysteine-rich secretory proteins, acetylcholinesterases, natriuretic peptides, ficolins, phospholipase A₂ (PLA₂) inhibitors, PLA₂s, vascular endothelial growth factors, Kunitz-type protease inhibitors, cystatins, phospholipase Bs, cobra venom factors, waprins, SVMP inhibitors, matrix metalloproteinases, and hyaluronidases were also identified, albeit, at very low abundances ranging from 0.05% to 1.7% of the transcriptome. The venom proteomes of neonate and adult B. irregularis were also both overwhelmingly (78 and 84%, respectively) dominated by monomeric and dimeric 3FTxs, followed by moderately abundant (21% (N) and 13% (A)) CRISPs, low abundance (1% (N), 3% (A)) PIII-SVMPs, and very low abundance (<0.01%) PLA₂ and SVMP inhibitors. The differences in relative toxin abundances identified between neonate and adult snakes likely correlates to shifts in prey preference between the two age classes, from nearly-exclusively lizards to lizards, birds and small mammals. Immunoaffinity antivenomics with experimentally designed rabbit anti-Brown Treesnake (anti-BTS) venom IgGs against homologous venom from adult snakes demonstrated that CRISPs, PIII-SVMPs, and 60-70% of 3FTxs were effectively immunocaptured. Western blot analysis showed that all venom proteins were recognized by anti-BTS IgGs, and cross-reactivity with other rear-fanged snake venoms was also observed. Incubation of anti-BTS venom IgGs with crude B. irregularis venom resulted in a significant decrease in proteolytic (SVMP) activity against azocasein. These results provide the first comparative venomic and anti-venomic analysis of neonate and adult B. irregularis from Guam, further highlighting evolutionary trends in venom composition among rear-fanged venomous snakes.

SIGNIFICANCE PARAGRAPH

The Brown Treesnake (Boiga irregularis) has caused extensive ecological and economic damage to the island of Guam where it has become a classic example of the negative impacts of invasive species. In the current study, we report the first combined transcriptomic and proteomic analysis of B. irregularis venom of Guam origin. The transcriptome of an adult snake contained toxin sequences belonging to 18 protein families, with three-finger toxin (3FTx) isoforms being the most abundant and representing 94% of all venom protein transcript reads. Our bottom-up and top-down venomic analyses confirmed that 3FTxs are the major components of B. irregularis venom, and a comparative analysis of neonate and adult venoms demonstrate a clear ontogenetic shift in toxin abundance, likely driven by dietary variation between the two age classes. Second-generation antivenomics and Western blot analysis using purified anti-Brown Treesnake rabbit serum IgGs (anti-BTS IgGs) showed strong immunoreactivity toward B. irregularis venom. Interestingly, our anti-BTS IgGs did not cross-react with 3FTxs found in several other rear-fanged snake venoms, or against 3FTxs in the venom of the elapid Ophiophagus hannah, indicating that epitopes in these 3FTx molecules are quite distinct.

Targeting Metastasis with Snake Toxins: Molecular Mechanisms

Metastasis involves the migration of cancer cells from a primary tumor to invade and establish secondary tumors in distant organs, and it is the main cause for cancer-related deaths. Currently, the conventional cytostatic drugs target the proliferation of malignant cells, being ineffective in metastatic disease. This highlights the need to find new anti-metastatic drugs. Toxins isolated from snake venoms are a natural source of potentially useful molecular scaffolds to obtain agents with anti-migratory and anti-invasive effects in cancer cells. While there is greater evidence concerning the mechanisms of cell death induction of several snake toxin classes on cancer cells; only a reduced number of toxin classes have been reported on (i.e., disintegrins/disintegrin-like proteins, C-type lectin-like proteins, C-type lectins, serinproteases, cardiotoxins, snake venom cystatins) as inhibitors of adhesion, migration, and invasion of cancer cells. Here, we discuss the anti-metastatic mechanisms of snake toxins, distinguishing three targets, which involve (1) inhibition of extracellular matrix components-dependent adhesion and migration, (2) inhibition of epithelial-mesenchymal transition, and (3) inhibition of migration by alterations in the actin/cytoskeleton network.

Venomics of Remipede Crustaceans Reveals Novel Peptide Diversity and Illuminates the Venom's Biological Role

We report the first integrated proteomic and transcriptomic investigation of a crustacean venom. Remipede crustaceans are the venomous sister group of hexapods, and the venom glands of the remipede Xibalbanus tulumensis express a considerably more complex cocktail of proteins and peptides than previously thought. We identified 32 venom protein families, including 13 novel peptide families that we name xibalbins, four of which lack similarities to any known structural class. Our proteomic data confirm the presence in the venom of 19 of the 32 families. The most highly expressed venom components are serine peptidases, chitinase and six of the xibalbins. The xibalbins represent Inhibitory Cystine Knot peptides (ICK), a double ICK peptide, peptides with a putative Cystine-stabilized α-helix/β-sheet motif, a peptide similar to hairpin-like β-sheet forming antimicrobial peptides, two peptides related to different hormone families, and four peptides with unique structural motifs. Remipede venom components represent the full range of evolutionary recruitment frequencies, from families that have been recruited into many animal venoms (serine peptidases, ICKs), to those having a very narrow taxonomic range (double ICKs), to those unique for remipedes. We discuss the most highly expressed venom components to shed light on their possible functional significance in the predatory and defensive use of remipede venom, and to provide testable ideas for any future bioactivity studies.

Comparative venom gland transcriptomics of Naja kaouthia (monocled cobra) from Malaysia and Thailand: elucidating geographical venom variation and insights into sequence novelty

Background

The monocled cobra (Naja kaouthia) is a medically important venomous snake in Southeast Asia. Its venom has been shown to vary geographically in relation to venom composition and neurotoxic activity, indicating vast diversity of the toxin genes within the species. To investigate the polygenic trait of the venom and its locale-specific variation, we profiled and compared the venom gland transcriptomes of N. kaouthia from Malaysia (NK-M) and Thailand (NK-T) applying next-generation sequencing (NGS) technology.

Methods

The transcriptomes were sequenced on the Illumina HiSeq platform, assembled and followed by transcript clustering and annotations for gene expression and function. Pairwise or multiple sequence alignments were conducted on the toxin genes expressed. Substitution rates were studied for the major toxins co-expressed in NK-M and NK-T.

Results and discussion

The toxin transcripts showed high redundancy (41–82% of the total mRNA expression) and comprised 23 gene families expressed in NK-M and NK-T, respectively (22 gene families were co-expressed). Among the venom genes, three-finger toxins (3FTxs) predominated in the expression, with multiple sequences noted. Comparative analysis and selection study revealed that 3FTxs are genetically conserved between the geographical specimens whilst demonstrating distinct differential expression patterns, implying gene up-regulation for selected principal toxins, or alternatively, enhanced transcript degradation or lack of transcription of certain traits. One of the striking features that elucidates the inter-geographical venom variation is the up-regulation of α-neurotoxins (constitutes ∼80.0% of toxin’s fragments per kilobase of exon model per million mapped reads (FPKM)), particularly the long-chain α-elapitoxin-Nk2a (48.3%) in NK-T but only 1.7% was noted in NK-M. Instead, short neurotoxin isoforms were up-regulated in NK-M (46.4%). Another distinct transcriptional pattern observed is the exclusively and abundantly expressed cytotoxin CTX-3 in NK-T. The findings suggested correlation with the geographical variation in proteome and toxicity of the venom, and support the call for optimising antivenom production and use in the region. Besides, the current study uncovered full and partial sequences of numerous toxin genes from N. kaouthia which have not been reported hitherto; these include N. kaouthia-specific l-amino acid oxidase (LAAO), snake venom serine protease (SVSP), cystatin, acetylcholinesterase (AChE), hyaluronidase (HYA), waprin, phospholipase B (PLB), aminopeptidase (AP), neprilysin, etc. Taken together, the findings further enrich the snake toxin database and provide deeper insights into the genetic diversity of cobra venom toxins.

Molecular and Cellular Mechanisms of Apoptosis during Dissociated Spermatogenesis

Apoptosis is a tightly controlled process by which tissues eliminate unwanted cells. Spontaneous germ cell apoptosis in testis has been broadly investigated in mammals that have an associated spermatogenesis pattern. However, the mechanism of germ cell apoptosis in seasonally breeding reptiles following a dissociated spermatogenesis has remained enigmatic. In the present study, morphological evidence has clearly confirmed the dissociated spermatogenesis pattern in Pelodiscus sinensis. TUNEL and TEM analyses presented dynamic changes and ultrastructural characteristics of apoptotic germ cells during seasonal spermatogenesis, implying that apoptosis might be one of the key mechanisms to clear degraded germ cells. Furthermore, using RNA-Seq and digital gene expression (DGE) profiling, a large number of apoptosis-related differentially expressed genes (DEGs) at different phases of spermatogenesis were identified and characterized in the testis. DGE and RT-qPCR analysis revealed that the critical anti-apoptosis genes, such as Bcl-2, BAG1, and BAG5, showed up-regulated patterns during intermediate and late spermatogenesis. Moreover, the increases in mitochondrial transmembrane potential in July and October were detected by JC-1 staining. Notably, the low protein levels of pro-apoptotic cleaved caspase-3 and CytC in cytoplasm were detected by immunohistochemistry and western blot analyses, indicating that the CytC-Caspase model might be responsible for the effects of germ cell apoptosis on seasonal spermatogenesis. These results facilitate understanding the regulatory mechanisms of apoptosis during spermatogenesis and uncovering the biological process of the dissociated spermatogenesis system in reptiles.

Translationally Controlled Tumor Protein (TCTP/HRF) in Animal Venoms

Proteins from TCTP/HRF family were identified as venom toxins of spiders from different genus. We have found a TCTP toxin in the venom gland of Loxosceles intermedia, a venomous spider very common in South Brazil. TCTP from L. intermedia, named LiTCTP, was cloned, produced in a heterologous prokaryotic system, and the recombinant toxin was biochemically characterized. Our results point that LiTCTP is involved in the inflammatory events of Loxocelism, the clinical signs triggered after Loxosceles sp. bite, which include intense inflammatory reaction at the bite site followed by local necrosis. TCTP toxins were also identified in spiders from different genus. There are very few articles about TCTP toxins in other venomous animals in the literature, although a NCBI database search on the protein sequences reveals TCTP on snake's venom glands transcriptomic and genomic studies. Studies on TCTP as a venom toxin are very few and its biological role as a venom component in prey capture is still unknown.

Fatal self-envenomation in a brown tree snake, Boiga irregularis, from south-east Queensland

The case history and clinical signs of a fatal self-envenomation event by a brown tree snake, Boiga irregularis, in South-east Queensland, Australia, are presented. Clinical signs began 20 minutes post-envenomation with muscle twitching, ataxia, and heat seeking behavior which progressed to partial paralysis by 6 hours, generalised paralysis and respiratory arrest at 10 hours and cardiac arrest by 12 hours post-envenomation. Clinical signs are suggestive of potent neurotoxicity for B. irregularis to its own venom.

Proteomic analysis of venom variability and ontogeny across the arboreal palm-pitvipers (genus Bothriechis)

Bothriechis is a genus of eleven currently recognized slender and arboreal venomous snakes, commonly called palm-pitvipers that range from southern Mexico to northern South America. Despite dietary studies suggesting that palm-pitvipers are generalists with an ontogenetic shift toward endothermic prey, venom proteomic analyses have revealed remarkable divergence between the venoms of the Costa Rican species, B. lateralis, B. schlegelii, B. supraciliaris, and B. nigroviridis. To achieve a more complete picture of the venomic landscape across Bothriechis, the venom proteomes of biodiversity of the northern Middle American highland palm-pitvipers, B. thalassinus, B. aurifer, and B. bicolor from Guatemala, B. marchi from Honduras, and neonate Costa Rican B. lateralis and B. schlegelii, were investigated. B. thalassinus and B. aurifer venoms are comprised by similar toxin arsenals dominated by SVMPs (33-39% of the venom proteome), CTLs (11-16%), BPP-like molecules (10-13%), and CRISPs (5-10%), and are characterized by the absence of PLA₂ proteins. Conversely, the predominant (35%) components of B. bicolor are D49-PLA₂ molecules. The venom proteome of B. marchi is similar to B. aurifer and B. thalassinus in that it is rich in SVMPs and BPPs, but also contains appreciable amounts (14.3%) of PLA₂s. The major toxin family found in the venoms of both neonate B. lateralis and B. schlegelii, is serine proteinase (SVSP), comprising about 20% of their toxin arsenals. The venom of neonate B. schlegelii is the only palm-pitviper venom where relative high amounts of Kunitz-type (6.3%) and γPLA₂ (5.2%) inhibitors have been identified. Despite notable differences between their proteomes, neonate venoms are more similar to each other than to adults of their respective species. However, the ontogenetic changes taking place in the venom of B. lateralis strongly differ from those that occur in the venom of B. schlegelii. Thus, the ontogenetic change in B. lateralis produces a SVMP-rich venom, whereas in B. schlegelii the age-dependent compositional shift generates a PLA₂-rich venom. Overall, genus-wide venomics illustrate the high evolvability of palm-pitviper venoms. The integration of the pattern of venom variation across Bothriechis into a phylogenetic and biogeographic framework may lay the foundation for assessing, in future studies, the evolutionary path that led to the present-day variability of the venoms of palm-pitvipers.

SIGNIFICANCE

Bothriechis represents a monophyletic basal genus of eleven arboreal palm-pitvipers that range from southern Mexico to northern South America. Despite palm-pitvipers' putative status as diet generalists, previous proteomic analyses have revealed remarkable divergence between the venoms of Costa Rican species, B. lateralis, B. schlegelii, B. supraciliaris, and B. nigroviridis. Our current proteomic study of Guatemalan species, B. thalassinus, B. aurifer, and B. bicolor, Honduran B. marchi, and neonate B. lateralis and B. schlegelii from Costa Rica was undertaken to deepen our understanding of the evolutionary pattern of venom proteome diversity across Bothriechis. Ancestral characters are often, but not always, preserved in an organism's development. Venoms of neonate B. lateralis and B. schlegelii are more similar to each other than to adults of their respective species, suggesting that the high evolvability of palm-pitviper venoms may represent an inherent feature of Bothriechis common ancestor. Our genus-wide data identified four nodes of venom phenotype differentiation across the phylogeny of Bothriechis. Integrated into a phylogenetic and biogeographic framework, the pattern of venom variation across Bothriechis may lay the groundwork to establish whether divergence was driven by selection for efficient resource exploitation in arboreal 'islands', thereby contributing to the ecological speciation of the genus.

RNA sequence analyses of r-Moj-DM treated cells: TXNIP is required to induce apoptosis of SK-Mel-28

RNA sequencing of untreated and r-Moj-DM treated SK-Mel-28 cells was performed after 6 h, to begin unraveling the apoptotic pathway induced by r-Moj-DM. Bioinformatic analyses of RNA sequencing data yielded 40 genes that were differentially expressed. Nine genes were upregulated and 31 were downregulated. qRT-PCR was used to validate differential expression of 13 genes with known survival or apoptotic-inducing activities. Expression of BNiP3, IGFBP3, PTPSF, Prune 2, TGF-ß, and TXNIP were compared from cells treated with r-Moj-DN (a strong apoptotic inducer) or r-Moj-DA (a non-apoptotic inducer) for 1 h, 2 h, 4 h, and 6 h after treatment. Our results demonstrate that significant differences in expression are only detected after 4 h of treatment. In addition, expression of TXNIP (an apoptotic inducer) remains elevated at 4 h and 6 h only in r-Moj-DN treated cells. Based on the consistency of elevated TXNIP expression, we further studied TXNIP as a novel target of disintegrin activation. Confocal microscopy of anti-TXNIP stained SK-Mel-28 cells suggests nuclear localization of TXNIP after r-Moj-DM treatment. A stable TXNIP knockdown SK-Mel-28 cell line was produced to test TXNIP' role in the apoptotic induction by r-Moj-DM. High cell viability (74.3% ±9.1) was obtained after r-Moj-DM treatment of TXNIP knocked down SK-Mel-28 cells, compared to 34% ±0.187 for untransduced cells. These results suggest that TXNIP is required early in the apoptotic-inducing pathway resulting from r-Moj-DM binding to the αv integrin subunit.

Trends in the Evolution of Snake Toxins Underscored by an Integrative Omics Approach to Profile the Venom of the Colubrid Phalotris mertensi

Only few studies on snake venoms were dedicated to deeply characterize the toxin secretion of animals from the Colubridae family, despite the fact that they represent the majority of snake diversity. As a consequence, some evolutionary trends observed in venom proteins that underpinned the evolutionary histories of snake toxins were based on data from a minor parcel of the clade. Here, we investigated the proteins of the totally unknown venom from Phalotris mertensi (Dipsadinae subfamily), in order to obtain a detailed profile of its toxins and to appreciate evolutionary tendencies occurring in colubrid venoms. By means of integrated omics and functional approaches, including RNAseq, Sanger sequencing, high-resolution proteomics, recombinant protein production, and enzymatic tests, we verified an active toxic secretion containing up to 21 types of proteins. A high content of Kunitz-type proteins and C-type lectins were observed, although several enzymatic components such as metalloproteinases and an L-amino acid oxidase were also present in the venom. Interestingly, an arguable venom component of other species was demonstrated as a true venom protein and named svLIPA (snake venom acid lipase). This finding indicates the importance of checking the actual protein occurrence across species before rejecting genes suggested to code for toxins, which are relevant for the discussion about the early evolution of reptile venoms. Moreover, trends in the evolution of some toxin classes, such as simplification of metalloproteinases and rearrangements of Kunitz and Wap domains, parallel similar phenomena observed in other venomous snake families and provide a broader picture of toxin evolution.

Colubrid Venom Composition: An -Omics Perspective

Snake venoms have been subjected to increasingly sensitive analyses for well over 100 years, but most research has been restricted to front-fanged snakes, which actually represent a relatively small proportion of extant species of advanced snakes. Because rear-fanged snakes are a diverse and distinct radiation of the advanced snakes, understanding venom composition among "colubrids" is critical to understanding the evolution of venom among snakes. Here we review the state of knowledge concerning rear-fanged snake venom composition, emphasizing those toxins for which protein or transcript sequences are available. We have also added new transcriptome-based data on venoms of three species of rear-fanged snakes. Based on this compilation, it is apparent that several components, including cysteine-rich secretory proteins (CRiSPs), C-type lectins (CTLs), CTLs-like proteins and snake venom metalloproteinases (SVMPs), are broadly distributed among "colubrid" venoms, while others, notably three-finger toxins (3FTxs), appear nearly restricted to the Colubridae (sensu stricto). Some putative new toxins, such as snake venom matrix metalloproteinases, are in fact present in several colubrid venoms, while others are only transcribed, at lower levels. This work provides insights into the evolution of these toxin classes, but because only a small number of species have been explored, generalizations are still rather limited. It is likely that new venom protein families await discovery, particularly among those species with highly specialized diets.

Full-Length Venom Protein cDNA Sequences from Venom-Derived mRNA: Exploring Compositional Variation and Adaptive Multigene Evolution

Envenomation of humans by snakes is a complex and continuously evolving medical emergency, and treatment is made that much more difficult by the diverse biochemical composition of many venoms. Venomous snakes and their venoms also provide models for the study of molecular evolutionary processes leading to adaptation and genotype-phenotype relationships. To compare venom complexity and protein sequences, venom gland transcriptomes are assembled, which usually requires the sacrifice of snakes for tissue. However, toxin transcripts are also present in venoms, offering the possibility of obtaining cDNA sequences directly from venom. This study provides evidence that unknown full-length venom protein transcripts can be obtained from the venoms of multiple species from all major venomous snake families. These unknown venom protein cDNAs are obtained by the use of primers designed from conserved signal peptide sequences within each venom protein superfamily. This technique was used to assemble a partial venom gland transcriptome for the Middle American Rattlesnake (Crotalus simus tzabcan) by amplifying sequences for phospholipases A₂, serine proteases, C-lectins, and metalloproteinases from within venom. Phospholipase A₂ sequences were also recovered from the venoms of several rattlesnakes and an elapid snake (Pseudechis porphyriacus), and three-finger toxin sequences were recovered from multiple rear-fanged snake species, demonstrating that the three major clades of advanced snakes (Elapidae, Viperidae, Colubridae) have stable mRNA present in their venoms. These cDNA sequences from venom were then used to explore potential activities derived from protein sequence similarities and evolutionary histories within these large multigene superfamilies. Venom-derived sequences can also be used to aid in characterizing venoms that lack proteomic profiles and identify sequence characteristics indicating specific envenomation profiles. This approach, requiring only venom, provides access to cDNA sequences in the absence of living specimens, even from commercial venom sources, to evaluate important regional differences in venom composition and to study snake venom protein evolution.

This work demonstrates that full-length venom protein messenger RNAs are present in secreted venoms and can be used to acquire full-length protein sequences of toxins from both front-fanged (Elapidae, Viperidae) and rear-fanged (Colubridae) snake venoms, eliminating the need to use venom glands. Full-length transcripts were obtained from venom samples that were fresh, newly lyophilized, old, field desiccated or commercially prepared, representing a significant advance over previous attempts which produced only partial sequence transcripts. Transcripts for all major venom protein families (metalloproteinases, serine proteases, C-type lectins, phospholipases A₂ and three-finger toxins) responsible for clinically significant snakebite symptoms were obtained from venoms. These sequences aid in the identification and characterization of venom proteome profiles, allowing for the identification of peptide sequences, specific isoforms, and novel venom proteins. The application of this technique will help to provide venom protein sequences for many snake species, including understudied rear-fanged snakes. Venom protein transcripts offer important insights into potential snakebite envenomation profiles and the molecular evolution of venom protein multigene families. By requiring only venom to obtain venom protein cDNAs, the approach detailed here will provide access to cDNA-based protein sequences from commercial and other venom sources, facilitating study of snake venom protein composition and evolution.

Stabilising the Integrity of Snake Venom mRNA Stored under Tropical Field Conditions Expands Research Horizons

BACKGROUND

Snake venoms contain many proteinaceous toxins that can cause severe pathology and mortality in snakebite victims. Interestingly, mRNA encoding such toxins can be recovered directly from venom, although yields are low and quality is unknown. It also remains unclear whether such RNA contains information about toxin isoforms and whether it is representative of mRNA recovered from conventional sources, such as the venom gland. Answering these questions will address the feasibility of using venom-derived RNA for future research relevant to biomedical and antivenom applications.

METHODOLOGY/PRINCIPAL FINDINGS

Venom was extracted from several species of snake, including both members of the Viperidae and Elapidae, and either lyophilized or immediately added to TRIzol reagent. TRIzol-treated venom was incubated at a range of temperatures (4-37°C) for a range of durations (0-48 hours), followed by subsequent RNA isolation and assessments of RNA quantity and quality. Subsequently, full-length toxin transcripts were targeted for PCR amplification and Sanger sequencing. TRIzol-treated venom yielded total RNA of greater quantity and quality than lyophilized venom, and with quality comparable to venom gland-derived RNA. Full-length sequences from multiple Viperidae and Elapidae toxin families were successfully PCR amplified from TRIzol-treated venom RNA. We demonstrated that venom can be stored in TRIzol for 48 hours at 4-19°C, and 8 hours at 37°C, at minimal cost to RNA quality, and found that venom RNA encoded multiple toxin isoforms that seemed homologous (98-99% identity) to those found in the venom gland.

CONCLUSIONS/SIGNIFICANCE

The non-invasive experimental modifications we propose will facilitate the future investigation of venom composition by using venom as an alternative source to venom gland tissue for RNA-based studies, thus obviating the undesirable need to sacrifice snakes for such research purposes. In addition, they expand research horizons to rare, endangered or protected snake species and provide more flexibility to performing fieldwork on venomous snakes in tropical conditions.

Functional characterizations of venom phenotypes in the eastern diamondback rattlesnake (Crotalus adamanteus) and evidence for expression-driven divergence in toxic activities among populations

Phenotypes frequently vary across and within species. The connection between specific phenotypic effects and function, however, is less understood despite being essential to our understanding of the adaptive process. Snake venoms are ideal for identifying functionally important phenotypic variation because venom variation is common, and venoms can be functionally characterized through simple assays and toxicity measurements. Previous work with the eastern diamondback rattlesnake (Crotalus adamanteus) used multivariate statistical approaches to identify six unique venom phenotypes. We functionally characterized hemolytic, gelatinase, fibrinogenolytic, and coagulant activity for all six phenotypes, as well as one additional venom, to determine if the statistically significant differences in toxin expression levels previously documented corresponded to differences in venom activity. In general, statistical differences in toxin expression predicted the identified functional differences, or lack thereof, in toxic activity, demonstrating that the statistical approach used to characterize C. adamanteus venoms was a fair representation of biologically meaningful differences. Minor differences in activity not accounted for by the statistical model may be the result of amino-acid differences and/or post-translational modifications, but overall we were able to link variation in protein expression levels to variation in function as predicted by multivariate statistical approaches.