Venom Composition of Neglected Bothropoid Snakes from the Amazon Rainforest: Ecological and Toxinological Implications

Snake venoms have evolved in several families of Caenophidae, and their toxins have been assumed to be biochemical weapons with a role as a trophic adaptation. However, it remains unclear how venom contributes to the success of venomous species for adaptation to different environments. Here we compared the venoms from Bothrocophias hyoprora, Bothrops taeniatus, Bothrops bilineatus smaragdinus, Bothrops brazili, and Bothrops atrox collected in the Amazon Rainforest, aiming to understand the ecological and toxinological consequences of venom composition. Transcriptomic and proteomic analyses indicated that the venoms presented the same toxin groups characteristic from bothropoids, but with distinct isoforms with variable qualitative and quantitative abundances, contributing to distinct enzymatic and toxic effects. Despite the particularities of each venom, commercial Bothrops antivenom recognized the venom components and neutralized the lethality of all species. No clear features could be observed between venoms from arboreal and terrestrial habitats, nor in the dispersion of the species throughout the Amazon habitats, supporting the notion that venom composition may not shape the ecological or toxinological characteristics of these snake species and that other factors influence their foraging or dispersal in different ecological niches.

ToxCodAn-Genome: an automated pipeline for toxin-gene annotation in genome assembly of venomous lineages

BACKGROUND

The rapid development of sequencing technologies resulted in a wide expansion of genomics studies using venomous lineages. This facilitated research focusing on understanding the evolution of adaptive traits and the search for novel compounds that can be applied in agriculture and medicine. However, the toxin annotation of genomes is a laborious and time-consuming task, and no consensus pipeline is currently available. No computational tool currently exists to address the challenges specific to toxin annotation and to ensure the reproducibility of the process.

RESULTS

Here, we present ToxCodAn-Genome, the first software designed to perform automated toxin annotation in genomes of venomous lineages. This pipeline was designed to retrieve the full-length coding sequences of toxins and to allow the detection of novel truncated paralogs and pseudogenes. We tested ToxCodAn-Genome using 12 genomes of venomous lineages and achieved high performance on recovering their current toxin annotations. This tool can be easily customized to allow improvements in the final toxin annotation set and can be expanded to virtually any venomous lineage. ToxCodAn-Genome is fast, allowing it to run on any personal computer, but it can also be executed in multicore mode, taking advantage of large high-performance servers. In addition, we provide a guide to direct future research in the venomics field to ensure a confident toxin annotation in the genome being studied. As a case study, we sequenced and annotated the toxin repertoire of Bothrops alternatus, which may facilitate future evolutionary and biomedical studies using vipers as models.

CONCLUSIONS

ToxCodAn-Genome is suitable to perform toxin annotation in the genome of venomous species and may help to improve the reproducibility of further studies. ToxCodAn-Genome and the guide are freely available at https://github.com/pedronachtigall/ToxCodAn-Genome.

Comparing morphological and secretory aspects of cephalic glands among the New World coral snakes brings novel insights on their biological roles

Oral and other cephalic glands have been surveyed by several studies with distinct purposes. Despite the wide diversity and medical relevance of the New World coral snakes, studies focusing on understanding the biological roles of the glands within this group are still scarce. Specifically, the venom glands of some coral snakes were previously investigated but all other cephalic glands remain uncharacterized. In this sense, performing morphological and molecular analysis of these glands may help better understand their biological role. Here, we studied the morphology of the venom, infralabial, rictal, and harderian glands of thirteen species of Micrurus and Micruroides euryxanthus. We also performed a molecular characterization of these glands from selected species of Micrurus using transcriptomic and proteomic approaches. We described substantial morphological variation in the cephalic glands of New World coral snakes and structural evidence for protein-secreting cells in the inferior rictal glands. Our molecular analysis revealed that the venom glands, as expected, are majorly devoted to toxin production, however, the infralabial and inferior rictal glands also expressed some toxin genes at low to medium levels, despite the marked morphological differences. On the other hand, the harderian glands were dominated by the expression of lipocalins, but do not produce toxins. Our integrative analysis, including the prediction of biological processes and pathways, helped decipher some important traits of cephalic glands and better understand their biology.

