Identification of novel proteins from the venom of a cryptic snake Drysdalia coronoides by a combined transcriptomics and proteomics approach

Three finger toxins of elapids: structure, function, clinical applications and its inhibitors

The WHO lists snakebite as a "neglected tropical disease". In tropical and subtropical areas, envenoming is an important public health issue. This review article describes the structure, function, chemical composition, natural inhibitors, and clinical applications of Elapids' Three Finger Toxins (3FTX) using scientific research data. The primary venomous substance belonging to Elapidae is 3FTX, that targets nAChR. Three parallel β-sheets combine to create 3FTX, which has four or five disulfide bonds. The three primary types of 3FTX are short-chain, long-chain, and nonconventional 3FTX. The functions of 3FTX depend on the specific toxin subtype and the target receptor or ion channel. The well-known effect of 3FTX is probably neurotoxicity because of the severe consequences of muscular paralysis and respiratory failure in snakebite victims. 3FTX have also been studied for their potential clinical applications. α-bungarotoxin has been used as a molecular probe to study the structure and function of nAChRs (Nicotinic Acetylcholine Receptors). Acid-sensing ion channel (ASIC) isoforms 1a and 1b are inhibited by Mambalgins, derived from Black mamba venom, which hinders their function and provide an analgesic effect. α- Cobra toxin is a neurotoxin purified from Chinese cobra (Naja atra) binds to nAChR at the neuronal junction and causes an analgesic effect for moderate to severe pain. Some of the plants and their compounds have been shown to inhibit the activity of 3FTX, and their mechanisms of action are discussed.

Molecular phylogeny reveals distinct evolutionary lineages of the banded krait, Bungarus fasciatus (Squamata, Elapidae) in Asia

The banded krait, Bungarus fasciatus is a widespread elapid snake, likely to comprise several distinct species in different geographic regions of Asia. Therefore, based on molecular phylogenetics and comparative morphology data, we present an overview of the systematic composition of the species to delimit potential biogeographic boundaries. Our phylogenetic analyses, based on four mitochondrial genes, reveal the existence of at least three evolutionary lineages within B. fasciatus, corresponding to Indo-Myanmar, Sundaic and eastern Asian lineages. We are convinced that there are at least three taxonomic entities within the nomen B. fasciatus and restrict the distribution of B. fasciatus sensu stricto to the Indo-Myanmar region. We also provide additional natural history data of the taxon from eastern India. Finally, we advocate further studies to establish the degree of reproductive isolation among these diverging evolutionary lineages and to reassess the systematic status of this species complex especially the Sundaic and eastern Asian lineages.

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

Simple Summary

Around the world, snake envenomation poses a serious health risk. Proteins with pharmacological effects are present in snake venom. Recent studies elaborate snake venom and its potential application, including as a cancer drug and antibacterial substances. Our study aimed to analyze the global profile of the literature in snake venom research from documents indexed in the Scopus database between 1933 and 2022. In total, 2999 documents were published with Brazil showing the highest productivity. Antivenom, proteomics, and transcriptomics are emerging as hot topics on a global scale. The present study offers a distinctive overview of snake venom research conducted worldwide.

Abstract

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.

Transcriptomic Analysis of the Venom Gland and Enzymatic Characterization of the Venom of Phoneutria depilata (Ctenidae) from Colombia

The transcriptome of the venom glands of the Phoneutria depilata spider was analyzed using RNA-seq with an Illumina protocol, which yielded 86,424 assembled transcripts. A total of 682 transcripts were identified as potentially coding for venom components. Most of the transcripts found were neurotoxins (156) that commonly act on sodium and calcium channels. Nevertheless, transcripts coding for some enzymes (239), growth factors (48), clotting factors (6), and a diuretic hormone (1) were found, which have not been described in this spider genus. Furthermore, an enzymatic characterization of the venom of P. depilata was performed, and the proteomic analysis showed a correlation between active protein bands and protein sequences found in the transcriptome. The transcriptomic analysis of P. depilata venom glands show a deeper description of its protein components, allowing the identification of novel molecules that could lead to the treatment of human diseases, or could be models for developing bioinsecticides.

Venom-Derived Neurotoxins Targeting Nicotinic Acetylcholine Receptors

Acetylcholine was the first neurotransmitter described. The receptors targeted by acetylcholine are found within organisms spanning different phyla and position themselves as very attractive targets for predation, as well as for defense. Venoms of snakes within the Elapidae family, as well as those of marine snails within the Conus genus, are particularly rich in proteins and peptides that target nicotinic acetylcholine receptors (nAChRs). Such compounds are invaluable tools for research seeking to understand the structure and function of the cholinergic system. Proteins and peptides of venomous origin targeting nAChR demonstrate high affinity and good selectivity. This review aims at providing an overview of the toxins targeting nAChRs found within venoms of different animals, as well as their activities and the structural determinants important for receptor binding.

The Sequence and Three-Dimensional Structure Characterization of Snake Venom Phospholipases B

