Evolutionary analysis of novel serine proteases in the venom gland transcriptome of Bitis gabonica rhinoceros

Background: Serine proteases are major components of viper venom and target various stages of the blood coagulation system in victims and prey. A better understanding of the diversity of serine proteases and other enzymes present in snake venom will help to understand how the complexity of snake venom has evolved and will aid the development of novel therapeutics for treating snake bites. Methodology and Principal Findings: Four serine protease-encoding genes from the venom gland transcriptome of Bitis gabonica rhinoceros were amplified and sequenced. Mass spectrometry suggests the four enzymes corresponding to these genes are present in the venom of B. g. rhinoceros. Two of the enzymes, rhinocerases 2 and 3 have substitutions to two of the serine protease catalytic triad residues and are thus unlikely to be catalytically active, though they may have evolved other toxic functions. The other two enzymes, rhinocerases 4 and 5, have classical serine protease catalytic triad residues and thus are likely to be catalytically active, however they have glycine rather than the more typical aspartic acid at the base of the primary specificity pocket (position 189). Based on a detailed analysis of these sequences we suggest that alternative splicing together with individual amino acid mutations may have been involved in their evolution. Changes within amino acid segments which were previously proposed to undergo accelerated change in venom serine proteases have also been observed. Conclusions and Significance: Our study provides further insight into the diversity of serine protease isoforms present within snake venom and discusses their possible functions and how they may have evolved. These multiple serine protease isoforms with different substrate specificities may enhance the envenomation effects and help the snake to adapt to new habitats and diets. Our findings have potential for helping the future development of improved therapeutics for snake bites.

Analysis of Bothrops jararacussu venomous gland transcriptome focusing on structural and functional aspects: I - gene expression profile of highly expressed phospholipases A(2)

Snake venom glands are a rich source of bioactive molecules such as peptides, proteins and enzymes that show important pharmacological activity leading to in local and systemic effects as pain, edema, bleeding and muscle necrosis. Most studies on pharmacologically active peptides and proteins from snake venoms have been concerned with isolation and structure elucidation through methods of classical biochemistry. As an attempt to examine the transcripts expressed in the venom gland of Bothrops jararacussu and to unveil the toxicological and pharmacological potential of its products at the molecular level, we generated 549 expressed sequence tags (ESTs) from a directional cDNA library. Sequences obtained from single-pass sequencing of randomly selected cDNA clones could be identified by similarities searches on existing databases, resulting in 197 sequences with significant similarity to phospholipase A(2) (PLA(2)), of which 83.2% were Lys49-PLA(2) homologs (BOJU-1), 0.1% were basic Asp49-PLA(2)s (BOJU-II) and 0.6% were acidic Asp49-PLA(2)s (BOJU-III). Adjoining this very abundant class of proteins we found 88 transcripts codifying for putative sequences of metalloproteases, which after clustering and assembling resulted in three full-length sequences: BOJUMET-I, BOJUMET-II and BOJUMET-III; as well as 25 transcripts related to C-type lectin like protein including a full-length cDNA of a putative galactose binding C-type lectin and a cluster of eight serine-proteases transcripts including a full-length cDNA of a putative serine protease. Among the full-length sequenced clones we identified a nerve growth factor (Bj-NGF) with 92% identity with a human NGF (NGHUBM) and an acidic phospholipase A2 (BthA-I-PLA(2)) displaying 85-93% identity with other snake venom toxins. Genetic distance among PLA(2)s from Bothrops species were evaluated by phylogenetic analysis. Furthermore, analysis of full-length putative Lys49-PLA(2) through molecular modeling showed conserved structural domains, allowing the characterization of those proteins as group II PLA(2)s. The constructed cDNA library provides molecular clones harboring sequences that can be used to probe directly the genetic material from gland venom of other snake species. Expression of complete cDNAs or their modified derivatives will be useful for elucidation of the structure-function relationships of these toxins and peptides of biotechnological interest. (C) 2004 Elsevier SAS. All rights reserved.

Elements of the granular gland peptidome and transcriptome persist in air-dried skin of the South American orange-legged leaf frog, Phyllomedusa hypocondrialis.