The Venom Composition of the Snake Tribe Philodryadini: 'Omic' Techniques Reveal Intergeneric Variability among South American Racers

Snakes of the Philodryadini tribe are included in the Dipsadidae family, which is a diverse group of rear-fanged snakes widespread in different ecological conditions, including habitats and diet. However, little is known about the composition and effects of their venoms despite their relevance for understanding the evolution of these snakes or even their impact on the occasional cases of human envenoming. In this study, we integrated venom gland transcriptomics, venom proteomics and functional assays to characterize the venoms from eight species of the Philodryadini tribe, which includes the genus Philodryas, Chlorosoma and Xenoxybelis. The most abundant components identified in the venoms were snake venom metalloproteinases (SVMPs), cysteine-rich secretory proteins (CRISPs), C-type lectins (CTLs), snake endogenous matrix metalloproteinases type 9 (seMMP-9) and snake venom serinoproteinases (SVSPs). These protein families showed a variable expression profile in each genus. SVMPs were the most abundant components in Philodryas, while seMMP-9 and CRISPs were the most expressed in Chlorosoma and Xenoxybelis, respectively. Lineage-specific differences in venom composition were also observed among Philodryas species, whereas P. olfersii presented the highest amount of SVSPs and P. agassizii was the only species to express significant amounts of 3FTx. The variability observed in venom composition was confirmed by the venom functional assays. Philodryas species presented the highest SVMP activity, whereas Chlorosoma species showed higher levels of gelatin activity, which may correlate to the seMMP-9 enzymes. The variability observed in the composition of these venoms may be related to the tribe phylogeny and influenced by their diets. In the presented study, we expanded the set of venomics studies of the Philodryadini tribe, which paves new roads for further studies on the evolution and ecology of Dipsadidae snakes.

Independent recruitment of different types of phospholipases A2 to the venoms of Caenophidian snakes: the rise of PLA2-IIE within Pseudoboini (Dipsadidae)

Snake venom harbors a wide and diverse array of enzymatic and nonenzymatic toxic components, allowing them to exert myriad effects on their prey. However, they appear to trend toward a few optimal compositional scaffolds, dominated by four major toxin classes: SVMPs, SVSPs, 3FTxs and PLA2s. Nevertheless, the latter appears to be restricted to vipers and elapids, as it has never been reported as a major venom component in rear-fanged species. Here, by investigating the original transcriptomes from 19 species distributed in eight genera from the Pseudoboini tribe (Dipsadidae: Xenodontinae) and screening among seven additional tribes of Dipsadidae and three additional families of advanced snakes, we discovered that a novel type of venom, PLA2, resembling a PLA2-IIE, has been recruited to the venom of some species of the Pseudoboini tribe, where it is a major component. Proteomic and functional analyses of these venoms further indicate that these PLA2s play a relevant role in the venoms from this tribe. Moreover, we reconstructed the phylogeny of PLA2s across different snake groups and show that different types of these toxins have been recruited in at least five independent events in caenophidian snakes. Additionally, we present the first compositional profiling of Pseudoboini venoms. Our results demonstrate how relevant phenotypic traits are convergently recruited by different means and from homologous and nonhomologous genes in phylogenetically and ecologically divergent snake groups, possibly optimizing venom composition to overcome diverse adaptative landscapes.

Venom phenotype conservation suggests integrated specialization in a lizard-eating snake

Biological specialization reduces the size of niche space while increasing efficiency in the use of available resources. Specialization often leads to phenotypic changes via natural selection aligning with niche space constraints. Commonly observed changes are in size, shape, behavior, and traits associated with feeding. One often selected trait for dietary specialization is venom, which, in snakes, often shows variation dependent on diet across and within species. The Neotropical Blunt-headed Treesnake (Imantodes cenchoa) is a highly specialized, rear-fanged, arboreal, lizard hunter that displays a long thin body, enlarged eyes, and a large Duvernoy's gland. However, toxin characterization of I. cenchoa has never been completed. Here, we use RNA-seq and mass spectrometry to assemble, annotate, and analyze the venom gland transcriptomes of four I. cenchoa from across their range. We find a lack of significant venom variation at the sequence and expression levels, suggesting venom conservation across the species. We propose this conservation provides evidence of a specialized venom repertoire, adapted to maximize efficiency of capturing and processing lizards. Importantly, this study provides the most complete venom gland transcriptomes of I. cenchoa and evidence of venom specialization in a rear-fanged snake, giving insight into selective pressures of venom across all snake species.