Snake venom phospholipases B (SVPLBs) are the least studied enzymes. They constitute about 1% of Bothrops crude venoms, however, in other snake venoms, it is present in less than 1%. These enzymes are considered the most potent hemolytic agent in the venom. Currently, no structural information is available about these enzymes from snake venom. To better understand its three-dimensional structure and mechanisms of envenomation, the current work describes the first model-based structure report of this enzyme from Bothrops moojeni venom named as B. moojeni phospholipase B (PLB_Bm). The structure model of PLB_Bm was generated using model building software like I-TESSER, MODELLER 9v19, and Swiss-Model. The build PLB_Bm model was validated using validation tools (PROCHECK, ERRAT, and Verif3D). The analysis of the PLB_Bm modeled structure indicates that it contains 491 amino acid residues that form a well-defined four-layer αββα sandwich core and has a typical fold of the N-terminal nucleophile aminohydrolase (Ntn-hydrolase). The overall structure of PLB_Bm contains 18 β-strands and 17 α-helices with many connecting loops. The structure divides into two chains (A and B) after maturation. The A chain is smaller and contains 207 amino acid residues, whereas the B chain is larger and contains 266 amino acid residues. The sequence and structural comparison among homologous snake venom, bacterial, and mammals PLBs indicate that differences in the length and sequence composition may confer variable substrate specificity to these enzymes. Moreover, the surface charge distribution, average volume, and depth of the active site cavity also vary in these enzymes. The present work will provide more information about the structure-function relationship and mechanism of action of these enzymes in snakebite envenomation.

The Urgent Need to Develop Novel Strategies for the Diagnosis and Treatment of Snakebites

Snakebite envenoming (SBE) is a priority neglected tropical disease, which kills in excess of 100,000 people per year. Additionally, many millions of survivors also suffer through disabilities and long-term health consequences. The only treatment for SBE, antivenom, has a number of major associated problems, not least, adverse reactions and limited availability. This emphasises the necessity for urgent improvements to the management of this disease. Administration of antivenom is too frequently based on symptomatology, which results in wasting crucial time. The majority of SBE-affected regions rely on broad-spectrum polyvalent antivenoms that have a low content of case-specific efficacious immunoglobulins. Research into small molecular therapeutics such as varespladib/methyl-varespladib (PLA2 inhibitors) and batimastat/marimastat (metalloprotease inhibitors) suggest that such adjunctive treatments could be hugely beneficial to victims. Progress into toxin-specific monoclonal antibodies as well as alternative binding scaffolds such as aptamers hold much promise for future treatment strategies. SBE is not implicit during snakebite, due to venom metering. Thus, the delay between bite and symptom presentation is critical and when symptoms appear it may often already be too late to effectively treat SBE. The development of reliable diagnostical tools could therefore initiate a paradigm shift in the treatment of SBE. While the complete eradication of SBE is an impossibility, mitigation is in the pipeline, with new treatments and diagnostics rapidly emerging. Here we critically review the urgent necessity for the development of diagnostic tools and improved therapeutics to mitigate the deaths and disabilities caused by SBE.

Drysdalin, an antagonist of nicotinic acetylcholine receptors highlights the importance of functional rather than structural conservation of amino acid residues

Snake venom neurotoxins are potent antagonists of nicotinic acetylcholine receptors (nAChRs). Here, we describe a novel member of class 3c long-chain neurotoxin drysdalin from the venom of Drysdalia coronoides. Drysdalin lacks three of the eight conserved classical functional residues critical for nAChRs interaction. Despite such a drastic alteration of the functional site, recombinant drysdalin showed irreversible postsynaptic neurotoxicity with nanomolar potency and selectively antagonizes the rodent muscle (α1)2β1δε, and human α7 and α9α10 nAChRs, but had no significant activity at the human α3β2, α3β4, α4β2, and α4β4 nAChRs. Substitution of Leu34 and Ala37 residues with the conserved Arg had minimal impact on the potency whereas conserved Phe replacement of residue Arg30 substantially reduced or abolished inhibitory activity. In contrast, truncation of the 24-residue long C-terminal tail leads to complete loss in (a) activity at α9α10 nAChR; and (b) irreversibility with reduced potency at the muscle and α7 nAChRs. Overall, the non-conserved Arg30 residue together with the uniquely long C-terminal tail contribute to the inhibitory activity of drysdalin at the nAChRs suggesting, at least for drysdalin, functional rather than sequence conservation plays a critical role in determining the activity of the toxin.

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.

Computational Studies of Snake Venom Toxins

Most snake venom toxins are proteins, and participate to envenomation through a diverse array of bioactivities, such as bleeding, inflammation, and pain, cytotoxic, cardiotoxic or neurotoxic effects. The venom of a single snake species contains hundreds of toxins, and the venoms of the 725 species of venomous snakes represent a large pool of potentially bioactive proteins. Despite considerable discovery efforts, most of the snake venom toxins are still uncharacterized. Modern bioinformatics tools have been recently developed to mine snake venoms, helping focus experimental research on the most potentially interesting toxins. Some computational techniques predict toxin molecular targets, and the binding mode to these targets. This review gives an overview of current knowledge on the ~2200 sequences, and more than 400 three-dimensional structures of snake toxins deposited in public repositories, as well as of molecular modeling studies of the interaction between these toxins and their molecular targets. We also describe how modern bioinformatics have been used to study the snake venom protein phospholipase A2, the small basic myotoxin Crotamine, and the three-finger peptide Mambalgin.