The defensive strategy of amphibians against predator attack relies heavily on the secretion of noxious/toxic chemical cocktails from specialized skin granular glands. Bioactive peptides constitute a major component of secretions in many species and the most complex are produced by neotropical leaf frogs of the sub-family Phyllomedusinae. We recently reported that these skin secretions contain elements of both the granular gland peptidome and transcriptome and that polyadenylated mRNAs constituting the latter are protected from degradation by interactions with endogenous amphipathic peptides. This thus permits parallel amino acid sequencing of peptides and nucleic acid sequencing of cloned precursor transcripts from single lyophilized samples of secretion. Here we report that the protection afforded is sufficiently robust to permit transcriptome studies by cloning of full-length polyadenylated peptide precursor encoding mRNAs from libraries constructed using ambient temperature air-dried skin from recently deceased specimens as source material. The technique was sufficiently sensitive to permit the identification of cDNAs encoding antimicrobial peptides constituted by six different isoforms of phylloseptin and two dermaseptins. Also, for the first time, establishment of the nucleic acid and amino acid sequence of the precursor encoding the phyllomedusine frog skin bradykinin-related peptide, phyllokinin, from cloned cDNA, was achieved. These data unequivocally demonstrate that the granular gland transcriptome persists in air-dried amphibian skin—a finding that may have fundamental implications in the study of archived materials but also in the wider field of molecular biology.

Combined snake venomics and venom gland transcriptomic analysis of Bothropoides pauloensis

Unraveling the repertoire of venom toxins of Bothropoides pauloensis was assessed by snake venomics and venom gland transcriptomic surveys. Both approaches yielded converging overall figures, pointing to metalloproteinases (similar to 37%), PLA(2)s (26-32%), and vasoactive (bradykinin-potentiating) peptides (12-17%) as the major toxin classes. The high occurrence of SVMPs, PLA(2) molecules, vasoactive peptides, along with serine proteinases, explains the local and systemic effects observed in envenomations by B. pauloensis. Minor (<3%) C-type lectin, serine proteinase, L-amino acid oxidase, nerve growth factor, and CRISP molecules were also identified in the transcriptome and the proteome. Low abundance (0.3%) EST singletons coding for vascular endothelial growth factor (svVEGF), ohanin, hyaluronidase, and 5' nucleotidase were found only in the venom gland cDNA library. At the molecular level, the transcriptomic and proteomic datasets display low compositional concordance. In particular, although there is good agreement between transcriptome and proteome in the identity of BPPs, PLA(2) molecules and L-amino acid oxidase, both datasets strongly depart in their C-type lectin and SVMP complements. These data support the view that venom composition is influenced by transcriptional and translational mechanisms and emphasize the value of combining proteomic and transcriptomic approaches to acquire a more complete understanding of the toxinological profile and natural history of the snake venom. (C) 2012 Elsevier B.V. All rights reserved.

Cloning from tissue surrogates: antimicrobial peptide (esculentin) cDNAs from the defensive skin secretions of Chinese ranid frogs

The defensive skin secretions of amphibians are a rich source of bioactive peptides. Here we describe a rapid technique for skin granular gland transcriptome cloning from a surrogate tissue-the secretion itself. cDNA libraries were constructed from lyophilized skin secretion from each of the Chinese frogs (Rana schmackeri, Rana versabilis, and Rana plancyi fukienensis) using magnetic oligo(dT) bead-captured polyadenylated mRNA as templates. Specific esculentin cDNAs were amplified by 3'-RACE using a degenerate primer designed for a consensus nucleotide sequence in the 5' untranslated region of previously characterized ranid frog peptide cDNAs. The cloned cDNAs were found to encode the antimicrobial peptides esculentins 1 and 2 from each of the species examined. The presence of predicted peptide structures in skin secretions was confirmed by MALDI-TOF mass spectrometry and automated Edman degradation. This experimental approach can thus rapidly expedite parallel transcriptome and peptidome analysis of amphibian granular gland secretions without harming or sacrificing donor animals.

Nigrocin-2 peptides from Chinese Odorrana frogs – integration of UPLC/MS/MS with molecular cloning in amphibian skin peptidome analysis

