From genome to "venome": molecular origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences and related body proteins

Venomic analyses of Scolopendra viridicornis nigra and Scolopendra angulata (Centipede, Scolopendromorpha): Shedding light on venoms from a neglected group

Centipedes are venomous arthropods responsible for a significant number of non-lethal human envenomations. Despite this, information about the composition and function of their venom contents is scarce. In this study, we have used a 'structure to function' proteomic approach combining two-dimensional chromatography (2D-LC), electrospray ionization quadrupole/time-of-flight mass spectrometry (ESI-Q-TOF/MS), N-terminal sequencing and similarity searching to better understand the complexities of the venoms from two Brazilian centipede species: Scolopendra viridicornis nigra and Scolopendra angulata. Comparisons between the LC profiles and the mass compositions of the venoms of the two species are provided. The observed molecular masses ranged from 3019.62 to 20996.94Da in S. viridicornis nigra (total: 62 molecular masses) and from 1304.73 to 22639.15Da in S. angulata (total: 65 molecular masses). Also, the N-termini of representatives of 10 protein/peptide families were successfully sequenced where nine of them showed no significant similarity to other protein sequences deposited in the Swiss-Prot database. A screening for insecto-toxic activities in fractions from S. viridicornis venom has also been performed. Six out of the 12 tested fractions were responsible for clear toxic effects in house flies. This work demonstrates that centipede venoms might be a neglected but important source of new bioactive compounds.

Computer modeling of binding of diverse weak toxins to nicotinic acetylcholine receptors

Weak toxins are the "three-fingered" snake venoms toxins grouped together by having an additional disulfide in the N-terminal loop I. In general, weak toxins have low toxicity, and biological targets have been identified for some of them only, recently by detecting the effects on the nicotinic acetylcholine receptors (nAChR). Here the methods of docking and molecular dynamics simulations are used for comparative modeling of the complexes between four weak toxins of known spatial structure (WTX, candoxin, bucandin, gamma-bungarotoxin) and nAChRs. WTX and candoxin are those toxins whose blocking of the neuronal alpha7- and muscle-type nAChR has been earlier shown in binding assays and electrophysiological experiments, while for the other two toxins no such activity has been reported. Only candoxin and WTX are found here to give stable solutions for the toxin-nAChR complexes. These toxins appear to approach the binding site similarly to short alpha-neurotoxins, but their final position resembles that of alpha-cobratoxin, a long alpha-neurotoxin, in the complex with the acetylcholine-binding protein. The final spatial structures of candoxin and WTX complexes with the alpha7 neuronal or muscle-type nAChR are very similar and do not provide immediate answer why candoxin has a much higher affinity than WTX, but both of them share a virtually irreversible mode of binding to one or both these nAChR subtypes. Possible explanation comes from docking and MD simulations which predict fast kinetics of candoxin association with nAChR, no gross changes in the toxin conformation (with smaller toxin flexibility on alpha7 nAChR), while slow WTX binding to nAChR is associated with slow irreversible rearrangement both of the tip of the toxin loop II and of the binding pocket residues locking finally the toxin molecule. Computer modeling showed that the additional disulfide in the loop I is not directly involved in receptor binding of WTX and candoxin, but it stabilizes the structure of loop I which plays an important role in toxin delivery to the binding site. In summary, computer modeling visualized possible modes of binding for those weak toxins which interact with the nAChR, provided no solutions for those weak toxins whose targets are not the nAChRs, and demonstrated that the additional disulfide in loop I cannot be a sound criteria for joining all weak toxins into one group; the conclusion about the diversity of weak toxins made from computer modeling is in accord with the earlier phylogenetic analysis.

Novel families of toxin-like peptides in insects and mammals: a computational approach

