Elongation Factor 1-Alpha Sequences Do Not Support an Early Divergence of the Acoela

Cédric Berney, Jan PawlowskiGo, and Louisette Zaninetti

Département de Zoologie et Biologie Animale, Université de Genève, Switzerland


    Abstract
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 
The phylogenetic position of the Acoela is a key problem in the understanding of metazoan evolution. Recent studies based on 18S ribosomal DNA (rDNA) sequences have placed the Acoela in an extremely basal position as the sister group to all other extant triploblastic animals, suggesting that the phylum Platyhelminthes is polyphyletic. In order to test the results obtained with 18S rDNA, we sequenced elongation factor 1-alpha (EF1a) for the acoel Convoluta roscoffensis and five species of Turbellaria (two polyclads, Leptoplana tremellaris, and Prostheceraeus vittatus, and three triclads, Crenobia alpina, Schmidtea polychroa, and Girardia tigrina). Phylogenetic analyses of EF1a sequences show that the acoel sequences branch within the Platyhelminthes, in opposition to the 18S rDNA data. Moreover, comparison of the central variable region of EF1a shows similar sequence signatures between C. roscoffensis and the three triclad species. Although EF1a sequences fail to prove the monophyly of the phylum Platyhelminthes, they do not confirm the early divergence of the Acoela.


    Introduction
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 
Acoel flatworms, ordinarily classified in the phylum Platyhelminthes as an order of the class Turbellaria, have recently been given great attention because of their particularly problematic phylogenetic position. Recent studies based on 18S ribosomal DNA (rDNA) sequences cast doubt on the monophyly of the phylum Platyhelminthes (Carranza, Baguña, and Riutort 1997Citation ) and suggested a possible primitive position of acoel flatworms, making them the earliest-diverging Bilateria and the sister group to all other extant triploblastic animals (Ruiz-Trillo et al. 1999Citation ).

The 18S rDNA gene is the most commonly used gene in phylogenetic analyses, and its database is the most complete for all living organisms (Hillis and Dixon 1991Citation ; Maidak et al. 1999Citation ). Like other ribosomal genes, the 18S gene is very easy to amplify because of its large number of copies per genome. Furthermore, the alternation in a single molecule of very conserved and much more variable regions makes 18S rDNA a powerful tool for molecular phylogeny at different taxonomic levels. However, it has been shown that in many phyla, such as Rhombozoa, Dicyemida, Nematoda, Gnathostomulida, and Chaetognatha, 18S rDNA genes do not evolve at a constant rate (Katayama et al. 1995Citation ; Winnepenninckx et al. 1995Citation ; Pawlowski et al. 1996Citation ; Littlewood et al. 1998Citation ). As a result, phylogenetic analyses lead to an artifactual basal clustering of the fast-evolving phyla, known as the long-branch attraction (LBA) phenomenon. The Acoela, most of which have been shown to be fast-clock organisms for the 18S rDNA gene (Carranza, Baguña, and Riutort 1997Citation ), represent a particularly good example of the LBA effect.

There are two main ways to overcome the LBA problem. The first one is to broadly sample problematic groups in order to find non-fast-clock species among each of them. This method was successfully used with the phylum Nematoda, shown to be closely related to the Arthropoda (Aguinaldo et al. 1997Citation ), and led to a better comprehension of general metazoan phylogeny, with the separation of the triploblastic animals into three main clades, the Lophotrochozoa, the Ecdysozoa, and the Deuterostomia (Halanych et al. 1995Citation ; Aguinaldo et al. 1997Citation ). Ruiz-Trillo et al. (1999)Citation applied this method to the Acoela and retained the unique non-fast-clock sequence of Paratomella rubra out of 18 acoel species. When only taxa with uniform rates of change were included in their analyses, the sequence of P. rubra appeared in a very basal position, as the first offshoot of the triploblastic animals. Therefore, the authors concluded that the acoel flatworms do not belong to the phylum Platyhelminthes and are probably the closest living relatives of the first Bilateria.

A second way to avoid the LBA effect is to find another gene for which rates of substitution are more constant between all metazoan phyla. The problem in finding such a gene is that no other database of metazoan sequences is as complete as the one for 18S rDNA. This is the main reason why so few studies of animal phylogeny are based on other genes, despite the general agreement that molecular phylogenetic inference requires congruent results from multiple gene sequences to be really conclusive (Baldauf and Palmer 1993Citation ; Regier and Schultz 1997Citation ; Abouheif, Zardoya, and Meyer 1998Citation ).

