Complementary DNA sequencing and identification of mRNAs from the venomous gland of Agkistrodon piscivorus leucostoma

Bibliometric Analysis of Literature in Snake Venom-Related Research Worldwide (1933-2022)

Snake envenomation is a severe economic and health concern affecting countries worldwide. Snake venom carries a wide variety of small peptides and proteins with various immunological and pharmacological properties. A few key research areas related to snake venom, including its applications in treating cancer and eradicating antibiotic-resistant bacteria, have been gaining significant attention in recent years. The goal of the current study was to analyze the global profile of literature in snake venom research. This study presents a bibliometric review of snake venom-related research documents indexed in the Scopus database between 1933 and 2022. The overall number of documents published on a global scale was 2999, with an average annual production of 34 documents. Brazil produced the highest number of documents (n = 729), followed by the United States (n = 548), Australia (n = 240), and Costa Rica (n = 235). Since 1963, the number of publications has been steadily increasing globally. At a worldwide level, antivenom, proteomics, and transcriptomics are growing hot issues for research in this field. The current research provides a unique overview of snake venom research at global level from 1933 through 2022, and it may be beneficial in guiding future research.

Structural, enzymatic and pharmacological profiles of AplTX-II - A basic sPLA2 (D49) isolated from the Agkistrodon piscivorus leucostoma snake venom

A basic sPLA₂ (D49) from the venom of snake Agkistrodon piscivorus leucostoma (AplTX-II) was isolated, purified and characterized. We determined the enzymatic and pharmacological profiles of this toxin. AplTX-II was isolated with a high level of purity through reverse phase chromatography and molecular exclusion. The enzyme showed pI 9.48 and molecular weight of 14,003 Da. The enzymatic activity of the AplTX-II depended on Ca²⁺ pH and temperature. The comparison of the primary structure with other sPLA₂s revealed that AplTX-II presented all the structural reasons expected for a basic sPLA₂s. Additionally, we have resolved its structure with the docked synthetic substrate NOBA (4-nitro-3-octanoyloxy benzoic acid) by homology modeling, and performed MD simulations with explicit solvent. Structural similarities were found between the enzyme's modeled structure and other snake sPLA₂ X-Ray structures, available in the PDB database. NOBA and active-site water molecules spontaneously adopted stable positions and established interactions in full agreement with the reaction mechanism, proposed for the physiological substrate, suggesting that NOBA hydrolysis is an excellent model to study phospholipid hydrolysis.

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.

Transcriptomics and venomics: implications for medicinal chemistry

Over the last three decades, transcriptomic studies of venom gland cells have continuously evolved, opening up new possibilities for exploring the molecular diversity of animal venoms, a prerequisite for the discovery of new drug candidates and molecular phylogenetics. The molecular complexity of animal venoms is much greater than initially thought. In this review, we describe the different technologies available for transcriptomic studies of venom, from the original individual cloning approaches to the more recent global Next Generation Sequencing strategies. Our understanding of animal venoms is evolving, with the discovery of complex and diverse bio-optimized cocktails of compounds, including mostly peptides and proteins, which are now beginning to be studied by academic and industrial researchers.

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

Background

Largely because of their direct, negative impacts on human health, the venoms of front-fanged snakes of the families Viperidae and Elapidae have been extensively characterized proteomically, transcriptomically, and pharmacologically. However, relatively little is known about the molecular complexity and evolution of the venoms of rear-fanged colubrid snakes, which are, with a few notable exceptions, regarded as harmless to humans. Many of these snakes have venoms with major effects on their preferred prey, and their venoms are probably as critical to their survival as those of front-fanged elapids and viperids.

Results

We sequenced the venom-gland transcriptomes from a specimen of Hypsiglena (Desert Night Snake; family Colubridae, subfamily Dipsadinae) and of Boiga irregularis (Brown Treesnake; family Colubridae, subfamily Colubrinae) and verified the transcriptomic results proteomically by means of high-definition mass spectrometry. We identified nearly 3,000 nontoxin genes for each species. For B. irregularis, we found 108 putative toxin transcripts in 46 clusters with <1% nucleotide divergence, and for Hypsiglena we identified 79 toxin sequences that were grouped into 33 clusters. Comparisons of the venoms revealed divergent venom types, with Hypsiglena possessing a viper-like venom dominated by metalloproteinases, and B. irregularis having a more elapid-like venom, consisting primarily of three-finger toxins.