Most animal toxins are short proteins that appear in venom and vary in sequence, structure and function. A common characteristic of many such toxins is their apparent structural stability. Sporadic instances of endogenous toxin-like proteins that function in non-venom context have been reported. We have utilized machine learning methodology, based on sequence-derived features and guided by the notion of structural stability, in order to conduct a large-scale search for toxin and toxin-like proteins. Application of the method to insect and mammalian sequences revealed novel families of toxin-like proteins. One of these proteins shows significant similarity to ion channel inhibitors that are expressed in cone snail and assassin bug venom, and is surprisingly expressed in the bee brain. A toxicity assay in which the protein was injected to fish induced a strong yet reversible paralytic effect. We suggest that the protein may function as an endogenous modulator of voltage-gated Ca(2+) channels. Additionally, we have identified a novel mammalian cluster of toxin-like proteins that are expressed in the testis. We suggest that these proteins might be involved in regulation of nicotinic acetylcholine receptors that affect the acrosome reaction and sperm motility. Finally, we highlight a possible evolutionary link between venom toxins and antibacterial proteins. We expect our methodology to enhance the discovery of additional novel protein families.

Antimicrobial activity of omwaprin, a new member of the waprin family of snake venom proteins

We have isolated and characterized omwaprin, a 50-amino-acid cationic protein from the venom of inland taipan (Oxyuranus microlepidotus). It is a new member of the waprin family of snake venom proteins. A synthetic gene was designed and constructed for expressing the recombinant protein in Escherichia coli. Recombinant omwaprin was used for carrying out functional analyses. The protein is non-toxic to Swiss albino mice at doses of up to 10 mg/kg when administered intraperitoneally. However, it shows selective and dose-dependant antibacterial activity against Gram-positive bacteria. The minimum inhibitory doses were in the range 2-10 microg for selected species of bacteria in radial diffusion assays. The antibacterial activity is salt-tolerant up to 350 mM NaCl. However, omwaprin lost its antibacterial activity upon reduction and alkylation of its cysteine residues, or upon deletion of six N-terminal amino acid residues, four of which are positively charged. These observations indicate that the three-dimensional structure constrained by four disulfide bonds and the N-terminal residues are essential for its activity. The mechanism of action is via membrane disruption, as shown by scanning electron microscopy. Importantly, omwaprin lacks haemolytic activity on human erythrocytes. This demonstrates the specificity of omwaprin for bacterial membranes. Unlike other reported WAP (whey acidic protein) domain-containing antibacterial proteins, including elafin, EPPIN (epididymal proteinase inhibitor), SWAM1 and SWAM2 [single WAP (whey acidic protein) motif proteins 1 and 2] and SLPI (secretory leucocyte proteinase inhibitor), omwaprin shows species-specific activity on the Gram-positive bacteria tested.

Neurotropic effects of venoms and other factors that promote prey acquisition

Mammals envenomed by either the Eastern diamondback rattlesnake (Crotalus adamanteus) or the cottonmouth (Agkistrodon piscivorus piscivorus) exhibit an immediate but transitory pupillar contraction, a parasympathomimetic effect mediated through the ciliary ganglion that can be prevented by a retrobulbar injection of anesthetic. The venom of the cottonmouth injected into the lymph spaces of the frog (Rana pipiens) produces an immediate and total collapse of the lung sacs. Applied locally to the surface, it produces a constriction that eventually collapses the entire sac. Tests of venoms and toxins from both anterior and posterior parts of the venom apparatus indicate that the lung-collapsing moiety originates in the accessory, not the main portion of the venom gland. This is the first example of a functional specialization within the whole structure. It seems that this factor is elaborated primarily in snakes that prey upon frogs, although insufficient data are available from this study to confirm this. In both reptile species, the predatory strike is accompanied by an immediate effect, perhaps mediated by the parasympathetic nervous system, designed to incapacitate the prey and facilitate capture. These effects cannot now be attributed to neurotoxins because the effect of the former is transitory (and not lethal) and neither has been purified sufficiently to determine potency or structure. Both take part in securing, but not killing, the prey, and both directly oppose the sympathetic nervous system "fright-fight/flight" response. Evidence is presented to support the possibility that known epigenetic mechanisms are capable of effecting heritable changes in gene expression that could allow for the development of factors that facilitate prey acquisition and promote rapid adaptation to environmental change.

Structure of two genes encoding parallel prothrombin activators in Tropidechis carinatus snake: gene duplication and recruitment of factor X gene to the venom gland

BACKGROUND

Proteins with new function originate through gene duplication followed by divergence. In nature, occurrence of structurally and functionally similar proteins performing highly diverse physiological roles within an organism is rare. Several Australian elapid snakes have two parallel prothrombin activating systems with distinct physiological roles. For example, in Tropidechis carinatus, trocarin D exists in the venom and acts as toxin, whereas coagulation factor (F) X exists in plasma and plays a role in hemostasis.

