Diversification of a single ancestral gene into a successful toxin superfamily in highly venomous Australian funnel-web spiders
- Sandy S Pineda†1,
- Brianna L Sollod†2, 7,
- David Wilson†1, 3, 8,
- Aaron Darling1, 9,
- Kartik Sunagar4, 5,
- Eivind A B Undheim1, 6,
- Laurence Kely6,
- Agostinho Antunes4, 5,
- Bryan G Fry1, 6Email author and
- Glenn F King1Email author
© Pineda et al.; licensee BioMed Central Ltd. 2014
Received: 3 December 2013
Accepted: 26 February 2014
Published: 5 March 2014
Spiders have evolved pharmacologically complex venoms that serve to rapidly subdue prey and deter predators. The major toxic factors in most spider venoms are small, disulfide-rich peptides. While there is abundant evidence that snake venoms evolved by recruitment of genes encoding normal body proteins followed by extensive gene duplication accompanied by explosive structural and functional diversification, the evolutionary trajectory of spider-venom peptides is less clear.
Here we present evidence of a spider-toxin superfamily encoding a high degree of sequence and functional diversity that has evolved via accelerated duplication and diversification of a single ancestral gene. The peptides within this toxin superfamily are translated as prepropeptides that are posttranslationally processed to yield the mature toxin. The N-terminal signal sequence, as well as the protease recognition site at the junction of the propeptide and mature toxin are conserved, whereas the remainder of the propeptide and mature toxin sequences are variable. All toxin transcripts within this superfamily exhibit a striking cysteine codon bias. We show that different pharmacological classes of toxins within this peptide superfamily evolved under different evolutionary selection pressures.
Overall, this study reinforces the hypothesis that spiders use a combinatorial peptide library strategy to evolve a complex cocktail of peptide toxins that target neuronal receptors and ion channels in prey and predators. We show that the ω-hexatoxins that target insect voltage-gated calcium channels evolved under the influence of positive Darwinian selection in an episodic fashion, whereas the κ-hexatoxins that target insect calcium-activated potassium channels appear to be under negative selection. A majority of the diversifying sites in the ω-hexatoxins are concentrated on the molecular surface of the toxins, thereby facilitating neofunctionalisation leading to new toxin pharmacology.
KeywordsSpider toxin Spider venom Hexatoxin ω-hexatoxin κ-hexatoxin Australian funnel-web spider Molecular evolution Gene duplication Positive selection Negative selection
Venoms have proven to be key evolutionary innovations for many divergent animal lineages [1, 2]. Although the most extensively studied venoms are from the medically important scorpions, snakes, and spiders, venom systems are present in many other lineages including cnidarians, echinoderms, molluscs, fish, lizards, and mammals [1, 2]. These venoms have evolved to serve a variety of purposes, including prey capture, competitor deterrence, and defense against predators. There has been considerable innovation both in the chemical composition of these venoms as well as the method of venom delivery, which includes barbs, beaks, fangs, harpoons, nematocysts, pinchers, proboscises, spines, spurs, and stingers [1, 2].
From a molecular evolutionary perspective, the venoms of snakes are the best understood. There is now abundant evidence that snake venoms evolved by recruitment of genes encoding normal body proteins followed by extensive duplication, neofunctionalization, and in some instances relegation to the status of pseudogene [1, 3–5]. In many cases, these genes have been explosively replicated to produce large multigene families. This process is analogous to the birth-and-death model of evolution proposed for multigene families involved in adaptive immunity, such as the major histocompatibility complex and immunoglobulin VH genes . However, the evolutionary trajectory is less clear for the venoms of spiders, scorpions, and molluscs, which are dominated by disulfide-rich peptides of mass 2–9 kDa [7–12]. These peptides typically possess high affinity and often-exquisite specificity for particular classes of ion channels and other nervous system targets [13–15]. These neurotoxic functions are perhaps not surprising given that the primary role of these venoms is to paralyse or kill envenomated prey [11, 16, 17].
Results and discussion
The ω-hexatoxins are expressed as prepropeptide precursors
All of the ω-HXTXs are expressed as prepropeptide precursors that are posttranslationally processed to yield the mature toxin sequence (Figure 2A). The highly hydrophobic 22-residue signal sequence is of similar length to that reported for peptide-toxin precursors from spiders, scorpions, and cone snails . The 20-residue propeptide sequence is highly acidic, with a net charge of -4, a feature that has been noted for numerous spider-toxin propeptide sequences [24–28] but which is not characteristic of toxin precursors from other venomous animals. Moreover, the presence of a propeptide region contrasts with most scorpion-toxin precursors in which the signal sequence is fused directly to the mature toxin without an intervening propeptide . The reason for the highly acidic propeptide region in spider toxin precursors remains to be determined, but it may be related to specific interactions between the toxin precursor and components of the secretory and/or protein folding pathway in spider venom glands.
The propeptide sequence terminates with a dibasic Arg-Arg signature; dibasic sequences are common recognition sites for proteolytic removal of propeptide segments in neuropeptide precursors from both vertebrates  and invertebrates . While Arg is the terminal residue in virtually all known spider-toxin propeptide sequences, the penultimate residue is variable, though it is commonly Asp, Glu, or Lys [25, 26, 31–33].
ω-hexatoxins belong to a large toxin-gene superfamily
Orthologs of ω-HXTX-Hv1a were identified in all five species of Australian funnel-web spider examined in this study. These species are distributed along the eastern seaboard of Australia with a geographic spread of more than 2000 km (Figure 1). The hexathelids are a group of approximately 40 species divided into three genera: Atrax, Hadronyche and Illawarra[34–36]. They are adapted to forest environments but can also be found in habitats that range from montane herblands and open woodland to closed forest . Conservation of the ω-hexatoxin family of toxins over this wide range of environments and differing prey distributions implies that there has been strong evolutionary pressure to maintain these peptides as part of the venom arsenal, which is perhaps not surprising given that they are broadly active against many different arthropods [11, 17, 37].
In addition to obvious homologs of ω-HXTX-Hv1a, the RACE and transcriptomic analyses revealed additional families of toxins that had almost identical signal sequences to the ω-HXTX-1 transcripts, but divergent propeptide and mature toxin sequences. We named one of these families the ω/κ-HXTX family (described as U-ACTX in ). The ω/κ-HXTX peptides appear to be distributed in two of the species examined (A. robustus and H. versuta); this reinforces the idea that these toxins most likely arose ancestrally by duplication of a ω-HXTX-1 gene followed by hypermutation of the propeptide and mature-toxin regions in order to create a new function (neofunctionalization). The conservation and radiation of these toxins across this family of spiders implies that they are not nonfunctional relics of an explosive radiation of this toxin-gene superfamily, and we confirmed this by showing that recombinant ω/κ-HXTX-Hv1a is highly insecticidal . The high insecticidal potency of this family of peptides is believed to result from a synergistic effect on insect voltage-gated calcium (CaV) channels and calcium-activated potassium (KCa) channels .
RACE analysis of the venom-gland cDNA library from H. versuta also led to amplification of transcripts encoding the insecticidal toxin κ-HXTX-Hv1c , and sequencing of the venom-gland transcriptome from H. modesta also uncovered an ortholog of this toxin (Figure 3B). This was entirely unexpected since this toxin has a vastly different primary structure to ω-HXTX-Hv1a . Moreover, in addition to the six conserved cysteine residues in ω-HXTX-Hv1a that form an inhibitor cystine knot (ICK) motif [40, 41], κ-HXTX-Hv1c contains two additional cysteine residues that form an extremely rare vicinal disulfide bond [42–44]. Furthermore, in contrast to ω-HXTX-Hv1a, which blocks insect CaV channels, κ-HXTX-Hv1c is a potent and specific blocker of KCa channels . Nevertheless, the near identity of the signal sequence in these two toxin families and the conservation of cysteine residues in the mature toxin indicate that they evolved from the same ancestral toxin gene and are members of the same gene superfamily.
We did not find orthologs of κ-HXTX-Hv1c in any of the other three species of Australian funnel-web spider (H. infensa, A. robustus, and H. venenata). However, κ-HXTX-Hv1a, κ-HXTX-Hv1b, and κ-HXTX-Hv1c are expressed at very low levels in H. versuta venom , and consequently we cannot rule out the possibility that these toxins are present in the venom of the other three spiders but the transcript levels are too low to be detected using the methods employed here.
The Shiva superfamily of peptide toxins
It has previously been suggested that superfamilies of spider-venom peptides evolved from a single ancestral gene via explosive gene duplication ; the work described here further supports this idea as it is clear that the ω-HXTXs, ω/κ-HXTXs, and κ-HXTXs belong to a large superfamily of toxins that arose via gene duplication (Figure 2C). We have chosen to name spider-toxin gene superfamilies after deities of death and destruction since the major biological role of these toxins is to paralyze and/or kill envenomated prey. Accordingly, we have named the ω/κ-HXTX/ω-HXTX/κ-HXTX gene superfamily after the Hindu deity Shiva, commonly known as the “destroyer”.
There are two striking differences between the Shiva superfamily precursors and transcripts encoding human neuropeptides and other secreted proteins. First, whereas precursors of human neuropeptides often encode multiple mature neuropeptide sequences [29, 52], we and others have not found any examples of spider-toxin transcripts that encode more than a single mature toxin sequence. Secondly, in direct contrast to the toxin precursors, the sequence of the mature human neuropeptide(s) is usually strongly conserved whereas there is significantly more variability in the signal sequence. This is perhaps not surprising given that human neuropeptides usually act on a single well-defined molecular target whereas spider toxins typically target a specific subtype of receptor or ion channel that nevertheless might vary significantly in primary structure between prey taxa (Note that most spiders are generalist predators that target a phylogenetically diverse range of prey). Thus, expressing a family of related toxins in the venom (essentially a mini-combinatorial peptide library) might ensure that the desired receptor/ion channel is targeted, regardless of prey taxa.
Position-specific cysteine codon bias
Mature ω-HXTXs contain three disulfide bonds with 1–4, 2–5, 3–6 connectivity. These disulfides form an ICK motif that provides these toxins with a high degree of chemical, thermal and biological stability . Although it is clear from a protein structure viewpoint why these six cysteine residues need to be strictly maintained in order to preserve the toxin’s 3D scaffold, one would not expect to find a preference for either one of the two possible cysteine codons (TGT and TGC). Intriguingly, however, previous analysis of ω-HXTX precursors revealed a strong bias for TGC at four of the six cysteine positions in the mature toxin region . An extended logo analysis  incorporating all of the newly discovered sequences reported in this study corroborated the previously observed codon bias (Figure 4B). We found an extreme TGC codon bias for the four cysteine residues that form the 1–4 and 3–6 disulfide bridges in the ω-HXTX family but not for the two cysteines that form the 2–5 disulfide bond (Figure 4B). The observed position-specific codon bias is not simply a manifestation of global codon bias in these spiders as we have observed a preference for TGT as opposed to TGC for cysteine residues in other hexatoxin superfamilies (data not shown). Moreover, we did not observe extreme codon bias for any other conserved residue in the mature hexatoxins.
Position-specific cysteine-codon bias has also been observed in superfamilies of cone snail toxins and it has been proposed that these codons might serve as attractants for a mutator complex that includes a poorly processive and highly mutagenic polymerase (e.g., DNA Pol V) that promotes radiation of the toxin superfamily by facilitating hypermutation of the mature toxin region [49, 50]. However, there is currently no direct evidence that cysteine-codon bias plays a part in directing the evolution of spider or cone snail toxins.
Molecular evolution analyses
Nucleotide and complementary protein analyses for ω toxins
Accessible surface areaf
0.99 ± 0.38
0.87 ± 0.42
1.41 ± 0.49
1.47 ± 0.37
M W , M V , V 0 , μ
8, 8, 8, 8
1.25 ± 0.37
1.28 ± 0.41
M W , M V , V 0 , μ
8, 8, 8, 8
1.62 ± 0.63
1.55 ± 0.36
1.39 ± 0.48
1.45 ± 0.38
0.80 ± 0.44
0.71 ± 0.41
0.53 ± 0.40
0.51 ± 0.30
1.34 ± 0.45
1.38 ± 0.40
1.10 ± 0.33
1.01 ± 0.42
Molecular evolution of ω and κ toxins from Australian funnel-web spiders
ω >1d: 3
ω <1e: 5
ω >1a: 1
ω <1b: 0
ω >1a: 3
ω <1b: 7
(0 + 1)
To derive further support for the positively selected sites detected by nucleotide analyses, we employed a complementary protein-level approach implemented in TreeSAAP (Table 1). TreeSAAP identified two positively selected sites in the ω-HXTXs that were in common with the sites identified by site-model 8 of PAML (Table 1). Evolutionary fingerprint analyses (Figure 5A) clearly revealed several residues in the ω-HXTXs and the combined toxin dataset that evolve under the influence of positive selection, while a majority of residues in the κ-HXTXs remained under evolutionary constraint (Figure 5A,B). Thus, evolution of the ω-HXTXs has been significantly influenced by short bursts of episodic adaptations, while the κ-HXTXs appear to be under negative selection.
Constraints on mutation of the mature toxin sequence
Figure 7B and C show schematic representations of the 3D structure of κ-HXTX-Hv1c and ω-HXTX-Hv1a, respectively. The two toxins can be considered to comprise four inter-cystine loops, which are labelled 1–4 from N- to C-terminus. Although there is an obvious similarity in the disposition of the three centrally located disulfide bridges that form the cystine-knot motif in each toxin, the overall topology of the toxins, as well as the size and relative orientation of the four inter-cystine loops, appears quite different. However, the structural overlay in Figure 7D, which was generated automatically by the DaliLite structural alignment program , reveals that the two structures are in fact remarkably similar.
The DaliLite alignment yields a root mean square deviation of 2.4 Å over the backbone atoms of the 28 aligned residues, indicating that the two toxins are indeed structural homologs. The three central disulfide bridges and loop 1 align remarkably well. Loop 4, which encompasses the β-hairpin present in both toxins, also aligns well except for the four-residue insertion in ω-HXTX-Hv1a (Figure 7A), which increases the size of the hairpin loop at the tip of Loop 4. The major structural differences between the two toxins are the very different orientations of loops 2 and 3. However, these structural variations cannot disguise the fact that the two toxins essentially conform to the same 3D scaffold despite their extraordinary sequence divergence (16% identity if the cysteine framework is excluded). This ability to maintain a consistent molecular architecture despite massive variation in the inter-cystine loop sequences has important implications for the mechanism by which this superfamily of peptide toxins has evolved.
Spiders and other venomous animals rely on the production of pharmacologically complex venoms for defense, prey capture, and competitor deterrence. The major components of most spider venoms are disulfide-rich peptides that have evolved to target a wide range of receptors and ion channels in the insect nervous system. The ω-HXTX and κ-HXTX families were the first peptides isolated from Australian funnel-web spiders that were shown to be insecticidal . Analysis of all transcripts encoding these peptides showed that they are initially expressed as prepropeptides that are proteolytically processed to yield a 36–37 residue mature peptide that contains three disulfide bridges that form an ICK motif plus a non-canonical vicinal disulfide bond in the κ-HXTXs.
The extreme diversity of primary structure within the Shiva toxin superfamily suggests that there have been few evolutionary restraints on sequence diversification outside of the disulfide bridges that direct the 3D fold of these peptides. The ω-HXTXs, in particular, seem to have evolved under the influence of positive Darwinian selection in an episodic fashion, whereas the κ-HXTXs appear to be constrained by negative selection pressures. Functional assessments of these toxins should shed further light on why they have adopted quite contrasting molecular evolutionary regimes. ω-HXTXs were also found to have adopted RAVER, where a large number of the episodically diversifying sites are concentrated on the molecular surface, facilitating the generation of novel pharmacological sites. These toxins may therefore be good candidates for in vitro evolution studies designed to produce modified peptides with desired therapeutic  or agrochemical  properties. Most importantly, this study reinforces the idea that the remarkable chemical and pharmacological complexity of spider venoms may be derived from a relatively small number of ancestral genes.
Identification of ω-HXTX-1 homologs via rapid amplification of cDNA ends
Venom-gland cDNA libraries were prepared from individual specimens of the following species of Australian funnel-web spider (Arthropoda: Chelicerata: Arachnida: Araneae: Opisthothelae: Mygalomorphae: Hexathelidae): Hadronyche infensa, H. versuta, H. venenata, and Atrax robustus, which were collected from geographically distinct regions of Australia (Figure 1A). Spiders were cooled to -20°C for 40–60 min, then paired venom glands were carefully dissected from each specimen (Figure 1C). Each pair of venom glands was combined, then polyA + mRNA was extracted using a QuickPrep Micro mRNA Purification Kit (Amersham Pharmacia Biotech/GE Healthcare Life Sciences, Rydalmere, NSW, Australia) and stored at -20°C until further use. cDNA libraries were constructed using a Marathon cDNA Amplification Kit (Clontech Laboratories, Mount View, CA, USA). From the mRNA template, single-stranded cDNA was constructed using Superscript III reverse transcriptase (Life Technologies, Grand Island, NY, USA) and a poly(dT) anchor primer (Echoclonanch-2, GGGCAGGT17). Second-strand synthesis was performed according to the kit specifications except the cDNA was purified using a Concert Rapid PCR Purification kit (Gibco/Life Technologies) instead of a phenol/chloroform extraction. A Marathon cDNA Amplification adaptor (Clontech Laboratories) was then ligated to the double-stranded cDNA. After overnight ligation, the sample was precipitated using 10 μl of 5% w/w glycogen, 10 μl of 3 M sodium acetate pH 5.2, and 100 μl of 100% ethanol at -20°C. The sample was subsequently washed with 80% ethanol and dried for 10 min prior to resuspension in Tris-EDTA buffer.
Transcripts encoding ω-HXTX-Hv1a and paralogs/orthologs thereof were subsequently obtained via rapid amplification of cDNA ends (RACE) (Figure 2B) . A redundant 3’ PCR primer based on residues 24–31 of ω-HXTX-Hv1a (ω-HV1 5’ RTTNCCRTTYTCRTTYTCYTTRAA 3’) was used in conjunction with a 5’ universal adaptor primer in a 5’ RACE experiment designed to extract information about the upstream region of the ω-HXTX-Hv1a transcript (EchoAP1: 5’ CACCCCTAATACGACTCACTATAGG 3’). A gene-specific primer for 3’ RACE was then designed based on the leader sequence obtained from the 5’ RACE experiment (3’ RACE primer: 5’ TGCTGCAATATGAATACCGC 3’. This primer was used in combination with the Echoclonanch-2 oligo(dT) primer (5’ GGGCAGGTTTTTTTTTTTTTTTTT 3’) to generate transcripts that encode a signal sequence homologous to that of ω-HXTX-Hv1a.
PCR products were extracted from agarose gels using a Gibco gel purification kit, precipitated using Pellet Paint Co-Precipitant kit (Novagen/EMD Millipore, Billerica, MA, USA), then phosphorylated with kinase in preparation for cloning. PCR products were then ligated into pSMART and transformed into E. cloni cells (Lucigen, Middleton, WI, USA) using the Lucigen CloneSmart Blunt Cloning kit. Transformed clones were cultured for one hour in Recovery Medium, then plated with 50 μg/mL ampicillin to allow for overnight growth. PCR screening was then used to select colonies with the expected insert size for DNA sequencing. DNA sequences (and the corresponding protein sequences) were collated and analysed using Geneious Pro, version 3.8.5  and signal sequence cleavage sites were predicted using SignalP, version 3.0 .
Identification of ω-HXTX-1 homologs via transcriptomics
Paired venom glands from Hadronyche modesta were dissected out and pooled. Total RNA was extracted using the standard TRIzol® Plus method (Invitrogen/Life Technologies) according to the manufacturer’s protocol. One microgram of Total RNA was used to construct a cDNA library using the CREATOR™ SMART™ cDNA library construction kit (Clontech Laboratories) following the manufacturer’s protocol. Briefly, total RNA was reversed transcribed using the SMART™ Moloney Murine Leukemia virus (MMLV) reverse transcriptase. Second-strand synthesis was completed using long distance polymerase chain reaction (PCR) as follows: 1 min at 95°C followed by 20 cycles of 1 min at 95°C and 6 min at 68°C. Products were then digested and size fractionated using a CHROMA SPIN-400 DEPC-H2O column (Clontech Laboratories) and then ligated into the pDNR-lib donor vector. Recombinant plasmids were electroporated into E-shot™ DHB10™-T1R electro competent cells (Invitrogen/Life Technologies). 384 clones were randomly selected and sequenced by capillary electrophoresis on an Applied Biosystems 3730×l DNA analyzer (Applied Biosystems/Life Technologies) at the Brisbane node of the Australian Genome Research Facility (AGRF). Sequences were processed so vector and polyA + tails were clipped using CLC Main Work Bench (CLC-Bio), and the Blast2GO bioinformatic suite [74, 75] was used to provide Gene Ontology, BLAST and domain/Interpro annotation. Signal sequence cleavage sites were predicted using SignalP, version 3.0 .
In accordance with the recently introduced systematic nomenclature for naming peptide toxins from venomous animals , ω-ACTX-Hv1a  and J-ACTX-Hv1c  have been renamed ω-HXTX-Hv1a and κ-HXTX-Hv1a, respectively, and the various paralogs and orthologs uncovered in this study have been named accordingly. Briefly, the Greek letter denotes the molecular target of the peptide, followed by the generic name indicating the family from which the toxin is derived; in this case the abbreviation is HXTX for hexatoxin. After the generic family name, a two-letter abbreviation is used to denote the genus and species, indicated by upper and lowercase letters respectively (i.e., Hv for H. versuta, Hi for H. infensa, etc.). The name of the species is immediately followed by a numeral that helps to distinguish different toxins with similar pharmacology and this number is followed by letter that denotes the paralog number (this is based on the number of different encoded mature toxin sequences).
Molecular evolution analyses
A total of 73 nucleotide and 90 peptide sequences were aligned using the default settings in Geneious Pro, version 3.8.5  then manually adjusted for optimal alignment prior to the following molecular evolution analyses (see Additional file 1: Figure S1 and S2).
Test for recombination
To overcome the effects of recombination on the phylogenetic and evolutionary interpretations , we employed Single Breakpoint algorithms implemented in the HyPhy package and assessed the effect of recombination on all the toxin forms examined in this study [78, 79]. When potential breakpoints were detected using the small sample Akaike information criterion (AIC), the sequences were compartmentalized or partitioned before conducting selection analyses to allow the recombining units to have distinct phylogenies (as described in [80, 81]).
We evaluated the influence of natural selection on the toxins using maximum-likelihood models [82, 83] implemented in CODEML of the PAML software . We employed site-specific models that estimate positive selection statistically as an ω value significantly greater than 1. We compared likelihood values for three pairs of models with different assumed ω distributions as no a priori expectation exists for ω: M0 (constant ω rates across all sites) versus M3 (allows ω to vary across sites within n discrete categories, where n ≥ 3); M1a (a model of neutral evolution) where all sites are assumed to be either under negative (ω <1) or neutral selection (ω = 1) versus M2a (a model of positive selection) which in addition to the site classes mentioned for M1a, assumes a third category of sites; sites with ω >1 (positive selection) and M7 (β) versus M8 (β and ω), and models that mirror the evolutionary constraints of M1 and M2 but assume that ω values are drawn from a β distribution . Only if the alternative models (M3, M2a and M8: allow sites with ω >1) show a better fit in a Likelihood Ratio Test (LRT) relative to their null models (M0, M1a and M7: do not allow sites ω >1), are their results considered significant. LRT is estimated as twice the difference in maximum likelihood values between nested models and compared with the χ2 distribution with the appropriate degree of freedom—the difference in the number of parameters between the two models. The Bayes empirical Bayes (BEB) approach  was used to identify amino acids under positive selection by calculating the posterior probabilities that a particular amino acid belongs to a given selection class (neutral, conserved, or highly variable). Sites with greater posterior probability (PP ≥ 95%) of belonging to the 'ω > 1 class’ were inferred to be positively selected.
FUBAR [55–57] implemented in HyPhy [78, 79] was employed to detect sites evolving under positive and negative selection. MEME , which is designed to overcome the drawbacks of site-specific assessments, was used to detect episodic diversifying selection. Mutations were also assessed via a complementary protein-level approach implemented in TreeSAAP . An evolutionary fingerprint analysis was carried out using the ESD algorithm implemented in datamonkey [78, 79, 88] in order to clearly depict the proportion of sites under different regimes of selection.
Logo plots  showing cysteine codon bias were constructed using Geneious software, version 5.4.
Shiva superfamily phylogenetic tree
The molecular evolutionary history of the Shiva superfamily toxins was reconstructed using Bayesian inference as implemented in MrBayes version 3.2.1 , using lset rates = invgamma with the prset aamodelpr = mixed command, which enables the program to optimize between the nine different amino acid substitution matrices implemented in MrBayes. WAG  was chosen as the best substitution matrix by the program. Tree searches were run using four Markov chains for a minimum of 10 million generations, sampling every 100th tree. The log likelihood score of each saved tree was plotted against the number of generations to establish the point at which the log-likelihood scores of the analyses reached their asymptote. 25% of the total trees sampled were discarded as burnin. The posterior probabilities for clades were established by constructing a majority rule consensus tree for all trees generated after completion of the burnin. The tree was rooted using the sequence of ω-actinopoditoxin-Mb1a (also known as ω-missulenatoxin-Mb1a) from the Eastern mouse spider Missulena bradleyi as an outgroup; this toxin also blocks CaV channels and it has sequence homology to the ω-HXTXs, including conservation of the ω-HXTX-Hv1a pharmacophore. However, although M. bradleyi and Australian funnel-web spiders (family Hexathelidae) both belong to the infraorder Mygalomorphae, M. bradleyi is a member of the Actinopodidae family.
Structural alignment of omega and kappa hexatoxins
Atomic coordinates for ω-HXTX-Hv1a  and κ-HXTX-Hv1c  were downloaded from the Protein DataBank (PDB accession codes 1AXH and 1DL0, respectively). A structural alignment of the toxins was automatically generated using DaliLite . All structure figures were generated using MacPyMOL . The Consurf webserver  was used for mapping evolutionary selection pressures on 3D homology models.
Availability of supporting data
Nucleic acid and protein sequence alignments and their respective accession numbers can be accessed from the supplementary material along with tables relevant to the molecular evolution analyses.
Rapid amplification of cDNA ends
Pico titre plate
Australian genome research facility
- H. infensa:
- H. modesta:
- H. versuta:
- H. venenata:
- A. robustus:
Mixed effects model evolution
Single likelihood ancestor counting
Random effects likelihood
Fast unbiased approximate bayesian
Phylogenetic analysis by maximum likelihood
Bayes empirical Bayes
Accessible surface area
Rapid Accumulation of Variations in Exposed Residues.
This work was supported by the Australian Research Council (Discovery Grants DP0774245 and Grant DP1095728 to GFK) and the U.S. National Science Foundation (MCB9983242 to GFK). BGF gratefully acknowledges financial support from the Australian Research Council (ARC) and The University of Queensland. KS was funded by a PhD grant (SFRH/BD/61959/2009) from F.C.T (Fundação para a Ciência e a Tecnologia). AA was funded by the project PTDC/AAC-AMB/121301/2010 (FCOMP-01-0124-FEDER-019490) and PesT-C/MAR/LA0015/2011 from F.C.T. SSP was supported by an International Postgraduate Research Scholarship (IPRS) from The University of Queensland. EABU was supported by The University of Queensland (IPRS, UQ Centennial Scholarship, and UQ Advantage Top-Up Scholarship) and the Norwegian State Education Loans Fund. The authors thank Mr Glen Gregg for collection of Hadronyche venenata specimens.
- Fry BG, Roelants K, Champagne DE, Scheib H, Tyndall JDA, King GF, Nevalainen TJ, Norman JA, Lewis RJ, Norton RS, Renjifo C, de la Vega RC R: The toxicogenomic multiverse: convergent recruitment of proteins into animal venoms. Annu Rev Genom Hum Genet. 2009, 10: 483-511. 10.1146/annurev.genom.9.081307.164356.View ArticleGoogle Scholar
- Casewell NR, Wüster W, Vonk FJ, Harrison RA, Fry BG: Complex cocktails: the evolutionary novelty of venoms. Trends Ecol Evol. 2013, 28: 219-229. 10.1016/j.tree.2012.10.020.PubMedView ArticleGoogle Scholar
- Fry BG, Wüster W, Kini RM, Brusic V, Khan A, Venkataraman D, Rooney AP: Molecular evolution and phylogeny of elapid snake venom three-finger toxins. J Mol Evol. 2003, 57: 110-129. 10.1007/s00239-003-2461-2.PubMedView ArticleGoogle Scholar
- Fry BG, Wüster W: Assembling an arsenal: origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences. Mol Biol Evol. 2004, 21: 870-883. 10.1093/molbev/msh091.PubMedView ArticleGoogle Scholar
- Fry BG: From genome to venome: molecular origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences and related body proteins. Genome Res. 2005, 15: 403-420. 10.1101/gr.3228405.PubMed CentralPubMedView ArticleGoogle Scholar
- Nei M, Gu X, Sitnikova T: Evolution by the birth-and-death process in multigene families of the vertebrate immune system. Proc Natl Acad Sci U S A. 1997, 94: 7799-7806. 10.1073/pnas.94.15.7799.PubMed CentralPubMedView ArticleGoogle Scholar
- Olivera BM, Hillyard DR, Marsh M, Yoshikami D: Combinatorial peptide libraries in drug design: lessons from venomous cone snails. Trends Biotechnol. 1995, 13: 422-426. 10.1016/S0167-7799(00)88996-9.PubMedView ArticleGoogle 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: 131-139. 10.1016/j.peptides.2004.07.016.PubMedView ArticleGoogle Scholar
- Escoubas P: Molecular diversification in spider venoms: a web of combinatorial peptide libraries. Mol Divers. 2006, 10: 545-554. 10.1007/s11030-006-9050-4.PubMedView ArticleGoogle Scholar
- de la Vega RC R, Schwartz EF, Possani LD: Mining on scorpion venom biodiversity. Toxicon. 2010, 56: 1155-1161. 10.1016/j.toxicon.2009.11.010.View ArticleGoogle Scholar
- King GF, Hardy MC: Spider-venom peptides: structure, pharmacology, and potential for control of insect pests. Annu Rev Entomol. 2013, 58: 475-496. 10.1146/annurev-ento-120811-153650.PubMedView ArticleGoogle Scholar
- Smith JJ, Herzig V, King GF, Alewood PF: The insecticidal potential of venom peptides. Cell Mol Life Sci. 2013, 70: 3665-3693. 10.1007/s00018-013-1315-3.PubMedView ArticleGoogle Scholar
- King GF, Escoubas P, Nicholson GM: Peptide toxins that selectively target insect NaV and CaV channels. Channels. 2008, 2: 100-116. 10.4161/chan.2.2.6022.PubMedView ArticleGoogle Scholar
- King GF: Venoms as a platform for human drugs: translating toxins into therapeutics. Expert Opin Biol Ther. 2011, 11: 1469-1484. 10.1517/14712598.2011.621940.PubMedView ArticleGoogle Scholar
- Lewis RJ, Dutertre S, Vetter I, Christie MJ: Conus venom peptide pharmacology. Pharmacol Rev. 2012, 64: 259-298. 10.1124/pr.111.005322.PubMedView ArticleGoogle Scholar
- Escoubas P, Rash L: Tarantulas: eight-legged pharmacists and combinatorial chemists. Toxicon. 2004, 43: 555-574. 10.1016/j.toxicon.2004.02.007.PubMedView ArticleGoogle Scholar
- Tedford HW, Gilles N, Menez A, Doering CJ, Zamponi GW, King GF: Scanning mutagenesis of ω-atracotoxin-Hv1a reveals a spatially restricted epitope that confers selective activity against insect calcium channels. J Biol Chem. 2004, 279: 44133-44140. 10.1074/jbc.M404006200.PubMedView ArticleGoogle Scholar
- Fletcher JI, Smith R, O’Donoghue SI, Nilges M, Connor M, Howden MEH, Christie MJ, King GF: The structure of a novel insecticidal neurotoxin, ω-atracotoxin-HV1, from the venom of an Australian funnel web spider. Nat Struct Biol. 1997, 4: 559-566. 10.1038/nsb0797-559.PubMedView ArticleGoogle Scholar
- Wang X, Smith R, Fletcher JI, Wilson H, Wood CJ, Howden ME, King GF: Structure-function studies of ω-atracotoxin, a potent antagonist of insect voltage-gated calcium channels. Eur J Biochem. 1999, 264: 488-494. 10.1046/j.1432-1327.1999.00646.x.PubMedView ArticleGoogle Scholar
- Chong Y, Hayes JL, Sollod B, Wen S, Wilson DT, Hains PG, Hodgson WC, Broady KW, King GF, Nicholson GM: The ω-atracotoxins: selective blockers of Insect M-LVA and HVA calcium channels. Biochem Pharmacol. 2007, 74: 623-638. 10.1016/j.bcp.2007.05.017.PubMedView ArticleGoogle Scholar
- Tedford HW, Fletcher JI, King GF: Functional significance of the β-hairpin in the insecticidal neurotoxin ω-atracotoxin-Hv1a. J Biol Chem. 2001, 276: 26568-26576. 10.1074/jbc.M102199200.PubMedView ArticleGoogle Scholar
- Tedford HW, Maggio F, Reenan RA, King GF: A model genetic system for testing the in vivo function of peptide toxins. Peptides. 2007, 28: 51-56. 10.1016/j.peptides.2006.08.026.PubMedView ArticleGoogle Scholar
- Khan SA, Zafar Y, Briddon RW, Malik KA, Mukhtar Z: Spider venom toxin protects plants from insect attack. Transgenic Res. 2006, 15: 349-357. 10.1007/s11248-006-0007-2.PubMedView ArticleGoogle Scholar
- Diniz MR, Paine MJ, Diniz CR, Theakston RD, Crampton JM: Sequence of the cDNA coding for the lethal neurotoxin Tx1 from the Brazilian “armed” spider Phoneutria nigriventer predicts the synthesis and processing of a preprotoxin. J Biol Chem. 1993, 268: 15340-15342.PubMedGoogle Scholar
- Krapcho KJ, Kral RM, Vanwagenen BC, Eppler KG, Morgan TK: Characterization and cloning of insecticidal peptides from the primitive weaving spider Diguetia canities. Insect Biochem Mol Biol. 1995, 25: 991-1000. 10.1016/0965-1748(95)00029-U.PubMedView ArticleGoogle Scholar
- Wang X-H, Connor M, Wilson D, Wilson HI, Nicholson GM, Smith R, Shaw D, Mackay JP, Alewood PF, Christie MJ, King GF: Discovery and structure of a potent and highly specific blocker of insect calcium channels. J Biol Chem. 2001, 276: 40806-40812.Google Scholar
- Ostrow KL, Mammoser A, Suchyna T, Sachs F, Oswald R, Kubo S, Chino N, Gottlieb PA: cDNA sequence and in vitro folding of GsMTx4, a specific peptide inhibitor of mechanosensitive channels. Toxicon. 2003, 42: 263-274. 10.1016/S0041-0101(03)00141-7.PubMedView ArticleGoogle Scholar
- Liang S: An overview of peptide toxins from the venom of the Chinese bird spider Selenocosmia huwena Wang [Ornithoctonus huwena (Wang)]. Toxicon. 2004, 43: 575-585. 10.1016/j.toxicon.2004.02.005.PubMedView ArticleGoogle Scholar
- Hook V, Funkelstein L, Lu D, Bark S, Wegrzyn J, Hwang S-R: Proteases for processing proneuropeptides into peptide neurotransmitters and hormones. Annu Rev Pharmacol Toxicol. 2008, 48: 393-423. 10.1146/annurev.pharmtox.48.113006.094812.PubMed CentralPubMedView ArticleGoogle Scholar
- Veenstra JA: Mono- and dibasic proteolytic cleavage sites in insect neuroendocrine peptide precursors. Arch Insect Biochem Physiol. 2000, 43: 49-63. 10.1002/(SICI)1520-6327(200002)43:2<49::AID-ARCH1>3.0.CO;2-M.PubMedView 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
- Chen J, Deng M, He Q, Meng E, Jiang L, Liao Z, Rong M, Liang S: Molecular diversity and evolution of cystine knot toxins of the tarantula Chilobrachys jingzhao. Cell Mol Life Sci. 2008, 65: 2431-2444. 10.1007/s00018-008-8135-x.PubMedView ArticleGoogle Scholar
- Wong ES, Hardy MC, Wood D, Bailey T, King GF: SVM-based prediction of propeptide cleavage sites in spider toxins identifies toxin innovation in an Australian tarantula. PLoS One. 2013, 8: e66279-10.1371/journal.pone.0066279.PubMed CentralPubMedView ArticleGoogle Scholar
- Raven RJ: The spider infraorder Mygalomorphae (Araneae): cladistics and systematics. Bull Am Mus Nat Hist. 1986, 182: 1-180.Google Scholar
- Tedford HW, Sollod BL, Maggio F, King GF: Australian funnel-web spiders: master insecticide chemists. Toxicon. 2004, 43: 601-618. 10.1016/j.toxicon.2004.02.010.PubMedView ArticleGoogle Scholar
- Gray MR: A revision of the Australian funnel-web spiders (Hexathelidae: Atracinae). Records Austr Mus. 2010, 62: 285-392. 10.3853/j.0067-1975.62.2010.1556.View ArticleGoogle Scholar
- Mukherjee AK, Sollod BL, Wikel SK, King GF: Orally active acaricidal peptide toxins from spider venom. Toxicon. 2006, 47: 182-187. 10.1016/j.toxicon.2005.10.011.PubMedView ArticleGoogle Scholar
- King GF, Sollod-Mcfarland B, Nicholson GM, Gunning SJ: Insecticidal Polypeptides and Methods of use Thereof U.S. Patent No. 7,354,993. 2008Google Scholar
- Wang X-H, Connor M, Smith R, Maciejewski MW, Howden MEH, Nicholson GM, Christie MJ, King GF: Discovery and characterization of a family of insecticidal neurotoxins with a rare vicinal disulfide bond. Nat Struct Biol. 2000, 7: 505-513. 10.1038/75921.PubMedView ArticleGoogle Scholar
- Pallaghy PK, Nielsen KJ, Craik DJ, Norton RS: A common structural motif incorporating a cystine knot and a triple-stranded β-sheet in toxic and inhibitory polypeptides. Protein Sci. 1994, 3: 1833-1839. 10.1002/pro.5560031022.PubMed CentralPubMedView ArticleGoogle Scholar
- King GF, Tedford HW, Maggio F: Structure and function of insecticidal neurotoxins from Australian funnel-web spiders. J Toxicol-Toxin Rev. 2002, 21: 359-389.View ArticleGoogle Scholar
- Wang G, Peterkofsky A, Clore GM: A novel membrane anchor function for the N-terminal amphipathic sequence of the signal-transducing protein IIAGlucose of the Escherichia coli phosphotransferase system. J Biol Chem. 2000, 275: 39811-39814. 10.1074/jbc.C000709200.PubMedView ArticleGoogle Scholar
- Maggio F, King GF: Scanning mutagenesis of a Janus-faced atracotoxin reveals a bipartite surface patch that is essential for neurotoxic function. J Biol Chem. 2002, 277: 22806-22813. 10.1074/jbc.M202297200.PubMedView ArticleGoogle Scholar
- Mobli M, de Araujo AD, Lambert LK, Pierens GK, Windley MJ, Nicholson GM, Alewood PF, King GF: Direct visualization of disulfide bonds through diselenide proxies using 77Se NMR spectroscopy. Angew Chem Int Ed. 2009, 48: 9312-9314. 10.1002/anie.200905206.View ArticleGoogle Scholar
- Gunning SJ, Maggio F, Windley MJ, Valenzuela SM, King GF, Nicholson GM: The Janus-faced atracotoxins are specific blockers of invertebrate KCa channels. FEBS J. 2008, 275: 4045-4059. 10.1111/j.1742-4658.2008.06545.x.PubMedView ArticleGoogle Scholar
- Woodward SR, Cruz LJ, Olivera BM, Hillyard DR: Constant and hypervariable regions in conotoxin propeptides. EMBO J. 1990, 9: 1015-1020.PubMed CentralPubMedGoogle 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 NY Acad Sci. 1999, 870: 223-237. 10.1111/j.1749-6632.1999.tb08883.x.PubMedView ArticleGoogle 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 U S A. 1999, 96: 6820-6823. 10.1073/pnas.96.12.6820.PubMed CentralPubMedView ArticleGoogle Scholar
- Conticello SG, Pilpel Y, Glusman G, Fainzilber M: Position-specific codon conservation in hypervariable gene families. Trends Genet. 2000, 16: 57-59. 10.1016/S0168-9525(99)01956-3.PubMedView ArticleGoogle 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: 120-131. 10.1093/oxfordjournals.molbev.a003786.PubMedView ArticleGoogle Scholar
- Espiritu DJD, Watkins M, Dia-Monje V, Cartier GE, Cruz LJ, Olivera BM: Venomous cone snails: molecular phylogeny and the generation of toxin diversity. Toxicon. 2001, 39: 1899-1916. 10.1016/S0041-0101(01)00175-1.PubMedView ArticleGoogle Scholar
- Fricker LD: Neuropeptide-processing enzymes: applications for drug discovery. AAPS J. 2005, 7: E449-E455. 10.1208/aapsj070244.PubMed CentralPubMedView ArticleGoogle Scholar
- Saez NJ, Senff S, Jensen JE, Er SY, Herzig V, Rash LD, King GF: Spider-venom peptides as therapeutics. Toxins. 2010, 2: 2851-2871. 10.3390/toxins2122851.PubMed CentralPubMedView ArticleGoogle Scholar
- Schneider TD, Stephens RM: Sequence logos: a new way to display consensus sequences. Nucl Acids Res. 1990, 18: 6097-6100. 10.1093/nar/18.20.6097.PubMed CentralPubMedView ArticleGoogle Scholar
- Kosakovsky Pond SL, Frost SDW: Not so different after all: a comparison of methods for detecting amino acid sites under selection. Mol Biol Evol. 2005, 22: 1208-1222. 10.1093/molbev/msi105.PubMedView ArticleGoogle Scholar
- Murrell B, Moola S, Mabona A, Weighill T, Sheward D, Kosakovsky Pond SL, Scheffler K: FUBAR: a Fast, Unconstrained Bayesian AppRoximation for inferring selection. Mol Biol Evol. 2013, 30: 1196-1205. 10.1093/molbev/mst030.PubMed CentralPubMedView ArticleGoogle Scholar
- Murrell B, Wertheim JO, Moola S, Weighill T, Scheffler K, Kosakovsky Pond SL: Detecting individual sites subject to episodic diversifying selection. PLoS Genet. 2012, 8: e1002764-10.1371/journal.pgen.1002764.PubMed CentralPubMedView ArticleGoogle Scholar
- Kozminsky-Atias A, Zilberberg N: Molding the business end of neurotoxins by diversifying evolution. FASEB J. 2012, 26: 576-586. 10.1096/fj.11-187179.PubMedView ArticleGoogle Scholar
- Sunagar K, Jackson TNW, Undheim EAB, Ali SA, Antunes A, Fry BG: Three-fingered RAVERs: rapid accumulation of variation in exposed residues of snake venom toxins. Toxins. 2013, 5: 2172-2208. 10.3390/toxins5112172.PubMed CentralPubMedView ArticleGoogle Scholar
- Kini RM, Chan YM: Accelerated evolution and molecular surface of venom phospholipase A2 Enzymes. J Mol Evol. 1999, 48: 125-132. 10.1007/PL00006450.PubMedView ArticleGoogle Scholar
- Brust A, Sunagar K, Undheim EAB, Vetter I, Yang DC, Casewell NR, Jackson TNW, Koludarov I, Alewood PF, Hodgson WC, Lewis RJ, King GF, Antunes A, Hendrikx I, Fry BG: Differential evolution and neofunctionalization of snake venom metalloprotease domains. Mol Cell Proteomics. 2013, 12: 651-663. 10.1074/mcp.M112.023135.PubMed CentralPubMedView ArticleGoogle Scholar
- Ruder T, Sunagar K, Undheim EAB, Ali SA, Wai T-C, Low DHW, Jackson TNW, King GF, Antunes A, Fry BG: Molecular phylogeny and evolution of the proteins encoded by coleoid (cuttlefish, octopus, and squid) posterior venom glands. J Mol Evol. 2013, 76: 192-204. 10.1007/s00239-013-9552-5.PubMedView ArticleGoogle Scholar
- Low DHW, Sunagar K, Undheim EAB, Ali SA, Alagon AC, Ruder T, Jackson TNW, Pineda Gonzalez S, King GF, Jones Antunes A, Fry BG: Dracula’s children: molecular evolution of vampire bat venom. J Proteomics. 2013, 89: 95-111.PubMedView ArticleGoogle Scholar
- Sunagar K, Fry BG, Jackson TNW, Casewell NR, Undheim EAB, Vidal N, Ali SA, King GF, Vasudevan K, Vasconcelos V, Antunes A: Molecular evolution of vertebrate neutrophins: co-option of the highly conserved nerve groth factor gene into the advance snake arsenal. PLoS One. 2013, 8: e81827-10.1371/journal.pone.0081827.PubMed CentralPubMedView ArticleGoogle Scholar
- Nisani Z, Dunbar SG, Hayes WK: Cost of venom regeneration in Parabuthus transvaalicus (Arachnida: Buthidae). Comp Biochem Physiol A-Mol Integr Physiol. 2007, 147: 509-513. 10.1016/j.cbpa.2007.01.027.PubMedView ArticleGoogle Scholar
- Nisani Z, Boskovic DS, Dunbar SG, Kelln W, Hayes WK: Investigating the chemical profile of regenerated scorpion (Parabuthus transvaalicus) venom in relation to metabolic cost and toxicity. Toxicon. 2012, 60: 315-323. 10.1016/j.toxicon.2012.04.343.PubMedView ArticleGoogle Scholar
- Morgenstern D, King GF: The venom optimization hypothesis revisited. Toxicon. 2013, 63: 120-128.PubMedView ArticleGoogle Scholar
- de Araujo AD, Herzig V, Windley MJ, Dziemborowicz S, Mobli M, Nicholson GM, Alewood PF, King GF: Do vicinal disulfide bridges mediate functionally important redox transformations in proteins?. Antioxid Redox Signal. 2013, 19: 1976-1980. 10.1089/ars.2013.5365.PubMed CentralPubMedView ArticleGoogle Scholar
- Maggio F, King GF: Role of the structurally disordered N- and C-terminal residues in the Janus-faced atracotoxins. Toxicon. 2002, 40: 1355-1361. 10.1016/S0041-0101(02)00154-X.PubMedView ArticleGoogle Scholar
- Holm L, Kääriäinen S, Rosenström P, Schenkel A: Searching protein structure databases with DaliLite v.3. Bioinformatics. 2008, 24: 2780-2781. 10.1093/bioinformatics/btn507.PubMed CentralPubMedView ArticleGoogle Scholar
- Frohman MA: Rapid amplification of complementary DNA ends for generation of full-length complementary DNAs: thermal RACE. Methods Enzymol. 1993, 218: 340-356.PubMedView ArticleGoogle Scholar
- Drummond AJ, Ashton B, Buxton S, Cheung M, Cooper A, Duran C, Field M, Heled J, Kearse M, Markowitz S, Moir R, Stones-Havas S, Sturrock S, Thierer T, Wilson A: Geneious v5.4. 2011. 2011, 2011: 2011-Available from http://www.geneious.comGoogle Scholar
- Bendtsen JD, Nielsen H, Von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340: 783-795. 10.1016/j.jmb.2004.05.028.PubMedView ArticleGoogle Scholar
- Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucl Acids Res. 2008, 36: 3420-3435. 10.1093/nar/gkn176.PubMed CentralPubMedView ArticleGoogle Scholar
- Götz S, Arnold R, Sebastián-León P, Martín-Rodríguez S, Tischler P, Jehl M-A, Dopazo J, Rattei T, Conesa A: B2G-FAR, a species-centered GO annotation repository. Bioinformatics. 2011, 27: 919-924. 10.1093/bioinformatics/btr059.PubMed CentralPubMedView ArticleGoogle Scholar
- King GF, Gentz MC, Escoubas P, Nicholson GM: A rational nomenclature for naming peptide toxins from spiders and other venomous animals. Toxicon. 2008, 52: 264-276. 10.1016/j.toxicon.2008.05.020.PubMedView ArticleGoogle Scholar
- Posada D, Crandall KA: The effect of recombination on the accuracy of phylogeny estimation. J Mol Evol. 2002, 54: 396-402. 10.1007/s00239-001-0034-9.PubMedView ArticleGoogle Scholar
- Kosakovsky Pond SL, Frost SDW: Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics. 2005, 21: 2531-2533. 10.1093/bioinformatics/bti320.View ArticleGoogle Scholar
- Delport W, Poon AFY, Frost SDW, Kosakovsky Pond SL: Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics. 2010, 26: 2455-2457. 10.1093/bioinformatics/btq429.PubMed CentralPubMedView ArticleGoogle Scholar
- Kosakovsky Pond SL, Posada D, Gravenor MB, Woelk CH, Frost SDW: Automated phylogenetic detection of recombination using a genetic algorithm. Mol Biol Evol. 2006, 23: 1891-1901. 10.1093/molbev/msl051.PubMedView ArticleGoogle Scholar
- Scheffler K, Martin DP, Seoighe C: Robust inference of positive selection from recombining coding sequences. Bioinformatics. 2006, 22: 2493-2499. 10.1093/bioinformatics/btl427.PubMedView ArticleGoogle Scholar
- Goldman N, Yang Z: A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol Biol Evol. 1994, 11: 725-736.PubMedGoogle Scholar
- Yang Z: Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998, 15: 568-573. 10.1093/oxfordjournals.molbev.a025957.PubMedView ArticleGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.PubMedView ArticleGoogle Scholar
- Nielsen R, Yang Z: Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics. 1998, 148: 929-936.PubMed CentralPubMedGoogle Scholar
- Yang Z, Wong WSW, Nielsen R: Bayes Empirical Bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005, 22: 1107-1118. 10.1093/molbev/msi097.PubMedView ArticleGoogle Scholar
- Woolley S, Johnson J, Smith MJ, Crandall KA, McClellan DA: TreeSAAP: selection on amino acid properties using phylogenetic trees. Bioinformatics. 2003, 19: 671-672. 10.1093/bioinformatics/btg043.PubMedView ArticleGoogle Scholar
- Kosakovsky Pond SL, Murrell B, Fourment M, Frost SDW, Delport W, Scheffler K: A random effects branch-site model for detecting episodic diversifying selection. Mol Biol Evol. 2011, 28: 3033-3043. 10.1093/molbev/msr125.PubMed CentralPubMedView ArticleGoogle Scholar
- Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP: MrBayes 3.2: efficient bayesian phylogenetic inference and model choice across a large model space. Systematic Biol. 2012, 61: 539-542. 10.1093/sysbio/sys029.View ArticleGoogle Scholar
- Whelan S, Goldman N: A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol Biol Evol. 2001, 18: 691-699. 10.1093/oxfordjournals.molbev.a003851.PubMedView ArticleGoogle Scholar
- DeLano WL: The PyMOL Molecular Graphics System. 2002, San Carlos, CA, USA: DeLano ScientificGoogle Scholar
- Glaser F, Pupko T, Paz I, Bell RE, Bechor-Shental D, Martz E, Ben-Tal N: ConSurf: identification of functional regions in proteins by surface-mapping of phylogenetic information. Bioinformatics. 2003, 19: 163-164. 10.1093/bioinformatics/19.1.163.PubMedView ArticleGoogle 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 credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.