A Review and Database of Snake Venom Proteomes

Advances in the last decade combining transcriptomics with established proteomics methods have made possible rapid identification and quantification of protein families in snake venoms. Although over 100 studies have been published, the value of this information is increased when it is collated, allowing rapid assimilation and evaluation of evolutionary trends, geographical variation, and possible medical implications. This review brings together all compositional studies of snake venom proteomes published in the last decade. Compositional studies were identified for 132 snake species: 42 from 360 (12%) Elapidae (elapids), 20 from 101 (20%) Viperinae (true vipers), 65 from 239 (27%) Crotalinae (pit vipers), and five species of non-front-fanged snakes. Approximately 90% of their total venom composition consisted of eight protein families for elapids, 11 protein families for viperines and ten protein families for crotalines. There were four dominant protein families: phospholipase A₂s (the most common across all front-fanged snakes), metalloproteases, serine proteases and three-finger toxins. There were six secondary protein families: cysteine-rich secretory proteins, l-amino acid oxidases, kunitz peptides, C-type lectins/snaclecs, disintegrins and natriuretic peptides. Elapid venoms contained mostly three-finger toxins and phospholipase A₂s and viper venoms metalloproteases, phospholipase A₂s and serine proteases. Although 63 protein families were identified, more than half were present in <5% of snake species studied and always in low abundance. The importance of these minor component proteins remains unknown.

The medical threat of mamba envenoming in sub-Saharan Africa revealed by genus-wide analysis of venom composition, toxicity and antivenomics profiling of available antivenoms

Mambas (genus Dendroaspis) are among the most feared of the medically important elapid snakes found in sub-Saharan Africa, but many facets of their biology, including the diversity of venom composition, remain relatively understudied. Here, we present a reconstruction of mamba phylogeny, alongside genus-wide venom gland transcriptomic and high-resolution top-down venomic analyses. Whereas the green mambas, D. viridis, D. angusticeps, D. j. jamesoni and D. j. kaimosae, express 3FTx-predominant venoms, black mamba (D. polylepis) venom is dominated by dendrotoxins I and K. The divergent terrestrial ecology of D. polylepis compared to the arboreal niche occupied by all other mambas makes it plausible that this major difference in venom composition is due to dietary variation. The pattern of intrageneric venom variability across Dendroaspis represented a valuable opportunity to investigate, in a genus-wide context, the variant toxicity of the venom, and the degree of paraspecific cross-reactivity between antivenoms and mamba venoms. To this end, the immunological profiles of the five mamba venoms were assessed against a panel of commercial antivenoms generated for the sub-Saharan Africa market. This study provides a genus-wide overview of which available antivenoms may be more efficacious in neutralising human envenomings caused by mambas, irrespective of the species responsible. The information gathered in this study lays the foundations for rationalising the notably different potency and pharmacological profiles of Dendroaspis venoms at locus resolution. This understanding will allow selection and design of toxin immunogens with a view to generating a safer and more efficacious pan-specific antivenom against any mamba envenomation.

BIOLOGICAL SIGNIFICANCE

The mambas (genus Dendroaspis) comprise five especially notorious medically important venomous snakes endemic to sub-Saharan Africa. Their highly potent venoms comprise a high diversity of pharmacologically active peptides, including extremely rapid-acting neurotoxins. Previous studies on mamba venoms have focused on the biochemical and pharmacological characterisation of their most relevant toxins to rationalize the common neurological and neuromuscular symptoms of envenomings caused by these species, but there has been little work on overall venom composition or comparisons between them. Only very recently an overview of the composition of the venom of two Dendroaspis species, D. angusticeps and D. polylepis, has been unveiled through venomics approaches. Here we present the first genus-wide transcriptomic-proteomic analysis of mamba venom composition. The transcriptomic analyses described in this paper have contributed 29 (D. polylepis), 23 (D. angusticeps), 40 (D. viridis), 25 (D. j. jamesoni) and 21 (D. j. kaimosae), novel full-length toxin sequences to the non-redundant Dendroaspis sequence database. The mamba genus-wide venomic analysis demonstrated that major D. polylepis venom components are Kunitz-fold family toxins. This feature is unique in relation to the relatively conserved three-finger toxin (3FTx)-dominated venom compositions of the green mambas. Venom variation was interpreted in the context of dietary variation due to the divergent terrestrial ecology of D. polylepis compared to the arboreal niche occupied by all other mambas. Additionally, the degree of cross-reactivity conservation of mamba venoms was assessed by antivenomics against a panel of commercial antivenoms generated for the sub-Saharan Africa market. This study provides a genus-wide overview to infer which available antivenoms may be capable of neutralising human envenomings caused by mambas, irrespective of the species responsible. The information gathered in this study lays the foundations for rationalising the pharmacological profiles of mamba venoms at locus resolution. This understanding will contribute to the generation of a safer and more efficacious pan-Dendroaspis therapeutic antivenom against any mamba envenomation.

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.

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.