RESULTS

Here, we show that FX and the trocarin D genes are expressed in a highly tissue-specific manner in T. carinatus. To understand their origin, recruitment and evolution, we determined the complete structure of their genes. Both genes have eight exons with identical exon-intron boundaries. All the introns are 92-99% identical with the exception of intron 1, indicating a recent gene duplication event. The first intron of the trocarin D gene is also nearly identical to that of the FX gene, except for two deletions (255 and 1406 bp) and three insertions (214, 1975, and 2174 bp). The third insertion has a potential scaffold/matrix attached region. The putative promoter of the snake FX gene shares similar cis-elements compared with those of human and murine FX genes. Interestingly, the trocarin D promoter has a 264-bp insertion with core promoter sequences and cis-elements that are known to induce high-level expression. This insertion might be responsible for switching from constitutive expression of the FX gene to inducible expression of trocarin D. Thus, we named this segment as VERSE (Venom Recruitment/Switch Element).

CONCLUSION

To our knowledge, this is the first molecular evidence for the recruitment of a duplicated gene for expression in venom glands by a simple insertion.

Loss of introns along the evolutionary diversification pathway of snake venom disintegrins evidenced by sequence analysis of genomic DNA from Macrovipera lebetina transmediterranea and Echis ocellatus

Analysis of cDNAs from Macrovipera lebetina transmediterranea (Mlt) and Echis ocellatus (Eo) venom gland libraries encoding disintegrins argued strongly for a common ancestry of the messengers of short disintegrins and those for precursors of dimeric disintegrin chains. We now report the sequence analysis of disintegrin-coding genes from these two vipers. Genomic DNAs for dimeric disintegrin subunits Ml_G1 and Ml_G2 (Mlt) and Eo_D3 (Eo) contain single 1-kb introns exhibiting the 5'-GTAAG (donor)/3'-AG (acceptor) consensus intron splicing signature. On the other hand, the short RTS-disintegrins Ml_G3 (Mlt) and Eo_RTS (Eo) and the short RGD-disintegrin ocellatusin (Eo) are transcribed from intronless genomic DNA sequences, indicating that the evolutionary pathway leading to the emergence of short disintegrins involved the removal of all intronic sequences. The insertion position of the intron within Ml_G1, Ml_G2, and Eo_D3 is conserved in the genes for vertebrate ADAM (A disintegrin and metalloproteinase) protein disintegrin-like domains and within the gene for the medium-size snake disintegrins halystatins 2 and 3. However, a comparative analysis of currently available disintegrin(-like) genes outlines the view that a minimization of both the gene organization and the protein structure underlies the evolution of the snake venom disintegrin family.

Identification of a human ortholog of the mouseDcppgene locus, encoding a novel member of the CSP-1/Dcpp salivary protein family

Salivary fluid, the collective product of numerous major and minor salivary glands, contains a range of secretory proteins that play key defensive, digestive, and gustatory roles in the oral cavity. To understand the distinct protein “signature” contributed by individual salivary glands to salivary secretions, we studied a family of proteins shown by in vitro mRNA translation to be abundantly expressed in mouse sublingual glands. Molecular cloning, Southern blotting, and restriction fragment length polymorphism analyses showed these to represent one known and two novel members of the common salivary protein (CSP-1)/Demilune cell and parotid protein (Dcpp) salivary protein family, the genes for which are closely linked in the T-complex region of mouse chromosome 17. Bioinformatic analysis identified a putative human CSP-1/Dcpp ortholog, HRPE773, expressed predominantly in human salivary tissue, that shows 31% amino acid identity and 45% amino acid similarity to the mouse Dcpp query sequence. The corresponding human gene displays a similar structure to the mouse Dcpp genes and is located on human chromosome 16 in a region known to be syntenic with the T-complex region of mouse chromosome 17. The predicted mouse and human proteins both display classical NH2-terminal signal sequences, putative jacalin-related lectin domains, and potential N-linked glycosylation sites, suggesting secretion via sublingual saliva into the oral cavity where they may display antimicrobial activity or provide a defensive coating to enamel. Identification of a human CSP-1/Dcpp ortholog therefore provides a key tool for investigation of salivary protein function in human oral health and disease.

