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

Structural studies with crotoxin B from Crotalus durissus collilineatus venom suggest a heterodimeric assembly formed by two new isoforms

In accidents involving Crotalus snakes, the crotoxin complex (CTX) plays lethal action due to its neurotoxic activity. On the other hand, CTX have potential biotechnological application due to its anti-tumoral, anti-inflammatory, antimicrobial, analgesic and immunomodulatory properties. CTX is a heterodimer composed of Crotoxin A (CA or crotapotin), the acidic nontoxic and non-enzymatic component and; Crotoxin B (CB), a basic, toxic and catalytic PLA2. Currently, there are two classes of CTX isoforms, whose differences in their biological activities have been attributed to features presented in CB isoforms. Here, we present the crystal structure of CB isolated from the Crotalus durissus collilineatus venom. It amino acid sequence was assigned using the SEQUENCE SLIDER software, which revealed that the crystal structure is a heterodimer composed of two new CB isoforms (colCB-A and colCB-B). Bioinformatic and biophysical analyses showed that the toxin forms a tetrameric assembly in solution similar to CB from Crotalus durissus terrificus venom, despite some differences observed at the dimeric interface. By the previously proposed classification, the colCB-B presents features of the class I isoforms while colCB-A cannot be classified into classes I and II based on its amino acid sequence. Due to similar features observed for other CB isoforms found in the NCBI database and the results obtained for colCB-A, we suggest that there are more than two classes of CTX and CB isoforms in crotalic venoms.

A Review of Rattlesnake Venoms

Venom components are invaluable in biomedical research owing to their specificity and potency. Many of these components exist in two genera of rattlesnakes, Crotalus and Sistrurus, with high toxicity and proteolytic activity variation. This review focuses on venom components within rattlesnakes, and offers a comparison and itemized list of factors dictating venom composition, as well as presenting their known characteristics, activities, and significant applications in biosciences. There are 64 families and subfamilies of proteins present in Crotalus and Sistrurus venom. Snake venom serine proteases (SVSP), snake venom metalloproteases (SVMP), and phospholipases A2 (PLA2) are the standard components in Crotalus and Sistrurus venom. Through this review, we highlight gaps in the knowledge of rattlesnake venom; there needs to be more information on the venom composition of three Crotalus species and one Sistrurus subspecies. We discuss the activity and importance of both major and minor components in biomedical research and drug development.

Comparative analysis of PacBio and ONT RNA sequencing methods for Nemopilema Nomurai venom identification

Recent studies on marine organisms have made use of third-generation sequencing technologies such as Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT). While these specialized bioinformatics tools have different algorithmic designs and performance capabilities, they offer scalability and can be applied to various datasets. We investigated the effectiveness of PacBio and ONT RNA sequencing methods in identifying the venom of the jellyfish species Nemopilema nomurai. We conducted a detailed analysis of the sequencing data from both methods, focusing on key characteristics such as CD, alternative splicing, long-chain noncoding RNA, simple sequence repeat, transcription factor, and functional transcript annotation. Our findings indicate that ONT generally produced higher raw data quality in the transcriptome analysis, while PacBio generated longer read lengths. PacBio was found to be superior in identifying CDs and long-chain noncoding RNA, whereas ONT was more cost-effective for predicting alternative splicing events, simple sequence repeats, and transcription factors. Based on these results, we conclude that PacBio is the most specific and sensitive method for identifying venom components, while ONT is the most cost-effective method for studying venogenesis, cnidocyst (venom gland) development, and transcription of virulence genes in jellyfish. Our study has implications for future sequencing technologies in marine jellyfish, and highlights the power of full-length transcriptome analysis in discovering potential therapeutic targets for jellyfish dermatitis.

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.

Comparison of Four Methods of RNA Extraction and cDNA Synthesis from The Venom of Peruvian Snakes of the Genus Bothrops of Clinical Importance