Peptidomics is a powerful set of tools for the identification, structural elucidation and discovery of novel regulatory peptides and for monitoring the degradation pathways of structurally and catalytically important proteins. Amphibian skin secretions, arising from specialized granular glands, often contain complex peptidomes containing many components of entirely novel structure and unique site-substituted analogues of known peptide families. Following the discovery that the granular gland transcriptome is present in such secretions in a PCR-amenable form, we designed a strategy for peptide structural characterization involving the integration of ‘shotgun’ cloning of cDNAs encoding peptide precursors, deduction of putative bioactive peptide structures, and confirmation of these structures using tandem MS/MS sequencing. Here, we illustrate this strategy by means of elucidation of the primary structures of nigrocin-2 homologues from the defensive skin secretions of four species of Chinese Odorrana frogs, O. schmackeri, O. livida, O. hejiangensis and O. versabilis. Synthetic replicates of the peptides were found to possess antimicrobial activity. Nigrocin-2 peptides occur widely in the skin secretions of Asian ranid frogs and in those of the Odorrana group, and are particularly well-represented and of diverse structure in some species. Integration of the molecular analytical technologies described provides a means for rapid structural characterization of novel peptides from complex natural libraries in the absence of systematic online database information.

Construction of cDNA libraries from trifluoroacetic acid-solvated amphibian skin secretions: molecular cloning of multiple bombinin-like peptide precursor transcripts from a library of yellow-bellied toad (Bombina variegata) secretion

Amphibian skin secretions are renowned as complex mixtures of bioactive peptides many of which are analogues of endogenous regulatory peptides. While skin secretions can be obtained non-invasively for peptidome analysis, parallel studies on the granular gland transcriptome required specimen sacrifice. The aim of the present study was to analyse archived skin secretions to determine the robustness of bioactive peptide precursor-encoding polyadenylated mRNAs in an attempt to extract maximum molecular information from rare samples. A range of solvated skin secretion samples were examined after lyophilisation for their potential to generate viable and comprehensive cDNA libraries based upon polyadenylated mRNA capture and amplification/cloning using appropriate commercial kits. Here we present unequivocal data that the granular gland transcriptome persists in a PCR amenable format even after storage for as long as 12 years in 0.1%(v/v) aqueous trifluoroacetic acid (TFA). We used a pooled skin secretion sample (2 ml) from the yellow-bellied toad, Bombina variegata (n = 14), containing the equivalent of 5 mg/ml of lyophilised skin secretion, that had been used in part for peptide isolation purposes in 1998 and had been stored at - 20 °C since that time. In the first cloning experiment, 12 different bombinin-like peptide precursor cDNAs were cloned encoding 17 different bombinins, the majority of which were novel. Subsequently, bombesin and bradykinin-related peptide precursor transcripts have been cloned successfully. These data illustrate the unexpected stability/longevity of the transcriptome in these secretions — a finding with implications for both this field of research and for the wider field of molecular biology.

Antimicrobial peptides and alytesin are co-secreted from the venom of the Midwife toad, Alytes maurus (Alytidae, Anura): Implications for the evolution of frog skin defensive secretions

The skin secretions of frogs and toads (Anura) have long been a known source of a vast abundance of bioactive substances. In the past decade, transcriptome data of the granular glands of anuran skin has given new impetus to investigations of the putative constituent peptides. Alytes obstetricans was recently investigated and novel peptides with antimicrobial activity were isolated and functionally characterised. However, genetic data for the evolutionarily ancient lineage to which Alytes belongs (midwife toads; Alytidae) remains unavailable.

Here we present the first such genetic data for Alytidae, derived via the granular gland transcriptome of a closely-related species of midwife toad, Alytes maurus. First, we present nucleotide sequences of the entire peptide precursors for four novel antimicrobial peptides (AMPs). The two precursors resemble those from Bombinatoridae in both their structural architecture and amino acid sequence. Each precursor comprises two AMPs as tandem repeats, with a member of the alyteserin-1 family (alyteserin-1Ma: GFKEVLKADLGSLVKGIAAHVAN-NH2 or alyteserin-1Mb: GFKEVLKAGLGSLVKGIPAHVAN-NH2) followed by its corresponding member from the alyteserin-2 family (alyteserin-2Ma: FIGKLISAASGLLSHL-NH2 or alyteserin-2Mb: ILGAIIPLVSGLLSHL-NH2). Synthetic replicates of the four AMPs possessed minimal inhibitory concentrations (MICs) ranging from 9.5 to 300 µM, with the most potent being alyteserin-2Ma. Second, we also cloned the cDNA encoding an alytesin precursor, with the active alytesin exhibiting high sequence identity to bombesin-related peptides from other frogs. All putative mature peptide sequences were confirmed to be present in the skin secretion via LC/MS.