Elongation factor 1-alpha (EF1a) is a protein found in all eukaryotic lineages. EF1a is very abundant in the cytosol, where it is involved in the GTP-dependent binding of aminoacyl-tRNAs to the A site of the ribosomes in the second step of translation from mRNAs to proteins. Its universal occurrence and good amino acid sequence conservation makes it an ideal phylogenetic tool with which to determine very ancient relationships, such as the relative branching order of the most primitive Eukaryotes (see, e.g., Hashimoto et al. 1994Citation ; Nordnes, Krauss, and Johansen 1994Citation ; Baldauf, Palmer, and Doolittle 1996Citation ). Furthermore, the nucleotide sequences of EF1a genes were used at a much lower level to study phylogenies of families and genera, using the degenerated third position of codons (see, e.g., Cho et al. 1995Citation ; Mitchell et al. 1997Citation ). EF1a has also been used for resolving the phylogenetic relationships between animal phyla and classes as a good means to confirm or cast doubt on results based on other genes (Kojima et al. 1993Citation ; Regier and Shultz 1997Citation ; Kojima 1998Citation ).

In this study, we use EF1a in the first attempt to infer the phylogenetic position of the Acoela using a protein-coding gene. A fragment of about 950 nt of EF1a was amplified and sequenced for the acoel flatworm Convoluta roscoffensis, as well as for three species of Tricladida (Crenobia alpina, Schmidtea polychroa, and Girardia tigrina) and two species of Polycladida (Prostheceraeus vittatus and Leptoplana tremellaris). We compared it to other metazoan sequences in order to establish the phylogenetic position of the Acoela and test the validity of the results obtained with 18S rDNA sequences.


    Materials and Methods
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 
Sample Collection and DNA Extractions
Living specimens of C. roscoffensis were collected on the west coast of Cotentin (Normandy, France) and kept alive for several months in the laboratory. Other specimens were collected in Roscoff (Brittany, France) and preserved in 70% alcohol. Two members of the order Polycladida (L. tremellaris and P. vittatus) were collected near the marine biology station of Luc-sur-Mer (Normandy, France) and also preserved in 70% alcohol. Living specimens of S. polychroa and G. tigrina (Tricladida, Dugesiidae) were collected in Geneva, Switzerland, and living specimens of C. alpina (Tricladida, Planariidae) were collected above Martigny, Switzerland. They were all kept alive in an aquarium.

DNA from one individual of each species was extracted in guanidine lysis buffer, precipitated with isopropanol, and dissolved in distilled water (Tkach and Pawlowski 1999Citation ).

PCR Amplification, Cloning, and Sequencing
PCR amplifications were all performed in a total volume of 50 µl with an amplification profile consisting of 40 cycles of 30 s at 94 °C, 30 s at 52 °C, and 120 s at 72 °C, followed by 5 min at 72 °C for final extension. To avoid contamination with the endosymbiotic algae of C. roscoffensis, a set of four metazoan-specific primers was designed for the PCR-amplification of a fragment of about 950 nt of EF1a (see table 1 for the sequences and positions of these primers). These metazoan-specific primers were used to amplify the same fragment of EF1a for the two polyclads and the three triclads. All PCR products were purified using the High Pure PCR purification kit (Roche) and then ligated in the p-GEM-T vector (Promega) and cloned in XL-2 Ultracompetent Cells (Stratagene), all according to the manufacturers' instructions. The clones were then sequenced on an ABI 377 Prism sequencer. The sequences reported in this paper have been deposited in the GenBank database under accession numbers AJ250908AJ250914.


View this table:
[in this window]
[in a new window]
 
Table 1 List of the Primers Used in this Study

 
Phylogenetic Analyses
The sequences of EF1a obtained in this study were manually aligned to other available metazoan and fungal sequences using the Genetic Data Environment software, version 2.2 (Larsen et al. 1993Citation ). All species names, taxonomic ranks, and GenBank accession numbers of EF1a sequences used in our analyses are listed in table 2 . A relative-rate test was performed using RRTree (Robinson-Rechavi and Huchon 2000).


View this table:
[in this window]
[in a new window]
 