RNA purification and cDNA synthesis represents the starting point for molecular analyses of snake venom proteins-enzymes. Usually, the sacrifice of snakes is necessary for venom gland extraction to identify protein-coding transcripts; however, the venom can be used as a source of transcripts. Although there are methods for obtaining RNA from venom, no comparative analysis has been conducted in the Bothrops genus. In the present study, we compared four commercial methods for RNA purification and cDNA synthesis from venom (liquid, lyophilized, or long-term storage) of four clinically relevant species of Peruvian Bothrops. Our results show that the TRIzol method presents the highest yield of RNA purified from venom (59 ± 11 ng/100 µL or 10 mg). The SuperScript First-Strand Synthesis System kit produced high amounts of cDNA (3.2 ± 1.2 ng cDNA/ng RNA), and the highest value was from combination with the Dynabeads mRNA DIRECT kit (4.8 ± 2.0 ng cDNA/ng RNA). The utility of cDNA was demonstrated with the amplification of six relevant toxins: thrombin-like enzymes, P-I and P-III metalloproteinases, acid and basic phospholipases A2, and disintegrins. To our knowledge, this is the first comparative study of RNA purification and cDNA synthesis methodologies from Bothrops genus venom.

Venom Proteomics of Trimeresurus gracilis, a Taiwan-Endemic Pitviper, and Comparison of Its Venom Proteome and VEGF and CRISP Sequences with Those of the Most Related Species

Trimeresurus gracilis is an endemic alpine pitviper in Taiwan with controversial phylogeny, and its venom proteome remains unknown. In this study, we conducted a proteomic analysis of T. gracilis venom using high-performance liquid chromatography-tandem mass spectrometry and identified 155 toxin proteoforms that belong to 13 viperid venom toxin families. By searching the sequences of trypsin-digested peptides of the separated HPLC fractions against the NCBI database, T. gracilis venom was found to contain 40.3% metalloproteases (SVMPs), 15.3% serine proteases, 6.6% phospholipases A2, 5.0% L-amino acid oxidase, 4.6% Cys-rich secretory proteins (CRISPs), 3.2% disintegrins, 2.9% vascular endothelial growth factors (VEGFs), 1.9% C-type lectin-like proteins, and 20.2% of minor toxins, nontoxins, and unidentified peptides or compounds. Sixteen of these proteoforms matched the toxins whose full amino-acid sequences have been deduced from T. gracilis venom gland cDNA sequences. The hemorrhagic venom of T. gracilis appears to be especially rich in PI-class SVMPs and lacks basic phospholipase A2. We also cloned and sequenced the cDNAs encoding two CRISP and three VEGF variants from T. gracilis venom glands. Sequence alignments and comparison revealed that the PI-SVMP, kallikrein-like proteases, CRISPs, and VEGF-F of T. gracilis and Ovophis okinavensis are structurally most similar, consistent with their close phylogenetic relationship. However, the expression levels of some of their toxins were rather different, possibly due to their distinct ecological and prey conditions.

Extracellular Vesicles from Bothrops jararaca Venom Are Diverse in Structure and Protein Composition and Interact with Mammalian Cells

Snake venoms are complex cocktails of non-toxic and toxic molecules that work synergistically for the envenoming outcome. Alongside the immediate consequences, chronic manifestations and long-term sequelae can occur. Recently, extracellular vesicles (EVs) were found in snake venom. EVs mediate cellular communication through long distances, delivering proteins and nucleic acids that modulate the recipient cell's function. However, the biological roles of snake venom EVs, including possible cross-organism communication, are still unknown. This knowledge may expand the understanding of envenoming mechanisms. In the present study, we isolated and characterized the EVs from Bothrops jararaca venom (Bj-EVs), giving insights into their biological roles. Fresh venom was submitted to differential centrifugation, resulting in two EV populations with typical morphology and size range. Several conserved EV markers and a subset of venom related EV markers, represented mainly by processing enzymes, were identified by proteomic analysis. The most abundant protein family observed in Bj-EVs was 5'-nucleotidase, known to be immunosuppressive and a low abundant and ubiquitous toxin in snake venoms. Additionally, we demonstrated that mammalian cells efficiently internalize Bj-EVs. The commercial antibothropic antivenom partially recognizes Bj-EVs and inhibits cellular EV uptake. Based on the proteomic results and the in vitro interaction assays using macrophages and muscle cells, we propose that Bj-EVs may be involved not only in venom production and processing but also in host immune modulation and long-term effects of envenoming.

Evidence that genomic incompatibilities and other multilocus processes impact hybrid fitness in a rattlesnake hybrid zone