Molecular evolution of PIII-SVMP and RGD disintegrin genes from the genus Crotalus

Several types of disintegrins have been isolated from Crotalus spp rattlesnakes, including RGD disintegrins, and PIII-SVMPs. We isolated six cDNAs from snake venom glands using RT-PCR. Three RGD disintegrins (atroxatin, mojastin, and viridistatin) and three PIII-SVMPs (catroriarin, scutiarin, and viristiarin) cDNAs were isolated from the rattlesnakes Crotalus atrox, Crotalus scutulatus scutulatus, and Crotalus viridis viridis, respectively. Atroxatin and Viridistatin shared 90% amino acid identity to each other, and 87% identity to Mojastin. Scutiarin and Viristiarin were identical. All PIII-SVMPs isolated in this study shared the highest amino acid identity with Catrocollastatin. cDNA and protein sequences for RGD disintegrins, one MVD disintegrin, and PIII-SVMPs of the genus Crotalus (present in the NCBI database), were used in phylogenetic analysis. Neighbor-joining analysis of PIII-SVMP and RGD/MVD disintegrin-coding DNA sequences showed that these groups of genes separate into separate clades. A Phi(ST) pairwise comparison and Analysis of Molecular Variance (AMOVA) between PIII-SVMPs and RGD/MVD disintegrins showed significant genetic differences. Mutations observed in ten of the cDNAs analyzed did not affect Cys-coding sequences. Our K(A)/K(S) data suggest that rapid evolution occurred between the genes coding for PIII-SVMPs resulting, in the production of RGD disintegrin-coding genes. However, once these genes diverged, mutations in the PIII-SVMP-coding genes were accumulated less frequently.

Some aspects of the venom proteome of the Colubridae snakePhilodryas olfersiirevealed from a Duvernoy's (venom) gland transcriptome

We investigated the putative toxins of Philodryas olfersii (Colubridae), a representative of a family of snakes neglected in venom studies despite their growing medical importance. Transcriptomic data of the venom gland complemented by proteomic analysis of the gland secretion revealed the presence of major toxin classes from the Viperidae family, including serine proteases, metalloproteases, C-type lectins, Crisps, and a C-type natriuretic peptide (CNP). Interestingly, the phylogenetic analysis of the CNP precursor showed it as a linker between two related precursors found in Viperidae and Elapidae snakes. We suggest that these precursors constitute a monophyletic group derived from the vertebrate CNPs.

Molecular Cloning of Echis ocellatus Disintegrins Reveals Non-Venom-Secreted Proteins and a Pathway for the Evolution of Ocellatusin

We report the cloning and sequence analysis of Echis ocellatus cDNAs coding for dimeric disintegrin subunits and for the short disintegrin ocellatusin. All the dimeric disintegrin subunit messengers belong to the short-coding class, indicating that short messengers may be more widely distributed than previously thought. Mass spectrometric analysis of the HPLC-separated venom proteins was performed to characterize the dimeric disintegrins expressed in the venom proteome. In addition to previously reported EO4 and EO5 heterodimers, a novel dimeric disintegrin containing RGD- and KGD-bearing subunits was identified. However, a WGD-containing polypeptide encoded by clone Eo1-1 was not detected in the venom, suggesting the occurrence of larger genomic than proteomic diversity, which could represent part of a non-venom-secreted reservoir of disintegrin that may eventually acquire physiological relevance for the snake upon changes of ecological niches and prey habits. On the other hand, the realization of the existence of two distinct messengers coding for the short disintegrin ocellatusin reveals key events of the evolutionary emergence of the short disintegrin ocellatusin from a short-coding dimeric disintegrin precursor by two nucleotide mutations.

Characterization of venom components from the scorpion Androctonus crassicauda of Turkey: Peptides and genes

