The venom-gland transcriptome of the eastern diamondback rattlesnake (Crotalus adamanteus)
© Rokyta et al.; licensee BioMed Central Ltd. 2012
Received: 14 March 2012
Accepted: 2 July 2012
Published: 16 July 2012
Snake venoms have significant impacts on human populations through the morbidity and mortality associated with snakebites and as sources of drugs, drug leads, and physiological research tools. Genes expressed by venom-gland tissue, including those encoding toxic proteins, have therefore been sequenced but only with relatively sparse coverage resulting from the low-throughput sequencing approaches available. High-throughput approaches based on 454 pyrosequencing have recently been applied to the study of snake venoms to give the most complete characterizations to date of the genes expressed in active venom glands, but such approaches are costly and still provide a far-from-complete characterization of the genes expressed during venom production.
We describe the de novo assembly and analysis of the venom-gland transcriptome of an eastern diamondback rattlesnake (Crotalus adamanteus) based on 95,643,958 pairs of quality-filtered, 100-base-pair Illumina reads. We identified 123 unique, full-length toxin-coding sequences, which cluster into 78 groups with less than 1% nucleotide divergence, and 2,879 unique, full-length nontoxin coding sequences. The toxin sequences accounted for 35.4% of the total reads, and the nontoxin sequences for an additional 27.5%. The most highly expressed toxin was a small myotoxin related to crotamine, which accounted for 5.9% of the total reads. Snake-venom metalloproteinases accounted for the highest percentage of reads mapping to a toxin class (24.4%), followed by C-type lectins (22.2%) and serine proteinases (20.0%). The most diverse toxin classes were the C-type lectins (21 clusters), the snake-venom metalloproteinases (16 clusters), and the serine proteinases (14 clusters). The high-abundance nontoxin transcripts were predominantly those involved in protein folding and translation, consistent with the protein-secretory function of the tissue.
We have provided the most complete characterization of the genes expressed in an active snake venom gland to date, producing insights into snakebite pathology and guidance for snakebite treatment for the largest rattlesnake species and arguably the most dangerous snake native to the United States of America, C. adamanteus. We have more than doubled the number of sequenced toxins for this species and created extensive genomic resources for snakes based entirely on de novo assembly of Illumina sequence data.
Human envenomation by snakes is a worldwide issue that claims more than 100,000 lives per year and exacts untold costs in the form of pain, disfigurement, and loss of limbs or limb function [1–3]. Despite the significance of snakebites, their treatments have remained largely unchanged for decades. The only treatments currently available are traditional antivenoms derived from antisera of animals, usually horses , innoculated with whole venoms [5, 6]; such an approach is the only readily available option for largely uncharacterized, complex mixtures of proteins such as snake venoms. Although often lifesaving and generally effective against systemic effects, these antivenoms have little or no effect on local hemorrhage or necrosis [7–9], which are major aspects of the pathology of viperid bites and can result in lifelong disability [4, 5]. These traditional treatments also sometimes lead to adverse reactions in patients . Advances in treatment approaches will depend on a complete knowledge of the nature of the offending toxins, but current estimates of the numbers of unique toxins present in snake venoms are in excess of 100 , a number not approached in even the most extensive venom-characterization efforts to date .
The significance of snake venoms extends well beyond the selective pressures they may directly impose upon human populations. Snake venoms have evolutionary consequences for those species that snakes prey upon [12, 13], as well as species that prey upon the snakes , and their study can therefore provide insights into predator-prey coevolution. Snake venom components have been leveraged as drugs and drug leads [15–17] and have been used directly as tools for studying physiological processes such as pain reception . In addition to the significance of the toxins, the nature of the extreme specialization of snake venom glands for the rapid but temporary production and export of large quantities of protein could provide insights into basic mechanisms of proteostasis, the breakdown of which is thought to contribute to neurodegenerative diseases such as Parkinson’s and Alzheimer’s .
The eastern diamondback rattlesnake (Crotalus adamanteus) is a pit viper native to the southeastern United States and is the largest member of the genus Crotalus, reaching lengths of up to 2.44 m . The diet of C. adamanteus consists primarily of small mammals (e.g., squirrels, rabbits, and mouse and rat species) and birds, particularly ground-nesting species such as quail . Because of its extreme size and consequent large venom yield, C. adamanteus is arguably the most dangerous snake species in the United States and is one of the major sources of snakebite mortality throughout its range . Crotalus adamanteus has recently become of interest from a conservation standpoint because of its declining range, which at one time included seven states along the southeastern Coastal Plain . This species has now apparently been extirpated from Louisiana and is listed as endangered in North Carolina [23, 24]. As a consequence of recent work by Rokyta et al.  based on 454 pyrosequencing, the venom of C. adamanteus is among the best-characterized snake venoms; 40 toxins have been identified.
Transcriptomic characterizations of venom glands of snakes [25–28] and other animals [29–32] have relied almost exclusively on low-throughput sequencing approaches. Sanger sequencing, with its relatively long, high-quality reads, has been the only method available until recently and has provided invaluable data on the identities of venom genes. Because venomous species are primarily nonmodel organisms, high-throughput sequencing approaches have been slow to pervade the field of venomics (but see Hu et al. ), despite becoming commonplace in other transcriptomic-based fields. Rokyta et al.  recently used 454 pyrosequencing to characterize venom genes for C. adamanteus. More recently, Durban et al.  used 454 sequencing to study the venom-gland transcriptomes of a mix of RNA from eight species of Costa Rican snakes. Whittington et al.  used a hybrid approach with both 454 and Illumina sequencing to characterize the platypus venom-gland transcriptome, although they had a reference genome sequence, making de novo assembly unnecessary. Pyrosequencing is expensive and low-throughput relative to Illumina sequencing, and the high error rate, particularly for homopolymer errors , significantly increases the difficulty of identifying coding sequences without reference sequences.
We sequenced the venom-gland transcriptome of the eastern diamondback rattlesnake with Illumina technology using a paired-end approach coupled with short insert sizes effectively to produce longer, high-quality reads on the order of approximately 150 nt to facilitate de novo assembly (an approach similar to that of Rodrigue et al.  for metagenomics). The difference in read length from that of 454 sequencing was compensated for by the increase of more than two orders of magnitude in the number of reads. We demonstrated de novo assembly and analysis of a venom-gland transcriptome using only Illumina sequences and provided a comprehensive characterization of both the toxin and nontoxin genes expressed in an actively producing snake venom gland.
Results and discussion
Venom-gland transcriptome sequencing and assembly
ABySS assembly summaries
> 200 nt
8.77 × 107
3.70 × 107
3.67 × 107
2.42 × 107
1.79 × 107
1.14 × 107
5.60 × 106
5.52 × 106
3.61 × 105
8.55 × 107
3.52 × 107
3.50 × 107
2.24 × 107
1.62 × 107
1.01 × 107
4.79 × 106
4.72 × 106
3.27 × 105
8.13 × 107
3.28 × 107
3.27 × 107
2.03 × 107
1.43 × 107
8.70 × 106
4.02 × 106
3.95 × 106
2.93 × 105
7.54 × 107
3.00 × 107
3.00 × 107
1.79 × 107
1.24 × 107
7.28 × 106
3.28 × 106
3.23 × 106
2.47 × 105
6.79 × 107
2.64 × 107
2.65 × 107
1.51 × 107
1.02 × 107
5.84 × 106
2.52 × 106
2.54 × 106
2.12 × 105
NGen assembly summaries
Total unique full-length toxins =
Expression levels of full-length toxin clusters
% total reads
% toxin reads
GenBank TSA accessions
9.68 × 10−1
9.38 × 10−1
9.15 × 10−1
9.14 × 10−1
9.13 × 10−1
8.97 × 10−1
8.38 × 10−1
8.30 × 10−1
7.82 × 10−1
7.74 × 10−1
7.69 × 10−1
6.98 × 10−1
6.73 × 10−1
6.31 × 10−1
6.00 × 10−1
5.96 × 10−1
5.41 × 10−1
5.41 × 10−1
5.40 × 10−1
5.33 × 10−1
5.15 × 10−1
5.04 × 10−1
4.80 × 10−1
4.66 × 10−1
4.31 × 10−1
4.24 × 10−1
4.15 × 10−1
3.59 × 10−1
3.59 × 10−1
3.58 × 10−1
3.49 × 10−1
3.36 × 10−1
3.15 × 10−1
2.90 × 10−1
2.88 × 10−1
2.57 × 10−1
2.41 × 10−1
2.23 × 10−1
2.16 × 10−1
1.90 × 10−1
1.71 × 10−1
1.71 × 10−1
1.59 × 10−1
1.34 × 10−1
1.19 × 10−1
1.16 × 10−1
1.12 × 10−1
1.10 × 10−1
1.09 × 10−1
1.05 × 10−1
9.53 × 10−2
9.40 × 10−2
7.35 × 10−2
5.83 × 10−2
5.62 × 10−2
4.86 × 10−2
3.22 × 10−2
3.09 × 10−2
2.22 × 10−2
2.10 × 10−2
1.14 × 10−2
7.83 × 10−3
5.11 × 10−3
4.99 × 10−3
4.46 × 10−3
3.65 × 10−3
3.36 × 10−3
3.14 × 10−3
2.68 × 10−3
2.24 × 10−3
1.70 × 10−3
5.60 × 10−4
3.80 × 10−4
We used the number or percentage of reads mapping to a particular transcript as a measure of its abundance. Although average coverage might be a more appropriate proxy for the number of copies of a given transcript present, because it accounts for differences in transcript lengths, we prefer read counts as a measure of the expression expenditure on a given transcript because they better reflect the energetic cost associated with producing the encoded protein and are consistent with previous work using low-throughput sequencing (see, e.g., Pahari et al. ). In addition, this measurement should more closely match proteomic-based measurements of the contents of venom components (see, e.g., Gibbs et al. ) which come in the form of the percentages of total peptide bonds in the sample.
Snake venom metalloproteinases
We identified 39 unique sequences and 16 clusters of snake-venom metalloproteinases (SVMPs) that accounted for 24.4% of the reads mapping to toxin sequences and 8.6% of the total reads (Figure 3A and Table 3). In terms of total reads, the SVMPs were the most abundant class of toxins in the C. adamanteus venom-gland transcriptome. SVMPs are the primary sources of the local and systemic hemorrhage associated with envenomation by viperids and are divided into a number of subclasses based on their domain structure [44, 45]. All SVMPs have a metalloproteinase domain characterized by a zinc-binding motif. All of the SVMPs identified for C. adamanteus belong to either the type II or the type III subclass. Type II SVMPs (SVMPIIs) have a disintegrin domain in addition to the metalloproteinase domain, which may be proteolytically cleaved posttranslationally to produce a free disintegrin. Type III SVMPs (SVMPIIIs) have a disintegrin-like and a cysteine-rich domain in addition to the metalloproteinase domain. We found 8 clusters of each of these two subclasses with 23 unique SVMPII sequences and 16 unique SVMPIII sequences. SVMPII and SVMPIII clusters comprise 16.4% and 8.0% of the reads mapping to toxins respectively (Figure 3). The sequences in both subclasses are diverse. The maximum pairwise nt divergence for the SVMPIIs was 10.0%, corresponding to a maximum amino-acid divergence of 18.1%. For the SVMPIIIs, the maximum pairwise nt divergence was 20.4% with a maximum amino-acid divergence of 42.3%. Although SVMPs were the dominant toxins as a class, the individual SVMP cluster with the highest abundance was SVMPII-5, which was only the eighth most abundant toxin cluster (Figure 2B and Table 3).
Mackessy  categorized rattlesnake venoms as type I or type II on the basis of their toxicities and metalloproteinase activities. These two measurements tend to be inversely related in rattlesnakes: species (or populations) with low LD50 values tend also to have low or undetectable hemorrhagic activities. SVMPs are the major hemorrhagic components of snake venoms, and high toxicity appears to be caused mostly by neurotoxic venom components. Low-toxicity venoms with high metalloproteinase activity are classified as type I, and high-toxicity venoms with low metalloproteinase activity are classified as type II. On the basis of the abundance of SVMPs in the venom-gland transcriptome, C. adamanteus clearly has type I venom, although the relatively low toxicity of its venom  is at least partially compensated for by its large size and venom yield.
The most diverse and the second most abundant toxin class in the C. adamanteus venom-gland transcriptome was the C-type lectin (CTL) class. We identified 37 unique sequences and 21 clusters of CTLs that accounted for 22.2% of the reads mapping to toxins and 7.8% of the total reads (Figure 3A and Table 3). CTLs generally either inhibit or activate components of plasma or blood-cell types, thereby interfering with hemostasis . Most known snake-venom CTLs function as heterodimers or even more complex arrangements , probably accounting in part for their diversity. The divergence among members of this class within the C. adamanteus genome was extreme, although all members preserved a CTL-like domain. Some pairs shared virtually no conserved amino-acid positions. Three of the CTL clusters provide evidence for the relevance of alternative splicing in the generation of toxin proteins. CTL-3f, CTL-4e, and CTL9b all have 48-nt insertions in the same region but are otherwise similar or identical to other members of their clusters.
Snake venom serine proteinases
The third most abundant toxin class for C. adamanteus was the snake-venom serine proteinases (SVSPs). We identified 18 unique sequences and 14 clusters in this toxin class, accounting for 20.0% of the toxin reads and 7.1% of the total reads (Figure 3A and Table 3). Three of the 10 most highly expressed individual toxins were SVSPs (Figure 2). SVSPs interfere with a wide array of reactions involving blood coagulation and hemostasis and belong to the trypsin family of serine proteases [49, 50]. Mackessy  detected significant thrombin-like and kallikrein-like activity in the venom of C. adamanteus, which are attributable to the action of SVSPs. The diversity of SVSPs within the C. adamanteus genome is high; maximum pairwise nt divergence is 20.6% and amino-acid divergence is 47.4%.
The members of two SVSP clusters differ in a way that should be noted. The lengths of SVSPs are generally well conserved throughout the class. SVSP-7a has a 27-nt insertion relative to the two other members of its cluster but is otherwise identical to SVSP-7b. This difference could reflect the presence of alternative splicing for this gene. SVSP-3a is unique among the C. adamanteus SVSPs or those known from other snake species in apparently having a 65-amino-acid extension of its C-terminal region. The other member of its cluster, SVSP-3b, has a single deletion of a C nt in a poly-C tract that terminates its coding sequence consistently with other known SVSPs. The reads generating the SVSP-3a form vastly outnumber those for the SVSP-3b form; more than 95% of the reads support the extended version of the protein. The effect, if any, of this C-terminal extension remains to be determined.
Previous work with C. adamanteus identified only a single phospholipase A2 (PLA2) sequence , but we identifed seven unique sequences in six clusters (Figure 2 and Table 3), accounting for 7.8% of the toxin reads and 2.8% of the total reads (Figure 3). PLA2s are among the most functionally diverse classes of snake-venom toxins and have pharmocological effects ranging from neurotoxicity (presynaptic or postsynaptic) to myotoxicity and cardiotoxicity. Anticoagulant and hemolytic effects due to PLA2s are also known [51, 52]. Compared to other toxin classes of C. adamanteus, the diversity of PLA2s is low. Five of the six clusters are all within 5% nt divergence of one another. PLA2-3 is the lone, high-divergence outlier, differing by more than 31% at the nt level from the other clusters. PLA2-3 is also expressed at the lowest level of any of the PLA2s (Table 3).
Other high-abundance toxins
The SVMPs, CTLs, SVSPs, and PLA2s account for 74% of the reads mapping to toxin sequences (Figure 3), 73% of the toxin clusters, and 82% of the unique toxin sequences. The remaining toxins belong to 16 different classes. Many of these are low-abundance transcripts (Figure 2 and Table 3) and may not actually function as significant toxins, whereas several others have high to moderate abundances and represent significant components of the venom.
The most abundant toxin transcript and the most abundant transcript overall (Figure 2) was a small basic myotoxin related to crotamine [53, 54]. The precursor protein is just 70 amino acids in length with a predicted 22-amino-acid signal peptide. This transcript was detected by Rokyta et al. , but the coding sequence was prematurely truncated in their sequence because of a single nt deletion. This toxin accounts for 16.8% of the toxin reads (Figure 3A) and 5.9% of the total reads. Crotamine, originally isolated from the venom of C. durissus, causes spastic paralysis in mice and is found in the venoms of many species of Crotalus. Muscle spasms, twitching, and paralysis of the legs have been reported for human envenomations by C. adamanteus. Interestingly, Straight et al.  noted that individuals of C. adamanteus from populations in southern and central Florida lack this toxin in their venoms. Given that this myotoxin is the most abundant transcript in the venom of our specimen, its absence in southern populations points to a dramatic difference in venoms within this species and the potential for significantly different pathological effects associated with bites from different C. adamanteus populations.
A single L-amino-acid oxidase (LAAO) transcript was the second most abundant toxin transcript (Figure 2B), consistent with the previously detected LAAO activity in the C. adamanteus venom . This single transcript accounted for 5.3% of the reads mapping to toxins and 1.9% of the total reads. LAAOs are flavoproteins, giving the venom its yellow color; can be edema- or apoptosis-inducing; and can induce or inhibit platelet aggregation . These effects are probably mediated by H2O2released during the oxidation reaction catalyzed by the enzyme. The 29th most abundant toxin transcript was a cysteine-rich secretory protein (CRISP) (Figure 2B and Table 3), accounting for 1.3% of the toxin reads (Figure 3A). Although CRISPs are widely found in snake venoms, their precise effects are not well established , but they appear to interfere with smooth-muscle contraction [58, 59]. A single transcript for a bradykinin-potentiating and C-type natriuretic peptide transcript (BPP) was found to account for 0.7% of the toxin reads (Figure 3A). The encoded protein is similar to a protein identified in Sistrurus catenatus (GenBank accession: DQ464265) that was hypothesized to reduce blood pressure in envenomated prey . A loss of blood pressure has been reported in human envenomations by C. adamanteus.
Other low-abundance toxins
The remaining 17 clusters are classified as “others” in Figure 3A. Because each has a relatively low expression level (Table 3), many of these should be considered putative toxins until their presence in the C. adamanteus venom is confirmed proteomically and pharmacological effects are associated with them.
Rokyta et al.  detected the presence of a transcript encoding a protein homologous to ohanin from Ophiophagus hannah[60, 61] and to a homologous protein from Lachesis muta; we found a transcript identical to that of Rokyta et al. . Pung et al. [60, 61] found the O. hannah version of this protein to increase pain sensitivity (hyperalgesia) and to induce temporary hypolocomotion in mice and proposed naming the class vespryns (VESP). Exceptionally intense pain has been reported after envenomation of humans by C. adamanteus, although whether such pain is due to a specific toxin is not clear.
We detected three different nucleotidases (NUCs) and five different phosphodiesterases (PDEs) in the venom-gland transcriptome of C. adamanteus. Only one of the NUCs and three of the PDEs had signal peptides, and we therefore only considered these as potential toxins: NUC, PDE, PDE-4, and PDE-6 (Table 3). The roles of these enzymes in venoms are uncertain, but their primary function may be to liberate toxic nucleosides [63–65]. Significant PDE activity has been detected previously in the venom of C. adamanteus.
The C. adamanteus venom-gland transcriptome contained three Kunitz-type protease inhibitors (KUNs). Two of these shared more than 75% animo-acid identity with a KUN from Austrelaps labialis (GenBank accession: B2BS84), an Australian elapid. All three KUNs have domains that place them in the superfamily of bovine pancreatic trypsin-like inhibitors, and snake toxins from this family are known to inhibit plasma serine proteinases. Although KUNs are commonly observed in snake venoms, their role in envenomation (if any) is not well defined . The three KUNs detected for C. adamanteus are all at relatively low abundances, suggesting that they are not major components of the venom.
We identified two transcripts, HYAL-1 and HYAL-2, encoding hyaluronidase-like proteins. Hyaluronidases are generally regarded as venom components that promote the dissemination of other venom components by degrading the extracellular matrix at the site of injection , although they may have more direct toxic effects . The coding sequences of our two transcripts differ only in the presence of a 765-nt deletion in HYAL-2 relative to HYAL-1. Truncated hyaluronidases such as HYAL-2 have been detected in the venoms of other viperid species  and may represent an example of alternative splicing. We also identified a transcript encoding a glutaminyl-peptide cyclotransferase (glutaminyl cyclase; GC). Many snake venom components have N termini blocked by pyroglutamate, and GCs catalyze the formation of this block. This component is related more to maturation and protection of other toxins and probably contributes only indirectly to toxicity .
We identified six growth-factor-related sequences in the venom-gland transcriptome of C. adamanteus: a nerve growth factor (NGF), a neurotrophic factor (NF), two vascular endothelial growth factors (VEGF) in a single cluster, and a cysteine-rich with EGF-like domain protein (CREGF). The NGF transcript encodes a 241 amino-acid precursor protein and shares 99% amino-acid identity with a NGF from C. durissus (GenBank accession: AAG30924). The NF transcript encodes a 180-amino-acid precursor that shares homology with mesencephalic astrocyte-derived neurotrophic factors. We found no close venom-related sequences for this NF in the available databases. The VEGF sequences appear to be alternatively spliced versions of one another. VEGF-1a encodes a 192-amino-acid precursor, and VEGF-1b encodes a 148-amino-acid precursor. Aside from the 132-nt deletion in VEGF-1b relative to VEGF-1a, their coding sequences are identical. Both forms have database matches of the same length with 99% amino-acid identity from Trimeresurus flavoviridis (GenBank accessions: AB154418 and AB154419). Finally, we detected the same cysteine-rich with EGF-like domain protein as described by Rokyta et al. .
The final two putative toxin transcripts are of questionable significance because of their low expression levels. A single sequence with 77% amino-acid identity to a waprin (WAP) sequence from Philodryas olfersii (GenBank accession: EU029742), a rear-fanged colubrid, was detected. Related sequences have been detected in a variety of other rear-fanged snake species, but such proteins are only known to exhibit antimicrobial activity . We detected a venom factor (VF) transcript that shares 87% animo-acid identity with a VF from Austrelaps superbus (GenBank accession: AY903291) . The C. adamanteus VF transcript encodes a 1,652-amino-acid precursor with a 22-amino-acid signal peptide. The best-studied member of this toxin family is cobra venom factor, which is known to activate the complement system . The extremely low expression levels of these transcripts may indicate that they represent the orthologous genes to the ancestors of the known toxic forms and may therefore have no toxic functions.
Comparison to previous work
Correspondence with the results of Rokyta et al.]
% nt divergence
454 version incomplete
454 version has 1-nt deletion that truncates the coding sequence prematurely
454 version incomplete
454 version has 123-nt insertion
454 version incomplete; no signal peptide; no longer considered toxin
No longer considered toxin
454 version incomplete
454 version incomplete
454 version incomplete
454 version incomplete
454 version has 1-nt deletion that truncates the coding sequence prematurely
454 version incomplete
454 version incomplete
454 version incomplete
The 20 most highly expressed nontoxin transcripts
Protein disulfide isomerase
Rearrange disulfide bonds (ER)
Cytochrome C oxidase subunit I
Electron transport chain
Electron transport chain
Translation elongation factor 1 α 1
Protein chaperone (ER)
Endoplasmin (HSP90 family)
Protein chaperone (ER)
78 kDa glucose-regulated protein
Protein chaperone (ER)
Heat shock protein 5 (GRP78 splice variant?)
Protein chaperone (ER)
NADH dehydrogenase subunit 5
Electron transport chain
Cytochrome C oxidase subunit III
Electron transport chain
Protein disulfide isomerase A6
Rearrange disulfide bonds (ER)
Protein disulfide isomerase A3
Rearrange disulfide bonds (ER)
NADH dehydrogenase subunit 1
Electron transport chain
NADH dehydrogenase subunit 4
Electron transport chain
Protein disulfide isomerase A4
Rearrange disulfide bonds (ER)
Translation elongation factor 2
mRNA stability and translation
Actin, cytoplasmic 2
The abundances of several major classes of nontoxins are provided in Figure 3B. We identified 57 sequences with functions related to protein folding [19, 78–80], including various classes of heat-shock proteins, protein-disulfide isomerases, peptidyl-prolyl cis-trans isomerases, dnaJ-complex components, and T-complex components. These sequences together accounted for 28.4% of the total reads mapping to nontoxins. Ribosomal-protein transcripts (cytoplasmic and mitochondrial) accounted for 9.5% of the nontoxin reads, and mitochondrial genes accounted for another 9.0%. Finally, we identified 110 sequences transcripts encoding proteins involved in protein degradation [81, 82], including proteins involved in the ubiquitin-proteasome system and the ER-associated protein-degradation system , which accounted for 2.6% of the nontoxin reads. Protein-quality control should be essential in a high-throughput protein-producing tissue such as a snake venom gland.
Toxin and protease inhibitors detected in the venom-gland transcripts
9.80 × 10−4
1.23 × 10−3
1.31 × 10−3
Metalloproteinase inhibitor 1
1.04 × 10−3
Metalloproteinase inhibitor 2
2.57 × 10−3
Metalloproteinase inhibitor 3
1.01 × 10−3
PLA2 inhibitor beta
1.54 × 10−3
PLA2 inhibitor gamma B 1
3.68 × 10−3
PLA2 inhibitor gamma B 2
7.70 × 10−4
PLA2 inhibitor B
1.07 × 10−3
9.68 × 10−3
9.40 × 10−4
Sequence accession numbers
The original, unmerged sequencing reads were submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive under accession number SRA050594. The annotated toxin and nontoxin sequences were submitted to the GenBank Transcriptome Shotgun Assembly (TSA) database under accession numbers JU173621–JU173743 (toxins) and JU173744–JU176622 (nontoxins).
We have described the most comprehensive venom-gland transcriptomic characterization of a snake species to date and provided full-length coding sequences for 123 unique toxin proteins and 2,879 unique nontoxin proteins. We have demonstrated the use of Illumina sequencing technology for the sequencing and de novo assembly of a tissue-specific transcriptome for a nonmodel species, C. adamanteus, for which genome-scale resources were previously unavailable. Because the nontoxin sequences in particular should be conserved across snake species, our results should greatly facilitate similar work with other venomous species, serving as an assembly template and reducing the number of reads for which de novo assembly will be necessary.
The expressed toxin genes in the venom gland of C. adamanteus provide a detailed portrait of a type I rattlesnake venom . The most abundant transcript expressed in the C. adamanteus venom gland encoded a myotoxin homologous to crotamine. Crotamine is known to induce spastic paralysis , a symptom that has been observed in human envenomations by C. adamanteus. Like those of most viperids, the bites of C. adamanteus result in significant tissue damage and necrosis, and we found that SVMPs, the major class of hemorrhagic toxins, dominated venom-gland gene expression. The second most abundant toxin transcript overall was an LAAO, which are also noted for causing local tissue damage . Coagulopathy is a common occurrence with pit-viper bites . The CTLs and SVSPs were also both diverse and abundant in the venom-gland transcriptome of C. adamanteus, and both classes primarily attack the hemostatic system. In terms of gene sequences of venom components, the venom of C. adamanteus is now the best-characterized snake venom, although a thorough proteomic analysis of the venom is still needed. The sequences we have generated will greatly facilitate such a proteomic characterization by serving as a database against which to query mass-spectrum results.
The expression patterns of the nontoxin genes in the venom gland of C. adamanteus reflect the protein-secretory function of the tissue and the high energetic demands of rapid venom production . The most highly expressed nontoxin genes were those involved in the production and processing of proteins and energy production to support these activities. Molecular chaperones and PDIs were particularly abundant. Though the expression patterns for nontoxins were not surprising, future comparisons with other snake species, especially those from other snake families, may be able to elucidate the origin and early stages of the evolution of the venom gland.
Venom-gland transcriptome sequencing
We sequenced the venom-gland transcriptome of a single animal from Florida (Wakulla County): an adult female weighing 393 g with a snout-to-vent length of 792 mm and a total length of 844 mm. To stimulate transcription in the venom glands, we anesthetized the snake by propofol injection (10 mg/kg) and extracted venom by electrostimulation under anesthesia . After venom extraction, the animal was allowed to recover for four days while transcription levels reached their maxima . The snake was euthanized by injection of sodium pentobarbitol (100 mg/kg), and its venom glands were subsequently removed. The above techniques were approved by the Florida State University Institutional Animal Care and Use Committee (IACUC) under protocol #0924.
Sequencing and nonnormalized cDNA library preparation were performed by the HudsonAlpha Institute for Biotechnology Genomic Services Laboratory (http://www.hudsonalpha.org/gsl/). Transcriptome sequencing was performed essentially as described by Mortazavi et al.  in a modification of the standard Illumina methods described in detail in Bentley et al. . Total RNA was reduced to poly-A+ RNA with oligo-dT beads. Two rounds of poly-A+ selection were performed. The purified mRNA was then subjected to a mild heat fragmentation followed by random priming for first-strand synthesis. Standard second-strand synthesis was followed by standard library preparation with the double-stranded cDNA as input material. This approach is similar to that of Illumina’s TruSeq RNA-seq library preparation kit. Sequencing was performed in one lane on the Illumina HiSeq 2000 with 100-base-pair paired-end reads.
Transcriptome assembly and analysis
assuming any of the four nucleotides is equally likely to be at any position. To be conservative, we only merged reads if the minimum probability was less than 10−10and the second smallest probability was at least 1000 times larger (Figure 1A). The latter condition was meant to help avoid merging reads that span highly repetitive regions. For cases in which the insert size was less than the read length, sequence data outside the overlap were assumed to represent adaptors and were deleted. We updated quality scores for the overlapping positions following the approach of Rodrigue et al. . For merged reads, quality scores for nonoverlapping bases were left unchanged (Figure 1B). The unmerged reads were typically those pairs from the longer end of the insert-size distribution.
Because of the inherent difficulty in de novo transcriptome assembly, we used a diverse array of assembly approaches and combined the results for a final data set. We performed assemblies using ABySS version 1.2.6 [37, 38] under a wide array of parameter values using both the merged and unmerged reads. In particular, we used k-mer values of 51, 61, 71, 81, and 91 and varied the coverage (c) and erode (e) parameters from 2 to 1,000. We set E = 0, m = 20, and s = 200 for all assemblies. Trans-ABySS  provided little or no improvement of our assemblies, primarily because assembly quality appeared to be more dependent on the coverage and erode parameters than on the k-mer length. We also conducted assemblies using both the merged and unmerged reads with Velvet version 1.1.02  and k-mer values of 71, 81, and 91. We selected the best of these assemblies on the basis of the N 50 values for further assembly into transcripts with Oases version 0.1.20 (http://www.ebi.ac.uk/∼zerbino/oases/) . For Oases, we set the minimum transcript length to 300 nt and the coverage cutoff to 10. We also followed the approach of Rokyta et al.  and used the NGen2.2 assembler from DNAStar (http://www.dnastar.com/). Because this assembler is limited to 20–30 million reads, we used only the merged reads. We performed four independent assemblies: three with 20 million merged reads each and one with the remaining 12,114,709 merged reads. Each assembly was performed with the default settings for high-stringency, de novo transcriptome assembly for long Illumina reads, including default quality trimming. The high-stringency setting corresponded to setting the minimum match percentage to 90%. We retained contigs comprising at least 100 reads.
In addition to the all-at-once assembly approaches above, we developed an iterative approach that was both more effective at generating full-length transcripts and more computationally efficient. The first step consisted of applying our Extender program (see below) as a de novo assembler starting from 1,000 reads. Full-length transcripts were identified with blastx searches (see below), then used as templates in a reference-based assembly in NGen3.1 with a 98% minimum match percentage to filter reads corresponding to identified transcripts. Ten million of the unassembled sequences were then used in a de novo transcriptome assembly in NGen3.1 with the same settings as described above for de novo assembly except that the minimum match percentage was increased to 93% and contigs comprising less than 200 sequences were discarded. The resulting sequences were identified, where possible, by means of blastx searches, and the identified full-length transcripts were used in another templated assembly to generate a further-reduced set of reads. This iterative process was repeated two additional times.
To provide transcriptional profiles of the venom gland, we performed GO annotation with Blast2GO . We ran full analyses on one of NGen assemblies of 20 million merged reads, including blastx searches, GO mapping, and annotation. We used the default Blast2GO parameters throughout. We converted the GO annotation to generic GO-slim terms. We ran the same analysis on the combined set of annotated nontoxin sequences.
For gene identification and annotation, we conducted blastx searches using mpiblast version 1.6.0 (http://www.mpiblast.org/) of the consensus sequences of contigs of our assemblies against the NCBI nonredundant protein database (nr; downloaded March 2011 and updated through November 2011). We used an E-value cut-off of 10−4, and only the top 10 matches were considered. For toxin identification, hit descriptions were searched for a set of keywords based on known snake-venom toxins and protein classes. Any sequence matching these keywords was checked for a full-length coding sequence. We generally only retained transcripts with full-length coding sequences (but see below). For the iterative assembly approach, the remaining, presumably nontoxin-encoding, contigs were screened for those whose match lengths were at least 90% of the length of at least one of their database matches. This step was intended to minimize the number of fragmented or partial sequences that were considered for annotation. In addition, we sorted the contigs of the three 20-million-sequence NGen assemblies from the all-at-once approach on the basis of the number of reads and attempted to annotate the top 500 contigs from one assembly and the top 100 from the other two.
We estimated transcript abundances using high-stringency reference-based assemblies in NGen3.1 with a minimum match percentage of 95. Ten million of the merged reads were mapped onto the full-length, annotated transcripts, and the percentage of reads mapping to each transcript was used as a proxy for abundance.
The purpose of Extender is to estimate quickly one or more full-length transcript sequences from a large number of high-quality sequence reads. The procedure begins with one or more seed sequences provided by the user. The seeds can be known sequences (e.g., partial transcripts from a previous assembly) or simply sequences of one or more of the reads. The Extender procedure begins by hashing the k-mers observed at the two ends of the seeds. If k is set to 50, for example, then the 50-base sequence present at the 5’ end of each seed is used as a key in a hash table, and the hash value is a pointer to the seed in the list of seeds. A second hash table is likewise used for k-mers from the 3’ ends of the seeds. Note that this method requires that all initial k-mers be unique (that no two sequence ends be identical). Once the seeds are hashed, the seeds are extended with the set of reads provided by the user as follows. The two k-mers from the ends of each read are looked up in each hash table. If the key is present in the hash table, the seed is extended by concatenation of the nonoverlapping bases from the read onto the appropriate end of the seed. If the key is absent, the reverse complement of the read is used to extend the seed if the end k-mers are found. After each extension, the k-mer key facilitating the extension is removed from the hash table and the new k-mer key is added (the reference to the seed remains the same). The procedure is repeated until the reads have been cycled through N times, where N is chosen by the user. Cycling is beneficial because the Extender does not reset to the beginning of the read list when an extension is made.
Extension of a seed typically terminates when the end of the full-length transcript is reached or when a sequencing error is encountered in the end of an incorporated read. The presence of low-frequency biological artifacts (e.g., unspliced introns) may also result in termination of the extension. In order to improve the accuracy of the consensus sequence prediction, Extender can create replicate seeds for a particular seed by sequentially trimming one base at a time from both ends. Using replicate seeds allows several independent sequences that represent the same target consensus sequence to be generated simultaneously, and these replicates are entirely independent because they begin with different keys. The user can obtain the final estimate of the sequence corresponding to each original seed by taking the consensus across replicates or by simply choosing the replicate producing the longest sequence. We took the former approach for all of our assembly efforts. Overall, Extender is highly inefficient with its use of data and requires many long, high-quality reads, but it is extremely computationally efficient, having short run times and low memory requirements.
We used Extender in two different ways: to complete partial toxin transcripts and as a de novo assembler. For the former, we used partial toxin transcripts from NGen assemblies that were found to have fragments of coding sequence homologous with known toxins. The partial transcripts were trimmed to just the partial coding sequence and used as seeds. To use Extender as a de novo assembler, we seeded it with 1,000 random reads. For both applications, we used a k-mer size of 100, 20 replicates, 10 cycles through the complete set of merged reads excluding all reads with any bases with quality scores less than 30.
Bradykinin potentiating and C-type natriuretic peptides
Cysteine-rich with EGF-like domain
Cysteine-rich secretory protein
Kunitz-type protease inhibitor
L amino-acid oxidase
Nerve growth factor
Protein disulfide isomerase
Snake venom metalloproteinase (types II and III)
Snake venom serine proteinase
Vascular endothelial growth factor
The authors thank Kenneth P. Wray for dissecting the venom glands and Darryl Heard for training DRR in the electrostimulation technique for venom extraction. Computational resources were provided by the Florida State University High-Performance Computing cluster, and the authors thank James C. Wilgenbusch for assistance in the use of these resources. Funding for this work was provided to DRR and ARL by Florida State University.
- Chippaux JP: Snake-bites: appraisal of the global situation. Bull WHO. 1998, 76: 515-524.PubMed CentralPubMed
- O’Neil ME, Mack KA, Gilchrist J, Wozniak EJ: Snakebite injuries treated in United States emergency departments, 2001–2004. Wilderness Env Med. 2007, 18: 281-287. 10.1580/06-WEME-OR-080R1.1.View Article
- Langley RL: Deaths from reptile bites in the United States, 1979–2004. Clin Toxicol. 2009, 47: 44-47. 10.1080/15563650801968313.View Article
- Theakston RDG, Warrell DA, Griffiths E: Report of a WHO workshop on the standardization and control of antivenoms. Toxicon. 2003, 41: 541-557. 10.1016/S0041-0101(02)00393-8.View ArticlePubMed
- Smith J, Bush S: Envenomations by reptiles in the United States. Handbook of Venoms and Toxins of Reptiles. Edited by: Mackessy SP. 2010, CRC Press, Boca Raton, Florida, 475-490.
- Neves-Ferreira AGC, Valente RH, Perales J, Domont GB: Natural inhibitors: innate immunity to snake venoms. Handbook of Venoms and Toxins of Reptiles. Edited by: Mackessy SP. 2010, CRC Press, Boca Raton, Florida, 259-284.
- Rucavado A, Lomonte B: Neutralization of myonecrosis, hemorrhage, and edema induced by Bothrops apser snake venom by homologous and heterologous pre-existing antibodies in mice. Toxicon. 1996, 34: 567-577. 10.1016/0041-0101(95)00162-X.View ArticlePubMed
- Huang KF, Chow LP, Chiou SH: Isolation and characterization of a novel proteinase inhibitor from the snake serum of Taiwan habu (Trimeresurus mucrosquamatus). Biochem Biophys Res Commun. 1999, 263: 610-616. 10.1006/bbrc.1999.1421.View ArticlePubMed
- Valente RH, Dragulev B, Perales J, Fox JW, Domont GB: BJ46a, a snake venom metalloproteinase inhibitor. Eur J Biochem. 2001, 268: 3042-3052. 10.1046/j.1432-1327.2001.02199.x.View ArticlePubMed
- Serrano SMT, Shannon JD, Wang D, Camargo ACM, Fox JW: A multifaceted analysis of viperid snake venoms by two-dimensional gel electrophoresis: An approach to understanding venom proteomics. Proteomics. 2005, 5: 501-510. 10.1002/pmic.200400931.View ArticlePubMed
- Rokyta DR, Wray KP, Lemmon AR, Lemmon EM, Caudle SB: A high-throughput venom-gland transcriptome for the eastern diamondback rattlesnake (Crotalus adamanteus) and evidence for pervasive positive selection across toxin classes. Toxicon. 2011, 57: 657-671. 10.1016/j.toxicon.2011.01.008.View ArticlePubMed
- Biardi JE, Chien DC, Coss RG: California ground squirrel (Spermophilus beecheyi) defenses against rattlesnake venom digestive and hemostatic toxins. J Chem Ecol. 2005, 31: 2501-2518. 10.1007/s10886-005-7610-1.View ArticlePubMed
- Biardi JE, Nguyen KT, Lander S, Whitley M, Nambiar KP: A rapid and sensitive fluorometric method for the quantitative analysis of snake venom metalloproteases and their inhibitors. Toxicon. 2011, 57: 342-347. 10.1016/j.toxicon.2010.12.014.PubMed CentralView ArticlePubMed
- Jansa SA, Voss RS: Adaptive evolution of the venom-targeted vWF protein in opossums that eat pitvipers. PLoS One. 2011, 6: e20997-10.1371/journal.pone.0020997.PubMed CentralView ArticlePubMed
- Harvey AL, Bradley KN, Cochran SA, Rowan EG, Pratt JA, Quillfeldt JA, Jerusalinsky DA: What can toxins tell us for drug discovery?. Toxicon. 1998, 36: 1635-1640. 10.1016/S0041-0101(98)00156-1.View ArticlePubMed
- Ménez A: Functional architectures of animal toxins: a clue to drug design?. Toxicon. 1998, 36: 1557-1572. 10.1016/S0041-0101(98)00148-2.View ArticlePubMed
- Escoubas P, King GF: Venomics as a drug discovery platform. Expert Rev Proteomics. 2009, 6: 221-224. 10.1586/epr.09.45.View ArticlePubMed
- Bohlen CJ, Chesler AT, Sharif-Naeini R, Medzihradszky KF, Zhou S, King D, Sánchez EE, Burlingame AL, Basbaum AI, Julius D: A heteromeric Texas coral snake toxin targets acid-sensing ion channels to produce pain. Nature. 2011, 479: 410-414. 10.1038/nature10607.PubMed CentralView ArticlePubMed
- Hartl FU, Bracher A, Hayer-Hartl M: Molecular chaperones in protein folding and proteostasis. Nature. 2011, 475: 324-332. 10.1038/nature10317.View ArticlePubMed
- Klauber LM: Rattlesnakes: Their Habits, Life Histories, and Influence on Mankind. 1997, University of California Press, Berkeley, California
- Gold BS, Dart RC, Barish RA: Bites of venomous snakes. N Engl J Med. 2002, 347: 347-356. 10.1056/NEJMra013477.View ArticlePubMed
- Conant R, Collins JT: A Fieldguide to Reptiles and Amphibians of Eastern and Central North America. 1998, Houghton Mifflin Harcourt, New York, New York
- Palmer WM, Braswell AL: Reptiles of North Carolina. 1995, University of North Carolina Press, Chapel Hill, North Carolina
- Dundee HA, Rossman DA: The Amphibians and Reptiles of Louisiana. 1996, Louisiana University Press, Baton Rouge, Louisiana
- Pahari S, Mackessy SP, Kini RM: The venom gland transcriptome of the desert 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.PubMed CentralView ArticlePubMed
- 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, 10: 564-10.1186/1471-2164-10-564.PubMed CentralView ArticlePubMed
- Leão LI, Ho PL, de L M Junqueira-de Azevedo I: Transcriptomic basis for an antiserum against Micrurus corallinus (coral snake) venom. BMC Genomics. 2009, 10: 112-10.1186/1471-2164-10-112.PubMed CentralView ArticlePubMed
- Jiang Y, Li Y, Lee W, Xu X, Zhang Y, Zhao R, Zhang Y, Wang W: Venom gland transcriptomes of two elapid snakes (Bungarus multicinctus and Naja atra) and evolution of toxin genes. BMC Genomics. 2011, 12: 1-10.1186/1471-2164-12-1.PubMed CentralView ArticlePubMed
- Morgenstern D, Rohde BH, King GF, Tal T, Sher D, Zlotkin E: The tale of a resting venom gland: transcriptome of a replete venom gland from the scorpion Hottentotta judaicus. Toxicon. 2011, 57: 695-703. 10.1016/j.toxicon.2011.02.001.View ArticlePubMed
- Whittington CM, Papenfuss AT, Locke DP, Mardis ER, Wilson RK, Abubucker S, Mitreva M, Wong ESW, Hsu AL, Kuchel PW, Belov K, Warren WC: Novel venom gene discovery in the platypus. Genome Biol. 2010, 11: R95-10.1186/gb-2010-11-9-r95.PubMed CentralView ArticlePubMed
- Gremski LH, Silveira RBD, Chaim OM, Probst CM, Ferrer VP, Nowatzki J, Weinschutz HC, Madeira HM, Gremski W, Nader HB, Senff-Ribeiro A, Veiga SS: A novel expression profile of the Loxosceles intermedia spider venomous gland revealed by transcriptome analysis. Mol BioSyst. 2010, 6: 2403-2416. 10.1039/c004118a.View ArticlePubMed
- Ruiming Z, Yibao M, Yawen H, Zhiyong D, Yingliang W, Zhijian C, Wenxin L: Comparative venom gland transcriptome analysis of the scorpion Lychas mucronatus reveals intraspecific toxic gene diversity and new venomous components. BMC Genomics. 2010, 11: 452-10.1186/1471-2164-11-452.PubMed CentralView ArticlePubMed
- Hu H, Bandyopadhyay PK, Olivera BM, Yandell M: Characterization of the Conus bullatus genome and its venon-duct transcriptome. BMC Genomics. 2011, 12: 60-10.1186/1471-2164-12-60.PubMed CentralView ArticlePubMed
- Durban J, Juárez P, Angulo Y, Lomonte B, Flores-Diaz M, Alape-Girón A, Sasa M, Sanz L, Gutiérrez JM, Dopazo J, Conesa A, Calvete JJ: Profiling the venom gland transcriptomes of Costa Rican snakes by 454 pyrosequencing. BMC Genomics. 2011, 12: 259-10.1186/1471-2164-12-259.PubMed CentralView ArticlePubMed
- Gilles A, Meglécz E, Pech M, Ferreira S, Malausa T, Martin JF: Accuracy and quality assessment of 454 GS-FLX Titanium pyrosequencing. BMC Genomics. 2011, 12: 245-10.1186/1471-2164-12-245.PubMed CentralView ArticlePubMed
- Rodrigue S, Materna AC, Timberlake SC, Blackburn MC, Malmstrom RR, Aim EJ, Chisholm SW: Unlocking short read sequencing for metagenomics. PLoS One. 2010, 5: e11840-10.1371/journal.pone.0011840.PubMed CentralView ArticlePubMed
- Birol I, Jackman SD, Nielsen CB, Qian JQ, Varhol R, Stazyk G, Morin RD, Zhao Y, Hirst M, Schein JE, Horsman DE, Connors JM, Gascoyne RD, Marra MA, Jones SJM: De novo transcriptome assembly with ABySS. Bioinformatics. 2009, 25: 2872-2877. 10.1093/bioinformatics/btp367.View ArticlePubMed
- Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJM, Birol I: ABySS: a parallel assembler for short read sequence data. Genome Res. 2009, 19: 1117-1123. 10.1101/gr.089532.108.PubMed CentralView ArticlePubMed
- Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18: 821-829. 10.1101/gr.074492.107.PubMed CentralView ArticlePubMed
- Schulz MH, Zerbino DR, Vingron M, Birney E: Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics. 2012, 28: 1086-1092. 10.1093/bioinformatics/bts094.PubMed CentralView ArticlePubMed
- Feldmeyer B, Wheat CW, Krezdorn N, Rotter B, Pfenninger M: Short read Illumina data for the de novo assembly of a non-model snail species transcriptome (Radix balthica, Basommatophora, Pulmonata), and a comparison of assembler performance. BMC Genomics. 2011, 12: 317-10.1186/1471-2164-12-317.PubMed CentralView ArticlePubMed
- Dohm JC, Lottaz C, Borodina T, Himmelbauer H: Substantial biases in ultra-short read data sets from high-throughput DNA sequencing. Nucleic Acids Res. 2008, 36: e105-10.1093/nar/gkn425.PubMed CentralView ArticlePubMed
- Gibbs HL, Sanz L, Calvete JJ: Snake population venomics: proteomics-based analyses of individual variation reveals significant gene regulation effects on venom protein expression in Sistrurus rattlesnakes. J Mol Evol. 2009, 68: 113-125. 10.1007/s00239-008-9186-1.View ArticlePubMed
- Fox JW, Serrano SMT: Structural considerations of the snake venom metalloproteinases, key members of the M12 reprolysin family of metalloproteinases. Toxicon. 2005, 45: 969-985. 10.1016/j.toxicon.2005.02.012.View ArticlePubMed
- Fox JW, Serrano SMT: Snake venom metalloproteinases. Handbook of Venoms and Toxins of Reptiles. Edited by: Mackessy SP. 2010, CRC Press, Boca Raton, Florida, 95-113.
- Mackessy SP: Venom composition in rattlesnakes: trends and biological significance. The Biology of Rattlesnakes. Edited by: Hayes WK, Beaman KR, Cardwell MD, Bush SP. 2008, Loma Linda University Press, Loma Linda, California, 495-510.
- Du XY, Clemetson KJ: Reptile C-type lectins. Handbook of Venoms and Toxins of Reptiles. Edited by: Mackessy SP. 2010, CRC Press, Boca Raton, Florida, 359-375.
- Walker JR, Nagar B, Young NM, Hirama T, Rini JM: X-ray crystal structure of a galactose-specific C-type lectin possessing a novel decameric quaternary structure. Biochemistry. 2004, 43: 3783-3792. 10.1021/bi035871a.View ArticlePubMed
- Serrano SMT, Maroun RC: Snake venom serine proteinases: sequence homology vs. substrate specificity, a paradox to be solved. Toxicon. 2005, 45: 1115-1132. 10.1016/j.toxicon.2005.02.020.View ArticlePubMed
- Phillips DJ, Swenson SD, Francis S, Markland J: Thrombin-like snake venom serine proteinases. Handbook of Venoms and Toxins of Reptiles. Edited by: Mackessy SP. 2010, CRC Press, Boca Raton, Florida, 139-154.
- Lynch VJ: Inventing an arsenal: adaptive evolution and neofunctionalization of snake venom phospholipase A2genes. BMC Evol Biol. 2007, 7: 2-10.1186/1471-2148-7-2.PubMed CentralView ArticlePubMed
- Doley R, Zhou X, Kini RM: Snake venom phospholipase A2 enzymes. Handbook of Venoms and Toxins of Reptiles. Edited by: Mackessy SP. 2010, CRC Press, Boca Raton, Florida, 173-205.
- Rádis-Baptista G, Oguiura N, Hayashi MAF, Camargo ME, Grego KF, Oliveira EB, Yamane T: Nucleotide sequence of crotamine isoform precursors from a single South American rattlesnake (Crotalus durissus terrificus). Toxicon. 1999, 37: 973-984. 10.1016/S0041-0101(98)00226-8.View ArticlePubMed
- Oguiura N, Boni-Mitake M, Rádis-Baptista G: New view on crotamine, a small basic polypeptide myotoxin from South American rattlesnake venom. Toxicon. 2005, 46: 363-370. 10.1016/j.toxicon.2005.06.009.View ArticlePubMed
- Straight RC, Glenn JL, Wolt TB, Wolfe MC: Regional differences in content of small basic peptide toxins in the venoms of Crotalus adamanteus and Crotalus horridus. Comp Biochem Physiol B. 1991, 100: 51-58. 10.1016/0305-0491(91)90083-P.PubMed
- Tan NH, Fung SY: Snake venom L-amino acid oxidases. Handbook of Venoms and Toxins of Reptiles. Edited by: Mackessy SP. 2010, CRC Press, Boca Raton, Florida, 221-235.
- Heyborne WH, Mackessy SP: Cysteine-rich secretory proteins in reptile venoms. Handbook of Venoms and Toxins of Reptiles. Edited by: Mackessy SP. 2010, CRC Press, Boca Raton, Florida, 325-336.
- Yamazaki Y, Hyodo F, Morita T: Wide distribution of cysteine-rich secretory proteins in snake venoms: isolation and cloning of novel snake venom cysteine-rich secretory proteins. Arch Biochem Biophys. 2003, 412: 133-141. 10.1016/S0003-9861(03)00028-6.View ArticlePubMed
- Yamazaki Y, Morita T: Structure and function of snake venom cysteine-rich secretory proteins. Toxicon. 2004, 44: 227-231. 10.1016/j.toxicon.2004.05.023.View ArticlePubMed
- Pung YF, Wong PTH, Kumar PP, Hodgson WC, Kini RM: Ohanin, a novel protein from king cobra venom, induces hypolocomotion and hyperalgesia in mice. J Biol Chem. 2005, 280: 13137-13147.View ArticlePubMed
- Pung YF, Kumar SV, Rajagopalan N, Fry BG, Kumar PP, Kini RM: Ohanin, a novel protein from king cobra venom: its cDNA and genomic organization. Gene. 2006, 371: 246-256. 10.1016/j.gene.2005.12.002.View ArticlePubMed
- Junqueira-de-Azevedo ILM, Ching ATC, Carvalho E, Faria F, Nishiyama Jr 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, 173: 877-889. 10.1534/genetics.106.056515.PubMed CentralView ArticlePubMed
- Aird SD: Ophidian envenomation strategies and the role of purines. Toxicon. 2002, 40: 335-393. 10.1016/S0041-0101(01)00232-X.View ArticlePubMed
- Aird SD: The role of purine and pyrimidine nucleosides in snake venoms. Handbook of Venoms and Toxins of Reptiles. Edited by: Mackessy SP. 2010, CRC Press, Boca Raton, Florida, 393-419.
- Dhananjaya BL, Vishwanath BS, D’Souza CJM: Snake venom nucleases, nucleotidases, and phosphomonoesterases. Handbook of Venoms and Toxins of Reptiles. Edited by: Mackessy SP. 2010, CRC Press, Boca Raton, Florida, 155-171.
- Shafqat J, Zaidi ZH, Jörnvall H: Purification and characterization of a chymotrypsin Kunitz inhibitor type of polypeptide from the venom of cobra (Naja naja naja). FEBS Lett. 1990, 275: 6-8. 10.1016/0014-5793(90)81426-O.View ArticlePubMed
- Harrison RA, Ibison F, Wilbraham D, Wagstaff SC: Identification of cDNAs encoding viper venom hyaluronidases: cross-generic sequence conservation of full-length and unusually short variant transcripts. Gene. 2007, 392: 22-33. 10.1016/j.gene.2006.10.026.View ArticlePubMed
- Kemparaju K, Girish KS, Nagaraju S: Hyaluronidases, a neglected class of glycosidases from snake venom: beyond a spreading factor. Handbook of Venoms and Toxins of Reptiles. Edited by: Mackessy SP. 2010, CRC Press, Boca Raton, Florida, 237-258.
- Pawlak J, Kini RM: Snake venom glutaminyl cyclase. Toxicon. 2006, 48: 278-286. 10.1016/j.toxicon.2006.05.013.View ArticlePubMed
- Fry BG, Scheib H, van der Weerd L, Young B, McNaughtan J, Ramjan SFR, Vidal N, Poelmann RE, Norman JA: Evolution of an arsenal. Mol Cell Proteomics. 2008, 7: 215-246.View ArticlePubMed
- Rehana S, Kini RM: Molecular isoforms of cobra venom factor-like proteins in the venom of Austrelaps superbus. Toxicon. 2007, 50: 32-52. 10.1016/j.toxicon.2007.02.016.View ArticlePubMed
- Eggertsen G, Lind P, Sjöquist J: Molecular characterization of the complement activating protein in the venom of the Indian cobra (Naja n. siamensis). Mol Immunol. 1981, 18: 125-133. 10.1016/0161-5890(81)90078-X.View ArticlePubMed
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMed
- McCue MD: Cost of producing venom in three North American pitviper species. Copeia. 2006, 2006: 818-825. 10.1643/0045-8511(2006)6[818:COPVIT]2.0.CO;2.View Article
- Mackessy SP, Baxter LM: Bioweapons synthesis and storage: the venom gland of front-fanged snakes. Zool Anz. 2006, 245: 147-159. 10.1016/j.jcz.2006.01.003.View Article
- Wang Q, Zhang Z, Blackwell K, Carmichael GG: Vigilins bind to promiscuously A-to-I-edited RNAs and are involved in the formation of heterochromatin. Curr Biol. 2005, 15: 384-391. 10.1016/j.cub.2005.01.046.View ArticlePubMed
- Nishikura K: Functions and regulation of RNA editing by ADAR deaminases. Annu Rev Biochem. 2010, 79: 321-349. 10.1146/annurev-biochem-060208-105251.PubMed CentralView ArticlePubMed
- Hartl FU: Molecular chaperones in cellular protein folding. Nature. 1996, 381: 571-580. 10.1038/381571a0.View ArticlePubMed
- Fink AL: Chaperone-mediated protein folding. Physiol Rev. 1999, 79: 425-449.PubMed
- Young JC, Agashe VR, Siegers K, Hartl FU: Pathways of chaperone-mediated protein folding in the cytosol. Nat Rev Mol Cell Biol. 2004, 5: 781-791. 10.1038/nrm1492.View ArticlePubMed
- Finley D: Recognition and processing of ubiquitin-protein conjugates by the proteasome. Annu Rev Biochem. 2009, 78: 477-513. 10.1146/annurev.biochem.78.081507.101607.PubMed CentralView ArticlePubMed
- Buchberger A, Bukau B, Sommer T: Protein quality control in the cytosol and the endoplasmic reticulum: brothers in arms. Mol Cell. 2010, 40: 238-252. 10.1016/j.molcel.2010.10.001.View ArticlePubMed
- Bagola K, Mehnert M, Jarosch E, Sommer T: Protein dislocation from the ER. Biochim Biophys Acta. 2011, 1808: 925-936. 10.1016/j.bbamem.2010.06.025.View ArticlePubMed
- Huang KF, Chiou SH, Ko TP, Wang AHJ: Determinants of the inhibition of a Taiwan habu venom metalloproteinase by its endogenous inhibitors by X-ray crystallography and synthetic inhibitor analogues. Eur J Biochem. 2002, 269: 3047-3056. 10.1046/j.1432-1033.2002.02982.x.View ArticlePubMed
- Richards R, Trabi M, Johnson LA, de Jersey J, Masci PP, Lavin MF, St Pierre L: Cloning and characterization of novel cystatins from elapid snake venom glands. Biochimie. 2011, 93: 659-668. 10.1016/j.biochi.2010.12.008.View ArticlePubMed
- McCleary RJR, Heard DJ: Venom extraction from anesthetized Florida cottonmouths, Agkistrodon piscivorus conanti, using a portable nerve stimulator. Toxicon. 2010, 55: 250-255. 10.1016/j.toxicon.2009.07.030.View ArticlePubMed
- Rotenberg D, Bamberger ES, Kochva E: Studies on ribonucleic acid synthesis in the venom glands of Vipera palaestinae (Ophidia, Reptilia). Biochem J. 1971, 121: 609-612.PubMed CentralView ArticlePubMed
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-628. 10.1038/nmeth.1226.View ArticlePubMed
- Bentley DR, Balasubramanian S, Swerdlow HP, Smith GP, Milton J, Brown CG, Hall KP, Evers DJ, Barnes CL, Bignell HR, Boutell JM, Bryant J, Carter RJ, Cheetham RK, Cox AJ, Ellis DJ, Flatbush MR, Gormley NA, Humphray SJ, Irving LJ, Karbelashvili MS, Kirk SM, Li H, Liu X, Maisinger KS, Murray LJ, Obradovic B, Ost T, Parkinson ML, Pratt MR, Rasolonjatovo IMJ, Reed MT, Rigatti R, Rodighiero C, Ross MT, Sabot A, Sankar SV, Scally A, Schroth GP, Smith ME, Smith VP, Spiridou A, Torrance PE, Tzonev SS, Vermaas EH, Walter K, Wu X, Zhang L, Alam MD, Anastasi C, Aniebo IC, Bailey DMD, Bancarz IR, Banerjee S, Barbour SG, Baybayan PA, Benoit VA, Benson KF, Bevis C, Black PJ, Boodhun A, Brennan JS, Bridgham JA, Brown RC, Brown AA, Buermann DH, Bundu AA, Burrows JC, Carter NP, Castillo N, Chiara E, Catenazzi M, Chang S, Neil Cooley R, Crake NR, Dada OO, Diakoumakos KD, Dominguez-Fernandez B, Earnshaw DJ, Egbujor UC, Elmore DW, Etchin SS, Ewan MR, Fedurco M, Fraser LJ, Fuentes Fajardo KV, Furey WS, George D, Gietzen KJ, Goddard CP, Golda GS, Granieri PA, Green DE, Gustafson DL, Hansen NF, Harnish K, Haudenschild CD, Heyer NI, Hims MM, Ho JT, Horgan AM, Hoschler K, Hurwitz S, Ivanov DV, Johnson MQ, James T, Jones TAH, Kang GD, Kerelska TH, Kersey AD, Khrebtukova I, Kindwall AP, Kingsbury Z, Kokko-Gonzales PI, Kumar A, Laurent MA, Lawley CT, Lee SE, Lee X, Liao AK, Loch JA, Lok M, Luo S, Mammen RM, Martin JW, McCauley PG, McNitt P, Mehta P, Moon KW, Mullens JW, Newington T, Ning Z, Ng BL, Novo SM, O’Neill MJ, Osborne MA, Osnowski A, Ostadan O, Paraschos LL, Pickering L, Pike AC, Pike AC, Pinkard DC, Pliskin DP, Podhasky J, Quijano VJ, Raczy C, Rae VH, Rawlings SR, Rodriguez AC, Roe PM, Rogers J, Bacigalupo MCR, Romanov N, Romieu A, Roth RK, Rourke NJ, Ruediger ST, Rusman E, Sanches-Kuiper RM, Schenker MR, Seoane JM, Shaw RJ, Shiver MK, Short SW, Sizto NL, Sluis JP, Smith MA, Sohna JES, Spence EJ, Stevens K, Sutton N, Szajkowski L, Tregidgo CL, Turcatti G, vandeVondele S, Verhovsky Y, Virk SM, Wakelin S, Walcott GC, Wang J, Worsley GJ, Yan J, Yau L, Zuerlein M, Rogers J, Mullikin JC, Hurles ME, McCooke NJ, West JS, Oaks FL, Lundberg PL, Klenerman D, Durbin R, Smith AJ: Accurate whole human genome sequencing using reversible terminator chemistry. Nature. 2008, 456: 53-59. 10.1038/nature07517.PubMed CentralView ArticlePubMed
- Robertson G, Schein J, Chiu R, Corbett R, Field M, Jackman SD, Mungall K, Lee S, Okada HM, Qian JQ, Griffith M, Raymond A, Thiessen N, Cezard T, Butterfield YS, Newsome R, Chan SK, She R, Varhol R, Kamoh B, Prabhu AL, Tam A, Zhao Y, Moore RA, Hirst M, Marra MA, Jones SJM, Hoodless PA, Birol I: De novo assembly and analysis of RNA-seq data. Nat Methods. 2010, 7: 909-912. 10.1038/nmeth.1517.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.