A transcriptomic analysis of gene expression in the venom gland of the snake Bothrops alternatus (urutu)
BMC Genomics volume 11, Article number: 605 (2010)
The genus Bothrops is widespread throughout Central and South America and is the principal cause of snakebite in these regions. Transcriptomic and proteomic studies have examined the venom composition of several species in this genus, but many others remain to be studied. In this work, we used a transcriptomic approach to examine the venom gland genes of Bothrops alternatus, a clinically important species found in southeastern and southern Brazil, Uruguay, northern Argentina and eastern Paraguay.
A cDNA library of 5,350 expressed sequence tags (ESTs) was produced and assembled into 838 contigs and 4512 singletons. BLAST searches of relevant databases showed 30% hits and 70% no-hits, with toxin-related transcripts accounting for 23% and 78% of the total transcripts and hits, respectively. Gene ontology analysis identified non-toxin genes related to general metabolism, transcription and translation, processing and sorting, (polypeptide) degradation, structural functions and cell regulation. The major groups of toxin transcripts identified were metalloproteinases (81%), bradykinin-potentiating peptides/C-type natriuretic peptides (8.8%), phospholipases A2 (5.6%), serine proteinases (1.9%) and C-type lectins (1.5%). Metalloproteinases were almost exclusively type PIII proteins, with few type PII and no type PI proteins. Phospholipases A2 were essentially acidic; no basic PLA2 were detected. Minor toxin transcripts were related to L-amino acid oxidase, cysteine-rich secretory proteins, dipeptidylpeptidase IV, hyaluronidase, three-finger toxins and ohanin. Two non-toxic proteins, thioredoxin and double-specificity phosphatase Dusp6, showed high sequence identity to similar proteins from other snakes. In addition to the above features, single-nucleotide polymorphisms, microsatellites, transposable elements and inverted repeats that could contribute to toxin diversity were observed.
Bothrops alternatus venom gland contains the major toxin classes described for other Bothrops venoms based on trancriptomic and proteomic studies. The predominance of type PIII metalloproteinases agrees with the well-known hemorrhagic activity of this venom, whereas the lower content of serine proteases and C-type lectins could contribute to less marked coagulopathy following envenoming by this species. The lack of basic PLA2 agrees with the lower myotoxicity of this venom compared to other Bothrops species with these toxins. Together, these results contribute to our understanding of the physiopathology of envenoming by this species.
Studies of the global composition of snake venoms have advanced dramatically in recent years with the introduction of transcriptomic and proteomic approaches for analyzing venom gland gene expression and venom composition [1–3]. The introduction of these two major fields in toxinology has created new perspectives for identifying novel molecules for drug discovery [4, 5] and for improving antivenom development and the clinical treatment of snakebite [6–8]. A combination of these two approaches can be particularly useful in providing a more comprehensive understanding of venom composition [3, 9–12].
Detailed transcriptomic studies of snake venom glands began with the report of Junqueira-de-Azevedo and Ho  who analyzed expressed sequence tags (ESTs) from the venom gland of the Brazilian pitviper Bothrops insularis, a species endemic to the island of Queimada Grande off the coast of the State of São Paulo in southeastern Brazil. Since then, venom gland transcriptomes have been reported for various other South American venomous snakes, including Bothrops species [11, 14–16], Crotalus durissus collilineatus (South American rattlesnake) , Lachesis muta (bushmaster) , Micrurus corallinus (coral snake)  and the colubrid Philodryas olfersii.
The genus Bothrops is responsible for most cases of snakebite throughout Latin America  and in Brazil accounts for 80-90% of bites by venomous snakes. Envenoming by Bothrops species results in extensive local effects, including pain, edema, inflammation, hemorrhage and necrosis [22–24], as well as systemic actions that include coagulopathy, internal hemorrhage, circulatory shock and acute renal failure [21, 25–28]. These actions are mediated by a variety of venom components, with the most extensively studied being (hemorrhagic) metalloproteinases  and myotoxic phospholipases A2 (PLA2) [23, 30, 31].
Bothrops alternatus occurs in southeastern and southern Brazil, Uruguay, central and northern Argentina and south-eastern and southern Paraguay . This species is an important cause of snakebite, although the prevalence of bites varies considerably throughout its geographic distribution depending on the regional human and snake population densities and the presence of more abundant species such as Bothrops jararaca, the major cause of snakebite in southeastern Brazil. Thus, for example, notifications of bites by Brazilian Bothrops spp. received by the Instituto Butantan from 1902-1945 (n = 4,881) suggest that B. alternatus accounts for ~6% of bites by this genus . In contrast, in a more recent series of 3,139 cases attended at the Hospital Vital Brazil, Instituto Butantan in São Paulo, from 1981-1990, only two (0.14%) of the 1,412 cases in which the Bothrops species was identified were caused by B. alternatus compared to 1,376 (97.5%) by B. jararaca. In regions where B. alternatus is more abundant or predominates (southern Brazil, Uruguay and parts of Argentina), this species contributes to much a greater proportion of snakebites, e.g., 18% in Argentina, but may be much higher locally (≥ 60% of snakebites in the Argentinian province of Corrientes) .
Bothrops alternatus produces clinical manifestations characteristic of this genus [21, 28]. In a series of 32 cases involving this species , most of the patients (74.2%) were 15-49 years old and were bitten in the lower limbs (84.3%). All patients developed local pain, edema and most (31; 96.9%) had incoagulable blood, 13 (40.6%) developed hemorrhage, 7 (21.9%) had blisters and 3 (9.4%) showed necrosis; there were no fatalities. In addition, renal failure has been observed after bites by this species in Uruguay . The administration of antivenom (produced by the Instituto Nacional de Producción de Biológicos - INPB in Argentina and the Instituto Butantan, Fundação Ezequiel Dias or Instituto Vital Brazil in Brazil) is the standard treatment for envenoming by B. alternatus. These antivenoms may be bivalent (raised against B. alternatus and B. neuwiedii in Argentina) or poplyvalent (raised against the clinically most relevant Bothrops species in Argentina and Brazil - B. alternatus, B. jararaca, B. jararacussu, B. neuwiedi and B. moojeni) [38, 39]. In the case series studied by Bauab et al. , a median of 40 ml of antivenom (four 10 ml ampoules) was administered.
Although B. alternatus venom contains a variety of enzymatic and biological activities, relatively few of the venom proteins involved have been purified and characterized, with most investigations having dealt with metalloproteinases and disintegrins [40–50], PLA2[51–55], coagulant enzymes [56–58], L-amino acid oxidase (LAO)  and phosphodiesterase . Despite these investigations, the venom of B. alternatus is still less understood than those of other Bothrops species such as B. jararaca and B. jararacussu. In order to obtain a more comprehensive understanding of the toxinology of this species, we have used a transcriptomic approach to examine venom gland gene expression in B. alternatus and compared the results with those for other members of this genus.
Venom glands and RNA isolation
Venom glands from three adult B. alternatus snakes were obtained three days after venom extraction. The choice of this interval was based essentially on other studies [16, 17, 61, 62] that used this period, although intervals of two  and four [11, 13] days post-milking have also been used for Bothrops species. The choice of a 2-4 day post-milking interval prior to gland removal is partly based on studies of venom mRNA and protein synthesis in snake venom glands that show maximal production 3-8 days post-milking [reviewed in ], and partly on histological and functional studies in B. jararacussu and B. jararaca, respectively, showing that changes in gland epithelial morphology and venom production peak around four days post-milking.
Each pair of venom glands was homogenized in liquid nitrogen and total RNA was extracted with TRIzol reagent (Invitrogen, UK), according to the manufacturer's instructions. Three cDNA libraries (one from each snake) were independently constructed with CloneMiner cDNA library construction kits (Invitrogen, UK), according to the manufacturer's instructions. The cDNA libraries were processed and analyzed using an in-house bioinformatics pipeline, with all annotations being done manually (contig by contig and singlet by singlet) for each library. The production of three independent libraries meant that it was possible to undertake detailed analyses of SNPs, microsatellites and inversions not generally done for Bothrops species since in most studies the venom glands of different snakes are usually combined into a single pool for mRNA extraction, with a subsequent loss of information.
This work was approved by the institutional Committee for Ethics in Animal Experimentation (CEEA/UNICAMP, protocol no. 864-1) and was done according to the ethical guidelines of the Brazilian Society of Laboratory Animal Science (SBCAL, formerly the Brazilian College for Animal Experimentation - COBEA).
The cDNA libraries were sequenced using BigDye terminator 3.0 kits and an automated DNA capillary sequencer (ABI PRISM 3700 DNA Analyzer, Applied Biosystems, Foster City, CA, USA). All of the cDNA sequences were 5' sequenced using the primer M13F (5'-TGTAAAACGACGGCCAGT-3').
Clusterization, assembly and identification of Bothrops alternatus expressed-sequence tags
The Phred program  was used to obtain sequences and quality files from chromatograms obtained from expressed-sequence tag (EST) sequencing. The EST cleaning pipeline described by Baudet and Dias  was then used to pre-process the ESTs and prepare the sequences for assembly. This pipeline removes sequences with plasmid similarity, polyA/polyT regions, low base quality and slippage signals. Sequences <100 bp long after cleaning were discarded. CAP3 software  was used to cluster and assemble the clean sequences into contigs and singlets (unisequences). For assembly, an overlap of 100 bp and an identity of at least 95% were used as criteria to detect pairwise similarities.
Annotation of Bothrops alternatus ESTs
After clustering and assembly, a BLAST search was done to identify similarities between the ESTs and sequences deposited in public databases. All of the sequences were aligned against the GenBank non-redundant (nr) protein database using BLASTX and BLASTN  with an E-value cut-off of 1e-5. The B. alternatus ESTs were also screened against two locally generated sequence databases, SerpP and SerpN, that included all snake protein and nucleotide sequences from GenBank, respectively. In addition, the ESTs were compared with the complete genome of the lizard Anolis carolinensis (http://genome.ucsc.edu/cgi-bin/hgGateway?db=anoCar1). Gene Ontology annotation was done with Blast2GO  using GO-slim terms. The uncharacterized ESTs were examined for the presence of a signal peptide by using SignalP 3.0 software (http://www.cbs.dtu.dk/services/SignalP/).
Sequence alignments for selected proteins were done with the program ClustalW .
Single nucleotide polymorphisms
The software QualitySNP  was used to identify single-nucleotide polymorphisms (SNPs). Non-synonymous and synonymous SNPs (nsSNPs and sSNPs, respectively) were identified by detecting open-reading frames (ORFs) of contigs with SNPs using the FASTA algorithm run against the 2008 version of UniProt . The possibility of SNPs arising from artifacts during DNA sequencing was minimized by the fact that the cDNA libraries were prepared independently from three snakes and that we used consensus sequences from contigs with at least three reads from separate sequencing plates for which the cDNA was prepared and the reactions run on different days. These procedures considerably reduced the possiblity of artifacts derived from DNA sequencing and strengthened our conclusions regarding the presence of SNPs.
Identification of transposable elements and long inverted repeats
Alignment of the unisequences to repetitive elements in RepBase release 14.08  was done with BLASTN that was automated using in-house Perl scripts (available upon request). The E-value cut-off was set at 1 × 10-10 and only alignments of at least 50 bp were considered for unisequences. In addition, the alignments with database sequences had to show >80% identity over at least 10% of their lengths. To identify possible functional integrations of transposable elements (TEs) into host genes, we searched for chimeric transcripts between TEs and protein-coding genes. Unisequences that aligned with a TE in the database over 10% to 90% of their sequence were used to match the sequence back to the corresponding host gene by masking the element in the sequence and querying the NCBI nr nucleotide database. For tBLASTX remote searches, a minimum of 50% of the sequence was required to be involved in the best hit, with at least 70% identity.
Additional repeated elements that were not present in RepBase were observed in some unisequences. These novel elements were identified using the palindrome program in the EMBOSS package . Unisequences that contained inverted repeats were identified by aligning each transcript against itself using BLASTN (E-value cut-off of 1e-5) and then inspected visually to check for alignments in opposite strands. The upper limit for mismatches between the two repeated segments was 10%.
Comparison of the B. alternatus EST library with other Bothrops species
The pattern of gene expression in B. alternatus venom gland was compared with EST data for other Bothrops species (B. atrox, B. insularis[11, 13], B. jararaca and B. jararacussu) available in public databases, and with the major toxin classes detected by proteomic analyses of Bothrops venoms. In this comparison and in the description of the different toxin groups discussed below, the relative percentages of the venom gland EST categories were calculated as a percentage of the overall (total) number of all ESTs identified, whereas the relative percentages of the ESTs related to specific toxin groups was calculated based on the total number of toxin ESTs identified.
Results and Discussion
Venom gland EST database
cDNA libraries were generated from three pairs of B. alternatus venom glands and analyzed using the pipeline shown in Figure 1. A total of 5,350 valid unisequences was obtained by single-pass sequencing and, after trimming, the lengths of the sequences ranged from 100 bp to 848 bp (mean: 562 bp) (Figure 2). The ESTs were assembled into 838 contigs and 4512 singletons. Thirty percent of the sequences (1,605 ESTs) were hits, of which 1,245 (23% of total ESTs and ~78% of hits) were related to toxins; the remaining 70% (3,745 ESTs) were no-hits. Figure 3 shows the proportion of hits and no-hits for selected databases. The proportion of hits was greater against the GenBank and SerpN nucleotide databases compared to the corresponding protein databases. There was a particularly low number of hits in SerpP, which likely reflects the limited number of non-venom snake protein sequences available in public databases.
SignalP analysis identified 147 ESTs with signal peptides, 62 of which were related to snake toxins (37 for PLA2, 7 for metalloproteinases, 6 for proteins related to platelet glycoprotein 1β binding protein, 6 for proteins related to ACF 1/2 β-chain, 4 for serine proteinases, and 1 each for disintegrins and bradykinin-potentiating peptides or BPPs), 50 were for other proteins and 35 were for unknown proteins.
All of the ESTs were deposited in the dbEST division of GenBank (accession numbers GW575430 to GW583300).
Gene ontology (GO)
A global GO analysis  of unisequences in relation to molecular function revealed that the largest number of transcripts was related to protein binding, followed by peptidase and transferase activities, calcium-binding, enzyme regulation, nucleic acid binding, structural protein and transporter activities, lipid binding, and transcriptional regulation (Figure 4). Analysis of the biological processes revealed two major groups arbitrarily defined as those with >100 ESTs and those with ≤100 ESTs. The former contained genes related to stress, responses to external stimuli, development, symbiosis, cell differentiation, DNA metabolism, cell death and signal transduction, whereas the latter consisted of genes related to catabolism, precursor biosynthesis, translation, protein modification, transcription, lipid metabolism, protein and ion transport, cytoskeletal organization, cell proliferation and cell cycle (Figure 5). The biological processes identified here generally agreed with the broad categories (general metabolism, transcription and translation, processing and sorting, (polypeptide) degradation, structural functions and cell regulation) reported for other Bothrops species [11, 13–16] and L. muta.
Selected venom gland proteins
Venom gland receptors
We detected one transcript coding for a protein similar to α1D-adrenoceptor from B. jararaca. Since α- and β-adrenergic mechanisms have an important role in venom production in B. jararaca[76–78], principally through receptors structurally related to α1D- and α1B-adrenoceptors , a similar regulatory mechanism mediated by these receptors probably also occurs in B. alternatus. Another receptor identified was vascular endothelin type A (ETA) receptor (one singlet). This receptor mediates the vasoconstrictor activity of the endothelin family of peptides in B. jararaca aorta  and could be important in the regulation of vascular tone in blood vessels irrigating the venom gland.
A contig was identified for a calglandulin-like EF-hand protein. Venom gland calglandulin, initially identified in B. insularis based on an EST analysis of this venom gland  and subsequently cloned and expressed from this species  (also detected in elapid venom glands [9, 81]), is a Ca2+-binding protein structurally similar to calmodulin and troponin C. This protein, which is venom gland-specific (not found in other snake organs or venom) , may have a role in the transport and secretion of venom toxins, in addition to possibly acting as an intracellular Ca2+-chelator to regulate toxin activity.
Dual specificity phosphatases (Dusp) have important functions in embryogenesis, cell growth and immune responses by acting as negative feedback regulators of mitogen-activated protein kinases (MAPK) . DUSP6 has an important role in snake embryogenesis, with high gene expression in the anterior part of the presomitic mesoderm [83, 84]. This protein has not previously been detected in venom gland transcriptomes or in analyses of venom composition. We detected a partial sequence for DUSP6 protein that shared similarity with the corresponding protein from Pantherophis guttatus (corn snake; accession number ABW82165) (Figure 6). The role of DUSP6 in the venom gland is unclear but may be related to development of the gland secretory epithelium and venom production.
Thioredoxin (Trx) is a protein with a variety of activities, including roles in DNA synthesis, protein disulfide bond reduction and the degradation of H2O2 that may be related to protection against oxidative stress and the induction of apoptosis . Thioredoxin participates in redox reactions through the reversible oxidation of its active center dithiol to a disulfide, thereby catalyzing dithiol-disulfide exchange reactions; the reduction of the active site disulfide of oxidized Trx to regenerate the dithiol of reduced Trx is catalyzed by thioredoxin reductase, a selenium-containing flavoprotein . Trx is also involved in the reversible S-nitrosylation of cysteine residues in target proteins, an important step in signaling by intracellular nitric oxide (NO). We found a protein sequence in the B. alternatus transcriptome that was related to Trx previously identified in venom gland cDNA from O. hannah (accession number AAK09384) (Figure 7). The function of Trx in venom glands is unknown, but it may be involved in protecting epithelial secretory cells of the gland from oxidative stress and death (by necrosis or apoptosis), particularly since venom components such as LAO can cause cell death through the formation of H2O2[87, 88].
Venom components identified by ESTs
Figure 8 summarizes the venom components detected based on EST analysis. The major classes detected were metalloproteinases/disintegrins, BPP/C-type natriuretic peptide (CNP) precursors, PLA2, serine proteinases and C-type lectins. Genes expressed in lower abundance included cysteine-rich secretory proteins (CRISPs), taicatoxin-like protein, prothrombin activator, a catrocollastatin precursor and dipeptidylpeptidase IV (DPP IV).
Major toxin classes
Bothrops venoms are rich in a variety of snake venom metalloproteinases (SVMPs) that contribute to local and systemic bleeding after envenoming [21, 28, 29]. SVMPs are currently classified in three major classes (PI, PII and PIII, with the latter two containing five and four subclasses, respectively) , with PIII metalloproteinases being particularly abundant and extensively studied in Bothrops venoms. Metalloproteinases were the most abundant venom components in B. alternatus venom gland (Figure 8) and were almost exclusively class PIII proteins. SVMPs accounted for 1010 ESTs that were grouped into 22 contigs, the most abundant of which were related to jararhagin from B. jararaca, Russell's viper (Vipera russelli) venom factor × heavy chain (RVV-X) and hemorrhagic factor 3 (HF3), also from B. jararaca (all class PIII proteins) (Figure 9). Hits were also obtained for Bothrops metalloproteinases such as berythactivase from B. erythromelas, bothropasin from B. jararaca and metalloproteinases II (Bojumet II) and III (Bojumet III) from B. jararacussu. There was a low number of transcripts for class PII proteins, including one transcript coding for contortrostatin and one contig related to insularinase. No class PI SVMPs were detected. The overwhelming abundance of transcripts for class PIII SVMPs (with few PII and no PI SVMPs) observed here agrees well with the results of a proteomic analysis of this venom in which only class PIII SVMPs were detected . Our finding of transcripts related to jararhagin and berythactivase also agrees with these authors who detected proteins related to these two SVMPs in B. alternatus venom.
In agreement with the abundance of SVMP transcripts, two PIII metalloproteinases have been isolated and characterized from this venom. Souza et al. reported the characterization of a 55 kDa PIII metalloproteinase (alternagin) from Brazilian B. alternatus venom that can undergo autolysis to release a 28 kDa ECD-disintegrin-like cysteine-rich domain (alternagin-C). Alternagin inhibits the binding of K562 cells to collagen by selectively blocking α2β1 integrin, in a manner similar to jararhagin from B. jararaca venom; this inhibitory action is apparently mediated by alternagin-C. Subsequently, Gay et al.[44, 45] also isolated alternagin (which they termed balteragin) from Argentine B. alternatus. This protein causes edema, hemorrhage and necrosis when administered intramuscularly in mice, and systemic hemorrhage, primarily in the lungs, kidneys and liver, when administered intravenously. Cominetti et al. described the characterization of a 130 kDa dimeric PIII metalloproteinase/RDG-disintegrin (BaG; monomeric form = 55 kDa) from Brazilian B. alternatus venom. This enzyme, which accounts for at least 0.2% of venom protein, inhibits ADP-induced platelet aggregation via a mechanism independent of its enzymatic activity. BaG also blocks the adhesion of K562 cells to fibronectin (and detaches cells previously adhered to this protein), a phenomenon mediated by binding to α5β1 integrin, but has no effect on the adhesion of these cells to collagen type I (mediated by α2β1 integrin).
The processing of PII and PIII SVMPs gives rise to RDG-disintegrins and ECD-disintegrin-like cysteine-rich (DC) domains, respectively, that exert a variety of biological activities through interaction with cell surface integrins . In the B. alternatus ESTs, there were four ESTs coding for disintegrins, and at least three disintegrins were identified in a proteomic analysis of this venom . As indicated above, autolysis of alternagin gives rise to alternagin-C, an ECD-disintegrin-like cysteine-rich domain . Alternagin-C exerts a variety of activities, including the induction of neutrophil migration via integrin signaling , stimulation of human umbilical vein endothelial cell migration [43, 46] and modulation of angiogenesis in vitro and in vivo, but has a limited beneficial effect on muscle regeneration . More recently, another disintegrin (DisBa-01) identified in a cDNA library from B. alternatus venom gland, has been shown to inhibit angiogenesis and melanoma mestastasis by interacting with αvβ3 integrin  and prevents platelet adhesion to fibrinogen, thereby prolonging bleeding time in mice ; the latter action is mediated by interaction with platelet αIIbβ3 integrin.
Bothrops venoms are rich in acidic and basic PLA2 with Asp or Lys at position 49 in their active site. Many of these PLA2 cause pain, edema, inflammation and myotoxicity [23, 30], as well as neuromuscular blockade in vitro. The content of myotoxic PLA2 varies considerably among Bothrops species, e.g., B. asper, B. colombiensis, B. fonsecai, B. jararacussu, B. moojeni, B. neuwiedi and B. pradoi have an elevated content of these enzymes whereas others, such as B. alternatus, B. atrox, B. cotiara, B. erythromelas and B. jararaca, have few or no myotoxins [6, 92–96]. In addition, some species, e.g., B. jararacussu[14, 97] contain a variety of PLA2 in addition to myotoxins. PLA2 accounted for 5.6% of B. alternatus venom components (1.28% of the total EST library) (Figure 8). This PLA2 content was similar to that of B. insularis (6.7%)  and B. jararaca (6.7%) , but considerably less than for B. atrox (14.6%; ) and B. jararacussu (58%; ). The B. alternatus PLA2 formed two major groups: those related to BOTIN PLA2 (PA2-BOTIN or BinTX-I; SwissProt/trEMBL accession number Q8QG87), one of two acidic Asp49 PLA2 isoforms (BinTX-I = pI 5.05; BinTX-II = pI 4.49) characterized from the venom of B. insularis, with 41 ESTs, and those related to an acidic (pI 4.5) hypotensive Asp49 PLA2 (Q8AXY1) from B. jararacussu venom , with 28 ESTs.
BinTX-I causes mouse hind-paw edema, as well as myonecrosis and partial neuromuscular blockade in chick biventer cervicis preparations . Hence, it is probable that the corresponding PLA2 in B. alternatus venom exert similar effects. In contrast, the acidic hypotensive PLA2 from B. jararacussu is not myotoxic, cytotoxic or lethal but causes edema and hypotension, and inhibits platelet aggregation . Nisenbom et al.[51, 52] reported the characterization of an acidic (pI 5.04-5.08) PLA2 from B. alternatus that accounted for most of the lethality of this venom in mice and produced cardiovascular alterations such as dyspnea, tachycardia, arrhythmia and circulatory shock, as well as tissue damage (hemorrhage, necrosis, etc.). This enzyme therefore shares properties with B. insularis and B. jararacussu PLA2.
Basic Lys49 PLA2s characteristic of Bothrops venoms were not detected in this venom gland library, in agreement with a proteomic analysis of this venom  and with EST analyses of B. insularis and the bushmaster L. muta. However, two ESTs coding for the basic PLA2 of crotoxin (the principal neurotoxin of C. d. terrificus; also identified in B. insularis ESTs; ) and one transcript for the PLA2 β-neurotoxin ammodytoxin (from Vipera ammodytes) were detected. The general absence of basic Lys49 PLA2 in the B. alternatus transcriptome agrees with the very low content of myotoxic PLA2 already reported for this venom [93, 95]. Indeed, to date, only one basic (pI 8.63) myotoxic Lys49 PLA2 has been isolated and characterized from this venom [54, 55]. This PLA2, which accounts for <1% of the venom content, causes mouse footpad edema, skeletal muscle myonecrosis, lysis of C2C12 skeletal muscle myoblasts and presynaptic neuromuscular blockade in chick biventer cervicis and mouse phrenic nerve-diaphragm preparations. However, given the very low level of this PLA2 in the venom, its overall contribution to the activities of the venom is likely to be minimal.
In addition to PLA2, inhibitors of this enzyme were also detected (17 ESTs with hits only in the 3'UTR; similar to PLA2 inhibitors of Trimeresurus flavoviridis). Venom gland PLA2 inhibitors show considerable sequence identity with γ-PLI (phospholipase inhibitors) identified in snake sera, including that of B. alternatus.
Coagulopathy is a major complication of systemic envenoming by Bothrops species that, in conjunction with the action of venom metalloproteinases, can contribute to systemic hemorrhage [21, 28]. Coagulation disturbances caused by Bothrops venoms are mediated by a variety of enzymes that act at various points of the coagulation cascade, including fibrinogen degradation and fibrin consumption. The clinical outcome of this combined action is incoagulable blood. Some of the serine proteinases may also have kinin-releasing activity that can contribute to local pain, vasodilation and systemic hypotension [101, 102]. In the B. alternatus transcriptome, serine proteinases were the fourth most abundant group of toxins (Figure 8) (1.9% and 0.45% of toxin and total transcripts, respectively), with ~42% (10 ESTs) being related to BthaTL, a serine proteinase identified in this venom by Vitorino-Cardoso et al. based on gene cloning and amino acid sequence analysis. This protein, which shares sequence identity (> 60%) with a variety of snake venom thrombin-like enzymes (TLE), is probably the same as the TLE balterobin previously purified and characterized by Smolka et al.. Other proteins in this group were related to serine proteinase isoform 5 (7 ESTs), venom serum proteinase-like protein 2 precursor (5 ESTs) and venom serine proteinase HS112 precursor (1 EST) and TLE (1 EST). In addition, there was one transcript with a hit in a UTR for a protein similar to KN-BJ2, a kinin-releasing, fibrinogen clotting enzyme from B. jararaca venom [15, 102, 103]. The range of serine proteinases identified here in ESTs agrees with the presence of coagulant (BthaTL and several TLEs homologous to ancrod, batroxobin, proteinase HS114 and pallase) and kinin-releasing enzymes identified in the proteome of this venom . The abundance of serine proteinases in the B. alternatus transcriptome was similar to the 2% reported for B. jararacussu but considerably less than in other species of this genus, e.g., ranging from 8.2% in B. atrox to 28.6% in B. jararaca. Differences in the relative proportion of serine proteinases present in the venoms could contribute to variations in the severity of coagulopathy among Bothrops species.
Proteins with C-type lectin-like domains are widespread in the animal kingdom, including venoms. Venom C-type lectins form two major groups: true C-type lectins with carbohydrate binding domains and C-type lectin-like proteins with non-carbohydrate binding domains, the latter consisting mainly of factor IX/X binding proteins and proteins that interact with platelet receptors. C-type lectins are involved in venom-induced coagulation disturbances and have been identified in various Bothrops ESTs [11, 13–16]. These proteins formed the fifth most abundant class of toxin transcripts in B. alternatus (18 ESTs with protein, 1.46%, or 0.33% of the total cDNA library) (Figure 8) and were related to the α subunit of platelet glycoprotein Ib-binding protein (8 ESTs), the ACF 1/2 β chain (7 ESTs), coagulation factor IX-binding protein (2 ESTs) and bothrojaracin (1 EST). There was also a transcript for convulxin (one singlet in the 3'UTR), a C-type lectin initially identified in the venom of C. d. terrificus (also identified in B. insularis ESTs; ). The abundance of C-type lectins was considerably less than in B. atrox (4.8%; ), B. insularis (14.5%; ), B. jararaca (8.3%; ) and B. jararacussu (5%; ); if these relative proportions are reflected in the venoms, then C-type lectins may contribute less to coagulopathy in bites by B. alternatus than by other Bothrops spp.
Few proteins of this class have been characterized from B. alternatus venom. Bothrojaracin, a 27 kDa heterodimer (13 kDa and 15 kDa) initially purified from B. jararaca venom, is homologous to several C-type lectin family proteins and inhibits thrombin activity, e.g., platelet aggregation and fibrinogen clotting, by interacting with exosites I and II of this enzyme to form a non-covalent complex [104, 105]. Castro et al. characterized bothroalternin, a bothrojaracin-like protein, from B. alternatus venom; this protein 27 kDa homodimer cross-reacts with antibodies to bothrojaracin, shares sequence identity (68-84%) with the α and β chains of factor IX/X binding protein and botrocetin from B. jararaca venom, and inhibits thrombin-induced platelet aggregation. Bothroalternin has also been detected in the proteome of B. alternatus.
Bradykinin-potentiating peptide and C-type natriuretic peptide precursors
Bothrops venoms have long been known to contain bradykinin-potentiating peptides (BPPs) that enhance venom-induced hypotension by inhibiting angiotensin-converting enzyme, a pivotal enzyme in the formation of angiotensin II (vasoconstrictor peptide) and degradation of bradykinin (vasodilator peptide) . Numerous such peptides have since been identified in Bothrops venoms ( and references therein), and cloning and sequence analysis studies in the late 1990s identified genes encoding BPPs in B. jararaca venom gland [108, 109]. These same studies also reported the presence of genes encoding for C-type natriuretic peptides (CNPs) that could contribute to venom-induced hypotension. Transcriptomic studies of Bothrops venom glands have confirmed these findings and shown that genes encoding these two peptides account for 6-20% of venom toxin ESTs [11, 13, 15, 16]. A similar percentage (~9%) was observed here for B. alternatus (Figure 8), for which 110 ESTs coded for BPP/CNP precursors (specifically the CNP region) while another 156 ESTs showed similarities only in the 3'UTR. In contrast to the abundance of BPP/CNP transcripts detected in Bothrops spp. by EST analysis, proteomic analyses have shown that these peptides account for <1% of the venom protein/peptide content [6, 110], although in B. insularis the venom BPP content is 10% . A similar situation applies to CNPs for which some peptidomic analyses have detected no such peptides in venom, despite the identification of transcripts in EST analyses, e.g., B. insularis and L. muta (; but cf. ) (see section on 'Inverted repeats' below for additional comment on CNP detection).
Minor toxin classes
Cysteine-rich secretory proteins
Snake venom cysteine-rich secretory proteins (CRISPs) are 23-26 kDa proteins with a high content of cysteine residues (16) that form eight disulfide bridges. CRISPs are widespread in snake venoms and have been detected in transcriptomic [11, 13–16] and proteomic [6, 7, 11, 96, 110, 112, 113] studies of Bothrops venoms, where they account for 0.5-2% and 0-3.6% of venom genes and toxins, respectively. In agreement with these studies, we detected five ESTs with protein (0.4% of toxin transcripts) for CRISPs in B. alternatus venom gland. Although the role of venom CRISPs remains poorly understood, several members of this family block a variety of ion channels, including L-type Ca2+ channels, cyclic-nucleotide-gated ion channels, voltage-activated K+ channels (Kv1.3), high-conductance Ca2+-activated K+ channels (BCa) and the ryanodine-sensitive intracellular Ca2+ channel. Peichoto et al. have recently shown that a CRISP from the venom of the South American colubrid Philodryas patagoniensis causes myonecrosis. Since myonecrosis is an important local response to envenoming by Bothrops spp., including B. alternatus, it is possible that CRISPs could contribute to this activity.
Dipeptidylpeptidase IV (DPP IV) (also known as CD26) has a wide distribution in snake venoms . DPP IV genes have been identified in the venom glands of B. jararaca, Gloydius blomhoffii brevicaudus, L. muta and various Australian elapids . Gasparello-Clemente and Silveira  reported DPP IV in Bothrops venoms, with the highest activity in B. alternatus. In agreement with this, we detected a transcript coding for DDP IV protein in B. alternatus. Venom DPP IV, which has enzymatic properties similar to those of other eukaryotic DPP IV , may contribute to venom-induced cardiovascular alterations, possibly by degrading endogenous peptides through association with exosome-like vesicles in fresh venom , and possibly by interfering with glucose homeostasis and the immune and neuroendocrine systems . DPP IV may also be important in processing polypeptide precursors of venom peptides, as suggested for wasp and bee venoms [120, 121].
Snake venom growth factors comprise primarily vascular endothelial growth factor (svVEGF) and nerve growth factor (NGF), both of which have been detected in Bothrops ESTs [11, 13–16]. We detected five ESTs for VEGF similar to B. insularis VEGF [13, 122], one EST similar to L. muta VEGF  and one EST similar to C. d. terrificus NGF. However, the proportion of growth factor transcripts in B. alternatus was lower than in other Bothrops ESTs. Junqueira-de-Azevedo et al. showed that svVEGF from B. insularis venom was able to increase vascular permeability and suggested that this protein could be involved in local and systemic vascular responses to envenoming.
Hyaluronidase, which has an important role in facilitating venom diffusion from the site of inoculation through its degradation of hyaluronic acid in the extracellular matrix, is widespread in Bothrops venoms . However, hyaluronidase has not generally been detected in transcriptomic and proteomic analyses of this genus, perhaps because of its low level of expression in venom glands. In agreement with this, we found only one EST similar to truncated hyaluronidase from Bitis arietans in the B. alternatus cDNA library.
L-Amino acid oxidase
L-Amino acid oxidase (LAO), which is widespread in Bothrops venoms , exerts a variety of biological activities, including interference with platelet aggregation, cytotoxicity and microbicidal activity [87, 88]. These deleterious effects are mediated largely via the production of H2O2 during the oxidation of α-keto amino acids. LAO accounts for 0.5-2.6% of toxin transcripts in Bothrops venom gland transcriptomes [13–16] and has been purified and characterized from these venoms . In the case of B. alternatus, we detected seven ESTs for this toxin in the nucleotide region: four of these were related to LAO from B. jararaca and three were ESTs in the 3'UTR related to Ophiophagus hannah LAO. In contrast to this low transcript abundance, LAO accounted for 6.9% of the venom proteins identified in a proteomic analysis of B. alternatus venom . LAO purified from B. alternatus venom is an acidic (pI ~5.37), homodimeric (123 kDa) glycoprotein that induces platelet aggregation, causes edema, is bactericidal and slightly hemorrhagic ; this enzyme may contribute to the cytotoxicity of B. alternatus venom in cultured Madin-Darby canine kidney (MDCK) cells .
Three-finger toxins (3-FTx)
Three-finger toxins (3-FTx) consist predominantly of elapid neurotoxins, including α-neurotoxins, cardiotoxins and fasciculins, and a variety of less well-characterized venom proteins [126, 127]. Originally thought to be restricted to elapids, 3-FTx have since been identified in crotalid [18, 128] and colubrid [127, 129, 130] venom glands. 3-FTx exert a variety of biological activities that include blockade of acetylcholine (nicotinic and muscarinic) receptors, β-adrenergic receptors, L-type calcium channels and integrins, inhibition of acetylcholinesterase, and cardiotoxicity mediated by interation with phospholipids.
We identified 30 ESTs for 3-FTx in the B. alternatus venom gland library (0.56% of total ESTs; 2.4% of toxin ESTs) that shared 90% similarity with the intron II region of the gene for 3FTx-3 from Sistrurus catenatus edwardsi, but no hits with colubrid, elapid or L. muta 3-FTx. Pahari et al. also noted that the 3-FTx nucleotide and protein sequences of S. c. edwardsi bore no relationship to those of L. muta, a finding confirmed by phylogenetic analysis that placed the L. muta toxins distant from those of S. c. edwardsi. Complete sequencing of the B. alternatus genes would be helpful in determining the precise relationship between these 3-FTx and those of S. c. edwardsi.
This is the first identification of 3-FTx genes in Bothrops and, together with other studies, suggests that this class of toxins may occur in a variety of New World pitvipers, i.e., Bothrops [this study], Sistrurus and Lachesis (based on transcriptomic analyses) and Atropoides mexicanus (nummifer) (but not in Atropoides picadoi) (based on proteomic analysis) . However, transcriptomic analyses have not detected these toxins in other New World pitvipers, e.g., Agkistrodon piscivorus leucostoma and C. d. collilineatus, or in Old World pitvipers (Agkistrodon (Deinagkistrodon) acutus[133, 134]) and vipers (Bitis gabonica and Echis species [12, 135]). In addition to this inter-generic variation, there is also intra-generic variation in the occurrence of 3-FTx. Thus, whereas 3-FTx have been detected in a transcriptomic analysis of S. c. edwardsi, these toxins have not been detected in a proteomic analysis of venoms from several members of this genus . Likewise, A. mexicanus (nummifer) venom, but not that of A. picadoi, contains 3-FTx . Such variation may reflect the low abundance of transcripts and proteins (making their detection difficult) and/or the non-uniform recruitment of these toxins into the venom proteome . This could explain why transcriptomic [11, 13–16] and proteomic [6, 7, 11, 96, 110, 112, 113, 137] analyses have generally not detected these genes and proteins in Bothrops species.
The physiological relevance of 3-FTx to envenoming by Bothrops species is unclear, particularly in view of the low transcript abundance indicated by our results. However, Bothrops venoms cause neuromuscular blockade in avian and mammalian nerve-muscle preparations in vitro. Although most of this activity has been attributed to the action of basic myotoxic PLA2, 3-FTx could also contribute to this response through their ability to interfere with neuromuscular transmission.
Envenoming by Bothrops species results in local pain at the bite site [21, 28]. Venom components implicated in this phenomenon include PLA2, metalloproteinases and serine proteinases that act via the release of endogenous mediators such as arachidonic acid metabolites (leucotrienes and prostaglandins), bradykinin, proinflammatory cytokines and neuronal nitric oxide [23, 24, 26]. In contrast, little is known of the involvement of other Bothrops venom proteins in this pain. The identification here of 24 ESTs with hits in the 3'UTR of the precursor for ohanin, an ~12 kDa hypolocomotion and hyperalgesic protein initially identified in the venom of the king cobra, O. hannah, and part of the vespryn family of proteins [139, 140], indicates that other venom proteins, in addition to those indicated above, may be involved in Bothrops venom-induced pain. Since its initial identification, ohanin (vespryn) has been found in Australian elapids [9, 141] and pitvipers such as A. acutus, C. d. collilineatus and L. muta. Our finding is the first report for ohanin in Bothrops venom glands. Together, these reports indicate that ohanins (and vespryns in general) are likely to be widespread in snake venoms where they may contribute to venom-induced hyperalgesia.
We detected two ESTs coding for a protein related to the serine proteinase inhibitor component of taicatoxin (TCX), a multimeric protein composed of an α-neurotoxin-like peptide (8 kDa), a neurotoxin PLA2 (16 kDa) and a serine proteinase inhibitor (7 kDa) (1:1:4 stoichiometry) initially isolated from the venom of the Australian taipan (Oxyuranus scutellatus scutellatus) [142, 143]. TCX blocks high threshold, voltage-dependent calcium channel currents in cardiac membranes and cultured ventricular myocytes, although there may also be a role for the PLA2 activity in the action of this toxin [143, 144]. The role of putative TCX-like proteins in B. alternatus is unclear but could contribute to the cardiovascular actions of this venom.
In addition to the foregoing toxins, we identified a low number of hits for several other toxins with BLAST results only in the nucleotide database. These hits included toxins related to cardiotoxins (12 reads), a muscarinic-like toxin (three reads) related to that of Bungarus multicinctus, precursors similar to post-synaptic α-neurotoxin NTX-2 from Naja sputatrix (one contig) and neurotoxin 6, and a molecule related to cobra-venom factor (one read).
Snake venoms contain a variety of nucleotidases (phosphodiesterases, 5'-nucleotidase, acid and alkaline phosphatases and ADP/ATPases) and nucleases (deoxyribonuclease - DNase and ribonuclease - RNase) that have a potentially important role in envenoming, particularly in affecting platelet aggregation and cardiovascular responses (hypotension, vascular permeability) [119, 145]. In agreement with this, we obtained transcripts coding for a variety of genes related to some of these proteins, including an acidic DNase similar to a mouse DNAse IIα (lysosomal DNase) precursor (3 ESTs), ecto-5'-nucleotidases (one EST coding for the enzyme in zebrafish, Danio rerio, and another transcript coding for the same enzyme in horse, Equus caballus), ectonucleotide pyrophosphatase/phosphodiesterase 3 (two ESTs, one each for monkey, Macaca mulatta, and mouse, Mus musculus) and an adenosine deaminase related to Xenopus laevis enzyme (one transcript). Although several of these genes have been detected in other transcriptomic analyses [13, 18, 61, 133, 135], it is currently unclear whether the corresponding proteins are secreted into the venom or simply part of normal intracellular metabolism in the venom gland. We have purified phosphodiesterase , 5'-nucleotidase and an acidic DNase (DNase II) (unpublished findings) from B. alternatus venom, but since the structure of these proteins is unknown it is unclear to what extent they share similarity with the proteins coded by the genes identified in the venom gland ESTs. Ogawa et al. have recently shown that G. blomhoffi ecto-5'-nucleotidase is structurally similar to this enzyme from various other vertebrates.
Other genes identified that could potentially be involved in envenoming included a cytokine-like protein (one EST similar to horse cytokine) (also identified in other snake venom ESTs [61, 133]) and tumor necrosis factor (two ESTs similar to chicken, Gallus gallus, protein), both of which could contribute to the venom-induced local inflammatory response.
Selected features of B. alternatus venom gland ESTs
Single nucleotide polymorphisms and simple repeats (microsatellites)
Polymorphisms in the nucleic acid sequences of snake venom proteins, particularly PLA2, have been identified in several species [147, 148]. These genetic modifications arise from rapid gene duplication followed by single nucleotide polymorphisms (SNPs), with an increase in non-synonymous nucleotide substitutions that alter the DNA sequence encoding the protein. Ohno et al. suggested that such alterations result in the rapid appearance of novel toxins with different biological activities. Our analysis revealed 132 putative SNPs in the B. alternatus transcriptome, of which 31 (16 non-synonymous and 15 synonymous substitutions) were located in ORFs, as determined based on alignment against the UniProt database. In addition, we identified 32 insertion-deletion polymorphisms (indels) (Table 1; Additional file 1). Although not extensively studied, detailed analysis of venom protein SNPs could be useful for population genetic studies and for assessing the importance of rapid sequence changes in generating the observed diversity of genes involved in venom production .
The screening of unisequences with the tandem repeats finder tool resulted in the identification of 321 sequences with possible microsatellite regions; when annotated sequences with transposable elements (TEs) were excluded, 244 sequences with microsatellite regions were identified (data not shown). These sequences could be potentially useful for the development of molecular markers for characterization of the genetic variability and population structure of B. alternatus throughout its geographic distribution.
Eukaryotic genomes contain a large number of repeated sequences, a high proportion of which may consist of transposable elements (TEs). In snakes, TEs have been previously identified in PLA2 genes from the venom gland of Vipera ammodytes[149, 150]. These TEs are ruminant retroposons corresponding to 5'-truncated Bov-B long interspersed repeated DNA (LINE) and were identified in ammodytin L (a natural mutant of a group II PLA2) and ammodytoxin C (similar structure to other mammalian group II PLA2) genes. Alignments meeting our criteria (see Methods) allowed the identification of repeated elements in 492 unisequences, which corresponded to 9% of the B. alternatus ESTs (Additional file 2); (retro)transposons have also been identified in B. insularis and L. muta. Several TE families were identified among the B. alternatus transcripts, the most abundant being BovB, which was identified in 120 unisequences (Table 2). The expression of these TE families indicated that the genome harbored a large number of potentially active elements that gave rise to a variety of ESTs.
Further examination of the 492 unisequences harboring TEs revealed that 193 of them were fused to a sequence not included in the repeats database. To determine the identity of these sequences, we masked the TEs and BLASTed the masked unisequences against the GeneBank nr database. This approach revealed that in 80 unisequences, the TE was fused to a protein-coding gene (Additional file 3). Most transposition events involving protein-coding regions are deleterious to host genomes because they disrupt mRNA translation, localization and stability. However, since TEs associated with PLA2 genes have been found to be transcribed [136, 137], we hypothesize that the TEs in B. alternatus may act as a source of genetic novelty by playing an important role in the origin and diversity of toxin genes.
The alignment of B. alternatus unisequences to corresponding reference genes indicated that some of the genes contained long inverted repeats (IR) in which a segment of the sequence was inverted from its original orientation in the reference gene (Additional file 4). An interesting example of such inversions involved genes for BPPs/CNPs, for which 16 unisequences with inverted repeats were observed (data not shown). Figure 10 shows the case in which an IR occurred in a portion of mRNA coding for part of a CNP. As mentioned above in the discussion of CNP genes, although EST studies have indicated the presence of genes for these peptides in Bothrops and Lachesis species, proteomic studies have not always been successful in isolating and identifying these peptides [11, 111]. The reason for this discrepancy between EST and proteomic findings has not been adequately addressed but one explanation could be related to the occurrence of IR in these genes that may interfere with adequate gene transcription and/or translation in vivo, thereby preventing toxin production.
The origin of such IR transcripts is unclear. Theoretically, they could be artifacts generated by end-to-end joining of non-contiguous cDNA sequences and template switching by reverse transcriptase during cDNA synthesis; template switching can artificially delete portions of cDNAs and be wrongly interpreted as alternative transcripts. During first strand cDNA synthesis, poly(A)+ RNAs are primed by an oligo(dT) and when the reverse transcriptase reaches the 5'-end of the junction site in the mRNA, the nascent cDNA switches the template to the IR in the antisense mRNA that will now be used as a new template. However, since the homologous recombination between two RNA templates promoted by reverse transcriptase requires RNase H activity, and since the PowerScript™ Reverse Transcriptase used here for cDNA synthesis lacked RNase H activity, it seems unlikely that homology-dependent template switching could have generated these hybrid cDNAs artificially.
Further reasons that the IR observed here are unlikely to be artifacts include the following: 1) end-to-end joining during library construction would randomly ligate molecules and one would expect not only intra-molecular hybrids but also ligation of independently transcribed mRNAs, 2) in the case of an artifact, the transcript of highest abundance would be expected to be frequently ligated to other transcripts to form inter-molecular hybrids, but this was not observed in our library, and 3) the finding that the same ESTs containing inversions were present in multiple sequence reads from the same library and in independent libraries. Together, these observations suggest that this phenomenon was not an artifact.
Comparison of the B. alternatus transcriptome with other Bothrops species
Following the initial report of gene expression analysis in the venom gland of B. insularis (, recently updated by Valente et al.), similar studies were reported for other Bothrops species, including B. atrox, B. jararaca and B. jararacussu. The availability of data for several Bothrops species provides an opportunity to compare related species and draw some general conclusions regarding venom composition in this genus. Few such transcriptomic comparisons are available for other genera .
To date, our study represents the largest EST database for Bothrops; the 5,350 ESTs reported here considerably exceed those reported for B. atrox (610) , B. insularis (initially 610, but updated to 2042) [11, 13], B. jararaca (2318)  and B. jararacussu (549) . We chose to generate a larger database than those previously reported for Bothrops species in an attempt to identify novel genes in addition to the main toxin groups already known for this genus. Indeed, this larger database allowed the detection of genes for toxins such as ohanin and 3-FTx, previously identified in other snake genera but not in Bothrops, in addition to non-venom proteins such as DUSP6 and thioredoxin.
The generation of a large database for B. alternatus did not significantly alter the profile of the most common toxin groups (metalloproteinases, BPPs/CNPs, PLA2, serine proteinases, C-type lectins, growth factors, etc.) when compared with other Bothrops species. This finding suggests that a large database is not essential for identification of the major toxin groups but can be useful for increasing the chances of detecting toxins with very low transcript abundance. This conclusion agrees with data for B. insularis in which increasing the EST database from 610  to 2042  ESTs did not markedly alter the relative proportions of metalloproteinases (41.7% vs. 43.2%, for 610 vs. 2042 ESTs), BPPs/CNPs (19.7% vs. 15.8%), C-type lectins (14.6% vs. 14.2%), serine proteinases (9.6% vs. 11.2%), PLA2 (6.7% vs. 5.4%), svVEGF (4.3% vs. 4.7%), LAO (2.6% vs. 3.5%), CRISPs (0.6% vs. 1.5%) and NGF (0.3% vs. 0.4%). In the present study, 78% of the hits corresponded to toxin transcripts, which compares favorably with values for other Bothrops species (54% - 78%) [11, 13–16] and indicates that the relative proportion of toxin transcripts is not directly related to the size of the cDNA library.
Conversely, the generation of a large database resulted in a considerably greater number of no-hits: 70% in this study compared to 13% - 25% for other Bothrops studies [11, 13–16] and 7% - 54% for other snake genera [12, 18–20, 61, 128, 132–135]. These no-hits reflect the limited amount of information available for Bothrops species and other snakes in venom gland EST databases, and represent a potentially rich source for the identification of novel toxins.
Figure 11 compares the relative abundance of the major toxin classes observed in Bothrops species based on EST analyses. In all cases, there was a predominance of metalloproteinases, BPPs, PLA2, serine proteinases and C-type lectins, with less abundant groups being LAO, CRISPs and growth factors (principally svVEGF and NGF). There was considerable inter-specific variation in the content of the major toxin. Thus, B. alternatus had the highest proportion of metalloproteinase transcripts among the five species, being more than three-fold more abundant than in B. jararacussu. PLA2 abundance was similar to B. insularis, greater than B. jararaca but less than B. atrox and B. jararacussu; the latter species was the only one in which PLA2 transcripts were more abundant than metalloproteinases (at least two-fold greater). The proportion of BPP/CNP transcripts in B. alternatus was similar to B. atrox and B. jararaca but about half that of B. insularis, while serine proteinases and C-type lectins were generally less abundant than in other Bothrops species. As indicated above, a lower content of serine proteinases and C-type lectins in the venom could account for the less severe coagulopathy observed clinically for envenoming by B. alternatus compared to other Bothrops species . Despite the inter-specific variation in the relative proportion of toxin classes, these findings confirm that most Bothrops venom components can be classified into a few major groups. This conclusion agrees with proteomic analyses of Bothrops venoms that have also identified these groups as the major toxin families [6, 7, 11, 90, 96, 97, 110, 112, 113, 137] (Figure 12). In addition to inter-specific variation, these proteomic studies have also reported individual, age-dependent and geographic variation in the toxin content of these major classes [110, 112, 137].
For five Bothrops species (B. alternatus, B. atrox, B. insularis, B. jararaca and B. jararacussu) there are transcriptomic and proteomic analyses that allow comparison of the toxin frequencies in the different classes. For metalloproteinases and PLA2, there is reasonably good agreement between the proportion of transcripts and the corresponding levels of these proteins detected in the venoms, whereas for other classes, e.g., BPPs, C-type lectins and serine proteinases, there are often marked discrepancies between the transcriptomic and proteomic data (cf. Figures 11 and 12). In the case of B. alternatus, there was good agreement between the proportion of ESTs (this study) and venom content  of PLA2 (5.5% vs 7.8%, respectively) and C-type lectins (1.4% vs. 1.7%), but considerable divergence between these two data sets in the case of metalloproteinases (81.4% vs. 43.1%), serine proteinases (1.9% vs. 24.1%) and LAO (0.6% vs. 6.9%) [this study and ref. ]. Divergent transcriptomic and proteomic results have also been observed for certain toxin groups in other snake genera, e.g., Echis species  and L. muta. The causes and implications of such discrepancies have been discussed elsewhere  and indicate the need for caution in interpreting transcriptomic data as being representative of the final venom composition.
Finally, it should be noted that the principal toxin classes identified by transcriptomic and proteomic analyses are those that have been shown to contribute to the major local and systemic effects produced by these venoms, i.e., potent hemorrhagic activity (mediated by metalloproteinases), extensive myonecrosis (myotoxic PLA2), inflammatory responses (metallo- and serine proteinases, PLA2), coagulopathy (C-type lectins, serine proteinases, disintegrins) and cardiovascular actions (metalloproteinases, PLA2, serine proteinases such as kallikrein-like enzymes, and peptides such as BPPs) seen experimentally and clinically [21, 22, 26–30].
An EST database containing 5,350 ESTs was produced for the venom gland of B. alternatus, the largest to date for the genus Bothrops. The venom gland expressed the major toxin groups (metalloproteinases, PLA2, serine proteases, C-type lectins and BPP/CNPs) identified in transcriptomic and proteomic analyses of other Bothrops species. In addition, genes for toxins (ohanin, 3-FTx and taicatoxin-like protein) and non-venom components (Dusp6 and thioredoxin) not previously detected in Bothrops venoms were also observed in this gland. The results described here expand our understanding of the composition of B. alternatus venom and help to explain the major symptoms associated with envenoming by this species.
C-type natriuretic peptide
cysteine-rich secretory protein
- DPP IV:
expressed sequence tags
nerve growth factor
open reading frame
single nucleotide polymorphism
snake venom metalloproteinase
snake venom vascular endothelial growth factor
Ménez A, Stöcklin R, Mebs D: 'Venomics' or: The venomous systems genome project. Toxicon. 2006, 47: 255-259. 10.1016/j.toxicon.2005.12.010.
Calvete JJ, Juárez P, Sanz L: Snake venomics. Strategy and applications. J Mass Spectrom. 2007, 42: 1405-1414. 10.1002/jms.1242.
Fox JW, Serrano SMT: Exploring snake venom proteomes: multifaceted analyses for complex toxin mixtures. Proteomics. 2008, 8: 909-920. 10.1002/pmic.200700777.
Georgieva D, Arni RK, Betzel C: Proteome analysis of snake venom toxins: pharmacological insights. Expert Rev Proteomics. 2008, 5: 787-797. 10.1586/147894184.108.40.2067.
Escoubas P, King GF: Venomics as a drug discovery platform. Expert Rev Proteomics. 2009, 6: 221-224. 10.1586/epr.09.45.
Calvete JJ, Borges A, Segura A, Flores-Díaz M, Alape-Girón A, Gutiérrez JM, Diez N, De Sousa L, Kiriakos D, Sánchez E, Faks JG, Escolano J, Sanz L: Snake venomics and antivenomics of Bothrops colombiensis, a medically important pitviper of the Bothrops atrox-asper complex endemic to Venezuela: contributing to its taxonomy and snakebite management. J Proteomics. 2009, 72: 227-240. 10.1016/j.jprot.2009.01.005.
Gutiérrez JM, Sanz L, Escolano J, Fernández J, Lomonte B, Angulo Y, Rucavado A, Warrell DA, Calvete JJ: Snake venomics of the Lesser Antillean pit vipers Bothrops caribbaeus and Bothrops lanceolatus: correlation with toxicological activities and immunoreactivity of a heterologous antivenom. J Proteome Res. 2008, 7: 4396-4408. 10.1021/pr8003826.
Gutiérrez JM, Lomonte B, León G, Alape-Girón A, Flores-Díaz M, Sanz L, Angulo Y, Calvete JJ: Snake venomics and antivenomics: proteomic tools in the design and control of antivenoms for the treatment of snakebite envenoming. J Proteomics. 2009, 72: 165-182. 10.1016/j.jprot.2009.01.008.
St. Pierre L, Birrell GW, Earl ST, Wallis TP, Gorman JJ, de Jersey J, Masci PP, Lavin MF: Diversity of toxic components from the venom of the evolutionarily distinct black whip snake, Demansia vestigiata. J Proteome Res. 2007, 6: 3093-3107. 10.1021/pr0701613.
Sanz L, Escolano J, Ferretti M, Biscoglio MJ, Rivera E, Crescenti EJ, Angulo Y, Lomonte B, Gutiérrez JM, Calvete JJ: Snake venomics of the South and Central American bushmasters. Comparison of the toxin composition of Lachesis muta gathered from proteomic versus transcriptomic analysis. J Proteomics. 2008, 71: 46-60. 10.1016/j.jprot.2007.10.004.
Valente RH, Guimarães PR, Junqueira M, Neves-Ferreira AG, Soares MR, Chapeaurouge A, Trugilho MR, León IR, Rocha SL, Oliveira-Carvalho AL, Wermelinger LS, Dutra DL, Leão LI, Junqueira-de-Azevedo IL, Ho PL, Zingali RB, Perales J, Domont GB: Bothrops insularis venomics: a proteomic analysis supported by transcriptomic-generated sequence data. J Proteomics. 2009, 72: 241-255. 10.1016/j.jprot.2009.01.001.
Wagstaff SC, Sanz L, Juárez P, Harrison RA, Calvete JJ: Combined snake venomics and venom gland transcriptomic analysis of the ocellated carpet viper, Echis ocellatus. J Proteomics. 2009, 71: 609-623. 10.1016/j.jprot.2008.10.003.
Junqueira-de-Azevedo ILM, Ho PL: A survey of gene expression and diversity in the venom glands of the pit viper snake Bothrops insularis through the generation of expressed sequence tags (ESTs). Gene. 2002, 299: 279-291. 10.1016/S0378-1119(02)01080-6.
Kashima S, Roberto PG, Soares AM, Astolfi-Filho S, Pereira JO, Giuliati S, Faria M, Xavier MAS, Fontes MRM, Giglio JR, França SC: Analysis of Bothrops jararacussu venomous gland transcriptome focusing on structural and functional aspects: 1 - gene expression profile of highly expressed phopholipases A2. Biochimie. 2004, 86: 211-239. 10.1016/j.biochi.2004.02.002.
Cidade DAP, Simão TA, Dávila AMR, Wagner G, Junqueira-de-Azevedo ILM, Ho PL, Bon C, Zingali RB, Albano RM: Bothrops jararaca venom gland transcriptome: analysis of the gene expression pattern. Toxicon. 2006, 48: 437-461. 10.1016/j.toxicon.2006.07.008.
Neiva M, Arraes FB, de Souza JV, Rádis-Baptista G, Prieto da Silva AR, Walter ME, de Macedo Brigido M, Yamane T, López-Lozano JL, Astolfi-Filho S: Transcriptome analysis of the Amazonian viper Bothrops atrox venom gland using expressed sequence tags (ESTs). Toxicon. 2009, 53: 427-436. 10.1016/j.toxicon.2009.01.006.
Boldrini-França J, Rodrigues RS, Fonseca FP, Menaldo DL, Ferreira FB, Henrique-Silva F, Soares AM, Hamaguchi A, Rodrigues VM, Otaviano AR, Homsi-Brandeburgo MI: Crotalus durissus collilineatus venom gland transcriptome: analysis of gene expression profile. Biochimie. 2009, 91: 586-595. 10.1016/j.biochi.2009.02.001.
Junqueira-de-Azevedo ILM, Ching ATC, Carvalho E, Faria F, Nishiyama MY, Ho PL, Diniz MRV: Lachesis muta (Viperidae) cDNAs reveal diverging pit viper molecules and scaffolds typical of cobra (Elapidae) venoms: implications for snake toxin repertoire evolution. Genetics. 2006, 172: 877-889. 10.1534/genetics.106.056515.
Leão LI, Ho PL, Junqueira-de-Azevedo ILM: Transcriptomic basis for an antiserum against Micrurus corallinus (coral snake) venom. BMC Genomics. 2009, 10: 112-10.1186/1471-2164-10-112.
Ching ATC, Rocha MMT, Paes Leme AF, Pimenta DC, Furtado MFD, Serrano SMT, Ho PL, Junqueira-de-Azevedo ILM: Some aspects of the venom proteome of the Colubirdae snake Philodryas olfersii revealed from a Duvernoy's (venom) gland transcriptome. FEBS Lett. 2006, 580: 4417-4422. 10.1016/j.febslet.2006.07.010.
Warrell DA: Snakebites in Central and South America: epidemiology, clinical features and clinical management. Venomous Reptiles of the Western Hemisphere. Edited by: Campbell JA, Lamar WW. 2004, Ithaca: Comstock Publishing Associates/Cornell University Press, 2: 709-761.
Gutiérrez JM, Lomonte B: Efectos locales en el envenenamiento ofídico en América Latina. Animais Peçonhentos no Brasil: Biologia, Clínica e Terapêutica dos Acidentes. Edited by: Cardoso JLC, França FOS, Wen FH, Málaque CMS, Haddad V Jr. 2003, São Paulo: Sarvier/FAPESP, 310-323.
Teixeira CF, Landucci EC, Antunes E, Chacur M, Cury Y: Inflammatory effects of snake venom myotoxic phospholipase A2. Toxicon. 2003, 42: 947-962. 10.1016/j.toxicon.2003.11.006.
Teixeira C, Cury Y, Moreira V, Picolo G, Chaves F: Inflammation induced by Bothrops asper venom. Toxicon. 2009, 54: 988-997. 10.1016/j.toxicon.2009.05.026.
Pinho FM, Yu L, Burdmann EA: Snakebite-induced acute kidney injury in Latin America. Semin Nephrol. 2008, 28: 354-362. 10.1016/j.semnephrol.2008.04.004.
Gutiérrez JM, Rucavado A, Chaves F, Díaz C, Escalante T: Experimental pathology of local tissue damage induced by Bothrops asper snake venom. Toxicon. 2009, 54: 958-975. 10.1016/j.toxicon.2009.01.038.
Gutiérrez JM, Escalante T, Rucavado A: Experimental pathophysiology of systemic alterations induced by Bothrops asper snake venom. Toxicon. 2009, 54: 976-987. 10.1016/j.toxicon.2009.01.039.
França FOS, Málaque CMS: Acidente botrópico. Animais Peçonhentos no Brasil: Biologia, Clínica e Terapêutica dos Acidentes. Edited by: Cardoso JLC, França FOS, Wen FH, Málaque CMS, Haddad V Jr. 2003, São Paulo: Sarvier/FAPESP, 72-86.
Gutiérrez JM, Rucavado A, Escalante T, Díaz C: Hemorrhage induced by snake venom metalloproteinases: biochemical and biophysical mechanisms involved in microvessel damage. Toxicon. 2005, 45: 997-1011. 10.1016/j.toxicon.2005.02.029.
Gutiérrez JM, Ownby CL: Skeletal muscle degeneration induced by venom phospholipases A2: insights into the mechanisms of local and systemic myotoxicity. Toxicon. 2003, 42: 915-931. 10.1016/j.toxicon.2003.11.005.
Montecucco C, Gutiérrez JM, Lomonte B: Cellular pathology induced by snake venom phospholipase A2 myotoxins and neurotoxins: common aspects of their mechanisms of action. Cell Mol Life Sci. 2008, 65: 2897-2912. 10.1007/s00018-008-8113-3.
Campbell JA, Lamar WW: Venomous Reptiles of the Western Hemisphere. 2004, Ithaca: Comstock Publishing Associates/Cornell University Press
Fonseca F: Animais Peçonhentos. 1949, São Paulo: Instituto Butantan
Ribeiro LA, Jorge MT: Acidente por serpentes do gênero Bothrops: série de 3.139 casos. Rev Soc Bras Med Trop. 1997, 30: 475-480.
Ruiz de Torrent RM, Ulón SN, Sario HM: Epidemiológia del accidente ofídico en la provincia de Corrientes. 2004, Corrientes, Argentina: Universidad Nacional del Nordeste - Comunicaciones Científicas y Tecnológicas, Abstract V-029
Bauab FA, Junqueira GR, Corradini MC, Silveira PV, Nishioka AS: Clinical and epidemiological aspects of the 'urutu' lance-headed viper (Bothrops alternatus) bite in a Brazilian hospital. Trop Med Parasitol. 1994, 45: 243-245.
Purtscher H, Burger M, Savio E, Juanicó LR, Domínguez BV, Caorsi E, Tambler M, Laborde A: Ofidismo y aracnidismo en el Uruguay. Diagnostico, complicaciones y tratamiento. Rev Med Uruguay. 1983, 7: 1-37.
Cardoso DF, Yamaguchi IK, Moura-da-Silva AM: Produção de soros antitoxinas e perspectivas de modernização por técnicas de biologia molecular. Animais Peçonhentos no Brasil: Biologia, Clínica e Terapêutica dos Acidentes. Edited by: Cardoso JLC, França FOS, Wen FH, Málaque CMS, Haddad V Jr. 2003, São Paulo: Sarvier/FAPESP, 367-379.
de Salud Ministerio: Resolución 34/2007 - Guia de prevención, diagnóstico, tratamiento y vigilancia epidemiológica de los envenenamientos ofídicos. 2007, Buenos Aires: Dirección de la Calidad de los Servicios de Salud, Ministerio da Salud
Souza DH, Iemma MR, Ferreira LL, Faria JP, Oliva ML, Zingali RB, Niewiarowski S, Selistre-de-Araujo HS: The disintegrin-like domain of the snake venom metalloprotease alternagin inhibits α2β1 integrin-mediated cell adhesion. Arch Biochem Biophys. 2000, 384: 341-350. 10.1006/abbi.2000.2120.
Mariano-Oliveira A, Coelho AL, Terruggi CH, Selistre-de-Araujo HS, Barja-Fidalgo C, de Freitas MS: Alternagin-C, a non-RGD-disintegrin, induces neutrophil migration via integrin signaling. Eur J Biochem. 2003, 270: 4799-4808. 10.1046/j.1432-1033.2003.03867.x.
Cominetti MR, Ribeiro JU, Fox JW, Selistre-de-Araujo HS: BaG, a new dimeric metalloproteinase/disintegrin from the Bothrops alternatus snake venom that interacts with α5β1 integrin. Arch Biochem Biophys. 2003, 416: 171-179. 10.1016/S0003-9861(03)00298-4.
Cominetti MR, Terruggi CH, Ramos OH, Fox JW, Mariano-Oliveira A, De Freitas MS, Figueiredo CC, Morandi V, Selistre-de-Araujo HS: Alternagin-C, a disintegrin-like protein, induces vascular endothelial cell growth factor (VEGF) expression and endothelial cell proliferation in vitro. J Biol Chem. 2004, 279: 18247-18255. 10.1074/jbc.M311771200.
Gay CC, Leiva LC, Maruñak S, Teibler P, Acosta de Pérez O: Proteolytic, edematogenic and myotoxic activities of a hemorrhagic metalloproteinase isolated from Bothrops alternatus venom. Toxicon. 2005, 46: 546-554. 10.1016/j.toxicon.2005.06.019.
Gay CC, Maruñak SL, Teibler P, Ruiz R, Acosta de Pérez OC, Leiva LC: Systemic alterations induced by a Bothrops alternatus hemorrhagic metalloproteinase (baltergin) in mice. Toxicon. 2009, 53: 53-59. 10.1016/j.toxicon.2008.10.010.
Selistre-de-Araujo HS, Cominetti MR, Terruggi CH, Mariano-Oliveira A, de Freitas MS, Crépin M, Figueiredo CC, Morandi V: Alternagin-C, a disintegrin-like protein from the venom of Bothrops alternatus, modulates α2β1 integrin-mediated cell adhesion, migration and proliferation. Braz J Med Biol Res. 2005, 38: 1505-1511. 10.1590/S0100-879X2005001000007.
Ramos OH, Terruggi CH, Ribeiro JU, Cominetti MR, Figueiredo CC, Bérard M, Crépin M, Morandi V, Selistre-de-Araujo HS: Modulation of in vitro and in vivo angiogenesis by alternagin-C, a disintegrin-like protein from Bothrops alternatus snake venom and by a peptide derived from its sequence. Arch Biochem Biophys. 2007, 46: 1-6. 10.1016/j.abb.2007.02.021.
Durigan JL, Peviani SM, Russo TL, Delfino GB, Ribeiro JU, Cominetti MR, Selistre-de-Araujo HS, Salvini TF: Effects of alternagin-C from Bothrops alternatus on gene expression and activity of metalloproteinases in regenerating skeletal muscle. Toxicon. 2008, 52: 687-694. 10.1016/j.toxicon.2008.07.018.
Ramos OH, Kauskot A, Cominetti MR, Bechyne I, Salla Pontes CL, Chareyre F, Manent J, Vassy R, Giovannini M, Legrand C, Selistre-de-Araujo HS, Crépin M, Bonnefoy A: A novel αvβ3-blocking disintegrin containing the RDG motive, DisBa-01, inhibits bFGF-induced angiogenesis and melanoma metastasis. Clin Exp Mestastasis. 2008, 25: 53-64. 10.1007/s10585-007-9101-y.
Kauskot A, Cominetti MR, Ramos OH, Bechyne I, Renard JM, Hoylaerts MF, Crépin M, Legrand C, Selistre-de-Araujo HS, Bonnefoy A: Hemostatic effects of recombinant DisBa-01, a disintegrin from Bothrops alternatus. Front Biosci. 2008, 13: 6604-6616. 10.2741/3176.
Nisenbom HE, Seki C, Vidal JC: Phospholipase A2 from Bothrops alternatus (víbora de la cruz) venom. Purification and some characteristic properties. Toxicon. 1986, 24: 259-272. 10.1016/0041-0101(86)90151-0.
Nisenbom HE, Perazzo JC, Monserrat AJ, Vidal JC: Contribution of phospholipase A2 to the lethal potency of Bothrops alternatus (víbora de la cruz) venom. Toxicon. 1986, 24: 807-817. 10.1016/0041-0101(86)90106-6.
Nisenbom HE, Perazzo JC, Monserrat AJ, Vidal JC: Effect of chemical modification with p-bromophenacyl bromide on the enzymatic and lethal properties of phospholipase A2 from Bothrops alternatus (víbora de la cruz) venom. Toxicon. 1988, 26: 1137-1144. 10.1016/0041-0101(88)90298-X.
Ponce-Soto LA, Lomonte B, Gutiérrez JM, Rodrigues-Simioni L, Novello JC, Marangoni S: Structural and functional properties of BaTX, a new Lys49 phospholipase A2 homologue isolated from the venom of the snake Bothrops alternatus. Biochim Biophys Acta. 2007, 1770: 585-593.
Ponce-Soto LA, Barros JC, Marangoni S, Hernandez S, Dal Belo CA, Corrado AP, Hyslop S, Rodrigues-Simioni L: Neuromuscular activity of BaTX, a presynaptic basic PLA2 isolated from Bothrops alternatus snake venom. Comp Biochem Physiol C Toxicol Pharmacol. 2009, 150: 291-297. 10.1016/j.cbpc.2009.05.007.
Castro HC, Dutra DL, Oliveira-Carvalho AL, Zingali RB: Bothroalternin, a thrombin inhibitor from the venom of Bothrops alternatus. Toxicon. 1998, 36: 1903-1912. 10.1016/S0041-0101(98)00111-1.
Smolka MB, Marangoni S, Oliveira B, Novello JC: Purification and partial characterization of a thrombin-like enzyme, balterobin, from the venom of Bothrops alternatus. Toxicon. 1998, 36: 1059-1063. 10.1016/S0041-0101(98)80008-1.
Vitorino-Cardoso AF, Pereira Ramos OH, Homsi-Brandeburgo MI, Selistre-de-Araujo HS: Insights into the substrate specificity of a novel snake venom serine peptidase by molecular modeling. Comp Biochem Physiol B Biochem Mol Biol. 2006, 144: 334-342. 10.1016/j.cbpb.2006.03.010.
Stábeli RG, Marcussi S, Carlos GB, Pietro RC, Selistre-de-Araujo HS, Giglio JR, Oliveira EB, Soares AM: Platelet aggregation and antibacterial effects of an L-amino acid oxidase purified from Bothrops alternatus snake venom. Bioorg Med Chem. 2004, 12: 2881-2886. 10.1016/j.bmc.2004.03.049.
Valério AA, Corradini AC, Panunto PC, Mello SM, Hyslop S: Purification and characterization of a phosphodiesterase from Bothrops alternatus snake venom. J Protein Chem. 2002, 21: 495-503. 10.1023/A:1022414503995.
Francischetti IMB, My-Pham V, Harrison J, Garfield MK, Ribeiro JMC: Bitis gabonica (gaboon viper) snake venom gland: toward a catalog of the full-length transcripts (cDNA) and proteins. Gene. 2004, 337: 55-69. 10.1016/j.gene.2004.03.024.
Wagstaff SC, Harrison RA: Venom gland EST analysis of the saw-scaled viper, Echis ocellatus, reveals novel α9β1 integrin-binding motifs in venom metalloproteinases and a new group of putative toxins, renin-like aspartic proteases. Gene. 2006, 377: 21-32. 10.1016/j.gene.2006.03.008.
Kochva E: Oral glands of the Reptilia. Biology of the Reptilia. Physiology B. Edited by: Gans C, Gans KA. 1978, London: Academic Press, 8: 43-161.
Carneiro SM, Pinto VR, Jared C, Lula LA, Faria FP, Sesso A: Morphometric studies on venom secretory cells from Bothrops jararacussu (jararacuçu) before and after venom extraction. Toxicon. 1991, 29: 569-580. 10.1016/0041-0101(91)90051-R.
Yamanouye N, Britto LR, Carneiro SM, Markus RP: Control of venom production and secretion by sympathetic outflow in the snake Bothrops jararaca. J Exp Biol. 1997, 200: 2457-2556.
Green P: PHRAP Documentation. 1999, University of Washington, Seattle, [http://bozeman.mbt.washington.edu/phrap.docs/phrap.html]
Baudet C, Dias Z: Analysis of slipped sequences in EST projects. Genet Mol Res. 2006, 5: 169-181.
Huang X, Madan A: CAP3: a DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nuc Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.
Thompson JD, Gibson TJ, Higgins DG: Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics. 2002, 2: 2.3.1-2.3.22.
Tang J, Vosman B, Voorrips RE, van der Linden CG, Leunissen JA: QualitySNP: a pipeline for detecting single nucleotide polymorphisms and insertions/deletions in EST data from diploid and polyploid species. BMC Bioinformatics. 2006, 7: 438-10.1186/1471-2105-7-438.
Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, Martin MJ, Natale DA, O'Donovan C, Redaschi N, Yeh LS: UniProt: the Universal Protein knowledgebase. Nucleic Acids Res. 2004, 32: D115-D119. 10.1093/nar/gkh131.
Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genom Res. 2005, 110: 462-467. 10.1159/000084979.
Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16: 276-277. 10.1016/S0168-9525(00)02024-2.
Yamanouye N, Carneiro SM, Scrivano CN, Markus RP: Characterization of β-adrenoceptors responsible for venom production in the venom gland of the snake Bothrops jararaca. Life Sci. 2000, 67: 217-226. 10.1016/S0024-3205(00)00626-3.
Kerchove CM, Carneiro SM, Markus RP, Yamanouye N: Stimulation of the α-adrenoceptor triggers the venom production cycle in the venom gland of Bothrops jararaca. J Exp Biol. 2004, 207: 411-416. 10.1242/jeb.00778.
Kerchove CM, Luna MS, Zablith MB, Lazari MF, Smaili SS, Yamanouye N: α1-Adrenoceptors trigger the snake venom production cycle in secretory cells by activating phosphatidylinositol 4,5-bisphosphate hydrolysis and ERK signaling pathway. Comp Biochem Physiol A Mol Integr Physiol. 2008, 150: 431-437. 10.1016/j.cbpa.2008.04.607.
Borgheresi RA, Leroy JM, Yogi A, Dos Santos RA, Breno MC, Tostes RC: Pharmacologic and molecular characterization of the vascular ETA receptor in the venomous snake Bothrops jararaca. Exp Biol Med (Maywood). 2006, 231: 729-735.
Junqueira-de-Azevedo IL, Pertinhez T, Spisni A, Carreño FR, Farah CS, Ho PL: Cloning and expression of calglandulin, a new EF-hand protein from the venom glands of Bothrops insularis snake in E. coli. Biochim Biophys Acta. 2003, 1648: 90-98.
St. Pierre L, Woods R, Earl S, Masci PP, Lavin MF: Identification and analysis of venom gland-specific genes from the coastal taipan (Oxyranus scutellatus) and related species. Cell Mol Life Sci. 2005, 62: 2679-2693. 10.1007/s00018-005-5384-9.
Jeffrey KL, Camps M, Rommel C, Mackay CR: Targeting dual-specificity phosphatases: manipulating MAP kinase signaling and immune responses. Nat Rev Drug Discov. 2007, 6: 391-403. 10.1038/nrd2289.
Goldbeter A, Pourquié O: Modeling the segmentation clock as a network of coupled oscillations in the Notch, Wnt and FGF signaling pathways. J Theor Biol. 2008, 252: 574-585. 10.1016/j.jtbi.2008.01.006.
Gomez C, Ozbudak EM, Wunderlich J, Baumann D, Lewis J, Pourquié O: Control of segment number in vertebrate embryos. Nature. 2008, 454: 335-339. 10.1038/nature07020.
Arnér ES, Holmgren A: Physiological functions of thioredoxin and thioredoxin reductase. Eur J Biochem. 2000, 267: 6102-6109. 10.1046/j.1432-1327.2000.01701.x.
Arnér ES: Focus on mammalian thioredoxin reductases - important selenoproteins with versatile functions. Biochim Biophys Acta. 2009, 1790: 495-526.
Du XY, Clemetson KJ: Snake venom L-amino acid oxidase. Toxicon. 2002, 40: 659-665. 10.1016/S0041-0101(02)00102-2.
Zuliani JP, Kayano AM, Zaqueo KD, Neto AC, Sampaio SV, Soares AM, Stábeli RG: Snake venom L-amino acid oxidases: some consideration about their functional characterization. Protein Pept Lett. 2009, 16: 908-912. 10.2174/092986609788923347.
Fox JW, Serrano SMT: Insights into and speculations about snake venom metalloproteinase (SVMP) synthesis, folding and disulfide bond formation and their contribution to venom complexity. FEBS J. 2008, 275: 3016-3030. 10.1111/j.1742-4658.2008.06466.x.
Öhler M, Georgieva D, Seifert J, von Bergen M, Arni RK, Genov N, Betzel C: The venomics of Bothrops alternatus is a pool of acidic proteins with predominant hemorrhagic and coagulopathic activities. J Proteome Res. 2010, 9: 2442-2437. 10.1021/pr901128x.
Gallacci M, Cavalcante WLG: Understanding the in vitro neuromuscular activity of snake venom Lys49 phospholipase A2 homologues. Toxicon. 2010, 55: 1-11. 10.1016/j.toxicon.2009.10.025.
Moura-da-Silva AM, Cardoso DF, Tanizaki MM: Differences in distribution of myotoxic proteins in venoms from different Bothrops species. Toxicon. 1990, 28: 1293-1301. 10.1016/0041-0101(90)90094-N.
Moura-da-Silva AM, Desmond H, Laing G, Theakston RDG: Isolation and comparison of myotoxins isolated from venoms of different species of Bothrops snakes. Toxicon. 1991, 29: 713-723. 10.1016/0041-0101(91)90063-W.
Moura-da-Silva AM, Cardoso DF, Tanizaki MM, Mota I: Neutralization of myotoxic activity of Bothrops venoms by antisera to purified myotoxins and to crude venoms. Toxicon. 1991, 29: 1471-1480. 10.1016/0041-0101(91)90003-A.
Andrião-Escarso SH, Soares AM, Rodrigues VM, Angulo Y, Dias C, Lomonte B, Gutiérrez JM, Giglio JR: Myotoxic phospholipases A2 in Bothrops snake venoms: effect of chemical modifications on the enzymatic and pharmacological properties of bothropstoxins from Bothrops jararacussu. Biochimie. 2000, 82: 755-763. 10.1016/S0300-9084(00)01150-0.
Tashima AK, Sanz L, Camargo ACM, Serrano SMT, Calvete JJ: Snake venomics of the Brazilian pitvipers Bothrops cotiara and Bothrops fonsecai. Identification of taxonomy markers. J Proteomics. 2008, 71: 473-485. 10.1016/j.jprot.2008.07.007.
Correa-Netto C, Teixeira-Araujo R, Aguiar AS, Melgarejo AR, De-Simone SG, Soares MR, Foguel D, Zingali RB: Immunome and venome of Bothrops jararacussu: a proteomic approach to study the molecular immunology of snake toxins. Toxicon. 2010, 55: 1222-1235. 10.1016/j.toxicon.2009.12.018.
Cogo JC, Lilla S, Souza GHMF, Hyslop S, de Nucci G: Purification, sequencing and structural analysis of two acidic phospholipases A2 from the venom of Bothrops insularis (jararaca ilhoa). Biochimie. 2006, 88: 1947-1959. 10.1016/j.biochi.2006.10.006.
Andrião-Escarso SH, Soares AM, Fontes MR, Fuly AL, Corrêa FM, Rosa JC, Greene LJ, Giglio JR: Structural and functional characterization of an acidic platelet aggregation inhibitor and hypotensive phospholipase A2 from Bothrops jararacussu snake venom. Biochem Pharmacol. 2002, 64: 723-732. 10.1016/S0006-2952(02)01210-8.
Estevão-Costa MI, Rocha BC, Mudado MA, Redondo R, Franco GR, Fortes-Dias CL: Prospection, structural analysis and phylogenetic relationships of endogenous γ-phospholipase A2 inhibitors in Brazilian Bothrops snakes (Viperidae, Crotalinae). Toxicon. 2008, 52: 122-129. 10.1016/j.toxicon.2008.04.167.
Zaganelli GL, Zaganelli MG, Magalhães A, Diniz CR, de Lima ME: Purification and characterization of a fibrinogen-clotting enzyme from the venom of jararacuçu (Bothrops jararacussu). Toxicon. 1996, 34: 807-819. 10.1016/0041-0101(96)00006-2.
Serrano SM, Hagiwara Y, Murayama N, Higuchi S, Mentele R, Sampaio CA, Camargo AC, Fink E: Purification and characterization of a kinin-releasing and fibrinogen-clotting serine proteinase (KN-BJ) from the venom of Bothrops jararaca, and molecular cloning and sequence analysis of its cDNA. Eur J Biochem. 1998, 251: 845-853. 10.1046/j.1432-1327.1998.2510845.x.
Fox JW, Ma L, Nelson K, Sherman NE, Serrano SMT: Comparison of indirect and direct approaches using ion-trap and Fourier transform ion cyclotron resonance mass spectrometry for exploring viperid venom proteomes. Toxicon. 2006, 47: 700-714. 10.1016/j.toxicon.2006.01.022.
Zingali RB, Jandrot-Perrus M, Guillin M, Bon C: Bothrojaracin, a new thrombin inhibitor isolated from Bothrops jararaca venom: characterization and mechanism of thrombin inhibition. Biochemistry. 1993, 32: 10794-10802. 10.1021/bi00091a034.
Arocas V, Zingali RB, Guillin MC, Bon C, Jandrot-Perrus M: Bothrojaracin: a potent two-site-directed thrombin inhibitor. Biochemistry. 1996, 35: 9083-9089. 10.1021/bi960043l.
Ferreira SH: A bradykinin-potentiating factor (BPF) present in the venom of Bothrops jararaca. Br J Pharmacol Chemother. 1965, 24: 163-169.
Menin L, Perchuć A, Favreau P, Perret F, Michalet S, Schöni R, Wilmer M, Stöcklin R: High throughput screening of bradykinin-potentiating peptides in Bothrops moojeni snake venom using precursor ion mass spectrometry. Toxicon. 2008, 51: 1288-1302. 10.1016/j.toxicon.2008.02.019.
Murayama N, Hayashi MA, Ohi H, Ferreira LA, Hermann VV, Saito H, Fujita Y, Higuchi S, Fernandes BL, Yamane T, de Camargo AC: Cloning and sequence analysis of a Bothrops jararaca cDNA encoding a precursor of seven bradykinin-potentiating peptides and a C-type natriuretic peptide. Proc Natl Acad Sci USA. 1997, 94: 1189-1193. 10.1073/pnas.94.4.1189.
Higuchi S, Murayama N, Saguchi K, Ohi H, Fujita Y, Camargo AC, Ogawa T, Deshimaru M, Ohno M: Bradykinin-potentiating peptides and C-type natriuretic peptides from snake venom. Immunopharmacology. 1999, 44: 129-135. 10.1016/S0162-3109(99)00119-8.
Núñez V, Cid P, Zanz L, de la Torre P, Angulo Y, Lomonte B, Gutiérrez JM, Calvete JJ: Snake venomics and antivenomics of Bothrops atrox venoms from Colombia and the Amazonian regions of Brazil, Perú and Ecuador suggest the occurrence of geographic variation of venom phenotype by a trend towards paedomorphism. J Proteomics. 2009, 73: 57-78. 10.1016/j.jprot.2009.07.013.
Soares MR, Oliveira-Carvalho AL, Wermelinger LS, Zingali RB, Ho PL, Junqueira-de-Azevedo IL, Diniz MRV: Identification of novel bradykinin-potentiating peptides and C-type natriuretic peptide from Lachesis muta venom. Toxicon. 2005, 46: 31-38. 10.1016/j.toxicon.2005.03.006.
Alape-Girón A, Sanz L, Escolano J, Flores-Díaz M, Madrigal M, Sasa M, Calvete JJ: Snake venomics of the lancehead pitviper Bothrops asper: geographic, individual, and ontogentic variation. J Proteome Res. 2008, 7: 3556-3571. 10.1021/pr800332p.
Alape-Girón A, Flores-Díaz M, Sanz L, Madrigal M, Escolano J, Sasa M, Calvete JJ: Studies on the venom proteome of Bothrops asper: perspectives and applications. Toxicon. 2009, 54: 938-948. 10.1016/j.toxicon.2009.06.011.
Peichoto ME, Mackessy SP, Teibler P, Tavares FL, Burckhardt PL, Breno MC, Acosta O, Santoro ML: Purification and characterization of a cysteine-rich protein from Philodryas patagoniensis snake venom. Comp Biochem Physiol C Toxicol Pharmacol. 2009, 150: 79-84. 10.1016/j.cbpc.2009.03.002.
Aird SD: Snake venom dipeptidyl peptidase IV: taxonomic distribution and quantitative variation. Comp Biochem Physiol B Biochem Mol Biol. 2008, 150: 222-228. 10.1016/j.cbpb.2008.03.005.
Ogawa Y, Mamura Y, Murayama N, Yanoshita R: Characterization and cDNA cloning of dipeptidyl peptidase IV from the venom of Gloydius blomhoffi brevicaudus. Comp Biochem Physiol B Biochem Mol Biol. 2006, 145: 35-42. 10.1016/j.cbpb.2006.05.013.
Gasparello-Clemente E, Silveira PF: Fluorometric assay using naphthylamide substrates for assessing novel venom peptidase activities. Toxicon. 2002, 40: 1617-1626. 10.1016/S0041-0101(02)00180-0.
Ogawa Y, Kanai-Azuma M, Akimoto Y, Kawakami H, Yanoshita R: Exosome-like vesicles in Gloydius blomhoffii blomhoffii venom. Toxicon. 2008, 51: 984-993. 10.1016/j.toxicon.2008.02.003.
Aird SD: Ophidian envenomation strategies and the role of purines. Toxicon. 2002, 40: 335-393. 10.1016/S0041-0101(01)00232-X.
Kreil G, Haiml L, Suchanek G: Stepwise cleavage of the pro part of promelittin by dipeptidylpeptidase IV. Evidence for a new type of precursor-product conversion. Eur J Biochem. 1980, 111: 49-58. 10.1111/j.1432-1033.1980.tb06073.x.
Lee VA, Tu WC, Jinn TR, Peng CC, Lin LJ, Tzen JT: Molecular cloning of the precursor polypeptide of mastoparan B and its putative processing enzyme, dipeptidyl peptidase IV, from the black-bellied hornet, Vespa basilis. Insect Mol Biol. 2007, 16: 231-237. 10.1111/j.1365-2583.2006.00718.x.
Junqueira-de-Azevedo IL, Farsky SH, Oliveira ML, Ho PL: Molecular cloning and expression of a functional snake venom vascular endothelium growth factor (VEGF) from the Bothrops insularis pit viper. A new member of the VEGF family of proteins. J Biol Chem. 2001, 276: 39836-39842. 10.1074/jbc.M106531200.
Queiroz GP, Pessoa LA, Portaro FC, Furtado MF, Tambourgi DV: Interspecific variation in venom composition and toxicity of Brazilian snakes from Bothrops genus. Toxicon. 2008, 52: 842-851. 10.1016/j.toxicon.2008.10.002.
Pessatti M, Fontana JD, Furtado MF, Guimãraes MF, Zanette LR, Costa WT, Baron M: Screening of Bothrops snake venoms for L-amino acid oxidase activity. Appl Biochem Biotechnol. 1995, 51-52: 197-210. 10.1007/BF02933424.
Nascimento JM, Franchi GC, Nowill AE, Collares-Buzato CB, Hyslop S: Cytoskeletal rearrangement and cell death induced by Bothrops alternatus snake venom in cultured Madin-Darby canine kidney cells. Biochem Cell Biol. 2007, 85: 591-605. 10.1139/O07-067.
Fry BG, Wüster W, Kini RM, Brusic V, Khan A, Venkataraman D, Rooney AP: Molecular evolution and phylogeny of elapid snake venom three-finger toxins. J Mol Evol. 2003, 57: 110-129. 10.1007/s00239-003-2461-2.
Fry BG, Scheib H, van der Weerd L, Young B, McNaughtan J, Ramjan SF, Vidal N, Poelmann RE, Norman JA: Evolution of an arsenal: structural and functional diversification of the venom system in the advanced snakes (Caenophidia). Mol Cell Proteomics. 2008, 7: 215-246.
Pahari S, Mackessy SP, Kini RM: The venom gland transcriptome of the dessert massasauga rattlesnake (Sistrurus catenatus edwardsii): towards an understanding of venom composition among advanced snakes (Superfamily Colubroidea). BMC Mol Biol. 2007, 8: 115-10.1186/1471-2199-8-115.
Pawlak J, Kini RM: Unique gene organization of colubrid three-finger toxins: complete cDNA and gene sequences of denmotoxin, a bird-specific toxin from colubrid snake Boiga dendrophila (mangrove catsnake). Biochimie. 2008, 90: 868-877. 10.1016/j.biochi.2008.02.016.
Pawlak J, Mackessy SP, Sixberry NM, Stura EA, Le Du MH, Ménez R, Foo CS, Ménez A, Nirthanan S, Kini RM: Irditoxin, a novel covalently linked heterodimeric three-finger toxin with high taxon-specific neurotoxicity. FASEB J. 2009, 23: 534-545. 10.1096/fj.08-113555.
Angulo Y, Escolano J, Lomonte B, Gutiérrez JM, Sanz L, Calvete JJ: Snake venomics of Central American pitvipers: clues for rationalizing the distinct envenomation profiles of Atropoides nummifer and Atropoides picadoi. J Proteome Res. 2008, 7: 708-719. 10.1021/pr700610z.
Jia Y, Cantu BA, Sánchez EE, Pérez JC: Complementary DNA sequencing and identification of mRNAs from the venomous gland of Agkistrodon piscivorus leucostoma. Toxicon. 2008, 51: 1457-1466. 10.1016/j.toxicon.2008.03.028.
Qinghua L, Xiaowei Z, Wei Y, Chenji L, Yijun H, Pengxin Q, Wingwen S, Songnian H, Guangmei Y: A catalog for transcripts in the venom gland of the Agkistrodon acutus: identification of the toxins potentially involved in coagulopathy. Biochem Biophys Res Commun. 2006, 341: 522-531. 10.1016/j.bbrc.2006.01.006.
Zhang B, Liu Q, Yin W, Zhang X, Huang Y, Luo Y, Qiu P, Su X, Yu J, Hu S, Yan G: Transcriptome analysis of Deinagkistrodon acutus venomous gland focusing on cellular structures and functional aspects using expressed sequence tags. BMC Genomics. 2006, 7: 152-10.1186/1471-2164-7-152.
Casewell NR, Harrison RA, Wüster W, Wagstaff SC: Comparative venom gland transcriptome surveys of the saw-scaled vipers (Viperidae: Echis) reveal substantial intra-family gene diversity and novel venom transcripts. BMC Genomics. 2009, 30: 564-10.1186/1471-2164-10-564.
Sanz L, Gibbs HL, Mackessy SP, Calvete JJ: Venom proteomes of closely related Sistrurus rattlesnakes with divergent diets. J Proteome Res. 2006, 5: 2098-2112. 10.1021/pr0602500.
Guércio AP, Shevchenko A, Shevchenko A, López-Lozano JL, Paba J, Sousa MV, Ricart CAO: Ontogenetic variations in the venom proteome of the Amazonian snake Bothrops atrox. Proteome Sci. 2006, 4: 11-10.1186/1477-5956-4-11.
Zamunér SR, Cruz-Höfling MA, Corrado AP, Hyslop S, Rodrigues-Simioni L: Comparison of the neurotoxic and myotoxic effects of Brazilian Bothrops venoms and their neutralization by commercial antivenom. Toxicon. 2004, 44: 259-271. 10.1016/j.toxicon.2004.05.029.
Pung YF, Wong PT, Kumar PP, Hodgson WE, Kini RM: Ohanin, a novel protein from king cobra venom, induces hypolocomotion and hyperalgesia in mice. J Biol Chem. 2005, 280: 13137-13147. 10.1074/jbc.M414137200.
Pung YF, Kumar SV, Rajagopalan N, Fry BG, Kumar PP, Kini RM: Ohanin, a novel protein from king cobra venom and its cDNA and genomic organization. Gene. 2006, 371: 246-256. 10.1016/j.gene.2005.12.002.
Birrell GW, Early STH, Wallis TP, Masci PP, de Jersey J, Gorman JJ, Lavin MF: The diversity of bioactive proteins in Australian snake venoms. Mol Cell Proteomics. 2007, 6: 973-986. 10.1074/mcp.M600419-MCP200.
Brown AM, Yatani A, Lacerda AE, Gurrola GB, Possani LD: Neurotoxins that act selectively on voltage-dependent cardiac calcium channels. Circ Res. 1987, 61: 16-19.
Possani LD, Martin BM, Yatani A, Mochca-Morales J, Zamudio FZ, Gurrola GB, Brown AM: Isolation and physiological characterization of taicatoxin, a complex toxin with specific effects on calcium channels. Toxicon. 1992, 30: 1343-1364. 10.1016/0041-0101(92)90511-3.
Fantini E, Athias P, Tirosh R, Pinson A: Effect of TaiCatoxin (TCX) on the electrophysiological, mechanical and biochemical characteristics of spontaneously beating ventricular cardiomyocytes. Mol Cell Biochem. 1996, 160-161: 61-66. 10.1007/BF00240032.
Sales PBV, Santoro ML: Nucleotidase and DNase activities in Brazilian snake venom. Comp Biochem Physiol C. 2008, 147: 85-95.
Ogawa Y, Murayama N, Yanoshita R: Molecular cloning and characterization of ecto-5'-nucleotidase from the venoms of Gloydius blomhoffi. Toxicon. 2009, 54: 408-412. 10.1016/j.toxicon.2009.05.004.
Ohno M, Chijiwa T, Oda-Ueda N, Ogawa T, Hattori S: Molecular evolution of myotoxic phospholipases A2 from snake venom. Toxicon. 2003, 42: 841-854. 10.1016/j.toxicon.2003.11.003.
Soto JG, Powell RL, Reyes SR, Wolana L, Swanson LJ, Sanchez EE, Perez JC: Genetic variation of a disintegrin gene found in the American copperhead snake (Agkistrodon contortrix). Gene. 2006, 373: 1-7. 10.1016/j.gene.2005.11.039.
Kordis D, Gubenšek F: Bov-B long interspersed repeated DNA (LINE) sequences are present in Vipera ammodytes phospholipase A2 genes and in genomes of Viperidae snakes. Eur J Biochem. 1997, 246: 772-779. 10.1111/j.1432-1033.1997.00772.x.
Kordis D, Gubenšek F: The Bov-B lines found in Vipera ammodytes toxic PLA2 genes are widespread in snake genomes. Toxicon. 1998, 36: 1585-1590. 10.1016/S0041-0101(98)00150-0.
The authors thank Daniela Stancato for technical assistance and Aline da Costa Lima for help in preparing the figures for the manuscript. K.C.C. is supported by a doctoral studentship from Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP). S.H. is supported by a research fellowship from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq). This work was supported by CNPq.
The authors declare that they have no competing interests.
KCC was responsible for most of the experimental work, sequence annotation, data analysis, and wrote the paper under the supervision of SH. MJDS did part of the experimental work, provided day-to-day supervision of the experimental work and analyzed the results. GGLC performed the bioinformatics analyses, annotated the sequences and helped draft the manuscript. TTT helped with the bioinformatics analyses and manuscript preparation. LEVDB contributed to the sequence annotation. RV helped with the bioinformatics analyses. MM was involved in delineating the experimental design in the initial stages of the work and contributed to discussions of early experiments. SH conceived and wrote the project (together with KCC), supervised the project, helped in analyzing the results, and wrote the paper with KCC. All authors read and commented on the manuscript before submission.
Electronic supplementary material
Additional file 1:Table showing the positions and descriptions of SNPs identified in B. alternatus ESTs using the program QualitySNP. (XLS 48 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Cardoso, K.C., Da Silva, M.J., Costa, G.G. et al. A transcriptomic analysis of gene expression in the venom gland of the snake Bothrops alternatus (urutu). BMC Genomics 11, 605 (2010). https://doi.org/10.1186/1471-2164-11-605