A combined proteomic and transcriptomic analysis of slime secreted by the southern bottletail squid, Sepiadarium austrinum (Cephalopoda)

Sepiadarium austrinum, the southern bottletail squid, is a small squid that inhabits soft sediments along Australia's south-east coast. When provoked, it rapidly secretes large volumes of slime, presumably as a form of chemical defense. We analyzed the proteomic composition of this slime using tandem mass spectrometry and transcriptomics and found that it was remarkably complex with 1735 identified protein groups (FDR:0.01). To investigate the chemical defense hypothesis we performed an Artemia toxicity assay and used sequence analysis to search for toxin-like molecules. Although the slime did not appear to be toxic to Artemia we found 13 proteins in slime with the hallmarks of toxins, namely cysteine richness, short length, a signal peptide and/or homology to known toxins. These included three short (80-130AA) cysteine rich secreted proteins with no homology to proteins on the NCBI or UniProt databases. Other protein families found included, CAP, phospholipase-B, ShKT-like peptides, peptidase S10, Kunitz BPTI and DNase II. Quantitative analysis using intensity based absolute quantification (iBAQ via MaxQuant) revealed 20 highly abundant proteins, accounting for 67% of iBAQ signal, and three of these were toxin-like. No mucin homologues were found suggesting that the structure of the slime gel may be formed by an unknown mechanism.

BIOLOGICAL SIGNIFICANCE

This study is the first known instance of a slime secretion from a cephalopod to be analyzed by proteomics methods and is the first investigation of a member of the family Sepiadariidae using proteomic methods. 1735 proteins were identified with 13 of these fitting criteria established for the identification of putative toxins. The slime is dominated by 20 highly abundant proteins with secreted, cysteine rich proteins. The study highlights the importance of 'omics approaches in understanding novel organisms.

Combined proteomics and transcriptomics identifies sting-related toxins of jellyfish Cyanea nozakii

Jellyfish sting has become a worldwide issue of critical concern to human health and safety in coastal areas in recent decades. Cyanea nozakii is one of the dominant blooming species and dangerous stingers in China. However, it remains unclear how many and what types of toxins are present in the venom. So, we used a combined transcriptomics and proteomics approach to investigate the venom composition of jellyfish C. nozakii. In total 4,608,524 Illumina valid reads were obtained to de novo assemble to 40,434 unigenes in the transcriptomics analysis. And, a total of 311,635 MS/MS spectra with 12,247 unique MS/MS spectra were generated to 1556 homologous proteins in the proteomics analysis. 174 potential toxin proteins were identified, with 27 proteins homology to the toxins from venomous animals, including phospholipase A2, zinc metalloproteinase-disintegrin agkistin, serine protease inhibitor, plancitoxin-1, alpha-latrocrustotoxin-Lt1a, etc. This study described the transcriptomics and venom proteomics of jellyfish C. nozakii for the first time. Our findings provide a comprehensive understanding of the venom composition of C. nozakii. Furthermore, the results may also be very helpful for the discovery of novel bioactive proteins, as well as the development of effective treatments for jellyfish sting in the future.

BIOLOGICAL SIGNIFICANCE

Jellyfish Cyanea nozakii is one of the most dangerous stingers in the coast of china. Hundreds of thousands of people would be stung every year and victims suffered a severe pain, itch, swelling, inflammation, wheal and even more serious consequence. However, it remains unclear how many and what types of toxins are present as well as the relationship between the clinical symptoms and toxins. Our combined transcriptomics and proteomics findings can provide a comprehensive understanding of the venom composition of C. nozakii and will also be helpful for the development of effective treatments for jellyfish sting in the future.

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.

Identification of hyaluronidase and phospholipase B in Lachesis muta rhombeata venom

Hyaluronidases contribute to local and systemic damages after envenoming, since they act as spreading factors cleaving the hyaluronan presents in the connective tissues of the victim, facilitating the diffusion of venom components. Although hyaluronidases are ubiquitous in snake venoms, they still have not been detected in transcriptomic analysis of the Lachesis venom gland and neither in the proteome of its venom performed previously. This work purified a hyaluronidase from Lachesis muta rhombeata venom whose molecular mass was estimated by SDS-PAGE to be 60 kDa. The hyaluronidase was more active at pH 6 and 37 °C when salt concentration was kept constant and more active in the presence of 0.15 M monovalent ions when the pH was kept at 6. Venom was fractionated by reversed-phase liquid chromatography (RPLC). Edman sequencing after RPLC failed to detect hyaluronidase, but identified a new serine proteinase isoform. The hyaluronidase was identified by mass spectrometry analysis of the protein bands in SDS-PAGE. Additionally, phospholipase B was identified for the first time in Lachesis genus venom. The discovery of new bioactive molecules might contribute to the design of novel drugs and biotechnology products as well as to development of more effective treatments against the envenoming.

