Concerted evolution of sea anemone neurotoxin genes is revealed through analysis of the Nematostella vectensis genome

Gene families, which encode toxins, are found in many poisonous animals, yet there is limited understanding of their evolution at the nucleotide level. The release of the genome draft sequence for the sea anemone Nematostella vectensis enabled a comprehensive study of a gene family whose neurotoxin products affect voltage-gated sodium channels. All gene family members are clustered in a highly repetitive approximately 30-kb genomic region and encode a single toxin, Nv1. These genes exhibit extreme conservation at the nucleotide level which cannot be explained by purifying selection. This conservation greatly differs from the toxin gene families of other animals (e.g., snakes, scorpions, and cone snails), whose evolution was driven by diversifying selection, thereby generating a high degree of genetic diversity. The low nucleotide diversity at the Nv1 genes is reminiscent of that reported for DNA encoding ribosomal RNA (rDNA) and 2 hsp70 genes from Drosophila, which have evolved via concerted evolution. This evolutionary pattern was experimentally demonstrated in yeast rDNA and was shown to involve unequal crossing-over. Through sequence analysis of toxin genes from multiple N. vectensis populations and 2 other anemone species, Anemonia viridis and Actinia equina, we observed that the toxin genes for each sea anemone species are more similar to one another than to those of other species, suggesting they evolved by manner of concerted evolution. Furthermore, in 2 of the species (A. viridis and A. equina) we found genes that evolved under diversifying selection, suggesting that concerted evolution and accelerated evolution may occur simultaneously.     

Purification, sequence, and pharmacological properties of sea anemone toxins from Radianthus paumotensis. A new class of sea anemone toxins acting on the sodium channel

Four new toxins have been isolated from the sea anemone Radianthus paumotensis: RpI, RpII, RpIII, and RpIV. They are polypeptides comprised of 48 or 49 amino acids; the sequence of RpII has been determined. Toxicities of these toxins in mice and crabs are similar to those of the other known sea anemone toxins, but they fall into a different immunochemically defined class. The sequence of RpII shows close similarities with the N-terminal end (up to residue 20) of the previously sequenced long sea anemone toxins, but most of the remaining part of the molecule is completely different. Like the other sea anemone toxins, Radianthus toxins are active on sodium channels; they slow down the inactivation process. Through their Na+ channel action, Radianthus toxins stimulate Na+ influx into tetrodotoxin-sensitive neuroblastoma cells and tetrodotoxin-resistant rat skeletal myoblasts. The efficiency of the toxins is similar in the two cellular systems. In that respect, Radianthus toxins behave much more like scorpion neurotoxins than sea anemone toxins from Anemonia sulcata or Anthopleura xanthogrammica. In binding experiments to synaptosomal Na+ channels, Radianthus toxins compete with toxin II from the scorpion Androctonus australis but not with toxins II and V from Anemonia sulcata.     

Amino acid sequence of a sea anemone toxin from Parasicyonis actinostoloides

Amino-acid sequence of a toxin from sea anemone, Parasicyonis actinostoloides, is determined. The toxin consists of 31 amino acid residues and is cross-linked with four disulphide bridges. The sequence has some similarity to that of toxin III and no similarity to those of toxin I and toxin II both from sea anemone, Anemonia sulcata, or to that of Anthopleurin A from Anthopleura xanthogrammica.     

Optimization of the expression of equistatin in Pichia pastoris.

To improve the expression of equistatin, a proteinase inhibitor from the sea anemone Actinia equina, in the yeast Pichia pastoris, we prepared gene variants with yeast-preferred codon usage and lower repetitive AT and GC content. The full gene optimization approximately doubled the level of steady-state mRNA and protein accumulated in the culture medium. The removal of a short stretch of 12 additional nucleotides from the multiple cloning site (MCS) sequence in the vector pPIC9 had an enhancement effect similar to full gene optimization (factor 1.5) at the mRNA level. However, at the protein level, this increase was 4- to 10-fold. The optimized gene without the MCS sequence yielded 1.66 g/L active protein in a bioreactor and was purified by a new two-step procedure with a recovery of activity that was >95%. This production level constitutes an overall improvement of about 20-fold relative to our previously published results. The characteristics of the MCS sequence element are discussed in the light of its apparent ability to act as negative expression regulator.     