Hybrid zones provide valuable opportunities to understand the genomic mechanisms that promote speciation by providing insight into factors involved in intermediate stages of speciation. Here, we investigate introgression in a hybrid zone between two rattlesnake species (Crotalus viridis and Crotalus oreganus concolor) that have undergone historical allopatric divergence and recent range expansion and secondary contact. We use Bayesian genomic cline models to characterize genomic patterns of introgression between these lineages and identify loci potentially subject to selection in hybrids. We find evidence for a large number of genomic regions with biased ancestry that deviate from the genomic background in hybrids (i.e., excess ancestry loci), which tend to be associated with genomic regions with higher recombination rates. We also identify suites of excess ancestry loci that show highly correlated allele frequencies (including conspecific and heterospecific combinations) across physically unlinked genomic regions in hybrids. Our findings provide evidence for multiple multilocus evolutionary processes impacting hybrid fitness in this system.

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.

Functional, immunological characterization, and anticancer activity of BaMtx: A new Lys49- PLA2 homologue isolated from the venom of Peruvian Bothrops atrox snake (Serpentes: Viperidae)

Bothorps atrox is responsible for most of the ophidism cases in Perú. As part of the envenoming, myotoxicity is one of the most recurrent and destructive effects. In this study, a myotoxin, named BaMtx, was purified from B. atrox venom to elucidate its biological, immunological, and molecular characteristics. BaMtx was purified using CM-Sephadex-C-25 ion-exchange resin and SDS-PAGE analysis showed a unique protein band of 13 kDa or 24 kDa under reducing or non-reducing conditions, respectively. cDNA sequence codified a 122-aa mature protein with high homology with other Lys49-PLA₂s; modeled structure showed a N-terminal helix, a β-wing region, and a C-terminal random coil. This protein has a poor phospholipase A2 enzymatic activity. BaMtx has myotoxic (DMM = 12.30 ± 0.95 μg) and edema-forming (DEM = 26.00 ± 1.15 μg) activities. Rabbit immunization with purified enzyme produced anti-BaMtx antibodies that reduced 50.28 ± 10.15% of myotoxic activity and showed significant cross-reactivity against B. brazili and B pictus venoms. On the other hand, BaMtx exhibits mild anti-proliferative and anti-migratory effects on breast cancer cells, affecting the ROS and NADH levels, which may reduce mitochondrial respiration. These results contribute to the understanding of B. atrox Lys49-PLA2 effects and establish the anticancer potential de BaMtx.

A non-lethal method for studying scorpion venom gland transcriptomes, with a review of potentially suitable taxa to which it can be applied

Scorpion venoms are mixtures of proteins, peptides and small molecular compounds with high specificity for ion channels and are therefore considered to be promising candidates in the venoms-to-drugs pipeline. Transcriptomes are important tools for studying the composition and expression of scorpion venom. Unfortunately, studying the venom gland transcriptome traditionally requires sacrificing the animal and therefore is always a single snapshot in time. This paper describes a new way of generating a scorpion venom gland transcriptome without sacrificing the animal, thereby allowing the study of the transcriptome at various time points within a single individual. By comparing these venom-derived transcriptomes to the traditional whole-telson transcriptomes we show that the relative expression levels of the major toxin classes are similar. We further performed a multi-day extraction using our proposed method to show the possibility of doing a multiple time point transcriptome analysis. This allows for the study of patterns of toxin gene activation over time a single individual, and allows assessment of the effects of diet, season and other factors that are known or likely to influence intraindividual venom composition. We discuss the gland characteristics that may allow this method to be successful in scorpions and provide a review of other venomous taxa to which this method may potentially be successfully applied.

Intraspecific venom variation of Mexican West Coast Rattlesnakes (Crotalus basiliscus) and its implications for antivenom production