Discovery of toxin-encoding genes from the false viper Macropisthodon rudis, a rear-fanged snake, by transcriptome analysis of venom gland

Although rear-fanged snakes are often considered as non-threatening to humans, some species are lethal or medically hazardous. The toxin components and bioactivities of front-fanged snakes have been extensively studied; however, only limited research has explored the venoms of rear-fanged snakes. The false viper, Macropisthodon rudis, is widespread in southern China, but little is known about the toxins that this snake produces. Here, we analyzed the transcriptome of the venom gland of M. rudis using high-throughput sequencing with an illumina HiSeq 2000. The raw data were assembled and annotated using public databases. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and gene ontology (GO) were analyzed. Using sequence comparisons, snake venom metalloproteinases (SVMPs) and a phosphodiesterase (PDE) were discovered in the venom gland of M. rudis.

Venom-gland transcriptome and venom proteome of the Malaysian king cobra (Ophiophagus hannah)

Background

The king cobra (Ophiophagus hannah) is widely distributed throughout many parts of Asia. This study aims to investigate the complexity of Malaysian Ophiophagus hannah (MOh) venom for a better understanding of king cobra venom variation and its envenoming pathophysiology. The venom gland transcriptome was investigated using the Illumina HiSeq™ platform, while the venom proteome was profiled by 1D-SDS-PAGE-nano-ESI-LCMS/MS.

Results

Transcriptomic results reveal high redundancy of toxin transcripts (3357.36 FPKM/transcript) despite small cluster numbers, implying gene duplication and diversification within restricted protein families. Among the 23 toxin families identified, three-finger toxins (3FTxs) and snake-venom metalloproteases (SVMPs) have the most diverse isoforms. These 2 toxin families are also the most abundantly transcribed, followed in descending order by phospholipases A₂ (PLA₂s), cysteine-rich secretory proteins (CRISPs), Kunitz-type inhibitors (KUNs), and L-amino acid oxidases (LAAOs). Seventeen toxin families exhibited low mRNA expression, including hyaluronidase, DPP-IV and 5’-nucleotidase that were not previously reported in the venom-gland transcriptome of a Balinese O. hannah. On the other hand, the MOh proteome includes 3FTxs, the most abundantly expressed proteins in the venom (43 % toxin sbundance). Within this toxin family, there are 6 long-chain, 5 short-chain and 2 non-conventional 3FTx. Neurotoxins comprise the major 3FTxs in the MOh venom, consistent with rapid neuromuscular paralysis reported in systemic envenoming. The presence of toxic enzymes such as LAAOs, SVMPs and PLA₂ would explain tissue inflammation and necrotising destruction in local envenoming. Dissimilarities in the subtypes and sequences between the neurotoxins of MOh and Naja kaouthia (monocled cobra) are in agreement with the poor cross-neutralization activity of N. kaouthia antivenom used against MOh venom. Besides, the presence of cobra venom factor, nerve growth factors, phosphodiesterase, 5’-nucleotidase, and DPP-IV in the venom proteome suggests its probable hypotensive action in subduing prey.

Conclusion

This study reports the diversity and abundance of toxins in the venom of the Malaysian king cobra (MOh). The results correlate with the pathophysiological actions of MOh venom, and dispute the use of Naja cobra antivenoms to treat MOh envenomation. The findings also provide a deeper insight into venom variations due to geography, which is crucial for the development of a useful pan-regional antivenom.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1828-2) contains supplementary material, which is available to authorized users.

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.

Interrogating the Venom of the Viperid Snake Sistrurus catenatus edwardsii by a Combined Approach of Electrospray and MALDI Mass Spectrometry

The complete sequence characterization of snake venom proteins by mass spectrometry is rather challenging due to the presence of multiple isoforms from different protein families. In the present study, we investigated the tryptic digest of the venom of the viperid snake Sistrurus catenatus edwardsii by a combined approach of liquid chromatography coupled to either electrospray (online) or MALDI (offline) mass spectrometry. These different ionization techniques proved to be complementary allowing the identification a great variety of isoforms of diverse snake venom protein families, as evidenced by the detection of the corresponding unique peptides. For example, ten out of eleven predicted isoforms of serine proteinases of the venom of S. c. edwardsii were distinguished using this approach. Moreover, snake venom protein families not encountered in a previous transcriptome study of the venom gland of this snake were identified. In essence, our results support the notion that complementary ionization techniques of mass spectrometry allow for the detection of even subtle sequence differences of snake venom proteins, which is fundamental for future structure-function relationship and possible drug design studies.

Testing the Toxicofera: comparative transcriptomics casts doubt on the single, early evolution of the reptile venom system