An integrative view of the toxic potential of Conophis lineatus (Dipsadidae: Xenodontinae), a medically relevant rear-fanged snake

Most traditional research on snake venoms has focused on front-fanged snake families (Viperidae, Elapidae, and Atractaspididae). However, venom is now generally accepted as being a much more broadly possessed trait within snakes, including species traditionally considered harmless. Unfortunately, due to historical inertia and methodological challenges, the toxin repertoires of non-front-fanged snake families (e.g., Colubridae, Dipsadidae, and Natricidae) have been heavily neglected despite the knowledge of numerous species capable of inflicting medically relevant envenomations. Integrating proteomic data for validation, we perform a de novo assembly and analysis of the Duvernoy's venom gland transcriptome of the Central American Road Guarder (Dipsadidae: Xenodontinae: Conophis lineatus), a species known for its potent bite. We identified 28 putative toxin transcripts from 13 toxin families in the Duvernoy's venom gland transcriptome, comprising 63.7% of total transcriptome expression. In addition to ubiquitous snake toxin families, we proteomically confirmed several atypical venom components. The most highly expressed toxins (55.6% of total toxin expression) were recently described snake venom matrix metalloproteases (svMMPs), with 48.0% of svMMP expression contributable to a novel svMMP isoform. We investigate the evolution of the new svMMP isoform in the context of rear-fanged snakes using phylogenetics. Finally, we examine the morphology of the venom apparatus using μCT and explore how the venom relates to autecology and the highly hemorrhagic effects seen in human envenomations. Importantly, we provide the most complete venom characterization of this medically relevant snake species to date, producing insights into the effects and evolution of its venom, and point to future research directions to better understand the venoms of 'harmless' non-front-fanged snakes.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

ToxCodAn: a new toxin annotator and guide to venom gland transcriptomics

MOTIVATION

Next-generation sequencing has become exceedingly common and has transformed our ability to explore nonmodel systems. In particular, transcriptomics has facilitated the study of venom and evolution of toxins in venomous lineages; however, many challenges remain. Primarily, annotation of toxins in the transcriptome is a laborious and time-consuming task. Current annotation software often fails to predict the correct coding sequence and overestimates the number of toxins present in the transcriptome. Here, we present ToxCodAn, a python script designed to perform precise annotation of snake venom gland transcriptomes. We test ToxCodAn with a set of previously curated transcriptomes and compare the results to other annotators. In addition, we provide a guide for venom gland transcriptomics to facilitate future research and use Bothrops alternatus as a case study for ToxCodAn and our guide.

RESULTS

Our analysis reveals that ToxCodAn provides precise annotation of toxins present in the transcriptome of venom glands of snakes. Comparison with other annotators demonstrates that ToxCodAn has better performance with regard to run time ($>20x$ faster), coding sequence prediction ($>3x$ more accurate) and the number of toxins predicted (generating $>4x$ less false positives). In this sense, ToxCodAn is a valuable resource for toxin annotation. The ToxCodAn framework can be expanded in the future to work with other venomous lineages and detect novel toxins.

MITGARD: an automated pipeline for mitochondrial genome assembly in eukaryotic species using RNA-seq data

MOTIVATION

Over the past decade, the field of next-generation sequencing (NGS) has seen dramatic advances in methods and a decrease in costs. Consequently, a large expansion of data has been generated by NGS, most of which have originated from RNA-sequencing (RNA-seq) experiments. Because mitochondrial genes are expressed in most eukaryotic cells, mitochondrial mRNA sequences are usually co-sequenced within the target transcriptome, generating data that are commonly underused or discarded. Here, we present MITGARD, an automated pipeline that reliably recovers the mitochondrial genome from RNA-seq data from various sources. The pipeline identifies mitochondrial sequence reads based on a phylogenetically related reference, assembles them into contigs, and extracts a complete mtDNA for the target species.