Conclusions

Despite the difficulty of procuring venom from rear-fanged species, we were able to complete all analyses from a single specimen of each species without pooling venom samples or glands, demonstrating the power of high-definition transcriptomic and proteomic approaches. We found a high level of divergence in the venom types of two colubrids. These two venoms reflected the hemorrhagic/neurotoxic venom dichotomy that broadly characterizes the difference in venom strategies between elapids and viperids.

Venom gland transcriptomics for identifying, cataloging, and characterizing venom proteins in snakes

Snake venoms are cocktails of protein toxins that play important roles in capture and digestion of prey. Significant qualitative and quantitative variation in snake venom composition has been observed among and within species. Understanding these variations in protein components is instrumental in interpreting clinical symptoms during human envenomation and in searching for novel venom proteins with potential therapeutic applications. In the last decade, transcriptomic analyses of venom glands have helped in understanding the composition of various snake venoms in great detail. Here we review transcriptomic analysis as a powerful tool for understanding venom profile, variation and evolution.

Genetic Basis for Variation of Metalloproteinase-Associated Biochemical Activity in Venom of the Mojave Rattlesnake (Crotalus scutulatus scutulatus)

The metalloproteinase composition and biochemical profiles of rattlesnake venom can be highly variable among rattlesnakes of the same species. We have previously shown that the neurotoxic properties of the Mojave rattlesnake (Crotalus scutulatus scutulatus) are associated with the presence of the Mojave toxin A subunit suggesting the existence of a genetic basis for rattlesnake venom composition. In this report, we hypothesized the existence of a genetic basis for intraspecies variation in metalloproteinase-associated biochemical properties of rattlesnake venom of the Mojave rattlesnake. To address this question, we PCR-amplified and compared the genomic DNA nucleotide sequences that code for the mature metalloproteinase domain of fourteen Mojave rattlesnakes captured from different geographical locations across the southwest region of the United States. In addition, the venoms from the same rattlesnakes were tested for their ability to hydrolyze fibrinogen, fibrin, casein, and hide powder azure and for induction of hemorrhage in mice. Overall, based on genomic sequencing and biochemical data, we classified Mojave rattlesnake venom into four distinct groups of metalloproteinases. These findings indicate that differences in nucleotide sequences encoding the mature proteinase domain and noncoding regions contribute to differences in venom metalloproteinase activities among rattlesnakes of the same species.

Cloning of cDNAs and molecular characterisation of C-type lectin-like proteins from snake venoms

C-type lectin-like proteins (CTLPs) isolated from snake venoms are the largest and most complex non-mammalian vertebrate C-type lectin-like domain family. In the present study, we simultaneously amplified four cDNAs encoding different types of CTLP subunits from the venoms of two different species of snakes by RT-PCR with a single sense primer and a nested universal primer - two CTLP subunit-encoding cDNAs were cloned from Deinagkistrodon acutus venom and two from Agkistrodon halys Pallas venom. All four cloned CTLP subunits exhibited typical motifs in their corresponding domain regions but with relatively-low sequence similarities to each other. Compared with previously-published CTLPs, the four cloned CTLPs subunits showed slight variations in the calcium-binding sites and the disulphide bonding patterns. To our knowledge, these data constitute the first example of co-expression of CTLP platelet glycoprotein Ib-binding subunits and coagulation factors in Agkistrodon halys Pallas venom.

A novel method for the purification of low soluble recombinant C-type lectin proteins

Snake venoms contain a complex mixture of many biological molecules including proteins. The purification of recombinant proteins is a key step in studying their function and structure with affinity chromatography as the common method used in their purification. In bacterial expression systems, hydrophobic recombinant proteins are usually precipitated into inclusion bodies, and contaminants are typically associated with tagged proteins after purification. The purpose of this study was to develop a procedure to purify hydrophobic recombinant proteins without an affinity tag. Snake venom mature C-type lectin-like proteins (CLPs) with a tag were cloned, expressed, and purified by repeated sonication and wash steps. The effects of the signal peptide on the expression and solubility of the recombinant protein were investigated. The CLPs in washed inclusion bodies were solubilized and refolded by dialysis. The CLPs without a tag were successfully purified with a yield 38 times higher than the traditional method, and inhibited blood platelet aggregation with an IC(50) of 100.57 μM in whole blood. This novel procedure is a rapid, and inexpensive method to purify functional recombinant hydrophobic CLPs from snake venoms useful in the development of drug therapies.