The identification of apparently conserved gene complements in the venom and salivary glands of a diverse set of reptiles led to the development of the Toxicofera hypothesis - the single, early evolution of the venom system in reptiles. However, this hypothesis is based largely on relatively small scale EST-based studies of only venom or salivary glands and toxic effects have been assigned to only some putative Toxicoferan toxins in some species. We set out to examine the distribution of these proposed venom toxin transcripts in order to investigate to what extent conservation of gene complements may reflect a bias in previous sampling efforts. Our quantitative transcriptomic analyses of venom and salivary glands and other body tissues in five species of reptile, together with the use of available RNA-Seq datasets for additional species, shows that the majority of genes used to support the establishment and expansion of the Toxicofera are in fact expressed in multiple body tissues and most likely represent general maintenance or "housekeeping" genes. The apparent conservation of gene complements across the Toxicofera therefore reflects an artefact of incomplete tissue sampling. We therefore conclude that venom has evolved multiple times in reptiles.

Venomics and antivenomics of Bothrops erythromelas from five geographic populations within the Caatinga ecoregion of northeastern Brazil

The Caatinga lancehead, Bothrops erythromelas, is a medically relevant species, responsible for most of the snakebite accidents in most parts of its distribution range in northeastern Brazil. The spectrum and geographic variability of its venom toxins were investigated applying a venomics approach to venom pools from five geographic areas within the Caatinga ecoregion. Despite its wide habitat, populations of B. erythromelas from Ceará, Pernambuco, Juazeiro, Paraiba, and Ilha de Itaparica exhibit highly conserved venom proteomes. Mirroring their compositional conservation, the five geographic venom pools also showed qualitatively and quantitatively overlapping antivenomic profiles against antivenoms generated in Vital Brazil (BR) and Clodomiro Picado (CR) Institutes, using different venoms in the immunization mixtures. The paraspecificity exhibited by the Brazilian SAB and the Costa Rican BCL antivenoms against venom toxins from B. erythromelas indicates large immunoreactive epitope conservation across genus Bothrops during the last ~14 million years, thus offering promise for the possibility of generating a broad-spectrum bothropic antivenom. Biological Significance Accidental snakebite envenomings represent an important public health hazard in Brazil. Ninety per cent of the yearly estimated 20-30,000 snakebite accidents are caused by species of the Bothrops genus. Bothrops erythromelas, a small, moderately stocky terrestrial venomous snake, is responsible for most of the snakebite accidents in its broad distribution range in the Caatinga, a large ecoregion in northeastern Brazil. To gain a deeper insight into the spectrum of medically important toxins present in the venom of the Caatinga lancehead, we applied a venomics approach to define the proteome and geographic variability of adult B. erythromelas venoms from five geographic regions. Although intraspecific compositional variation between venoms among specimens from different geographic regions has long been appreciated by herpetologists and toxinologists as a general feature of highly adaptable and widely distributed snake species, the five B. erythromelas populations investigated exhibit highly conserved venom proteomes. The overall toxin profile of the Caatinga lancehead's venom explains the local and systemic effects observed in envenomations by B. erythromelas. The five geographic venom pools sampled also showed qualitatively and quantitatively overlapping antivenomic profiles against antivenoms generated using different bothropic venoms in the immunization mixtures. The large immunoreactive epitope conservation across genus Bothrops offers promise for the generation of a broad-spectrum bothropic antivenom.

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.

Isolation and structural and pharmacological characterization of α-elapitoxin-Dpp2d, an amidated three finger toxin from black mamba venom

We isolated a novel, atypical long-chain three-finger toxin (TFT), α-elapitoxin-Dpp2d (α-EPTX-Dpp2d), from black mamba (Dendroaspis polylepis polylepis) venom. Proteolytic digestion with trypsin and V8 protease, together with MS/MS de novo sequencing, indicated that the mature toxin has an amidated C-terminal arginine, a posttranslational modification rarely observed for snake TFTs. α-EPTX-Dpp2d was found to potently inhibit α7 neuronal nicotinic acetylcholine receptors (nAChR; IC₅₀, 58 ± 24 nM) and muscle-type nAChR (IC₅₀, 114 ± 37 nM) but did not affect α3β2 and α3β4 nAChR isoforms at 1 μM concentrations. Competitive radioligand binding assays demonstrated that α-EPTX-Dpp2d competes with epibatidine binding to the Lymnea stagnalis acetylcholine-binding protein (Ls-AChBP; IC₅₀, 4.9 ± 2.3 nM). The activity profile and binding data are reminiscent of classical long-chain TFTs with a free carboxyl termini, suggesting that amidation does not significantly affect toxin selectivity. The crystal structure of α-EPTX-Dpp2d was determined at 1.7 Å resolution and displayed a dimeric toxin assembly with each monomer positioned in an antiparallel orientation. The dimeric structure is stabilized by extensive intermolecular hydrogen bonds and electrostatic interactions, which raised the possibility that the toxin may exist as a noncovalent homodimer in solution. However, chemical cross-linking and size-exclusion chromatography coupled with multiangle laser light scattering (MALLS) data indicated that the toxin is predominantly monomeric under physiological conditions. Because of its high potency and selectivity, we expect this toxin to be a valuable pharmacological tool for studying the structure and function of nAChRs.