Table 2 List of the Species Used in this Study, Along with GenBank Accession Numbers

 
An evolutionary tree was inferred from the deduced amino acid sequences using the neighbor-joining (NJ) method (Saitou and Nei 1987Citation ) applied to Dayhoff's PAM distance matrix (Dayhoff, Schwartz, and Orcutt 1978Citation ). The reliability of internal branches was assessed using the bootstrap method (Felsenstein 1988Citation ) with 1,000 bootstrap replicates. The Phylo-Win program (Galtier and Gouy 1996Citation ) was used for distance computations and NJ building and bootstrapping. All maximum-parsimony (MP) analyses were performed with PAUP* (Swofford 1998Citation ). The most parsimonious trees were determined using a heuristic search procedure with 100 random-addition-sequence replicates and tree bisection-reconnection branch swapping. Reliability of internal branches was evaluated with 500 bootstrap replicates. All maximum-likelihood (ML) analyses were carried out with the PUZZLE software, version 4.02 (Strimmer and von Haeseler 1996Citation ), except for the exhaustive tree searches using constrained tree topologies, which were performed with the protml program of the MOLPHY package, version 2.3 (Adachi and Hasegawa 1996Citation ).


    Results
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 
The amplified fragment of EF1a varies in size, from 947 nt (316 amino acids) in the two polyclads to 962 nt (321 amino acids) in C. roscoffensis to 965 nt (322 amino acids) in the three triclads. The mean length of vertebrate EF1a is 1386 nucleotides (462 amino acids). The amplified EF1a of C. roscoffensis and the five turbellarians corresponds to amino acid positions 45–374 of human EF1a. The two polyclad sequences share an intron of about 110 nt in the codon corresponding to amino acid 343. The two sequences of C. roscoffensis obtained from two different DNA extractions, the first (1) from an individual from Roscoff and the second (2) from an individual from the west coast of Cotentin, differ in five positions, of which four are silent transitions in the third codon positions and one is a transition in the second position which changes glycine to glutamic acid at position 230, as shown in figure 1 .



View larger version (60K):
[in this window]
[in a new window]
 
Fig. 1.—Alignment of the 12-amino-acid insertion shared by fungi and metazoans in the central region of EF1a, showing the five additional amino acids of Convoluta roscoffensis and the six additional amino acids of Crenobia alpina, Schmidtea polychroa, and Girardia tigrina. The insertion begins with a conserved tryptophan and ends with a conserved glycine. Convoluta roscoffensis and the three triclads share the common sequence signature KK(ED)E. Numbers at bottom indicate the positions of the corresponding amino acid residues in human EF1a (GenBank accession number J04617). The only difference between the two amino acid sequences of C. roscoffensis is visible at position 230

 
The seven sequences obtained were aligned to diverse eukaryotic sequences (see table 2 ). The alignment shows that the variation in size between C. roscoffensis and the different members of the Turbellaria occurs in the central variable region of EF1a. This region consists of an insertion of about 12 amino acids (positions 214–225) shared only by metazoans and fungi, as shown by Baldauf and Palmer (1993)Citation . Whereas the two polyclad sequences are of the same length as those of fungi and most animals, the sequences of C. roscoffensis and the three triclads are five and six amino acids longer, respectively. The six additional amino acids of the triclad sequences can be quite readily aligned to the five additional amino acids of the sequences of C. roscoffensis, because they share the clear amino acid sequence signature KK(ED)E (fig. 1 ).