The soluble venom from the scorpion Androctonus crassicauda was fractionated by high performance liquid chromatography. At least 44 different sub-fractions were resolved and collected for finger print mass analysis using an electrospray mass spectrometer. This analysis revealed the presence of 80 distinct molecular mass components, from which five were further characterized. A peptide, named Acra1 was fully sequenced. It contains 58 amino acid residues cross-bridged by six cysteines forming three disulfide pairs, with a molecular mass of 6497 Da. A second purified peptide named Acra2 was partially sequenced with a molecular mass of 7849 Da. Acra1 is toxic and Acra2 is lethal to mice, at the dose assayed. Additionally, a cDNA library of the venomous gland of one specimen was prepared and several clones were obtained among which is one that codes for Acra1. Three analog gene sequences were found with point mutations either in the section that corresponds to the mature peptide or to the signal peptide. The signal peptide is 22 amino acid residues long. Several other gene sequences obtained suggest the presence in this venom of three distinct groups of peptides, among which are peptides similar to known Na(+)-channel specific toxins of other scorpions. A new type of peptide was identified with odd number of cysteines (seven), allowing the formation of heterodimers with molecular masses in the range of 16,000 atomic mass units (a.m.u.).

Is it possible to develop pan-arthropod vaccines?

Hematophagous arthropods that transmit the etiological agents of arthropod-borne diseases have become the focus of anti-vector vaccines, targeted mainly at components of their saliva and midgut. These efforts have been directed mostly towards developing species-specific vaccines. An alternative is to target cross-reactive epitopes that have been preserved during evolution of the arthropods. The N- and O-linked glycans that are attached to arthropod glycoproteins are one of the potential targets of this pan-arthropod vaccine approach. Here, we discuss how genetically modified Drosophila melanogaster cells can be used to synthesize and to deliver these arthropod glycans to vertebrate hosts.

Molecular Cloning of Disintegrin-like Transcript BA-5A from a Bitis arietans Venom Gland cDNA Library: A Putative Intermediate in the Evolution of the Long-Chain Disintegrin Bitistatin

We report the cloning and sequence analysis of BA-5A from a venom gland cDNA library of the puff adder, Bitis arietans, that encodes a novel ECD-disintegrin-like domain. BA-5A is a unique PII disintegrin. It contains the 16 cysteine residues that are conserved in all known disintegrin-like domains of ADAM proteins and snake venom metalloproteinases but lacks the cysteine-rich domain. These features suggest that BA-5A may represent an intermediate in the evolutionary pathway of the long disintegrin bitistatin and that removal of the cysteine-rich domain and loss of the PIII-specific disulfide bond were separate events along the structural diversification pathway of disintegrins, the former predating the latter. The protein family composition of the Bitis arietans venom, as determined by combination of reversed-phase HPLC and proteomic analysis, was as follows: Zn(2+)-metalloproteinase (38.5%), serine proteinase (19.5%), disintegrin (17.8%), C-type lectin-like (13.2%), PLA(2) (4.3%), Kunitz-type inhibitor (4.1%), cystatin (1.7%), and unknown (0.9%). BA-5A could not be detected in the venom proteome of Bitis arietans. The occurrence of this very low-abundance (< 0.05%) or nonexpressed disintegrin transcript indicates a hitherto unrecognized structural diversity of this protein family. Whether BA-5A plays a physiological role or represents an orphan protein which could eventually evolve a role in the adaptation of snakes to changing ecological niches and prey habits deserves further investigation.

Transcriptome analysis of Deinagkistrodon acutus venomous gland focusing on cellular structure and functional aspects using expressed sequence tags

Background

The snake venom gland is a specialized organ, which synthesizes and secretes the complex and abundant toxin proteins. Though gene expression in the snake venom gland has been extensively studied, the focus has been on the components of the venom. As far as the molecular mechanism of toxin secretion and metabolism is concerned, we still knew a little. Therefore, a fundamental question being arisen is what genes are expressed in the snake venom glands besides many toxin components?

Results

To examine extensively the transcripts expressed in the venom gland of Deinagkistrodon acutus and unveil the potential of its products on cellular structure and functional aspects, we generated 8696 expressed sequence tags (ESTs) from a non-normalized cDNA library. All ESTs were clustered into 3416 clusters, of which 40.16% of total ESTs belong to recognized toxin-coding sequences; 39.85% are similar to cellular transcripts; and 20.00% have no significant similarity to any known sequences. By analyzing cellular functional transcripts, we found high expression of some venom related genes and gland-specific genes, such as calglandulin EF-hand protein gene and protein disulfide isomerase gene. The transcripts of creatine kinase and NADH dehydrogenase were also identified at high level. Moreover, abundant cellular structural proteins similar to mammalian muscle tissues were also identified. The phylogenetic analysis of two snake venom toxin families of group III metalloproteinase and serine protease in suborder Colubroidea showed an early single recruitment event in the viperids evolutionary process.