Quantitative high-throughput profiling of snake venom gland transcriptomes and proteomes (Ovophis okinavensis and Protobothrops flavoviridis)

Background

Advances in DNA sequencing and proteomics have facilitated quantitative comparisons of snake venom composition. Most studies have employed one approach or the other. Here, both Illumina cDNA sequencing and LC/MS were used to compare the transcriptomes and proteomes of two pit vipers, Protobothrops flavoviridis and Ovophis okinavensis, which differ greatly in their biology.

Results

Sequencing of venom gland cDNA produced 104,830 transcripts. The Protobothrops transcriptome contained transcripts for 103 venom-related proteins, while the Ovophis transcriptome contained 95. In both, transcript abundances spanned six orders of magnitude. Mass spectrometry identified peptides from 100% of transcripts that occurred at higher than contaminant (e.g. human keratin) levels, including a number of proteins never before sequenced from snakes. These transcriptomes reveal fundamentally different envenomation strategies. Adult Protobothrops venom promotes hemorrhage, hypotension, incoagulable blood, and prey digestion, consistent with mammalian predation. Ovophis venom composition is less readily interpreted, owing to insufficient pharmacological data for venom serine and metalloproteases, which comprise more than 97.3% of Ovophis transcripts, but only 38.0% of Protobothrops transcripts. Ovophis venom apparently represents a hybrid strategy optimized for frogs and small mammals.

Conclusions

This study illustrates the power of cDNA sequencing combined with MS profiling. The former quantifies transcript composition, allowing detection of novel proteins, but cannot indicate which proteins are actually secreted, as does MS. We show, for the first time, that transcript and peptide abundances are correlated. This means that MS can be used for quantitative, non-invasive venom profiling, which will be beneficial for studies of endangered species.

Predicting function from sequence in a large multifunctional toxin family

Venoms contain active substances with highly specific physiological effects and are increasingly being used as sources of novel diagnostic, research and treatment tools for human disease. Experimental characterisation of individual toxin activities is a severe rate-limiting step in the discovery process, and in-silico tools which allow function to be predicted from sequence information are essential. Toxins are typically members of large multifunctional families of structurally similar proteins that can have different biological activities, and minor sequence divergence can have significant consequences. Thus, existing predictive tools tend to have low accuracy. We investigated a classification model based on physico-chemical attributes that can easily be calculated from amino-acid sequences, using over 250 (mostly novel) viperid phospholipase A₂ toxins. We also clustered proteins by sequence profiles, and carried out in-vitro tests for four major activities on a selection of isolated novel toxins, or crude venoms known to contain them. The majority of detected activities were consistent with predictions, in contrast to poor performance of a number of tested existing predictive methods. Our results provide a framework for comparison of active sites among different functional sub-groups of toxins that will allow a more targeted approach for identification of potential drug leads in the future.

Protease inhibitors from marine venomous animals and their counterparts in terrestrial venomous animals

The Kunitz-type protease inhibitors are the best-characterized family of serine protease inhibitors, probably due to their abundance in several organisms. These inhibitors consist of a chain of ~60 amino acid residues stabilized by three disulfide bridges, and was first observed in the bovine pancreatic trypsin inhibitor (BPTI)-like protease inhibitors, which strongly inhibit trypsin and chymotrypsin. In this review we present the protease inhibitors (PIs) described to date from marine venomous animals, such as from sea anemone extracts and Conus venom, as well as their counterparts in terrestrial venomous animals, such as snakes, scorpions, spiders, Anurans, and Hymenopterans. More emphasis was given to the Kunitz-type inhibitors, once they are found in all these organisms. Their biological sources, specificity against different proteases, and other molecular blanks (being also K⁺ channel blockers) are presented, followed by their molecular diversity. Whereas sea anemone, snakes and other venomous animals present mainly Kunitz-type inhibitors, PIs from Anurans present the major variety in structure length and number of Cys residues, with at least six distinguishable classes. A representative alignment of PIs from these venomous animals shows that, despite eventual differences in Cys assignment, the key-residues for the protease inhibitory activity in all of them occupy similar positions in primary sequence. The key-residues for the K⁺ channel blocking activity was also compared.

Characterization of monomeric and multimeric snake neurotoxins and other bioactive proteins from the venom of the lethal Australian common copperhead (Austrelaps superbus)