The close structural resemblance of the alyteserin genes that we isolated for A. maurus with those of Bombina provide independent molecular evidence for a close evolutionary relationship between these genera as well as more support for the convergent evolution of the AMP system within anurans. In contrast to the more evolutionarily conserved nature of neuropeptides (including alytesin, which we also isolated), the more variable nature of the AMP system together with the sporadic distribution of AMPs among anuran amphibians fuels in part our hypothesis that the latter system was co-opted secondarily to fulfil a function in the innate immune system, having originally evolved for defence against potential macropredators.

Snake Venomic of Crotalus durissus terrificus-Correlation with Pharmacological Activities

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Venomics Profiling of Thamnodynastes strigatus Unveils Matrix Metalloproteinases and Other Novel Proteins Recruited to the Toxin Arsenal of Rear-Fanged Snakes

Rear-fanged and aglyphous snakes are usually considered not dangerous to humans because of their limited capacity of injecting venom. Therefore, only a few studies have been dedicated to characterizing the venom of the largest parcel of snake fauna. Here, we investigated the venom proteome of the rear-fanged snake Thamnodynastes strigatus, in combination with a transcriptomic evaluation of the venom gland. About 60% of all transcripts code for putative venom components. A striking finding is that the most abundant type of transcript (similar to 47%) and also the major protein type in the venom correspond to a new kind of matrix metalloproteinase (MMP) that is unrelated to the classical snake venom metalloproteinases found in all snake families. These enzymes were recently suggested as possible venom components, and we show here that they are proteolytically active and probably recruited to venom from a MMP-9 ancestor. Other unusual proteins were suggested to be venom components: a protein related to lactadherin and an EGF repeat-containing transcript. Despite these unusual molecules, seven toxin classes commonly found in typical venomous snakes are also present in the venom. These results support the evidence that the arsenals of these snakes are very diverse and harbor new types of biologically important molecules.

Isolation, N-glycosylations and Function of a Hyaluronidase-Like Enzyme from the Venom of the Spider Cupiennius salei.

STRUCTURE OF CUPIENNIUS SALEI VENOM HYALURONIDASE Hyaluronidases are important venom components acting as spreading factor of toxic compounds. In several studies this spreading effect was tested on vertebrate tissue. However, data about the spreading activity on invertebrates, the main prey organisms of spiders, are lacking. Here, a hyaluronidase-like enzyme was isolated from the venom of the spider Cupiennius salei. The amino acid sequence of the enzyme was determined by cDNA analysis of the venom gland transcriptome and confirmed by protein analysis. Two complex N-linked glycans akin to honey bee hyaluronidase glycosylations, were identified by tandem mass spectrometry. A C-terminal EGF-like domain was identified in spider hyaluronidase using InterPro. The spider hyaluronidase-like enzyme showed maximal activity at acidic pH, between 40-60°C, and 0.2 M KCl. Divalent ions did not enhance HA degradation activity, indicating that they are not recruited for catalysis. FUNCTION OF VENOM HYALURONIDASES Besides hyaluronan, the enzyme degrades chondroitin sulfate A, whereas heparan sulfate and dermatan sulfate are not affected. The end products of hyaluronan degradation are tetramers, whereas chondroitin sulfate A is mainly degraded to hexamers. Identification of terminal N-acetylglucosamine or N-acetylgalactosamine at the reducing end of the oligomers identified the enzyme as an endo-β-N-acetyl-D-hexosaminidase hydrolase. The spreading effect of the hyaluronidase-like enzyme on invertebrate tissue was studied by coinjection of the enzyme with the Cupiennius salei main neurotoxin CsTx-1 into Drosophila flies. The enzyme significantly enhances the neurotoxic activity of CsTx-1. Comparative substrate degradation tests with hyaluronan, chondroitin sulfate A, dermatan sulfate, and heparan sulfate with venoms from 39 spider species from 21 families identified some spider families (Atypidae, Eresidae, Araneidae and Nephilidae) without activity of hyaluronidase-like enzymes. This is interpreted as a loss of this enzyme and fits quite well the current phylogenetic idea on a more isolated position of these families and can perhaps be explained by specialized prey catching techniques.