Phylogenetic analysis of serine proteases from Russell's viper (Daboia russelli siamensis) and Agkistrodon piscivorus leucostoma venom

Serine proteases are widely found in snake venoms. They have variety of functions including contributions to hemostasis. In this study, five serine proteases were cloned and characterized from two different cDNA libraries: factor V activator (RVV-V), alpha fibrinogenase (RVAF) and beta fibrinogenase (RVBF) from Russell's viper (Daboia russelli siamensis), and plasminogen activator (APL-PA) and protein C activator (APL-C) from Agkistrodon piscivorus leucostoma. The snake venom serine proteases were clustered in phylogenetic tree according to their functions. K(A)/K(S) values suggested that accelerated evolution has occurred in the mature protein coding regions in cDNAs of snake venom serine proteases.

A high-throughput venom-gland transcriptome for the Eastern Diamondback Rattlesnake (Crotalus adamanteus) and evidence for pervasive positive selection across toxin classes

Despite causing considerable human mortality and morbidity, animal toxins represent a valuable source of pharmacologically active macromolecules, a unique system for studying molecular adaptation, and a powerful framework for examining structure-function relationships in proteins. Snake venoms are particularly useful in the latter regard as they consist primarily of a moderate number of proteins and peptides that have been found to belong to just a handful of protein families. As these proteins and peptides are produced in dedicated glands, transcriptome sequencing has proven to be an effective approach to identifying the expressed toxin genes. We generated a venom-gland transcriptome for the Eastern Diamondback Rattlesnake (Crotalus adamanteus) using Roche 454 sequencing technology. In the current work, we focus on transcripts encoding toxins. We identified 40 unique toxin transcripts, 30 of which have full-length coding sequences, and 10 have only partial coding sequences. These toxins account for 24% of the total sequencing reads. We found toxins from 11 previously described families of snake-venom toxins and have discovered two putative, previously undescribed toxin classes. The most diverse and highly expressed toxin classes in the C. adamanteus venom-gland transcriptome are the serine proteinases, metalloproteinases, and C-type lectins. The serine proteinases are the most abundant class, accounting for 35% of the toxin sequencing reads. Metalloproteinases are the most diverse; 11 different forms have been identified. Using our sequences and those available in public databases, we detected positive selection in seven of the eight toxin families for which sufficient sequences were available for the analysis. We find that the vast majority of the genes that contribute directly to this vertebrate trait show evidence for a role for positive selection in their evolutionary history.

A transcriptomic analysis of gene expression in the venom gland of the snake Bothrops alternatus (urutu)

Background

The genus Bothrops is widespread throughout Central and South America and is the principal cause of snakebite in these regions. Transcriptomic and proteomic studies have examined the venom composition of several species in this genus, but many others remain to be studied. In this work, we used a transcriptomic approach to examine the venom gland genes of Bothrops alternatus, a clinically important species found in southeastern and southern Brazil, Uruguay, northern Argentina and eastern Paraguay.

Results

A cDNA library of 5,350 expressed sequence tags (ESTs) was produced and assembled into 838 contigs and 4512 singletons. BLAST searches of relevant databases showed 30% hits and 70% no-hits, with toxin-related transcripts accounting for 23% and 78% of the total transcripts and hits, respectively. Gene ontology analysis identified non-toxin genes related to general metabolism, transcription and translation, processing and sorting, (polypeptide) degradation, structural functions and cell regulation. The major groups of toxin transcripts identified were metalloproteinases (81%), bradykinin-potentiating peptides/C-type natriuretic peptides (8.8%), phospholipases A₂ (5.6%), serine proteinases (1.9%) and C-type lectins (1.5%). Metalloproteinases were almost exclusively type PIII proteins, with few type PII and no type PI proteins. Phospholipases A₂ were essentially acidic; no basic PLA₂ were detected. Minor toxin transcripts were related to L-amino acid oxidase, cysteine-rich secretory proteins, dipeptidylpeptidase IV, hyaluronidase, three-finger toxins and ohanin. Two non-toxic proteins, thioredoxin and double-specificity phosphatase Dusp6, showed high sequence identity to similar proteins from other snakes. In addition to the above features, single-nucleotide polymorphisms, microsatellites, transposable elements and inverted repeats that could contribute to toxin diversity were observed.

Conclusions