Conclusion

Gene cataloguing and profiling of the venom gland of Deinagkistrodon acutus is an essential requisite to provide molecular reagents for functional genomic studies needed for elucidating mechanisms of action of toxins and surveying physiological events taking place in the very specialized secretory tissue. So this study provides a first global view of the genetic programs for the venom gland of Deinagkistrodon acutus described so far and an insight into molecular mechanism of toxin secreting.

All sequences data reported in this paper have been submitted into the public database [GenBank: DV556511-DV565206].

Snake venom glutaminyl cyclase

Glutaminyl cyclase (QC) catalyzes N-terminal glutamine cyclization of many endocrine peptides and is typically abundant in brain tissue. As three-finger toxins in the venoms of colubrid snakes Boiga dendrophila and Boiga irregularis contain N-terminal pyroglutamate, we searched for QC in venom glands of both snakes. Here we report cDNA sequences of QC from brain and venom gland tissues of Boiga species. We propose that QC expressed in snake venom gland tissue plays a role in the N-terminal pyroglutamate formation of several snake venom toxins, indirectly contributing to venom potency.

Molecular cloning of disintegrins from Cerastes vipera and Macrovipera lebetina transmediterranea venom gland cDNA libraries: insight into the evolution of the snake venom integrin-inhibition system

We report the cloning and sequence analysis of Cerastes vipera and Macrovipera lebetina transmediterranea cDNAs coding for short non-RGD (Arg-Gly-Asp) disintegrins and for dimeric disintegrin subunits. The mRNAs belong to the short-coding class, suggesting that these disintegrin mRNAs may be more widely distributed than previously thought. Our data also argue for a common ancestry of the mRNAs of short disintegrins and those coding for precursors of dimeric disintegrin chains. The Macrovipera lebetina transmediterranea dimeric disintegrin reported to inhibit the laminin-binding integrins alpha3beta1, alpha6beta1 and alpha7beta1 was analysed using a proteomic approach and was shown to bear MLD (Met-Leu-Asp) and VGD (Val-Gly-Asp) motifs. The results highlight the fact that disintegrins have evolved a restricted panel of integrin-blocking sequences that segregate with defined branches of the phylogenetic tree of the integrin alpha-chains, providing novel insights into the evolutionary adaptation of the snake venom antagonists to the ligand-binding sites of their target integrin receptors.

Lachesis muta (Viperidae) cDNAs reveal diverging pit viper molecules and scaffolds typical of cobra (Elapidae) venoms: implications for snake toxin repertoire evolution

Efforts to describe toxins from the two major families of venomous snakes (Viperidae and Elapidae) usually reveal proteins belonging to few structural types, particular of each family. Here we carried on an effort to determine uncommon cDNAs that represent possible new toxins from Lachesis muta (Viperidae). In addition to nine classes of typical toxins, atypical molecules never observed in the hundreds of Viperidae snakes studied so far are highly expressed: a diverging C-type lectin that is related to Viperidae toxins but appears to be independently originated; an ohanin-like toxin, which would be the third member of the most recently described class of Elapidae toxins, related to human butyrophilin and B30.2 proteins; and a 3FTx-like toxin, a new member of the widely studied three-finger family of proteins, which includes major Elapidae neurotoxins and CD59 antigen. The presence of these common and uncommon molecules suggests that the repertoire of toxins could be more conserved between families than has been considered, and their features indicate a dynamic process of venom evolution through molecular mechanisms, such as multiple recruitments of important scaffolds and domain exchange between paralogs, always keeping a minimalist nature in most toxin structures in opposition to their nontoxin counterparts.

Snake venom metalloproteases--structure and function of catalytic and disintegrin domains