Intraspecific variation in snake venoms has been widely documented worldwide. However, there are few studies on this subject in Mexico. Venom characterization studies provide important data used to predict clinical syndromes, to evaluate the efficacy of antivenoms and, in some cases, to improve immunogenic mixtures in the production of antivenoms. In the present work, we evaluated the intraspecific venom variation of Crotalus basiliscus, a rattlesnake of medical importance and whose venom is used in the immunization of horses to produce one of the Mexican antivenoms. Our results demonstrate that there is variation in biological and biochemical activities among adult venoms and that there is an ontogenetic change from juvenile to adult venoms. Juvenile venoms were more lethal and had higher percentages of crotamine and crotoxin, while adult venoms had higher percentages of snake venom metalloproteases (SVMPs). Additionally, we documented crotoxin-like PLA₂ variation in which specimens from Zacatecas, Sinaloa and Michoacán (except 1) lacked the neurotoxin, while the rest of the venoms had it. Finally, we evaluated the efficacy of three lots of Birmex antivenom and all three were able to neutralize the lethality of four representative venoms but were not able to neutralize crotamine. We also observed significant differences in the LD₅₀ values neutralized per vial among the different lots. Based on these results, we recommend including venoms containing crotamine in the production of antivenom for a better immunogenic mixture and to improve the homogeneity of lots.

A Meta-Analysis of the Protein Components in Rattlesnake Venom

The specificity and potency of venom components give them a unique advantage in developing various pharmaceutical drugs. Though venom is a cocktail of proteins, rarely are the synergy and association between various venom components studied. Understanding the relationship between various components of venom is critical in medical research. Using meta-analysis, we observed underlying patterns and associations in the appearance of the toxin families. For Crotalus, Dis has the most associations with the following toxins: PDE; BPP; CRL; CRiSP; LAAO; SVMP P-I and LAAO; SVMP P-III and LAAO. In Sistrurus venom, CTL and NGF have the most associations. These associations can predict the presence of proteins in novel venom and understand synergies between venom components for enhanced bioactivity. Using this approach, the need to revisit the classification of proteins as major components or minor components is highlighted. The revised classification of venom components is based on ubiquity, bioactivity, the number of associations, and synergies. The revised classification can be expected to trigger increased research on venom components, such as NGF, which have high biomedical significance. Using hierarchical clustering, we observed that the genera's venom compositions were similar, based on functional characteristics rather than phylogenetic relationships.

Rapid Identification of Phospholipase A₂ Transcripts from Snake Venoms

Phospholipase A₂ (PLA₂) is a major component in snake venoms and it is found in many different isoforms. To identify transcripts encoding different PLA₂ isoforms, we developed a simple, rapid procedure. Total RNA was extracted from the venoms of three cottonmouth snakes and two diamondback rattlesnakes, and further reverse-transcribed into complementary DNA (cDNA). Using one pair of cottonmouth PLA₂-specific primers and a Reverse Transcription Polymerase Chain Reaction (RT-PCR) technique, we identified 27 unique full-length PLA₂ transcripts, including nine sequences identical to the previously documented ones in the database and one novel GIII-like PLA₂. Two common transcripts respectively encoding Asp49 and Lys49 PLA₂ isoforms were identified in all three cottonmouth venoms, that contain more PLA₂ transcripts than the diamondback rattlesnake venoms. The placement of cloned PLA₂ transcripts in snake venom PLA₂s was further discussed by phylogenetic analysis. The procedure developed in this study paves the way for accelerated acquisition of transcriptome data on any other venom toxin families. The results obtained are crucial for insight into the structure and function of PLA₂ isoforms for scientific and potential therapeutic purposes.

DNA barcodes from snake venom: a broadly applicable method for extraction of DNA from snake venoms

DNA barcoding is a simple technique used to develop a large-scale system of classification that is broadly applicable across a wide variety of taxa. DNA-based analysis of snake venoms can provide a system of classification independent of currently accepted taxonomic relationships by generating DNA barcodes specific to each venom sample. DNA purification from dried snake venoms has previously required large amounts of starting material, has resulted in low yields and inconsistent amplification, and was possible with front-fanged snakes only. Here, we present a modified DNA extraction protocol applied to venoms of both front- and rear-fanged snakes that requires significantly less starting material (1 mg) and yields sufficient amounts of DNA for successful PCR amplification of regions commonly used for DNA barcoding. [Formula: see text].

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.

Engineering varied serine protease inhibitors by converting P1 site of BF9, a weakly active Kunitz-type animal toxin