Envenomation by Australian copperheads results mainly in muscle paralysis largely attributed to the presence of postsynaptic α-neurotoxins. However, poorly reversible neurotoxic effects suggest that these venoms may contain snake presynaptic phospholipase A2 neurotoxins (SPANs) that irreversibly inhibit neurotransmitter release. Using size-exclusion liquid chromatography, the present study isolated the first multimeric SPAN complex from the venom of the Australian common copperhead, Austrelaps superbus. The multimeric SPAN P-elapitoxin-As1a (P-EPTX-As1a) along with two novel monomeric SPANs and a new postsynaptic α-neurotoxin were then pharmacologically characterized using the chick biventer cervicis nerve-muscle preparation. All SPANs inhibited nerve-evoked twitch contractions at the neuromuscular junction without inhibiting contractile responses to cholinergic agonists or KCl. These actions are consistent with a prejunctional action to inhibit neurotransmitter release, without direct myotoxicity. Furthermore, the multimeric P-EPTX-As1a caused tetanic 'fade' in muscle tension under high frequency nerve stimulation, and produced a triphasic alteration to neurotransmitter release. These actions have been previously noted with other multimeric SPAN complexes such as taipoxin. Moreover, the neurotoxic α-subunit of P-EPTX-As1a shows high homology to taipoxin α-chain. Several other coagulopathic and myotoxic high mass proteins including a class PIII snake venom metalloproteinase, C-type lectin, l-amino acid oxidase, acetylcholinesterase and phospholipase B were also identified that may contribute to the overall toxicity of A. superbus venom. In conclusion, clinicians should be aware that early antivenom intervention might be necessary to prevent the onset of irreversible presynaptic neurotoxicity caused by multimeric and monomeric SPANs and that A. superbus venom is potentially capable of producing coagulopathic and myotoxic effects.

A brief review of the scientific history of several lesser-known snake venom proteins: l-amino acid oxidases, hyaluronidases and phosphodiesterases

When considering the proteins and toxins in snake venom one's thoughts generally migrate to the proteinases, neurotoxins and phospholipases since these families of proteins are comprised by many of the toxins found in venom. However as modern proteomic and transcriptomic venom research has abundantly shown snake venoms are complex and containing numerous families of protein beyond the "big three". In this brief review we will discuss three of the lesser discussed proteins typically found in snake venoms: l-amino acid oxidases (LAAO); hyaluronidases and phosphodiesterases. These proteins have long been known to be part of many venoms' proteomes with reports appearing in the literature as early as 1944 for LAAO, 1947 for hyaluronidase (spreading factor), and 1932 for venom phosphodiesterase. These are more or less contemporary with the first reports (circa 1950) on snake venom proteases. Thus, the relatively modest literature on these snake venom proteins stems not from lack of early discovery but rather more likely to their ostensibly minor role in snake venom pathophysiology. In this review we will provide an overview of the experimental history of these venom proteins, their biochemical and structural features and their role in snake venom toxinology with the aim of bringing a fuller, more comprehensive, understanding of the history of laboratory research on snake venoms. In addition, there are some comments on these proteins from investigators who were actively engaged in their investigation.

Snake venomics of two poorly known Hydrophiinae: Comparative proteomics of the venoms of terrestrial Toxicocalamus longissimus and marine Hydrophis cyanocinctus

The venom proteomes of Toxicocalamus longissimus and Hydrophis cyanocinctus, a fossorial and a marine species, respectively, of the Hydrophiinae genus of Elapidae, were investigated by Edman degradation of RP-HPLC isolated proteins, and de novo MS/MS sequencing of in-gel derived tryptic peptide ions. The toxin arsenal of T. longissimus is made up of 1-2 type-I PLA(2) molecules, which account for 6.5% of the venom proteins, a minor PIII-SVMP (1.4% of the venom toxins), and ~20 members of the 3FTx family comprising 92% of the venom proteome. Seventeen proteins (5 type-I PLA(2)s and 12 3FTxs) were found in the venom of H. cyanocinctus. Three-finger toxins and type-I PLA(2) proteins comprise, respectively, 81% and 19% of its venom proteome. The simplicity of the H. cyanocinctus venom proteome is highlighted by the fact that only 6 venom components (3 short-chain neurotoxins, two long-chain neurotoxins, and one PLA(2) molecule) exhibit relative abundances >5%. As expected from its high neurotoxin abundance, the LD(50) for mice of H. cyanocinctus venom was fairly low, 0.132μg/g (intravenous) and 0.172μg/g (intraperitoneal). Our data indicate that specialization towards a lethal cocktail of 3FTx and type-I PLA(2) molecules may represent a widely adopted trophic solution throughout the evolution of Elapidae. Our results also points to a minimization of the molecular diversity of the toxin arsenal of the marine snake Hydrophis cyanocinctus in comparison to the venom proteome of its terrestrial relatives, and highlight that the same evolutionary solution, economy of the toxin arsenal, has been convergently adopted by different taxa in response to opposite selective pressures, loss and gain of neurotoxicity.

The structural and functional diversification of the Toxicofera reptile venom system

The evolutionary origin and diversification of the reptilian venom system is described. The resolution of higher-order molecular phylogenetics has clearly established that a venom system is ancestral to snakes. The diversification of the venom system within lizards is discussed, as is the role of venom delivery in the behavioural ecology of these taxa (particularly Varanus komodoensis). The more extensive diversification of the venom system in snakes is summarised, including its loss in some clades. Finally, we discuss the contentious issue of a definition for "venom", supporting an evolutionary definition that recognises the homology of both the venom delivery systems and the toxins themselves.