RESULTS

We demonstrate that MITGARD can reconstruct the mitochondrial genomes of several species throughout the tree of life. We noticed that MITGARD can recover the mitogenomes in different sequencing schemes and even in a scenario of low-sequencing depth. Moreover, we showed that the use of references from congeneric species diverging up to 30 million years ago (MYA) from the target species is sufficient to recover the entire mitogenome, whereas the use of species diverging between 30 and 60 MYA allows the recovery of most mitochondrial genes. Additionally, we provide a case study with original data in which we estimate a phylogenetic tree of snakes from the genus Bothrops, further demonstrating that MITGARD is suitable for use on biodiversity projects. MITGARD is then a valuable tool to obtain high-quality information for studies focusing on the phylogenetic and evolutionary aspects of eukaryotes and provides data for easily identifying a sample using barcoding, and to check for cross-contamination using third-party tools.

Replacement and Parallel Simplification of Nonhomologous Proteinases Maintain Venom Phenotypes in Rear-Fanged Snakes

Novel phenotypes are commonly associated with gene duplications and neofunctionalization, less documented are the cases of phenotypic maintenance through the recruitment of novel genes. Proteolysis is the primary toxic character of many snake venoms, and ADAM metalloproteinases, named snake venom metalloproteinases (SVMPs), are largely recognized as the major effectors of this phenotype. However, by investigating original transcriptomes from 58 species of advanced snakes (Caenophidia) across their phylogeny, we discovered that a different enzyme, matrix metalloproteinase (MMP), is actually the dominant venom component in three tribes (Tachymenini, Xenodontini, and Conophiini) of rear-fanged snakes (Dipsadidae). Proteomic and functional analyses of these venoms further indicate that MMPs are likely playing an "SVMP-like" function in the proteolytic phenotype. A detailed look into the venom-specific sequences revealed a new highly expressed MMP subtype, named snake venom MMP (svMMP), which originated independently on at least three occasions from an endogenous MMP-9. We further show that by losing ancillary noncatalytic domains present in its ancestors, svMMPs followed an evolutionary path toward a simplified structure during their expansion in the genomes, thus paralleling what has been proposed for the evolution of their Viperidae counterparts, the SVMPs. Moreover, we inferred an inverse relationship between the expression of svMMPs and SVMPs along the evolutionary history of Xenodontinae, pointing out that one type of enzyme may be substituting for the other, whereas the general (metallo)proteolytic phenotype is maintained. These results provide rare evidence on how relevant phenotypic traits can be optimized via natural selection on nonhomologous genes, yielding alternate biochemical components.

Size Matters: An Evaluation of the Molecular Basis of Ontogenetic Modifications in the Composition of Bothrops jararacussu Snake Venom

Ontogenetic changes in venom composition have been described in Bothrops snakes, but only a few studies have attempted to identify the targeted paralogues or the molecular mechanisms involved in modifications of gene expression during ontogeny. In this study, we decoded B. jararacussu venom gland transcripts from six specimens of varying sizes and analyzed the variability in the composition of independent venom proteomes from 19 individuals. We identified 125 distinct putative toxin transcripts, and of these, 73 were detected in venom proteomes and only 10 were involved in the ontogenetic changes. Ontogenetic variability was linearly related to snake size and did not correspond to the maturation of the reproductive stage. Changes in the transcriptome were highly predictive of changes in the venom proteome. The basic myotoxic phospholipases A₂ (PLA₂s) were the most abundant components in larger snakes, while in venoms from smaller snakes, PIII-class SVMPs were the major components. The snake venom metalloproteinases (SVMPs) identified corresponded to novel sequences and conferred higher pro-coagulant and hemorrhagic functions to the venom of small snakes. The mechanisms modulating venom variability are predominantly related to transcriptional events and may consist of an advantage of higher hematotoxicity and more efficient predatory function in the venom from small snakes.

