Assembling an arsenal, the scorpion way
© Kozminsky-Atias et al; licensee BioMed Central Ltd. 2008
Received: 28 August 2008
Accepted: 16 December 2008
Published: 16 December 2008
For survival, scorpions depend on a wide array of short neurotoxic polypeptides. The venoms of scorpions from the most studied group, the Buthida, are a rich source of small, 23–78 amino acid-long peptides, well packed by either three or four disulfide bridges that affect ion channel function in excitable and non-excitable cells.
In this work, by constructing a toxin transcripts data set from the venom gland of the scorpion Buthus occitanus israelis, we were able to follow the evolutionary path leading to mature toxin diversification and suggest a mechanism for leader peptide hyper-conservation. Toxins from each family were more closely related to one another than to toxins from other species, implying that fixation of duplicated genes followed speciation, suggesting early gene conversion events. Upon fixation, the mature toxin-coding domain was subjected to diversifying selection resulting in a significantly higher substitution rate that can be explained solely by diversifying selection. In contrast to the mature peptide, the leader peptide sequence was hyper-conserved and characterized by an atypical sub-neutral synonymous substitution rate. We interpret this as resulting from purifying selection acting on both the peptide and, as reported here for the first time, the DNA sequence, to create a toxin family-specific codon bias.
We thus propose that scorpion toxin genes were shaped by selective forces acting at three levels, namely (1) diversifying the mature toxin, (2) conserving the leader peptide amino acid sequence and intriguingly, (3) conserving the leader DNA sequences.
The order Scorpiones constitutes one of the most ancient groups of animals on earth (more than 400 million years of evolution), represented today by approximately 1500 different species [1, 2] which can be divided into four groups. The venoms of scorpions from the most studied group, the Buthida, are a rich source of small, 23–78 amino acid-long peptides, well packed by either three or four disulfide bridges, that affect ion channel function in excitable and non-excitable cells [1, 2]. Typically, toxin peptides contain a structural signature defined by the presence of a cysteine-stabilized α/β motif, in which two disulfide bridges covalently link a α-helical segment with one strand of a β-sheet structure [3–5]. The dense core typically contains the sequence motifs, CXXXC and CXC (where C stands for cysteine and X stands for any amino acid residue), shown to be essentially conserved [1, 2, 6]. Toxins are mainly neurotoxic peptides, interacting specifically with various ionic channels (i.e. Na+, K+, Cl- and Ca2+ channels). The Na+ channel toxins are gating modifiers that can be divided into two families, according to their pharmacological effect: α- and β-toxins. α-toxins bind in a voltage dependent mode and slow channel inactivation [7, 8]. β-toxins bind to a separate site, independently of membrane potential, and shift the channel activation to the hyperpolarizing direction [7, 8]. β-toxins can be further divided, according to their amino acid sequence and selectivity, into "true" β, excitatory and depressant toxins. As opposed to all others, the excitatory and depressant toxins are not toxic to mammals. Short chain toxins form a large family of peptides that block K+ or Cl- channels [5, 9].
Within each toxin family, the amino acid sequence of the leader peptide is conserved, while the mature toxin-coding domain displays hyper-diversification, especially in those residues exposed to the surface of the molecule and which affect toxin activity [2, 10]. These characteristics (i.e. conserved cysteine scaffold, conserved leader peptide and hyper-variable mature toxin) are shared with toxic peptides from other animal orders, such as snakes and cone snails [6, 11, 12]. In all these animals, evolution shows an apparent preference for toxin-diversifying processes [11, 13].
Scorpions feed on insects, spiders and other small invertebrates. Consequently, scorpion venom contains toxins which effectively paralyze this type of prey. The venom, however, also serves a defensive function, and is, therefore, highly lethal to vertebrates, including mammals. This has been suggested as the driving force for the apparently accelerated rate of evolution displayed by mature toxins . Although a number of studies have considered conopeptides and snake toxins from an evolutionary perspective [11, 13], thus far, the mechanisms leading to accelerated evolution of scorpion toxins have not been thoroughly addressed.
Accordingly, to study the evolution of scorpion toxin gene families, in the present report we have unraveled the venom gland content of the Israeli scorpion, Buthus occitanus israelis (Boi). In all identified toxin families, the leader peptide sequences displayed remarkable conservation, owing to purifying selection forces aimed at conserving both the amino acid and nucleotide sequences, while the mature toxin-coding sequences have undergone diversifying selection. These observations provided the unique opportunity to study the evolutionary fate of gene families simultaneously undergoing opposite selective pressures in separate domains. The data obtained provide insight into the evolutionary pathways leading to the diversification of the mature toxin and imply the existence of a novel mechanism leading to family-specific leader peptide hyper-conservation.
The venom gland contains a vast variety of neurotoxins
Striking differences in the transcription levels of different toxins could also be observed (Fig. 1C). Order-of-magnitude differences were observed between the relative expression levels of different transcripts from the same family, such as α-toxins, excitatory toxins as well as Cl-- and K+-channel toxins. Seventy five percent of the transcripts had relatively low transcript representation (≤ 2 independently isolated clones). Whereas, on the other hand, 5% of the transcripts had relatively high transcript representation (≥ 10 clones). Hence, related scorpion toxins may be expressed at vastly different levels in the venom gland.
Are similar cDNA transcripts distinct genes or allelic variants?
Since the cDNA library was constructed from the pooled RNA from 10 different individuals, it was imperative to assess whether the observed transcript variety constituted the presence of multiple genes or was an outcome of the sample containing a much smaller number of multi-allele genes. To this end, all homologous transcript pairs with a percent identity of greater than 90% were studied. These correspond to four homologous cDNA pairs, namely: BoiTx1 and BoiTx260, putative K+ channel toxins exhibiting 97% identity; BoiTx492 and BoiTx475, putative depressant toxins exhibiting 92.8% identity; BoiTx611 and BoiTx357, putative Cl- channel toxins exhibiting 99% identity and BoiTx651 and BoiTx458, putative β-like toxins exhibiting 94.4% identity. The homologous pairs were subjected to PCR amplification using primers designed to identify both tested sequences. The templates for these experiments were the genomic DNA of all 10 scorpions from whence the cDNA library was generated. Individual cDNA clones served as control. In addition, the homologous cDNA pairs were digested with restriction enzymes directed against a recognition site found on only one of the clones in each pair (i.e. NcoI, HindIII, and LweI for pairs 1, 2 and 3, 4, respectively). Analysis of the amplification and restriction digestion products on agarose gel revealed that gene pairs 2, 3 and 4 exhibited digestion of only one partner of the pair in all 10 scorpions. On the other hand, clone pair 1 (97% identity) exhibited a single digestion pattern (i.e. heterozygous) in some individuals or no digestion at all in others. This digestion pattern implies that the BoiTx1 and BoiTx260 transcripts (i.e. pair 1) represent allelic variants of the same transcript, while all other sequences most likely correspond to separate genes.
Toxin divergence occurred after speciation
Leader peptide conservation and mature toxin hyper-diversification
To examine whether other toxin families (e.g. α-, β-like and Cl- channel toxins) possess the same pattern of conservation versus hyper-diversification, an unrooted neighbor-joining tree diagram of the different domains was plotted (Fig. 3B). While segregation of the different gene families is clearly observed in the leader domain (Fig. 3B, left plot), most transcripts diverge directly from the origin of the mature domain tree (Fig. 3B, right plot). The tight cluster of the leader domain tree, versus the star-like shape of the mature domain tree, demonstrates that the same pattern of conserved leader and diverse mature peptide sequence exists in all of the distinctive toxin families.
These data thus demonstrate that scorpion toxins gene families can be defined primarily on the basis of their highly conserved leader domains, in combination with distinctly conserved cysteine patterns, as was previously suggested for other scorpion species .
The toxin gene domains experience opposite selective patterns
Apparent position-specific codon conservation in strong negatively-selected sites
Nucleotide diversity is greater within the mature toxin-coding domain than in other gene domains
Can this difference in substitution rates between the leader and the mature domains be attributed solely to the influence of selective forces? To assess this possibility, the nucleotide sequences of all Boi depressant toxin introns were isolated and determined. One of the depressant toxins genes lacked an intron. All other genes possess a 298 to 327 bp-long intron, located near the 3' end of the leader-coding sequence, at a position similar to other scorpion toxins genes [26, 27]. This intron domain is thought to reflect a neutral nucleotide mutation-fixation rate. Indeed, Boi toxin intron nucleotide diversity exceeds that of the leader-coding domain but is much lower than that of the mature toxin coding domain (Fig. 6A, B). To assess the neutral substitution rates affecting the coding domains while eliminating the effects of protein-shaping selective forces, the substitution rates per site in the wobble (3rd codon) position of the negatively-selected sites were determined (Fig. 6A). The neutral substitution rate within the mature toxin coding region was similar to that of the intron. On the other hand, the neutral substitution rate within the leader-coding region was significantly (~2-fold) lower (Fig. 6A). This can be interpreted as an indication of a decrease in the synonymous substitution rate and of position-specific codon bias in the leader-coding domain.
Transition preference in the leader toxin coding domain
We further characterized Boi's toxin genes dataset by examining the transversion/transition (Tv/Ts) ratios in depressant and α-toxins. We found a two-fold excess of transitions in the leader sequence in comparison with that of the mature-coding domain (Fig. 6C). To ascertain the neutrally selected Tv/Ts ratio, we examine this ratio within the intron region. The Tv/Ts ratio of the mature-coding domain was very similar to that of the intron while the leader was found to be under a strong transition bias.
Scorpion venoms have evolved during the last 400 million years to constitute an arsenal of high affinity ion channel modulators. Here, we investigated evolutionary pathways leading to the establishment of this pharmacological wealth by constructing and analyzing a Boi toxin cDNA library. The cDNA library transcripts exemplify the tremendous variety of toxin genes prevalent in an old world scorpion. The apparent hyper-diversity within each toxin family could be explained, in part, by the presence of multiple alleles per gene. Our comparison of genomic DNA versus cDNA sequences in several scorpion individuals indicates that in all but one of the cDNA clone pairs, diversity is not likely to reflect allelic variants but rather the existence of distinct genes. The degree of diversification within the different toxin families was versatile. Boi putative K+ channel blockers are much more diverse than Na+ channel modifiers, within the same activity sub-type (i.e. depressant, excitatory, α- and β-toxins) (Fig. 1A). We thus speculate that the high divergence of K+ channel toxins is due to the higher variance in possible targets (e.g. mammals have ~10-fold more K+ channel sub-types than Na+ channels).
Gene duplications upon speciation
The number of toxin genes in a single scorpion species indicates that several gene duplications occurred during scorpion evolution. Gene duplication is typically a consequence of unequal crossing-over, retro-positioning or chromosomal duplication . In another species, the sea anemone Nematostella vectensis, it was recently suggested that toxin gene duplication was achieved by unequal crossing-over . Although most scorpion toxin genes could have arisen via gene duplication, we have identified here, for the first time, a Boi depressant gene lacking the conserved intron located within the leader sequence coding domain. This molecular feature is indicative of the retro-positioning which occurs when an mRNA is retro-transcribed to cDNA and then inserted into the genome.
Following examination of toxin diversity in several venomous animals (snakes and cone snails [11, 17, 18, 20]) it was suggested that duplication of the toxin genes followed speciation. Our inter-species examination of depressant and α-toxin families in Boi and BmK scorpions revealed that most of the toxin subfamilies are unique to each species, suggesting again that speciation preceded gene duplication.
Here we suggest a possible evolutionary mechanism enabling fixation of toxin gene duplicates upon speciation. Similar to point mutations, duplications occur in an individual, which can be either fixed or lost in the population. If the new allele is selectively neutral, as compared with pre-existing alleles, that allele has only the small probability of 1/2 N of being fixed in a diploid population , where N is the effective population size. This suggests that the vast majority of the duplicated genes will be lost. The speciation process entails a vast decrease in the new effective population size, thus opening a short window of time enabling an increase in the fixation rate of gene duplication. Therefore, the association of duplication with speciation could be due to a rise in the duplication fixation rate upon speciation, caused by a decrease in the effective population size. For those duplicated genes that do become fixed, fixation is time-consuming. On average, it takes 4 N generations for a neutral allele to become fixed in a population of N individuals . Furthermore, the number of toxins in each family is similar in the two scorpions, as is the overall number of toxins in the venom. The latter was shown to also be the case in other scorpions , as well as with other species [11, 32]. Such observations imply that upon a dramatic change in the environment that promotes speciation, rapid fixation rate of duplicated toxin genes occurs so as to allow the introduction of novel toxic phenotypes. At the same time, a similar number of toxin genes are lost. This may allow for adaptation to novel ecological niches and pray, while maintaining a minimal concentration of each of the toxins within the limited venom pool, as was recently suggested for the PLA2 genes of snakes .
Could gene conversion cause the intra-species toxin similarity? In the sea anemone N. vectensis , as well as in two species of sea snakes [33, 34], a very high degree of conservation has been reported within certain toxin families, suggesting the involvement of concerted evolution processes . Here, the degree of diversity among paralogous genes roles out the possibility of recent gene conversion, although the higher degree of similarity between paralogous genes, as compared with orthologous genes, can not rule out the occurrence of early gene conversion events .
Mature toxin gene diversification
Protein neofunctionalism after gene duplication is believed to have played a major role in evolution , although the mechanisms by which this process occurs remains controversial. When comparing the amino acid sequences of protoxins from all families, hyper-diversification of the mature toxin and conservation of the leader peptide were observed (Fig. 3). In this study, we examined the Boi depressant toxin gene family to evaluate the evolutionary processes leading to this phenomenon. Within the mature-coding region, diversifying forces (Fig. 4) are sufficient to explain the high observed substitution rate (Fig. 6A, B) given that the neutral substitution rate within the synonymous sites (i.e. substitution rate at the wobble position at negatively-selected sites) is the same as that of the mostly neutral intron (Fig. 6A). The high substitution rate of the positively-selected sites gave rise to alleged position-specific codon conservation at adjacent negatively-selected sites (Fig. 5).
Leader toxin hyper-conservation
The strong inter-species leader conservation, as indicated by their defined segregation between Boi and BmK scorpions within depressant (Fig. 2B) and α-toxin families (not shown), supports our suggestion for a codon usage-dependant translational regulatory mechanism. This implies that a transcriptional regulatory mechanism may be at work that differently affects the toxin families tested and is species-specific (Fig. 2B), as expected. Although, we can not rule out other mechanisms for non-neutral evolution at synonymous sites, such as regulation of mRNA stability and/or intron splicing , here we propose that the expression levels of toxin families are differentially controlled via regulation at the translation level. Our suggested mechanism exploits the relative abundance of different tRNA species to either restrict or enhance translation rates.
In addition, a clear bias for transitions over transversions was observed in the leader-coding domain in comparison to the Tv/Ts ratio within the intron and the mature toxin-coding domain (Fig. 6C). This could be attributed to strong negative selection forces inflicted upon this region, as at the most variable site, i.e., the wobble position, a transition would be preferred since in 94% of cases, this would result in a silent synonymous mutation.
In summary, positive selection upon mature toxin residues was reported in the venoms of other venomous organisms, such as snakes [12, 13, 39], cone snails [11, 17] and spiders , suggesting a common diversifying evolutionary mechanism. Indeed, hyper-variability in gene families is apparent in systems evolved to recognize foreign molecules, whether those genes encode venom-derived toxins, gene families of the immune system or antigenic parasite surface proteins . In cone snails, a hyper-variability-generating molecular mechanism relying upon an error prone-like DNA polymerase was suggested .
In this work, by constructing and analyzing a dataset of toxin transcripts from the venom gland of the scorpion Boi, we were able to follow the evolutionary path leading to scorpion toxin diversification. Duplicated genes encoding toxin peptides were fixed within the genome mainly following speciation. Upon such fixation, the mature toxin coding domain was subjected to diversifying selection to increase species fitness. This resulted in alleged position-specific codon conservation at adjacent negatively-selected sites. On the other hand, the leader peptide is characterized by a conservation of the amino acid sequence, a sub-neutral synonymous substitution rate and a strong transition bias. This was explained by purifying selection forces acting to conserve both the peptide and DNA sequences. This suggests the presence of a family-specific, codon-sensitive translational regulatory mechanism for toxin genes.
cDNA library construction and screening
Buthus occitanus israelis (Boi) scorpions were collected at Sde Boker, Israel. Total RNA was extracted from the venom glands of 10 scorpions using an EZ RNA extraction kit (Promega), 2 days after electrical 'milking' of the venom. mRNA was purified with a PolyATract mRNA Isolation System (Promega). Double-stranded cDNA was synthesized from 2 μg of mRNA using the Universal Riboclone cDNA Synthesis System (Promega). Blunt end cDNA clones were cloned into the pBluescript KS+ (pBS) plasmid digested with SmaI and transfected into Escherichia coli DH5α cells. Electro-transformation of the cDNA ligation yielded a library comprising approximately 5 × 104 primary clones. To determine insert size, individual clones were amplified by PCR using T7 and T3 primers. To avoid sequencing through the poly A tail, insert orientation was determined by PCR using the T3, T7 and poly T primers. Four hundred and twenty randomly chosen cDNA clones, ranging in length from 250 to 600 bp, were then sequenced to obtain a reliable representation of the toxin content in the venom gland.
Plasmid DNA was purified using a Wizard plus SV Miniprep kit (Promega). Restriction enzyme digestions, DNA ligations and phosphorylation and dephosphorylation reactions were performed according to manufacturer instructions (Fermentas or New England Biolabs). Competent E. coli DH5α cells were transformed by either heat shock or electroporation procedures. Polymerase chain reactions (PCR) were performed with a PTC-2000 apparatus (MJ Research), using Pfu (Fermentas) or Taq (Promega) enzymes, as indicated.
Differentiation between genes and alleles
Genomic DNA was isolated from the same 10 scorpions used for the cDNA library construction. Four pairs of homologous genes exhibiting over 90% identity were selected for analysis. As such, four PCR primer pairs were designed, according to the cDNA sequence, to amplify each gene pair. To distinguish between the clones, digestion by a restriction enzyme direct against a recognition site found on only one of the genes in each pair was preformed.
Sequencing the introns of depressant toxins genes
Introns from the depressant toxin family were amplified using Pfu DNA polymerase with Boi genomic DNA as template. PCR forward primers were constructed to match a highly conserved sequence of the 5' UTR and leader regions. Reverse primers were designed to match a hyper-variable region within the mature toxin-coding domain to ensure specificity. PCR amplification products were cloned into a SmaI-digested pBS vector. Positive clones were then sequenced to retrieve the intron sequence.
Toxin classification and homology identification was achieved using the BLASTN and BLASTP programs . Individual transcripts were aligned using MUSCLE algorithm . Alignments were refined manually by the Jalview program . Unrooted phylogenetic trees were constructed using the neighbor-joining method and maximum likelihood algorithms  using MEGA4 software  and then visualized with TreeView . Synonymous versus non-synonymous substitution rates were analyzed using the DNAsp  and Selecton  programs. Transitions versus transversion rates, as well as codon usage profiles, were drawn using MEGA4 software. Nucleotide diversity rates were estimated by the DNAsp program.
We are grateful to Moran Gershoni for helpful discussions, Ofer Ovadia for help with statistical analyses and Jerry Eichler for advice on the manuscript. This work was supported by grants from the Israel Science Foundation (431/03) and the Zlotowski Center for Neuroscience to N.Z.
- Goyffon M, Landon C: Scorpion toxins and defensins. Comptes rendus des seances de la Societe de biologie et de ses filiales. 1998, 192 (3): 445-462.PubMedGoogle Scholar
- Rochat H, Bernard P, Couraud F: Scorpion toxins: chemistry and mode of action. Advances in cytopharmacology. 1979, 3: 325-334.PubMedGoogle Scholar
- Garcia ML, Hanner M, Kaczorowski GJ: Scorpion toxins: tools for studying K+ channels. Toxicon. 1998, 36 (11): 1641-1650.View ArticlePubMedGoogle Scholar
- Garcia ML, Hanner M, Knaus HG, Slaughter R, Kaczorowski GJ: Scorpion toxins as tools for studying potassium channels. Methods Enzymol. 1999, 294: 624-639.View ArticlePubMedGoogle Scholar
- Rodriguez de la Vega RC, Possani LD: Current views on scorpion toxins specific for K+-channels. Toxicon. 2004, 43 (8): 865-875.View ArticlePubMedGoogle Scholar
- Possani LD, Merino E, Corona M, Bolivar F, Becerril B: Peptides and genes coding for scorpion toxins that affect ion-channels. Biochimie. 2000, 82 (9–10): 861-868.View ArticlePubMedGoogle Scholar
- Jover E, Couraud F, Rochat H: Two types of scorpion neurotoxins characterized by their binding to two separate receptor sites on rat brain synaptosomes. Biochemical and biophysical research communications. 1980, 95 (4): 1607-1614.View ArticlePubMedGoogle Scholar
- Jover E, Martin-Moutot N, Couraud F, Rochat H: Binding of scorpion toxins to rat brain synaptosomal fraction. Effects of membrane potential, ions, and other neurotoxins. Biochemistry. 1980, 19 (3): 463-467.View ArticlePubMedGoogle Scholar
- DeBin JA, Maggio JE, Strichartz GR: Purification and characterization of chlorotoxin, a chloride channel ligand from the venom of the scorpion. Am J Physiol. 1993, 264 (2 Pt 1): C361-369.PubMedGoogle Scholar
- Mebs D: Toxicity in animals. Trends in evolution?. Toxicon. 2001, 39 (1): 87-96.View ArticlePubMedGoogle Scholar
- Conticello SG, Gilad Y, Avidan N, Ben-Asher E, Levy Z, Fainzilber M: Mechanisms for evolving hypervariability: the case of conopeptides. Mol Biol Evol. 2001, 18 (2): 120-131.View ArticlePubMedGoogle Scholar
- Fujimi TJ, Nakajyo T, Nishimura E, Ogura E, Tsuchiya T, Tamiya T: Molecular evolution and diversification of snake toxin genes, revealed by analysis of intron sequences. Gene. 2003, 313: 111-118.View ArticlePubMedGoogle Scholar
- Ohno M, Menez R, Ogawa T, Danse JM, Shimohigashi Y, Fromen C, Ducancel F, Zinn-Justin S, Le Du MH, Boulain JC, et al: Molecular evolution of snake toxins: is the functional diversity of snake toxins associated with a mechanism of accelerated evolution?. Prog Nucleic Acid Res Mol Biol. 1998, 59: 307-364.View ArticlePubMedGoogle Scholar
- Zhijian C, Feng L, Yingliang W, Xin M, Wenxin L: Genetic mechanisms of scorpion venom peptide diversification. Toxicon. 2006, 47 (3): 348-355.View ArticlePubMedGoogle Scholar
- Rodriguez de la Vega RC, Possani LD: Overview of scorpion toxins specific for Na+ channels and related peptides: biodiversity, structure-function relationships and evolution. Toxicon. 2005, 46 (8): 831-844.View ArticlePubMedGoogle Scholar
- Goudet C, Chi CW, Tytgat J: An overview of toxins and genes from the venom of the Asian scorpion Buthus martensi Karsch. Toxicon. 2002, 40 (9): 1239-1258.View ArticlePubMedGoogle Scholar
- Duda TF, Palumbi SR: Molecular genetics of ecological diversification: duplication and rapid evolution of toxin genes of the venomous gastropod Conus. Proc Natl Acad Sci USA. 1999, 96 (12): 6820-6823.PubMed CentralView ArticlePubMedGoogle Scholar
- Olivera BM, Walker C, Cartier GE, Hooper D, Santos AD, Schoenfeld R, Shetty R, Watkins M, Bandyopadhyay P, Hillyard DR: Speciation of cone snails and interspecific hyperdivergence of their venom peptides. Potential evolutionary significance of introns. Ann N Y Acad Sci. 1999, 870: 223-237.View ArticlePubMedGoogle Scholar
- Zhu S, Bosmans F, Tytgat J: Adaptive evolution of scorpion sodium channel toxins. J Mol Evol. 2004, 58 (2): 145-153.View ArticlePubMedGoogle Scholar
- Lynch VJ: Inventing an arsenal: adaptive evolution and neofunctionalization of snake venom phospholipase A2 genes. BMC Evol Biol. 2007, 7 (2): 2-PubMed CentralView ArticlePubMedGoogle Scholar
- Tian C, Yuan Y, Zhu S: Positively selected sites of scorpion depressant toxins: possible roles in toxin functional divergence. Toxicon. 2008, 51 (4): 555-562.View ArticlePubMedGoogle Scholar
- Kini RM, Chan YM: Accelerated evolution and molecular surface of venom phospholipase A2 enzymes. J Mol Evol. 1999, 48 (2): 125-132.View ArticlePubMedGoogle Scholar
- Gibbs HL, Rossiter W: Rapid Evolution by Positive Selection and Gene Gain and Loss: PLA2 Venom Genes in Closely Related Sistrurus Rattlesnakes with Divergent Diets. J Mol Evol. 2008, 66 (2): 151-166.View ArticlePubMedGoogle Scholar
- Stern A, Doron-Faigenboim A, Erez E, Martz E, Bacharach E, Pupko T: Selecton 2007: advanced models for detecting positive and purifying selection using a Bayesian inference approach. Nucleic Acids Res. 2007, W506-511. 35 Web ServerGoogle Scholar
- Karon JM: The covarion model for the evolution of proteins: parameter estimates and comparison with Holmquist, Cantor, and Jukes' stochastic model. J Mol Evol. 1979, 12 (3): 197-218.View ArticlePubMedGoogle Scholar
- Kozminsky-Atias A, Somech E, Zilberberg N: Isolation of the first toxin from the scorpion Buthus occitanus israelis showing preference for Shaker potassium channels. FEBS Lett. 2007, 581 (13): 2478-2484.View ArticlePubMedGoogle Scholar
- Froy O, Sagiv T, Poreh M, Urbach D, Zilberberg N, Gurevitz M: Dynamic diversification from a putative common ancestor of scorpion toxins affecting sodium, potassium, and chloride channels. J Mol Evol. 1999, 48 (2): 187-196.View ArticlePubMedGoogle Scholar
- Zhang J: Evolution by gene duplication: an update. Trends Ecol Evol. 2003, 18 (6): 292-298.View ArticleGoogle Scholar
- Moran Y, Weinberger H, Sullivan JC, Reitzel AM, Finnerty JR, Gurevitz M: Concerted evolution of sea anemone neurotoxin genes is revealed through analysis of the Nematostella vectensis genome. Mol Biol Evol. 2008, 25 (4): 737-747.View ArticlePubMedGoogle Scholar
- Takahata N: Neutral theory of molecular evolution. Curr Opin Genet Dev. 1996, 6 (6): 767-772.View ArticlePubMedGoogle Scholar
- Possani LD, Becerril B, Delepierre M, Tytgat J: Scorpion toxins specific for Na+-channels. Eur J Biochem. 1999, 264 (2): 287-300.View ArticlePubMedGoogle Scholar
- Sollod BL, Wilson D, Zhaxybayeva O, Gogarten JP, Drinkwater R, King GF: Were arachnids the first to use combinatorial peptide libraries?. Peptides. 2005, 26 (1): 131-139.View ArticlePubMedGoogle Scholar
- Li M, Fry BG, Kini RM: Putting the brakes on snake venom evolution: the unique molecular evolutionary patterns of Aipysurus eydouxii (Marbled sea snake) phospholipase A2 toxins. Mol Biol Evol. 2005, 22 (4): 934-941.View ArticlePubMedGoogle Scholar
- Tamiya T, Fujimi TJ: Molecular evolution of toxin genes in Elapidae snakes. Mol Divers. 2006, 10 (4): 529-543.View ArticlePubMedGoogle Scholar
- Nei M, Rooney AP: Concerted and birth-and-death evolution of multigene families. Annu Rev Genet. 2005, 39: 121-152.PubMed CentralView ArticlePubMedGoogle Scholar
- Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. Science. 2000, 290 (5494): 1151-1155.View ArticlePubMedGoogle Scholar
- Claros MG, Brunak S, von Heijne G: Prediction of N-terminal protein sorting signals. Curr Opin Struct Biol. 1997, 7 (3): 394-398.View ArticlePubMedGoogle Scholar
- Chamary JV, Parmley JL, Hurst LD: Hearing silence: non-neutral evolution at synonymous sites in mammals. Nat Rev Genet. 2006, 7 (2): 98-108.View ArticlePubMedGoogle Scholar
- Hseu TH, Jou ED, Wang C, Yang CC: Molecular evolution of snake venom toxins. J Mol Evol. 1977, 10 (2): 167-182.View ArticlePubMedGoogle Scholar
- Field MC, Boothroyd JC: Sequence divergence in a family of variant surface glycoprotein genes from trypanosomes: coding region hypervariability and downstream recombinogenic repeats. J Mol Evol. 1996, 42 (5): 500-511.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402.PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797.PubMed CentralView ArticlePubMedGoogle Scholar
- Clamp M, Cuff J, Searle SM, Barton GJ: The Jalview Java alignment editor. Bioinformatics. 2004, 20 (3): 426-427.View ArticlePubMedGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.PubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599.View ArticlePubMedGoogle Scholar
- Page RD: TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12 (4): 357-358.PubMedGoogle Scholar
- Rozas J, Rozas R: DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis. Bioinformatics. 1999, 15 (2): 174-175.View ArticlePubMedGoogle Scholar
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.