Although there were a lot of weakly active animal toxins in the venoms, their values and applications are still mysterious, such as BF9, which is a Kunitz-type toxin isolated from the venom of the elapid snake Bungarus fasciatus. Here, we used BF9 to be a molecular scaffold, and engineered eight BF9-derived peptides by changing P1 site Asn17 of BF9, such as BF9-N17Y and BF9-N17T designed from the polar subfamily, BF9-N17L and BF9-N17G designed from the Non-polar subfamily, BF9-N17D designed from acidic subfamily, and BF9-N17H, BF9-N17K and BF9-N17R designed from basic subfamily. Through enzyme inhibitor experiment assays, we found a potent and selective chymotrypsin inhibitor BF9-N17Y, a potent and selective coagulation factor XIa inhibitor BF9-N17H, and two highly potent coagulation factor XIa inhibitors BF9-N17K and BF9-N17. APTT and PT assays further showed that BF9-N17H, BF9-N17K and BF9-N17R were three novel anticoagulants with selectively intrinsic coagulation pathway inhibitory activity. Considering that natural weakly active animal toxins are also a huge peptide resource, our present work might open a new window about pharmacological applications of weakly active animal toxins, which might be good templates for potent and selective molecular probe and lead drug designs.

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.

Modern trends in animal venom research - omics and nanomaterials

Animal venom research is a specialized investigation field, in which a number of different methods are used and this array is constantly expanding. Thus, recently emerged omics and nanotechnologies have already been successfully applied to venom research. Animal venoms have been studied for quite a long time. The traditional reductionist approach has been to isolate individual toxins and then study their structure and function. Unfortunately, the characterization of the venom as a whole system and its multiple effects on an entire organism were not possible until recent times. The development of new methods in mass spectrometry and sequencing have allowed such characterizations of venom, encompassing the identification of new toxins present in venoms at extremely low concentrations to changes in metabolism of prey organisms after envenomation. In particular, this type of comprehensive research has become possible due to the development of the various omics technologies: Proteomics, peptidomics, transcriptomics, genomics and metabolomics. As in other research fields, these omics technologies ushered in a revolution for venom studies, which is now entering the era of big data. Nanotechnology is a very new branch of technology and developing at an extremely rapid pace. It has found application in many spheres and has not bypassed the venom studies. Nanomaterials are quite promising in medicine, and most studies combining venoms and nanomaterials are dedicated to medical applications. Conjugates of nanoparticles with venom components have been proposed for use as drugs or diagnostics. For example, nanoparticles conjugated with chlorotoxin - a toxin in scorpion venom, which has been shown to bind specifically to glioma cells - are considered as potential glioma-targeted drugs, and conjugates of neurotoxins with fluorescent semiconductor nanoparticles or quantum dots may be used to detect endogenous targets expressed in live cells. The data on application of omics and nanotechnologies in venom research are systematized concisely in this paper.

Understanding Biological Roles of Venoms Among the Caenophidia: The Importance of Rear-Fanged Snakes

Snake venoms represent an adaptive trophic response to the challenges confronting a limbless predator for overcoming combative prey, and this chemical means of subduing prey shows several dominant phenotypes. Many front-fanged snakes, particularly vipers, feed on various vertebrate and invertebrate prey species, and some of their venom components (e.g., metalloproteinases, cobratoxin) appear to have been selected for "broad-brush" incapacitation of different prey taxa. Using proteomic and genomic techniques, the compositional diversity of front-fanged snakes is becoming well characterized; however, this is not the case for most rear-fanged colubroid snakes. Because these species consume a high diversity of prey, and because venoms are primarily a trophic adaptation, important clues for understanding specific selective pressures favoring venom component composition will be found among rear-fanged snake venoms. Rear-fanged snakes typically (but not always) produce venoms with lower complexity than front-fanged snakes, and there are even fewer dominant (and, arguably, biologically most relevant) venom protein families. We have demonstrated taxon-specific toxic effects, where lizards and birds show high susceptibility while mammals are largely unaffected, for both Old World and New World rear-fanged snakes, strongly indicating a causal link between toxin evolution and prey preference. New data are presented on myotoxin a, showing that the extremely rapid paralysis induced by this rattlesnake toxin is specific for rodents, and that myotoxin a is ineffectual against lizards. Relatively few rear-fanged snake venoms have been characterized, and basic natural history data are largely lacking, but directed sampling of specialized species indicates that novel compounds are likely among these specialists, particularly among those species feeding on invertebrate prey such as scorpions and centipedes. Because many of the more than 2200 species of colubroid snakes are rear-fanged, and many possess a Duvernoy's venom gland, understanding the nature of their venoms is foundational to understanding venom evolution in advanced snakes.

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.