Bothrops alternatus venom gland contains the major toxin classes described for other Bothrops venoms based on trancriptomic and proteomic studies. The predominance of type PIII metalloproteinases agrees with the well-known hemorrhagic activity of this venom, whereas the lower content of serine proteases and C-type lectins could contribute to less marked coagulopathy following envenoming by this species. The lack of basic PLA₂ agrees with the lower myotoxicity of this venom compared to other Bothrops species with these toxins. Together, these results contribute to our understanding of the physiopathology of envenoming by this species.

Genome evolution in Reptilia, the sister group of mammals

The genomes of birds and nonavian reptiles (Reptilia) are critical for understanding genome evolution in mammals and amniotes generally. Despite decades of study at the chromosomal and single-gene levels, and the evidence for great diversity in genome size, karyotype, and sex chromosome diversity, reptile genomes are virtually unknown in the comparative genomics era. The recent sequencing of the chicken and zebra finch genomes, in conjunction with genome scans and the online publication of the Anolis lizard genome, has begun to clarify the events leading from an ancestral amniote genome--predicted to be large and to possess a diverse repeat landscape on par with mammals and a birdlike sex chromosome system--to the small and highly streamlined genomes of birds. Reptilia exhibit a wide range of evolutionary rates of different subgenomes and, from isochores to mitochondrial DNA, provide a critical contrast to the genomic paradigms established in mammals.

Molecular cloning and characterization of cDNAs encoding metalloproteinases from snake venom glands

Snake venom metalloproteinases (SVMPs) are a superfamily of zinc-dependent proteases and participate in a number of important biological, physiological and pathophysiological processes. In this work, we simultaneously amplified nine cDNAs encoding different classes of metalloproteinases from glands of four different snake species (Agkistrodon contortrix laticinctus, Crotalus atrox, Crotalus viridis viridis and Agkistrodon piscivorus leucostoma) by RT-PCR with a pair of primers. Among the encoded metalloproteinases, two enzymes (AclVMP-I and AplVMP-I), three enzymes (CaVMP-II, CvvVMP-II and AplVMP-II) and four enzymes (AclVMP-III, CaVMP-III, CvvVMP-III and AplVMP-III) with the characteristic motif (HEXXHXXGXXH) of metalloproteinase belong to type P-I, P-II and P-III enzymes, respectively. Disintegrin domains of CaVMP-II and CvvVMP-II from two Crotatus snakes contain RGD-motif whereas AplVMP-II from Agkistrodon snake has KGD-motif. Instead of R/KGD-motif within disintegrin domain of SVMP-II enzyme, CaVMP-III, CvvVMP-III and AplVMP-III enzymes contain SECD-motif, while AclVMP-III has DDCD-motif in their corresponding position of disintegrin-like domains. There are 12 Cys amino acids in cysterin-rich domains of each P-III enzyme. Moreover, a disintegrin precursor (AplDis) with RGD-motif also simultaneously amplified from the glands of A.p. leucostoma while amplifying AplVMP-II and AplVMP-III, which indicated that different types of SVMPs and related genes are present in a single species of snake and share a consensus sequence at the 3' and 5' untranslated regions. RT-PCR result also showed that P-III is highly expressed in Crotalus snakes than in Agkistrodon snakes. Aligning the deduced amino acid sequence of these enzymes with other SVMPs from GenBank database indicated that this is the first report on the isolation of cDNAs encoding P-II and P-III enzymes from C.v. viridis and A.p. leucostoma snakes. The availability of these SVMP sequences directly facilitated further studies of structure characterization and diversified function analysis.

cDNA cloning, expression and fibrin(ogen)olytic activity of two low-molecular weight snake venom metalloproteinases