Myriapod haemocyanin: the first three-dimensional reconstruction of Scolopendra subspinipes and preliminary structural analysis of S. viridicornis

Haemocyanins (Hcs) are copper-containing, respiratory proteins that occur in the haemolymph of many arthropod species. Here, we report the presence of Hcs in the chilopode Myriapoda, demonstrating that these proteins are more widespread among the Arthropoda than previously thought. The analysis of transcriptome of S. subspinipes subpinipes reveals the presence of two distinct subunits of Hc, where the signal peptide is present, and six of prophenoloxidase (PPO), where the signal peptide is absent, in the 75 kDa range. Size exclusion chromatography profiles indicate different quaternary organization for Hc of both species, which was corroborated by TEM analysis: S. viridicornis Hc is a 6 × 6-mer and S. subspinipes Hc is a 3 × 6-mer, which resembles the half-structure of the 6 × 6-mer but also includes the presence of phenoloxidases, since the 1 × 6-mer quaternary organization is commonly associated with hexamers of PPO. Studies with Chelicerata showed that PPO activity are exclusively associated with the Hcs. This study indicates that Scolopendra may have different proteins playing oxygen transport (Hc) and PO function, both following the hexameric oligomerization observed in Hcs.

An integrated analysis of mRNA and sRNA transcriptional profiles in Coffea arabica L. roots: insights on nitrogen starvation responses

Coffea arabica L. is an important agricultural commodity, accounting for 60% of traded coffee worldwide. Nitrogen (N) is a macronutrient that is usually limiting to plant yield; however, molecular mechanisms of plant acclimation to N limitation remain largely unknown in tropical woody crops. In this study, we investigated the transcriptome of coffee roots under N starvation, analyzing poly-A+ libraries and small RNAs. We also evaluated the concentration of selected amino acids and N-source preferences in roots. Ammonium was preferentially taken up over nitrate, and asparagine and glutamate were the most abundant amino acids observed in coffee roots. We obtained 34,654 assembled contigs by mRNA sequencing, and validated the transcriptional profile of 12 genes by RT-qPCR. Illumina small RNA sequencing yielded 8,524,332 non-redundant reads, resulting in the identification of 86 microRNA families targeting 253 genes. The transcriptional pattern of eight miRNA families was also validated. To our knowledge, this is the first catalog of differentially regulated amino acids, N sources, mRNAs, and sRNAs in Arabica coffee roots.

The complete mitochondrial genome of Bothrops jararaca (Reptilia, Serpentes, Viperidae)

The complete mitochondrial genome, containing 17,526 bp, was determined from the pitviper Bothrops jararaca. It is the first mitogenome for the most medically important genus of snake in Latin America. This mitogenome has common snake mitochondrial features such as a duplicated control region that has nearly identical sequences at two different locations of the mitogenome and a translocation of tRNA-Leu (UUR). Besides, we found a translocation of the tRNA-Pro compared to Colubridae snakes. Finally, an unusual possible duplication containing a tRNA-Phe was observed for the first time and may represent a marker of the genus.

Colubrid Venom Composition: An -Omics Perspective

Snake venoms have been subjected to increasingly sensitive analyses for well over 100 years, but most research has been restricted to front-fanged snakes, which actually represent a relatively small proportion of extant species of advanced snakes. Because rear-fanged snakes are a diverse and distinct radiation of the advanced snakes, understanding venom composition among "colubrids" is critical to understanding the evolution of venom among snakes. Here we review the state of knowledge concerning rear-fanged snake venom composition, emphasizing those toxins for which protein or transcript sequences are available. We have also added new transcriptome-based data on venoms of three species of rear-fanged snakes. Based on this compilation, it is apparent that several components, including cysteine-rich secretory proteins (CRiSPs), C-type lectins (CTLs), CTLs-like proteins and snake venom metalloproteinases (SVMPs), are broadly distributed among "colubrid" venoms, while others, notably three-finger toxins (3FTxs), appear nearly restricted to the Colubridae (sensu stricto). Some putative new toxins, such as snake venom matrix metalloproteinases, are in fact present in several colubrid venoms, while others are only transcribed, at lower levels. This work provides insights into the evolution of these toxin classes, but because only a small number of species have been explored, generalizations are still rather limited. It is likely that new venom protein families await discovery, particularly among those species with highly specialized diets.