Snake venoms are relevant sources of toxins that have evolved towards the engineering of highly active compounds. In the last years, research efforts have produced great advance in their understanding and uses. Metalloproteases with disintegrin domains are among the most abundant toxins in many Viperidae snake venoms. This review will focus on the structure, function and possible applications of the metalloprotease and disintegrin domains.

Ohanin, a novel protein from king cobra venom: its cDNA and genomic organization

Ohanin, from king cobra venom, is a novel protein which induces hypolocomotion and hyperalgesia in mice [Pung, Y.F., Wong, P.T.H., Kumar, P.P., Hodgson W.C., Kini, R.M., 2005. Ohanin, a novel protein from king cobra venom induces hypolocomotion and hyperalgesia in mice. J. Biol. Chem. 280, 13137-13147.]. It is weakly similar to PRY-SPRY domains (B30.2-like domain). Here we report the complete cDNA and genomic organization of ohanin. Interestingly, cDNA sequence does not show significant sequence similarity to any known sequences, including those of B30.2-like domain-containing proteins. Its full-length cDNA sequence of 1558 bp encodes for prepro-ohanin with a propeptide segment at the C-terminal. Ohanin is the first member of a new subfamily of proteins containing B30.2-like domain with short N-terminal segment. We named this subfamily as vespryns. There are two mRNA subtypes differing in their 5'-untranslated regions. Southern hybridization study shows that ohanin is encoded by a single gene. Its genomic sequence is 7086 bp with five exons and four introns, and the two types of mRNAs are generated by alternative splicing of exon 2. Our results indicate that ohanin and vespryns may have evolved from the same ancestral gene as B30.2 domain.

Snake fangs from the Lower Miocene of Germany: evolutionary stability of perfect weapons

There is a general consensus that most of today's nonvenomous snakes are descendants of venomous snakes that lost their venomous capabilities secondarily. This implies that the evolutionary history of venomous snakes and their venom apparatus should be older than the current evidence from the fossil record. We compared some of the oldest-known fossil snake fangs from the Lower Miocene of Germany with those of modern viperids and elapids and found their morphology to be indistinguishable from the modern forms. The primary function of recent elapid and viperid snake fangs is to facilitate the extremely rapid, stab-like application of highly toxic venoms. Our findings therefore indicate that the other components of the venom-delivery system of Early Miocene vipers and elapids were also highly developed, and that these snakes used their venom in the same way as their modern relatives. Thus, the fossil record supports the view that snakes used their venoms to rapidly subdue prey long before the mid-Tertiary onset of the global environmental changes that seem to have supported the successful radiation of venomous snakes.

Genomic organization of four novel nondisulfide-bridged peptides from scorpion Mesobuthus martensii Karsch: gaining insight into evolutionary mechanism

At least 25 nondisulfide-bridged peptides (NDBPs) have been identified and characterized from scorpions. However, the genomic organization of the genes that encode these peptides have not been reported yet. BmKa1, BmKa2 and BmKb1 are three novel genes that code for NDBPs identified by our group from Mesobuthus martensii Karsch. Based on their cDNA sequences, the genomic DNA sequences encoding these peptides were obtained using the PCR method. Sequence analysis showed that three distinct genomic structural patterns are used to encode these three peptides. The BmKa1 gene is not interrupted by any introns. However, the BmKa2 gene is composed of two exons, interrupted by a 67 bp intron that is located in the DNA region encoding the mature peptide. Two genomic homologues of the BmKb1 cDNA sequence, named BmKb1' and BmKb2, respectively, were obtained. The BmKb1' gene contains one intron of 593 bp, inserted into the DNA region that encodes the signal peptide. Similarly, the BmKb2 gene also contains an intron that interrupts the exon that encodes the NDBP signal peptide. The amino acid sequences deduced for BmKb2 and BmKb1' differ only at one position. The data suggest that the genomic organizational pattern of NDBPs displays more divergence than that exhibited by the genes that encode disulfide-bridged peptides from scorpions.

Early evolution of the venom system in lizards and snakes

