- Research article
- Open Access
Transcriptome analysis of the venom gland of the Mexican scorpion Hadrurus gertschi (Arachnida: Scorpiones)
© Schwartz et al; licensee BioMed Central Ltd. 2007
Received: 17 March 2007
Accepted: 16 May 2007
Published: 16 May 2007
Scorpions like other venomous animals posses a highly specialized organ that produces, secretes and disposes the venom components. In these animals, the last postabdominal segment, named telson, contains a pair of venomous glands connected to the stinger. The isolation of numerous scorpion toxins, along with cDNA-based gene cloning and, more recently, proteomic analyses have provided us with a large collection of venom components sequences. However, all of them are secreted, or at least are predicted to be secretable gene products. Therefore very little is known about the cellular processes that normally take place inside the glands for production of the venom mixture. To gain insights into the scorpion venom gland biology, we have decided to perform a transcriptomic analysis by constructing a cDNA library and conducting a random sequencing screening of the transcripts.
From the cDNA library prepared from a single venom gland of the scorpion Hadrurus gertschi, 160 expressed sequence tags (ESTs) were analyzed. These transcripts were further clustered into 68 unique sequences (20 contigs and 48 singlets), with an average length of 919 bp. Half of the ESTs can be confidentially assigned as homologues of annotated gene products. Annotation of these ESTs, with the aid of Gene Ontology terms and homology to eukaryotic orthologous groups, reveals some cellular processes important for venom gland function; including high protein synthesis, tuned posttranslational processing and trafficking. Nonetheless, the main group of the identified gene products includes ESTs similar to known scorpion toxins or other previously characterized scorpion venom components, which account for nearly 60% of the identified proteins.
To the best of our knowledge this report contains the first transcriptome analysis of genes transcribed by the venomous gland of a scorpion. The data were obtained for the species Hadrurus gertschi, belonging to the family Caraboctonidae. One hundred and sixty ESTs were analyzed, showing enrichment in genes that encode for products similar to known venom components, but also provides the first sketch of cellular components, molecular functions, biological processes and some unique sequences of the scorpion venom gland.
Scorpion venoms are very complex mixtures with hundreds of different components produced by the highly specialized venom glands. The most prominent components of scorpion venoms are the peptides responsible for the neurotoxic effects associated with their sting, for which more than 350 different have been described (extensive databases can be found in Tox-Prot  and SCORPION ). Most of these toxins are structurally related disulphide-rich short proteins (23–75 amino acid residues long), which affect cellular communication by modulating Na+ or K+ ion-channels permeability . Due to their importance in scorpion envenomation and their usefulness as molecular and pharmacological probes for studying ion-channels, most of the work performed to date are focused at these neurotoxins, with relative few other components ever described; among which are heterodimeric phospholipases A2 (v.gr. [4–6]), non-disulphide short peptides with cytolytic activity and a few other functions [7, 8]. Recent proteomic analyses [9–16] have documented the overall composition for nine scorpion species, all of them from the family Buthidae and most of them belonging to the Tityus genus. These analyses confirmed the gross estimation of an average of one hundred different proteins in each one of the venoms . Approximately half of them comprehend components with molecular masses in the range of commonly found scorpion toxins (2,000–8,000 Da). These numbers contrast heavily with the known universe of protein components (near four hundreds) described to exist in scorpion venoms, from which only about 12% are not classified within the known scorpion toxin families.
Further insights into scorpion venom compositions have been achieved by gene cloning by PCR-based methods conducted with cDNA libraries. For example, almost one hundred toxin precursors have been sequenced from venom gland libraries of the buthid scorpion Mesobuthus martensii (v.gr. [18–20]). Unfortunately the spectrum of sequences obtained through PCR-based approach is limited by the specificity of the PCR primers used. It is worth noticing that although PCR-based methods along with the abundant isolation and characterization of scorpion toxins and, more recently, proteomic profiling of whole venoms, have provided us with a large number of sequences, all these components are secreted from the venom glands. Little is known about the biological processes that are taking place inside the venom gland cells. Therefore, we elected to use a transcriptome approach to improve the understanding of the composition of Hadrurus gertschi venom gland.
The scorpion H. gertschi Soleglad (1976) belongs to the family Caraboctonidae  and is considered no dangerous to humans. H. gertschi is endemic to Mexico, occurring exclusively in the State of Guerrero, and lives underground in tunnels excavated in the soil. From the venom of this scorpion few components have been isolated and studied: hadrurin, an antimicrobial and cytolytic peptide ; HgeTx1, a K+ channel blocker ; hadrucalcine, a peptide capable of activating skeletal Ryanodine receptors [Schwartz et al., in preparation], and; the precursors HgeScplp and HgeβKTx, which encode for long-chain peptides similar to Scorpine and βKTx's, respectively . Although hadrurin was reported as component of H. aztecus venom , the specimens used in that work were not taxonomically identified and latter it was realized that scorpions from that geographical region should be named H. gertschi; this species assignment was confirmed by identification of relevant taxonomic keys in the specimens.
In the present work we randomly generated and analyzed 160 expressed sequence tags (ESTs) from a cDNA library of the venom gland of H. gertschi. These 160 ESTs corresponded to 0.15% of the whole non-amplified cDNA library and were generated from a non-normalized cDNA library. After clustering the resulting dataset, we identified transcripts possibly associated with different cellular functions. The possible roles of some of the transcripts are discussed, although many have unknown functions. Furthermore, we present 8 full length sequences of new toxins.
Single-pass gene sequencing from cDNA libraries is an affordable strategy to mine the transcript profile of a given tissue . This strategy has been used to analyze the transcripts profiles for few other venomous organisms, such as cnidarians [26, 27], cone snails , fishes , snakes (v.gr. [30–32]) and spiders . To the best of our knowledge, this is the first report of an ESTs strategy conducted with any scorpion venom gland. Moreover, this is the first comprehensive molecular study of a non-buthid scorpion, which could serve for comparative purposes when studying the details of the process by which buthid scorpions have been assembling their neurotoxic arsenal.
cDNA library and EST analysis
Identification of the transcripts predicted to be involved in common cellular processes and those similar to known venom components. The putative identity corresponds to the eukaryotic orthologous group (KOG), as detailed in Methods.
Gene products predicted to be involved in common cellular processes
KOG1376 Alpha tubulin
KOG3412 60S ribosomal protein L28
KOG0279 G protein beta subunit-like protein
KOG1954 Endocytosis/signaling protein EHD1
KOG3418 60S ribosomal protein L27
KOG3449 60S acidic ribosomal protein P2
KOG0714 Molecular chaperone (DnaJ superfamily)
KOG3458 NADH:ubiquinone oxidoreductase, NDUFA8/PGIV/19 kDa subunit
KOG0863 20S proteasome, regulatory subunit alpha type PSMA1/PRE5
KOG3311 Ribosomal protein S18 (40S)
KOG1629 Bax-mediated apoptosis inhibitor TEGT/BI-1
KOG0898 40S ribosomal protein S15
KOG2597 Predicted aminopeptidase of the M17 family
KOG3752 Ribonuclease H
KOG0876 Manganese superoxide dismutase
KOG4075 Cytochrome c oxidase, subunit IV/COX5b
KOG2667 COPII vesicle protein
KOG2403 Succinate dehydrogenase, flavoprotein subunit
KOG2486 Predicted GTPase
KOG4604 Uncharacterized conserved protein
Gene products similar to known venom components
α-KTx 6 subfamily
Novel NDBP group 5
Novel NDBP group 5
Novel NDBP group 3
Novel group III heterodimeric phospholipase
KOG4295 Serine proteinase inhibitor (KU family)
Similarity searches and sequence annotation
All sequences were submitted for blastn and blastx searches against nr database; an e-value < 10-5 was used as cut-off for confidential homologue detection. In addition to the ESTs for which no clear ORF have been identified, nearly 19% of the new sequences (in 3 contigs and 10 singlets) do not match with any entry in the database. Altogether, unassigned ESTs account for close to 48% of the total dataset, a value similar to other transcriptome studies which show values varying from 13% to 56% of non matched sequences [26–33, 35]. Noteworthy, these putative gene products signify a source of new information about scorpion venom gland specific genes. In addition to ESTs with no database match, 2 reads presented identity with sequences that have been already described but with no functional assessment, hereby named unknown proteins. One of these, HGE032|2273, is the clone which codes for a protein similar to the CG9896-like proteins identified in the scorpions Mesobuthus gibbosus and M. cyprius. The other one, HGE033|2208, encodes for a scorpion homologue of short proteins conserved in eukaryotic organisms (KOG4604, pfam04418.6), but whose function is not known yet.
It is worth noticing that in the group of identified proteins, toxins account for 59.7% of the transcripts (30% of the unique sequences). The distribution of all ESTs is depicted in Figure 2, it can be observed that ESTs coding for toxins are well represented in the H. gertschi venom gland transcriptome. Further, considering the non-matched ESTs as a possible source for new toxins, it can be assumed that these molecules are preferentially expressed over proteins related with the other cellular functions.
GO-sorted annotated sequences
Transcriptome analysis suggest cellular processes relevant for scorpion venom glands function
Although our sampling of the venom glands library is still incomplete, the diversity and nature of the annotated transcripts provide the first glimpse about molecular processes taking part in the scorpion venom gland cells. Since we constructed a non-amplified library, it could be expected that clone number reflects the actual prevalence of a given transcript. Moreover, by extension, different transcripts belonging to the same – confidentially assigned – GO category might suggest this category as important within the biological processes of scorpion venom glands.
For example, intuitively, the venom glands should support high protein synthesis and secretion in order to produce the large amounts of, secreted and renewable, venom proteins. In concordance, 8.2% of the total transcripts and 16.2% of the unique sequences match with either ribosomal components (1 contig and 4 singlets) or proteins involved in cellular trafficking (2 contigs and 1 singlet). Both processes are energetically costly and, consistently, 4.1% of whole ESTs and 7.4% of identified protein precursors match with components of the energy-producer organelle mitochondrion, whereas 2.7% and 4.4%, respectively, are putative homologues of proteins directly involved in the energy-producer oxidative phosphorilation or tricarboxilyc acid cycle. Indeed, scorpions whose venoms were artificially depleted shows increased oxygen consumption .
The importance of correct protein processing in the context of scorpion venom gland is emphasized by the presence of transcripts encoding for proteins involved in correct folding (HGE008|2233), posttranslational processing (HGE014|2453 and HGE018|2232) or proteasome-dependent degradation of proteins (HGE010|2258). One of these (HGE014|2453), match with aminopeptidases of the M17 family (KOG2597), which are exopeptidases involved in the processing and regular turnover of intracellular proteins, although their precise role in cellular metabolism is unclear. In particular aminopeptidases of the M17 family cleave leucine residues from the N-terminal of polypeptide chains, but substantial rates are evident for all amino acids . We predicted, by SignalP 3.0, its signal peptide-processing sequence (VAS-LK) suggesting that this transcript is secreted as a venom component, and indeed it could be important for posttranslational modifications of venomous gland components. Moreover, HGE018|2232 match with the glycosylating enzymes β-1,4-mannosyltransferases (KOG2941); which is consistent with the presence of glycosylated proteins in scorpion venoms (v.gr. ).
An interesting finding in our database was a transcript encoding the ribonuclease H domain of a non-LTR retrotransposon (nLTRrt) of clade R1 (sensu ). This nLTRrt clade is usually found inserted into telomeres and has been identified in several arthropods (including one arachnid) and some fungi . It is worth noticing that mobile elements and their remnants account for large proportions of most eukaryotic genomes, in which they have had central roles in genome evolution and hypervariation. The expression of transposases indicates that mobile elements might contribute to the diversification of venom toxins. Recently, Glushkov et al.,  reported the search for nLTRrt in 22 scorpion species, in which degenerate oligonucleotides based on consensus sequences of seven clades of nLTRrt were used to a PCR-based fishing approach. Unfortunately, even though these authors reported that PCR products of the expected size where obtained with R1-based degenerate primers, they only presented data for CR1, I and Jockey nLTRrt clades.
Toxins and other venom components
Our H. gertschi scorpion venom gland library is clearly enriched on toxin-like sequences, with more than 17% (in 4 contigs and 1 singlet) of the sequenced ESTs being similar to known families of scorpion toxins. Another 14% of the total ESTs (in 4 contigs) encodes for precursors which are homologues of previously characterized non-toxin scorpion venom components (see below). Considering that all these sequences contain putative signal peptides – identified by SignalP 3.0  – and their relative abundance, we suggest that these ESTs may encode secreted venom components. In fact, 3 of these clusters encode peptides already found in the venom of H. gertschi, which are currently being studied by our group [24 and unpublished].
β-KTx and scorpine-like peptides
One cluster (2 reads) coding the HgeβKTx and another (4 reads) coding the Hge scorpine like were identified in the transcriptome of H. gertschi and their sequence have already been reported . Here we present a distinct EST, homologous to the scorpine-like group of long-chain three disulphide-bridged scorpion venom peptides, named Hge scorpine like 2 (Figure 4b and Supplementary Figure 1C in Additional file 1).
Bradykinin-potentiating peptide like
One cluster (HGE029|NDBP_3.7, 2 reads; see Supplementary Figure 2B in Additional file 1) encodes a homologue of the bradykinin-potentiating peptide precursor (BmK3 or BmKpbb) from the scorpion Mesobuthus martensii, classified as NDBP-3.3 . The angiotensin-bradykinin system is a central hormonal system for the regulation of blood pressure. The angiotensin-converting enzyme (ACE) converts angiotensin I to angiotensin II and degrades bradykinin. Bradykinin potentiating peptides have been isolated from Tityus serrulatus (peptide T ) and Buthus occitanus (K12 ). Peptide T, a 13-amino acid linear peptide, potentiates the contractile activity of bradykinin on isolated smooth muscle, inhibits the hydrolysis of bradykinin by ACE, and enhances the depressor effect of bradykinin on arterial blood pressure in the anesthetized rats . Peptide K12 displays similar bradykinin potentiating activities . BmKbpp was identified from B. martensii Karsch by cDNA cloning based on the peptide K12 amino acid sequence . The last 21 residues of C-terminal region of BmKbpp showed 57% similarity with peptide K12. Based on the fact that BmKbpp also exhibits high similarity with parabutoporin and others antimicrobial peptides from scorpions, it was suggested that BmKbpp may be a molecule with a dual-function, and that the BmKbpp precursor may be processed in two alternative ways to produces two different mature molecules: BmKbpp and a peptide with only the C-terminal 21 residues of BmKbpp . Figure 5b shows the bradykinin-potentiating peptide like from H. gertschi, here named NDBP-3.7.
Kunitz-type carboxypeptidase inhibitor
Gene cloning of animal toxins has been extensively performed by PCR method, using primers deduced from direct protein sequencing, usually by Edman degradation or mass spectrometry analysis. These studies are aimed at the isolation of specific active components. However, this approach is not entirely suitable for search of unforeseen components that could be present in the venomous gland under study. The strategy is biased by the fact that only those genes that are sharing sequence similarities are usually discovered by this technique. For this reason, we adopted the molecular approach of generating and analyzing ESTs from the H. gertschi venom gland as the strategy to produce a general overview of the venom gland transcriptome. This strategy confirms the highly specialized nature of scorpion venom glands as toxin-producer, allowing the description, for the first time, of putative proteins that certainly are involved in cellular processes relevant for the venom glands' function. Additionally, the unguided mining also reveals novel predicted venom components, highlighting the usefulness of the transcriptome approach to improve venom profiling.
cDNA library construction
A cDNA library was constructed from total RNA extracted from a single telson of a H. gertschi scorpion. The scorpion was milked 5 days before RNA extraction. For RNA isolation the 'Total RNA Isolation System' of Promega (Madison, WI) was used. With this material a full-length cDNA phagemid library was prepared using the SMART cDNA Library Construction Kit (CLONTECH Lab., Palo Alto, CA). The titre of the non-amplified cDNA library obtained was 2.8 × 105 cfu/mL with 99% recombinant clones. For the PCR the oligonucleotides TriplEx2-5' (from CLONTECH Lab, Palo Alto, CA) and CD3/3' (5'-AAT CTA GAG GCC GAG GCG GCC GAC ATG-3') , designed on the basis of the CDSIII/3' sequence tag used for library construction, were used as primers.
DNA sequencing and bioinformatic analyses
Selected plasmids were isolated according to a standard alkaline lysis protocol, and single-pass sequencing of the 5'-termini was conducted with the primer TriplEx2-5' (CLONTECH Lab, Palo Alto, CA) using an automatic machine (Model 3100, Applied Biosystems, Foster city, CA) according to the manufacturer's instructions. The nucleotide sequences obtained in this work are deposited in dbEST (GenBank: EL698878–EL698945). To extract the high quality sequence region, the ESTs were subjected to the Phred program  with the window length set to 100 and the standard quality to 20. The CrossMatch program was used to remove vector and E. coli DNA sequences. ESTs that shared an identity of > 95 out of 100 nucleotides were assembled in contiguous sequences with the CAP3 program . All these bioinformatics analysis were simultaneously run at  using default setup. H. gertschi venom gland ESTs (clusters and singlets) were searched against nr database  using blastx and blastn algorithms  with an e-value cutoff set to < 10-5 to identify putative functions of the new ESTs. The signal peptide was predicted with the SignalP 3.0 program . Multiple sequences alignments were obtained using T-COFFEE  or CLUSTAL_X . The pairwise identities were calculated with BioEdit . In order to evaluate the most expressed ESTs, a less stringent alignment using a shared identity of > 65 out of 100 nucleotides was used.
With the aim of providing more useful sequence annotations for comparative studies, we selected to identify the eukaryotic orthologous group  to which each sequence would belong, instead of reporting the – most often used – highest scored BLAST hit. The best match for each sequence with KOG database was identified by Kognitor – a BLAST-based web server that search within KOGs only – and validated by independent PSI-BLAST searches . The rationale followed was that PSI-BLAST should return all members of the given KOG within the first two iterations. Furthermore, in order to improve sequence annotation and to gain insights into the cellular processes where each sequence could be involved, the ESTs were submitted to GOblet  and AmiGO  web servers. The aim of this procedure was to identify relevant Gene Ontology terms. Both programs perform blastx searches against UniProt (GOblet) and all GO-annotated protein databases (AmiGO). Cross-validation with identified KOGs was performed by inspection of the retrieved hits in order to identify homologues of the corresponding KOGs. Subsequently the GenBank entries of the given KOG were used as query for AmiGO. GO terms were assigned only when the novel sequence and at least one sequence of each KOGs retrieved coincident hits (usually five or more out of the ten highest-scoring hits).
The authors are indebted to Dr Blanca I. García and M.Sc. Timoteo Olamendi-Portugal for their technical assistance. This work was partially financed by grants from CNPq – Brazil number 201294/2003-6 (to EFS), CONACyT-Mexico and CNPq-Brazil (to LDP and EFS) and grants from Dirección General de Asuntos del Personal Académico (Universidad Nacional Autónoma de México) number IN227507 and CONACyT-Fondo Salud number 14109 (Mexican Government) to LDP.
- Jungo F, Bairoch A: Tox-Prot, the toxin protein annotation program of the Swiss-Prot protein knowledgebase. Toxicon. 2005, 45: 293-301. 10.1016/j.toxicon.2004.10.018.PubMedView ArticleGoogle Scholar
- Tan PT, Veeramani A, Srinivasan KN, Ranganathan S, Brusic V: SCORPION2: a database for structure-function analysis of scorpion toxins. Toxicon. 2006, 47: 356-363. 10.1016/j.toxicon.2005.12.001.PubMedView ArticleGoogle Scholar
- Possani LD, Rodríguez de la Vega RC: Scorpion venom peptides. Handbook of Biologically Active Peptides. Edited by: Kastin AJ. 2006, San Diego, Academic Press, 339-354.View ArticleGoogle Scholar
- Ramanaiah M, Parthasarathy PR, Venkaiah B: Purification and properties of phospholipase A2 from the venom of scorpion, (Heterometrus fulvipes). Biochem Int. 1990, 20: 931-940.PubMedGoogle Scholar
- Zamudio FZ, Conde R, Arevalo C, Becerril B, Martin BM, Valdivia HH, Possani LD: The mechanism of inhibition of ryanodine receptor channels by imperatoxin I, a heterodimeric protein from the scorpion Pandinus imperator. J Biol Chem. 1997, 272: 11886-11894. 10.1074/jbc.272.18.11886.PubMedView ArticleGoogle Scholar
- Valdez-Cruz NA, Batista CV, Possani LD: Phaiodactylipin, a glycosylated heterodimeric phospholipase A from the venom of the scorpion Anuroctonus phaiodactylus. Eur J Biochem. 2004, 271: 1453-1464. 10.1111/j.1432-1033.2004.04047.x.PubMedView ArticleGoogle Scholar
- Zeng XC, Corzo G, Hahin R: Scorpion venom peptides without disulfide bridges. IUBMB Life. 2005, 57: 13-21.PubMedView ArticleGoogle Scholar
- Luo F, Zeng XC, Hahin R, Cao ZJ, Liu H, Li WX: Genomic organization of four novel nondisulfide-bridged peptides from scorpion Mesobuthus martensii Karsch: gaining insight into evolutionary mechanism. Peptides. 2005, 26: 2427-2433. 10.1016/j.peptides.2005.06.008.PubMedView ArticleGoogle Scholar
- Pimenta AM, Stocklin R, Favreau P, Bougis PE, Martin-Eauclaire MF: Moving pieces in a proteomic puzzle: mass fingerprinting of toxic fractions from the venom of Tityus serrulatus (Scorpiones, Buthidae). Rapid Commun Mass Spectrom. 2001, 15: 1562-1572. 10.1002/rcm.415.PubMedView ArticleGoogle Scholar
- Batista CV, del Pozo L, Zamudio FZ, Contreras S, Becerril B, Wanke E, Possani LD: Proteomics of the venom from the Amazonian scorpion Tityus cambridgei and the role of prolines on mass spectrometry analysis of toxins. J Chromatogr B. 2004, 803: 55-66. 10.1016/j.jchromb.2003.09.002.View ArticleGoogle Scholar
- Diego-Garcia E, Batista CV, Garcia-Gomez BI, Lucas S, Candido DM, Gomez-Lagunas F, Possani LD: The Brazilian scorpion Tityus costatus Karsch: genes, peptides and function. Toxicon. 2005, 45: 273-283. 10.1016/j.toxicon.2004.10.014.PubMedView ArticleGoogle Scholar
- Barona J, Batista CV, Zamudio FZ, Gomez-Lagunas F, Wanke E, Otero R, Possani LD: Proteomic analysis of the venom and characterization of toxins specific for Na+- and K+ -channels from the Colombian scorpion Tityus pachyurus. Biochim Biophys Acta. 2006, 1764: 76-84.PubMedView ArticleGoogle Scholar
- Batista CV, D'Suze G, Gomez-Lagunas F, Zamudio FZ, Encarnacion S, Sevcik C, Possani LD: Proteomic analysis of Tityus discrepans scorpion venom and amino acid sequence of novel toxins. Proteomics. 2006, 6: 3718-3727. 10.1002/pmic.200500525.PubMedView ArticleGoogle Scholar
- Borges A, Garcia CC, Lugo E, Alfonzo MJ, Jowers MJ, Op den Camp HJ: Diversity of long-chain toxins in Tityus zulianus and Tityus discrepans venoms (Scorpiones, Buthidae): molecular, immunological, and mass spectral analyses. Comp Biochem Physiol C. 2006, 142: 240-252.Google Scholar
- Caliskan F, Garcia BI, Coronas FI, Batista CV, Zamudio FZ, Possani LD: Characterization of venom components from the scorpion Androctonus crassicauda of Turkey: peptides and genes. Toxicon. 2006, 48: 12-22. 10.1016/j.toxicon.2006.04.003.PubMedView ArticleGoogle Scholar
- Nascimento DG, Rates B, Santos DM, Verano-Braga T, Barbosa-Silva A, Dutra AA, Biondi I, Martin-Eauclaire MF, De Lima ME, Pimenta AM: Moving pieces in a taxonomic puzzle: venom 2D-LC/MS and data clustering analyses to infer phylogenetic relationships in some scorpions from the Buthidae family (Scorpiones). Toxicon. 2006, 47: 628-639. 10.1016/j.toxicon.2006.01.015.PubMedView ArticleGoogle Scholar
- Possani LD, Becerril B, Delepierre M, Tytgat J: Scorpion toxins specific for Na+-channels. Eur J Biochem. 1999, 264: 287-300. 10.1046/j.1432-1327.1999.00625.x.PubMedView ArticleGoogle Scholar
- Goudet C, Chi CW, Tytgat J: An overview of toxins and genes from the venom of the Asian scorpion Buthus martensi Karsch. Toxicon. 1999Google Scholar
- Zeng XC, Wang SX, Zhu Y, Zhu SY, Li WX: Identification and functional characterization of novel scorpion venom peptides with no disulfide bridge from Buthus martensii Karsch. Peptides. 2004, 25: 143-150. 10.1016/j.peptides.2003.12.003.PubMedView ArticleGoogle Scholar
- Zeng XC, Luo F, Li WX: Molecular dissection of venom from Chinese scorpion Mesobuthus martensii: identification and characterization of four novel disulfide-bridged venom peptides. Peptides. 2006, 27: 1745-1754. 10.1016/j.peptides.2006.01.012.PubMedView ArticleGoogle Scholar
- Soleglad ME, Fet V: High-level systematics and phylogeny of the extant scorpions (Scorpiones: Orthosterni). Euscorpius. 2003, 11-Google Scholar
- Torres-Larios A, Gurrola GB, Zamudio FZ, Possani LD: Hadrurin, a new antimicrobial peptide from the venom of the scorpion Hadrurus aztecus. Eur J Biochem. 2000, 267: 5023-31. 10.1046/j.1432-1327.2000.01556.x.PubMedView ArticleGoogle Scholar
- Schwartz EF, Schwartz CA, Gomez-Lagunas F, Zamudio FZ, Possani LD: HgeTx1, the first K+-channel specific toxin characterized from the venom of the scorpion Hadrurus gertschi Soleglad. Toxicon. 2006, 48: 1046-1053. 10.1016/j.toxicon.2006.08.009.PubMedView ArticleGoogle Scholar
- Diego-Garcia E, Schwartz EF, D'Suze G, Gonzalez SA, Batista CV, Garcia BI, Rodríguez de la Vega RC, Possani LD: Wide phylogenetic distribution of Scorpine and long-chain beta-KTx-like peptides in scorpion venoms: identification of "orphan" components. Peptides. 2007, 28: 31-37. 10.1016/j.peptides.2006.06.012.PubMedView ArticleGoogle Scholar
- Whitton C, Daub J, Thompson M, Blaxter M: Expressed sequence tags: medium-throughput protocols. Methods Mol Biol. 2004, 270: 75-92.PubMedGoogle Scholar
- Yang Y, Cun S, Xie X, Lin J, Wei J, Yang W, Mou C, Yu C, Ye L, Lu Y, Fu Z, Xu A: EST analysis of gene expression in the tentacle of Cyanea capillata. FEBS Lett. 2003, 538: 183-191. 10.1016/S0014-5793(03)00164-9.PubMedView ArticleGoogle Scholar
- Hydra EST database. [http://www.hydrabase.org]
- Pi C, Liu J, Peng C, Liu Y, Jiang X, Zhao Y, Tang S, Wang L, Dong M, Chen S, Xu A: Diversity and evolution of conotoxins based on gene expression profiling of Conus litteratus. Genomics. 2006, 88: 809-819. 10.1016/j.ygeno.2006.06.014.PubMedView ArticleGoogle Scholar
- Magalhaes GS, Junqueira-de-Azevedo ILM, Lopes-Ferreira M, Lornzini DM, Ho PL, Moura-da-Silva AM: Transcriptome analysis of expressed sequence tags from the venom glands of the fish Thalassophryne nattereri. Biochimie. 2006, 88: 693-699. 10.1016/j.biochi.2005.12.008.PubMedView ArticleGoogle Scholar
- Cidade DA, Simao TA, Davila AM, Wagner G, Junqueira-de-Azevedo Ide L, 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.PubMedView ArticleGoogle Scholar
- Wagstaff SC, Harrison RA: Venom gland EST analysis of the saw-scaled viper, Echis ocellatus, reveals novel alpha9beta1 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.PubMedView ArticleGoogle Scholar
- 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 structure and functional aspects using expressed sequence tags. BMC Genomics. 2006, 7: 152-10.1186/1471-2164-7-152.PubMed CentralPubMedView ArticleGoogle Scholar
- Kozlov S, Malyavka A, McCutchen B, Lu A, Schepers E, Herrmann R, Grishin E: A novel strategy for the identification of toxinlike structures in spider venom. Proteins. 2005, 59: 131-140. 10.1002/prot.20390.PubMedView ArticleGoogle Scholar
- Audic S, Claverie J-M: The significance of digital gene expression profiles. Genome Res. 1997, 7: 986-995.PubMedGoogle Scholar
- Dana AN, Hillenmeyer ME, Lobo NF, Kern MK, Romans PA, Collins FH: Differential gene expression in abdomens of the malaria vector mosquito, Anopheles gambiae, after sugar feeding, blood feeding and Plasmodium berghei infection. BMC Genomics. 7: 119-10.1186/1471-2164-7-119.Google Scholar
- Schmitt AO, Specht T, Beckmann G, Dahl E, Pilarsky CP, Hinzmann B, Rosenthal A: Exhaustive mining of EST libraries for genes differentially expressed in normal and tumour tissues. Nucleic Acids Res. 1999, 27: 4251-4260. 10.1093/nar/27.21.4251.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu D, Graber JH: Quantitative comparison of EST libraries requires compensation for systematic biases in cDNA generation. BMC Genomics. 2006, 7: 77-10.1186/1471-2164-7-77.View ArticleGoogle Scholar
- Tatusov RL, Fedorova ND, Jackson JW, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, Rao BS, Smirnov S, Sverdlov AV, Vasudevan S, Wolf YI, Yin JJ, Natale DA: The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003, 4: 41-10.1186/1471-2105-4-41.PubMed CentralPubMedView ArticleGoogle Scholar
- The Gene Ontology Consortium: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralView ArticleGoogle Scholar
- Groth D, Lehrach H, Hennig S: GOblet: a platform for Gene Ontology annotation of anonymous sequence data. Nucleic Acids Res. 2004, 32: W313-W319. 10.1093/nar/gkh406.PubMed CentralPubMedView ArticleGoogle Scholar
- Nisani Z, Dunbar SG, Hayes WK: Cost of venom regeneration in Parabuthus transvaalicus (Arachnida: Buthidae). Comp Biochem Physiol A. 2007, 147: 509-513. 10.1016/j.cbpa.2007.01.027.View ArticleGoogle Scholar
- Matsui M, Fowler JH, Walling LL: Leucine aminopeptidases: diversity in structure and function. Biol Chem. 2006, 387: 1535-1544. 10.1515/BC.2006.191.PubMedView ArticleGoogle Scholar
- Malik HS, Burke WD, Eickbush TH: The age and evolution of non-LTR retrotransposable elements. Mol Biol Evol. 1999, 16: 793-805.PubMedView ArticleGoogle Scholar
- Malik HS: Ribonuclease H evolution in retrotransposable elements. Cytogenet Genome Res. 2005, 110: 392-401. 10.1159/000084971.PubMedView ArticleGoogle Scholar
- Glushkov S, Novikova O, Blinov A, Fet V: Divergent non-LTR retrotransposon lineages from the genomes of scorpions (Arachnida: Scorpiones). Mol Genet Genomics. 2006, 275: 88-296. 10.1007/s00438-005-0079-3.View ArticleGoogle Scholar
- SignalP 3.0 Server. [http://www.cbs.dtu.dk/services/SignalP/]
- Tytgat J, Chandy KG, Garcia ML, Gutman GA, Martin-Eauclaire MF, van der Walt JJ, Possani LD: A unified nomenclature for short-chain peptides isolated from scorpion venoms: alpha-KTx molecular subfamilies. Trends Pharmacol Sci. 1999, 20: 444-447. 10.1016/S0165-6147(99)01398-X.PubMedView ArticleGoogle Scholar
- Rodríguez de la Vega RC, Possani LD: Current views on scorpion toxins specific for K+-channels. Toxicon. 2004, 43: 865-875. 10.1016/j.toxicon.2004.03.022.PubMedView ArticleGoogle Scholar
- UniProt Consortium: The Universal Protein Resource (UniProt). Nucleic Acids Res. 2007, 35: D193-197. 10.1093/nar/gkl929.View ArticleGoogle Scholar
- Zhu S, Huys I, Dyason K, Verdonck F, Tytgat J: Evolutionary trace analysis of scorpion toxins specific for K-channels. Proteins. 2004, 54: 361-370. 10.1002/prot.10588.PubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman 50DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralPubMedView ArticleGoogle Scholar
- Bagdany M, Batista CV, Valdez-Cruz NA, Somodi S, Rodríguez de la Vega RC, Licea AF, Varga Z, Gaspar R, Possani LD, Panyi G: Anuroctoxin, a new scorpion toxin of the alpha-KTx 6 subfamily, is highly selective for Kv1.3 over IKCa1 ion channels of human T lymphocytes. Mol Pharmacol. 2005, 67: 1034-1044. 10.1124/mol.104.007187.PubMedView ArticleGoogle Scholar
- Dai L, Corzo G, Naoki H, Andriantsiferana M, Nakajima T: Purification, structure-function analysis, and molecular characterization of novel linear peptides from scorpion Opisthacanthus madagascariensis. Biochem Biophys Res Commun. 2002, 293: 1514-1522. 10.1016/S0006-291X(02)00423-0.PubMedView ArticleGoogle Scholar
- Ferreira LA, Alves EW, Henriques OB: Peptide T, a novel bradykinin potentiator isolated from Tityus serrulatus scorpion venom. Toxicon. 1993, 31: 941-947. 10.1016/0041-0101(93)90253-F.PubMedView ArticleGoogle Scholar
- Meki AR, Nassar AY, Rochat H: A bradykinin-potentiating peptide (peptide K12) isolated from the venom of Egyptian scorpion Buthus occitanus. Peptides. 1995, 16: 1359-1365. 10.1016/0196-9781(95)02036-5.PubMedView ArticleGoogle Scholar
- Zeng XC, Li WX, Peng F, Zhu ZH: Cloning and characterization of a novel cDNA sequence encoding the precursor of a novel venom peptide (BmKbpp) related to a bradykinin-potentiating peptide from Chinese scorpion Buthus martensii Karsch. IUBMB Life. 2000, 49: 207-210. 10.1080/152165400306214.PubMedView ArticleGoogle Scholar
- Valdez-Cruz NA, Segovia L, Corona M, Possani LD: Sequence analysis and phylogenetic relationships of gene encoding heterodimeric Phospholipases A2 from the venom of the scorpion Anuroctonus phaiodactylus. Gene. 2007, 396: 149-158. 10.1016/j.gene.2007.03.007.PubMedView ArticleGoogle Scholar
- Harvey AL, Robertson B: Dendrotoxins: structure-activity relationships and effects on potassium ion channels. Curr Med Chem. 2004, 11: 3065-3072.PubMedView ArticleGoogle Scholar
- Bruhn HD, Beress L: Polyvalent proteinase inhibitor from the sea anemone Anemonia sulcata: effect on coagulation and fibrinolysis. Thromb Haemost. 1978, 39: 552-554.PubMedGoogle Scholar
- Bayrhuber M, Vijayan V, Ferber M, Graf R, Korukottu J, Imperial J, Garrett JE, Olivera BM, Terlau H, Zweckstetter M, Becker S: Conkunitzin-S1 is the first member of a new Kunitz-type neurotoxin family. Structural and functional characterization. J Biol Chem. 2005, 280: 23766-23770. 10.1074/jbc.C500064200.PubMedView ArticleGoogle Scholar
- Peng K, Lin Y, Liang SP: Nuclear magnetic resonance studies on huwentoxin-XI from the Chinese bird spider Ornithoctonus huwena: 15N labeling and sequence-specific 1H, 15N nuclear magnetic resonance assignments. Acta Biochim Biophys Sin (Shanghai). 2006, 38: 457-466. 10.1111/j.1745-7270.2006.00191.x.View ArticleGoogle Scholar
- Lazo GR, Tong J, Miller R, Hsia C, Rausch CY, Kang Y, Anderson OD: Software scripts for quality checking of high throughput nucleic acid sequencers. Biotechniques. 2001, 30: 1300-1305.PubMedGoogle Scholar
- Huang X, Madan A: CAP3: a DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMed CentralPubMedView ArticleGoogle Scholar
- Laboratório de Biologia Molecular – Universidade de Brasilia. [https://www.biomol.unb.br/]
- NCBI BLAST. [http://www.ncbi.nlm.nih.gov/blast]
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.PubMedView ArticleGoogle Scholar
- Notredame C, Higgins DG, Heringa J: T-Coffee: a novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000, 302: 205-217. 10.1006/jmbi.2000.4042.PubMedView ArticleGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralPubMedView ArticleGoogle Scholar
- BioEdit. [http://www.mbio.ncsu.edu/BioEdit/bioedit.html]
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.