Venom-related transcripts from Bothrops jararaca tissues provide novel molecular insights into the production and evolution of snake venom

Attempts to reconstruct the evolutionary history of snake toxins in the context of their co-option to the venom gland rarely account for nonvenom snake genes that are paralogous to toxins, and which therefore represent important connectors to ancestral genes. In order to reevaluate this process, we conducted a comparative transcriptomic survey on body tissues from a venomous snake. A nonredundant set of 33,000 unigenes (assembled transcripts of reference genes) was independently assembled from six organs of the medically important viperid snake Bothrops jararaca, providing a reference list of 82 full-length toxins from the venom gland and specific products from other tissues, such as pancreatic digestive enzymes. Unigenes were then screened for nontoxin transcripts paralogous to toxins revealing 1) low level coexpression of approximately 20% of toxin genes (e.g., bradykinin-potentiating peptide, C-type lectin, snake venom metalloproteinase, snake venom nerve growth factor) in body tissues, 2) the identity of the closest paralogs to toxin genes in eight classes of toxins, 3) the location and level of paralog expression, indicating that, in general, co-expression occurs in a higher number of tissues and at lower levels than observed for toxin genes, and 4) strong evidence of a toxin gene reverting back to selective expression in a body tissue. In addition, our differential gene expression analyses identify specific cellular processes that make the venom gland a highly specialized secretory tissue. Our results demonstrate that the evolution and production of venom in snakes is a complex process that can only be understood in the context of comparative data from other snake tissues, including the identification of genes paralogous to venom toxins.

A transcriptomic view of the proteome variability of newborn and adult Bothrops jararaca snake venoms

Background

Snake bite is a neglected public health problem in communities in rural areas of several countries. Bothrops jararaca causes many snake bites in Brazil and previous studies have demonstrated that the pharmacological activities displayed by its venom undergo a significant ontogenetic shift. Similarly, the venom proteome of B. jararaca exhibits a considerable variation upon neonate to adult transition, which is associated with changes in diet from ectothermic prey in early life to endothermic prey in adulthood. Moreover, it has been shown that the Brazilian commercial antibothropic antivenom, which is produced by immunization with adult venom, is less effective in neutralizing newborn venom effects. On the other hand, venom gland transcripts of newborn snakes are poorly known since all transcriptomic studies have been carried out using mRNA from adult specimens.

Methods/Principal Findings

Here we analyzed venom gland cDNA libraries of newborn and adult B. jararaca in order to evaluate whether the variability demonstrated for its venom proteome and pharmacological activities was correlated with differences in the structure of toxin transcripts. The analysis revealed that the variability in B. jararaca venom gland transcriptomes is quantitative, as illustrated by the very high content of metalloproteinases in the newborn venom glands. Moreover, the variability is also characterized by the structural diversity of SVMP precursors found in newborn and adult transcriptomes. In the adult transcriptome, however, the content of metalloproteinase precursors considerably diminishes and the number of transcripts of serine proteinases, C-type lectins and bradykinin-potentiating peptides increase. Moreover, the comparison of the content of ESTs encoding toxins in adult male and female venom glands showed some gender-related differences.

Conclusions/Significance

We demonstrate a substantial shift in toxin transcripts upon snake development and a marked decrease in the metalloproteinase P-III/P-I class ratio which are correlated with changes in the venom proteome complexity and pharmacological activities.