Calitoxin, a neurotoxic peptide from the sea anemone Calliactis parasitica: amino acid sequence and electrophysiological properties

We have isolated a new toxin, calitoxin (CLX), from the sea anemone Calliactis parasitica whose amino acid sequence differs greatly from that of other sea anemone toxins. The polypeptide chain contains 46 amino acid residues, with a molecular mass of 4886 Da and an isoelectric point at pH 5.4. The amino acid sequence determined by Edman degradation of the reduced, S-carboxymethylated polypeptide chain and tryptic and chymotryptic peptides is Ile-Glu-Cys-Lys-Cys-Glu-Gly-Asp-Ala-Pro-Asp-Leu-Ser-His-Met-Thr-Gly-Thr- Val-Tyr - Phe-Ser-Cys-Lys-Gly-Gly-Asp-Gly-Ser-Trp-Ser-Lys-Cys-Asn-Thr-Tyr-Thr-Ala- Val-Ala - Asp-Cys-Cys-His-Glu-Ala. No cysteine residues were present in the peptide. Similarly to other sea anemone toxins, calitoxin interacts, in crustacean nerve muscle preparations, with axonal and not with muscle membranes, inducing a massive release of neurotransmitter that causes a strong muscle contraction. The low homology of CLX with RP II and ATX II toxins has implications regarding the role played by particular amino acid residues.     

The amino acid sequence of toxin RpIII from the sea anemone, Radianthus paumotensis

The amino acid sequence of the sodium channel toxin RpIII from the sea anemone Radianthus paumotensis has been determined. The protein is homologous with five analogous toxins from three anemone species, and is most similar to a less toxic protein, RpII, from the same organism. Twelve residues are conserved in all six toxins, one of which is an arginine residue thought to be essential for toxicity. The others (Cys, Gly, Pro and Trp) tend to be conserved in other sets of homologous proteins to maintain functional folds. Comparisons of the sequences suggest the existence of two separate but related classes of toxins cumon the three species of anemone.     

Isolation, characterization, and amino acid sequence of a polypeptide neurotoxin occurring in the sea anemone Stichodactyla helianthus

An aqueous exudate collected from frozen and thawed bodies of a Caribbean sea anemone, Stichodactyla (formerly Stoichactis) helianthus, contained a polypeptide neurotoxin (Sh I) selectively toxic to crustaceans. The polypeptide was purified by G-50 Sephadex, phosphocellulose, and sulfopropyl-Sephadex chromatography and shown to have a molecular size of 5200 daltons and a pI of 8.3. The amino acid sequence determined by automatic Edman degradations of whole RCM Sh I and of its clostripain, staphylococcal protease, and cyanogen bromide digest peptides is A1ACKC5DDEGP10DIRTA15PLTGT20VDLGS25CNAGW30EKCAS35YYTII40ADCCR45KKK . Only 33% of this sequence is identical with the sequence of Anemonia sulcata toxin II, a sea anemone toxin isolated from the taxonomic family Actiniidae. The six half-cystines are located in equivalent positions to those of the actiniid toxins and account for nearly half of the residues common to all of the toxins. However, 69% of the Sh I sequence is identical with that of toxin II from Heteractis paumotensis, another sea anemone belonging to the family Stichodactylidae. Stichodactylid toxins lack the initial N-terminal residue of actiniid toxins and possess three consecutive acidic residues at positions 6-8, a single tryptophan at position 30, and four consecutive basic residues at positions 45-48 (C-terminus). A rabbit IgG prepared by Sh I immunization bound Sh I with a K0.5 of 4.7 nM but failed to bind homologous actiniid (Anemonia sulcata II, Condylactis gigantea III) or bolocerid (Bolocera tuedae II) polypeptide neurotoxins.(ABSTRACT TRUNCATED AT 250 WORDS)     