Thirty-eight sequences of diverse metazoan and fungal species were retained for the phylogenetic analyses, with the latter being used as a designated outgroup in all analyses. The choice of Fungi as an outgroup is congruent with previous studies based on several proteins (Baldauf and Palmer 1993Citation ; Nikoh et al. 1994Citation ) which confirmed that they are the closest eukaryotic relatives of the Metazoa. Sequences of some groups (nematodes, molluscs, echinoderms) were not selected because of their artifactual branching at the base of the Metazoa. Three hundred two amino acid residues, corresponding to the amplified nucleotide fragments excluding the primers and the short insert in C. roscoffensis and the triclads, were used in the analyses (the complete alignment is available at our web site, http://www.unige.ch/sciences/biologie/biani/msg). Of these 302 amino acids, 149 are constant (48.2%). Twenty-three variable positions are parsimony-uninformative, and 130 are parsimony-informative.

Prior to phylogenetic analyses, we performed two tests. First, we ran a relative-rate test (Robinson-Rechavi and Huchon 2000) on both nucleotide and amino acid sequences of EF1a to determine if this gene shows a relatively constant rate of substitution among lineages. The test confirmed at a 5% level that the differences in substitution rates among all lineages—including the two acoel sequences—were nonsignificant; in all lineages, the EF1a's evolve at similar rates. Then, we performed a likelihood-mapping analysis (Strimmer and von Haeseler 1997Citation ) to evaluate the phylogenetic content of the data. Likelihoods of all 148,995 possible quartets of sequences were mapped on a triangle, which was partitioned in seven regions. Of these quartets, 84.3% had completely resolved phylogenies, and only 7.8% were in the star-tree region (fig. 2 ), indicating that the 302 amino acid residues of EF1a used in our analyses contain a satisfying degree of phylogenetic information. NJ, MP, and ML analyses were then performed.



View larger version (15K):
[in this window]
[in a new window]
 
Fig. 2.—Likelihood-mapping analysis of the data set, represented as a triangle. Values at the corners indicate the percentages of well-resolved phylogenies for all 148,995 possible quartets of sequences, and values at the central and lateral regions are the percentages of unresolved phylogenies. The cumulatively high percentage (84.3 %) from the corner values indicates that the data set is phylogenetically informative

 
The phylogenetic tree inferred from the amino acid sequence of EF1a with the NJ method, using Dayhoff's PAM distance matrix, is shown in figure 3A. In the tree, the three species of diploblasts are grouped together and form the sister group to the triploblastic animals. Vertebrates and arthropods form monophyletic groups. Annelids plus pogonophorans also appear as a strongly supported monophyletic unit, with the two pogonophorans and the three clitellates forming two groups independently derived from polychaetes, as shown by McHugh (1997)Citation and Kojima (1998). The two sequences of C. roscoffensis and the five turbellarian species are grouped together. Convoluta roscoffensis fall within the sequences of Turbellaria, as the sister group to the polyclad species, with a bootstrap support value of 56%. The three triclad species are placed at the base of the group.



View larger version (87K):
[in this window]
[in a new window]
 
Fig. 3.—Phylogenetic position of the acoel flatworm Convoluta roscoffensis inferred from amino acid sequences of EF1a using (A) the neighbor-joining (NJ) method and (B) the maximum-parsimony method. In both trees, the two polyclads, the three triclads, and Convoluta roscoffensis appear as a monophyletic unit, a sister group to members of the Lophotrochozoa. The trees are rooted using five fungal sequences as a designated outgroup. A, NJ tree constructed using Dayhoff's PAM distance matrix, with all branch lengths drawn to scale. Numbers at internal nodes indicate percentages of 1,000 bootstrap replicates that support the branch; values less than 50% are omitted. B, Fifty percent majority-rule consensus of the eight most- parsimonious trees (percentiles less than 100% are indicated below the branches) obtained using a heuristic search procedure with 100 random-addition-sequence replicates and tree bisection-reconnection branch swapping. Trees length is 1,025 steps, with a consistency index of 0.447 and a retention index of 0.583. Bootstrap values of 500 replicates are indicated above the nodes; values less than 50% are omitted. The tree differs from the NJ tree mainly in the position of C. roscoffensis as the sister group to the three triclad species

 
The 50% majority-rule consensus tree of the eight most- parsimonious trees found with the MP method is shown in figure 3B. It has a length of 1,025 steps, with a consistency index (CI) of 0.447 and a retention index (RI) of 0.583. Four clades are present in both MP and NJ trees—the diploblastid clade, the vertebrate clade, the arthropod clade, and a lophotrochozoan clade formed by platyhelminths, annelids, and vestimentiferans. Their relative branching order is not resolved. All groups supported by bootstrap values greater than 70% in the NJ tree are also well supported in the MP analysis, except for the arthropod clade, which is supported by only a 47% bootstrap value. As in the NJ tree shown in figure 3A, the Platyhelminthes cluster together, but the main difference between the two trees is that the two sequences from C. roscoffensis are placed in the MP tree as the sister group to the three triclads with a 50% bootstrap value, with the two polyclads branching off at the base of the group.

Next, an ML tree was constructed with the quartet puzzling (QP) algorithm (Strimmer and von Haeseler 1996Citation ), which automatically assigns estimations of support to each internal branch. Dayhoff's PAM substitution model was used. The frequency of each amino acid was estimated from the data set, and a 5% level chi-square test was performed to confirm that the amino acid composition of each sequence was identical to the average amino acid composition of the whole alignment. The chosen model of rate heterogeneity was a discrete gamma distribution with eight categories, with all necessary parameters being estimated automatically from the data set. The QP search was conducted with 10,000 puzzling steps. As shown by the likelihood-mapping analysis, 11,592 of the 148,995 quartets analyzed (7.8%) were unresolved, leading to a multifurcating tree with a basal polytomy in the metazoan clade. The tree has a log likelihood value of -6,507.07. Six groups of metazoans are defined: the diploblastid clade, the vertebrate clade, the arthropod clade, the clade formed by annelids plus pogonophorans, and two platyhelminth clades. As in the NJ tree, the two sequences of C. roscoffensis appear as the sister group to the two polyclads, with a QP reliability of 65%. However, the monophyly of the Platyhelminthes was recovered in only 28% of all intermediate trees generated in the 10,000 puzzling steps, and the three triclad species cluster in a distinct group (data not shown).


    Discussion
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 
The results of our study are in opposition to 18S rDNA analyses, which place the acoel flatworms as the sister group to all other extant triploblastic animals, well apart from other Platyhelminthes (Ruiz-Trillo et al. 1999Citation ). Analyses of EF1a provide two independent pieces of evidence favoring a close relationship between the Acoela and the Turbellaria. First, C. roscoffensis shares an amino acid sequence signature with the three triclad species, which can be interpreted as a single insertion event in a common ancestor of both lineages (Baldauf and Palmer 1993Citation ). Second, in all types of phylogenetic analyses, Convoluta always branches as the sister group to members of the Turbellaria. This finding is in agreement with morphological data. The simplicity of acoel morphological features can be interpreted as secondary derived reductions (Ehlers 1985Citation ; Smith and Tyler 1985Citation ; Smith, Tyler, and Rieger 1986Citation ), and there are many rigorous morphological synapomorphies that support a sister group relationship between the Acoela and some members of the Turbellaria, i.e., the Nemertodermatida (Ehlers 1985Citation ; Smith and Tyler 1985Citation ; Ax 1996Citation ). Furthermore, the phylogenetic position of the group formed by C. roscoffensis and the Turbellaria in both NJ and MP analyses, placed as the sister group to other members of the Lophotrochozoa (i.e., annelids and pogonophorans), is in agreement with previous studies on the phylogenetic position of the phylum Platyhelminthes (Aguinaldo et al. 1997Citation ; Balavoine 1997Citation ).

However, the monophyly of the Platyhelminthes was never well supported in our analyses, as shown by very low bootstrap values (8% in the NJ analysis and 21% in the MP analysis). Interestingly, this low support for the platyhelminth clade is apparently mainly due to the important differences between polyclad and triclad sequences. Indeed, if the acoel sequences are not included in the analysis, the five Turbellaria no longer appear monophyletic with the NJ method (data not shown), with the triclads being the sister group to the annelid plus pogonophoran clade. On the other hand, when the three triclad sequences are not included in the analysis, the grouping of the acoel sequences with the two polyclads is strengthened, with a bootstrap support value of 74% (data not shown). Inversely, when the two polyclad sequences are removed from the analysis, the two acoel sequences group with the triclads with a bootstrap support value of 61% (data not shown). Similarly, with the MP method, removing the two polyclad sequences from the analysis strengthens the support value for the acoel plus triclad clade (71%, data not shown). In addition, if the triclad sequences are removed from the analysis, the two sequences of C. roscoffensis cluster with the polyclad sequences with a bootstrap support value of 62%. In all cases, the Acoela branch with one of the lineages of Turbellaria, and the platyhelminth clade always remains the sister group to the other lophotrochozoan species. In fact, the two sequences of C. roscoffensis share several common amino acid substitutions with the polyclads on the one hand and with the triclads on the other hand, with very few of them being common to the five Turbellaria. This may explain why the monophyly of the platyhelminth clade is never well supported, whereas the monophyly of the acoel plus triclad clade alone and that of the acoel plus polyclad clade alone are supported by quite confident bootstrap values.

Given the weak support of EF1a data for the monophyly of the Platyhelminthes, we were interested in testing their support for the hypothesis of an early branching of the Acoela, suggested by 18S rDNA data. With this aim, we used the exhaustive search mode of the protml program in MOLPHY, version 2.3 (Adachi and Hasegawa 1996Citation ), which allows constrained tree topologies, to determine the likelihood values of the totally resolved trees compatible with the six metazoan groups present in the QP tree. We retained only the 20 best tree topologies out of the 945 that were compatible with the QP analysis. Nine of them show the lophotrochozoan clade to be monophyletic, and in four, the Platyhelminthes are monophyletic and a sister group to the annelid plus pogonophoran clade. Moreover, the three diploblastid sequences are the first offshoot of the Metazoa in 18 of the 20 best trees. The Kishino-Hasegawa test (Kishino and Hasegawa 1989Citation ) showed at a 5% level that none of the 20 best topologies could be rejected as significantly worse than the other. Next, we performed the same test on two constrained topologies, defined as follows: we used the general topology obtained by Ruiz-Trillo et al. (1999)Citation as a constrained topology, first with the Platyhelminthes monophyletic, then with the two acoel sequences at the base of all other triploblastic species. The first constrained topology was included in the 20 best tree topologies; it has a log likelihood value of -6,373.49. However the topology which shows the Acoela as the first offshoot after the diploblasts, which has a log likelihood value of -6,392.86, was rejected as significantly worse than the best tree at a 5% level. This test shows that even if EF1a sequences cannot resolve the problem of the monophyly of the phylum Platyhelminthes, they give good evidence in favor of rejecting the hypothesis that acoel flatworms represent the first offshoot of the Bilateria.

A problematic feature of our analyses is the fact that the two sequences of C. roscoffensis grouped with the polyclad sequences in both NJ and ML analyses, whereas they grouped with the triclad sequences in the MP analysis. This can be explained by the fact that in both NJ and ML analyses, amino acid substitutions were weighted using Dayhoff's PAM substitution model, whereas in the MP analysis, all amino acid substitutions were equally weighted. As stated above, the two acoel sequences share a few more amino acid substitutions with the three triclads than with the two polyclads, leading to the grouping of C. roscoffensis with the triclads in the MP analysis. However, when Dayhoff's PAM substitution model is taken into account, these substitutions are weighted differently, and C. roscoffensis groups with the polyclads. Nevertheless, the sequence signature shared by C. roscoffensis and the three triclads suggests that the Acoela might be more closely related to the Tricladida than to the Polycladida. Interestingly, Campos et al. (1998)Citation also found a close relationship between Tricladida and Acoela in their molecular phylogeny of the phylum Platyhelminthes, based on partial 18S rDNA sequences. In the analyzed fragment of EF1a, about 15 specific amino acid substitutions define the triclad clade, whereas only about 10 define the polyclad clade. NJ and ML analyses might thus be misleading because of the numerous amino acid substitutions that define the triclad clade, which tend to group the acoel and polyclad sequences.

More sequences from representatives of all turbellarian orders and of members of the Neodermata would perhaps help to clarify this situation. It has been shown (Lecointre et al. 1993Citation ; Giribet and Ribera 1998Citation ) that increasing the number of taxa in a clade can help to strengthen its monophyly, although it also very rapidly increases the duration of the analyses. However, even more EF1a sequences may not definitively resolve the relationships within the Platyhelminthes. As shown by our analyses, the phylogeny of the phyla for which more EF1a sequences are available is not always well resolved. Whereas the internal topology of the arthropod clade is always well defined in all trees, the one of the annelid plus pogonophoran clade is much vaguer. This shows that the phylogenetic information contained in EF1a sequences is not the same for all metazoan groups. EF1a is thus probably not the perfect gene for inferring good metazoan phylogenies, but neither is 18S rDNA, as shown by Abouheif, Zardoya, and Meyer (1998)Citation . Both genes give conflicting results, which prompts us to search for another molecular marker to conclusively resolve the phylogenetic position of Acoela. Our study stresses the need for the use of several different genes to properly determine the phylogenetic position of any group of organisms.


View this table:
[in this window]
[in a new window]
 
Table 2 Continued

 

    Acknowledgements
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 
We would like to thank José Fahrni for sequencing assistance and Maria Holzmann, Juan Montoya, Colomban de Vargas, and Olivier Jousson for helpful comments on our work.


    Footnotes
 
Manolo Gouy, Reviewing Editor

1 Abbreviations: CI, consistency index; EF1a, elongation factor 1-alpha; LBA, long-branch attraction; NJ, neighbor joining; ML, maximum likelihood; MP, maximum parsimony; QP, quartet puzzling; rDNA, ribosomal DNA; RI, retention index. Back

2 Keywords: elongation factor 1-alpha molecular phylogeny Convoluta roscoffensis, Acoela Platyhelminthes Back

3 Address for correspondence and reprints: Jan Pawlowski, Université de Genève, Station de Zoologie, 154, route de Malagnou, 1224 Chêne-Bougeries, Genève, Switzerland. E-mail: jan.pawlowski{at}zoo.unige.ch Back


    literature cited
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 literature cited
 

    Abouheif, E., R. Zardoya, and A. Meyer. 1998. Limitations of metazoan 18S rRNA sequence data: implications for reconstructing a phylogeny of the animal kingdom and inferring the reliability of the Cambrian explosion. J. Mol. Evol. 47:394–405.[ISI][Medline]

    Adachi, J., and M. Hasegawa. 1996. MOLPHY version 2.3: programs for molecular phylogenetics based on maximum likelihood. Comput. Sci. Monogr. 28:1--150.

    Aguinaldo, A. M. A., J. M. Turbeville, L. S. Linford, M. C. Rivera, J. R. Garey, R. A. Raff, and J. A. Lake. 1997. Evidence for a clade of Nematodes, Arthropods and other moulting animals. Nature 387:489–493.

    Ax, P. 1996. Multicellular animals. A new approach to the phylogenetic order in nature. Springer Verlag, Berlin.

    Balavoine, G. 1997. The early emergence of platyhelminths is contradicted by the agreement between 18S rRNA and Hox data. C. R. Acad. Sci. Paris Life Sci. 320:83–94.

    Baldauf, S. L., and J. D. Palmer. 1993. Animals and Fungi are each other's closest relatives: congruent evidence from multiple proteins. Proc. Natl. Acad. Sci. USA 90:11558–11562.

    Baldauf, S. L., J. D. Palmer, and W. F. Doolittle. 1996. The root of the universal tree and the origin of eukaryotes based on elongation factor phylogeny. Proc. Natl. Acad. Sci. USA 93:7749–7754.

    Campos, A., M. P. Cummings, J. L. Reyes, and J. P. Laclette. 1998. Phylogenetic relationships of Platyhelminthes based on 18S ribosomal gene sequences. Mol. Phylogenet. Evol. 10:1–10.[ISI][Medline]

    Carranza, S., J. Baguña, and M. Riutort. 1997. Are the Platyhelminthes a monophyletic primitive group? An assessment using rRNA sequences. Mol. Biol. Evol. 14:485–497.[Abstract]

    Cho, S., A. Mitchell, J. C. Regier, C. Mitter, R. W. Poole, T. P. Friedlander, and S. Zhao. 1995. A highly conserved nuclear gene for low-level phylogenetics: elongation factor-1alpha recovers morphology-based tree for Heliothine moths. Mol. Biol. Evol. 12:650–656.[Abstract]

    Dayhoff, M. O., R. M. Schwartz, and B. C. Orcutt. 1978. A model of evolutionary change in proteins. Pp. 345–352 in M. O. Dayhoff, ed. Atlas of protein sequence Structure. Vol. 5, Suppl. 3. National Biomedical Research Foundation, Washington, D.C.

    Ehlers, U. 1985. Das phylogenetische System der Plathelminthes. Fischer Verlag, Stuttgart.

    Felsenstein, J. 1988. Phylogenies from molecular sequences: inference and reliability. Annu. Rev. Genet. 22:521–565.[ISI][Medline]

    Galtier, N., and M. Gouy. 1996. SEAVIEW and PHYLO_WIN: two graphic tools for sequence alignment and molecular phylogeny. Comput. Appl. Biosci. 12:543–548.[Abstract]

    Giribet, G., and C. Ribera. 1998. The position of Arthropods in the animal kingdom: a search for a reliable outgroup for internal arthropod phylogeny. Mol. Phylogenet. Evol. 9:481–488.[ISI][Medline]

    Halanych, K. M., J. D. Bacheller, A. M. A. Aguinaldo, S. M. Liva, D. M. Hillis, and J. A. Lake. 1995. Evidence from 18S ribosomal DNA that the Lophophorates are protostome animals. Science 267:1641–1643.

    Hashimoto, T., Y. Nakamura, F. Nakamura, T. Shirakura, J. Adachi, N. Goto, K. Okamoto, and M. Hasegawa. 1994. Protein phylogeny gives a robust estimation for early divergences of Eukaryotes: phylogenetic place of a mitochondrion-lacking Protozoan, Giardia lamblia. Mol. Biol. Evol. 11:65–71.[Abstract]

    Hillis, D. M., and M. T. Dixon. 1991. Ribosomal DNA: molecular evolution and phylogenetic inference. Q. Rev. Biol. 66:411–453.[Medline]

    Katayama, T., H. Wada, H. Furuya, N. Satoh, and M. Yamamoto. 1995. Phylogenetic position of the Dicyemid Mesozoa inferred from 18S rDNA sequences. Biol. Bull. 189:81–90.[Abstract/Free Full Text]

    Kishino, H., and M. Hasegawa. 1989. Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order in Hominoidea. J. Mol. Evol. 29:170–179.[ISI][Medline]

    Kojima, S. 1998. Paraphyletic status of Polychaeta suggested by phylogenetic analysis based on the amino acid sequences of elongation factor-1{alpha}. Mol. Phylogenet. Evol. 9:255–261.[ISI][Medline]

    Kojima, S., T. Hashimoto, M. Hasegawa, S. Murata, S. Ohta, H. Seki, and N. Okada. 1993. Close phylogenetic relationships between Vestimentifera (tube worms) and Annelida revealed by the amino acid sequence of elongation factor-1alpha. J. Mol. Evol. 37:66–70.[ISI][Medline]

    Larsen, N., G. J. Olsen, B. L. Maidak, M. J. McCaughey, R. Overbeek, T. J. Macke, T. L. Marsh, and C. R. Woese. 1993. The ribosomal database project. Nucleic Acids Res. 21:3021–3023.[Abstract]

    Lecointre, G., H. Philippe, H. L. Van Le, and H. Le Guyader. 1993. Species sampling has a major impact on phylogenetic inference. Mol. Phylogenet. Evol. 2:205–224.[Medline]

    Littlewood, D. T. J., M. J. Telford, K. A. Clough, and K. Rohde. 1998. Gnathostomulida—an enigmatic metazoan phylum from both morphological and molecular perspectives. Mol. Phylogenet. Evol. 9:72–79.[ISI][Medline]

    McHugh, D. 1997. Molecular evidence that echiurans and pogonophorans are derived annelids. Proc. Natl. Acad. Sci. USA 94:8006–8009.

    Maidak, B. L., J. R. Cole, C. T. Parker Jr. et al. (14 co- authors). 1999. A new version of the RDP (Ribosomal Database Project). Nucleic Acids Res. 27:171–173.[Abstract/Free Full Text]

    Mitchell, A., S. Cho, J. C. Regier, C. Mitter, R. W. Poole, and M. Matthews. 1997. Phylogenetic utility of elongation factor-1alpha in Noctuoidea (Insecta: Lepidoptera): the limits of synonymous substitution. Mol. Biol. Evol. 14:381–390.[Abstract]

    Nikoh, N., N. Hayase, N. Iwabe, K. Kuma, and T. Miyata. 1994. Phylogenetic relationships of the kingdoms Animalia, Plantae and Fungi, inferred from 23 different protein species. Mol. Biol. Evol. 11:762–768.[Abstract/Free Full Text]

    Nordnes, S., S. Krauss, and T. Johansen. 1994. cDNA sequence of zebrafish (Brachydanio rerio) translation elongation factor 1-alpha: molecular phylogeny of eukaryotes based on elongation factor 1-alpha protein sequences. Biochim. Biophys. Acta 1219:529–532.

    Pawlowski, J., J.-I. Montoya-Burgos, J. F. Fahrni, J. Wüest, and L. Zaninetti. 1996. Origin of the Mesozoa inferred from 18S rRNA gene sequences. Mol. Biol. Evol. 14:498–505.[Abstract]

    Regier, J. C., and J. W. Shultz. 1997. Molecular phylogeny of the major Arthropod groups indicates polyphyly of Crustaceans and a new hypothesis for the origin of Hexapods. Mol. Biol. Evol. 14:902–913.[Abstract]

    Robinson-Rechavi, M., and D. Huchon. 2000. RRTree: relative rate tests between groups of sequences on a phylogenetic tree. Bioinformatics (in press).

    Ruiz-Trillo, I., M. Riutort, D. T. J. Littlewood, E. A. Herniou, and J. Baguña. 1999. Acoel flatworms: earliest extant bilaterian Metazoans, not members of Platyhelminthes. Science 283:1919–1923.

    Saitou, N., and M. Nei. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4:406–425.[Abstract]

    Smith, J. P. S. III, and S. Tyler. 1985. The acoel turbellarians: kingpins of metazoan evolution or specialized offshot? Pp. 123–142 in S. C. Morris, J. D. George, R. Gibson, and H. M. Platt, eds. The origins and relationships of lower invertebrates. Oxford University Press, Oxford, England.

    Smith, J. P. S. III, S. Tyler, and R. M. Rieger. 1986. Is the Turbellaria polyphyletic? Hydrobiologia 132:13–21.

    Strimmer, K., and A. von Haeseler. 1996. Quartet puzzling: a quartet maximum likelihood method for reconstructing tree topologies. Mol. Biol. Evol. 13:964–969.[Free Full Text]

    ———. 1997. Likelihood-mapping: a simple method to visualize phylogenetic content of a sequence alignment. Proc. Natl. Acad. Sci. USA 94:6815–6819.

    Swofford, D. L. 1998. PAUP*, phylogenetic analysis using parsimony (*and other methods). Version 4.0b2a. Sinauer, Sunderland, Mass.

    Tkach, V., and J. Pawlowski. 1999. A new method of DNA extraction from the ethanol-fixed parasitic worms. Acta Parasitol. 44:147–148.[ISI]

    Winnepenninckx, B., T. Backeljau, L. Y. Mackey, J. M. Brooks, R. De Wachter, S. Kumar, and J. R. Garey. 1995. 18S rRNA data indicate that Aschelminthes are polyphyletic in origin and consist of at least three distinct clades. Mol. Biol. Evol. 12:1132–1137.[Abstract]

Accepted for publication March 2, 2000.