Bothrops jararaca is one of the most abundant venomous snake species in Brazil. It is primarily a nocturnal and generalist animal, however, it exhibits a notable ontogenetic shift in diet, feeding mainly on arthropods, lizards, and amphibians (ectothermic prey) through its juvenile phase and on small mammals (endothermic animals) during adult life. Due to its broad geographical distribution, this species is responsible for the majority of the accidents by Bothrops genus in Brazil. Studies on envenomation cases with newborn and adult B. jararaca snakes have shown distinct patterns, mainly related to blood coagulation disorders, which seems to be prominent in accidents with newborn specimens. Moreover, it has been demonstrated that the Brazilian commercial antibothropic antivenom, which is produced by immunization with adult venom, is less effective in neutralizing newborn venom effects. In this study we analyzed the venom gland transcriptome of newborn snake specimens and compared the content of toxin transcripts with that of adult specimens. We demonstrate that upon B. jararaca development, its repertoire of mRNAs encoding toxins changes both qualitatively and quantitatively and these alterations are associated with the venom proteome profiles and pharmacological activities displayed by newborn and adult specimens.

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 (∼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.

A new Factor Xa inhibitor from Amblyomma cajennense with a unique domain composition

Bioactive compounds of great interest are found in the saliva of hematophagous organisms. While exploring a cDNA library derived from the salivary glands of the tick Amblyomma cajennense, a transcript that codes for a protein with unique structure (containing an N-terminal Kunitz-type domain and a C-terminus with no homology to any annotated sequences) was found. The recombinant mature form of this protein ( approximately 13.5kDa) was produced in Escherichia coli BL21 (DE3), and it was able to inhibit Factor Xa (FXa) and extend global blood clotting times in vitro and ex vivo. Static and dynamic predictions of its tertiary structure indicate regions that may be related to its FXa inhibitor function.

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.

Transcriptome analysis of expressed sequence tags from the venom glands of the fish Thalassophryne nattereri

Thalassophryne nattereri (niquim) is a venomous fish found on the northern and northeastern coasts of Brazil. Every year, hundreds of humans are affected by the poison, which causes excruciating local pain, edema, and necrosis, and can lead to permanent disabilities. In experimental models, T. nattereri venom induces edema and nociception, which are correlated to human symptoms and dependent on venom kininogenase activity; myotoxicity; impairment of blood flow; platelet lysis and cytotoxicity on endothelial cells. These effects were observed with minute amounts of venom. To characterize the primary structure of T. nattereri venom toxins, a list of transcripts within the venom gland was made using the expressed sequence tag (EST) strategy. Here we report the analysis of 775 ESTs that were obtained from a directional cDNA library of T. nattereri venom gland. Of these ESTs, 527 (68%) were related to sequences previously described. These were categorized into 10 groups according to their biological functions. Sequences involved in gene and protein expression accounted for 14.3% of the ESTs, reflecting the important role of protein synthesis in this gland. Other groups included proteins engaged in the assembly of disulfide bonds (0.5%), chaperones involved in the folding of nascent proteins (1.4%), and sequences related to clusterin (1.5%), as well as transcripts related to calcium binding proteins (1.0%). We detected a large cluster (1.3%) related to cocaine- and amphetamine-regulated transcript (CART), a peptide involved in the regulation of food intake. Surprisingly, several retrotransposon-like sequences (1.0%) were found in the library. It may be that their presence accounts for some of the variation in venom toxins. The toxin category (18.8%) included natterins (18%), which are a new group of kininogenases recently described by our group, and a group of C-type lectins (0.8%). In addition, a considerable number of sequences (32%) was not related to sequences in the databases, which indicates that a great number of new toxins and proteins are still to be discovered from this fish venom gland.

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.

Natterins, a new class of proteins with kininogenase activity characterized from Thalassophryne nattereri fish venom

A novel family of proteins with kininogenase activity and unique primary structure was characterized using combined pharmacological, proteomic and transcriptomic approaches of Thalassophryne nattereri fish venom. The major venom components were isolated and submitted to bioassays corresponding to its main effects: nociception and edema. These activities were mostly located in one fraction (MS3), which was further fractionated. The isolated protein, named natterin, was able to induce edema, nociception and cleave human kininogen and kininogen-derived synthetic peptides, releasing kallidin (Lys-bradykinin). The enzymatic digestion was inhibited by kallikrein inhibitors as Trasylol and TKI. Natterin N-terminal peptide showed no similarity with already known proteins present in databanks. Primary structure of natterin was obtained by a transcriptomic approach using a representative cDNA library constructed from T. nattereri venom glands. Several expressed sequence tags (ESTs) were obtained and processed by bioinformatics revealing a major group (18%) of related sequences unknown to gene or protein sequence databases. This group included sequences showing the N-terminus of isolated natterin and was named Natterin family. Analysis of this family allowed us to identify five related sequences, which we called natterin 1-4 and P. Natterin 1 and 2 sequences include the N-terminus of the isolated natterin. Furthermore, internal peptides of natterin 1-3 were found in major spots of whole venom submitted to mass spectrometry/2DGE. Similarly to the ESTs, the complete sequences of natterins did not show any significant similarity with already described tissue kallikreins, kininogenases or any proteinase, all being entirely new. These data present a new task for the knowledge of the action of kininogenases and may help in understanding the mechanisms of T. nattereri fish envenoming, which is an important medical problem in North and Northeast of Brazil.

Comparative genomics of two Leptospira interrogans serovars reveals novel insights into physiology and pathogenesis

Leptospira species colonize a significant proportion of rodent populations worldwide and produce life-threatening infections in accidental hosts, including humans. Complete genome sequencing of Leptospira interrogans serovar Copenhageni and comparative analysis with the available Leptospira interrogans serovar Lai genome reveal that despite overall genetic similarity there are significant structural differences, including a large chromosomal inversion and extensive variation in the number and distribution of insertion sequence elements. Genome sequence analysis elucidates many of the novel aspects of leptospiral physiology relating to energy metabolism, oxygen tolerance, two-component signal transduction systems, and mechanisms of pathogenesis. A broad array of transcriptional regulation proteins and two new families of afimbrial adhesins which contribute to host tissue colonization in the early steps of infection were identified. Differences in genes involved in the biosynthesis of lipopolysaccharide O side chains between the Copenhageni and Lai serovars were identified, offering an important starting point for the elucidation of the organism's complex polysaccharide surface antigens. Differences in adhesins and in lipopolysaccharide might be associated with the adaptation of serovars Copenhageni and Lai to different animal hosts. Hundreds of genes encoding surface-exposed lipoproteins and transmembrane outer membrane proteins were identified as candidates for development of vaccines for the prevention of leptospirosis.

A prothrombin activator from Bothrops erythromelas (jararaca-da-seca) snake venom: characterization and molecular cloning

A novel prothrombin activator enzyme, which we have named 'berythractivase', was isolated from Bothrops erythromelas (jararaca-da-seca) snake venom. Berythractivase was purified by a single cation-exchange-chromatography step on a Resource S (Amersham Biosciences) column. The overall purification (31-fold) indicates that berythractivase comprises about 5% of the crude venom. It is a single-chain protein with a molecular mass of 78 kDa. SDS/PAGE of prothrombin after activation by berythractivase showed fragment patterns similar to those generated by group A prothrombin activators, which convert prothrombin into meizothrombin, independent of the prothrombinase complex. Chelating agents, such as EDTA and o -phenanthroline, rapidly inhibited the enzymic activity of berythractivase, like a typical metalloproteinase. Human fibrinogen A alpha-chain was slowly digested only after longer incubation with berythractivase, and no effect on the beta- or gamma-chains was observed. Berythractivase was also capable of triggering endothelial proinflammatory and procoagulant cell responses. von Willebrand factor was released, and the surface expression of both intracellular adhesion molecule-1 and E-selectin was up-regulated by berythractivase in cultured human umbilical-vein endothelial cells. The complete berythractivase cDNA was cloned from a B. erythromelas venom-gland cDNA library. The cDNA sequence possesses 2330 bp and encodes a preproprotein with significant sequence similarity to many other mature metalloproteinases reported from snake venoms. Berythractivase contains metalloproteinase, desintegrin-like and cysteine-rich domains. However, berythractivase did not elicit any haemorrhagic response. These results show that, although the primary structure of berythractivase is related to that of snake-venom haemorrhagic metalloproteinases and functionally similar to group A prothrombin activators, it is a prothrombin activator devoid of haemorrhagic activity. This is a feature not observed for most of the snake venom metalloproteinases, including the group A prothrombin activators.