Among extant reptiles only two lineages are known to have evolved venom delivery systems, the advanced snakes and helodermatid lizards (Gila Monster and Beaded Lizard). Evolution of the venom system is thought to underlie the impressive radiation of the advanced snakes (2,500 of 3,000 snake species). In contrast, the lizard venom system is thought to be restricted to just two species and to have evolved independently from the snake venom system. Here we report the presence of venom toxins in two additional lizard lineages (Monitor Lizards and Iguania) and show that all lineages possessing toxin-secreting oral glands form a clade, demonstrating a single early origin of the venom system in lizards and snakes. Construction of gland complementary-DNA libraries and phylogenetic analysis of transcripts revealed that nine toxin types are shared between lizards and snakes. Toxinological analyses of venom components from the Lace Monitor Varanus varius showed potent effects on blood pressure and clotting ability, bioactivities associated with a rapid loss of consciousness and extensive bleeding in prey. The iguanian lizard Pogona barbata retains characteristics of the ancestral venom system, namely serial, lobular non-compound venom-secreting glands on both the upper and lower jaws, whereas the advanced snakes and anguimorph lizards (including Monitor Lizards, Gila Monster and Beaded Lizard) have more derived venom systems characterized by the loss of the mandibular (lower) or maxillary (upper) glands. Demonstration that the snakes, iguanians and anguimorphs form a single clade provides overwhelming support for a single, early origin of the venom system in lizards and snakes. These results provide new insights into the evolution of the venom system in squamate reptiles and open new avenues for biomedical research and drug design using hitherto unexplored venom proteins.

Overview of scorpion toxins specific for Na+ channels and related peptides: biodiversity, structure-function relationships and evolution

Scorpion venoms contain a large number of bioactive components. Several of the long-chain peptides were shown to be responsible for neurotoxic effects, due to their ability to recognize Na(+) channels and to cause impairment of channel functions. Here, we revisited the basic paradigms in the study of these peptides in the light of recent data concerning their structure-function relationships, their functional divergence and extant biodiversity. The reviewed topics include: the criteria for classification of long-chain peptides according to their function, and a revision of the state-of-the-art knowledge concerning the surface areas of contact of these peptides with known Na(+) channels. Additionally, we compiled a comprehensive list encompassing 191 different amino acid sequences from long-chain peptides purified from scorpion venoms. With this dataset, a phylogenetic tree was constructed and discussed taking into consideration their documented functional divergence. A critical view on problems associated with the study of these scorpion peptides is presented, drawing special attention to the points that need revision and to the subjects under intensive research at this moment, regarding scorpion toxins specific for Na(+) channels and the other related long-chain peptides recently described.

The transcriptome of the salivary glands of the female western black-legged tick Ixodes pacificus (Acari: Ixodidae)

Sequencing of an Ixodes pacificus salivary gland cDNA library yielded 1068 sequences with an average undetermined nucleotide of 1.9% and an average length of 487 base pairs. Assembly of the expressed sequence tags yielded 557 contigs, 138 of which appear to code for secreted peptides or proteins based on translation of a putative signal peptide. Based on the BLASTX similarity of these contigs to 66 matches of Ixodes scapularis peptide sequences, only 58% sequence identity was found, indicating a rapid divergence of salivary proteins as observed previously for mosquito and triatomine bug salivary proteins. Here we report 106 mostly full-length sequences that clustered in 16 different families: Basic-tail proteins rich in lysine in the carboxy-terminal, Kunitz-containing proteins (monolaris, ixolaris and penthalaris families), proline-rich peptides, 5-, 9.4- and 18.7-kDa proteins of unknown functions, in addition to metalloproteases (class PIII-like) similar to reprolysins. We also have found a family of disintegrins, named ixodegrins that display homology to variabilin, a GPIIb/IIIa antagonist from the tick Dermacentor variabilis. In addition, we describe peptides (here named ixostatins) that display remarkable similarities to the cysteine-rich domain of ADAMST-4 (aggrecanase). Many molecules were assigned in the lipocalin family (histamine-binding proteins); others appear to be involved in oxidant metabolism, and still others were similar to ixodid proteins such as the anticomplement ISAC. We also identified for the first time a neuropeptide-like protein (nlp-31) with GGY repeats that may have antimicrobial activity. In addition, 16 novel proteins without significant similarities to other tick proteins and 37 housekeeping proteins that may be useful for phylogenetic studies are described. Some of these proteins may be useful for studying vascular biology or the immune system, for vaccine development, or as immunoreagents to detect prior exposure to ticks. Electronic version of the manuscript can be found at.