The 18S ribosomal RNA sequence of the sea anemone Anemonia sulcata and its evolutionary position among other eukaryotes

Evolutionary trees based on partial small ribosomal subunit RNA sequences of 22 metazoa species have been published [(1988) Science 239, 748-753]. In these trees, cnidarians (Radiata) seemed to have evolved independently from the Bilateria, which is in contradiction with the general evolutionary view. In order to further investigate this problem, the complete srRNA sequence of the sea anemone Anemonia sulcata was determined and evolutionary trees were constructed using a matrix optimization method. In the tree thus obtained the sea anemone and Bilateria together form a monophyletic cluster, with the sea anemone forming the first line of the metazoan group.     

Analysis of rare amino acid replacements supports the Coelomata clade

The recent analysis of a novel class of rare genomic changes, RGC_CAMs (after conserved amino acids-multiple substitutions), supported the Coelomata clade of animals as opposed to the Ecdysozoa clade (Rogozin et al. 2007). A subsequent reanalysis, with the sequences from the sea anemone Nematostella vectensis included in the set of outgroup species, suggested that this result was an artifact caused by reverse amino replacements and claimed support for Ecdysozoa (Irimia et al. 2007). We show that the internal branch connecting the sea anemone to the bilaterian animals is extremely short, resulting in a weak statistical support for the Coelomata clade. Direct estimation of the level of homoplasy, combined with taxon sampling with different sets of outgroup species, reinforces the support for Coelomata, whereas the effect of reversals is shown to be relatively minor.     

Identification of a new member of the GLWamide peptide family: physiological activity and cellular localization in cnidarian polyps

KPNAYKGKLPIGLWamide, a novel member of the GLWamide peptide family, was isolated from Hydra magnipapillata. The purification was monitored with a bioassay: contraction of the retractor muscle of a sea anemone, Anthopleura fuscoviridis. The new peptide, termed Hym-370, is longer than the other GLWamides previously isolated from H. magnipapillata and another sea anemone, A. elegantissima. The amino acid sequence of Hym-370 is six residues longer at its N-terminal than a putative sequence previously deduced from the cDNA encoding the precursor protein. The new longer isoform, like the shorter GLWamides, evoked concentration-dependent muscle contractions in both H. magnipapillata and A. fuscoviridis. In contrast, Hym-248, one of the shorter GLWamide peptides, specifically induced contraction of the endodermal muscles in H. magnipapillata. This is the first case in which a member of the hydra GLWamide family (Hym-GLWamides) has exhibited an activity not shared by the others. Polyclonal antibodies were raised to the common C-terminal tripeptide GLWamide and were used in immunohistochemistry to localize the GLWamides in the tissue of two species of hydra, H. magnipapillata and H. oligactis, and one species of sea anemone, A. fuscoviridis. In each case, nerve cells were specifically labeled. These results suggest that the GLWamides are ubiquitous among cnidarians and are involved in regulating the excitability of specific muscles.     

A case for a Glossina genome project

Given the medical and agricultural significance of Glossina, knowledge of the genomic aspects of the vector and vector-pathogen interactions are a high priority. In preparation for a full genome sequence initiative, an extensive set of expressed sequence tags (ESTs) has been generated from tissue-specific normalized libraries. In addition, bacterial artificial chromosome (BAC) libraries are being constructed, and information on the genome structure and size from different species has been obtained. An international consortium is now in place to further efforts to lead to a full genome project.     

Computational prediction of RNA editing sites