Sexual selection in true fruit flies (Diptera: Tephritidae) : transcriptome and experimental evidences for phytochemicals increasing male competitive ability

In male tephritid fruit flies of the genus Bactrocera, feeding on secondary plant compounds (sensu lato male lures = methyl eugenol, raspberry ketone and zingerone) increases male mating success. Ingested male lures alter the male pheromonal blend, normally making it more attractive to females and this is considered the primary mechanism for the enhanced mating success. However, the male lures raspberry ketone and zingerone are known, across a diverse range of other organisms, to be involved in increasing energy metabolism. If this also occurs in Bactrocera, then this may represent an additional benefit to males as courtship is metabolically expensive and lure feeding may increase a fly's short-term energy. We tested this hypothesis by performing comparative RNA-seq analysis between zingerone-fed and unfed males of Bactrocera tryoni. We also carried out behavioural assays with zingerone- and cuelure-fed males to test whether they became more active. RNA-seq analysis revealed, in zingerone-fed flies, up-regulation of 3183 genes with homologues transcripts to those known to regulate intermale aggression, pheromone synthesis, mating and accessory gland proteins, along with significant enrichment of several energy metabolic pathways and gene ontology terms. Behavioural assays show significant increases in locomotor activity, weight reduction and successful mating after mounting; all direct/indirect measures of increased activity. These results suggest that feeding on lures leads to complex physiological changes, which result in more competitive males. These results do not negate the pheromone effect, but do strongly suggest that the phytochemical-induced sexual selection is governed by both female preference and male competitive mechanisms.

Sexual selection in true fruit flies (Diptera: Tephritidae): transcriptome and experimental evidences for phytochemicals increasing male competitive ability

In male tephritid fruit flies of the genus Bactrocera, feeding on secondary plant compounds (sensu lato male lures = methyl eugenol, raspberry ketone and zingerone) increases male mating success. Ingested male lures alter the male pheromonal blend, normally making it more attractive to females and this is considered the primary mechanism for the enhanced mating success. However, the male lures raspberry ketone and zingerone are known, across a diverse range of other organisms, to be involved in increasing energy metabolism. If this also occurs in Bactrocera, then this may represent an additional benefit to males as courtship is metabolically expensive and lure feeding may increase a fly's short-term energy. We tested this hypothesis by performing comparative RNA-seq analysis between zingerone-fed and unfed males of Bactrocera tryoni. We also carried out behavioural assays with zingerone- and cuelure-fed males to test whether they became more active. RNA-seq analysis revealed, in zingerone-fed flies, up-regulation of 3183 genes with homologues transcripts to those known to regulate intermale aggression, pheromone synthesis, mating and accessory gland proteins, along with significant enrichment of several energy metabolic pathways and gene ontology terms. Behavioural assays show significant increases in locomotor activity, weight reduction and successful mating after mounting; all direct/indirect measures of increased activity. These results suggest that feeding on lures leads to complex physiological changes, which result in more competitive males. These results do not negate the pheromone effect, but do strongly suggest that the phytochemical-induced sexual selection is governed by both female preference and male competitive mechanisms.

Granular gland transcriptomes in stimulated amphibian skin secretions.

Amphibian defensive skin secretions are complex, species-specific cocktails of biologically active molecules, including many uncharacterized peptides. The study of such secretions for novel peptide discovery is time-limited, as amphibians are in rapid global decline. While secretion proteome analysis is non-lethal, transcriptome analysis has until now required killing of specimens prior to skin dissection for cDNA library construction. Here we present the discovery that polyadenylated mRNAs encoding dermal granular gland peptides are present in defensive skin secretions, stabilized by endogenous nucleic acid-binding amphipathic peptides. Thus parallel secretory proteome and transcriptome analyses can be performed without killing the specimen in this model amphibian system--a finding that has important implications in conservation of biodiversity within this threatened vertebrate taxon and whose mechanistics may have broader implications in biomolecular science.