Two cDNA clones, AplVMP1 and AplVMP2, were isolated from a snake (Agkistrodon piscivorus leucostoma) venom gland cDNA library. The full-length cDNA sequence of AplVMP1 with a calculated molecular mass of 46.61 kDa is 1233 bp in length. AplVMP1 encodes PI class metalloproteinase with an open reading frame of 411 amino acid residues that includes signal peptide, pro-domain and metalloproteinase domains. The full-length cDNA of the AplVMP2 (1371 bp) has a calculated molecular mass of 51.16 kDa and encodes PII class metalloproteinase. The open reading frame of AplVMP2 with a 457 amino acid residues is composed of signal peptide, pro-domain, metalloproteinase and disintegrin domains. AplVMP1 and AplVMP2 showed 85% and 93% amino acid identical to PI class enzyme Agkistrodon contortrix laticinctus ACLPREF and PII class enzyme Agkistrodon piscivorus piscivorus piscivostatin, respectively. When expressed in Escherichia coli, most of recombinant proteins of AplVMP1 and AplVMP2 were in insoluble inclusion bodies, with soluble yields of 0.7 mg/l and 0.4 mg/l bacterial culture, respectively. Both affinity purified recombinant proteins show proteolytic activity on fibrinogen, although having an activity lower than that of crude A. p. leucostoma venom. Proteolytic activities of AplVMP1 and AplVMP2 were completely abolished after incubation with a final concentration of 100 microM of EDTA or 1,10-phenanthroline. Both AplVMP1 and AplVMP2 were active in a fibrin-agarose plate but devoid of hemorrhagic activity when injected (up to 50 microg) subcutaneously into mice, and had no capacity to inhibit platelet aggregation.

Transcriptomic basis for an antiserum against Micrurus corallinus (coral snake) venom

BACKGROUND

Micrurus corallinus (coral snake) is a tropical forest snake belonging to the family Elapidae. Its venom shows a high neurotoxicity associated with pre- and post-synaptic toxins, causing diaphragm paralysis, which may result in death. In spite of a relatively small incidence of accidents, serum therapy is crucial for those bitten. However, the adequate production of antiserum is hampered by the difficulty in obtaining sufficient amounts of venom from a small snake with demanding breeding conditions. In order to elucidate the molecular basis of this venom and to uncover possible immunogens for an antiserum, we generated expressed sequences tags (ESTs) from its venom glands and analyzed the transcriptomic profile. In addition, their immunogenicity was tested using DNA immunization.

RESULTS

A total of 1438 ESTs were generated and grouped into 611 clusters. Toxin transcripts represented 46% of the total ESTs. The two main toxin classes consisted of three-finger toxins (3FTx) (24%) and phospholipases A(2) (PLA(2)s) (15%). However, 8 other classes of toxins were present, including C-type lectins, natriuretic peptide precursors and even high-molecular mass components such as metalloproteases and L-amino acid oxidases. Each class included an assortment of isoforms, some showing evidence of alternative splicing and domain deletions. Five antigenic candidates were selected (four 3FTx and one PLA(2)) and used for a preliminary study of DNA immunization. The immunological response showed that the sera from the immunized animals were able to recognize the recombinant antigens.

CONCLUSION

Besides an improvement in our knowledge of the composition of coral snake venoms, which are very poorly known when compared to Old World elapids, the expression profile suggests abundant and diversified components that may be used in future antiserum formulation. As recombinant production of venom antigens frequently fails due to complex disulfide arrangements, DNA immunization may be a viable alternative. In fact, the selected candidates provided an initial evidence of the feasibility of this approach, which is less costly and not dependent on the availability of the venom.

Recombinant expression and affinity purification of snake venom gland parvalbumin in Escherichia coli

Parvalbumins (PV) are small, acidic, water soluble and calcium-binding proteins generally present in muscular and nervous tissues. In the present study, we identified and characterized a cDNA clone encoding PV, named AplPV, from a snake (Agkistrodon piscivorus leucostoma) venom gland cDNA library. AplPV belongs to EF-hand proteins with six alpha-helices constituting three EF-hand domains. The deduced amino acid sequence of AplPV is 91% and 68% identical to the previously characterized PVs of Boa constrictor and Cyprinus carpio, respectively. The full-length cDNA was subcloned into the expression vector pGEX and transformed into Escherichia coli (E.coli) to produce recombinant protein. The bacterially expressed GST-AplPV fusion protein was highly expressed, and effectively purified by Glutathione-Sepharose affinity chromatography. A high concentration of thrombin protease specifically cleaved and removed the GST tag from fusion protein, and further purified by Benzamidine column for removal of thrombin protease. As a result, the 12 kDa AplPV recombinant protein alone was purified. To investigate the tissue-specific biological occurrence of AplPV, a polyclonal antibody (anti-AplPV-antibody) was raised against GST-AplPV fusion protein in rabbit. Western blot analysis revealed that immunoreactive bands were exhibited in both recombinant protein and samples of venom glands, but not in any crude venom. This specific occurrence indicates a specialized function of AplPV in snake venom glands.