MOTIVATION: Some organisms edit their messenger RNA resulting in differences between the genomic sequence for a gene and the corresponding messenger RNA sequence. This difference complicates experimental and computational attempts to find and study genes in organisms with RNA editing even if the full genomic sequence is known. Nevertheless, knowledge of these editing sites is crucial for understanding the editing machinery of these organisms. RESULTS: We present a computational technique that predicts the position of editing sites in the genomic sequence. It uses a statistical approach drawing on the protein sequences of related genes and general features of editing sites of the organism. We apply the method to the mitochondrion of the slime mold Physarum polycephalum. It correctly predicts over 90% of the amino acids and over 70% of the editing sites.     

Short reads, circular genome: skimming solid sequence to construct the bighorn sheep mitochondrial genome

As sequencing technology improves, an increasing number of projects aim to generate full genome sequence, even for nonmodel taxa. These projects may be feasibly conducted at lower read depths if the alignment can be aided by previously developed genomic resources from a closely related species. We investigated the feasibility of constructing a complete mitochondrial (mt) genome without preamplification or other targeting of the sequence. Here we present a full mt genome sequence (16,463 nucleotides) for the bighorn sheep (Ovis canadensis) generated though alignment of SOLiD short-read sequences to a reference genome. Average read depth was 1240, and each base was covered by at least 36 reads. We then conducted a phylogenomic analysis with 27 other bovid mitogenomes, which placed bighorn sheep firmly in the Ovis clade. These results show that it is possible to generate a complete mitogenome by skimming a low-coverage genomic sequencing library. This technique will become increasingly applicable as the number of taxa with some level of genome sequence rises.     

The acclimation of anemone fishes to sea anemones: Protection by changes in the fish''s mucous coat

The anemone fish Amphiprion clarkii (Bennett) lives unharmed among the tentacles of sea anemones following acclimation. What occurs during this process is unknown; however, two hypotheses have been formulated recently. One suggests that the fish coats itself with anemone mucus to mask its chemical stimulus for cnida discharge. The second suggests that the fish is protected by some alteration in its own mucous coating.

In order to test these alternative hypotheses, a surrogate anemone was constructed and an Amphiprion allowed to associate with it prior to being placed with a real anemone. If the subsequent acclimation times to a real anemone are reduced significantly after exposure to the surrogate anemone, one would conclude that the fish does not need to coat itself with anemone mucus to become protected.

After association with a surrogate anemone, it was found that the acclimation times of isolated A. clarkii to the sea anemone Macrodactyla doreensis (Quoy & Gaimard) decreased seven-fold. These results support the hypothesis that the fish is responsible for manufacturing its own protection from sea anemones, presumably by independently altering its own mucous coat during acclimation.     

Complete amino acid sequence of tenebrosin-C, a cardiac stimulatory and haemolytic protein from the sea anemone Actinia tenebrosa

The complete amino acid sequence of the cardiac stimulatory and haemolytic protein tenebrosin-C, from the Australian sea anemone Actinia tenebrosa, has been determined by Edman degradation of the intact molecule and fragments produced by treatment of the polypeptide chain with cyanogen bromide and enzymatic cleavage with endoproteinase Asp-N, thermolysin and trypsin. The molecule is a single-chain polypeptide consisting of 179 amino acid residues with a calculated molecular mass of 19,797 Da. Tenebrosin-C shows a high degree of amino acid sequence similarity (63%) with Stoichactis helianthus cytolysin III [Blumenthal, K. M. and Kem, W. R. (1983) J. Biol. Chem. 258, 5574-5581] and is identical to a partial sequence (90 residues) reported for equinatoxin, a cardiostimulatory and haemolytic protein isolated from the European sea anemone Actinia equina [Ferlan, I. and Jackson, K. (1983) Toxicon Suppl. 3, 141-144]. No amino acid sequence similarity was detected between tenebrosin-C and other protein sequences stored in available databases. The predicted secondary structure of tenebrosin-C suggests that it is a compact, highly structured molecule.